the simplest way to describe kinetic phenomena is by the boltzmann equation ( be ) .
there are however cases when the be can not be used , and one has to use quantum kinetic equations ( qkes ) that are in the simplest case equations for the one particle green functions . if the keldysh technique @xcite is used , the one particle green functions that participate in the description of kinetic processes are the retarded green function ( and its complex conjugate ) , and the `` statistical '' green function , which is simply the field correlator . for free particles these green functions are diagonal in the basis of the free particle states , and are described by the free particle energies @xmath0 , and their occupation numbers @xmath1 .
these are the quantities that appear in the be for the _ interacting _ particles _ in external fields_.    the description of kinetic phenomena by qkes instead of bes results in two effects .
the first effect is state mixing - due to interactions and external fields the green functions are nondiagonal in the basis of the free particle states .
this effect is responsible for coherence phenomena , of which a well known example is interband polarization described by semiconductor bloch equations @xcite .
the second effect is level shift and broadening - the green functions are nonzero off the energy shell within some width defined by the interactions and the external fields .
level broadening effects in qkes for phonons have been studied in refs .
@xcite .
the effects of level broadening are of importance in two cases .
the first case is that in which level broadening allows scattering processes that are forbidden by energy and momentum conservation .
an example of such a case is that of electron - electron scattering in a 2deg , in the presence of a quantizing magnetic field @xcite . if one tries to calculate the auger transition rates responsible for anti - stokes luminescence , that was observed in ref .
@xcite , using fermi s golden rule , one encounters delta functions with an argument equal to zero , due to singularities in the density of states , see ref .
@xcite .
the second case in which level broadening effects should be considered is that of narrow energy distributions : when there are small energy scales in the problem that are narrower than the quantum width of the state .
a simple example of such a case , that has one small energy scale , is that of an electron distribution , excited by a narrow band - width laser . the problem becomes more complex when another small energy scale is added , for example : the electrons are excited to the close vicinity of some threshold .
a problem of this nature will be dealt with in this paper .    in the time domain ,
level shift and broadening lead to non - markovian effects in the scattering processes @xcite .
these effects are important in the case of short excitation pulses @xcite .
if state mixing induced by coherent excitation can be neglected , and only shift and broadening effects are considered , the qkes can be presented in a form similar to the bes
. we will give a general recipe how to write such equations , and will use this recipe when considering the following problem : hot luminescence from a quantum wire ( qwr ) due to electrons that are excited by a narrow band , noncoherent laser , to the close vicinity of the threshold for optical lo phonon emission .
optical and transport phenomena in qwrs have been of great interest lately .
the great interest in transport phenomena @xcite is due to the possibility that a 1deg in a qwr will have the properties of a luttinger liquid .
the fermi edge singularity @xcite and the possibility of stimulated emission @xcite enlarged interest in optics of qwrs .
the important role played by lo phonons in the trapping of electrons @xcite from three - dimensional extended states , into one - dimensional localized states , has been reported in ref .
@xcite .
theoretical studies of qwrs were concentrated around elementary excitations of a 1deg @xcite , fermi edge singularity @xcite , excitonic effects @xcite , electron - phonon scattering rates @xcite and interactions @xcite , relaxation of photoexcited carriers @xcite , calculation of envelope states @xcite , and conductance of qwrs with self - consistent broadening effects @xcite .
the retarded green function and the `` statistical '' green function which appear in the kinetic equations are defined as @xmath2_{\pm }    \rangle,\ ] ] @xmath3_{\mp }        \rangle,\ ] ] where @xmath4 are the field operators in the heisenberg representation .
the square brackets with a plus ( minus ) sign signify an anticommutator ( commutator ) .
the upper signs refer to fermions and the lower signs refer to bosons ( this convention is used throughout the paper ) .
the spin degrees of freedom @xmath5 , will be suppressed from now on , and one can think of them as included in the space coordinates .    in case of interacting particles the green functions obey the dyson equation @xmath6 where @xmath7 is a matrix of green functions , and @xmath8 is a matrix of self energy functions @xcite : @xmath9 the superscript @xmath10 indicates free particle green functions .
the advanced green function @xmath11 .    in the case of free particles we can use the free particle eigenstates @xmath12 , and the corresponding annihilation operators @xmath13 , to write @xmath14 . in this case @xmath15 depends only on the time difference @xmath16 , and can be fourier transformed in time into @xmath17 ( the index @xmath18 stands for r , a or s ) . for an orthonormal set of states @xmath19 we can define @xmath20 such that @xmath21    from the definitions ( [ gr - def ] ) , ( [ gs - def ] ) we obtain ( near the singularities ) @xmath22,\ ] ] and @xmath23,\ ] ] where @xmath0 is the energy of state @xmath19 , and @xmath24 is the principal part function .
the upper ( lower ) term in the column corresponds to fermions ( bosons ) .
the occupation number of state @xmath19 , @xmath25 , is not necessarily the equilibrium occupation number .    in the case of interacting particles in external fields , if one neglects coherence effects , the green functions are still diagonal in @xmath19 , but the energy levels are broadened and shifted .
let us first look at the time invariant case ( systems under d.c .
conditions ) .
the level shift @xmath26 , and the level broadening @xmath27 , are defined as @xmath28       { \rm i m } \sigma^{\rm r}_{\alpha}(\epsilon),\ ] ] where @xmath29 is the retarded self energy of particles in state @xmath19 .
the column has the same meaning as in expression ( [ gr0-ae ] ) .
the retarded green function is an outcome of the dyson equation ( [ dyson ] ) , that is written as @xmath30 in the representation of the @xmath19 states . using the defintions ( [ s - w - def ] )
we obtain from equation ( [ dyson - gr ] ) ( for positive frequencies ) @xmath31 where @xmath32 is the `` smeared '' and `` shifted '' principal value , and @xmath33 is the `` smeared '' and `` shifted '' delta function .    comparing expression ( [ gr - ae ] ) with expression ( [ gr0-ae ] ) we see that @xmath34 is in fact a generalization of @xmath35 , where the level shift and level width have been introduced into the principal value and delta functions .
we therefore write @xmath36 as a generalization of @xmath37 , thus defining the occupation functions @xmath38 @xmath39    expressions ( [ gr - ae ] ) and ( [ gs - ae ] ) were written for the case of positive frequencies , which will be assumed from now on . in case of negative boson frequencies one has to multiply the retarded green function and @xmath40 in these expressions by minus one .
the definitions above agree with the relation between the statistical green function and the retarded green function at thermal equilibrium @xmath41 where @xmath42 is the equilibrium distribution function ( fermi or bose - einstein ) .
the smeared delta function obeys the following normalization @xmath43 and the occupation numbers are obtained from the occupation functions in the following manner @xmath44    relations ( [ equil ] ) - ( [ on - of ] ) can be obtained from the lehmann representation of the green functions .
these relations are valid only if the one particle green functions can be considered as diagonal in @xmath19 , in other words , in the absence of coherence effects .    as a result of definitions ( [ s - w - def ] ) and ( [ gs - ae ] ) ,
the three unknown functions @xmath45 and @xmath46 have been replaced by a different set of three unknown functions : @xmath47 and @xmath38 .
if the system is not under d.c .  conditions , one can define a new pair of time variables : the relative time @xmath48 , and the `` center of mass '' time @xmath49 .
only slow processes will be considered so that all quantities vary in @xmath50 with some time scale that is large compared to the inverse characteristic energy of the particles . in such a case one can still use the formalism given above for the green functions , but with @xmath50 as a parameter of the problem , in addition to @xmath51 ( see for example ref .
@xcite ) .    for slow processes @xmath52
again @xmath53 can be written as a generalization of @xmath35 , in the same manner as before , but now @xmath50 is a parameter which appears in all functions .
thus , expressions ( [ gr - ae ] ) and ( [ gs - ae ] ) for @xmath54 and @xmath55 , can be used with the three unknown functions @xmath56 and @xmath57 that depend on @xmath50 as well as on @xmath51 .
the known procedure for obtaining a kinetic equation for @xmath57 is to apply the operator @xmath58 to the dyson equation for @xmath59 @xcite , where the operator @xmath60 is equal to @xmath61 for electrons , and @xmath62 for photons and phonons ( we use the convention @xmath63 throughout the paper ) .
the slow variation in @xmath50 should be considered and then the equation should be fourier transformed in the time difference coordinates .
this procedure leads to the qke @xmath64 from now on the `` bar '' over @xmath65 will be dropped , with the understanding that all time variables are in fact center of mass time variables .
equation ( [ qke ] ) for photons and phonons includes a factor of @xmath66 that multiplies the time derivative on the left hand side .
since we are interested in off - shell energies that are close to the on - shell energy , this factor can be taken to be one .    in order to obtain equations for @xmath67 and @xmath68 one
has to substitute an expression for the retarded self energy function , in terms of the one particle green functions , into the definitions of @xmath27 and @xmath26 ( [ s - w - def ] ) .
these equations are in fact the imaginary and the real parts of equation ( [ dyson - grte ] ) .
this yields coupled , self - consistent equations for @xmath27 and @xmath26 , that are also @xmath69 dependent .
it is important to stress that the occupation function depends on time explicitly , while the dependence of @xmath27 and @xmath26 on time is only through their dependence on @xmath69 .
therefore in order to find @xmath69 one has to write a kinetic equation , that will include time evolution through time derivatives . the equations for the level shift and width will include their dependence on time only through the appearance of @xmath69 , and will not include time derivatives .
if one uses the self - consistent - born approximation for the self energy functions , one finds that the equations obtained for @xmath70 , and @xmath68 could have been obtained easily by applying a few generalizations to the be for @xmath1 . in the transfer from the be to the qke the number of unknowns increases from one ( @xmath69 ) , to three ( @xmath71 and @xmath26 ) , therefore apart from the kinetic equation we require two more equations for the unknown functions @xmath67 and @xmath68 , that appear in the equation for @xmath72 .
we present here a recipe which allows one to go from the be for @xmath73 , to the qke for @xmath72 , and from the decay term in the boltzmann equation ( which is defined below ) to an equation for @xmath67 .
an equation for @xmath68 can be obtained from the equation for @xmath67 using kramers - kronig relations .
the be , in its most general form , can be written as @xmath74     2 \pi    \delta ( \epsilon_{\alpha } - \epsilon_{\beta } + \epsilon_{\gamma} ... )\equiv s.\ ] ] this equation describes the evolution in time of the occupation number of particles in state @xmath19 , due to interactions with particles in states @xmath75 , with matrix elements @xmath76 .
different states can also mean different particles .
the right hand side of the equation includes delta functions that are responsible for energy conservation .
there are many cases in which it is convenient to think of the collision integral @xmath77 as composed of `` scattering in '' events and `` scattering out '' events : @xmath78 .
the first term on the right hand side of equation ( [ boltzmann ] ) is the scattering in term : it describes processes in which particles enter state @xmath19 due to the interaction with other particles .
this term can be written as @xmath79 , where @xmath80 includes all the terms on the right hand side of the equation that do not multiply @xmath1 .
we refer to @xmath80 as the `` generation '' term .
the second term is the scattering out term : it describes processes in which particles leave state @xmath19 .
this term is usually written as @xmath81 , where @xmath82 is the scattering _ out _ time from state @xmath19 .    in the context of the generalization of the be to a qke , it is more convenient to think of the collision integral @xmath77 as made up of `` generation '' and `` decay '' terms : @xmath83 , where the decay term @xmath84 includes all the terms containing @xmath1
. thus @xmath85 is the total decay rate of particles in state @xmath19 .
it is the _ total _ decay rate , and not only the scattering _ out _ rate , that is related to the retarded self energy @xcite .
the general form for the decay rate @xmath86 , as it is deduced from equation ( [ boltzmann ] ) , is @xmath87    2 \pi     \delta ( \epsilon_{\alpha } - \epsilon_{\beta } + \epsilon_{\gamma} ... ).\ ] ] in a d.c .
situation one obtains for the occupation numbers @xmath88    in order to transform the be into a qke one should apply the following rules . to the left hand side of equation ( [ boltzmann ] )
apply    * @xmath89 , where @xmath51 is the off - shell energy of particles in state @xmath19 .    to the right hand side of the equation apply the following steps    * the occupation numbers for particles in states @xmath90 are replaced with the occupation functions , that depend on the off - shell energies @xmath91 , corresponding to the on - shell energies @xmath92 respectively : @xmath93 , @xmath94 ... * the summand under the sum over all states @xmath75 is multiplied by the product of the smeared delta functions for particles in these states with the corresponding off - shell energies @xmath95 , and integrated over the off - shell energies : @xmath96 * the energy conservation delta function of the on - shell energies @xmath92 is replaced with an energy conservation delta function of their corresponding off - shell energies @xmath91 : @xmath97 . *
multiply the entire right hand side by @xmath98 .    in order to obtain an equation for @xmath67
one should apply the first three steps of the generalization rules for the right hand side of the equation , to the expression for @xmath86 ( [ decay ] ) .
the qke can be written generally in the form @xmath99 the function @xmath100 contains all the terms on the right hand side of the qke , that do not multiply @xmath72 .
it is the generating term - the term in the kinetic equation that is responsible for the particle generation rate .
equation ( [ qke - gg ] ) shows clearly that the particle level width @xmath67 is in fact the total particle decay rate .
note that due to ( [ norm ] ) the integration of equation ( [ qke - gg ] ) over @xmath51 results in the disappearance of the @xmath101 factor from the left hand side .
the occupation function in the time independent case is given by @xmath102    the representation of the qke in terms of @xmath71 and @xmath26 is more convenient than the representation in terms of green functions , since one can clearly see all physical processes , and use the physical intuition that one gained from the be in order to simplify the qke .
( 36,32)(0,0 ) ( 0,0 ) ( 3,1)(0,0)@xmath103 ( 7,29)(0,0)@xmath104 ( -1,27)(0,0)@xmath105 ( -1.2,19)(0,0)@xmath106 ( 0,35)(0,0)[b]@xmath107 ( -1.8,12)(0,0)@xmath108 ( -1.5,9)(0,0)@xmath109 ( 16.5,23)(0,0)[r]@xmath105 ( 20,17)(0,0)[br]@xmath110 ( 22.5,17)(0,0)[bl]@xmath111 ( 30,19)(0,0)@xmath112    we now use the formalism given in section [ qke ] in order to deal with the specific problem of hot luminescence from a qwr . the case in which electrons are excited from the highest valence subband to the lowest conduction subband via photons , and then relax via lo phonon emission ( see fig .
[ qwr ] ) is considered .
all other conduction and valence subbands are neglected ( in contrast to ref .
we consider a situation in which electrons are excited by a narrow band , noncoherent laser , just above the threshold for optical phonon emission @xmath113 , and describe the hot luminescence due to the recombination of electrons with holes in an impurity level .
it is assumed there is some mechanism due to which the electrons leave the bottom of the conduction subband , so that there is no electron accumulation there , but this recombination can be neglected in the balance equations near the threshold .
the spin degree of freedom will be completely disregarded .    unless specified otherwise we will assume that the electrons are excited above the threshold , since we are mainly interested in luminescence due to these electrons .
below the threshold the electrons relax to the bottom of the band by emitting acoustical phonons .
the relaxation rate due to acoustical phonons @xmath114 , is much smaller than that due to optical phonons , and can be neglected above the threshold .
the qwr runs along the @xmath115 direction so that the electron wave functions are localized in the @xmath116 plane , and free waves in the @xmath115 direction .
the wave function in the conduction band can be written as @xmath117 , where @xmath118 is the wire normalization length , and @xmath119 is a vector in the @xmath116 plane .
the electron energy is @xmath120 .
all energies are measured from the bottom of the conduction band .
the holes in the valence band have wave functions of the form @xmath121 and energy @xmath122 , with @xmath123 .
the exciting photons are taken to be plane waves and the photon frequency is @xmath124 ( @xmath125 is the photon momentum and @xmath126 is the light velocity ) .
the lo phonons are three - dimensional and have a flat dispersion law , @xmath127 .
the impurities are distributed randomly in the bulk .
the wave function of a hole in impurity @xmath128 is of the form @xmath129 , where @xmath130 are the coordinates of the impurity position .    in order to analyze the hot luminescence we neglect light polarization effects for simplicity , and assume that the wire and the crystal are cylindrically symmetric ( an assumption that would not fit v - grooved wires nor the assumption of ref .
thus , the wave functions of the emitted photons are @xmath131 where @xmath132 is the normalizing radius and @xmath118 the normalizing length of the crystal , and @xmath133 and @xmath134 are the cylindrical polar coordinates .
@xmath135 are the bessel functions , and @xmath136 is the @xmath137 zero of @xmath135 .
the luminescence photon frequency is given by @xmath138 ( where the prime signifies this frequency is of luminescence , and not of excitation ) .
we neglect all excitonic effects ( these are treated for the case of optical absorption in one - dimensional semiconductors in refs .
@xcite ) .
we begin by writing down the bes for the electrons in the conduction band , and for the luminescence photons .
it is shown why this description fails and one has to turn to the qkes .
we will then employ the generalization rules in order to obtain the qkes from the bes , and use the qkes in the analysis of the luminescence .
we assume low temperature and a weak excitation level .
at equilibrium , when there is no excitation , all the electrons are in the valence band , and there are no phonons ( @xmath139 ) , due to the low temperature .
in other words , there are no electrons in the conduction band ( @xmath140 ) , no holes in the valence band ( @xmath141 ) , and the impurity level is fully occupied by holes ( @xmath142 ) .    the excitation creates electrons in the conduction band , and holes in the valence band .
most of the excited electrons relax to the bottom of the conduction band , emitting lo phonons , but a small fraction of them recombine with holes in the impurity level , producing hot luminescence . by the assumption of weak excitation we mean that @xmath143 , @xmath144 and @xmath145 are small
. the luminescence would be even weaker than the excitation , therefore @xmath146 remains close to one . due to the weak excitation the luminescence is spontaneous , and the occupation numbers for the luminescence photons are also small .    in the derivation of the bes ,
terms that are nonlinear in the electron occupation numbers will be neglected .
we will also neglect nonequilibrium contributions to the phonon and the hole ( both in the valence band and in the impurity level ) occupation numbers .
the interaction of the electrons with the luminescence field is neglected in the balance equation for the electrons , in comparison with their interaction with the excitation field , due to the weakness of the excitation . in the case of c.w .
excitation the balance equation for the electrons is @xmath147 where @xmath148 is the matrix element for the electron - exciting - photon interaction , and @xmath149 is the matrix element of the electron - lo phonon interaction .
let us first look at the electron - exciting - photon matrix element . in this interaction
one photon of wave vector @xmath150 is absorbed , and a hole of momentum @xmath151 and an electron of momentum @xmath152 are created . the interaction , in second quantized form ,
is given by @xmath153 the mass @xmath154 is the bare electron mass and the constant @xmath155 is the `` bare '' electron momentum operator ( in the direction of the photon polarization , assumed linear for simplicity ) , sandwiched between the valence band top and the conduction band bottom bloch wave functions . @xmath156 and @xmath157 are the electron and hole field operators , respectively .
the photon field operator is given by @xmath158 where @xmath159 , @xmath160 are the components of @xmath150 in the @xmath116 plane , and in the @xmath115 direction , respectively , @xmath161 is the annihilation operator , and @xmath162 is the normalization volume of the crystal in which the quantum wire is embedded .    in order to find @xmath163 one should sandwich interaction ( [ int - exc ] ) between the states @xmath164 and @xmath165 .
this is equal to @xmath166 the photon momentum is small and so is @xmath119 , since it is limited by the wire cross section , therefore @xmath167 .
integrating over the @xmath115 coordinates the square of the matrix element is @xmath168    in the electron - lo phonon interaction one phonon of momentum @xmath169 is emitted , an electron @xmath152 is annihilated , and an electron @xmath170 is created .
the interaction of the electrons with the polarization created by the lo phonons is given by @xmath171 where @xmath172 is some interaction constant ( see for example ref .
@xcite , @xcite ) , and the phonon field operator is @xmath173 the components of @xmath169 are defined in the same manner as those of @xmath150 , @xmath174 is the annihilation operator , and @xmath175 is the reduced mass per unit cell in the lattice .
the matrix element @xmath149 can be found by sandwiching interaction ( [ int - lo ] ) between the states @xmath176 and @xmath177 .
carrying out the integration over the @xmath115 coordinates , the square of the matrix element is @xmath178 where we exchanged @xmath172 for the known frlich constant @xmath19 , using the relation @xmath179 ( see ref .
@xcite ) .    substituting @xmath180 and @xmath181 into equation ( [ e - be ] )
, we can sum over @xmath151 and @xmath169 . taking the normalizing volumes to infinity and thus exchanging the sums with integrals
we obtain @xmath182 where @xmath183 @xmath184 and @xmath185 the zeroth order bessel function @xmath186 is a result of the integration over @xmath187 .
since the matrix elements ( [ |m - e - ep|^21 ] ) and ( [ |m - e - ph|^21 ] ) are smooth functions of their arguments , and we are interested only in processes which involve electrons that were excited close to the threshold , we may take them at the threshold values .
the threshold value of @xmath188 is @xmath189 , where @xmath190 .
this is the frequency that will excite electrons from the valence band to the threshold in the conduction band .
the threshold value of @xmath191 is @xmath192 .
the electron wave functions that appear in @xmath5 limit the spatial integration to @xmath193 , where @xmath194 is the wire width , therefore we can estimate that @xmath195 for @xmath196 and @xmath197 for @xmath198 .
for gaas @xmath199 , thus for most wire widths we are in the regime of @xmath198 , in which @xmath5 is expected to decrease with increasing wire cross section .
the constant @xmath200 is the effective frhlich constant in the one - dimensional case .
according to ref .
@xcite @xmath201 ( which is of the same order as the bulk frhlich constant ) for a wire of cross section @xmath202 , while according to ref .
@xmath203 for a wire of cross section @xmath204 .
we now perform the integrations which appear in equation ( [ e - be1 ] ) .
we begin with the generating term @xmath205 , the free term ( one that does nt include @xmath144 ) in the balance equation ( [ e - be1 ] ) for the electron occupation number .
this term is physically the generation rate of electrons in the conduction band . since @xmath206 we can approximate @xmath207 . using @xmath208 , where @xmath209 stands for the solid angle increment
, we obtain @xmath210 where @xmath211 , is the photon density of states , and @xmath212 is the angular average of @xmath213 , for photons of frequency @xmath214 , that create electrons with energy @xmath215 .
we consider a narrow band photon excitation field such that the angular average of the occupation number is given by @xmath216 where @xmath217 is the excitation field energy density , @xmath218 is the spectral width , and @xmath219 is the central frequency .
we are interested in a narrow band excitation close to the threshold , therefore @xmath220 , and the detuning ( the difference between the central excitation frequency and the threshold frequency ) is small : @xmath221 .
the integrations that appear in the expression for the decay rate @xmath222 , the term that includes all coefficients of @xmath223 in the kinetic equation ( [ e - be1 ] ) , will now be carried out . the first contribution to @xmath222
is that of the photons , and is equal to the generating term ( [ e - be1-ep ] ) .
the second contribution to @xmath222 is due to lo phonons , @xmath224    the photon contribution to the decay rate is negligible compared to the phonon contribution since radiative processes are slow compared to non - radiative processes . therefore @xmath222 can be written as @xmath225 where @xmath226 the meaning of the energy scale @xmath227 will be explained at the end of this section . for @xmath228 @xcite one obtains @xmath229 , thus @xmath230 .    below the threshold , where @xmath231 , the acoustical phonon contribution to @xmath27 is important . since
the latter is a smooth function near the threshold , a good approximation for @xmath27 is @xmath232 at @xmath233 .
the occupation number of the electrons that were excited above the threshold is given by expression ( [ n - a ] ) , where @xmath234 and @xmath27 are given by expressions ( [ e - be1-ep ] ) and ( [ decay - ek ] ) respectively .
it is clear from the expressions above that @xmath144 depends on @xmath152 through @xmath215 only .
we now turn our attention to the kinetic equation for the luminescence photons , @xmath235 , \end{aligned}\ ] ] which for @xmath142 and in case of spontaneous luminescence is @xmath236 @xmath237 is the matrix element for the luminescence photon - electron interaction .
the time derivative in the equation above is kept in order to clarify what is the luminescence source .
let us first look at the matrix element that describes the recombination of an electron with a hole from a specific impurity .
we will then have to sum over all impurities and average over all impurity configurations .
in this process a photon of quantum numbers @xmath238 is emitted , and an electron @xmath152 and a hole in impurity @xmath128 are destroyed .
the interaction is given by @xmath239 where @xmath240 is the field operator of the luminescence photons .
the matrix element of the interaction above between the states @xmath241 and @xmath242 for holes in impurity @xmath128 is @xmath243    we now sum the matrix element squared over all impurities and average over all impurity configurations .
this yields @xmath244 carrying out the integrations over the @xmath115 coordinates we obtain @xmath245 where @xmath246 . the second argument of @xmath247 ,
@xmath248 was replaced by @xmath249 because @xmath206 .
since the impurities are randomly distributed in the bulk , the average over all impurity configurations of the product of the hole wave functions , that appear in ( [ |m - lp - e|^21 ] ) , can depend only on the difference between their coordinates .
thus we define @xmath250 the summation over impurities divided by the normalizing volume renders a factor of @xmath251 - the impurity spatial density .
the averaged matrix element squared can then be written as @xmath252    the matrix element squared is much larger for @xmath253 than for other values of @xmath254 , since @xmath134 and @xmath255 are constrained to a small region in the @xmath116 plane ( due to the electron wave functions ) , and @xmath256 is the only bessel function which is finite at @xmath257 , @xmath258 .
thus the matrix element squared is proportional to @xmath259 .
from now on the index @xmath254 is omitted with the understanding that we are dealing only with @xmath253 .
the matrix element squared is a smooth function of @xmath260 and @xmath152 , therefore these can be substituted by their threshold values .
the threshold values of @xmath260 and @xmath152 are @xmath261 and @xmath262 respectively . the kinetic equation for the luminescence photons can then be written as @xmath263 where the matrix element squared is @xmath264    substituting @xmath265 into equation ( [ lp - be1 ] ) and integrating over @xmath152 we obtain @xmath266    the luminescence source can be characterized by the spectral dependence of @xmath267 , the energy of the emitted field , of all spectral modes within the interval @xmath268 , per unit length of the wire , per unit time . multiplying the generating term of the luminescence photons ( the right hand side of equation ( [ lp - source ] ) ) by @xmath269 and summing over @xmath160 and @xmath69 , we obtain @xmath270    the variable @xmath260 in the square root can be replaced with its threshold value @xmath271 , due to the smoothness of this function . taking the normalization volume ( that appears in the matrix element ) to infinity the sums over @xmath160 and @xmath69
are transformed into integrals .
performing the integration we obtain the final result @xmath272 where @xmath273 is a constant , and @xmath274 is essentially the product of the electron generation rate and the electron life - time at the corresponding energy .
since @xmath275 is a smooth function of @xmath276 , the spectral distribution of the luminescence is given by @xmath277 .
let us first consider a narrow band excitation that does not overlap the threshold , i.e. @xmath278 . in this case
@xmath277 is a peaked function of @xmath276 , reproducing the shape of the excitation .
it is of a lorentzian shape of width @xmath279 , centered at @xmath280 , where @xmath281 ( see fig .
[ xi1 ] ) @xmath282 the rescaling of @xmath276 compared to @xmath283 follows from the obvious relations ( see fig . [ qwr ] ) @xmath284 and @xmath285 .
( 10,7)(0,0 ) ( 0,0 ) ( 6,-0.1)(0,0)[t]@xmath286 ( 0,6.5)(0,0)[b]@xmath287 ( 10.3,0)(0,0)[lb]@xmath288 ( 0,0)(0,0)[tr]@xmath106 ( 5.6,2.75)(0,0)[r]@xmath289 ( 6.4,2.75)(0,0)[l]@xmath290 ( 7.2,2.75)(0,0)[l]@xmath291 ( 7.3,4.7 ) ( 7.3,6.5)(0,0)[b]@xmath287 ( 10,4.7)(0,0)[lb]@xmath288 ( 7.3,4.7)(0,0)[tr]@xmath106    only close to the threshold , when @xmath292 , the spectral behavior of the luminescence differs from that of the excitation and is given by ( see the inset of fig . [ xi1 ] ) @xmath293    ( 4,3)(0,0 ) ( 0,0 ) ( 1.25,0)(0,0)[t]@xmath294 ( 0,2.6)(0,0)[b]@xmath287 ( 4.1,0.1)(0,0)[l]@xmath288 ( 0,0)(0,0)[tr]@xmath106    when the narrow band excitation overlaps the threshold , i.e.   @xmath295 , the spectral distribution of the luminescence is very different from that of the excitation ( see fig . [ xi3 ] ) : @xmath296 and @xmath297    if one is interested in luminescence due to electrons that were excited close to but _ below _ the threshold , via an excitation that is centered _ above _ the threshold , i.e.  @xmath298 and @xmath299 , then @xmath300 in expression ( [ xi - be ] ) should be replaced by @xmath114 . in this case
@xmath277 has the same lorentzian form as ( [ xi - lor ] ) with the amplitude @xmath301 replaced by @xmath302 .
one can expect that the luminescence distribution @xmath277 , derived from the be , is not correct for frequencies @xmath276 close to the threshold @xmath303 .
luminescence of such frequencies is due to electrons with energies @xmath215 near the threshold , that have a large width @xmath304 ) , ( see equation ( [ decay - ek ] ) ) , while the be assumes that the width of the states is small compared to other energy scales of the problem . in our case one of these energy scales is @xmath305 . equating the width with this energy scale , @xmath306 ,
one finds an additional energy scale of the problem ( and a critical decay rate ) , that is given by @xmath227 that has been defined before ( [ gamma - c ] ) .
when @xmath307 , we have @xmath308 , and the assumption of the be breaks down .
one would therefore expect that for @xmath307 in the description of the excited electrons , and for @xmath309 in the description of the hot luminescence , the results predicted by the boltzmann equation would fail .
we will learn from the qkes that in fact the situation is more complicated .      as we have seen in the previous section
, the boltzmann description may fail for electrons excited to the vicinity of the threshold and hence for luminescence photons emitted by these electrons .
therefore , close to the threshold , the electrons and the luminescence photons have to be treated using qkes .
these qkes will be written by applying the generalization rules to equations ( [ e - be1 ] ) and ( [ lp - be1 ] ) .
it is assumed that the energies of the holes , the phonons and the excitation photons are not broadened .    in the generalization to the qkes
only those terms that were kept in the boltzmann description will be retained .
threshold values of the arguments of the matrix elements will be substituted as was done in the boltzmann description , due to the smoothness of the matrix elements .
the qke for the electrons is @xmath310 where @xmath51 is the electron off - shell energy , @xmath69 is the electron occupation function , and @xmath27 and @xmath311 are the electron energy level width and shift , respectively .    the coefficients of @xmath312 in equation ( [ e - qke ] ) can be recognized as the electron level width due to lo phonons @xmath313 and the total electron width is given by @xmath314 the photon contribution to the level width was neglected ( as it was neglected in the boltzmann description of the decay term ) .
one can see that the generating term in equation ( [ e - qke ] ) is equal to expression ( [ e - be1-ep ] ) , the generating term of the electrons be , with the on - shell energy @xmath215 replaced by the off - shell energy @xmath51 ( in the exciting photon occupation number ) @xmath315    the occupation function @xmath316 is found by substituting @xmath234 and @xmath27 into ( [ n - ae ] ) .
the generating term @xmath317 depends on @xmath152 through @xmath318 only , therefore @xmath319 .
it follows from ( [ gamma - e ] ) that @xmath27 is independent of @xmath152 ( and so is the level shift ) .
therefore , @xmath320 .    equation ( [ gamma - e ] ) is a self consistent equation for the level width . integrating the right hand side of the equation
we find that @xmath321^{-1/2 } +     \right .
\\    \left . \left
[ \epsilon - \omega_{o } -           \delta \epsilon ( \epsilon - \omega_{o } ) -            \frac { i}{2 }            \gamma ( \epsilon - \omega_{o } ) \right]^{-1/2 }     \right\}.\end{aligned}\ ] ] when @xmath322 the arguments of @xmath323 and @xmath324 are close to the bottom of the band . in this case
the level shift can be neglected since it is just a small renormalization .
the level width close to the bottom of the band @xmath325 , is due to scattering with acoustic phonons and due to thermal recombination .
the latter contribution to @xmath325 can be neglected , while the acoustic phonon contribution can be considered as a constant ( though not the same as the width due to acoustic phonons at the threshold ) .
we assume that @xmath325 is small compared to all other energy scales of the problem .
thus close to the threshold @xmath326^{-1/2 } +            \left [ \epsilon - \omega_{o } -            \frac { i}{2 }             \gamma_{o } \right]^{-1/2 }     \right\}.\ ] ] the behavior of @xmath327 when @xmath51 is close to @xmath328 is shown in fig .
[ gamma ] .
( 22,15 ) ( 0,0 ) ( 1,-0.5)(0,0)[t]@xmath329 ( 10,-0.5)(0,0)[t]@xmath330 ( 12.3,-0.5)(0,0)[t]@xmath331 ( 21.3,-0.5)(0,0)[t]@xmath332 ( 11,14.5)(0,0)[b]@xmath333 ( 22.5,0)(0,0)[lb]@xmath334 ( 11.9,1.2)(0,0)[l]@xmath335    from expression ( [ gamma - th ] ) we see that above the threshold for @xmath336 , @xmath337 can be written as @xmath338 comparing expression ( [ gamma- < ] ) to expression ( [ decay - ek ] ) , it is evident that above the threshold @xmath337 is equal to @xmath222 when the on - shell energy @xmath215 is replaced by the off - shell energy @xmath51 . since @xmath325 is the smallest energy scale in the problem , expression ( [ gamma- < ] ) can be used for the electron level width above the threshold .
contrary to the boltzmann description , the electron - lo phonon scattering contributes to the electron level width below the threshold as well . when @xmath339 the electron width due to optical phonons decays like @xmath340 , much faster than the decay above the threshold ( see fig . [ gamma ] ) . the contribution of the acoustic phonons to the electron width below the threshold @xmath114 , will be negligible compared to @xmath341 for @xmath342 . since @xmath343 and @xmath344 the right hand side of the inequality above is much smaller than @xmath227 , and much larger than @xmath325
( see the crossing point of @xmath345 and @xmath114 in fig .
[ gamma ] ) .
the qke that is obtained by applying the generalization rules to the be for the luminescence photons ( [ lp - be1 ] ) is @xmath346 where @xmath276 is the photon off - shell energy , @xmath347 is the photon occupation function , and @xmath172 and @xmath348 are the photon energy level width and shift , respectively .
all the photon functions are written for @xmath253 .    by applying the generalization rules to the right hand side of
the full be ( [ lp - be ] ) one can find the photon level width @xmath349 note that @xmath350 corresponds to photon generation .
it is evident from expression ( [ gamma ] ) that @xmath172 does not depend on @xmath160 , and therefore neither does @xmath348 .    using ( [ gamma ] ) equation ( [ lp - qke ] ) can be written as @xmath351 the right hand side of equation ( [ lp - qke1 ] ) is the photon generating term , the term from which the luminescence spectral distribution will be obtained .
as in the classical case , the luminescence source can be characterized by the spectral dependence of @xmath352 , that is obtained by multiplying the right hand side of equation ( [ lp - qke1 ] ) by @xmath353 and summing over @xmath160 and @xmath69 @xmath354 the main contribution to @xmath355 comes from @xmath160 and @xmath69 such that @xmath260 is close to @xmath276 .
it is convenient to write @xmath356 in the following form @xmath357 where @xmath358 is given by ( [ c ] ) , and @xmath359 when the electrons energy is taken to be on shell , that is the electron smeared delta function is reduced to a singular delta function , @xmath277 is reduced to expression ( [ xi - be ] ) . in writing @xmath172 in the form ( [ gamma-1 ] ) , the normalization volume has been separated from the rest of the function , that is independent of the index @xmath69 , and can therefore be taken out of the summation in @xmath360 .    taking the normalization volume to infinity the sums over @xmath160 and @xmath69
are transformed into integrals . from equation ( [ gamma-1 ] ) it is evident that @xmath361 .
this is due to the fact that the photons are emitted into an `` infinite '' space , while their interaction with the electrons is confined to the finite volume of the wire . in this case @xmath172 and @xmath348 are negligible compared to the other widths , and the smeared photon delta function that appears in the expression for @xmath360 becomes a singular delta function , giving @xmath362 since @xmath363 is a smooth function it is clear at this point that in order to understand the dependence of @xmath360 on @xmath276 , the behavior of @xmath364 should be analyzed .    performing the integration over @xmath51 in ( [ bar - gam ] ) , and substituting the expression we got for @xmath365
, we get @xmath366 since @xmath367 , we substitute @xmath328 for @xmath215 in the one - dimensional density of states .
this integral contains in fact two lorentzians .
the first is of width @xmath368 : the electron level width at energy @xmath369 .
the second lorentzian is of width @xmath370 : a width that is proportional to the spectral width of the exciting photon source .
since the excitation is such that @xmath215 is close to @xmath328 , and all the energy scales that can characterize the level width are much smaller than @xmath328 , the integrand goes to zero when @xmath215 approaches zero from above .
thus , the lower boundary of the integral can be taken to minus infinity . as a result
the spectral distribution of the luminescence is given by the function @xmath371^{2 } + \delta \nu_{o}^{2}/4 } \right).\end{aligned}\ ] ]    the region of validity of the be is evident from equation ( [ bar - gam2 ] ) . in order to obtain the result of the be ( [ xi - be ] ) , one has to replace the first lorentzian with a delta function .
this can be done only when its width is smaller than that of the second lorentzian , i.e. @xmath372 .
it follows from this inequality that the boltzmann description of the luminescence spectra is correct only far from the threshold , when @xmath373 .
the quantum interval @xmath374 differs from the naive estimate @xmath227 ( see the discussion at the end of subsection [ b - des ] ) : it is larger for `` narrow '' band excitation ( @xmath375 ) , and smaller for `` wide '' band excitation ( @xmath376 ) .    for an excitation that is mostly of frequencies above @xmath377
, the spectral dependence of the energy of the emitted photons will behave according to the predictions of the be ( see figs .
[ xi1]-[xi3 ] ) .
however , as long as the excitation is above the threshold , very close to the threshold @xmath277 will increase linearly with @xmath276 , and not as a square root , see ( [ sqrt ] ) .
non - boltzmann behavior is obtained when @xmath378 , that is @xmath379 , and the second lorentzian can be treated as a delta function . in this case one
obtains @xmath380^{-1}.\ ] ] we bring here two specific examples of extreme non - boltzmann behavior of @xmath277 . in both examples the excitation is centered within the quantum interval , i.e. @xmath381 , and `` narrow '' band , i.e. @xmath375 . due to the latter inequality @xmath382 .    in the first case
the detuning is large , @xmath383 , therefore the energy width of an electron excited by the central excitation frequency , @xmath384 , that is of the order of @xmath385 , is small compared to its distance from the threshold .
this is true for most of the electrons .
in such a case the electron states are `` well defined '' , but the prediction of the be for the electron distribution is wrong , since its width @xmath386 is much smaller than the width of the states . as a result it follows from ( [ xi - q ] ) that the luminescence spectra is symmetric and centered at the classical position @xmath280 , with @xmath387 , but its width is given by the quantum - mechanical width of the electron states ( see fig .
[ xi - q1 ] ) @xmath388^{1/2}.\ ] ]    ( 20,10)(0,0 ) ( 0,0 ) ( 0,11.5)(0,0)[b]@xmath287 ( 18.5,0)(0,0)[lb]@xmath288 ( 11.3,-0.1)(0,0)[t]@xmath389 ( 0,0)(0,0)[tr]@xmath106 ( 10.5,4.3)(0,0)[r]@xmath289 ( 12,4.3)(0,0)[l]@xmath290 ( 13.5,4.3)(0,0)[l]@xmath291 ( 13,8 ) ( 13,11.5)(0,0)@xmath287 ( 13,8)(0,0)[tr]@xmath106 ( 18.5,8)(0,0)[lb]@xmath288    ( 30,16)(0,0 ) ( 0,0 ) ( 0,16)(0,0)[b]@xmath287 ( 25.5,0)(0,0)[lb]@xmath288 ( 15,0)(0,-0.5)[t]@xmath332 ( 0,0)(0,0)[tr]@xmath106    in the second case the detuning is small , @xmath390 , and the electrons are excited to `` badly defined '' states , since the broadening of these states is larger than their distance from the threshold . as a result
the luminescence spectral distribution differs greatly from that predicted by the be .
one finds from ( [ xi - q ] ) that ( see fig .
[ xi - q3 ] ) @xmath391 and @xmath392    in conclusion let us note that the hole dispersion plays an essential role in determining the quantum behavior of the luminescence spectra .
the width of the second lorentzian in ( [ bar - gam2 ] ) @xmath370 , is not equal to the width of the classical luminescence line @xmath393 , and the integral is not a simple convolution of the classical luminescence profile with the spectral function of an electron state .
one can see from ( [ bar - gam2 ] ) that for a flat hole band ( @xmath394 ) the second lorenzian does not depend on the integration variable and the first lorenzian is integrated to one , restoring the be result , unexpectedly .
the luminescence due to electrons that were excited _ below _
the threshold ( @xmath395 ) by an excitation centered _ above _ the threshold ( @xmath299 ) is not obtained by the simple exchange of @xmath341 by @xmath114 , as in the classical case . as long as @xmath396 one should take @xmath397 from expression ( [ gamma - th ] ) since @xmath398 , and therefore although the electrons are below the threshold , it is the optical phonon contribution to the electron width that is dominant .
we have presented a recipe that allows one to generalize the be to a qke .
the qke obtained by employing this recipe is the same as that obtained by the keldysh green function technique in the self - consistent born approximation .
the advantage of this method for writing the qke is that it provides a physical understanding of the terms in the qke , and it allows one to neglect those terms that were negligible in the boltzmann description .
this is due to the fact that the equations are written for quantities that are similar to the quantities described by the boltzmann equation .
we considered the specific example of hot luminescence from a qwr .
the recipe described above was used in order to generalize the set of bes that describe the problem to a set of qkes . solving these equations we were able to describe the luminescence spectral distribution .
we have shown that there is a domain of luminescence frequencies , that correspond to a domain of phtoexcited electron energies , for which the quantum description of the luminescence spectral distribution leads to a different behavior than that given by the boltzmann description .
this quantum domain could not be easily guessed from level width considerations .
two other nontrivial conclusions were obtained .
the first is the role played by the hole mass in the definition of the quantum domain .
when the hole dispersion relation is flat the quantum domain shrinks to zero and the results of the boltzmann description are retrieved .
the second nontrivial conclusion is that there is an energy domain below the threshold for lo phonon emission , in which the lo phonon contribution to the electron level width is dominant .