Next Article in Journal
Chemoresistive Sensors for Cellular Type Discrimination Based on Their Exhalations
Next Article in Special Issue
Promoted Mid-Infrared Photodetection of PbSe Film by Iodine Sensitization Based on Chemical Bath Deposition
Previous Article in Journal
Observation of Spin-Glass-like Behavior over a Wide Temperature Range in Single-Domain Nickel-Substituted Cobalt Ferrite Nanoparticles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretical Analysis of Terahertz Frequency Multiplier Based on Semiconductor Superlattices

1
School of Physics and Electronic Engineering, Jiangsu University, 301 Xuefu Road, Zhenjiang 212013, China
2
Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
3
Key Laboratory of Terahertz Solid-State Technology, Shanghai Institute of Microsystem and Information Technology, Chinese Academy of Sciences, 865 Changning Road, Shanghai 200050, China
*
Authors to whom correspondence should be addressed.
Nanomaterials 2022, 12(7), 1114; https://doi.org/10.3390/nano12071114
Submission received: 1 March 2022 / Revised: 22 March 2022 / Accepted: 24 March 2022 / Published: 28 March 2022

Abstract

:
We propose a terahertz frequency multiplier based on high order harmonic generation in a GaAs-based miniband superlattice driven by an electric field. The performance of the frequency multiplier is analyzed using the balance equation approach, which incorporates momentum and energy relaxation processes at different lattice temperatures. It is found that the generated high-order harmonic power is sensitive to temperature changes. The peak power appears around resonance between driving terahertz frequency and intrinsic Bloch frequency. In the presence of the magnetic field, the peak power shifts towards a stronger static electric field region. The simulated results about the dependence of the second and third harmonic powers on a DC electric field are in qualitative consistence with the experiments. The proposed terahertz frequency multiplier based on semiconductor superlattice, being compact and efficient, is provided as a good candidate for terahertz wave generation.

1. Introduction

Semiconductor superlattice [1] is a periodical structure with layers of two or more semiconductor materials with similar lattice constants. During the past few decades, in semiconductor superlattice, many interesting phenomena have been observed, including DC current suppression with high frequency voltage bias [2,3], coherent Hall effect under the crossed electric and magnetic fields [4], periodic and chaotic electron transport under D C + A C electric field [5,6], etc. These phenomena are relevant to the nonlinearity of an electronic system and have provided abundant opportunities to develop compact sources and a spectroscopic scheme in the ‘terahertz gap’. The underlying mechanism of them is negative differential mobility (NDM), resulting from Bragg scattering of electrons at the boundary of the Brillouin zone.
Over the past twenty years, terahertz technology has attracted lots of interest due to its wide and important applications, such as nondestructive testing [7], imaging [8,9,10], high speed communication [11,12,13], and so on. Terahertz sources play a crucial role in these applications. The approaches for developing terahertz sources include optical terahertz generation, terahertz quantum cascade laser (THz QCL), and solid-state electronic devices [14]. Tunesi et al. [15] have reported that ultrafast metallization excited by a dual pump optical setup can generate a terahertz wave and, through this framework, a complex time-dependent photocarrier dynamics can be observed in black silicon. These newly developed optical terahertz generation schemes can create impact in complex nonlinear imaging and in the optical control of electron–electron interactions. THz QCL [16,17,18,19,20] operating at 2–5 THz has proved to be effective, compact, and coherent. However, the shortcoming is that THz QCL still needs to operate at liquid helium temperature to obtain excellent working performance. By comparison, the superlattice oscillator [21] can radiate microwaves at room temperature due to negative differential conductance (NDC) [22]. The radiation frequency can extend to the millimeter wave region by reducing the period of superlattice. In addition, within a certain range of the applied voltage, the current oscillations within the superlattice can be self-synchronized leading to a dramatic rise of the output power [23].
Recently, the superlattice harmonic multiplier has been proposed to generate millimeter and terahertz waves. Many theoretical [2,24] and experimental works [25,26,27,28] have been devoted to investigating the current–voltage (I–V) characteristics and the output power of the harmonic multiplier. Winnerl et al. [25] have reported frequency doubling and tripling of terahertz radiation in a voltage-biased GaAs/AlAs superlattice by taking into account elastic scattering and strong losses due to impedance mismatch among the mesa elements in their experiment. Apostolakis et al. [27] found that asymmetric current flow caused by static bias or interface roughness can generate even harmonics, and high order nonlinear output is not proportionate to the increasing input power. In this case, up to the 6th harmonic from the superlattice frequency multiplier can be realized at room temperature [28]. Moreover, strong terahertz gain in the superlattice oscillators can be achieved by avoiding space-charge instabilities if driving frequency reaches terahertz range [29,30].
Despite the rapid progress in fundamental sciences of superlattices and its application in terahertz sources, some key issues regarding the nature and performance of frequency multiplier based on miniband superlattice driven by electric and magnetic fields remain less understood. In this paper, we theoretically study the performance of the superlattice frequency multiplier by applying the balance equation approach. It is found that increasing lattice temperature can suppress the output power. The distribution of output power peak for frequency doubling and tripling is contributed from resonance between AC frequency and Bloch frequency. The paper is organized as follows: in Section 2, we give the semiclassical transport model describing electron dynamics in the superlattice frequency multiplier. The numerical results of high order harmonic coefficient and output power are presented in Section 3. The effects of temperature, external electric, and magnetic fields on the emission power of frequency multiplier are carefully analyzed. In Section 4, we summarize our work.

2. Balance Equation Approach for Electron Response of a Superlattice Frequency Multiplier

In this section, we consider the motion of electrons along the growth direction (z-axis) in the miniband of a GaAs-based superlattice. The superlattice is driven by an external terahertz radiation and a magnetic field B. The motion of electrons is confined in the lowest miniband and the inter-miniband tunneling effect is ignored. Within the tight-binding approximation, the electron energy dispersion of semiconductor superlattice can be written in the form
ε ( k ) = 2 ( k x 2 + k y 2 ) 2 m * + Δ 2 1 cos k z d ,
where d is the period of superlattice, k = ( k x , k y , k z ) is the electron wave vector, π / d < k z π / d is the longitudinal (along the superlattice axis) components of wave vector, m * is the electron effective mass of GaAs, Δ is miniband width, and is the Plank constant.
Now, we apply an external electric field parallel to the growth direction of superlattice (z axis) and a magnetic field perpendicular to it, i.e., along the x-axis. The electrons confined in the miniband are accelerated by the applied fields and scattered by disorders and phonons, leading to an overall drift motion and heating of the electron system. The macroscopic average states of the superlattice can be described by the momentum and energy balance equations, in which the average drift velocity and average electron energy are both functions of the external fields. The balance equations for electron motions in the miniband are in the form [31,32,33,34]
d υ dr d t = e E ( t ) B υ y m z * ( h e z ) υ dr τ υ ,
d υ y d t = e B υ dr m * υ y τ υ ,
d h e d t = e E ( t ) υ dr h e h e 0 τ ε ,
where υ y and υ dr are electron average velocities along the y- and z-axis, respectively. h e = h e z + 1 2 m * υ y 2 is electron average energy from semiclassical picture with h e z the average energy along the z-axis. e is the electron charge, τ υ is the momentum relaxation time, and τ ε is the energy relaxation time. m z * ( h e z ) is the energy-dependent effective mass of electrons in the miniband and expressed as
m z * ( h e z ) = M 0 1 2 h e z / Δ ,
with M 0 = 2 2 / Δ d 2 being the longitudinal effective mass at the bottom of the miniband.
Within the framework of the balance equation approach, the microscopic transport state of electron system is described by the center-of-mass momentum p d and the relative electron temperature T e . The ensemble-averaged quantities υ dr , υ y , and h e z are functions of p d and T e . The average energy is written by h e z = Δ 2 1 α 1 ( T e ) cos p d d / , where the expression of α 1 ( T e ) can be found in Ref. [34]. The thermal equilibrium energy is h e 0 = Δ 2 1 α 1 ( T ) with T the lattice temperature. In the Boltzmann’s distribution limit, we have α 1 ( T ) = I 1 ( Δ / 2 k B T ) / I 0 ( Δ / 2 k B T ) with k B the Boltzmann’s constant, I 0 and I 1 the modified Bessel functions of the zeroth and first order, respectively. Thus, the thermal equilibrium energy h e 0 will be
h e 0 = Δ 2 1 I 1 ( Δ / 2 k B T ) I 0 ( Δ / 2 k B T ) .
The applied time-dependent electric field along the growth direction of the superlattice is E ( t ) = E dc + E ac cos ( 2 π f ac t ) , where E dc is the DC electric field, E ac and f ac are the amplitude and frequency of an AC electric field, respectively. f ac is in the region of terahertz frequency. By solving Equations (2)–(4), the drift velocity υ dr of electrons and average energy h e can be obtained, and we find that they get into steady time-periodic state after a short transient process. In steady state, the drift velocity υ dr and average energy h e are both periodic functions with T ac = 1 / f ac . To analyse DC and harmonic characteristics of the superlattice, we make a Fourier transform of drift velocity υ dr , expressed by
υ dr ( t ) = υ 0 + n = 1 υ n s sin ( n ω t ) + υ n c cos ( n ω t ) ,
where υ 0 is the DC component, n = 1 , 2 , 3 is a positive integer, and υ n s and υ n c are the coefficients of sine component and cosine component for the nth harmonics, respectively. They are given by
υ 0 = 1 T ac 0 T ac υ dr ( t ) d t ,
υ n s = 2 T ac 0 T ac υ dr ( t ) sin ( n ω t ) d t ,
υ n c = 2 T ac 0 T ac υ dr ( t ) cos ( n ω t ) d t .
The current density of the miniband superlattice under terahertz electric field can be written as J ( t ) = e N υ dr , with N being the electron concentration. It is proportional to the average velocity υ dr . Therefore, the current of high-order harmonics, defined as j n ( t ) = e N υ n s sin ( n ω t ) + υ n c cos ( n ω t ) , indicates that the miniband superlattice can operate as a high frequency multiplier.

3. Performance Analysis of the Superlattice Terahertz Frequency Multiplier

In this section, with the above balance equation approach, we study the effect of electromagnetic field on the performance of the superlattice frequency multiplier. Here, a GaAs-based n-type superlattice with period d = 4.8 nm , miniband width = 70 meV , and the density of carriers N = 1.07× 10 18 cm 3 is considered. We calculate the output power of superlattice at lattice temperature 4.2 K , 77 K , and 300 K , respectively. In the calculations, the momentum relaxation time τ υ and energy relaxation time τ ε are not constant. They are functions of average electron energy h e z [34].
First, we investigate the high order harmonic generation and emitted harmonic power of semiconductor superlattice under the terahertz electric field and in the absence of magnetic field B. Figure 1 shows the DC component υ 0 , the second and third harmonic coefficients of sine component and cosine component υ 2 s , and υ 2 c , υ 3 s , υ 3 c of drift velocity in the superlattice. The lattice temperature is T = 300 K , and the frequency of the applied electric field is 0.8 THz . From the figure, we can find that second harmonic coefficients υ 2 s , υ 2 c vanish with DC bias E dc = 0 kV / cm . In other words, the superlattice frequency multiplier can only achieve odd harmonic generation with DC bias E dc = 0 kV / cm . Although the fluctuation of second harmonic coefficients is much stronger compared to the third harmonic with E dc increasing, the oscillation amplitudes of the second and third harmonics are relatively large. This implies that semiconductor superlattice can be used as a terahertz frequency multiplier with considerable power emission.
By Electromagnetic Theory, the power of the emitted nth harmonic P n ( ω ) is defined as [27]
P n ( ω ) = A μ 0 c L 2 8 n r I n ( ω ) ,
where A = 100 μ m 2 is the mesa area of a superlattice, and μ 0 and c are the permeability and the speed of light, respectively. L is effective path length through the superlattice, and n r = 13 is the refractive index of GaAs. I n ( ω ) is the root-mean-square value of the nth component of the Fourier expansion of the induced current density written as
I n ( ω ) = j n ( t ) cos ( n ω t ) t 2 + j n ( t ) sin ( n ω t ) t 2 ,
where t is performed as an integration over the period T ac .
Now, we consider the emitted power of the second and third harmonics at different lattice temperatures T = 4.2 K , 77 K , and 300 K for the same driving frequency f ac = 0.8 THz . Figure 2a–c display the second harmonic power P 2 ( ω ) generated by the superlattice frequency multiplier as functions of E dc and E ac . Figure 2d–f reveals the third harmonic power P 3 ( ω ) as functions of E dc and E ac . The lattice temperatures are as follows: Figure 2a,d 4.2 K , Figure 2b,e 77 K , and Figure 2c,f 300 K . The red and orange areas indicate higher output power, while the blue zone represents a lower one. From Figure 2a–c, we can see that, at T = 4.2 K , 77 K , and 300 K , the power maxima of the second harmonic are 0.485 mW , 0.556 mW , and 0.275 mW , respectively. With lattice temperature increasing to room temperature T = 300 K , the maximum power of the second harmonic is reduced, while the zone of output power is only slightly expanded. This tendency results from the decrease of electron drift velocity with the increase of the temperature. It is interesting that these maxima all appear at the same E dc = 6.86 kV / cm and E ac = 14.48 kV / cm regardless of the temperature. With the Bloch frequency f B = e E dc d / ( 2 π ) , we further find that these peaks locate at exactly where f B = f ac and f B = 2 f ac . The maximum power originates from resonance between driving terahertz AC frequency and intrinsic Bloch frequency. Figure 2d–f shows emitted power P 3 ( ω ) of the third harmonic at 4.2 K , 77 K , and 300 K , respectively. The corresponding maximum values of third harmonic power are 0.363 mW , 0.344 mW , and 0.169 mW . Similar to the second harmonic, the maximum power decreases as the temperature increases. However, different from the second harmonic, in the power diagram of P 3 ( ω ) , locations of the main peaks are at E dc = 0 and E ac = 27.24 kV / cm . The appearance at E dc = 0 indicates that odd harmonic coefficients play a dominant role in the frequency tripler in the absence of a static electric field. Furthermore, by comparison of Figure 2a and Figure 2b, Figure 2b and Figure 2e, or Figure 2c and Figure 2f, we can conclude that the higher the harmonic order is, the weaker the output power exhibits. Thus, it is reasonable that we will not take the higher order harmonic radiation with n > 3 into consideration in our study.
Since working at room temperature is an important feature of terahertz devices, we now focus on the harmonic generation of semiconductor superlattice at T = 300 K . Further insight on how driving frequency can affect the location of harmonic peak power is given by the color contour plot in Figure 3. Here, the driving frequencies are set as f ac = 0.35 THz , 0.55 THz , 0.66 THz , and 0.8 THz . Figure 3 shows that, by increasing the parameter f ac , the emitted maximum power of P 2 ( ω ) and P 3 ( ω ) increase slowly. The power distribution of P 2 ( ω ) is shown in Figure 3a–d. The maximum of emitted power also occurs at the resonance between AC field f ac and Bloch frequency f B . For example, Figure 3a displays that peaks of the second harmonic power appear at E dc 3.01 kV / cm , 6.03 kV / cm , 9.04 kV / cm , and 12.05 kV / cm , etc. With the increase of DC bias E dc or AC field E ac , the magnitude of the peaks drops progressively. In particular, the distribution of the peaks shows a linear function dependence on E dc and E ac . Figure 3e–h depict the power distribution of P 3 ( ω ) , whose output power peak distribution is more complex and irregular than that of P 2 ( ω ) . The maximum power peak of P 3 ( ω ) appears at E dc = 0 . As driving frequency increases, we observe that the power maximum of P 3 ( ω ) appears at a stronger AC electric field, ranging from E ac = 11.82 kV / cm to E ac = 26.70 kV / cm in Figure 3e–h. The other peak distributions of P 3 ( ω ) are not directly related to f B = i f ac ( i = 1 , 2 , 3 ) , which is different from that of the second harmonic. Figure 3e,f manifest that the magnitude of the secondary peak appearing at a small value of ( E dc , E ac ) is half of the maximum peak. When driving frequency becomes large, as shown in Figure 3h for f ac = 0.8 THz , the secondary peak splits into two with peak powers appearing at E ac = 14.28 kV / cm and E ac = 17.25 kV / cm , respectively.
To focus on the characteristic of P n ( ω ) versus static electric field, we show the emitted power profile of high order harmonic in Figure 4. Three different AC electric fields E ac = 5 kV / cm , 6 kV / cm , and 7 kV / cm are applied at room temperature with the same frequency f ac = 0.8 THz . Despite the AC electric fields, the peak power appears at the same value of E dc , where driving AC frequency is in resonance with Bloch frequency. The experimental demonstration of frequency doubling and tripling with semiconductor superlattice frequency multiplier has been reported in Ref. [25]. The superlattice consisting of 100 periods is grown by molecular beam epitaxy on a GaAs substrate. A corner cube antenna system can be fabricated to guide the terahertz radiation from a free-electron laser into the mesa element of frequency multiplier. The generated second and third harmonics are coupled out of the mesa by the corner cube antenna. The doping density of superlattice layers is N = 5 × 10 17 cm 3 . The experimental results of the second and third harmonic powers are reproduced in the insets of Figure 4. Our calculations of P 2 ( ω ) and P 3 ( ω ) are in qualitative consistence with the experiments [25]. It is noted that the experimental power is smaller by several orders and this may be attributed to light doping density and device loss in measurement of the experiment.
Finally, we demonstrate that the high order harmonic output power can be tuned by introducing a magnetic field perpendicular to the growth axis. Figure 5a,b exhibit P 2 ( ω ) and P 3 ( ω ) , respectively, as functions of E dc . The applied AC electric field is set to be constant as E ac = 3 kV / cm while the magnetic fields are different as B = 0 , 0.8 T , 1 T , 1.5 T , 2 T , and 2.3 T . The superlattice miniband width is Δ = 70 meV and lattice temperature is still T = 300 K . The AC frequency is f ac = 0.8 THz . Illustrated from Figure 5, the harmonic output power peak shifts to a stronger DC electric field with the magnetic field increasing. It can be understood by the fact that electrons perform cyclotron oscillation with frequency f c = e B m * , which increases with magnetic field and corresponds to a stronger DC electric field. Interestingly, the harmonic power profile reveals that the magnetic field has no effect on the lineshape distribution, while magnitude of harmonic power is sensitive to the magnetic field. Thus, two-dimensional images rather than contour plots are used here. It is worth noting that, as shown in Figure 5b, output power enhances along with the increase of the magnetic field at E dc = 0 kV / cm not being universal, and this rule is broken at other E ac values. On the whole, the magnetic field has an inhibitory effect on emitted harmonic power.

4. Conclusions

In summary, using the balance equation approach with relaxation time depending on lattice temperature and miniband width, we theoretically analyze harmonic coefficients of drift velocity and emitted harmonic power of a semiconductor superlattice frequency multiplier in terahertz frequency range. It is found that high lattice temperature can suppress the peak of P n ( ω ) and result in enhancement of the peak width. Meanwhile, the peak of the second harmonic output power P 2 ( ω ) appears at the resonance between driving frequency and Bloch frequency. The peak distribution of P 3 ( ω ) is much more complicated except that the emission power maximum is dominant at E dc = 0 kV / cm . With the presence of the magnetic field, the peak position of emitted power can be tuned due to cyclotron oscillation. Our simulations are in good agreement with the experiments, which verifies the validity of the balance equation approach. The compact and efficient terahertz frequency multiplier based on a semiconductor superlattice is proposed to be a good candidate for terahertz wave generation.

Author Contributions

Methodology, C.W. and J.C.; Investigation, W.F., S.W., Y.Z. and C.W.; Formal analysis, W.F., S.W. and C.W.; Writing, W.F., S.W., Y.Z., C.W. and J.C.; Funding acquisition, W.F., C.W. and J.C. All authors have read and agreed to the submitted version of the manuscript.

Funding

The APC was funded by the Science and Technology Commission of Shanghai Municipal (Grant No. 21DZ1101102).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

This work was supported by the National Key R & D Program of China (Grant No. 2017YFA0701005), the National Natural Science Foundation of China (Grant Nos. 61927813, 61975225, 62005304, 62035014), and the Science and Technology Commission of Shanghai Municipal (Grant Nos. 21DZ1101102, 21ZR1474600).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Esaki, L.; Tsu, R. Superlattice and negative differential conductivity in semiconductors. J. Res. Dev. 1970, 14, 61–65. [Google Scholar] [CrossRef]
  2. Lei, X.L. Current suppression and harmonic generation by intense terahertz fields in semiconductor superlattices. J. Appl. Phys. 1997, 82, 718–721. [Google Scholar] [CrossRef]
  3. Schomburg, E.; Ignatov, A.A.; Grenzer, J.; Renk, K.F.; Pavel’ev, D.G.; Koschurinov, Y.; Melzer, B.J.; Ivanov, S.; Schaposchnikov, S.; Kop’ev, P.S. Suppression of current through an Esaki–Tsu GaAs/AlAs superlattice by millimeter wave irradiation. Appl. Phys. Lett. 1996, 68, 1096–1098. [Google Scholar] [CrossRef]
  4. Hummel, A.B.; Blöser, C.; Bauer, T.; Roskos, H.G.; Kosevich, Y.A.; Kxoxhler, K. Electro-optic investigation of the Coherent Hall Effect in semiconductor superlattices. Phys. Status Solidi b 2005, 242, 1175–1178. [Google Scholar] [CrossRef]
  5. Soskin, S.M.; Khovanov, I.A.; McClintock, P.V.E. Regular rather than chaotic origin of the resonant transport in superlattices. Phys. Rev. Lett. 2015, 114, 166802. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Ying, L.; Huang, D.; Lai, Y.-C. Multistability, chaos, and random signal generation in semiconductor superlattices. Phys. Rev. E 2016, 93, 062204. [Google Scholar] [CrossRef] [Green Version]
  7. Shur, M. Recent developments in terahertz sensing technology. Proc. SPIE 2016, 9836, 98362Q. [Google Scholar] [CrossRef]
  8. Mekkalaki, A.; Kulkarni, M. Design of 16*16 high speed MLM using KOM for image filter application. In Proceedings of the 2016 International Conference on Electrical, Electronics, and Optimization Techniques (ICEEOT), Chennai, India, 3–5 March 2016; p. 3504. [Google Scholar]
  9. Sung, S.; Brown, E.R.; Grundfest, W.S.; Taylor, Z.D.; Selvin, S.; Bajwa, N.; Chantra, S.; Nowroozi, B.; Garritano, J.; Goell, J.; et al. THz Imaging System for in vivo Human Cornea. IEEE Trans. Terahertz Sci. Technol. 2018, 8, 27. [Google Scholar] [CrossRef]
  10. Qiu, F.C.; Tan, Z.Y.; Fu, Z.L.; Wan, W.J.; Li, M.Q.; Wang, C.; Cao, J.C. Reflective scanning imaging based on a fast terahertz photodetector. Opt. Commun. 2018, 427, 170–174. [Google Scholar] [CrossRef]
  11. Tessmann, A.; Schlechtweg, M.; Bruch, D.; Lewark, U.J.; Leuther, A.; Massler, H.; Wagner, S.; Seelmann-Eggebert, M.; Hurm, V.; Aidam, R.; et al. Terahertz monolithic integrated circuits based on metamorphic HEMT technology for sensors and communication. In Proceedings of the 2013 Asia-Pacific Microwave Conference Proceedings (APMC), Seoul, Korea, 5–8 November 2013. [Google Scholar]
  12. Wang, C.; Lu, B.; Miao, L.; Deng, X. 0.34 THz T/R front-end for wireless communication. High Power Laser Part. Beams 2013, 25, 1530–1534. [Google Scholar] [CrossRef]
  13. Zeng, H.; Liang, H.; Zhang, Y.; Wang, L.; Liang, S.; Gong, S.; Li, Z.; Yang, Z.; Zhang, X.; Lan, F.; et al. High-precision digital terahertz phase manipulation within a multichannel field perturbation coding chip. Nat. Photonics 2021, 15, 751–757. [Google Scholar] [CrossRef]
  14. Tonouchi, M. Cutting-edge terahertz technology. Nat. Photonics 2007, 1, 97–105. [Google Scholar] [CrossRef]
  15. Tunesi, J.; Peters, L.; Gongora, J.S.T.; Olivieri, L.; Fratalocchi, A.; Pasquazi, A.; Peccianti, M. Terahertz emission mediated by ultrafast time-varying metasurfaces. Phys. Rev. Res. 2021, 3, L042006. [Google Scholar] [CrossRef]
  16. Razavipour, S.G.; Dupont, E.; Fathololoumi, S.; Chan, C.W.I.; Lindskog, M.; Wasilewski, Z.R.; Aers, G.; Laframboise, S.R.; Wacker, A.; Hu, Q.; et al. An indirectly pumped terahertz quantum cascade laser with low injection coupling strength operating above 150 K. J. Appl. Phys. 2013, 113, 203107. [Google Scholar] [CrossRef]
  17. Agnew, G.; Grier, A.; Taimre, T.; Lim, Y.L.; Nikolić, M.; Valavanis, A.; Cooper, J.; Dean, P.; Khanna, S.P.; Lachab, M.; et al. Efficient prediction of terahertz quantum cascade laser dynamics from steady-state simulations. Appl. Phys. Lett. 2015, 16, 161105. [Google Scholar] [CrossRef] [Green Version]
  18. Wan, W.J.; Li, H.; Cao, J.C. Homogeneous spectral broadening of pulsed terahertz quantum cascade lasers by radio frequency modulation. Opt. Express 2018, 26, 980–989. [Google Scholar] [CrossRef]
  19. Qi, X.; Bertling, K.; Taimre, T.; Agnew, G.; Lim, Y.L.; Gillespie, T.; Robinson, A.; Brünig, M.; Demić, A.; Dean, P.; et al. Observation of optical feedback dynamics in single-mode terahertz quantum cascade lasers: Transient instabilities. Phys. Rev. A 2021, 103, 033504. [Google Scholar] [CrossRef]
  20. Khalatpour, A.; Paulsen, A.K.; Deimert, C.; Wasilewski, Z.R.; Hu, Q. High-power portable terahertz laser systems. Nat. Photonics 2021, 15, 16–20. [Google Scholar] [CrossRef]
  21. Grenzer, J.; Ignatov, A.A.; Schomburg, E.; Renk, K.F.; Pavel’ev, D.G.; Koschurinov, Y.; Melzer, B.; Ivanov, S.; Schaposchnikov, S.; Kop’ev, P.S. Microwave oscillator based on Bloch oscillations of electrons in a semiconductor superlattice. Ann. Phys. 1995, 507, 184–190. [Google Scholar] [CrossRef]
  22. Telenkov, M.P.; Mityagin, Y.A. Sequential resonant tunneling between landau levels in GaAs/AlGaAs superlattices in strong tilted magnetic and electric fields. Int. J. Mod. Phys. B 2007, 21, 1594–1599. [Google Scholar] [CrossRef]
  23. Gaifullin, M.B.; Alexeeva, N.V.; Hramov, A.E.; Makarov, V.V.; Maksimenko, V.A.; Koronovskii, A.A.; Greenaway, M.T.; Fromhold, T.M.; Patanè, A.; Mellor, C.J.; et al. Microwave generation in synchronized semiconductor superlattices. Phys. Rev. A 2017, 7, 044024. [Google Scholar] [CrossRef] [Green Version]
  24. Xia, J.B. Theory of terahertz-photocurrent resonances in miniband superlattices. Phys. Rev. B 1998, 58, 3565. [Google Scholar] [CrossRef] [Green Version]
  25. Winnerl, S.; Schomburg, E.; Brandl, S.; Kus, O.; Renk, K.F.; Wanke, M.C.; Allen, S.J.; Ignatov, A.A.; Ustinov, V.; Zhukov, A.; et al. Frequency doubling and tripling of terahertz radiation in a GaAs/AlAs superlattice due to frequency modulation of Bloch oscillations. Appl. Phys. Lett. 2000, 77, 1259–1261. [Google Scholar] [CrossRef]
  26. Pereira, M.F.; Zubelli, J.P.; Winge, D.; Wacker, A.; Rodrigues, A.S.; Anfertev, V.; Vaks, V. Theory and measurements of harmonic generation in semiconductor superlattices with applications in the 100 GHz to 1 THz range. Phys. Rev. B 2017, 96, 045306. [Google Scholar] [CrossRef] [Green Version]
  27. Apostolakis, A.; Pereira, M.F. Controlling the harmonic conversion efficiency in semiconductor superlattices by interface roughness design. AIP Adv. 2019, 9, 015022. [Google Scholar] [CrossRef] [Green Version]
  28. Pereira, M.F.; Anfertev, V.; Shevchenko, Y.; Vaks, V. Giant controllable gigahertz to terahertz nonlinearities in superlattices. Sci. Rep. 2020, 10, 15950. [Google Scholar] [CrossRef]
  29. Hyart, T.; Alekseev, K.N.; Leppänen, A.; Thuneberg, E.V. Terahertz Bloch oscillator with suppressed electric domains: Effect of elastic scattering. Int. J. Mod. Phys. B 2009, 23, 4459–4472. [Google Scholar] [CrossRef] [Green Version]
  30. Kroemer, H. Wave Packet Dynamics in a Biased Finite-Length Superlattice. arXiv 2003, arXiv:cond-mat/0310019. [Google Scholar]
  31. Ignatov, A.A.; Renk, K.F.; Dodin, E.P. Esaki-Tsu superlattice oscillator: Josephson-like dynamics of carriers. Phys. Rev. Lett. 1996, 70, 1996. [Google Scholar] [CrossRef]
  32. Lei, X.L.; Horing, N.J.M.; Cui, H.L. Microwave conductivity of a confined superlattice. Superlattices Microstruct. 1993, 14, 243. [Google Scholar] [CrossRef]
  33. Kosevich, Y.A. Anomalous Hall velocity, transient weak supercurrent, and coherent Meissner effect in semiconductor superlattices. Phys. Rev. B 2001, 63, 205313. [Google Scholar] [CrossRef]
  34. Cao, J.C.; Liu, H.C.; Lei, X.L. Chaotic dynamics in quantum-dot miniband superlattices. Phys. Rev. B 2000, 61, 5546. [Google Scholar] [CrossRef]
Figure 1. DC component υ 0 and harmonic coefficients υ 2 s , υ 2 c , υ 3 s , and υ 3 c of drifty velocity υ dr are plotted as a function of E dc in GaAs-based superlattice driven by E ac = 10 kV / cm at T = 300 K . Other superlattice parameters are d = 4.8 nm and Δ = 70 meV .
Figure 1. DC component υ 0 and harmonic coefficients υ 2 s , υ 2 c , υ 3 s , and υ 3 c of drifty velocity υ dr are plotted as a function of E dc in GaAs-based superlattice driven by E ac = 10 kV / cm at T = 300 K . Other superlattice parameters are d = 4.8 nm and Δ = 70 meV .
Nanomaterials 12 01114 g001
Figure 2. The calculated second harmonic output powers (ac) and third harmonic ones (df) as functions of E dc and E ac . The harmonic power is generated by a GaAs-based superlattice with d = 4.8 nm driven by an AC electric field with frequency f ac = 0.8 THz . The lattice temperatures are as follows: (a,d) 4.2 K , (b,e) 77 K , and (c,f) 300 K .
Figure 2. The calculated second harmonic output powers (ac) and third harmonic ones (df) as functions of E dc and E ac . The harmonic power is generated by a GaAs-based superlattice with d = 4.8 nm driven by an AC electric field with frequency f ac = 0.8 THz . The lattice temperatures are as follows: (a,d) 4.2 K , (b,e) 77 K , and (c,f) 300 K .
Nanomaterials 12 01114 g002
Figure 3. The second and third harmonic output powers as functions of E dc and E ac . The period of superlattice is d = 4.8 nm and temperature is T = 300 K . The AC electric field frequencies f ac presented in figure are defined as follows: (a,e) 0.35 THz , (b,f) 0.55 THz , (c,g) 0.66 THz , and (d,h) 0.8 THz .
Figure 3. The second and third harmonic output powers as functions of E dc and E ac . The period of superlattice is d = 4.8 nm and temperature is T = 300 K . The AC electric field frequencies f ac presented in figure are defined as follows: (a,e) 0.35 THz , (b,f) 0.55 THz , (c,g) 0.66 THz , and (d,h) 0.8 THz .
Nanomaterials 12 01114 g003
Figure 4. (a) Second and (b) third harmonic output powers are plotted as functions of the DC electric field E dc . The applied AC electric field is E ac = 5 kV / cm (black line), 6 kV / cm (red line), and 7 kV / cm (blue line), respectively. The lattice temperature is T = 300 K and the AC frequency is f ac = 0.8 THz . The insets of (a,b) describe the experimental results of the second and third harmonic power, respectively, versus normalized voltage for different excitation power that is proportional to E ac 2 . Reprinted with permission from Ref. [25]. Copyright 2000 American Institute of Physics.
Figure 4. (a) Second and (b) third harmonic output powers are plotted as functions of the DC electric field E dc . The applied AC electric field is E ac = 5 kV / cm (black line), 6 kV / cm (red line), and 7 kV / cm (blue line), respectively. The lattice temperature is T = 300 K and the AC frequency is f ac = 0.8 THz . The insets of (a,b) describe the experimental results of the second and third harmonic power, respectively, versus normalized voltage for different excitation power that is proportional to E ac 2 . Reprinted with permission from Ref. [25]. Copyright 2000 American Institute of Physics.
Nanomaterials 12 01114 g004
Figure 5. (a) Second and (b) third harmonic output powers are plotted as functions of the DC electric field E dc . The applied perpendicular magnetic fields are B = 0 T (black line), 0.8 T (red line), 1 T (blue line), 1.5 T (green line), 2 T (purple line), and 2.3 T (yellow line), respectively. The lattice temperature is T = 300 K and the AC frequency is f ac = 0.8 THz .
Figure 5. (a) Second and (b) third harmonic output powers are plotted as functions of the DC electric field E dc . The applied perpendicular magnetic fields are B = 0 T (black line), 0.8 T (red line), 1 T (blue line), 1.5 T (green line), 2 T (purple line), and 2.3 T (yellow line), respectively. The lattice temperature is T = 300 K and the AC frequency is f ac = 0.8 THz .
Nanomaterials 12 01114 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Feng, W.; Wei, S.; Zheng, Y.; Wang, C.; Cao, J. Theoretical Analysis of Terahertz Frequency Multiplier Based on Semiconductor Superlattices. Nanomaterials 2022, 12, 1114. https://doi.org/10.3390/nano12071114

AMA Style

Feng W, Wei S, Zheng Y, Wang C, Cao J. Theoretical Analysis of Terahertz Frequency Multiplier Based on Semiconductor Superlattices. Nanomaterials. 2022; 12(7):1114. https://doi.org/10.3390/nano12071114

Chicago/Turabian Style

Feng, Wei, Shuting Wei, Yonghui Zheng, Chang Wang, and Juncheng Cao. 2022. "Theoretical Analysis of Terahertz Frequency Multiplier Based on Semiconductor Superlattices" Nanomaterials 12, no. 7: 1114. https://doi.org/10.3390/nano12071114

APA Style

Feng, W., Wei, S., Zheng, Y., Wang, C., & Cao, J. (2022). Theoretical Analysis of Terahertz Frequency Multiplier Based on Semiconductor Superlattices. Nanomaterials, 12(7), 1114. https://doi.org/10.3390/nano12071114

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop