Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the observable baryon asymmetry of the universe ( bau ) and the nonzero neutrino masses  @xcite provide strong evidences of physics beyond the standard model ( sm ) . the possibility of unification of the gauge interactions is one of the most intriguing guiding principles for building theories beyond the sm , in this case called the grand unified theories ( gut ) . interestingly , one of the most economical realistic guts may explain both bau and neutrino masses . it is well known that @xmath0 is the minimal group that can provide for grand unification . the first and simplest gut , proposed in @xcite , is based on this group . its scalar sector contains only two representations of @xmath4 and @xmath5 , while each generation of matter fields is embedded in two representations of @xmath6 and * 10*. the @xmath0 symmetry is broken down to the sm one by the vacuum expectation value ( vev ) of the scalar singlet in @xmath5 , and the sm higgs belongs to @xmath4 . in spite of its elegance , this model is ruled out since it can not :    * unify the measured sm gauge couplings ; * explain nonzero neutrino masses ; * unify the sm yukawa couplings of charged leptons and down quarks at high scale in the renormalizable level ; * explain present experimental lower bound on the proton lifetime of order @xmath7  years . the first two of these problems were solved by adding a new fermionic * 24 * representation  @xcite . however a whole set of high dimensional operators was needed to address the third problem . the renormalizable generalization of the model of @xcite with the minimal number of scalar fields contains only one new scalar @xmath8 representation  @xcite , so that the electroweak symmetry is broken by both @xmath9 and @xmath10 vevs , providing for realistic yukawas . ultimately , this nonsupersymmetric renormalizable adjoint @xmath0 theory contains scalars in @xmath4 , @xmath5 and @xmath8 , and matter in @xmath11 , @xmath12 ( @xmath13 ) and @xmath14 . ( see appendix  [ app : fields ] for the details . ) it is well known  @xcite that the neutrino masses may be generated through type i  @xcite and type iii  @xcite seesaw mechanisms in the renormalizable adjoint @xmath0 model . in addition , the bau may be explained  @xcite by generating the lepton asymmetry in the out - of - equilibrium decays of fermion triplets responsible for the type iii seesaw , and converting it to a baryon asymmetry by sphaleron transitions  @xcite in the usual baryogenesis  @xcite via leptogenesis ( lg ) scenario  @xcite . in this paper we investigate a new range of the parameters where neutrino masses get a sizable contribution from the one - loop coloured seesaw besides the tree level type i and iii seesaw mechanisms . however , since the loop term is proportional to a linear combination of type i and iii seesaw yukawas , the lightest neutrino still remains massless . with the new contribution to neutrino mass , the allowed range in the parameter space for successful lg is enlarged . one of the interesting results of this model of the neutrino masses and lg in the framework of adjoint @xmath0 theory is the upper bound on the proton lifetime of order @xmath15  years  @xcite , which can be tested in future . we remark that there is a parameter range in adjoint @xmath0 , which is hard to test at the proton decay experiments  @xcite . however within this range the observable neutrino masses and bau can not be generated .    in the next section we investigate generation of the neutrino masses together with the unification of the couplings in the considered model . in section  [ section : nu_masses_exp ] we discuss parametrization of the yukawa couplings in terms of the neutrino experimental data . in section  [ section : lg ] we give formulae for lg . we discuss the numerical results in section  [ sec : numerics ] and conclude in section  [ sec : summary ] . the interactions relevant for the generation of masses of the adjoint fermions are given by @xmath16 with @xmath17 . the masses of the fermions responsible for the neutrino masses are given by  @xcite @xmath18 where the scale of grand unification is @xmath19 . the regime of successful unification determines the approximate mass relations @xmath20 with @xmath21 .    for the mass spectrum of adjoint su(5 ) , which is allowed by unification , proton decay and lg , see fig . [ fg : spectrum ] with the description in section  [ sec : numerics ] . the yukawa interactions needed to generate the neutrino masses are given by @xmath22 which contains @xmath23 where @xmath24 and the definitions used for the representations and multiplets are given in the appendices  [ app : fields ] and [ app : multiplets ] . the scalar doublets @xmath25 and @xmath26 are related to the sm higgs @xmath27 and the heavy scalar @xmath28 by the orthogonal transformation @xmath29 with  @xcite @xmath30 the vevs of physical fields are @xmath31  gev and @xmath32 . hence the terms relevant for the neutrino masses generation at tree level , shown in fig . [ fig:1:loop : nu : masses ]  @xmath33 and @xmath34 , can be written as @xmath35 and interactions relevant for the one - loop contribution to the neutrino masses , shown in fig . [ fig:1:loop : nu : masses ]  @xmath36 , can be written as @xmath37 where in and we have defined @xmath38 notice that in eqs . ( [ h3 ] ) and ( [ h0 ] ) we corrected the result of  @xcite . the contribution of @xmath39 to the @xmath40 parameter is zero due to mass degeneracy of @xmath41 and @xmath42 , because the contributions to neutral and charged goldstones @xcite cancel each other .    only electrically neutral components of @xmath43 enter the mechanisms illustrated in fig . [ fig:1:loop : nu : masses ] . the quadratic coupling @xmath44 arises from several terms in the adjoint @xmath0 lagrangian . here we take it to be arbitrary , bounded only by perturbativity . the resulting neutrino mass matrix is @xmath45   \label{eq : nu_masses}\end{aligned}\ ] ] with @xmath46 the first term in eq . ( [ eq : nu_masses ] ) is generated in type iii , the second term in i seesaw , as in  @xcite . the third term comes from one - loop coloured seesaw  @xcite in the limit of small mass splitting @xmath47 , where @xmath48 is the mass of @xmath49 ( @xmath50 )  @xcite .    in @xcite the neutrino masses and lg were studied for @xmath51 and hence the one - loop contribution to neutrino masses was neglected . in this paper we are interested in the case of @xmath52 , in which the fermionic octet is relevant to the generation of neutrino masses .  show that the number of light massive neutrinos in the considered model is two , which is equal to the maximal number of nonzero independent vectors in the plane defined by the vectors @xmath53 and @xmath54 . and @xmath55 for the couplings are not allowed , because in this case there would be only one massive neutrino . ] however none of @xmath1 , @xmath2 and @xmath3 is decoupled from the generation of the neutrino masses due to its heavy mass or small couplings , as in the two right - handed neutrino models  @xcite . on the other hand , some of @xmath43 may not enter particular physical processes , e.g. , @xmath3 decouples from lg , see section  [ section : lg ] . such decoupling results in new physics effects , which are dependent on the yukawa couplings @xmath56 , in particular , the allowed parameter range of successful lg is enlarged , see section  [ sec : numerics ] . we remark that similar relation among the neutrino masses and lg may be derived in theory with two heavy majorana neutrinos and one light sterile neutrino , which has specific mass and mixing terms with other neutrinos parity and contribute to the neutrino masses due to the interaction with inert scalar doublet , see  @xcite . ] @xcite .    using the linear relation @xmath57 with @xmath58 eq . ( [ eq : nu_masses ] ) can be rewritten as @xmath59 where @xmath60 and @xmath61 with @xmath62 in the limit @xmath63 matrix @xmath64 goes to @xmath65 with @xmath66 .    the neutrino mass matrix in eq .  ( [ eq : nu_masses2 ] ) can be rewritten as @xmath67 by using the orthogonal transformation @xmath68 where @xmath69 is a real orthogonal matrix and @xmath70 , where @xmath71 are the eigenvalues of @xmath72 . the three neutrino mass spectra allowed by the oscillation data are    * normal hierarchical ( nh ) @xmath73 * inverted hierarchical ( ih ) @xmath74 * degenerate , @xmath75  ev ( @xmath76 ) ;    where @xmath77  ev@xmath78 and @xmath79  ev@xmath78 are the mass - squared differences of solar and atmospheric neutrino oscillations  @xcite . we use the standard casas - ibarra  @xcite parametrization of the yukawa couplings @xmath80 as @xmath81 where @xmath82 ; the pmns lepton mixing matrix @xmath83 can be written as @xmath84 where @xmath85 , @xmath86 , with @xmath87 , @xmath88 and @xmath89  @xcite , while the dirac @xmath90 and majorana @xmath91 @xmath92 violation phases are not fixed by the present experiments . finally , the @xmath93 complex matrix can be written as  @xcite @xmath94 in the normal and inverted hierarchy , respectively ; @xmath95 . we define the @xmath92 asymmetry as @xmath96}.\ ] ] in the considered model it is given by  @xcite @xmath97\ , f(m_{\rho_0}^2/m_{\rho_3}^2),\end{aligned}\ ] ] where @xmath98 = -\frac{1}{2\sqrt{x } } + \mathcal{o}(x^{-3/2 } ) \quad \text{for } \quad x\gg1\end{aligned}\ ] ] since the only non - vanishing contribution comes from the vertex correction  @xcite . we assume that the couplings of @xmath99 and @xmath100 interactions are small comparing to @xmath44 , so that two - loop effects in the @xmath92 asymmetry can be neglected .    using the total decay rate @xmath101 where @xmath102 ( signs of the components of @xmath2 ) , the decay parameter @xmath103 , given by the ratio of the decay width of a single @xmath2 component to the expansion rate when @xmath2 starts to become non - relativistic at @xmath104 , can be written as @xmath105 and rewritten as @xmath106 in terms of the rescaled decay rate  @xcite @xmath107 and the rescaled hubble expansion rate ( equilibrium @xmath2 mass ) @xmath108 for nh ( ih ) the strong washout regime requires @xmath109    using eqs . , , and , we have the rescaled decay rate @xmath110,\end{aligned}\ ] ] and the total @xmath92 asymmetry @xmath111   \sqrt{m_{1}m_{2}}\ , { \rm im}(s^*c)\right\ }   \nonumber\\\end{aligned}\ ] ] with @xmath112 , @xmath113 , and @xmath114 , @xmath115 for nh ( ih ) neutrinos .    in the limit @xmath63 we get the standard expressions @xmath116 where in eq . we corrected the result given in  @xcite by the factor of 1/2 . boltzmann equations in the unflavoured regime ( @xmath117  gev ) can be written as ( for more details see  @xcite and refs . therein ) @xmath118 where @xmath119 , and @xmath120 ( @xmath121 ) is the number density of @xmath122 calculated in a co - moving volume containing one @xmath2 ( all of its three components ) in ultrarelativistic thermal equilibrium : @xmath123 . initially , @xmath124 . the abundance at equilibrium is given by @xmath125 where @xmath126 is a modified bessel function . corresponds to the ` besselk[n , z ] ` function in mathematica . ] the decay factor is @xmath127 where @xmath27 is the hubble rate , and @xmath128 is the thermally averaged dilation factor . the sum @xmath129 , where @xmath130 is the contribution from higgs - mediated scatterings , is given by @xmath131,\ ] ] where @xmath132 is higgs mass and @xmath133    for the gauge scattering of @xmath2  @xcite we use the fit  @xcite @xmath134 where @xmath135  gev is the planck mass . the washout can be written as @xmath136 where @xmath137.\ ] ] the contribution of the non - resonant @xmath138 processes to the washout is @xmath139 where @xmath140  ev for nh ( ih ) neutrinos . after solving the boltzmann equations - , we obtain @xmath141 ( in our calculation we use final @xmath142 , where the fit in eq .   still applicable  @xcite ) , included in the final baryon asymmetry @xmath143 where 3 is the number of @xmath2 components , the dilution factor 86/2387 is calculated assuming standard photon production from the onset of lg till recombination  @xcite , and 12/37 is the fraction of @xmath144 asymmetry converted into a @xmath145 asymmetry assuming that sphalerons remain in equilibrium until after the electroweak phase transition  @xcite . the result in eq .   should be compared with the allowed values @xmath146 which come from the nucleosynthesis predictions and observed abundances of light elements  @xcite . gev . dashed line shows the unflavoured and dot - dashed line the flavoured lg bound.,scaledwidth=46.0% ]    in the considered model the unification scale and the masses of @xmath147 and @xmath148 are constrained by the proton decay : @xmath149  gev , @xmath150  gev and @xmath151  gev  @xcite . the parameter space allowed by the unification at one - loop level for @xmath152  gev is shown in fig . [ fg : unification ] , where the dashed line shows the bound from eq .   and the unflavoured lg constraint of @xmath117  gev , and the dot - dashed line shows the flavoured lg bound of @xmath153  gev .        the mass spectrum of fermions and scalars in renormalizable adjoint @xmath0 is shown in fig .  [ fg : spectrum ] , where the lower light gray belt shows the allowed @xmath148 mass range from unification and proton decay constraints , and the upper light gray belt shows the allowed range for @xmath43 masses from lg . the dark gray belt shows the allowed range for the scale of unification . all the sm particle masses lie near or below @xmath27 mass @xmath132 , while the rest of new particle masses in the theory can be set at the gut scale , besides the mass of @xmath154 field , which favours unification and unflavoured ( flavoured ) lg within the wide region @xmath155  gev , where the lower bound comes from the collider searches . the black lines for @xmath148 and @xmath43 in fig . [ fg : spectrum ] show characteristic values of their masses .     vs. @xmath103 in the case of nh ( ih ) . see text for details.,title="fig:",scaledwidth=48.0% ]   vs. @xmath103 in the case of nh ( ih ) . see text for details.,title="fig:",scaledwidth=48.0% ]    fig . [ fg : be ] shows the allowed ranges for @xmath156 versus @xmath103 for the particular parameter sets of @xmath157 , @xmath158  gev , and @xmath159 , @xmath160 and @xmath51 ( which reproduces the standard case in eqs .   ) , where the lighter region corresponds to larger @xmath161 . the results are weakly sensitive to the higgs mass within its expected range , and to @xmath162 due to only logarithmic dependence in eq .  . we use @xmath163  gev and @xmath164  gev . the upper bound for @xmath2 mass of @xmath165  gev is derived using eq .   and fig . [ fg : unification ] . to get the lower bounds on @xmath156 for every given value of @xmath103 , we maximize the solution   of the boltzmann equations   with the @xmath92 asymmetry from eq .   over the complex angle @xmath166 , and find the minimal value of @xmath156 required to produce the observed baryon asymmetry in eq . ( its lower bound ) . the resulting curves in fig . [ fg : be ] are smoothened and given to few percent accuracy . the allowed region for lg in the case of vanishing loop contribution to the neutrino masses is shown by the darkest regions , which agrees with the result of  @xcite . the lighter areas show the new allowed regions for lg for the larger values of @xmath161 up to 1 , when the loop contribution to the neutrino masses is significant . the dashed lines show the lower bounds on @xmath156 for the parameter set @xmath157 , @xmath167  gev and @xmath159 . the vertical lines show the lower bounds on @xmath103 in the strong washout regime . the dot - dashed and dotted curves show limits @xmath168 and @xmath169 , which are derived from eq .   using @xmath170 and 4 , respectively . below these curves the couplings @xmath171 remain well perturbative . we notice that the considered unflavored regime of lg is applicable for @xmath172  gev . below this level , shown in fig . [ fg : be ] ( left ) , the flavor effects should be taken into account . we remark that the maximal allowed value of @xmath173 , which is related to the mass of @xmath148 , can be tested at the future proton decay experiments  @xcite . we have presented a detailed study of the neutrino masses and the baryogenesis via lg in adjoint @xmath0 , which is a well motivated model since testable by the proton decay experiments . @xmath144 asymmetry is generated by the decays of @xmath174 triplet @xmath2 , with singlet @xmath1 running in the loop , while the neutrino masses are generated at the tree level by @xmath2 and @xmath1 exchanges and in one - loop coloured seesaw with propagating @xmath175 octets @xmath3 and @xmath148 . for the values of the octet scalar @xmath148 to sm higgs @xmath27 coupling @xmath44 of order 1 , the one - loop coloured seesaw is on the same order as the type i and iii seesaw . this new contribution to the neutrino masses relaxes the yukawa couplings for @xmath2 and @xmath1 , and makes the allowed range for unflavored lg significant not only for the normal hierarchy of the neutrino masses , but also for the inverted one , in contrast to the case of small @xmath44 . the lightest neutrino , however , remains massless , because the coloured seesaw loop diagram is proportional to the yukawa couplings of the type i and iii seesaw . for @xmath176 the @xmath92 asymmetry explicitly depends on several more parameters of adjoint @xmath0 , in comparing with the case of @xmath51 , such as : @xmath44 , the vev @xmath177 of the representation @xmath178 , and the ratio @xmath179 of the masses of adjoint fermions @xmath180 and @xmath181 . this increases the predictive power and testability of the model . we remark that some of the new particles , e.g. , @xmath154 and @xmath148 , may be light enough in the considered scenario to be tested at the lhc and the next generation of colliders  @xcite . under @xmath182 , the scalar fields are decomposed as @xmath183 and the matter fields are decomposed as @xmath184 where @xmath13 is the generation index , @xmath185 is the sm lepton doublet , @xmath186 are scalar @xmath174 doublets , and we denote @xmath187 , @xmath188 , @xmath189 and @xmath190 , where @xmath191 . the vevs are @xmath192_j^i , \qquad \langle(45_h)\rangle_j^{in}=0,\end{aligned}\ ] ] where @xmath193 , @xmath194 , and relation @xmath195  gev  @xcite provides proper masses of @xmath174 gauge bosons . more details on the representations are given in the appendices a of @xcite . the scalar potential was studied in @xcite , and the yukawa couplings were analyzed in @xcite . for the @xmath196 doublets and triplets we use the explicit forms @xmath197 respectively , where @xmath198 , @xmath199 , and @xmath200 are the pauli matrices . the adjoint representation @xmath201 ( @xmath193 ; @xmath202 ) can be written as @xmath203 where @xmath204 are the gell - mann matrices . * 45*@xmath205 is defined by the conditions : @xmath206 and @xmath207 . we omit index @xmath27 of scalar representations in the following . the explicit decomposition into sm multiplets can be deduced by considering the multiplet as the tensor product @xmath208 . below we will write the colour indices with letters from the beginning of the alphabet @xmath209 and the weak indices as @xmath210 . we have @xmath211_{dc } + \epsilon^{abd } [ s_{(\bar{3},1)}]_{dc } , \\    { 45}^{ab}_{r } & = \epsilon^{abc } [ s_{(\bar{3},2)}]_{cr } , \\    { 45}^{ar}_{b } & = \frac{1}{\sqrt{2 } } s_{8}^{ar } [ \lambda^{a}]^{a}_{b } + h_{2}^{r } \delta^{a}_{b } , \\    { 45}^{ar}_{s } & = \frac{1}{\sqrt{2 } } s_{(3,3)}^{ai } [ \sigma^{i}]^{r}_{s } + s_{(3,1)}^{a } \delta^{r}_{s } , \\    { 45}^{rt}_{s } & = -\frac{3}{2 } h_{2}^{t } \delta^{r}_{s},\end{aligned}\ ] ] where @xmath212_{ab } = -[s_{(\bar{3},1)}]_{ba}$ ] . we thank riccardo barbieri , anatoly borisov , konstantin stepanyantz , ilja dorsner , enrico bertuzzo and david straub for useful suggestions and comments . this work was supported by the esf mobilitas 2 grant mjd140 , mobilitas 3 grant mtt8 , esf grants 8090 and 8943 , and sf0690030s09 ( kk ) ; and by the eu itn `` unification in the lhc era '' , contract pitn - ga-2009 - 237920 ( unilhc ) and by miur under contract 2006022501 ( dz ) . p.  minkowski , _ phys.lett._*b67 * , 421 ( 1977 ) ; t.  yanagida , in _ proceedings of the workshop on the unified theory and the baryon number in the universe _ , eds . o.  sawada and a.  sugamoto ( tsukuba , japan , 1979 ) 95 ; m.  gell - mann , p.  ramond , and r.  slansky , in supergravity , eds . d.  z. freedman and p.van nieuwenhuizen ( north holland , amsterdam , 1979 ) 315 ; r.  n.  mohapatra and g.  senjanovic , _ phys.rev.lett._*44 * , 912 ( 1980 ) . r.  foot , _ et al . _ , z.phys.*c44 * , 441 ( 1989 ) ; e.  ma , phys.rev.lett.*81 * , 1171 - 1174 ( 1998 ) . s. blanchet and p. fileviez perez , jcap*0808 * , 037 ( 2008 ) [ 0807.3740 [ hep - ph ] ] ; mod.phys.lett.*a24 * , 1399 - 1409 ( 2009 ) [ 0810.1301 [ hep - ph ] ] . r.  barbieri , l.  j.  hall and v.  s.  rychkov , phys.rev.*d74 * , 015007 ( 2006 ) [ arxiv : hep - ph/0603188 ] . p. fileviez perez and m. b. wise , phys.rev.*d80 * , 053006 ( 2009 ) [ 0906.2950 [ hep - ph ] ] ; p. fileviez perez , aip conf.proc.*1222 * , 3 - 9 ( 2010 ) [ 0909.2698 [ hep - ph ] ] . m.  a.  luty , phys.rev.*d45 * , 455 ( 1992 ) ; l.  covi , e.  roulet and f.  vissani , phys.lett.*384 * , 169 ( 1996 ) [ peh - ph/9605319 ] . s. davidson , e. nardi and y.  nir , phys.rept.*466 * , 105 - 177 ( 2008 ) [ 0802.2962 [ hep - ph ] ] .
Abstract:

Answer: investigate baryogenesis via leptogenesis and generation of neutrino masses and mixings through the type i plus type iii seesaw plus an one - loop mechanism in the context of renormalizable adjoint @xmath0 theory . one light neutrino remains massless , because the contributions from three heavy majorana fermions @xmath1 , @xmath2 and @xmath3 to the neutrino mass matrix are not linearly independent . however none of these heavy fermions is decoupled from the generation of neutrino masses . this opens a new range in parameter space for successful leptogenesis , in particular , allows for inverted hierarchy of the neutrino masses .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: modern statistical physics encompasses a wide range of scientific problems from elementary particles to the evolution of the universe . its methods facilitate the determination of the macroscopic characteristics of complex nonequilibrium systems without having to exactly analyze the highly complicated individual microscopic dynamics of each constituent part . in many cases the system of interest can be separated into two distinct subsystems . one of them , which we will term the  _ dynamical subsystem _  , has only few degrees of freedom , whereas the other , the  _ bath _  , has , in an ideal case , the infinite number of them . the bath is a source of random perturbations ( fluctuations ) for the dynamical subsystem , also functioning as an absorber of dissipative energy . this separation is at the core of the theory of open quantum systems ( toqs ) based on papers of schwinger @xcite , developed in refs . @xcite and successfully applied to radiation damping @xcite and some problems in solid state physics @xcite .    in the present paper we demonstrate the applicability of toqs to the problem of spin dynamics in semiconductors , which has been extensively studied over the last decade both experimentally @xcite and theoretically @xcite in conjunction with numerous proposals for a variety of spin - based devices @xcite . the most promising materials for device purposes , the iii - v and ii - vi compounds , have been shown to have spin relaxation rates dominated by the dyakonov - perel ( dp ) mechanism @xcite at moderate temperatures and low hole concentrations . in contrast to earlier treatments , we develop a fully microscopic stochastic theory of spin dynamics in the presence of an external magnetic field , taking account of pertinent scattering mechanisms . our analysis is based on an innovative implementation of the two step relaxation process corresponding to the relaxation time hierarchy involved in ( a ) electron thermalization , and ( b ) spin relaxation . in the first stage of solution , we determine the relaxation rates and fluctuation characteristics of electron orbital motion due to coupling to optical and acoustic phonons and random impurities as a bath . spin relaxation dynamics ( the slowest process in the system ) can be neglected in stage ( a ) . the second stage , ( b ) , proceeds with the analysis of the spin relaxation process due to spin - orbit interaction , wherein the orbital degrees of freedom are considered as an  effective heat bath  , having the characteristics determined in stage ( a ) . the orbital motion can be considered as an intermediary transmitting fluctuations from the phonon / impurity bath to electron spin , and , also , transferring the dissipated energy flow in the opposite direction . a set of bloch equations with two distinct relaxation times ( longitudinal relaxation time , @xmath0 , responsible for spin magnetic moment relaxation , and transverse relaxation time , @xmath1 , responsible for decoherence ) is derived in this second stage . in both stages of our analysis we employ the toqs method . the structure of the paper is as follows . in section ii we present a brief outline of toqs and the principal relations used in subsequent sections . the microscopic formulation of the spin relaxation problem is given in section iii . in section iv we analyze the orbital electron dynamics in the presence of an external magnetic field with coupling to optic and acoustic phonons and random impurities . the spin - orbit interaction is examined in section v , where we derive the bloch equations on a microscopic basis . a summary of this work is given in section vi . the full hamiltonian of the system , separated into a  dynamical subsystem  and a  bath  , can be written as @xmath2 where @xmath3 is the hamiltonian of the dynamical subsystem , @xmath4 is the hamiltonian of the bath and @xmath5 describes their interaction . we can write the last term in a product form @xmath6 where @xmath7 is , in the general case , a nonlinear function of the dynamical subsystem variables , @xmath8 is a function of the bath variables , and the summation index @xmath9 can refer to either coordinate projections or a mode index . one can see from @xmath5 that both the dynamical subsystem and the bath influence each other . moreover , @xmath10 plays a role of a generalized force conjugate to the generalized bath coordinate , @xmath11 and vice versa ( figure1 ) .    in the case of gaussian statistics of the unperturbed bath variables , or when the coupling of the dynamical subsystem to the bath is weak , the full heisenberg operator of the bath is given by @xcite @xmath12 where @xmath13_{-}\rangle \eta ( t - t_{1 } )   \label{fi}\ ] ] is a linear response function , or the retarded green s function of the unperturbed bath variables . @xmath14 is the heaviside unit step function . another important function is the correlation function of the unperturbed bath variables @xmath15_{+}\rangle = m_{\alpha \beta } ( \tau ) . \label{em}\ ] ] here and below we use the notation @xmath16 _ { -}$ ] for a commutator , and @xmath16 _ { + } $ ] for an anticommutator . the linear response function of eq.([fi ] ) , and the correlation function of eq.([em ] ) , are related by means of the fluctuation - dissipation theorem @xmath17 where @xmath18 and @xmath19 to derive the fluctuation - dissipation theorem we employed the relation @xmath20_{+}\rangle = \langle \lbrack q_{\beta } ^{0}(0),q_{\alpha } ^{0}(\omega ) ] _ { -}\rangle coth{% \frac{\hbar \omega } { 2k_{b}t } , }   \label{coth}\ ] ] based on the gibbs distribution and , therefore , the fluctuation - dissipation theorem , in the present form , is valid only for systems in equilibrium . the heisenberg equation of motion for an arbitrary operator of the dynamical subsystem , @xmath21,is given by @xmath22_{-}-{\frac{1}{i\hbar } } % [ a(t),f_{\alpha } ( t)]_{-}q_{\alpha } ^{h}(t),\ ] ] where @xmath23 is a heisenberg bath operator including the influence of the dynamical subsystem . substituting the expression of eq . ( [ response ] ) for this operator , and taking account of the fact that the operator of the dynamical subsystem is commutative with @xmath23 at the same moment of time ( only ) , we obtain @xmath22_{-}-{\frac{1}{2}}[q_{\alpha } ^{0}(t),y_{\alpha } ( t)]_{+}-{\frac{1}{2}}\int dt_{1}\varphi _ { \alpha \beta } ( t , t_{1})[y_{\alpha } ( t),f_{\beta } ( t_{1})]_{+},\ ] ] where @xmath24_{-}.\ ] ] considering @xmath25 as a function of the unperturbed bath variables and employing the quantum analog of the furutsu - novikov theorem @xcite , we find that @xmath26_{+}\rangle = \int dt_{1}m_{\alpha \beta } ( t , t_{1})\langle { \frac{\delta y_{\alpha } ( t)}{\delta q_{\beta } ^{0}(t_{1})}}\rangle , \ ] ] and eliminate the appearance of the bath variables using functional derivatives as described by efremov and smirnov @xcite @xmath27_{-}\rangle \eta ( t - t_{1}).\ ] ] introducing the fluctuation source , @xmath28 , @xmath29_{+}+\int_{-\infty } ^{t}dt_{1}m_{\alpha \beta } ( t , t_{1}){\frac{i}{\hbar } } [ y_{\alpha } ( t),f_{\beta } ( t_{1})]_{-},\ ] ] with zero mean value @xmath30 , we find that the equation for the arbitrary operator @xmath21 of the dynamical subsystem has the form @xmath31_{-}-{\frac{1}{2}}\int_{-\infty } ^{\infty } dt_{1}\varphi _ { \alpha \beta } ( t , t_{1})[y_{\alpha } ( t),f_{\beta } ( t_{1})]_{+}-   \nonumber \\ & & -\int_{-\infty } ^{t}dt_{1}m_{\alpha \beta } ( t , t_{1}){\frac{i}{\hbar } } % [ y_{\alpha } ( t),f_{\beta } ( t_{1})]_{-}+\xi ( t).\end{aligned}\ ] ] in the case of weak coupling between the dynamical subsystem and the bath , the correlator of the fluctuation force is given by @xmath32 _ { + } \right\rangle = m_{\alpha \beta } ( t , t_{1})\left\langle { \frac{1}{2}}[y_{\alpha } ( t),y_{\beta } ( t_{1})]_{+}\right\rangle + r_{\alpha \beta } ( t , t_{1})\left\langle { \frac{1}{2}}[y_{\alpha } ( t),y_{\beta } ( t_{1})]_{-}\right\rangle , \ ] ] where @xmath33_{-}\rangle = \frac{\hbar } { 2i}\left ( \varphi _ { \alpha \beta } ( t , t_{1})-\varphi _ { \beta \alpha } ( t_{1},t)\right ) .\ ] ] we start from the model hamiltonian @xmath34 describing an electron with spin in the presence of a magnetic field directed along the @xmath35-axis , where the magnetic field and its vector potential are given by @xmath36 the first term in the hamiltonian ( [ ham ] ) is responsible for kinetic electron orbital motion , @xmath37 where the velocity component operators of the electron in a magnetic field can be written as @xmath38 _ { - } & = & -\frac{i\hbar \omega _ { c}}{m}.   \nonumber\end{aligned}\ ] ] the next term , describing spin motion in the presence of a magnetic field ( zeeman term ) , is given by @xmath39 where @xmath40 is the crystal @xmath40-factor ( in particular , @xmath41 for gaas ) , @xmath42 is the bohr magneton ( @xmath43 is the mass of a free electron , the effective mass , @xmath44 , is @xmath45 for gaas ) , @xmath46 is the frequency of the spin precession induced by the magnetic field , and @xmath47 is the standard set of pauli matrices . the free phonon hamiltonian is given by @xmath48 where @xmath49 is the phonon energy and @xmath50 and @xmath51 are the creation and annihilation operators of a phonon in the @xmath52-mode , respectively . the interaction of an electron with random impurities is described by the potential @xmath53 where @xmath54 are the spatial fourier components of the impurity potential , and @xmath55 denotes the volume of the crystal . the electron - phonon interaction is given by @xmath56 where @xmath57 presents the strength of the electron - phonon coupling . as mentioned above , the principal mechanism of spin - orbit interaction is the dyakonov - perel term , @xmath58 . in semiconductor crystals lacking inversion symmetry , the effective mass hamiltonian includes terms cubic in electron quasimomentum @xcite . the presence of these terms induces a random magnetic field @xmath59 given by @xmath60 where @xmath61 @xmath62 is the energy gap , and @xmath9 is a coefficient representing the strength of the spin - orbit coupling ( @xmath63 in gaas ) . the random magnetic field induced by the orbital motion interacts with electron spin through the interaction hamiltonian @xmath64    it is evident from the hamiltonian ( [ ham ] ) that there is no _ direct _ coupling of electron spin degrees of freedom to phonons and impurities . therefore , to analyze spin relaxation dynamics , we consider the orbital degrees of freedom as an intermediary that ( a ) transfers fluctuations from the phonon / impurity bath to electron spin , and ( b ) transfers energy from spin degrees of freedom to heat bath . this two - step process is illustrated in figure 2 . correspondingly , we address the analysis in two stages , and in the first stage we determine the statistical characteristics of orbital motion dictated by the phonon / impurity bath , neglecting the presence of the electron spin . the pertinent first - stage hamiltonian , @xmath65 , is given by @xmath66 where the first - stage interaction hamiltonian , @xmath67 , has the form @xmath68 comparing this expression to eq . ( [ hint ] ) , we take the summation index @xmath69 of eq . ( [ hint ] ) , as the mode number , @xmath52 , and the function @xmath70 plays the role of  @xmath7  of eq . ( [ hint ] ) , whereas the bath variable ,  @xmath8 , is @xmath71 .    in the second stage of analysis we employ orbital motion as an  effective heat bath  having the characteristics determined in the first stage . the second - stage hamiltonian , @xmath72 , can be written as @xmath73 where @xmath74 and @xmath75  may be taken as the remainder of the hamiltonian ( [ ham ] ) . to compare @xmath58 to eq . ( [ hint ] ) , we rewrite @xmath58 as    @xmath76    where the  effective heat bath  variables are given by @xmath77 in this , there is summation over coordinate projections and the  dynamic subsystem  variables are the electron spin projections . in this step we employ toqs to determine the statistical characteristics of electron orbital motion described by the operator equations @xcite @xmath78 and @xmath79 where @xmath80 are electron velocity operator components , and @xmath81 is the cyclotron frequency . the electron - bath interaction determines the relaxation rates , @xmath82,@xmath83 , the frequency shift , @xmath84 , and the fluctuation sources , @xmath85 , involved in eq . ( [ velocity ] ) , @xcite . the fourier transforms of the velocity correlation functions are given by @xmath86 _ { + } \right\rangle = \left\langle \frac{1}{2}\left [ v_{y}(\omega ) ; v_{y}\right ] _ { + } \right\rangle = \frac{k_{\bot } ( \omega ) } { 2}\left ( \frac{1}{\left ( \omega -\omega _ { c}\right ) ^{2}+\gamma _ { 0}^{2}}+\frac{1}{\left ( \omega + \omega _ { c}\right ) ^{2}+\gamma _ { 0}^{2}}\right ) ,   \label{xx}\ ] ] @xmath87 _ { + } \right\rangle = -\left\langle \frac{1}{2}\left [ v_{y}(\omega ) ; v_{x}\right ] _ { + } \right\rangle = \frac{k_{\bot } ( \omega ) } { 2i}\left ( \frac{1}{\left ( \omega -\omega _ { c}\right ) ^{2}+\gamma _ { 0}^{2}}-\frac{1}{\left ( \omega + \omega _ { c}\right ) ^{2}+\gamma _ { 0}^{2}}\right ) ,   \label{xy}\ ] ] and @xmath88 _ { + } \right\rangle = \frac{k_{z}(\omega ) } { \omega ^{2}+\gamma _ { z}^{2 } } ,   \label{zz}\ ] ] where @xmath89 _ { + } \right\rangle = \int d\left ( t - t_{1}\right ) e^{i\omega \left ( t - t_{1}\right ) } \left\langle \frac{1}{2}\left [ \xi _ { y}(t),\xi _ { y}(t_{1})% \right ] _ { + } \right\rangle , \ ] ] and @xmath90 _ { + } \right\rangle .\ ] ]    to be specific we take account of the contributions of polar optical phonons , deformational acoustic phonons and random impurities to electron orbital dynamics . in a first approximation in their ( weak ) coupling strengths to orbital motion , these contributions can be treated as mutually independent and the overall values of the damping rates and the time - fourier transforms of the correlation functions are given by ( superscripts @xmath91 and @xmath92  denote the contributions of optical and acoustic phonons , respectively , and @xmath93 denotes that of random impurities ) : @xmath94 the microscopic expressions of all these quantities are determined using the method of ref . @xcite , and are given below :    the response and correlation functions for polar optical phonons are @xmath95 resulting in @xmath96 and @xmath97 where @xmath98 + \right . \\ & & \left . + \left ( \coth ( \frac{\hbar \omega _ { 0}}{2k_{b}t})-1\right ) \left ( \omega _ { k}-\omega _ { 0}\right ) \exp \left [ -\frac{1}{2}\left ( \omega _ { k}-\omega _ { 0}\right ) ^{2}\tau _ { c}^{2}(k_{\perp } , k_{z})\right ] \right\ } . \nonumber\end{aligned}\ ] ] furthermore , we obtain @xmath99 and @xmath100 where @xmath101 + \right . \\ & & \left . + \exp \left [ -\frac{1}{2}\left ( \omega -\omega _ { k}-\omega _ { 0}\right ) ^{2}\tau _ { c}^{2}(k_{\perp } , k_{z})\right ] \right ) +   \nonumber \\ & & + \left ( \coth ( \frac{\hbar \omega _ { 0}}{2k_{b}t})-1\right ) \left ( \exp % \left [ -\frac{1}{2}\left ( \omega + \omega _ { k}-\omega _ { 0}\right ) ^{2}\tau _ { c}^{2}(k_{\perp } , k_{z})\right ] + \right . \nonumber \\ & & \left . + \exp \left [ -\frac{1}{2}\left ( \omega -\omega _ { k}+\omega _ { 0}\right ) ^{2}\tau _ { c}^{2}(k_{\perp } , k_{z})\right ] \right ) \right\ } ,   \nonumber\end{aligned}\ ] ] @xmath102 is the optical phonon frequency , @xmath103 , ( @xmath104 and @xmath105 are the hf and static permittivities of the crystal , respectively ) . the response and correlation functions of the acoustic phonons are given by @xmath106 which provide the contribution to the damping rate as @xmath107 and @xmath108 where @xmath109 + \right . \\ & & \left . + \left ( \coth ( \frac{\hbar uk}{2k_{b}t})-1\right ) \left ( \omega _ { k}-uk\right ) \exp \left [ -\frac{1}{2}\left ( \omega _ { k}-uk\right ) ^{2}\tau _ { c}^{2}(k_{\perp } , k_{z})\right ] \right\ } . \nonumber\end{aligned}\ ] ] accordingly , @xmath110 and @xmath111 where @xmath112 + \right . \\ & & \left . + \exp \left [ -\frac{1}{2}\left ( \omega -\omega _ { k}-uk\right ) ^{2}\tau _ { c}^{2}(k_{\perp } , k_{z})\right ] \right ) +   \nonumber \\ & & + \left ( \coth ( \frac{\hbar uk}{2k_{b}t})-1\right ) \left ( \exp \left [ -% \frac{1}{2}\left ( \omega + \omega _ { k}-uk\right ) ^{2}\tau _ { c}^{2}(k_{\perp } , k_{z})\right ] + \right .   \nonumber \\ & & \left . + \exp \left [ -\frac{1}{2}\left ( \omega -\omega _ { k}+uk\right ) ^{2}\tau _ { c}^{2}(k_{\perp } , k_{z})\right ] \right ) \right\ } . \nonumber\end{aligned}\ ] ] here , @xmath113 is the deformation potential , @xmath114 is the crystal density , and @xmath115 is the sound velocity .    for the case of static random impurities , there is not any response function , due to the lack of dynamics . the correlation function is given by @xmath116 where @xmath117 is the screening radius , @xmath118 , @xmath119 is the impurity concentration for species @xmath9 , and @xmath120 is their charge number . the relaxation rates associated with impurities have the forms @xmath121 and @xmath122 finally , the contributions to the fourier transforms of the fluctuation force  correlators can be written as @xmath123 + \exp \left [ -\frac{1}{2}\left ( \omega -\omega _ { k}\right ) ^{2}\tau _ { c}^{2}(k_{\perp } , k_{z})\right ] \right ) ,   \nonumber\end{aligned}\ ] ] and @xmath124 + \exp \left [ -\frac{1}{2}\left ( \omega -\omega _ { k}\right ) ^{2}\tau _ { c}^{2}(k_{\perp } , k_{z})\right ] \right ) ,   \nonumber\end{aligned}\ ] ]    in all these formulae we have used the notation @xmath125 @xmath126 @xmath127 . @xmath128 and @xmath129 can be determined self - consistently using @xmath130 and @xmath131 @xcite . all these results are obtained for relatively weak magnetic fields , @xmath132 . however , the energy shift , @xmath84 , due to electron - bath coupling is even smaller , @xmath133 , and will be neglected . in this section we examine the spin relaxation process in the presence of the  effective bath  with variables given by eq . ( [ q ] ) . their response functions and correlation functions have the form ( @xmath134 ) @xmath135_{+}\rangle\ ] ] @xmath136_{-}\rangle \eta ( t - t_{1})\ ] ] it is evident that these functions are of the sixth order of the electron velocity projections . using the quantum analog of the furutsu - novikov theorem for gaussian variables @xcite , they can be represented in terms of pair correlators . the resulting expressions for the spectral functions are given by @xmath137 @xmath138 _ { + } \right\rangle \left [ \left ( \left\langle v_{z}^{2}\right\rangle -\left\langle v_{x}^{2}\right\rangle \right ) ^{2}+4\left\langle v_{x}v_{y}\right\rangle \left\langle v_{x}v_{y}\right\rangle \right ] +   \label{sxx } \\ & & + \lambda ^{2}\int \frac{d\omega _ { 1}}{2\pi } \int \frac{d\omega _ { 2}}{2\pi } % \xi \left ( \omega _ { 1},\omega _ { 2},\omega -\omega _ { 1}-\omega _ { 2}\right ) \cdot    \nonumber \\ & & \cdot \left\ { 2\left\langle \frac{1}{2}\left [ v_{x}(\omega _ { 1}),v_{x}% \right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{x}(\omega _ { 2}),v_{x}\right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{x}(\omega -\omega _ { 1}-\omega _ { 2}),v_{x}\right ] _ { + } \right\rangle + \right . \nonumber \\ & & 2\left\langle \frac{1}{2}\left [ v_{x}(\omega _ { 1}),v_{x}\right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{z}(\omega _ { 2}),v_{z}% \right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{z}(\omega -\omega _ { 1}-\omega _ { 2}),v_{z}\right ] _ { + } \right\rangle -   \nonumber \\ & & \left . -4\left\langle \frac{1}{2}\left [ v_{x}(\omega _ { 1}),v_{x}\right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{x}(\omega _ { 2}),v_{y}% \right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{x}(\omega -\omega _ { 1}-\omega _ { 2}),v_{y}\right ] _ { + } \right\rangle \right\ } ,   \nonumber\end{aligned}\ ] ] @xmath139 _ { + } \right\rangle \left [ \left ( \left\langle v_{z}^{2}\right\rangle -\left\langle v_{x}^{2}\right\rangle \right ) ^{2}+4\left\langle v_{x}v_{y}\right\rangle \left\langle v_{x}v_{y}\right\rangle \right ] +   \label{sxy } \\ & & + \lambda ^{2}\int \frac{d\omega _ { 1}}{2\pi } \int \frac{d\omega _ { 2}}{2\pi } % \xi \left ( \omega _ { 1},\omega _ { 2},\omega -\omega _ { 1}-\omega _ { 2}\right ) \cdot    \nonumber \\ & & \cdot \left\ { 2\left\langle \frac{1}{2}\left [ v_{y}(\omega _ { 1}),v_{x}% \right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{y}(\omega _ { 2}),v_{x}\right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{y}(\omega -\omega _ { 1}-\omega _ { 2}),v_{x}\right ] _ { + } \right\rangle + \right . \nonumber \\ & & + 2\left\langle \frac{1}{2}\left [ v_{y}(\omega _ { 1}),v_{x}\right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{z}(\omega _ { 2}),v_{z}% \right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{z}(\omega -\omega _ { 1}-\omega _ { 2}),v_{z}\right ] _ { + } \right\rangle -   \nonumber \\ & & \left . -4\left\langle \frac{1}{2}\left [ v_{y}(\omega _ { 1}),v_{x}\right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{x}(\omega _ { 2}),v_{x}% \right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{x}(\omega -\omega _ { 1}-\omega _ { 2}),v_{x}\right ] _ { + } \right\rangle \right\ } ,   \nonumber\end{aligned}\ ] ] and @xmath140 _ { + } \right\rangle \left\ { \left\langle \frac{1}{2}\left [ v_{x}(\omega _ { 2}),v_{x}\right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{x}(\omega -\omega _ { 1}-\omega _ { 2}),v_{x}\right ] _ { + } \right\rangle -\right .    \nonumber \\ & & \left . -\left\langle \frac{1}{2}\left [ v_{x}(\omega _ { 2}),v_{y}\right ] _ { + } \right\rangle \left\langle \frac{1}{2}\left [ v_{x}(\omega -\omega _ { 1}-\omega _ { 2}),v_{y}\right ] _ { + } \right\rangle \right\ } ,   \nonumber\end{aligned}\ ] ] where we defined an auxiliary function @xmath141 which appears due to the replacement of the velocity commutators by their anticommutators in accordance with eq . ( [ coth ] ) . in eqs . ( [ sxx]-[szz ] ) we use the shortened notation @xmath142 _ { + } \right\rangle = \int \frac{d\omega } { 2\pi } \left\langle \frac{1}{2}\left [ v_{j}(\omega ) , v_{j}\right ] _ { + } \right\rangle , $ ] and @xmath143 _ { + } \right\rangle = \int \frac{d\omega } { 2\pi } \left\langle \frac{1}{2% } \left [ v_{x}(\omega ) , v_{y}\right ] _ { + } \right\rangle $ ] .employing a second application of the toqs method , we obtain equations for the average spin projections as follows : @xmath144 where the coefficients are given by @xmath145 + m_{yx}\left ( t , t_{1}\right ) \sin \left [ \omega _ { b}\left ( t - t_{1}\right ) \right ] + m_{zz}\left ( t , t_{1}\right ) \right\ } , \\ \gamma _ { xy } & = & \frac{4}{\hbar ^{2}}\int_{-\infty } ^{t}dt_{1}\left\ { m_{yy}\left ( t , t_{1}\right ) \sin \left [ \omega _ { b}\left ( t - t_{1}\right ) % \right ] -m_{yx}\left ( t , t_{1}\right ) \cos \left [ \omega _ { b}\left ( t - t_{1}\right ) \right ] \right\ } ,   \nonumber \\ \gamma _ { xz } & = & \frac{4}{\hbar ^{2}}\int_{-\infty } ^{t}dt_{1}\left\ { m_{zx}\left ( t , t_{1}\right ) \cos \left [ \omega _ { b}\left ( t - t_{1}\right ) % \right ] -m_{zy}\left ( t , t_{1}\right ) \sin \left [ \omega _ { b}\left ( t - t_{1}\right ) \right ] \right\ } ,   \nonumber \\ \gamma _ { x}^{0 } & = & \frac{2}{\hbar } \int_{-\infty } ^{t}dt_{1}\left\ { \varphi _ { zx}\left ( t , t_{1}\right ) \sin \left [ \omega _ { b}\left ( t - t_{1}\right ) % \right ] + \varphi _ { zy}\left ( t , t_{1}\right ) \cos \left [ \omega _ { b}\left ( t - t_{1}\right ) \right ] -\varphi _ { yz}\left ( t , t_{1}\right ) \right\ } ,   \nonumber\end{aligned}\ ] ] @xmath146 + m_{xy}\left ( t , t_{1}\right ) \cos \left [ \omega _ { b}\left ( t - t_{1}\right ) \right ] \right\ } , \\ \gamma _ { yy } & = & \frac{4}{\hbar ^{2}}\int_{-\infty } ^{t}dt_{1}\left\ { m_{xx}\left ( t , t_{1}\right ) \cos \left [ \omega _ { b}\left ( t - t_{1}\right ) % \right ] -m_{xy}\left ( t , t_{1}\right ) \sin \left [ \omega _ { b}\left ( t - t_{1}\right ) \right ] + m_{zz}\left ( t , t_{1}\right ) \right\ } ,   \nonumber \\ \gamma _ { yz } & = & \frac{4}{\hbar ^{2}}\int_{-\infty } ^{t}dt_{1}\left\ { m_{zx}\left ( t , t_{1}\right ) \sin \left [ \omega _ { b}\left ( t - t_{1}\right ) % \right ] + m_{zy}\left ( t , t_{1}\right ) \cos \left [ \omega _ { b}\left ( t - t_{1}\right ) \right ] \right\ } ,   \nonumber \\ \gamma _ { y}^{0 } & = & \frac{2}{\hbar } \int_{-\infty } ^{t}dt_{1}\left\ { \varphi _ { zy}\left ( t , t_{1}\right ) \sin \left [ \omega _ { b}\left ( t - t_{1}\right ) % \right ] -\varphi _ { zx}\left ( t , t_{1}\right ) \cos \left [ \omega _ { b}\left ( t - t_{1}\right ) \right ] + \varphi _ { xz}\left ( t , t_{1}\right ) \right\ } ,   \nonumber\end{aligned}\ ] ] and @xmath147 + \right . \nonumber \\ & & \left . + \left ( m_{yx}\left ( t , t_{1}\right ) -m_{xy}\left ( t , t_{1}\right ) \right ) \sin \left [ \omega _ { b}\left ( t - t_{1}\right ) \right ] \right\ } ,   \nonumber \\ \gamma _ { z}^{0 } & = & \frac{2}{\hbar } \int_{-\infty } ^{t}dt_{1}\left\ { \left ( \varphi _ { yx}\left ( t , t_{1}\right ) + \varphi _ { xy}\left ( t , t_{1}\right ) \right ) \cos \left [ \omega _ { b}\left ( t - t_{1}\right ) \right ] -\right .    \nonumber \\ & & \left . -\left ( \varphi _ { xx}\left ( t , t_{1}\right ) + \varphi _ { yy}\left ( t , t_{1}\right ) \right ) \sin \left [ \omega _ { b}\left ( t - t_{1}\right ) \right ] \right\ } . \nonumber\end{aligned}\ ] ] only six of the above twelve coefficients , @xmath148 , do not vanish , and we may identify relaxation times in terms of the nonvanishing @xmath149 as follows : @xmath150 where @xmath151 are the longitudinal and transverse relaxation times , responsible for the relaxation of the magnetic moment and decoherence respectively . furthermore , direct calculation shows that @xmath152 where @xmath153 may be identified as a frequency shift due to  effective bath  fluctuations , and @xmath154 is the equilibrium population difference between spin up and spin down states . with these calculated results and identifications , the spin equations take the usual bloch form , @xmath155    in terms of the spectral functions , @xmath156 , the relaxation rates ( inverse relaxation times ) , may be expressed as follows : @xmath157 and @xmath158 although eqs . ( [ t1 ] ) , ( [ t2 ] ) are relatively simple and yield some interesting qualitative conclusions , a quantitative analysis of them is difficult , because of the complexity of the expressions for the spectral functions , @xmath156 . this complexity is relieved under the prevailing assumption of weak coupling , which permits the replacement of the lorenzians involved in the integrands of eqs . ( [ sxx]-[szz ] ) by @xmath159 where @xmath160 is the dirac delta - function , and , consequently , @xmath161 ( in this , we note that in eqs . ( [ xx]),([xy ] ) and ( [ zz ] ) , we have @xmath162 all proportional to coupling strength @xmath163 ) . having employed eqs . ( [ lorenz ] ) , ( [ 1 ] ) in the integrals of eqs . ( [ sxx]-[szz ] ) representing @xmath164 @xmath165 and @xmath166 , we obtain the results of integration as : @xmath167 _ { + } \right\rangle \left [ \left ( \left\langle v_{z}^{2}\right\rangle -\left\langle v_{x}^{2}\right\rangle \right ) ^{2}+4\left\langle v_{x}v_{y}\right\rangle \left\langle v_{x}v_{y}\right\rangle \right ] +   \label{2 } \\ & & + \frac{\alpha ^{2}m^{3}}{8e_{g}}\left [ r_{1}(\omega ; \omega _ { c})+y_{1}(\omega ; \omega _ { c})\right ] ,   \nonumber\end{aligned}\ ] ] @xmath168 _ { + } \right\rangle \left [ \left ( \left\langle v_{z}^{2}\right\rangle -\left\langle v_{x}^{2}\right\rangle \right ) ^{2}+4\left\langle v_{x}v_{y}\right\rangle \left\langle v_{x}v_{y}\right\rangle \right ] + \\ & & + \frac{\alpha ^{2}m^{3}}{8e_{g}}\left [ r_{2}(\omega ; \omega _ { c})+y_{2}(\omega ; \omega _ { c})\right ] ,   \nonumber\end{aligned}\ ] ] @xmath169 where we have introduced the notation @xmath170 @xmath171 @xmath172 and @xmath173 these results are very general , providing explicit expressions for the spin relaxation times in terms of the material constants and coupling strengths . the magnetic field has been treated as relatively weak ( @xmath174 ) , but our estimates show that this range of fields includes those used in most experiments .    to make contact with earlier theories that did not take account of the magnetic field and used a phenomenological momentum scattering time @xcite , we put @xmath175 and @xmath176 in eqs . ( [ 2]-[3 ] ) , which immediately yields @xmath177 ^{2}}{\gamma _ { 0}^{3}}k(\omega ) \lambda ( \omega ; 3\alpha \gamma _ { 0}),\ ] ] where @xmath178 and @xmath179 in this zero field limit the relaxation times are equal and determined by @xmath180 in equilibrium , the rate of velocity fluctuations is defined only by temperature , and we obtain the result @xmath181 where @xmath182 is the kelvin temperature . this yields the common relaxation rate as @xmath183 which may be brought into coincidence with the expression usually used for the dp mechanism @xcite by the change of notation @xmath184 , with @xmath185 a phenomenological momentum scattering time . it should be noted that the choice of coefficient @xmath186 is unclear , with various authors assigning values in the range @xmath187 , with @xmath188 for scattering by deformational acoustic phonons and @xmath189 for scattering by polar optical phonons @xcite . on the other hand , our microscopic analysis provides a clearly defined value with account of all scattering mechanisms . moreover , we provide results for finite  magnetic field derived on a microscopic basis , rather than relying on the ansatz of ref . @xcite , which has @xmath190 , @xmath191 in contradistinction to this , our analysis in the presence of a magnetic field shows that the transverse and longitudinal spin relaxation times are not equal ( eqs . ( [ t1],[t2 ] ) ) . in summary , we have analyzed electron spin relaxation dynamics and decoherence in bulk semiconductors with an applied external magnetic field . in the absence of direct coupling of spin degrees of freedom to the phonon / impurity heat bath , we have used the orbital motion as an intermediary , transferring fluctuations from the bath to electron spin and , transmitting the dissipated energy from spin to the bath . to accomplish this , the two - stage procedure of solution has been employed : in the first stage , the fluctuation characteristics of the orbital motion were determined and the orbital degrees of freedom were used in the second stage as an  effective heat bath  for spin dynamics . in both stages of the analysis , we employed the method of the theory of open quantum systems , obtaining a set of bloch equations having longitudinal and transverse spin relaxation times , @xmath0 and @xmath1 , determined on a microscopic basis . our results provide analytic definition of the phenomenological parameters employed in earlier zero - field theories , in terms of material constants and coupling strengths . moreover , this analysis of electron spin dynamics in finite magnetic field yields explicit formulae for the longitudinal and transverse relaxation times . kikkawa and d.d . awschalom , phys . lett . * 80 * , 4313(1998 ) ; t.f . bogges , j.t . olesberg , c. yu , m.e . flatt , and w.h . lau , appl . lett . * 77 * , 1333 ( 2000 ) ; b. beschoten , e. johnston - halperin , d.k . young , m. poggio , j.e . grimaldi , s. keller , s.p . denbaars , u.k . mishra , e.l . hu , and d.d . awschalom , phys . b * 63 * 121202(r ) ( 2001 ) . lau , j.t . olesberg , and m.e . flatt , phys . b * 64 * , 161301(r ) ( 2001 ) ; a. bournel , p. dolfus , e. cassan , and p. hesto , appl . . lett . * 77 * 2346 ( 2000 ) ; d.m . frenkel , phys . b * 43 * 14228 ( 1991 ) ; p.h . song and k.w . kim , arxiv : cond - mat/0111076 ( 2001 ) . wolf , d.d . awschalom , r.a . buhrman , j.m . daughton , s. von molnr , m.l . roukes , a.y . chtchelkanova , d.m . treger , science * 294 * , 1488 ( 2001 ) ; g.a . prinz , science * 282 * , 1660 ( 1998 ) ; y. nishikawa , a. tackeuchi , s. nakamura , s. muto , and n. yokoyama , appl . . lett . * 66 * , 839 ( 1995 ) . m.i . dyakonov and v.i . perel , sov . jetp * 33 * , 1053 ( 1977 ) ; _ optical orientation , modern problems in condensed matter science _ , edited by f. meier and b.p . zakharchenya ( north - holland , amsterdam , 1984 ) , vol . 8 .
Abstract:

Answer: we employ the method of the theory of open quantum systems to analyze spin relaxation and decoherence in semiconductors in the presence of a magnetic field . we derive a set of bloch equations for electron spin with a fully microscopic determination of longitudinal and transverse relaxation times . electron scattering from optical and acoustic phonons and random impurities is taken into account . we obtain explicit expressions for the spin relaxation times in terms of material constants and coupling strengths , exhibiting formal agreement with earlier treatments in the zero magnetic field limit with microscopic specification of their phenomenological parameters .


Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: localization of a photonic mode by disorder in a statistically uniform and isotropic medium has been studied for more than two decades  @xcite , both in the optical and microwave spectral ranges . in this paper we propose and demonstrate a new mechanism for disorder induced strong three dimensional ( 3d ) localization of light . we achieved this in planar microcavities ( pmcs ) , containing a layer of strain - induced self - assembled quantum dots ( qds )  @xcite . the dispersion of light in the pmc plane is modified by the distributed bragg reflectors ( dbrs ) that trap the modes in perpendicular direction to the cavity plane . the qds constitute an active material which , under proper excitation , emits photons , due to recombination of qd confined excitons  @xcite . randomness in the qds in - plane position , composition and sizes produces fluctuations in the system s dielectric constant . this randomness provides the necessary `` attractive potential '' , which , together with the modified dispersion , brings about strong photon localization in the pmc plane . our novel mechanism differs from previous studies , since it provides genuine disorder - induced localization rather than leaky ( resonant ) modes  @xcite . it also differs from the `` transverse localization ''  @xcite , where the wave freely propagates in one direction while being trapped , by disorder , in the transverse direction  @xcite . we directly demonstrate the 3d localization of the photonic modes by measuring their photoluminescence ( pl ) intensity distribution . it is worthwhile to emphasize from the start that the comparatively small size of the excitation spot ( @xmath0 , or smaller ) plays a crucial role in the interpretation of the experimental data . the point is that _ in the absence of disorder _ the shape of the pl spectrum emerging from the cavity is quite sensitive to the size of the excitation spot . if the size is sufficiently small , the pl spectrum becomes very narrow , in marked difference with the broad shoulder predicted by the standard basic cavity physics. if there had been no disorder , the experimental conditions of our paper would have resulted in an extremely narrow pl peak . the significant broadening of the spectrum , as observed in our experiment , is due entirely to disorder . the broadening occurs both towards frequencies below and above the cavity frequency ( localized and resonant photonic modes , respectively ) . all this will be discussed in detail in the theoretical part of the paper . our sample was grown by molecular beam epitaxy on a [ 1,0,0 ] oriented gaas substrate . the sample consists of a pmc formed by a gaas one - wavelength resonator , sandwiched between dbrs made from 23(26 ) top ( bottom ) alternating gaas and alas quarter - wavelength - layer pairs . as an active material in the resonator s anti - node , a single strain - induced self - assembled in@xmath1ga@xmath2as qds layer with average density of 5@xmath3 qds per @xmath4 is used  @xcite . the use of such composition results in large and asymmetric qds with typical lengths and widths of about 100 and 30 nm respectively , as can be seen in the inset to fig . [ fig : one]a . we used two methods to study the pmc . a ) a diffraction limited confocal optical scanning microscopy ( a @xmath560 objective with numerical aperture of 0.85 , obtaining spatial resolution of @xmath6 ) . b ) a near field scanning optical microscopy ( nsom- nanonics cryoview @xmath7 with spatial resolution of @xmath8 ) . by applying either method we measure the lateral distribution of the electromagnetic field above the sample surface . two modes of excitation were used . in the first , the excitation was focused to a diffraction limit by the collecting objective  @xcite . in the second ( the only one used for nsom ) , the excitation was focused at an oblique angle to a spot diameter of about 30 @xmath9 . due to the high energy used for the excitation ( 632.8 nm light of a hene laser ) the images produced by the two excitation modes were almost identical . in fig . [ fig : one]a we present pl spectrum ( upper , black solid line ) from a single point on the surface of the pmc sample together with a pl spectrum which was taken from the same sample ( lower , solid blue line ) after the upper dbr mirror was completely etched . the spectrally broad blue pl line is centered at 1.33 ev and has full width at half maximum ( fwhm ) of approximately 10 mev . this spectral line results from s - shell electron - hole recombination within the inhomogeneous population of qds of various sizes and compositions . the radiative linewidth of a qd is @xmath10  @xcite . = .52    the rather symmetrical shape of the spectral line indicates that charge carriers are not distributed thermally between the qds . therefore , it is quite safe to assume that the spectrum accurately represents the actual energy distribution of the emission from individual qds . the spectrally sharp pl line in the figure is centered at 1.319 ev and has fwhm of 80 @xmath11ev . as we explain below , in a disorder free pmc this line width should have been about an order of magnitude narrower due to the finite excitation area and diffraction limited collection spot . the asymmetric spectral shape of the line and its actual linewidth which represents a q - factor of @xmath1215000 , in similarity to earlier measurements  @xcite , are due to the disorder in the dielectric constant as we quantitatively show below .    in the left inset to fig . [ fig : one]a we present on a semi logarithmic scale a spatially integrated ( ` far field ' ) pl spectrum from a 135@xmath9 wide @xmath13 long area on the sample s surface . the spectrum was obtained by summing up large number of individual pl spectra , each obtained from a diffraction limited areal spot . the dominating spectral feature is a large , asymmetrically broadened line around the energy of the pmc mode . the line decays exponentially towards lower energies with a characteristic energy of @xmath14 . towards higher energies the line decays with a characteristic energy of @xmath15 . the ratio between these characteristic energies agrees well with the prediction of our theoretical model for the disorder in the dielectric constant ( see below ) .    in figs . [ fig : one]b-[fig : one]d selective wavelength images of the scanned area are presented . the spectral domains which are used for each image are indicated by the green zones imposed on the far field spectrum in fig . [ fig : one]a . the image in fig . [ fig : one]b is obtained within a 0.5 mev energy window containing the pmc mode . it shows almost evenly distributed emission from the surface . in contrast , the image in fig . [ fig : one]c , which is obtained within a 0.5 mev window , 1 mev below the pmc mode , shows emission emerging from randomly distributed spots on the surface . likewise , the image in fig . [ fig : one]d , which is obtained from energy window of same width , located 3 mev lower , shows only 2 bright centers , from which all the emission results . = .46    in fig . [ fig : two ] we turn our attention to one of the low energy electromagnetic modes seen in fig .  [ fig : one]d . such localizing centers were observed in various mode energies and various locations on the sample surface . in the right panels of figs . [ fig : two]a-[fig : two]e , we display selective wavelength images obtained by scanning the aluminum coated fiber tip of our nsom microscope in tapping mode over the sample surface . the corresponding spectra in the left panels , are obtained from the brightest area pixels indicated by the + signs on the images . the green areas on the spectra indicate the spectral domains used for generating these images . these spectra and corresponding images reveal the following observations : four different localized photonic modes are observed at energies below that of the pmc mode . each mode is composed of two cross linearly polarized equally intense components ( not shown ) . the components are @xmath16 apart and the lower energy one is polarized along the [ 0,-1,1 ] direction  @xcite . the lowest energy mode is about 5 mev below the energy of the pmc mode . the in - plane spatial intensity distribution of this mode , as depicted by the image in fig . [ fig : two]a , is node - less , and has an elongated shape oriented along the [ 0,-1,1 ] crystallographic direction , along which the qds are also elongated ( see inset to fig .  [ fig : one ] ) . the mode dimensions are roughly 3@xmath52@xmath17 . the peak intensity is about 40 times stronger than that of the pmc mode . the next three images in fig . [ fig : two]b , c and d , present elongated spatial intensity distributions composed of one , two and three nodes , respectively . each mode s energy is @xmath121.5 - 2 mev higher than the energy of the preceding mode . the selective wavelength image in fig . [ fig : two]e is obtained at an energy window , which resides on the higher energy side of the pmc mode . it displays evenly distributed emission from everywhere except from the actual locus of the localized modes , where extended modes are forbidden . the magnitude and shape of the spot from which emission at this wavelength is missing is an estimate for the spatial extent of the region which localizes the optical modes . the region is comparable in shape and size to the most localized , lowest energy mode . the localized modes are insensitive to temperature and excitation density variations , hence they do not result from thermal or non - linear effects . the emission is linear with the excitation density between 1-@xmath18 watt/@xmath19 . increasing the temperature from 10 to 80k results in a small , rigid , red spectral shift , as expected from the temperature dependence of the material bandgap . since the sizes of the spots in fig . 2b-2d significantly exceed the size of an individual qd , the only way to account for the experimental observations is to assume that each spot corresponds to a _ localized _ electromagnetic mode . therefore the low frequency tail of the pl spectrum in fig . 1a should be interpreted as a tail of localized modes with energies below that of the pmc mode , @xmath20 , similar to the tail of localized electron states formed below the band - edge due to disorder  @xcite . to demonstrate how the presence of a cavity leads to the in - plane trapping , we consider the scalar wave equation corresponding to the frequency @xmath21 : [ eq : one1 ] d^2dz^2 + ^2 _ + ^2 c^2 = 0 , where @xmath22 is a uniform dielectric constant inside the cavity , @xmath23 , and @xmath24 is the perturbation due to the randomly positioned qds that reside in a narrow layer , of width @xmath25 nm , in the middle of the cavity . let us emphasize that an individual qd is much too small to be able to trap a photon . actually , in the continuum description employed in eq . ( [ eq : one1 ] ) individual qds are not resolved . what matters are the fluctuations in the areal density of the qds , which manifest themselves in fluctuations of the dielectric constant @xmath26 . since the electromagnetic modes are confined in the @xmath27-direction to a much shorter length than in the plane , the function @xmath28 can be factorized as @xmath29 , where @xmath30 is the cavity width , and @xmath31 satisfies the equation [ eq : two1 ] - c^2 ^2 _ - ^2 ( ) = - _ 0 ( _ c^2 - ^2 ) . in this equation , @xmath32 is the effective in - plane fluctuation of the dielectric constant obtained from @xmath33 by averaging over the @xmath27-distribution of the cavity mode intensity , namely , @xmath34 . the analogy between eq . ( [ eq : two1 ] ) and the schrdinger equation makes it clear that , for positive @xmath35 , eq .  ( [ eq : two1 ] ) admits localized solutions , with frequency @xmath21 smaller than the cavity frequency @xmath36 . thus , even a small enhancement in @xmath35 acts as an attractive potential localizing the mode laterally . to substantiate that this mechanism of trapping is relevant to our experimental observations , we note that the spatial intensity distributions and the energies of the modes in fig . [ fig : two ] are quite accurately fitted by solving eq . ( [ eq : two1 ] ) using a one dimensional parabolic model @xmath37 $ ] , @xmath38 describing a local fluctuation in the dielectric constant ( see fig .  [ fig : two]e , right panel , where the best fitted values are @xmath39 and @xmath40 ) . this agreement , and the fact that even in the case of the most localized modes , differently polarized components have almost the same energy , justifies the use of the scalar approximation .    to find the shape of the tail of the localized modes distribution , we have to specify the disorder `` potential '' , @xmath41 . the simplest possible model is the `` white - noise potential '' which amounts to treating the qds as point - like resonant units , randomly distributed in the plane ( this simple model does not describe the anisotropy observed in the inset of fig . [ fig : one ] as discussed below ) . fluctuations in the local dielectric constant are given as [ eq : three1 ] ( ) = a , where @xmath42 is the deviation of the number of qds ( per unit area near @xmath43 ) , with resonance frequency in the interval @xmath44 , from the average value @xmath45 . here @xmath46 is the total areal density of the qds , regardless of their frequency and @xmath47 is the probability distribution of the resonance frequencies ( the solid , blue line in fig . [ fig : one ] ) .    to estimate the constant @xmath48 , we consider the @xmath49 limit and imagine that the entire layer consists of qds , i.e. of in@xmath1ga@xmath50as . in this case , the enhancement of the dielectric constant in the layer , as compared to the surrounding material of gaas , is @xmath51 . this value should be multiplied by @xmath52 , which is the ratio of the qds layer width ( 5 nm ) to the cavity width ( @xmath12300 nm ) . there is an extra factor of 2 , since the qds are located at the anti - node of the pmc mode . thus , the effective enhancement , due to the presence of inas , is @xmath53 . the constant @xmath48 can now be estimated from eq . ( [ eq : three1 ] ) as @xmath54 , where @xmath55 is the typical area of a qd and @xmath56 is the central frequency of the distribution @xmath47 . note that the employed averaging in the @xmath27-direction , as well as our use of the scalar approximation is justified by the smallness of the parameter @xmath57 . indeed , the effect of disorder - induced coupling between different cavity modes is proportional to @xmath58 , whereas the effect of the qd layer on a given mode is proportional to @xmath57 . the function @xmath59 is a complicated , random function of position @xmath43 and frequency @xmath60 . using the white - noise model for @xmath61 and eq . ( [ eq : three1 ] ) , we obtain the probability @xmath62 for a particular realization @xmath63 of the fluctuating part of the dielectric constant @xmath64 ^ 2 \right\}.\ ] ] formally , the integral over @xmath60 in eq . ( [ eq : four1 ] ) diverges at @xmath65 . in this regard , it should be noted that the integral over @xmath60 is actually an approximation for a sum over discrete values , @xmath66 , of the resonant frequency of the dots in the relevent area ( of order 1  @xmath17 ) . the meaning of this area can be traced back to eq . ( [ eq : three1 ] ) which defines a macroscopic ( albeit still fluctuating in space ) dielectric constant . the spatial scale , over which @xmath63 is changing , should be at least few times larger than the  microscopic " distances ( size of the dots and their separation , of the order of 0.1 @xmath67 . thus , the minimal scale at which @xmath68 can be meaningfully defined should be in the sub micron range . we are interested in @xmath21 below @xmath20 , i.e. in frequency which belongs to the tail of the distribution @xmath69 . a simple estimate shows that for such @xmath21 the probability to encounter a dot with @xmath66 close to @xmath21 is extremely small . moreover , for such @xmath21 the integral over @xmath60 in eq . ( [ eq : four1 ] ) can be , roughly , replaced by @xmath70 ( note that such a replacement would not be possible for @xmath21 in the  bulk " of the probability distribution @xmath69 ) . the two equations , ( [ eq : two1 ] ) and ( [ eq : four1 ] ) , define the statistical problem of finding the probability that a localized mode , in a given frequency interval , will be created . this problem is completely identical to the two - dimensional version  @xcite of the problem of the tails of electron states in a random potential  @xcite . the final result for the probability @xmath71 to find a trap , in an area @xmath72 , which can support a localized mode in the frequency range @xmath73 can be presented as [ eq : five1 ] p ( ) , where the pre - exponential factor corresponds to the density of the pmc modes at frequencies @xmath74 . ( [ eq : five1 ] ) can be viewed as a qualitative interpolation between the maximal value of the modes density at @xmath75 , and the low - density `` tail '' . a more detailed crossover behavior between the tail at @xmath76 and @xmath77 at @xmath78 can be found in ref .  . important is that this crossover is determined by a _ single _ scale of frequencies @xmath79 which follows from eq . ( [ eq : five1 ] ) , so that in the tail region @xmath80 can be rewritten as @xmath81 $ ] . in eq . ( [ omegat ] ) we replaced @xmath82 by @xmath83 , since , experimentally , the difference @xmath83 is much bigger than the tail width , @xmath84 .    with the above estimate for the constant @xmath48 , eq . ( [ eq : five1 ] ) yields a characteristic energy of @xmath85 for the initial decay of the pl intensity below @xmath86 . it agrees well with the measured value ( dashed red line in the left inset to fig . [ fig : one]a ) . the relatively large uncertainty results from the uncertainties in the qds density and their average area .    up to now we have focused on explaining the pl spectrum tail for @xmath87 . in this domain , the pl intensity , @xmath88 , outside the cavity is simply proportional to the density of trapped modes , @xmath89 . in other words , the @xmath21-dependence of the light intensities inside and outside the cavity is the same for @xmath90 . this is _ not the case _ for the domain of frequencies @xmath91 to which we now turn . to calculate @xmath88 in this domain , we start by recapitulating the basic microcavity physics , in the absence of disorder . in this case the emission spectrum is controlled by the airy factor  @xcite @xmath92\vert^2},\ ] ] where @xmath93 is the angle at which the radiation exits the cavity , @xmath94 is the radiation frequency , @xmath95 ( @xmath96 ) is the reflection amplitude from the upper ( lower ) bragg mirror and @xmath97 is the transmission coefficient ( @xmath95 and @xmath96 are taken to be real ) . the phase @xmath98 is acquired in course of propagation between the boundaries @xmath99 . at frequencies , @xmath21 , close to @xmath20 eq . ( [ airy ] ) has a sharp maximum at @xmath100 , where the angle @xmath101 is defined as @xmath102^{1/2},\ ] ] so that @xmath103 can be simplified to @xmath104 note that condition @xmath100 coincides with the dispersion law of the waveguided mode . this follows from eq . ( [ eq : two1 ] ) in the absence of disorder . indeed , for in - plane wave vector , @xmath105 , eq . ( [ eq : two1 ] ) yields @xmath106 , leading to @xmath100 for @xmath107 . it follows from eq . ( [ simplified ] ) that the emission spectrum , @xmath88 , has a form of a plateau , which corresponds to the basic microcavity physics . indeed , integration over @xmath93 yields @xmath108.\ ] ] the second term describes a slight smearing of the plateau near @xmath109 . let us consider now photons which exit at an angle @xmath93 , from some small collection spot . it is important to realize that there is a correlation between the angle @xmath93 and the average distance @xmath110 ( measured from the collection spot ) at which the corresponding photons have been created . indeed , a typical photon will bounce @xmath111 times before escaping the cavity . therefore photons that bounce at angle @xmath93 and escape at the collection spot must be created at a distance [ l_theta ] l ( ) ~t^-1 d . thus , the finite size of an _ excitation _ area can serve as an efficient cutoff for a maximal escape angle @xmath112 . assuming a circular excitation area , of radius @xmath113 , and choosing a small collection spot at the center of the circle , one finds a maximal escape angle [ theta_m ] _ m = . of course , this relation is valid only if @xmath112 ( and therefore @xmath113 ) is not too large , namely , it should be smaller than the maximal escape angle for a fully excited cavity . since the angular cutoff translates into a frequency cutoff @xmath114 ( see eq . ( [ theta ] ) ) , we arrive at the important conclusion that the frequency range of the collected photons can become very narrow , if the excitation area is small enough . in our experiment @xmath115 , @xmath116 , so that the collected radiation constitutes a narrow peak at near @xmath20 , of relative width [ peak_width ] = = ( ) ^2 110 ^ -5 . obviously , this peak can be replaced by a @xmath57-function [ i0 ] i()=i_0(- _ c ) . qualitatively , the meaning of eq . ( [ i0 ] ) is that , due to the finite excitation spot , only the photons emitted normally to the mirror emerge outside the cavity . quantitatively we note that the relative width of the experimentally measured pmc mode ( @xmath117 ) is about an order of magnitude wider than the one calculated by eq . ( [ peak_width ] ) . as we show below , this broadening and the actual spectral shape of the pmc mode , results from the disorder . the replacement of the step - like spectrum of the _ clean _ cavity ( for infinite excitation spot ) by the @xmath57-function eq . ( [ i0 ] ) ( for finite excitation spot ) does not affect the emission spectrum of _ disordered _ cavity in the domain @xmath118 . in this domain the dependence , @xmath88 , is still determined by eq . ( [ eq : five1 ] ) . however , in the domain @xmath119 , the @xmath57-peak for the clean cavity gets broadened due to scattering of _ extended modes _ by the disorder .    in the presence of disorder , the above reasoning , which is based on the ray picture , is not applicable since _ all _ the solutions , @xmath120 , of eq . ( [ eq : two1 ] ) corresponding to the eigenstates , @xmath121 , have components with @xmath122 ( in the ray picture they correspond to rays with normal incidence to the microcavity plane ) . the amplitude of this component is is given by @xmath123 . therefore , the shape of the spectrum acquires the form @xmath124 where the function @xmath125 is defined as @xmath126 we note that the function @xmath125 as defined by eq . ( [ i1 ] ) is normalized @xmath127 thus , the disorder broadens the @xmath57-peak ( eq . ( [ i0 ] ) ) without changing its area . note also , that the calculation of the spectral shape , @xmath125 , is formally reduced to the calculation of the spectrum of two - dimensional _ exciton _ absorption in the presence of a white - noise random potential . this calculation was done earlier in ref . , where the zero - momentum components of the _ exciton _ states were selected by the matrix element for light absorption  @xcite . it is also important to note that the spectral shape , @xmath128 , is a function of the ratio @xmath129 , where @xmath84 is the same characteristic frequency , ( eq . ( [ omegat ] ) ) , which determines the spectral shape of the tail of the localized photon states . indeed , deep in the tail , @xmath130 , all the trapped solutions , @xmath120 , at a given @xmath131 are identical to each other , and eq . ( [ i1 ] ) is essentially reduced to the density equation of the trapped modes @xmath132 $ ] ( eq . ( [ eq : five1 ] ) ) .    in the opposite limit @xmath133 the density of photon states is constant . here , therefore , the decay of @xmath128 is determined by the fall - off of the integral @xmath123 . this integral can be estimated to first order by perturbation expansion in @xmath134 . thus it is proportional to @xmath135 . accurate calculation@xcite yields @xmath136 . this way , the width of the emission spectrum was found to be @xmath137 . the overall spectral shape of @xmath125 was found in ref . by combining the normalized ( eq . ( [ i2 ] ) ) low - frequency and high - frequency asymptotes . the calculated full width at half - maximum ( fwhm ) of the emission spectrum is @xmath137 , in excellent agreement with the experiment , where it is measured to be @xmath138 . the calculated ratio between the high energy and low energy slopes is 2:1 in very good agreement with the measured ratio of  2.5 . while the white noise model is successful in explaining the initial drop of the density of localized modes , it can not account for the experimental data in the deeper tail , starting around @xmath139 ev and below . it does not explain the non - monotonic decay of the probability as seen in fig . [ fig : one]a . it fails in predicting the observed multiplicity of confined modes which belong to the same localization region ( fig .  [ fig : two ] ) , and it does not explain the anisotropy in the electromagnetic mode shapes . the failure of the white noise model , for @xmath118 , is not surprising . indeed , when the mode frequency becomes smaller , localization of the mode becomes tighter . thus , various local features of the random potential - such as correlations in @xmath134 , the `` granular '' structure of the disorder and the inherent anisotropy of the system - become more important . obviously , all these features are not captured by the entirely featureless , universal white noise model .    proper account of these specific features of randomness might explain some of the data for the deeper @xmath21-tail . first , once correlations are introduced , the most probable deep fluctuation is no longer a fluctuation which captures one mode only . second , the probability of finding a localized mode common to two localizing centers , rather than to one , may prevail at a given energy below the energy of the pmc mode . indeed , localized modes between 1.317 ev and 1.318 ev always exhibit a spatial structure consisting of two bright spots separated by a dark region . this observation suggests existence of traps of a double - well shape . the origin of such traps , as well as the multimode traps in the deep tail ( below 1.316 ev ) , may be traced to the granular nature of the disorder . although , on the average , the qds are uniformly distributed , there exist configurations with several dots coming close to each other , within a distance of a dot size . such rare configurations constitute more efficient traps than the white - noise fluctuations and therefore dominate in the deeper tail . the clusters of qds are akin to clusters of potential wells in the lifshitz model of a disordered electronic system  @xcite . thus , anisotropic qd clustering , can account for the shape and multiplicity of the strongly localized modes in the deeper tail . we are unable to make quantitative predictions for the density and spatial structure of these modes , since no comprehensive theory exists . yet , we show that dielectric constant enhancement of order @xmath140 over distance of order @xmath141 , does account for the experimental observations ( fig . [ fig : two]e ) .    in conclusion , we demonstrate disorder - induced trapped photonic modes in a microcavity with an embedded layer of qds . in the lateral in - plane direction the modes are localized by spatial fluctuations of the dielectric constant due to randomness in the location and composition of the qds . our theory emphasizes the universal features of the trapping mechanism : the necessary combination of mirrors in the vertical direction and the lateral disorder . the theory quantitatively estimates the characteristic energy by which the pl intensity decays below the pmc mode . since the non universal , system - specific details of the randomness are not accurately known , the spatial shapes and energies of the modes deeper in the energy tail , are qualitatively discussed only . in addition we show both experimentally and theoretically that if the size of the excitation spot is sufficiently small , as it is in confocal and near field scanning microscopy , the microcavity emission spectrum becomes very narrow , in marked difference with the prediction of the standard cavity considerations , which apply to uniform excitation . the disorder broadens this spectrum , as well . thus , we quantitatively measure and calculate the disorder induced broadening below and above the cavity mode frequency .
Abstract:

Answer: we present the first experimental evidence for the existence of strongly localized photonic modes due to random two dimensional fluctuations in the dielectric constant . in one direction , the modes are trapped by ordered bragg reflecting mirrors of a planar , one wavelength long , microcavity . in the cavity plane , they are localized by disorder , which is due to randomness in the position , composition and sizes of quantum dots located in the anti - node of the cavity . we extend the theory of disorder induced strong localization of electron states to optical modes and obtain quantitative agreement with the main experimental observations .    ,    ,