1. Introduction
In general, a typical rectenna, i.e., rectifying antenna, consists of an antenna and rectifying circuit. The rectenna was first developed by Brown in 1963 [
1]. The rectifying circuit is used to convert electromagnetic (EM) waves into direct current (DC) electricity. The conventional rectifying circuits work properly with radio-frequency (RF) signals of a range from 3 kHz to 500 kHz. [
2] Certain carefully-designed rectifying circuits could be operating at up to 300 MHz with bipolar junction transistors and 2 GHz with integration of multiple MOS devices [
3,
4]. As the telecommunication is entering into the 5G and soon 6G era, both fast response and large bandwidth are needed, pushing the limit of frequency toward the high end of RF range (>100 GHz) [
5]. For example, millimeter waves (30-300 GHz) and terahertz signal will be used in 6G technology [
6,
7]. The 100 GHz signal demands 10 ps in response time, approaching the limit in transit time for an electron to diffuse across a silicon PN junction diode. Such a long diffusion time places a strict limit on the high-frequency application in converting alternative current (AC) into DC.
On the other hand, optical detection, as an important method for non-ionizing medical detections [
8,
9], has been increasingly used in wearable electronics, e.g., smart watches, to measure heartbeat, blood sugar, blood oxygen and so on. This technology requires the lights, a kind of high-frequency electromagnetic waves, to penetrate through the skin [
10]. The current wearable electronics are often equipped with visible lights, e.g., green and red lights, to perform the biomedical measurement [
11]. However, the medical information obtained by such wearable electronics is quite limited and not precise. THz propagation through the biological tissues would be accompanied by much less scattering loss because of its much longer wavelength compared to those of visible or infrared radiation [
12]. It is well known that EM waves at 10
12 Hz (THz) are able to image human tissues without harm and provide much more information at a higher resolution in comparison with other methods. For example, ultrasound, another method for non-ionizing detection, with frequencies ranging from 2 MHz to 12 MHz, has a resolution at 1 mm, much lower than that of optical and THz detection. Thus, THz electronics offer excellent safe evaluation, high resolution and precision, all of which are very attractive for future medical detections [
13,
14]. Since the optical and THz signal frequencies are too high for conventional semiconductor diodes, new diodes with a switching time faster than 1 ps are needed to rectify THz signal for medical detection.
In addition, rectenna has been studied for energy harvesting of 2.4 GHz EM waves for internet-of-things technology [
15]. To receive and rectify EM waves at optical frequencies, e.g., solar energy, optical rectennas [
16,
17] become attractive high-efficiency and low-cost solar cells if fast-switching diodes can be available. The conventional solar cells are based on a semiconductor PN junction whereas electron-hole pairs are generated by photons with energy larger or close to the energy bandgap (E
g) of the semiconductors. Silicon-based solar cells in the most commercially dominated products have a limited efficiency, about 16% for polycrystalline Si, 30% for monocrystalline Si and 55% for multi-junction cells [
18,
19,
20]. In comparison, the rectenna solar cells has a much higher efficiency (> 85%) and can be fabricated at low cost [
16].
Therefore, the common technical challenge for the above-mentioned emerging technologies (high-speed telecommunication, THz electronics and rectenna solar cells) hinges on new rectifier diodes. As shown in
Figure 1(a), in the conventional semiconductor PN junction diode, the excess carriers are injected to each side of the junction as minority carriers and diffuse through the diode. The diffusion time (τ
s) is close to the smallest value of electron lifetime (τ
e) or hole lifetime (τ
h). The lifetime of minority carrier depends on Auger lifetime due to Auger recombination [
21]. In typical, for Si with a doping of 10
18 cm
-3, the diffusion time is within a range from 10 ns to 1 μs, which is corresponding to switching frequencies from 1 MHz to 100 MHz. Optimization of device structures and diode materials, such as Schottky junction diodes and GaN diodes, has been studied to improve the switching frequency to GHz [
22]. However, THz and optical frequencies are still too high for the conventional PN junction diodes.
On the other hand, carrier transport based on tunneling can easily exceed the limit of diffusion in PN junctions, reaching the scale of femtoseconds (10
-15 s). The metal-insulator-metal (MIM) junction with an atomically thin insulator between two different metals is a simple yet very effective tunneling diode. A schematic of MIM tunneling diode and its energy barrier for tunneling are shown in Fig. 1(b). The electron tunneling time is determined by the energy barrier profile and applied bias across the barrier. Studies showed that the tunneling time across a tunneling diode could reach femtoseconds [
23] and even attoseconds (10
-18 s) in an atomic hydrogen. [
24]
The MIM heterostructures have recently been studied as a current rectifying diode [
25,
26,
27,
28]. Research interest of MIM diodes was focused on application in rectifying high-frequency electromagnetic waves, THz electronics, infrared light and even conversion of visible light into electricity [
16,
29,
30]. The MIM diodes based on bulk insulators, such as TiO
2, ZnO and NiO, could be applied in 28.3 THz rectenna switches [
31]. Amorphous InGaZn oxide was shown to improve turn-on voltage control in the diodes [
32]. Amorphous metal electrodes could also improve the current rectifying performance [
33]. Recently, MIM diodes based on 2D insulator materials have been reported [
34]. MIM diodes with single and multiple insulators have been studied and analyzed for the conversion of infrared light into electricity [
35,
36]. Despite the rapid progress, there still lacks a comprehensive study that combine first-principle calculation for designing 2D insulator materials, modeling of tunneling diodes, and simulation of rectifier circuits.
In this work, we applied density function theory (DFT) method to design and model MIM tunneling diodes with two-dimensional (2D) materials as the insulators. In the modeling, 1, 2, 3-layer TiO2, TaO2 and SnO2 were designed by oxidizing the corresponding 2D metal-sulfur materials. The monolayer or few-layer 2D materials were placed between two different metals to form a tunneling junction diode. In the device simulation, it is found that the Au/monolayer TiO2/Al tunneling junction diode exhibits the best rectifying performance. The simulated tunneling current-voltage (I-V) characteristics were fed through a rectifier in PSpice modeling, showing excellent signal rectifying from THz to optical frequencies.
2. Modeling and Simulation Method for Materials, Devices and Circuits
In this study, the atomistic simulations of 2D materials and MIM tunneling diodes were performed in QuantumWise Atomistix Toolkit (ATK) simulation platform [
37,
38,
39] from Synopsys based on density functional theory (DFT). The DFT calculations employed the norm-conserving Pseudo-Dojo pseudopotentials [
40] with Perdew-Burke-Ernzerhof (PBE) parametrization of exchange-correlation function [
41] for all materials and diode structures. To bring all atoms to the ground state, a 12×12×1 Monkhorst-Pack k-point grid mesh was used for the sampling of the Brillouin zone. The forces were kept to the least at 0.01 eV/Å. For the calculations of carrier transport in the tunneling diodes based on metal/2D insulators/metal heterostructures, we used first-principle calculations integrated with density functional non-equilibrium Green's function (NEGF) method. In addition, the current-voltage (I-V) characteristics was calculated by following the Landauer-Büttiker approach [
42]. The design of 2D materials and computation of carrier transport were implemented by using the above-mentioned QuantumWise ATK code. [
37]
In the simulation of 2D materials, 1, 2, 3-layer metal-sulfur materials are designed and modeled with structural relaxation. These ultra-thin 2D metal-sulfur materials are usually semiconductors with an energy bandgap. [
34] In the structural relaxation in simulation, the position of all atoms and the three-dimensional lattice parameters were optimized with the least Hellmann-Feynman force of 0.01 eV/Å on the metal-sulfur layer materials. In the optimization, Pulay-mixer algorithm was used as a fully self-consistent field (SCF) iteration control with a tolerance value of 10
-5 eV. The maximum number of iteration steps was set at a limit of 100. The self-consistent field computations were strictly tracked to guarantee full convergence within the iteration steps. Periodic boundary conditions were employed along the three directions (in-plane and out-of-plane) with a large vacuum region of 7 Å to minimize the interlayer interaction. [
43] Energy band structures and electronic properties of the metal-sulfur layer materials could be calculated once the structure optimization simulation was converged. In the design and simulation of metal-oxide layer materials, we proposed a simulated oxidation process: (1) the sulfur atoms were replaced with oxygen atoms; and (2) the new metal oxide materials were fully optimized by following the similar procedure in obtaining relaxed metal-sulfur layer materials.
The 2D layer oxide insulators were placed between two metal electrodes to construct metal/2D insulators/metal tunneling diodes. There is a 7-Å gap between the 2D insulators and electrodes at each side to protect the layer oxides from chemical reaction with the metal electrodes during simulation. As the layer insulators were surrounded by vacuum, our simulation was able to neglect the effect of substrate and electrodes. The simulation of tunneling current in the metal/2D insulators/metal diodes under a bias voltage was implemented by using DFT method in the Virtual Nanolab ATK package. [
37] Generalized Gradient Approximation (GGA) [
41] was adopted with PBE exchange correlation to describe and include the electron correlation and exchange energies in the calculation of carrier transport under bias voltage.
The design and modeling of a full-bridge AC/DC power converter built on the proposed MIM diodes was done in OrCAD PSpice simulation tool. [
44] In the circuit simulation, a PSpice device model was built for the MIM tunneling diodes by implementing the I-V characteristics and tunneling time obtained in above-mentioned DFT calculations. An AC voltage input of 1.0 V at 100 kHz was applied in the simulation and a DC output was analyzed to evaluate the performance of proposed MIM tunneling diodes. The output voltage as a function of the frequency of received signals from 1 Hz to 1000 THz was simulated based on the calculated tunneling time to analyze the high-frequency performance.