# Vertical cavity surface emitting terahertz laser

###### Abstract

Vertical cavity surface emitting terahertz lasers can be realized in conventional semiconductor microcavities with embedded quantum wells in the strong coupling regime. The cavity is to be pumped optically at half the frequency of the exciton state. Once a threshold population of excitons is achieved, a stimulated terahertz transition populates the lower exciton-polariton branch, and the cavity starts emitting laser light both in the optical and terahertz ranges. The lasing threshold is sensitive to the statistics of photons of the pumping light.

###### pacs:

78.67.Pt,78.66.Fd,78.45.+hCreation of efficient sources of terahertz (THz) radiation is of high importance for various fields of modern technology, including information transfer, biosensing, security and othersDragoman (). These applications are currently limited due to the lack of compact and reliable solid state sources of THz radiation. The main obstacle preventing the creation of such a source is the low rate of spontaneous emission of THz photons: according to Fermi’s golden rule this rate is proportional to the cube of the frequency and for THz transitions should be roughly tens of inverse milliseconds, while lifetimes of crystal excitations typically lie in the picosecond rangeDuc (); Doan (). The attempted strategies to improve this ratio include using the Purcell effectGerard (); Todorov () by embedding the sample inside a THz cavity, or using the cascade effect in quantum cascade lasersFaist () (QCLs). Nevertheless, until now QCLs in the spectral region around 1 THz remain costly, have macroscopic dimensions, are short-lived and still show a quantum efficiency of less than 1%.

Recently, it has been proposed that the emission rate of THz photons may be increased by bosonic stimulation if the THz transition feeds a condensate of exciton-polaritons. In Refs.Kavokin1, ; SavenkoPRL, ; delValle, the authors suggest using the transition between the upper and lower polariton branches in a semiconductor microcavity in the polariton lasing regime. The radiative trasition between two polariton modes, accompanied by emission of a THz photon, may become allowed if an electric field mixing polariton and dark exciton states is applied to the cavity. THz photons would be emitted in the plane of the cavity, and a supplementary lateral THz cavity would be needed to provide the positive feedback.

Here we propose an alternative model of a microcavity based THz laser, which has significant advantages compared to one based on the transition between two polariton branches. We propose using two-photon pumping of a exciton state, as has been realised already in GaAs based quantum well structuresCingolani (); Schemla (). The direct transition to or from the exciton state with emission or absorption of a single photon is forbidden by optical selection rules. Instead, a exciton can radiatively decay to the lower exciton-polariton mode formed by the exciton and cavity photon. This transition is accompanied by emission of a THz photonGarro (). The inverse process (THz absorption by a lower polariton mode with excitation of a exciton) has been recently observed experimentally KochCondmatt (). The THz transition from the state pumps the lowest energy exciton-polariton state, which eventually leads to the polariton lasing effect, widely discussed in the literature Christopoulos (). A macroscopic occupation of the lowest energy polariton state stimulates emission of THz photons, so that, in the polariton lasing regime, the cavity would emit one THz photon for each optical photon emitted by the polariton laser, ideally.

The design of the laser and the involved energy levels are illustrated schematically in Figure 1. This design has two crucial advantages with respect to that previously considered Kavokin1 (): it allows for operation with an optically allowed THz transition, and it provides vertical emission of THz photons. The whole structure is microscopic; no waveguides or THz cavities are needed. This latter point is also a significant advantage with respect to the quantum cascade laser, which operates in the wave-guide geometry. Moreover, the suggested scheme is very interesting from a fundamental point of view, since, as shown below, the threshold of the proposed laser appears to be sensitive to the statistics of photons of the pumping light. Here we present the quantum model of this laser based on the Liouville equation for the density matrix describing the exciton state, the lower polariton mode, and optical and THz photons.

The system under study consists of 2 excitons and 1 exciton-polaritons interacting with an external electromagnetic field and phonons in a semiconductor microcavity. Hence, our consideration involves both coherent and incoherent processes. Consequently, it is convenient to split the total Hamiltonian of the system into two parts,

(1) |

where the part includes all types of coherent processes in our system and part describes the decoherence in the system due to exciton interactions with acoustic phonons, modelled as a classical reservoir Magnusson (); Savenko ().

The equation for the density matrix can be written in the following form:

where is a Lindblad superoperator defined below. The incoherent part of the Hamiltonian can be divided into two parts, , where creates an excitation in the quantum system (and thus annihilates an excitation in the classical reservoir), and , conversely, annihilates an excitation in the quantum system and creates an excitation in the classical reservoir. The subscript index numerates different types of decoherence processes in the system.

The Lindblad terms can be then written as Savenko (); Ostatnicky ()

where accounts for energy conservation, and in realistic calculations should be taken as an average inverse broadening of the states in our system, . In the particular system which we consider, the coherent and incoherent parts of the Hamiltonian can be written as

(3) | |||||

(4) |

where

(5) | |||||

(6) | |||||

(7) | |||||

(8) |

Here and denote annihilation operators of the exciton and the lowest energy 1 polariton states, respectively, is the annihilation operator for laser photons exciting the exciton state, and is an annihilation operator for THz photons produced by the transition. , , and are the energies of the exciton, exciton-polariton, pump photon () and THz photon, , respectively. Constants and define the strengths of the two-photon absorption and radiative transitions, respectively. While the transition strength is determined by the dipole matrix elements, the main contribution to the two-photon excitation comes from processes with intermediate states in the valence or conduction bands ivchenkopikus (). We assume that the upper polariton state is higher in energy than the 2 exciton state, and thus can be excluded from consideration. The classical reservoir consists of the photons of the external laser light and THz photons. Note, that the coherent part of our Hamiltonian describes only the free modes of the system and thus is irrelevant to the kinetic equations.

Now, for the occupancy of the 2 exciton mode one can obtain after straightforward algebra:

(9) | |||||

where denotes averaging with the appropriate density matrix , and . The equations for the polariton mode occupancy , and the terahertz mode occupancy are analogous to those for . The occupancy of the pumping mode is defined by the intensity of the external pump and we do not need to write an independent dynamic equation for it. The same holds true for the higher order correlators involving pump operators, e.g. . It follows from Eq. (9) that the dynamic equations for the occupancies of the modes contain quantum correlators of fourth order, such as . For them one can also write the dynamic equations analogous to Eq. (9), which would contain correlators of sixth order. Proceeding further, one would obtain an infinite chain of coupled equations for the hierarchy of correlators.

In order to solve this chain of equations one needs to truncate the correlators at some stage. Here we use the mean-field approximation, which consists in truncation of the fourth-order correlators into products of second-order ones. One can approximate etc.

A particular analysis is needed for the truncation of the correlator containing four operators corresponding to the pumping mode. Using the definition of the second order coherence , it can be represented as

(10) |

Finally, the closed set of equations of motion describing the dynamics of our system reads

where we have introduced the non-radiative lifetime of the 2p exciton state , and lifetimes of lower polaritons and THz photons, and , respectively.

Let us analyze in more detail the kinetic equation for the state occupation . The terms describing incoming and outcoming scattering rates for this state can be written as

As we already mentioned in Eq. (10), the incoming rate is proportional to . The outgoing rate can be recast as . Hence, all terms containing cancel each other and Eq. (Vertical cavity surface emitting terahertz laser) holds. Note, that for processes of non-degenerate two-photon absorption or emission neglected here, two-photon emission processes where photon frequencies are different can be included in , the decay rate of the state. The scattering rates have standard form , independent of the photon field statistics enaki (). In what follows we assume that making it possible to neglect all other modes.

The key feature of Eq. (Vertical cavity surface emitting terahertz laser) is that the pumping term contains the second order coherence of the pump, i.e. strongly depends on its statistics. For a laser pump corresponding to the coherent state this statistics is Poissonian and . On the other hand different radiation sources may provide different statistics of pump photons and different values of . For example, for a thermal pump . Different values of lead to the different lasing thresholds in the system. Note that dependence of two-photon optical processes on the second order coherence of light has been discussed by Ivchenko in Ref. Ivchenko, . In Figure 2, we plot the steady state solutions for and , assuming that the terahertz mode occupation is zero, for the cases of coherent and thermal pumps. For a qualitative understanding, we note that the threshold is reached if , or equivalently when , since this signals the onset of Bose stimulation of the transition from the -state to the lower polariton, accompanied by terahertz emission. The terahertz generation rate is given by , and therefore shows the same threshold behaviour with pumping intensity as , as can also be seen in Figure 2. We observe that the threshold is higher in the case of a coherent pump, which may be understood in terms of the enhanced losses from the -state due to stimulated two-photon emission.

The dynamics of switching the system to polariton and terahertz lasing is essentially governed by the dynamics of the correlator , i.e. by the product of squared intensity and second order coherence of the pumping light. The efficiency of two-photon absorption and the polariton lifetime , dependent on the quality factor of the optical cavity, are two main factors which govern the threshold to lasing. The polarisation of pumping defines the polarisation of emission, and the frequency of THz radiation may be tuned by changing the exciton-photon detuning in the microcavity. Interestingly, THz lasing in our system is possible in the absence of a THz cavity. In order to demonstrate this, we have taken the terahertz photon lifetime in the calculation shown in Figure 2. One can see that THz lasing begins simultaneously with polariton lasing. The stimulated emission of THz photons is achieved due to the bosonic stimulation of this transition by occupation of the final ( polariton) state. The system emits THz radiation through the Bragg mirrors of the microcavity, which are transparent at THz frequencies.

In conclusion, we have developed a quantum theory of vertical cavity surface emitting terahertz lasers. The vertical terahertz lasing coexists with polariton lasing and, consequently, is characterised by a low pumping threshold. Interestingly, the value of the threshold is strongly dependent on the statistics of photons of the pumping light. For practical realization of compact terahertz lasers operating at room temperature, microcavities based on wide band gap semiconductors (GaN or ZnO) would seem to be the most advantageous. The resonant two photon pumping of 2 exciton states in such a cavity may be assured by a conventional vertical cavity light emitting diode (VCLED) emitting in red. This prospective structure would consist of a GaN or ZnO microcavity grown on the top of a GaAs based VCLED structure. As emission of coherent terahertz light coexists with the polariton lasing in the optical frequency range in this scheme, one can use the visible blue or green light produced by the polariton laser as a marker for the terahertz light beam, which may be important for applications in medicine, security control, etc.

AVK thanks Alberto Bramati, Elisabeth Giacobino, Elena Del Valle and Fabrice Laussy for fruitful discussions. This work has been supported by the visitors program of the International Institute of Physics (Natal, Brazil), and EU IRSES project ”POLAPHEN”. IAS acknowledges the support from Rannis “Center of Excellence in Polaritonics”. MMG was partically supported by RFBR.

## References

- (1) D. Dragoman, M. Dragoman, Progr. in Quant. Electronics 28, 1 (2004).
- (2) H. T. Duc et al, Phys. Rev. B 74, 165328 (2006).
- (3) T. D. Doan et al, Phys. Rev. B 72, 085301 (2005).
- (4) J.-M. Gerard and B. Gayral, Journ. Lightwave Technol. 17, 2089 (1999).
- (5) Y. Todorov et al, Phys. Rev. Lett. 99, 223603 (2007).
- (6) R.F. Kazarinov, and R.A. Suris, Sov. Phys. Semiconductors 5, 707 (1971); J. Faist et al, Science, 264, 553 (1994); E. Normand et al, Laser Focus World, 43, 90 (2007).
- (7) K. V. Kavokin et al, Appl. Phys. Lett. 97, 201111 (2010).
- (8) I.G. Savenko, I.A. Shelykh and M.A. Kaliteevski, Phys. Rev. Lett. 107, 027401 (2011).
- (9) E. del Valle and A. V. Kavokin, Phys. Rev. B 83, 193303 (2011).
- (10) I.M. Catalano et al, Phys. Rev B 40, 1312 (1989).
- (11) R. A. Kaindl, D. Hagele, M. A. Carnahan, and D. S. Chemla, Phys. Rev. B 79, 045320 (2009).
- (12) N. Garro, S.P. Kennedy, A.P. Heberle, R.T. Phillips, Physica Status Solidi (b), 221, 385 (2000).
- (13) J.L. Tomaino et al, arXiv:1112.2185 (2011).
- (14) see e.g. S. Christopolous et al, Phys. Rev. Lett. 98, 126405 (2007); A. Das et al, Phys. Rev. Lett. 107, 066405 (2011).
- (15) E. B. Magnusson et al, Phys. Rev. B 82, 195312 (2010).
- (16) I. G. Savenko et al, Phys. Rev. B 83, 165316 (2011).
- (17) T. Ostatnicky, I.A. Shelykh and A.V. Kavokin, Phys. Rev. B 81, 125319 (2010).
- (18) E.L. Ivchenko and G.E. Pikus, Superlattices and other heterostructures (Springer-Verlag, Berlin, 1997).
- (19) N. A. Enaki and O. B. Prepelitsa, Theoretical and Mathematical Physics 88, 967 (1991).
- (20) E.L. Ivchenko, Optical Spectroscopy of Semiconductor Nanostructures (Alpha Science, Harrow UK, 2005).