# Terahertz wave generation from hyper-Raman lines in two-level quantum systems driven by two-color lasers

###### Abstract

Based on spatial-temporal symmetry breaking mechanism, we propose a novel scheme for terahertz (THz) wave generation from hyper-Raman lines associated with the 0th harmonic (a particular even harmonic) in a two-level quantum system driven by two-color laser fields. With the help of analysis of quasi-energy, the frequency of THz wave can be tuned by changing the field amplitude of the driving laser. By optimizing the parameters of the laser fields, we are able to obtain arbitrary frequency radiation in the THz regime with appreciable strength (as strong as the typical harmonics). Our proposal can be realized in experiment in view of the recent experimental progress of even-harmonics generation by two-color laser fields.

###### pacs:

42.65.Ky, 42.50.Hz, 78.20.BhIntroduction Terahertz (THz) radiation, electromagnetic radiation with typical frequency from 0.1 THz to 10 THz-lying in the spectrum gap between the infrared and microwaves, has varieties of applications in information and communication technology, biology and medical sciences, homeland security, and global environmental monitoring etc. 1Tonouchi2007 (). The generation of THz wave is a challenge problem due to the lack of appropriate materials with small bandgaps in the usual optical approach. Great effort has been made to the design of THz sources 17Otsuji2006 (); 18Sekine2005 (); 19Orihashi2005 (); 20Crowe2005 (); 16Ito2005 (); 8Kawase2001 (); fwmix (); 10Faist1994 (); 11Kohler2002 (); 14williams2007 (); hu2004 (); 9Belkin2007 (); 21Williams2002 (); 22Bergner2005 (). THz wave can be obtained by both electronic and optical methods. The uni-travelling-carrier photodiode 16Ito2005 () and the quantum cascade laser10Faist1994 (); 11Kohler2002 (); 14williams2007 (); hu2004 (); 9Belkin2007 () are two examples. The electronic approach is limited in the low frequency end of the THz regime, and the optical approach usually focuses on the cases of small energy gap. Technological innovation in photonics and nanotechnology has provided us more new ways for THz radiation generation. Recently an interesting approach based on semiconductor nanostructure driven by acoustic wave was suggested 22a Ahn2007 (), and electrically pumped photonic-crystal THz laser was developed 22b Chassagneux2009 (). THz wave generation by up conversion method through high-order harmonics generation(HHG) in semiconductor quantum dot(QD) was also suggested in our previous work 3zhang2009PRB ().

One challenge of THz wave generation in optical approach in usual system is to generate low frequency radiation from a quantum system with a large energy gap. Usually, in the emission spectrum of a quantum system, there are harmonics, as well as the associated hyper-Raman lines. Hyper-Raman lines are caused by the transitions between the dressed bound states. Nth hyper-Raman lines are defined as those associated with the nth harmonics (they appear with the corresponding harmonics simultaneously due to the same symmetry reason, and locate near , where is the fundamental frequency of the incident laser) 23Millack (). Here is one observation that the shift of the hyper-Raman line to the associated harmonic could be small and could be tuned by the external field, even though the fundamental frequency is large. Then one may use the 0th hyper-Raman line (that associated with the 0th harmonic) to generate the low-frequency (THz) wave. There have been some studies on the harmonics and hyper-Raman lines. In many cases, there are only odd harmonics due to the particular symmetry and the associated hyper-Raman line is of high frequency and/or weak intensity. Therefore one has to solve these two problems to generate THz wave by using the hyper-Raman lines in the nonlinear optical processes.

In this article, we adopt an effective optical method to obtain THz wave in typical two-level quantum systems, which solves the two problems mentioned above. Based on the generalized spatial-temporal symmetry principle we developed recently 30sf2011NJP (), we propose a novel mechanism to obtain the 0th hyper-Raman lines in THz regime by introducing additional laser with frequency (). Using analysis of quasienergy and optimization of laser fields, we are able to obtain considerably intense THz radiation with desired frequency. With the help of our optimization method, it is quite likely that our proposal can be realized in experiment considering the recent experimental progress of even-harmonics generation by two-color laser.

Theoretical formulism Our two-level quantum system driven by laser fields [see the schematic diagram Fig. 1(a)] is described by the Hamiltonian

(1) |

where , is the Rabi frequency caused by periodic laser field , , e a unit vector. is the dipole between state and state . The energy spacing between the two states [with energies ] is set as .

The dynamics of our system is described by the equation of motion of the density matrix 30Narducci1990PRA ()

(2) |

where is the Hamiltonian, the last term describes possible dissipative effects (such as spontaneous phonon emission) and we set in the following. We numerically calculate the photon emission spectra of the two-level quantum systems by solving the density matrix in Eq.(2) through Runge-Kutta method with the time step of , total steps of 1200000 and the electron initially setting at the lower level. The average dipole can be calculated as . We use Fourier transformation to obtain the emission spectrum . In the calculation we set ( the driving field frequency) as the unit of energy. We choose as an example, yet in general one can also use other frequency (much larger than THz) laser as will be discussed later. The two-level quantum systems can be realized in many systems, such as atoms, molecules, and semiconductor quantum dots. Here we use the typical parameters of a quantum dot: the energy spacing , the dipole moment , phonon emission coefficient .

We explore the emission spectra using above numerical approach and the theory on the selection rule in high-order harmonic generation based on the generalized spatial-temporal symmetry. As described in our previous work30sf2011NJP (), for a quantum system described by Hamiltonian (1), if there exists one symmetric operation Q, which is the time shift () combined with another operation in spatial/spectrum domain, i.e. , such that the initial condition and the Hamiltonian (up to a sign) are invariant, and the dipole operator has a definite parity, then the emission spectrum contains no odd/even component if operator is even/odd. For our two-level quantum system driven by a monochromatic laser, i.e. , there exists one symmetric operation , which is the time shift combined with spatial operation : (here and in the following , refers to the annihilation operator for state ), and the Hamiltonian is invariant and the dipole operator is odd. Then odd harmonics alone are generated because of the spatial-temporal symmetry. Other types of spatial-temporal symmetry may lead to many interesting emission patterns 30sf2011NJP ().

Here we give some analysis on the hyper-Raman lines and their relation to the harmonic components. In our systems driven by a periodic field, the quasienergy state has the form , where the Floquet state can be written in the form . We consider a state . Then we have

(3) |

If the system posses a symmetry generated by the operator and the Floquet state has a definite parity under Q, i.e., , then we have . Then we have () for even/odd number and states , of same/different parity. Quite often the harmonics and the associated hyper-Raman lines appear simultaneously due to the same symmetry properties of the Floquet states. In the case with a monochromatic laser, we have odd harmonics and the associated hyper-Raman lines. While in other cases with symmetry broken, more harmonics and the associated hyper-Raman lines are generated.

Emission patterns

Let’s first look at the emission spectrum of a system driven by a monochromatic incident laser with frequency . As seen in Fig.1(b), there are odd harmonics as well as the associated hyper-Raman lines due to the spatial-temporal symmetry properties of the Floquet states as analyzed above. It is natural that, for the incident laser with frequency , the emission spectrum contains components with frequencies of , and …(odd orders of incident frequency ) as shown in Fig.1(c). One notices that there is no 0th hyper-Raman line in this case. If we use two-color lasers with frequencies and , interesting phenomena appear. As seen from Fig.1(d), even harmonics and the associated hyper-Raman lines are generated by introducing the second laser field. Especially a low frequency radiation, 0th hyper-Raman radiation is generated. We would like to point out that the emission spectrum for the case with driving field is not the supposition of those driven by and . For example, the harmonic in Fig. 1(d) neither appears in Fig.1(b) nor in Fig.1(c). It is also not the frequency summation or difference of the harmonics of systems driven by monochromatic incident laser with frequency and . In fact, the appearance of all even components (and the associated hyper-Raman lines) in Fig. 1(d) is the consequence of symmetry breaking as discussed above.

According to our theory, the symmetry of the quantum system generated by is broken by introducing the second laser with frequency (), and it is not broken by introducing the second laser with frequency (). These predictions are verified by our numerical results shown in Fig.1(e) (the second laser of frequency ) and Fig.1(f) (the second laser of frequency ). It is clear that the 0th hyper-Raman line does not appear in the cases with symmetry (see Figs.1(b) (c) (e)). Interestingly, the second laser field with frequency leads to the appearance of th harmonics (even for ) and the associated hyper-Raman line (in particular the 0th hyper-Raman line) as predicted by our theory, since the second laser with frequency breaks the spatial-temporal symmetry generated by . Using this mechanism, we can explain the experimental results of observing even harmonics in helium or plasma plumes (containing nanoparticles, carbon nanotubes, etc) driven by a two-color laser,29kim () where the second-harmonic driving term, despite being very small, breaks the symmetry and thus allows additional strong even harmonic components. Our theory also naturally explains the generation of even harmonics and associated hyper-Raman lines by breaking the spatial symmetry 6Kibis2009PRL ().

As seen from equation (3), the frequency of hyper-Raman line is determined by the quasienergy, which is tunable. One may ask whether one can use the hyper-Raman line associated with 1th harmonic to generate the low frequency/THz wave? Actually, as seen from Fig. 2(a) the intensity of the hyper-Raman line associated with the 1th harmonic decreases dramatically as the frequency decreasing. To make this point clearer, we compare the low frequency components associated the 1th harmonic and 0th harmonic for the case with two-color laser shown in Fig. 2(b). One sees that in the every low frequency regime, the intensity of the 0th hyper-Raman line (associated with the 0th harmonic) is much larger than that of the 1th hyper-Raman line. Therefore one should use the 0th hyper-Raman line to generate the THz radiation effectively.

Tuning of the frequency and intensity of THz wave

Now we discuss the optimization scheme of THz wave generation with driving field in our proposal. From equation (3), we see that the difference of the quasienergy is related to the shift of the hyper-Raman line form the corresponding harmonics, or the frequency of the 0th hyper-Raman line. This relation is verified by our numerical calculation as shown in Fig. 3. Moreover from Fig. 3, one sees that the 0th hyper-Raman line with very small frequency can be obtained in some parameter regimes, which can be used for THz wave generation. By applying this relation ( the frequency of the 0th hyper-Raman line), any designated frequency THz wave can be obtained by tuning the magnitudes of the two external fields and . For example, THz radiation of 2THz () can be obtained by tuning the external field as shown in Fig. 3 [intersection points of the solid straight-line () and the line for quasienergy]. Fig. 4(a) shows the values for and under which the low frequency is . There are disconnected parameter regimes of driving fields for generation of THz radiation. Furthermore, the intensity of the THz wave can be optimized. One can find the proper parameters for the maximal intensity with desired frequency [the triangle in Fig. 4(a)]. The corresponding emission spectrum is shown in Fig.4(b). It is seen that the intensity of THz radiation has been increased 200 times after optimization and we are able to obtain appreciably intense THz wave (as strong as typical harmonics). Here we estimate the emission power based on a system of arrays of semiconductor quantum quantum dots. The emission power of each dot is . There are about dots in the regime of size of wavelength which emit wave in phase. Therefore, the total emission power from a sample of submillimeter size is around . 6Kibis2009PRL ()

In our approach, the typical driving field intensity is in the order of , which is lower than the typical driving field intensity (in the order of used in other methods for THz generation and/or HHG by two-color laser fwmix (); 29kim (); even-two (). One should also use a long driving laser pulse [much longer than its period ()] as that used in Ref. [29kim, ]. The emitted hyper-Raman line frequency can not only be tuned by the driving field intensity as seen in Fig. 3, but also by the frequency of the incident laser. If we use the incident laser with frequency in the regime , we may obtain the 0th hyper-Raman line with frequency . In our method, there is another tunable parameter, the phase difference between two incident lasers. The main physical picture remains the same for different phases.

Unlike previous studies on non-linear driving two-level systems two-nonlinear (), the hyper-Raman lines and low frequency generation 23Millack (); 25Zhou2008 ref24 (); 24Antonino (); 23a Dakhnovskii ref21 (); HR (), our scheme of generating low frequency (THz) wave is based on the symmetry principle in typical two-level quantum systems of large energy spacing (in the order of ). It can be easily obtained from natural systems (such as atoms, molecules, and semiconductors) and artificial structures (such as quantum dots). The large energy spacing makes it robust against thermal fluctuation and also insensitive to the initial state 30sf2011NJP (). Our mechanism based on the symmetry principle is different from those (based on four-wave mixing fwmix () or transit current current ()) used before. Our non-perturbative analysis and calculation show that one only needs to tune the intensity of the two-color laser for typical two-level quantum systems, which is very convenient in experiments. We may be able to obtain tunable THz radiation due to the availability of stable, tunable driving sources above the THz regime 1Tonouchi2007 (); 29kim (). Moreover, people have already observed even harmonics in experiments with two-color laser 29kim (); even-two (). Therefore it is quite likely that the 0th hyper-Raman line (in THz regime) associated with 0th harmonic can be obtained in experiment if our optimization method is used.

Summary A novel mechanism based on spatial-temporal symmetry breaking is proposed to obtain THz wave from 0th hyper-Raman line in two-color pumped two-level quantum system. Quasienergy is calculated to determine the parameters of the incident laser to obtain radiation of arbitrary frequency from 0.1THz to 10THz. Upon optimization of the driving fields, we are able to obtain THz wave with desired frequency and appreciable intensity (as large as that of the typical harmonics).

Acknowledgments This work was partially supported by the National Science Foundation of China under Grants No.10874020, 11174042 and by the National Basic Research Program of China (973 Program) under Grants No. 2011CB922204, CAEP under Grant No. 2011B0102024, and the Project-sponsored by SRF for ROCS, SEM.

## References

- (1) M. Tonouchi, nature photonics 1, 97 (2007).
- (2) T. Otsuji, M. Hanabe, T. Nishimura, E. Sano, Opt. Express 14, 4815 (2006).
- (3) N. Sekine, K. Hirakawa, Phys. Rev. Lett. 94, 057408 (2005).
- (4) N. Orihashi, S. Suzuki, M. Asada, Appl. Phys. Lett. 87, 233501 (2005).
- (5) T. W. Crowe, W. L. Bishop, D. W. Perterfi eld, J. L. Hesler, R. M. Weikle, IEEE J. Solid-State Circuits 40, 2104 (2005).
- (6) H. Ito, F. Nakajima, T. Furuta, T. Ishibashi, Semicond. Sci. Technol. 20, S191 (2005).
- (7) K. Kawase, J. Shikata, I. Ito, J. Phys. D 34, R1 (2001).
- (8) X. Xie, J. Dai, and X.-C. Zhang, Phys. Rev. Lett. 96, 075005 (2006); T.-J. Wang, J.-F. Daigle, S. Yuan, F. Theberge, M. Chateauneuf, J. Dubois, G. Roy, H. Zeng, and S. L. Chin, Phys. Rev. A 83, 053801 (2011); T.-J. Wang, C. Marceau, Y. Chen, S. Yuan, F. Theberge, M. Chateauneuf, J. Dubois, and S. L. Chin, Appl. Phys. Lett. 96, 211113 (2010); J. Penano, P. Sprangle, B. Hafizi, D. Gordon, and P. Serafim, Phys. Rev. E 81, 026407 (2010).
- (9) J. Faist, F. Capasso, D. L. Sivco, C. Sirtori, A. L. Hutchinson, and A. Y. Cho, Science 264, 553 (1994).
- (10) R. Kohler, A. Tredicucci, F. Beltram, H. E. Beere, E. H. Linfeld, A. G. Davies, D. A. Ritchie, R. C. Iotti, and F. Rossi, Nature 417, 156 (2002).
- (11) B. S. Williams, Nat. Photonics 1, 517 (2007).
- (12) B. S. Williams, S. Kumar, Q. Hu, and J. L. Reno, Electron. Lett. 40, 431 (2004).
- (13) M. A. Belkin, F. Capasso, A. Belyanin, D. L. Sivco, A. Y. Cho, D. C. Oakley, C. J. Vineis, and G. W. Turner, Nat. Photonics 1, 288 (2007).
- (14) G. P. Williams, Rev. Sci.Instr. 73, 1461 (2002).
- (15) A. Bergner, U. Heugen, E. Brundermann, G. Schwaab, M. Havenith, D. R. Chamberlin, and E. E. Haller, Rev. Sci. Instr. 76, 063110 (2005).
- (16) K. J. Ahn, F. Milde, and A. Knorr, Phys. Rev. Lett. 98, 027401 (2007).
- (17) Y. Chassagneux, R. Colombelli, W. Maineult, S. Barbieri, H. E. Beere, D. A. Ritchie, S. P. Khanna, E. H. Linfield and A. G. Davies, Nature(London) 457, 174 (2009).
- (18) S.Q. Duan, W. Zhang, Y. Xie, W.D. Chu, and X.G. Zhao, Phys. Rev. B 80, 161304(R) (2009).
- (19) T. Millack and A. Maquet, J. Mod. Opt. 40, 2161 (1993).
- (20) S. F. Guo, S. Q. Duan, Y. Xie, W. D. Chu and W. Zhang, New J. Phys. 13, 053005 (2011).
- (21) L. M. Narducci, M. O. Scully, G.-L. Oppo, P. Ru, and J. R. Tredicce Phys. Rev. A 42, 1630 (1990).
- (22) R. A. Ganeev, H. Singhal, P. A. Naik, J. A. Chakera, H. S. Vora, R. A. Khan and P. D. Gupta, Phys. Rev. A 82, 053831 (2010); R. A. Ganeev, H. Singhal , P. A. Naik, I. A. Kulagin, P. V. Redkin, J. A. Chakera, M. Tayyab, R. A. Khan and P. D. Gupta, Phys. Rev. A 80, 033845 (2009).
- (23) O. V. Kibis, G. Ya. Slepyan, S. A. Maksimenko, and A. Hoffmann, Phys. Rev. Lett. 102, 023601 (2009).
- (24) I. J. Kim , C. M. Kim, H. T. Kim , G. H. Lee , Y. S. Lee, J. Y. Park, D. J. Cho and C. H. Nam, Phys. Rev. Lett. 94, 243901 (2005); I. J. Kim, G. H. Lee, S. B. Park, Y. S. Lee, T. K. Kim, C. H. Namb, T. Mocek and K. Jakubczak, Appl. Phys. Lett. 92, 021125 (2008); N. Ishii, A. Kosuge, T. Hayashi, T. Kanai, J. Itatani, S. Adachi, and S. Watanabe,Opt. express 16, 20876 (2008).
- (25) K. B. Nordstrom, K. Johnsen, S. J. Allen, A.-P. Jauho, B. Birnir, J. Kono, T. Noda, H. Akiyama, and H. Sakaki, Phys. Rev. Lett. 81, 457 (1998); L. Plaja and L. Roso, J. Mod. Opt. 40, 793 (1993).
- (26) A. D. Piazza and E. Fiordilino, Phys. Rev. A 64, 013802 (2001).
- (27) Z. Y. Zhou and J. M. Yuan, Phys. Rev. A 77, 063411 (2008).
- (28) Y. Dakhnovskii and H. Metiu, Phys. Rev. A 48, 2342 (1993).
- (29) C. Liu, S. Gong, R. Li, and Z. Xu, Phys. Rev. A 69, 023406 (2004); F. I. Gauthey, C. H. Keitel, P. L. Knight, and A. Maquet, Phys. Rev. A 52, 525 (1995); M. Frasca, Phys. Rev. A 60, 573 (1995); A. D. Piazza, E. Fiordilino and M. H. Mittleman, Phys. Rev. A 64, 013414 (1995); M. L. Pons, R. Taieb, and A. Maquet, Phys. Rev. A 54, 3634 (1996); H. Wang and X.-G. Zhao, J. Phys.: Condens. Matter 8, L285 (1996).
- (30) K. Y. Kim, J. H. Glownia, A. J. Taylor, and G. Rodriguez, Opt. Express 15, 4577 (2007); K. Y. Kim, A. J. Taylor, J. H. Glownia, and G. Rodriguez, Nat. Photon. 2, 605 (2008); Y. Minami, M. Nakajima, and T. Suemoto, Phys. Rev. A 83, 023828 (2011).