Abstract
The aim of this work is to improve the study of a phonon laser (saser) proposed by us several years ago[1]. This is a device capable to generate an intense coherent beam of acoustical phonons. Our acoustic laser consists in a double barrier heterostructure tailored such the energy difference between the ground and the first excited state in the well is close to the energy of the LO phonon. The electrons are directly injected into the excited level. Therefore they decay producing a high rate of LO phonons. These phonons are confined inside the well and decay into a pair of phonons[2]: LO -> <img src="http:/img/fbpe/bjp/v29n4/lotil.gif">+ TA. The TA phonons escape the well in the [111] direction constituting an intense coherent beam. Recently were studied (and sometimes realized experimentally) several kinds of phonon lasers. Up to our knowledge our saser is the only that has a very short wavelength (smaller than 25 Å) and a very long range (greater than 1000 mum). Because of that, such beam could have applications to acoustic nanoscopy, acoustic nanolithography and phonoelectronics. In early articles[1, 3, 4, 5, 6] we get the kinetic equations for the averaged electron and phonon populations. Quantum fluctuations were not taken into account. The system Hamiltonian is H = He + Hph + He<FONT FACE=Symbol>-</FONT>ph + Hph<FONT FACE=Symbol>-</FONT>ph + He<FONT FACE=Symbol>-</FONT>e. To solve this Hamiltonian we expand their eigenfunctions in the basis of the eigenstates¦jn1n2n3<FONT FACE=Symbol>ñ</FONT> of the single particle part of it. We obtain a set of coupled equations for the expansion coefficients that can be solved with some approximations. The results are qualitatively similar to those obtained previously.
A Quantum Formalism for a Terahertz Acoustic Laser
Ihosvany Camps Rodríguez,
Physics Faculty, University of Havana, Cuba and
Instituto de Fí sica, Universidade Federal Fluminense,
Niterói - RJ, Brazil.
Sergio Saul Makler,
Instituto de Física, Universidade Federal Fluminense,
Niteroi-RJ, Brazil
and Enrique Victoriano Anda
Instituto de Física, Universidade Federal Fluminense, Niteroi-RJ, Brazil and
Departamento de Física, Pontificia Universidade Católica do Rio de Janeiro,
Caixa Postal 38071, 22452-970, Rio de Janeiro - RJ, Brazil Received February 6, 1999 2
Received February 6, 1999
The aim of this work is to improve the study of a phonon laser (saser) proposed by us several years ago[1]. This is a device capable to generate an intense coherent beam of acoustical phonons. Our acoustic laser consists in a double barrier heterostructure tailored such the energy difference between the ground and the first excited state in the well is close to the energy of the LO phonon. The electrons are directly injected into the excited level. Therefore they decay producing a high rate of LO phonons. These phonons are confined inside the well and decay into a pair of phonons[2]: LO ®
+ TA. The TA phonons escape the well in the [111] direction constituting an intense coherent beam. Recently were studied (and sometimes realized experimentally) several kinds of phonon lasers. Up to our knowledge our saser is the only that has a very short wavelength (smaller than 25 Å) and a very long range (greater than 1000 mm). Because of that, such beam could have applications to acoustic nanoscopy, acoustic nanolithography and phonoelectronics. In early articles[1, 3, 4, 5, 6] we get the kinetic equations for the averaged electron and phonon populations. Quantum fluctuations were not taken into account. The system Hamiltonian is H = He + Hph + He-ph + Hph-ph + He-e. To solve this Hamiltonian we expand their eigenfunctions in the basis of the eigenstates ½jn1n2n3ñ of the single particle part of it. We obtain a set of coupled equations for the expansion coefficients that can be solved with some approximations. The results are qualitatively similar to those obtained previously.
I Introduction
In the recent years, the study of mesoscopic system has steadily grown. New theoretical approaches appeared at the time when new technologies were developed to produce systems of sizes ranging from 5 nm to 500 nm. In particular, double barrier heterostructures (DBH) have attracted a large attention due to its potential applications as ultra-high-frequency electronic oscillators, diodes, transistors and other electronic devices [7].
After the work of Goldman, Tsui and Cunningham [8, 9], it was realized the importance of electron-phonon interactions on the electronic properties of these structures. However, little importance has been given to the study of the phonons generated in this process, the way they propagate, their decay processes, etc.
This paper is dedicated mainly to the study of phonon generation in a DBH under the effect of an external applied bias. In previous works [3, 4, 6] this study was done through a set of phenomenological kinetic equations for the average populations. In this paper a set of quantum kinetic equations is obtained to describe the electron and phonon population using a full quantum treatment.
II The device
The device studied here consist on a DBH made of GaAs-AlGaAs. The sizes of the system are calculated in order to get De = e1-e0£ 0 for an applied bias V such that the first excited level is above the Fermi level
in the emitter. Here e0 is the energy of the ground state in the well, e1 is the energy of the first excited level and w0 is the G-point LO phonon frequency. By increasing V the levels in the well are lowered in relation to
. If
is low enough then the ground state falls below the bottom of the conduction band and the current is almost suppressed until the level e1 reaches
. By further augmenting V the current begins to flow through the excited level, but since De remains much less than
0phonon emission is inhibited. For a given bias V the resonant condition De »
0 is achieved and the electrons begin to decay to the ground state by emitting LO phonons. The potential profile and the level positions at this resonant condition are shown in figure 1.
As our system has wide barriers, the electrons remain confined in the well for long time. For an Al concentration greater than 0.25 [10] or 0.3 [11] the LO phonons are also confined inside the well. They can be also absorbed by exciting electrons from e0 to e1. This process acts in parallel with the decay of primary phonons due to anharmonicity. One of the products of this decay are the TA phonons. The device proposed here could produce a continuous beam of TA phonons in the range of 2 THz. These secondary TA phonons are coherent and form a beam that we call saser by analogy with a laser. The ultra short wavelength of these phonons permits potential applications that were discussed in a previous paper[5].
III The Hamiltonian
Our system is described by the Hamiltonian:
where
o is the single particle Hamiltonian given by an electronic part![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/fa26aab6fbb6eefccdc167109ce204d1b1d22b75.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/fa26aab6fbb6eefccdc167109ce204d1b1d22b75.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/446c39d524de2e8dbfe1d08aa197dc3bcd5fd938.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/fa26aab6fbb6eefccdc167109ce204d1b1d22b75.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/18c7baa136d39222d08eeba422acb7577f1f806f.gif)
int is the Hamiltonian that describes the electron-phonon (
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/fa26aab6fbb6eefccdc167109ce204d1b1d22b75.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/fa26aab6fbb6eefccdc167109ce204d1b1d22b75.gif)
The electronic system is described by a tight binding Hamiltonian with hoppings v between nearest neighbors. The Hamiltonian can be written in a Wannier basis, in terms of the electronic creation and annihilation operators cjlm+, cjlm at sites jlm,
with j, l, m = - , ¼,
.
Because the system has translational symmetry in the directions perpendicular to the grow direction z, the Hamiltonian can be uncoupled. We expand the operators cjlm in plane waves in the xy direction,
By doing that, we can treat the system as a sum over 1D Hamiltonians for each wave vector k^ perpendicular to the current in the z direction
where =
- 4v. The energies
are measured from the bottom of the conduction band of the emitter and
= ej +
. The energies ej are chosen to describe the energy profile of the DBH.
For the sake of simplicity we will left implicit the dependence with k^.
For the z direction we separate the space in three regions: the dispersion region and two semi-infinity unidimensional chains. To the left we will have planes with energy ej= 0 for j £ 0 and to the right ej = V for j ³ L+1, L is the length of the DBH. The corresponding eigenstates of these two regions are planes waves. If we disconnect the DBH from the left and right chains we get for the dispersion region the profile of an infinity (even not rectangular) well as it is shown in figure 2.
For the dispersion region we diagonalize a three-diagonal matrix of order L corresponding to the profile showed above getting the eigenvalues em and the eigenvectors |mñ of equation
Here,
ez is the part of![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/fa26aab6fbb6eefccdc167109ce204d1b1d22b75.gif)
After the diagonalization we get L discrete levels. We label these levels through the index m as m = 0, ¼. For the right chain we rename the planes, from L+1, L+2, ¼ to 1, 2, ¼; and for the left one, from 0, -1, -2,¼, -
to
,
, ¼,
. Therefore we get a new picture where we have two semi-infinity chains, the first one corresponding to the planes j =
, ¼,
and the other one corresponding to j = 1, ¼,
. Between the L levels obtained from the diagonalization of the matrix (8), only the two states with lower energies will participate significantly in the electronic transport. The others are far over the Fermi level, thus they do not contribute. Therefore the dispersion region is represented by j = 0 for the fundamental level and j =
for the excited level, with energies e0 and
respectively.
To connect the DBH with the left and right chains we computed the matrix elements
Only four values are relevant in our calculations: ,
,
, and
. Finally, the electronic Hamiltonian is
As it is well know [12], the dominant dispersion process in polar semiconductors is due to the coupling between electrons and LO phonons. In this work we also considered the phonons resulting from the decay of the LO phonons, i.e., the and the TA phonons. The Hamiltonians
LO,
,
TA used here are the simplest for each kind of phonons considering just one mode. They are
where
1,![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/c830c429f88a0389c5c2413ea9b53c32a6773253.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/c830c429f88a0389c5c2413ea9b53c32a6773253.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/18c7baa136d39222d08eeba422acb7577f1f806f.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/6ed860d74fd9d8cdcd448a79c95d00eac3ec8310.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/3e8e361a16ea4aa16ca81c35e255401e83db1eb7.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/3e8e361a16ea4aa16ca81c35e255401e83db1eb7.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/18c7baa136d39222d08eeba422acb7577f1f806f.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/fa26aab6fbb6eefccdc167109ce204d1b1d22b75.gif)
The electron-phonon interaction Hamiltonian has the form,
where g is a constant that measure the strength of the interaction.
The phonon-phonon interaction is given by[2]
This Hamiltonian describes the decay LO® + TA and its inverse process (recombination).
Now we define the operators:
The operators creates the eigenstates of
o:
formed by an electron at plane j and n1LO phonons, n2 phonons, and n3TA phonons.
Using the operators (16) and the Hamiltonian (1) we can calculate their equations of motion. In order to simplify the notation we call n for n1,n2,n3:
For the stationary solutions (t) =
(0)eiwt, the previous equations can be written as:
If we replace the Hamiltonian (1) into equation (19) we get the following equations:
where = ej+n1
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/c830c429f88a0389c5c2413ea9b53c32a6773253.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/c830c429f88a0389c5c2413ea9b53c32a6773253.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/c830c429f88a0389c5c2413ea9b53c32a6773253.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/c830c429f88a0389c5c2413ea9b53c32a6773253.gif)
We expand the eigenstates |Fñ of in the eigenstates of
o
Due to the orthogonality of |jnñ the amplitudes can be calculated as = á jn|Fñ = á0|
|Fñ.
IV The kinetic equations
To obtain the kinetic equations that describe the population of electrons and phonons we need to obtain the equations of motion for the electron number operators
cj and the phonon number operators![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/6ed860d74fd9d8cdcd448a79c95d00eac3ec8310.gif)
Defining the averaged populations for electrons in the excited and ground states as:
and the averaged populations for the LO, , and TA phonons as:
and using the equations (23) we get the set of kinetic quantum equations:
The population rates in (24,25,26, 27,28) are: the electronic input rate from the emitter
and the output rate from the excited level to the collector,
The net balance between the emission and absorption of LO phonons via the electron transitions is given by
The input rate from the emitter directly to the ground state,
is equal to zero near the saser resonant condition because the products () are real numbers, when the ground state merges below the bottom of the conduction band at the emitter. This is due to the fact that the modes connected to the ground state at left are evanescent modes.
The output rate from the ground state to the collector is
The decay of one phonon LO by the emission of a pair of and TA phonons, and the inverse process (i.e., recombination) are given by
The decay rate of the population turns out to be
Finally the escape rate of the TA phonons is given by
A diagram for the input and output rates is shown in figure 3.
To put the kinetic equations (24,25,26,27,28) as a function of the populations we have to replace the amplitudes from (22). By doing that, we get a new set of kinetic equations something different form the previous one[5]. For the sake of comparison we show the equations derived in reference [5],
where G is the input current, and R0 are the escape rate through the right barrier, w is the transition rate due to LO phonon emission, go is the decay rate of the LO phonons and t-1 = k2 is the escape rate of the TA phonons. Comparing (37) with (24) we see that the escape terms are coincident. By other hand, the second term appearing in (24) contain the second one in (37) plus some other cross terms that take into account that our system, due the escape of the electron from the dispersion region to the left chain, does not have levels but slightly hybridized peaks.
V The calculation of the current
In order to calculate the electronic current and the saser intensity we have to solve the system (22). In previous papers [4,6] we solved a restrict version of these equations in which the peaks corresponding to the excited level and the satellite of the ground state mixed in two polaronic branches. In that paper was not considered
ph-ph in the calculation of the electronic current. In this paper the phonon-phonon term is taken into account explicitely. In figure 4 we can see the two polaronic branches when the lifetime of the![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/18c7baa136d39222d08eeba422acb7577f1f806f.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/18c7baa136d39222d08eeba422acb7577f1f806f.gif)
![](https://minio.scielo.br/documentstore/1678-4448/FVN3WtBm3Xn8j7SddsyWhcn/18c7baa136d39222d08eeba422acb7577f1f806f.gif)
VI Conclusions
In this work, a set of kinetic equations was obtained using a full quantum treatment. This set is different from those presented in previous works where the kinetic equations were derived phenomenologically. Also, we consider for the first time all the relevant phonons and their interactions in the calculation of the current.
It can be seen that, if the decay of the phonons is very rapid, the current is the same as in previous calculations. Even if we consider a very long lifetime for these phonons, only small modifications are observed for the current.
In summary, an accurate full quantum calculation for a DBH phonon laser confirms the results obtained using rough phenomenological methods. The Hamiltonian presented here is suitable to study the synergetics of our laser.
Acknowledgments
One of us (IC) wants to acknowledge Latin American agency CLAF for financial support.
- [1] E. V. Anda, S. S. Makler, H. M. Pastawski, R. G. Barrera, Braz. J. of Phys. 24, 330 (1994).
- [2] F. Vallée, Phys. Rev. B 49, 2460 (1994).
- [3] S. S. Makler, D. E. Tuyarot, E. V. Anda and M. I. Vasilevskiy, Surface Science 361, 239 (1996).
- [4] D. E. Tuyarot, S. S. Makler, E. V. Anda and M. I. Vasilevskiy, Superlattices and Microstructures 22 (3) 427 (1997).
- [5] S. S. Makler, M. I. Vasilevskiy, E. V. Anda, D. E. Tuyarot, J. Weberszpil and H. M. Pastawski, J. of Phys. Cond. Matt. 10 (26), 5905 (1998).
- [6] J. Weberszpil, S. S. Makler, E. V. Anda and M. I. Vasilevskiy, Microelectronic Engineering 43 - 44, 471 (1998).
- [7] F. Capasso and S. Datta, Phys. Today 43, 74 (1990).
- [8] V. J. Goldman, D. C. Tsui and J. E. Cunningham, Phys. Rev. Lett. 58, 1256 (1987).
- [9] V. J. Goldman, D. C. Tsui and J. E. Cunningham, Phys. Rev. B 36, 7635 (1987).
- [10] B. Jusserand, F. Mallot, J.-M. Moison and G. Le Roux, Appl. Phys. Lett. 57, 560 (1990).
- [11] J. M. Jacob, D. M. Kim, A. Bouchalkha, J. J. Sony, J. F. Klem, H. Hou, C. W. Tu, H. Morkoç, Sol. State Comm. 91, 721 (1994).
- [12] H. Haken Quantum Field Theory of Solids, an Introduction, North-Holland, Amsterdam p. 218 (1996).
Publication Dates
-
Publication in this collection
23 Feb 2001 -
Date of issue
Dec 1999
History
-
Received
06 Feb 1999