the broadening of exciton emission lines due to electron - exciton collisions was first observed in bulk gaas @xcite . an experimental comparison between exciton - electron and exciton - exciton scattering mechanisms in bulk gaas @xcite and in gaas quantum wells ( qw s )
@xcite showed that the exciton - electron scattering efficiency is an order of magnitude larger than that of the exciton - exciton process .
both scattering processes are enhanced for the 2d excitons as compared to bulk excitons .
theoretical calculations of elastic and inelastic exciton - electron scattering were reported for bulk semiconductors @xcite and for qw s @xcite .
the semiclassical treatment of the latter neglected the exchange term which was shown recently to be the dominant term in qw s @xcite .
the existence of the negatively charged exciton @xmath1 ( trion ) was first proposed by lampert @xcite as the semiconductor analogue of the hydrogen ion .
due to their confinement in the growth direction , 2d trions have a binding energy which is an order - of - magnitude larger than that of the bulk trions @xcite .
this fact facilitated the observation of trions in cdte @xcite and in gaas @xcite qw s containing a low - density two - dimensional electron gas ( 2deg ) .
theoretical calculations of the binding energy of both negatively ( @xmath1 ) and positively ( @xmath4 ) charged excitons were performed using various trial functions , showing that only the singlet state is bound @xcite .
particular attention was devoted to the modification of the trion s properties in the presence of strong magnetic fields . in this limit ,
the @xmath1 binding energy is increased , and the electron spin triplet state becomes bound as well @xcite .
it was also observed that the dependence on the magnetic field of the @xmath1 and @xmath4 binding energies , as well as their zeeman splitting , differ drastically ( while being nearly identical at zero magnetic field ) @xcite .    while these aspects of the trions where extensively studied , there are very few reports on their broadening mechanisms @xcite . in this paper
we present a theoretical model for the scattering of both neutral and charged excitons with electrons in qw s .
this paper continues our previous work which introduced neutral excitons elastic scattering with electrons in the context of cavity polaritons @xcite .
here we extend our model , to incorporate both excitons and trions , taking into account elastic and inelastic scattering processes .
the model is based on calculating the exciton - electron direct and exchange interaction matrix elements , from which we derive the exciton scattering rates .
these are integrated over all final excitonic states , resulting in the exciton linewidth due to scattering , as a function of its initial momentum ( or energy ) .
a major difference between the exciton and trion scattering lies in the charge of the latter , which results in a divergence of its matrix elements in the limit of zero transferred momentum .
this divergence , originating from the infinite range of the coulomb potential , is treated by applying the lindhard model for the potential screening .
the screening action complicates the trion linewidth dependence on the electron density ; thus , while for electron densities larger than @xmath5 the trion linewidth due to electron scattering in a gaas qw at @xmath3 is comparable to that of the neutral exciton , for a very dilute 2deg ( @xmath6 ) the trion linewidth becomes roughly an order - of - magnitude larger than that of the neutral exciton .
inelastic scattering is also much more efficient in the trion s case , due to its smaller binding energy .
the calculated reflection spectra are obtained by considering the three electronic excitations of the system , namely , the heavy - hole and light - hole 1s neutral excitons , and the heavy - hole 1s charged exciton ( trion ) , with the appropriate oscillator strengths . a qualitative validation of our calculations is given by considering photoluminesence ( pl ) measurements that were done on a mixed type i  type ii gaas / alas qw ( mtqw ) structure @xcite .
the paper is organized as follows . in section [ ex - el ] ,
we present in detail the model for the neutral exciton - electron scattering . in section [ trion ]
we consider the trion - electron scattering .
first we construct the trion s wave function , which is used in the matrix elements calculations .
we then treat the divergences that originate from the infinite range of the coulomb potential , by applying the lindhard model for the potential screening .
the effect of the screening is discussed in the context of the trion linewidth dependence on the electron density .
in section [ bare - ex ] we use our model to calculate qw reflection spectra and discuss their relevance to the available experimental data . a summary and conclusions
are given in section [ conc ] . in appendices
a and b , we provide some details of the direct and exchange integrals calculations for the neutral exciton - electron scattering .
finally , a derivation of the trion binding energy , using a chandrasekhar - type trial function , is given in appendix c.
in this section we present a detailed description of the exciton - electron scattering , considering separately the elastic process , for which the exciton remains bound , and the inelastic process , where the exciton breaks into an unbound electron - hole pair .
we first construct the state of a single exciton with an in - plane center - of - mass ( cm ) momentum @xmath7 in the fermionic hilbert space of electron - hole pairs , using the notations of tassone and yamamoto @xcite .
it is a superposition of wave functions with different electron momenta @xmath8 and electron and hole @xmath9 coordinates , given by @xmath10 where @xmath11 is the electron ( hole ) creation operator with in - plane momentum @xmath7 and @xmath12 ( @xmath13 ) coordinate , @xmath14 , @xmath15 are the electron and exciton in - plane effective masses , respectively , @xmath16 , and @xmath17^{-3/2}\chi_e(z_e ) \chi_h(z_h)$ ] is the in - plane fourier transform of the exciton wave function ( see appendix a ) .
a state comprising an exciton and an unbound electron having @xmath18 and @xmath19 will be written as @xmath20 where @xmath21 is the electron wave function . applying the coulomb interaction operators @xmath22 , @xmath23 to this state , we have for the electron - hole and electron - electron interactions :    [ vehvee ] @xmath24    the first term in eq .
( [ veh ] ) represents the coulomb interaction between the constituents of the exciton , thus contributing to its self energy , and can be discarded in the calculation of the scattering matrix elements . in eqs.([veh ] ) and ( [ vee ] ) , @xmath25 is the 2-dimensional fourier transform of the coulomb interaction , where @xmath26 is the qw area , and the plus ( minus ) sign is used in @xmath22 ( @xmath23 ) . using the anticommutation relations for the fermion operators , one can easily find the scattering matrix elements to be    [ matel ] @xmath27\end{aligned}\ ] ] @xmath28,\end{aligned}\ ] ]    where we denote @xmath29 . the first term in each of eqs .
( [ matel ] ) contributes to the direct ( classical ) coulomb interaction and the second contributes to the exchange matrix element .
the direct term is calculated analytically ( the details can be found in appendix a ) , resulting in : @xmath30 , \label{vdir}\ ] ] where @xmath31 is the qw width , @xmath32 is a variational parameter associated with the exciton bohr radius in the qw , and we have defined the dimensionless functions @xmath33^{-3/2 } \nonumber \\
e}^{-lq}- 1+lq + \frac{5(lq)^3}{8\pi^2 } + \frac{3(lq)^5}{32 \pi^4}}{(\lambda q)^3 \left [ 1+(lq/2 \pi)^2 \right]^2}. \label{gh}\end{aligned}\ ] ] it is evident from eq .
( [ vdir ] ) that the direct term is identically zero for equal electron and hole masses ( @xmath34 ) .
similarly , we have for the exchange term @xmath35 h(\lambda k_1 ) , \label{vexch}\end{aligned}\ ] ] where we have defined @xmath36 . the exchange term ( [ vexch ] )
is computed numerically , however its angular part can be calculated analytically ( see appendix b ) .
it is convenient to transform to dimensionless direct and exchange integrals given by @xmath37 the direct and exchange integrals are plotted in fig .
[ dir_exc ] as a function of the transferred momentum @xmath38 , for the case @xmath39 , where the angular dependence of @xmath40 disappears . as @xmath41 the direct integral approaches zero ,
while the exchange integral has its maximum ( this is also the case for exciton - exciton interaction @xcite ) . in the general case ,
the exchange integral is a function of the transferred momentum @xmath38 , the momentum difference @xmath42 ( which can be regarded , for convenience , as the in - plane momentum of the colliding electron in the rest frame of the exciton ) , and of the angle @xmath43 .
the exchange interaction term has the following features : ( i ) the interaction favors the case @xmath44 .
physically , this means that the electron is inclined to transfer as much momentum as possible to the exciton , preferably in the same direction .
( ii ) the interaction retains its strength for quite large values of @xmath38 ( or @xmath42 ) even though the excitonic wave function vanishes much more rapidly with momentum .
( iii ) the differential cross section is largest for @xmath45 and decreases to a minimum for back scattering ( @xmath46 ) .
= 0.55    next we calculate the scattering rate of an exciton , with an initial in - plane momentum @xmath7 , to an excitonic state with an in - plane momentum @xmath47 .
the exchange matrix element ( [ vexch ] ) depends on five variables : @xmath48 , where we have defined @xmath49 .
it is convenient to transform from the angles @xmath50 and @xmath51 to the angles @xmath52 and @xmath53 @xcite : @xmath54 it can be verified that for the special case of @xmath55 , we have @xmath56 and @xmath57 .
conservation of energy and momentum for the exciton - electron scattering process reads @xmath58 where @xmath59 ( @xmath60 ) is the exciton ( electron ) kinetic energy .
this equation is satisfied with : @xmath61 \ , \ ] ] the scattering rate is then found using fermi s golden rule    the scattering rate is then found using fermi s golden rule @xmath62 \nonumber \\ & = & \frac{a m_e}{\pi \hbar^3 } \int_0^{2 \pi } \frac{d \delta } { \left| \cos \delta \right|}\frac{k_{\rm e}^*}{q } \left|v_{\rm dir } + v_{\rm exc } \right|^2   f_{\rm fd } ( k_{\rm e}^*)\left[1-f_{\rm fd } ( |{\bf k}_{\rm e}^*-{\bf q}| ) \right ] , \label{rate_ex}\end{aligned}\ ] ]    where @xmath63^{-1}$ ] is the electron fermi - dirac distribution function .
in the two dimensional case that we are considering , the chemical potential is @xmath64 , and we assume that the fermi energy is given by the free 2d electrons value : @xmath65 .
we note that in the high temperature limit , @xmath66 , the fermi - dirac distribution is practically classical , and that for low electron densities , the chemical potential becomes negative .    using the relations ( [ gamdel ] ) , the integrand in ( [ rate_ex ] )
can be expressed as a function of the variables : @xmath67 .
the exciton linewidth due to electron scattering is calculated by integrating over all final excitonic states @xmath68    figures [ fgam_ex ] a , b show the heavy - hole exciton ( @xmath69 ) and light - hole exciton ( @xmath70 ) linewidths at @xmath71 and @xmath3 for @xmath72 .
the large linewidths obtained for relatively low @xmath73 reflect the high efficiency of the electron scattering mechanism .
this should be compared to an exciton linewidth of @xmath74 0.2mev for acoustic phonon scattering at @xmath71 .
the larger linewidth obtained at @xmath3 ( fig .
[ fgam_ex]b ) is explained by noticing that the exciton - electron interaction matrix elements favor small energy transfer transitions .
thus , at high temperatures , the electrons have too high energy to be effective scatterers . decreasing
@xmath75 thus increases the scattering rate until the electron gas becomes degenerate , and fewer final states are available for the scattered electrons . in figure
[ fgam_ex]c the heavy hole exciton linewidth at @xmath55 is plotted as a function of @xmath73 for @xmath3 and @xmath76 ( the light hole exciton linewidth exhibits a similar behaviour ) .
it is seen that while at @xmath71 the linewidth amplitudes scale practically linearly with @xmath77 for a large range of electron densities ( up to @xmath78 ) , the linearity region at @xmath3 is much smaller .
this linearity threshold seems to be in accord with the onset of the phase - space filling effect which becomes noticeable at higher densities and effectively enlarges the exciton bohr radius @xcite . as the temperature decreases , the effect of phase - space filling becomes important at a much lower electron density .
we note that the functional dependence of the linewidths on the initial exciton momentum is almost unchanged within the linear regime . increasing @xmath73 further results in a shift of the maximum linewidth from @xmath55 to higher momenta .
= 0.55      we now consider the case when the electron - exciton scattering results in dissociation . this scattering process is depicted schematically in fig .
[ xion_scat ] .
( 150,150 ) ( 0,0)(150,150 )    ( 13,116)x ( 25,118 ) ( 30,115)(1,-1)30    ( 20,135)@xmath7    ( 15,30)e ( 25,32 ) ( 30,35)(1,1)30    ( 20,10)@xmath79    ( 80,90)h ( 90,92 ) ( 95,95)(1,1)30    ( 87,135)@xmath80    ( 80,73)e ( 90,75 ) ( 95,75)(1,0)40    ( 110,80)@xmath81    ( 80,55)e ( 90,58 ) ( 95,55)(1,-1)30    ( 100,10)@xmath82    the initial state of the system is given again by eq .
( [ initial ] ) while the final state is that of three free particles , symbolically written as @xmath83 , where @xmath84 is the momentum transferred from the free electron to the electron - hole pair , and @xmath81 is the second electron momentum . calculating the coulomb interaction matrix elements between the initial and final states of the system ,
we again obtain direct and exchange interaction terms :    @xmath85 \\
v_{{\rm exc}}^{\rm dis}\!\!\!&(&\!\!\!{\bf k}_{\rm x},{\bf q},{\bf k}',{\bf k}_{\rm e } ) = -\frac{4 e^2 \lambda^4}{\epsilon_0 l^2}\left(\frac{2 \pi}{a}\right)^{3/2 } h(\lambda \left|{\bf k}'-{\bf k}_{\rm e}\right| ) \times \nonumber \\ & &   \left [ g(\lambda|\alpha { \bf k}_{\rm x}+{\bf q}-{\bf k}'| ) -g(\lambda|\alpha { \bf k}_{\rm x}+{\bf q}-{\bf k}_{\rm e}| ) \right]\end{aligned}\ ] ]    where @xmath86 are given in ( [ gh ] ) .
we note that contrary to the elastic scattering case , here the direct and exchange matrix elements are comparable . in order to calculate the scattering rate we first introduce energy - momentum conservation @xmath87 where @xmath88 is the exciton binding energy .
denoting the angle @xmath89 together with @xmath90 which were defined previously , this equation is satisfied with @xmath91 \right\}\end{aligned}\ ] ] the scattering rate is then calculated to be : @xmath92 \times \\ & & \left [ 1-f_{\rm fd}(k ' ) \right ] \nonumber , \label{xrate_ion}\end{aligned}\ ] ] where the + ( - ) sign between the direct and exchange matrix elements corresponds to the singlet ( triplet ) electron spin configuration . for an unpolarized electron gas , an averaging over the two spin configurations
must be performed , in order to take into account the contributions of all the electrons to the exciton linewidth due to dissociating scattering : @xmath93 @xmath94 ( @xmath95 ) denotes the exciton linewidth contribution from the singlet ( triplet ) spin configuration , and both are calculated using eq .
( [ gam_ex ] ) .
we note that the spin configuration averaging is irrelevant in the elastic scattering case , since the direct term there is much smaller than the exchange term .
this results in similar contributions from both singlet and triplet configurations .
= 0.65    figure [ lw_ion]a shows the heavy hole exciton linewidth due to dissociating scattering at @xmath96 , and @xmath71 . although the magnitude of @xmath97 is of the same order as the elastic scattering linewidth ( compare to fig .
[ fgam_ex]a ) , its functional dependence on the exciton in - plane momentum is very different . in particular , we note that the maximal linewidth is obtained at a very large momentum ( @xmath98 ) .
this is due to the fact that in order for an exciton with initially small @xmath99 to be ionized , it must scatter on an electron with energy large enough to overcome its binding energy .
this is less likely as the temperature decreases , as is evident from fig .
[ lw_ion]b , where @xmath97 is plotted at @xmath3 .
the insets in figure [ lw_ion ] show the maximum values of @xmath97 vs. @xmath73 . as in the elastic scattering case
, the linear dependence on @xmath73 holds for much larger electron densities at @xmath71 .
at low temperatures and low @xmath73 , charged excitons are formed . in this section
we calculate the trion linewidth due to scattering with free electrons .
three scattering processes are considered : elastic scattering , where the trion remains bound , capturing scattering , where the heavy - hole exciton captures a free electron to form a trion , and dissociating scattering where the trion dissociates into a heavy - hole exciton and an extra free electron .
the latter process is important due to the small binding energy of the trion with respect to the exciton .
we neglect the scattering in which both electrons become unbound . in order to facilitate the matrix elements calculation , we consider a simple two parameter chandrasekhar type @xcite trial function for the trion @xmath100 , \label{phitr}\end{aligned}\ ] ] where the + ( - ) sign applies to the singlet ( triplet ) spin configuration . in eq .
( [ phitr ] ) , @xmath101 are the in - plane coordinates of the two electrons with respect to the hole in - plane coordinate , @xmath102 are the electrons @xmath9 coordinates , and @xmath103 , @xmath104 are the confinement functions in the @xmath9 direction , taken to be the same for all constituents ( eq .  ( [ chi ] ) ) .
the two electron orbitals are given by the normalized functions    @xmath105    and the trion wavefunction normalization factor is given by @xmath106 where @xmath107 the trial function ( [ phitr ] ) was used in @xcite for the case of a two dimensional negative - donor center ( i.e. , @xmath108 limit ) .
in addition , these authors performed the binding energy calculations using a more complicated version of ( [ phitr ] ) which includes a correlation term ( with an additional variational parameter ) .
the variational parameters @xmath109 are calculated by maximizing the trion binding energy for the two spin configurations .
the details of the binding energy calculation are given in appendix c , resulting in @xmath110 , @xmath111(@xmath112 , @xmath113 ) for the singlet ( triplet ) spin configuration , where @xmath114 is calculated with respect to the heavy hole 1s neutral exciton energy .
we note that the triplet configuration is barely bound , in accordance with the experimental observations , thus it can be disregarded . comparing with the much more elaborate treatment of stb _ et al _ @xcite ,
our results are fairly accurate , giving us the confidence to proceed in calculating the trion - electron scattering matrix elements using the wave function ( [ phitr ] ) .
we construct the state of a single trion with an in - plane cm momentum @xmath115 , similarly to that of a single exciton ( [ ex_st ] ) : @xmath116 where we have added the spin index to the electron creation operators , and denoted @xmath117 , and @xmath118 as the projection of a generic spin configuration of two electrons on the singlet spin configuration .
the in - plane fourier transformed trion wave function in eq .
( [ tr_st ] ) is given by @xmath119 where ,    @xmath120^{3/2 } } \\
\phi'_{{\bf k}}&=&\sqrt{\frac{8 \pi \lambda^{\prime 2}}{a } } \frac{1 } { \left[1+(\lambda ' k)^2 \right]^{3/2}}.\end{aligned}\ ] ]    a state comprising a trion and a free electron having @xmath121 , @xmath12 , and @xmath122 will be written as @xmath123 in the following we shall use these wave functions to calculate the scattering matrix elements for the various scattering processes .      applying the coulomb interaction operators @xmath124 , @xmath125 , and @xmath23 to the initial state of the system ( [ initr ] ) , and discarding the self energy contributions , we can write a generic expression for the three scattering matrix elements :    @xmath126    where we put the @xmath127 prefactor to indicate averaging over the initial electron spin states , and used a numeric index for the electron operators , indicating both spatial and spin degrees of freedom .
note that we have omitted the hole operators since they anticommute with the electron operators and do not contribute additional constraints . using the electrons anticommutation relations ,
the fermi vacuum expectation value of the operators reads : @xmath128 where we have symbolically written @xmath129 .
for the 1st electron interaction term we have : @xmath130 the second electron interaction term is found by making the substitutions : @xmath131 , @xmath132 in eqs .
( [ ops ] ) .
similarly , the electron - hole interaction term is found by substituting : @xmath131 in ( [ ops ] ) . applying the delta functions of ( [ delspin ] ) to eq .
( [ vei ] ) , it is evident that the contributions of the @xmath9-dependent parts are identical for all three terms , since we have taken identical confinement functions for the electron and the hole ( see eq .
[ chi ] ) . assuming that the trion remains in the singlet electron spin configuration
, we sum over the spin degrees of freedom , resulting in selection rules for the three electron final spin states , which determine the relevant signs of the various terms . performing the in - plane momentum integrations , and collecting the matrix elements into one direct and two exchange contributions
, we finally find :    @xmath133- \right .
\nonumber \\ & & \left . \left [ g(\lambda \alpha_t q/2 ) g(\lambda ' \alpha_t q/2)+\kappa^2g^2(\tilde{\lambda } \alpha_t q ) \right ] \right\ } \label{trion_dir}\end{aligned}\ ] ]    and    [ trion_exc ] @xmath134 \times \right .
\nonumber \\ & & \left .
\left [ g(\lambda ' |{\bf k}_1-\delta { \bf k}_t+{\bf q}|)- g(\lambda ' |\delta { \bf k}_t- { \bf q}| ) \right ] + [ \lambda \leftrightarrow \lambda']\times [ \lambda \leftrightarrow \lambda ' ] \right\ } \\
v^{(2)}_{\rm exc}({\bf q},\delta{\bf k}_t)&= & -\frac{{\cal n}_{\rm es}}{\pi } \int d^2 k_1 h(\lambda k_1)\left\ { \left [ \lambda \lambda ' \kappa g(\tilde{\lambda}|{\bf k}_1-\alpha_t q| ) g(\lambda ' |{\bf q}-\delta { \bf k}_t| ) + \lambda^2 g(\lambda '    k}_t| ) \right ] \times \right .
\nonumber \\ & & \left .
g(\lambda |{\bf k}_1-\delta { \bf k}_t+\alpha_t{\bf q}| ) + [ \lambda \leftrightarrow \lambda']\times[\lambda \leftrightarrow \lambda ' ] \right\},\end{aligned}\ ] ]    where we have denoted @xmath135 , @xmath136 , and @xmath137 . in eqs.([trion_dir])-([trion_exc ] ) , the prefactor , multiplying the dimensionless quantities , is given by @xmath138 we note that the direct term ( [ trion_dir ] ) diverges at @xmath139 .
this divergence originates from the infinite range of the coulomb potential , which is manifest whenever a scattering event between two charged particles occurs ( note that for the neutral exciton this divergence cancelled out , due to the equal contributions with opposite signs from the electron and hole constituting the exciton ) . in practice
the coulomb potential is screened by the presence of the 2deg .
since we are considering a low density 2deg , its fermi wavevector , @xmath140 is small , thus the semiclassical ( thomas - fermi ) approximation , which assumes that @xmath141 , is inadequate @xcite . in particular , the two - dimensional semiclassical analysis results ( at @xmath142 ) in a density - independent screening wavevector @xmath143 , implying that the presence of a single electron is sufficient to screen the external potential , which is clearly unrealistic .
we adopt instead the lindhard approach which utilizes a perturbative scheme to evaluate the induced charge density in first order of the total potential ( by total we mean the sum of the external potential , i.e. , that of the trion , and the 2deg potential that is induced by the trion s presence ) . the in - plane fourier transformed dielectric function which results from the effect of the screening is given by @xcite @xmath144 where @xmath145 and @xmath146 is the electrons effective bulk bohr radius .
two factors appear in eq .
( [ screene ] ) , which influence the screening .
the first is the screening form factor given by @xmath147 which arises from the finite qw width . in the infinite barrier limit we are considering , the confinement functions , @xmath148 ,
are given by eq .
( [ chi ] ) , and @xmath149 can be calculated analytically , resulting in @xmath150 ^ 2}.\ ] ] this screening form factor is smaller than 1 , thus it reduces the effect of the screening , compared with the strictly two - dimensional case ( @xmath151 ) .
the second factor which appears in eq .
( [ screene ] ) is given by @xmath152 we note that @xmath153 does not depend on the particular details of the qw , but only on the two - dimensional nature of the 2deg motion . substituting the screening dielectric function into the coulomb interaction , we find the screened potential to be @xmath154 thus the use of the screened potential implies replacing @xmath155 in eq .
( [ trion_dir ] ) with @xmath156 in general , the effect of the quasi two - dimensional screening saturates in larger electron densities , in contradistinction to the bulk case , resulting in a complicated dependence of the trion linewidth on @xmath73 .
we note that since there are no divergences in the exchange terms , replacing @xmath157 with @xmath158 in eqs .
( [ trion_exc ] ) amounts to a very small effect ( less than @xmath159 in the relevant parameter range ) .
= 0.6    figure [ matel_fig]a shows the direct and exchange dimensionless matrix elements as a function of transferred momentum @xmath38 , for the case @xmath160 , where the angular dependence of @xmath161 , @xmath162 disappears . the effect of the screening on the total dimensionless matrix element is shown in fig .
[ matel_fig]b , where various electron densities are considered , giving rise to a change in the screening action . as expected , the effect of the screening is noticeable in the region of small momenta of the direct matrix element .
= 0.57    the elastic scattering rate of the trion , @xmath163 , is calculated by writing the energy - momentum conservation equation and finding its solution @xmath164 , similarly to eqs.([econs])-([rate_ex ] ) .
integrating over all final trionic states , we calculate the trion linewidth due to elastic electron scattering , @xmath165 , plotted for two electron densities in fig .
[ trion_lw ] .
the linewidth dependence on @xmath73 is quite complicated .
first , as we are considering @xmath3 , the density range for which the linewidth exhibits a linear dependence on @xmath73 is quite small , similar to the neutral exciton scattering ( cf .
[ fgam_ex]c ) .
second , the screening potential becomes larger as the density decreases , thus for very small densities ( @xmath166 ) decreasing @xmath73 results in a larger trion linewidth .
we note that for such low densities , it is probable that other screening mechanisms become appreciable ( e.g. lattice impurities ) , thus such a behavior will probably not be observed experimentally .
we now consider the case in which a heavy - hole exciton captures a free electron thereby forming a trion .
the initial and final states of the system are given by eqs .
( [ initial ] ) and ( [ tr_st ] ) , respectively .
the calculation of the coulomb interaction matrix elements between these states essentially follows the procedure described in the previous section . in this process
there is no exchange term . instead , the two electrons which constitute the trion form a singlet spin configuration .
defining @xmath167 , the calculation yields : @xmath168 g\left(\lambda \left|{\bf k}_1+\frac{\delta { \bf k}}{\alpha_{\rm x}}\right| \right ) + \right .
\nonumber \\ & & \left .
\frac{\lambda \lambda'}{(\lambda_0+\lambda)^2 } \left [ g\left(\tilde{\lambda } |{\bf k}_1+\delta { \bf k}|\right)-g(\tilde{\lambda } \delta k ) \right ] \times \right .
\\ & & \left .
g\left ( \lambda ' \left|{\bf k}_1+\frac{\delta { \bf k}}{\alpha_{\rm x}}\right| \right ) \right\ } \nonumber\end{aligned}\ ] ] where we have defined the reduced bohr radii : @xmath169 and @xmath170 , and the prefactor : @xmath171    fig .
[ fig_extr]a shows the dependence of the dimensionless matrix element on @xmath42 . proceeding as before , we write conservation of energy : @xmath172 which is satisfied for @xmath173 where @xmath174 .
we note that for the energy conserving solutions , @xmath175 , the interaction matrix element is constant since @xmath176 .
thus , the dependence of the heavy hole exciton linewidth on its initial momentum , shown in fig .
[ fig_extr]b , comes solely from energy momentum conservation and the fermi - dirac distribution of the scattered electron : @xmath177 .
lastly , we consider the case in which the scattered trion dissociates into a neutral exciton and a free electron .
this scattering process is shown schematically in fig.[trex_scat ] .
the initial state of the system is given by eq .
( [ initr ] ) while the final state is that of three free particles , symbolically written as @xmath178 , where @xmath179 is the momentum transferred from the free electron to the trion , and @xmath7 is the neutral exciton cm momentum .
( 160,150 ) ( -5,0)(175,150 )    ( 15,115)t ( 25,118 ) ( 30,115)(1,-1)30    ( 0,135)@xmath180    ( 15,30)e ( 25,32 ) ( 30,35)(1,1)30    ( 10,10)@xmath181    ( 80,90)x ( 90,92 ) ( 95,95)(1,1)30    ( 87,135)@xmath182    ( 80,73)e ( 90,75 ) ( 95,75)(1,0)40    ( 100,83)@xmath183    ( 80,55)e ( 90,58 ) ( 95,55)(1,-1)30    ( 80,10)@xmath184    proceeding as before , we calculate the interaction matrix elements between the initial and final states of the system . to this end
we denote the final spin state of the three electrons as @xmath185 ( see fig .  [ trex_scat ] ) . assuming the interaction is spin independent
, this final state will be one of the two states : @xmath186 , where the state with @xmath187 ( @xmath188 ) corresponds to the initial electron spin @xmath122 being up ( down ) . actually , we need not concern ourselves with the exact form of the final electron spin state , since the projection of the initial spin state ( which is a tensor product of a singlet state with the additional free electron spin ) on a generic final spin state will extract the relevant contributions .
the proper ( anti)symmetrization of this final spin state is performed inherently by the anticommutation relations of the electron operators , as in eq .
( [ delspin ] ) . as in the elastic scattering case
, we find one direct and two exchange contributions . denoting @xmath189 , and @xmath190 as before ,
these terms are expressed as :    @xmath191 \right\ } + \right .
\nonumber \\ & & \left .
\frac{\lambda \lambda'}{(\lambda_0+\lambda')^2 } \left\ { { \renewcommand{\arraystretch}{0.5 } \begin{array}{l } \lambda \leftrightarrow \lambda ' \\ \tilde{\lambda } \rightarrow \tilde{\lambda } ' \end{array } } \right\ } \right ) \label{trex_dir}\end{aligned}\ ] ]    and    [ trex_exc ] @xmath192 \times   \right .
\nonumber \\ & & \left .
\left . g(\lambda'|
{ \bf q}-\delta { \bf k}_{\rm t}| ) + g(\lambda |{\bf k}_1-\delta { \bf k}_{\rm t}+{\bf
q}| ) g\left(\lambda ' \left| \frac{\delta { \bf k}}{\alpha_{\rm x}}+{\bf q } \right| \right ) \right\}+\ { \lambda \leftrightarrow \lambda ' \ } \right ) \\ v^{(2)}_{\rm exc}({\bf q},\delta{\bf k}_{\rm t},\delta{\bf k}_{\rm x})&= & -\frac{1}{2}{\cal n}_{\rm dis } h_s\left(\lambda_0 \left| \frac{\delta { \bf k}_{\rm x}}{\alpha_{\rm x}}-\frac{\delta { \bf k}_{\rm t}}{\alpha_{\rm t}}+{\bf q}\right| \right ) \left ( \frac{\lambda \lambda'}{(\lambda_0+\lambda)^2 } \left\ { g(\lambda '    \frac{\delta { \bf k}_{\rm t}-\beta_{\rm x}\delta { \bf k}_{\rm x}}{\alpha_{\rm x}}-{\bf q } \right| \right ) + \right .
\nonumber \\ & & \left .
g(\tilde{\lambda } \delta k ) \left [ g\left(\lambda ' \frac{\delta k}{\alpha_{\rm x}}\right ) - g(\lambda'| \delta { \bf k}_{\rm t}-{\bf
q}| ) \right ] \right\ } + \frac{\lambda \lambda'}{(\lambda_0+\lambda')^2 } \left\ { { \renewcommand{\arraystretch}{0.5 } \begin{array}{l } \lambda \leftrightarrow \lambda ' \\ \tilde{\lambda } \rightarrow \tilde{\lambda } ' \end{array } } \right\ } \right ) , \label{trex_exc2}\end{aligned}\ ] ]    where @xmath193 , @xmath194 and @xmath42 were defined in the previous subsection and the prefactor multiplying the dimensionless quantities in eqs.([trex_dir])-([trex_exc ] ) , is given by @xmath195    the direct term ( [ trex_dir ] ) and the 2nd exchange term ( [ trex_exc2 ] ) diverge at @xmath196 and at @xmath197 , respectively , implying the need of using a screened coulomb potential , as was done in the trion elastic scattering .
note that the direct term ( [ trex_dir ] ) depends only on @xmath198 , and on the angles @xmath199 and @xmath200 , whereas the exchange terms depend also on the initial electron momentum @xmath201 and on the angle @xmath202 . in practice
, we express everything in terms of @xmath203 , and the angles : @xmath204 , @xmath205 , and @xmath206 , given by :    @xmath207    the calculation of the scattering rate is performed as before , where conservation of energy now reads : @xmath208 equation ( [ cons_trex ] ) is satisfied for @xmath209 \right\}.\end{aligned}\ ] ]    = 0.37    using fermi s golden rule as in eq .
( [ xrate_ion ] ) , we find the scattering rate @xmath210 .
integrating over all final trionic states , we find the trion linewidth due to inelastic scattering as a function of its initial momentum @xmath211 , @xmath212 , shown in fig .
[ lw_trex ] for @xmath213 .
in this section we apply the linewidth calculations for the various scattering processes , presented in the previous sections , in order to calculate the line shapes of neutral excitons and trions .
these line shapes can be directly compared with measured reflection spectra , taken from mtqw structures containing a variable density 2deg . without electrons in the qw and for @xmath214 ,
the two dominant scattering processes are due to acoustic phonons and static disorder .
we note that the radiative lifetime of a free qw exciton is @xmath215 , resulting in a radiative homogeneous broadening of @xmath216 .
the exciton - acoustic phonon interaction leads to a homogeneous broadening whose lorentzian line shape , @xmath217 , is characterized by @xmath218 for the @xmath219 and @xmath220 for the @xmath221 at @xmath71 , using the parameters of gaas / alas qw s @xcite .
the 1s exciton energy , @xmath222 , has a gaussian distribution ( characterized by @xmath223 ) due to the roughness of the qw surface .
convolving these two distributions results in the familiar voigt function @xcite @xmath224 where we take @xmath225 .
when a low density 2deg is present , @xmath226 is convolved with a lorentzian lineshape that is associated with the homogeneous broadening that is due to the electron scattering : @xmath227 where @xmath228 is given by eq .
( [ gam_ex ] ) , with the exciton energy , @xmath229 , corresponding to @xmath99 . a pictorial description of the convolution in eq .  ( [ ls ] )
is given in fig .
[ conv ] , for both elastic and dissociating exciton - electron scattering processes ( the trion - electron scattering exhibits a similar behavior ) .
it is seen that the electron scattering admixes excitons having @xmath230 with the @xmath55 state .
the degree of admixture is determined by the value of the lorentzian that is due to electron scattering , given in the figure by the dashed line .
note that since @xmath231 is a decreasing function of @xmath229 , the lorentzian peak is shifted to higher exciton energies .
( this is observed experimentally as an increased shift with increasing @xmath77 @xcite ) .
= 0.45    in the calculation , the gaussian width , associated with the inhomogeneity , was taken as a fitting parameter , yielding the value @xmath232 for both @xmath219 and @xmath221 ( the same value was taken also for the trion line ) .
the calculated line shapes are identified with the imaginary part of the dielectric function in the qw that is related to the exciton ( trion ) resonance .
use of the kramers - kronig relations yields the real part of the dielectric function .
examples of calculated exciton line shapes together with the associated real part of the dielectric function are depicted in fig  [ line_sh ] for two electron densities .
= 0.4    combining the three exciton line shapes ( heavy - hole and light - hole neutral excitons , and heavy - hole trion ) , by using the appropriate oscillator strengths , the reflection spectra of a qw can be obtained .
for high temperature ( @xmath71 ) , where the trion is dissociated , these calculated spectra and their dependence on @xmath73 were shown to be in good agreement with the experimental data @xcite .    in the low temperature regime ( @xmath233 ) , the vindication of our results with experimental measurements is harder to perform . since the heavy hole exciton and trion lines are very close in energy , and the trion oscillator strength is much lower than that of the exciton , it is practically impossible to discriminate between them by reflection measurements .
pl experiments were performed in mtqw with 2deg , at low temperatures , facilitating the observation of the trion line @xcite , but there the detailed balance between the exciton and the trion line intensities is determined by the 2deg density , and the system dynamics should be addressed by solving , e.g. , rate equations .
= 0.6    in fig .
[ amnon]a the trion linewidth due to elastic electron scattering is shown as a function of @xmath73 for three temperatures .
the combined effect of both elastic and dissociating electron scattering on the trion linewidth in shown in fig.[amnon]b .
qualitatively , the trion line observed in @xcite indeed broadens considerably at much lower @xmath234 than the exciton line , as indicated by our results .
the electron densities stated in @xcite are probably an overestimate , as implied by the fact that the neutral exciton linewidth remains unchanged up to @xmath235 , in clear contradiction to both our calculations and other experimental data ( e.g. , @xcite ) .
it should be noted that calibrating the photoexcitation intensity @xmath236 to @xmath234 in mtqws is not an easy task .
even if one assumes a linear dependence of @xmath73 on @xmath236 ( which is valid for a limited range ) , the slope depends on the temperature , since the holes tunneling times decrease exponentially with @xmath75 , so a different calibration should be performed not only for different experiments , but also for each temperature . finally , the strong enhancement of the dissociating scattering with increasing temperature , evident in fig .
[ amnon]b , implies that the trion line should practically vanish for @xmath237 at very low electron densities .
in this paper we presented a theoretical study of the various exciton - electron scattering processes that take place in a qw which contains a low density 2deg .
we have demonstrated that even for small densities , of the order of @xmath238 in the case of a gaas / alas qw , exciton - electron scattering is much more efficient than both exciton - phonon and exciton - exciton scattering , thus playing a central role in the system dynamics .
the dependence of the exciton linewidth on @xmath73 was shown to be drastically different for the cases of a classical ( @xmath71 ) and a degenerate ( @xmath3 ) 2deg .
in general , the trion - electron scattering poses a much harder problem , one of the major difficulties being the evaluation of the screening effects .
we showed that the screening model has a crucial effect on the trion - electron scattering matrix elements , and hence on the resulting trion linewidth .
this is particularly true for very low densities ( @xmath239 ) , where the 2deg screening is less effective , making other possible mechanisms more dominant . using the calculated scattering rates we showed that dissociating scattering can produce a sizable broadening of the excitonic lines , in particular for the trion , whose binding energy is considerably smaller than that of the neutral exciton .
a simple method to produce excitonic line shapes from the calculated linewidths was devised , allowing our results to be easily compared with experimental reflection spectra .
a reliable source of experimental data is achievable from reflection measurements taken from a mtqw structure , embedded in a microcavity ( mc ) . in the limit of large mc mode energy ,
the exciton linewidths are restored , since their coupling with the mc mode is barely noticeable .
a major advantage of mc experiments is that the trion line is easily resolved by reflection @xcite , allowing a direct comparison with theory .
the calculated interaction matrix elements can be readily used to evaluate polariton - electron scattering rates in mcs with 2deg . in a recent experiment
we have demonstrated polariton final - state stimulation , assisted by these scattering mechanisms , which was hitherto observed using polariton - polariton scattering mechanism @xcite .
incorporation of the various scattering processes which were considered in this paper , in a consistent explanation for these nonlinear effects , will be the topic of a forthcoming publication .
this work was supported by the united states - israel binational science foundation ( bsf ) , jerusalem , israel , by the fund for promotion of research at the technion , and by the technion vpr fund - jewish communities of germany research fund .
g.r . gratefully acknowledges the financial help from the technion .
in this appendix we provide the details of the direct term calculation . separating the coordinates in the qw plane ( @xmath240 ) from the perpendicular coordinate ( @xmath9 ) , and denoting the electron and hole in - plane momenta by @xmath8 and @xmath241 , respectively
, we write the in - plane fourier transform of the exciton wave function @xmath242 $ ] @xmath243 \times \nonumber \\ & & { \rm e}^{-{\rm i}({\bf k}_1 \cdot { \bf r}_{e\| } + { \bf k}_2 \cdot { \bf r}_{h\|})},\end{aligned}\ ] ] where @xmath26 denotes the qw surface area . transforming to center - of - mass ( cm ) and relative coordinates in the qw plane ( @xmath244 , where @xmath245 , @xmath29 , and @xmath246 , @xmath15 are the electron and exciton in - plane effective masses , respectively ) we can decompose the exciton wave function into an envelope function , @xmath247 , and a free motion part related to the in - plane cm coordinate . denoting @xmath7 as the in - plane cm momentum we find @xmath248 } \nonumber \\ & = & \delta({\bf k}_{\rm x}- { \bf k}_1-{\bf k}_2 ) \phi_{{\bf
k}_1 -\alpha { \bf k}_{\rm x } } ( z_e , z_h)\end{aligned}\ ] ]    we use the simplest exciton wave function for the exciton ground state : @xmath249 , where @xmath250 is a normalization factor and @xmath32 is a variational parameter associated with the exciton bohr radius in the qw , which is fixed by maximizing the binding energy of the exciton @xcite .
we note that the use of a wave function separable in @xmath9 and @xmath251 , facilitates considerably the calculation of the scattering matrix elements , although it is strictly justifiable only for narrow well structures . assuming perfect confinement of electrons and holes in the qw , and taking the @xmath9 axis origin in the center of the qw , the confinement functions are : @xmath252 this approximation is quite accurate for almost all practical thicknesses @xmath31 . calculating the binding energy of the 1s exciton in a gaas qw of width @xmath253
, we find @xmath254(@xmath255 ) for the heavy - hole ( light - hole ) exciton .
these results , obtained with our single parameter trial function , are quite close to those obtained in @xcite , which took into account finite qw barriers , and anisotropic masses .    combining the first terms from eqs .
( [ matel ] ) , the direct term reads @xmath256
g(\lambda k_1 ) \times \nonumber \\ & & \cos^2(\pi z_e
/ l ) \cos^2(\pi
z_h / l ) \cos^2(\pi z_c / l),\end{aligned}\ ] ] where @xmath257 was defined in eq .
( [ gh ] ) .
we note that the electron wave functions in the qw plane are absent , as they contribute , together with the in - plane center of mass part of the exciton wave function , a fixed phase factor which is unimportant for our purposes .
the momentum integral is a simple convolution and can be readily evaluated @xmath258
g(\lambda k_1 ) \nonumber \\ & = & \frac{1}{4 \pi^2 } \int d^2 \rho \tilde{g}^2(\rho/\lambda ) \left [ \exp ( { \rm i } \mbox{\boldmath $ \rho$ } \cdot \beta { \bf q } ) - \exp ( { \rm i } \mbox{\boldmath $ \rho$ } \cdot \alpha { \bf q } ) \right ] \nonumber \\ & = & \frac{\pi}{2 \lambda^2 } \left [ g(\lambda \beta q/2 ) - g(\lambda \alpha q/2 ) \right]\end{aligned}\ ] ] where @xmath259 is the exciton spatial wave function in the qw plane . in order to evaluate the @xmath9 integrals
, we change to the coordinates @xmath260 which implies the following change in the integration limits : @xmath261 and performing the integrations results in the direct term , given by eq .
( [ vdir ] )
in this appendix we present an analytical calculation of the angular part of the exchange integral .
this calculation , while quite complicated , is worth the trouble for two reasons .
the first is , obviously , a major reduction in the time consumed by the numerical computation , making it feasible . as a subsidiary benefit ,
the partially analytical calculation removes some convergence problems encountered in the previously used numerical algorithm .    using dimensionless momenta
( e.g. , by writing @xmath38 we mean @xmath262 ) the exchange integral reads    @xmath263^{-3/2}-\left [ 1+({\bf q}- \delta { \bf k})^2 \right]^{-3/2 } \right\ } \times \nonumber \\ & & \left[1+({\bf k}_1- \delta { \bf k}+\alpha { \bf q})^2 \right]^{-3/2 }   \equiv -4 \pi \int dk_1 k_1 h(k_1 ) \left\ { ( i ) + ( ii ) \right\ } ,   \label{geni}\end{aligned}\ ] ]    where @xmath50 is the angle between @xmath179 and @xmath264 and @xmath51 is the angle between @xmath8 and @xmath264 .
we shall work out separately the two parts of eq .
( [ geni ] ) .
the first part is given by @xmath265^{-3/2},\end{aligned}\ ] ] where we have defined @xmath266 transforming @xmath267 and then to @xmath268 we have @xmath269^{-3/2}.\end{aligned}\ ] ] finding the roots of this 4th order polynomial we have @xmath270^{-3/2 } \label{i1}\end{aligned}\ ] ] where @xmath271 we note for later use that @xmath272 obey the inequalities @xmath273 one can reduce the polynomial in eq .
( [ i1 ] ) to a product of two quadratic binomial forms , by making use of the following transformation @xmath274 with @xmath275 , @xmath38 given by @xmath276 ( note that we assume @xmath277 , deferring for later the treatment in this special case ) .
noticing that @xmath275 and @xmath38 always have opposite signs , we choose @xmath278 and divide the integral into @xmath279 where the limit @xmath280 corresponds to @xmath281 and the limit @xmath282 corresponds to @xmath281 with @xmath38 and @xmath275 interchanged .
the integral now reads    @xmath283 ^ 2 } { \left [ ( y^2+s^2)(y^2+t^2 ) \right]^{3/2 } } - \frac{1}{(st)^3 } \int^{-q / p}_{-1 } dy \frac { \left [ y^2(1+p^2)+2y(1+pq)+1+q^2 \right]^2 } { \left [ ( y^2 + 1/s^2)(y^2 + 1/t^2 ) \right]^{3/2 } } \right\}\ ] ] where @xmath284 and @xmath285^{3/2}}.\ ] ] making the transformation @xmath286 in the first integral and @xmath287 in the second , we find @xmath288 ^ 2}{(1-k^2 \sin^2 \psi)^{3/2 } } - \right .
\\ & & \left .
st^3 \int_{\tan^{-1 } ( -s)}^{\tan^{-1 } ( -qs / p ) } d \psi \frac{\left [ ( p / s \sin \psi + q
\cos \psi)^2+(1/s \sin \psi + \cos \psi)^2 \right]^2}{(1-\tilde{k}^2 \sin^2 \psi)^{3/2 } } \right\},\end{aligned}\ ] ] where we have defined @xmath289 changing @xmath290 in the second integral , and using @xmath291 we find that the second integral exactly matches the first except that its upper limit is larger by @xmath292 ; thus we are left with @xmath293 ^ 2}{(1-\tilde{k}^2 \sin^2 \psi)^{3/2}}.\ ] ] this integral can be evaluated using incomplete elliptic integrals ( see p.  201 in @xcite ) .
furthermore , one can exploit the functional relations ( p.  911 in @xcite ) @xmath294 where @xmath295 denote incomplete elliptic integrals of the first and second kind , respectively , and @xmath296 denote complete elliptic integrals of the first and second kind , respectively .
the final result is @xmath297\left [ \frac{1}{t^2}(1+p^2)^2+(1+q^2)^2s^2 \right ] - \right .
\nonumber \\ & & \left .
2(3 + 3p^2q^2 + 4pq+p^2+q^2 ) \left [ 2s^2 { \bf e}(\tilde{k})-(s^2+t^2 ) { \bf k}(\tilde{k } ) \right ] \right\ } \ . \label{ires}\end{aligned}\ ] ]    for imaginary values of @xmath298 , the result is the same as eq .
( [ ires ] ) with the interchanges @xmath299 and @xmath300 , @xmath301 being defined in ( [ k ] ) .
we now return to examine two special cases which can not be treated using ( [ ires ] ) .
the first case corresponds to @xmath302 .
here we have @xmath303 although @xmath304 , @xmath298 is still finite , taking the value @xmath305 thus , enabling the use of ( [ ires ] ) with some modifications . remembering the inequalities ( [ ineq ] ) we note that the sign of @xmath306 is determined by the sign of @xmath307 .
we obtain , for @xmath308 @xmath309^{3/2 } ( d - b)^2 } \times \label{ires1 } \\ & & \!\!\ !
\left\ { \left[1 + \frac{(1+a^2)^2}{(b - a^2)(d - a^2 ) } \right ] \left [ ( 1+d - a^2){\bf e}(\tilde{k})- \right . \right .
\nonumber \\ & & \!\!\ ! \left .
2 { \bf k}(\tilde{k } ) \right ] \!\!- 2(3a^2 + 1 ) \!\ ! \left
[ 2{\bf e}(\tilde{k})-\!\ ! \left ( 1+\frac{b - a^2}{d - a^2 } \right){\bf k}(\tilde{k } ) \right ] \right\}. \nonumber\end{aligned}\ ] ] for @xmath310 the result is the same as eq .
( [ ires1 ] ) with the interchanges @xmath311 and @xmath312 .
the second special case corresponds to @xmath313 . here
the integral ( [ i1 ] ) reduces to @xmath314^{3/2 } } \int_{-\infty}^{\infty } dz \frac{(1+z^2)^2}{(z^2 + 2az+b)^3}.\ ] ] this integral can be easily evaluated ( see p.  81 in @xcite ) and the result is @xmath315^{3/2}(b - a^2)^{5/2}}.\ ] ]    finally , we evaluate the second part of the exchange integral , appearing in eq .
( [ geni ] ) .
we have @xmath316 where @xmath317 where defined in ( [ abcd ] ) . using the transformation @xmath318 , where @xmath319 , we have @xmath320 where @xmath321 .
this integral can be evaluated using the elliptic integral of second kind , resulting in @xmath322
in this appendix we present the calculation of the trion binding energy , based on the trial wave function given in eq.([phitr ] ) .
it is useful to work with the fourier transformed trion wave function ( [ tr_k ] ) , where we have for the two spin configurations , explicitly @xmath323.\end{aligned}\ ] ] the effective mass hamiltonian is given by              and we employed natural units of length and energy , namely , bulk effective bohr radius and rydberg .
the expectation value of the trion energy is given by @xmath330 where @xmath331 was defined in eq .
( [ s ] ) , and the upper ( lower ) signs apply to the singlet ( triplet ) spin configuration .
the various terms appearing in eq .
( [ etr ] ) are calculated in the in - plane fourier space , resulting in simple 2d convolutions , similar to those in appendix a. the results are :    @xmath332= \left ( \frac{\pi}{l } \right)^2(2+\sigma_z ) + \nonumber \\ & & \left[\frac{1}{\lambda^2}(1+\sigma)-\lambda^3\left(\frac{2}{l } \right)^2 \int_0^{\infty}dq q h(\lambda q)g(\lambda q/2 ) \right]+[\lambda \rightarrow \lambda ' ] \\ { \cal i}_2&=&\left(\frac{2}{l}\right)^2 \langle \phi \chi_1 \chi_h |h_0| \phi ' \chi_1 \chi_h \rangle = \frac{\kappa}{2}\left ( \frac{\pi}{l}\right)^2 ( 2+\sigma_z)+ \frac{4(1+\sigma)}{(\lambda+\lambda')^2}-(\lambda+\lambda ' ) \left(\frac{\tilde{\lambda}\kappa}{l}\right)^2 \int_0^{\infty}dq q h(\tilde{\lambda}q ) g(\tilde{\lambda}q )
\\ { \cal j}_1&=&\left(\frac{2}{l}\right)^2 \langle \phi_1 \chi_1 \phi'_2 \chi_2|h_{12}|\phi_1 \chi_1 \phi'_2 \chi_2 \rangle = 2\lambda^3 \left(\frac{2}{l}\right)^2 \int_0^{\infty}dq q h(2\lambda q ) g(\lambda q/2)g(\lambda'q/2 ) \\ { \cal j}_2&=&\left(\frac{2}{l}\right)^2 \langle \phi_1 \chi_1 \phi'_2 \chi_2 |h_{12}|\phi'_1 \chi_1   \phi_2 \chi_2   \rangle = 2\kappa^2\tilde{\lambda}^3 \left(\frac{2}{l}\right)^2 \int_0^{\infty}dq q h(2\tilde{\lambda}q ) g^2(\tilde{\lambda}q).\end{aligned}\ ] ]      in eqs .  ( [ ij ] ) , @xmath333 are the functions defined in eq .
( [ gh ] ) , and we have denoted @xmath334 .
we note that in the 2d limit ( @xmath335 ) , eliminating the kinetic energy term in the @xmath9 direction , eqs .
( [ ij ] ) reduce to    @xmath336 , & \lambda<\lambda ' \\ \frac{\lambda ' \kappa}{(\lambda'-\lambda)^2}\left [ \left(1 + \frac{\lambda^2}{\lambda^{\prime 2}}\right){\bf e}(1-\lambda^{\prime 2}/ \lambda^2)-2 { \bf k}(1-\lambda^{\prime 2}/ \lambda^2)\right ] , & \lambda>\lambda ' \\
\frac{3 \pi}{4 \lambda } , & \lambda=\lambda ' \end{array } \right . \\ { \cal j}_2&=&\frac{3 \pi \kappa}{2 ( \lambda+\lambda')}\end{aligned}\ ] ]    where @xmath337 , @xmath338 are the complete elliptic integrals of the first and second kind , respectively .
these closed form expressions coincide with the integral expressions of sandler and proetto @xcite , for @xmath339 and vanishing magnetic field ( note a misprint in eq .
a6 of this paper ) .    in order to maximize the trion binding energy , we must first write the expectation value of the exciton energy : @xmath340 resulting in the 1s exciton bohr radius and binding energy values stated in appendix a. note that in the 2d limit we have @xmath341 minimizing eq .
( [ ex0 ] ) in the @xmath339 limit , results in @xmath342 , @xmath343 , as expected .
the trion binding energy with respect to the exciton is related to the trion and exciton energies through @xmath344 where @xmath345 is the confined electron energy , in our infinite barrier model .
we note that the kinetic energy terms in the @xmath9 direction cancel out in the expression for the trion binding energy , so they can be removed from the calculation .
maximizing eq .
( [ ebtr ] ) with respect to the two variational parameters @xmath109 , we finally find the results given in section [ trion ] , where we used the gaas bulk values of the effective bohr radius ( @xmath346 ) , and rydberg ( 5.7mev ) .
99 r.  c.  c.  leite , j.  shah , and j.  p.  gordon , phys .
* 23 * , 1332 ( 1969 ) .
l.  schultheis , j.  kuhl , a.  honold and c.  w.  tu , phys .
* 57 * , 1635 ( 1986 ) .
a.  honold , l.  schultheis , j.  kuhl and c.  w.  tu , phys .
b * 40 * , 6442 ( 1989 ) .
v.  capozi , l.  pavesi , and j.  l.  staehli , phys .
b * 47 * , 6340 ( 1993 ) .
k.  matsuda , m.  hirooka , and s.  sunakawa , prog .
theor . phys . * 54 * , 79 ( 1975 ) ; s.  g.  elkomoss and g.  munschy , j. phys .
solids * 38 * , 557 ( 1977 ) ; * 40 * , 431 ( 1979 ) ; * 42 * , 1 ( 1981 ) . y .-
feng and h.  n.  spector , j. phys .
solids * 48 * , 593 ( 1987 ) ; superlattices and microstructures * 3 * , 459 ( 1987 ) ; y .- p .
feng z.  huang , and h.  n.  spector , j. phys .
solids * 50 * , 117 ( 1989 ) ; t.  s.  koh , y .- p .
feng and h.  n.  spector , j. appl .
phys . * 81 * , 2236 ( 1997 ) .
g.  ramon , r.  rapaport , a.  qarry , e.  cohen , a.  mann , arza ron , and l.  n.  pfeiffer , phys .
b * 65 * , 085323 ( 2002 ) .
m.  a.  lampert , phys .
lett . * 1 * , 450 ( 1958 ) . b.  stb and a.  ainane , superlatt . microstruct . * 5 * , 545 ( 1989 ) .
k.  kheng , r.  t.  cox , y.  merle daubign , f.  bassani , k.  saminadayar , and s.  tatarenko , phys .
* 71 * , 1752 ( 1993 ) .
g.  finkelstein , h.  shtrikman , and i.  bar - joseph , phys .
74 * , 976 ( 1995 ) ; a.  j.  shields , m.  pepper , d.  a.  ritchie , m.  y.  simmons , and g.  a.  c.  jones , phys .
b * 51 * , 18049 ( 1995 ) ; arza ron , h.  w.  yoon , m.  d.  sturge , a.  manassen , e.  cohen , and l.  n.  pfeiffer , sol . state comm .
* 97 * , 741 ( 1996 ) .
a.  i.  bobrysheva , m.  v.  grodetskii , and v.  t.  z.  zyukov , j. phys .
c * 16 * , 5723 ( 1983 ) ; a.  thilagam , phys . rev .
b * 55 * , 7804 ( 1997 ) . b.  stb , g.  munschy , l.  stauffer , f.  dujardin , and j.  murat , phys . rev .
b * 56 * , 12454 ( 1997 )
. j.  r.  chapman , n.  f.  johnson , and v.  n.  nicopoulos , phys .
b * 55 * , r10221 ( 1997 ) ; d.  m.  whittaker and a.  j.  shields ,
_ ibid _ * 56 * , 15185 ( 1997 ) ; a.  wjs , j.  j.  quinn , and p.  hawrylak , _ ibid _ * 62 * , 4630 ( 2000 ) . s.  glasberg , g.  finkelstein , h.  shtrikman , and i.  bar - joseph , phys .
b * 59 * , r10425 ( 1999 ) .
y.  yayon , m.  rappaport , v.  umansky , and i.  bar - joseph , phys .
b * 64 * , 081308(r ) ( 2001 ) .
a.  manassen , e.  cohen , arza ron , e.  linder , and l.  n.  pfeiffer , phys .
b * 54 * , 10609 ( 1996 ) .
f.  tassone and y.  yamamoto phys .
b * 59 * , 10830 ( 1999 ) . c.  ciuti , v.  savona , c.  piermarocchi , a.  quattropani and p.  schwendimann , phys .
b * 58 * , 7926 ( 1998 ) .
the angle between @xmath7 and @xmath347 is thus given by @xmath348 , where it is noted that since one can make the transformation @xmath349 , the result of the integration over @xmath350 is indifferent to the sign choice in @xmath351 .
d.  huang , j .-
chyi , and h.  morko , phys .
b * 42 * , 5147 ( 1990 ) .
s.  chandrasekhar , astrophys .
j. * 100 * , 176 ( 1944 ) ; h.  a.  bethe and e.  e.  salpeter , _ quantum mechanics of one- and two- electron atoms _ ( plenum , new york , 1957 ) .
n.  p.  sandler and c.  r.  proetto , phys .
b * 46 * , 7707 ( 1992 )
. g.  bastard , _ wave mechanics applied to semiconductor heterostructures _
( les editions de physique cnrs , paris , 1988 ) .
r.  rapaport , a.  qarry e.  cohen , arza ron , e.  linder , and l.  n.   pfeiffer , in _ proceedings of the 25th international conference on the physics of semiconductors _ , edited by n.  miura and t.  ando ( springer - verlag , 2001 ) , p.  679 .
j.  christen , and d.  bimberg , phys .
b * 42 * , 7213 ( 1990 ) .
v.  savona , l.  c.  andreani , p.  schwendimann and a.  quattropani , solid state comm . * 93 * , 733 ( 1995 ) . m.  grundmann and d.  bimberg , phys .
b * 38 * , 13486 ( 1988 ) . i.  s.  gradshteyn and i.  m.  ryzhik , _ tables of integrals , series , and products _
( academic press , san diego , 1994 ) .
r.  rapaport , r.  harel , e.  cohen , arza ron , e.  linder and l.  n.  pfeiffer , phys .
lett . * 84 * , 1607 ( 2000 ) ; phys .
a * 178 * , 481 ( 2000 ) .
g.  savvidis , j.  j.  baumberg , r.  m.  stevenson , m.  s.  skolnick , d.  m.  whittaker , and j.  s.  roberts , phys .
lett . * 84 * , 1547 ( 2000 ) ; j.  j.  baumberg , p.  g.  savvidis , r.  m.  stevenson , a.  i.  tartakovskii , m.  s.  skolnick , d.  m.  whittaker , and j.  s.  roberts , phys .
b * 62 * r16247 ( 2000 ) ; p.  g.  savvidis , c.  ciuti , j.  j.  baumberg , d.  m.  whittaker , m.  s.  skolnick , and j.  s.  roberts , _ ibid . _ * 64 * 075311 ( 2001 ) .
le si dang , d.  heger , r.  andr , f.  buf , and r.  romestain , phys .
lett . * 81 * , 3920 ( 1998 ) ; p.  senellart , j.  bloch , b.  sermage , and j.  y.  marzin , _ ibid . _ * 62 * , r16263 ( 2000 ) ; r.  huang , f.  tassone , and y.  yamamoto , _ ibid . _
* 61 * , r7854 ( 2000 ) .