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: edge localised modes ( elms ) are bursty , quasi - periodic expulsions of energy and particles from the plasma edge in the high - confinement state of toroidal magnetised plasmas @xcite . elms originate from the steep pressure gradient pedestal in the outer closed field line region of a tokamak . a substantial release of magnetic and potential energy , caused by a rapid rise of the rate of momentum transfer through the flux surfaces , leads to enhanced transport and loss of heat and particles into the open scrape - off layer ( sol ) field line region . the filamented peak heat fluxes during elms on the bounding divertor plates present serious restrictions on the performance of future fusion experiments like iter @xcite and necessitate the development of suitable control techniques @xcite . elms generically have been observed in all divertor tokamaks since the initial discovery of the high - confinement `` h - mode '' state in the early 1980s @xcite , while characteristics like frequency and intensity are found to vary widely depending on experimental conditions @xcite . stellarator experiments also are reported to show similar bursty quasi - periodic edge activity in the presence of edge transport barriers @xcite . indication that equilibrium drifts are part of the process has been found by experimental diagnosis of an asymmetry of the energy and particle fluxes on each of the divertor plates and evidence that this is sensitive to the directions of the toroidal current and magnetic field @xcite . phenomenologically , edge localised mode events in tokamaks have been compared to solar flare eruptions @xcite , and , in the cataclysmic variability of observed radiation emanating from the plasma , they may seem to bear also some remote resemblance to the outbursts of pressure - driven dwarf novae @xcite , although the specific instability mechanisms behind these phenomena are clearly of a substantially different nature . on the other hand , the nature of elms has also not yet been completely clarified , in the sense that there is no first - principles based predictive theory available . a number of characteristically different variations of elms have been observed in experiments ( e.g. type i , ii or iii elms ) which may actually be caused by different instability mechanisms @xcite . `` type iii '' elm observations show a remarkable similarity to global bursts found in computations of drift - wave turbulence , where the turbulence generated flows and fluxes are closely tied to the self - consistent equilibrium evolution in the plasma edge @xcite . `` type i '' elms in the experiment specifically show stronger magnetic activity during the burst and are usually associated with the onset of a magneto - hydrodynamic ( mhd ) ideal ballooning mode ( ibm ) instability when a threshold pressure gradient is reached in the steep edge pedestal of an h - mode plasma @xcite . the possible role of an additional current driven instability ( `` peeling - ballooning mode '' ) for `` type i '' elms has also received recent interest @xcite . the ideal ballooning instability with nonlinear phases or aftermath has provided a paradigm for recent and current study of the elm phenomenon theoretically . explosive mhd instability and associated critical phenomena have been advanced analytically @xcite . early braginskii fluid simulations were given in support of this @xcite . large , nonlinear mhd codes have been studying this in more detail , resolving several numerical problems and paying specific attention to the details of the magnetic geometry , with mixed results on the phenomenology of the instability beyond its linear phase @xcite . whether or not the ideal ballooning instability or its peeling / ballooning variant is actually responsible for elms occuring in experiments ( here noting the absence of a demonstrably well - resolved l - to - h mode transition in computations from first principles ) , the nonlinear phase of this `` ideal mhd blowout '' or `` ibm blowout '' phenomenon is of physical interest . the ability of the instability to destroy the edge pedestal layer is undisputed , and the term `` blowout '' is germane . typical physical parameters lead to a situation wherein the ideal mhd interchange growth rate is comparable to the parallel alfvn transit time , one and the same with the basic ideal ballooning criterion . this is fast with respect to mhd but not to microturbulence . provided they can simultaneously treat global mhd ( i.e. , the spectrum covers both global and ion gyroradius scales ) , edge - turbulence computations can also address the phenomenon and are even in a position to treat parts of it which fall outside the paradigm of the mhd model . here we present nonlinear gyrofluid computations of ibm blowout events localised to the edge / sol region of a tokamak plasma , resolved to below the ion gyroradius @xmath0 . as in the large mhd models the self - consistent evolution of the equilibrium is not only included but is also an integral part of the dependent variables . an electromagnetic 6-moment gyro - fluid model ( `` gemr '' ) for both electrons and ions is used , including an energy - conserving treatment of finite ion larmor radius effects and higher - moment terms which occur in the toroidal drift @xcite . the energetic consistency is a key feature since the reaction of parts of the system such as flows and currents which have low energy content can play a central role in the indirect nature of nonlinear dynamics by serving as transfer channels . although the gyrofluid model has a different formulation of nonlinear polarisation ( the generalisation of `` vorticity '' to a two - fluid setting ) than a fluid or mhd model , full correspondence in the regime of validity of the latter models has been shown @xcite . in this sense the model is a superset of global reduced mhd and electromagnetic microturbulence . the global dynamics is in the reduced mhd regime @xcite due to low absolute edge beta values ( @xmath1 ) but the small scales depend on treating the ion gyroradius to arbitrary order . the model does not assume instabilities to occur at any particular scale , but when they occur the resulting spread of the spectrum is found to reach @xmath0 within a few eddy turnover times of the onset of nonlinearity in every case computed . the work reported herein addressed the following points . the relationship of the mhd dynamics to the ideal threshold is much discussed but in the gyrofluid model several mechanisms of microinstability are also present . the threshold issue is obscured by the existence of ion temperature gradient ( itg ) drift instabilities at all beta values , with a smooth transition between them . nevertheless , varying character in the saturation phases is observed at different values . we also address the issue of possible `` rho - star '' dependence ( here the local rho - star is given by @xmath2 where @xmath3 is @xmath0 evaluated at @xmath4 and @xmath5 is the scale length of the @xmath6 profile , both evaluated at the mid - pedestal location in minor radius ) which would indicate a role for diamagnetic drifts . the gradient is sufficiently steep that @xmath7 so that none of these effects can be ordered small . we examine the necessary resolution to obtain converged cases , which gives the spectral range actively involved in nonlinear saturation . the main instability at experimentally relevant beta is the ideal ballooning instability but the main saturation process is energy exchange with broadband electromagnetic drift wave / itg turbulence which the instability itself generates . due mostly to this , the results turn out to be outside of both the magneto - hydrodynamic ( mhd ) and collisional braginskii paradigms , on which most previous approaches have been based . on the other hand , the gemr model still uses delta - f equations , so the actual profile phenomenology in the sol region is not well represented , and as in other work with gemr we concentrate on the properties of the nonlinear dynamics @xcite . following sections give the details of the gemr model , a discussion of the physical difference between itg and mhd dynamics and their relative roles , the methods used to pre - set the axisymmetric equilibrium state so that the computations are done in the absence of axisymmetric oscillations whose decay times can be slower than the ibm rise time , the computations of the blowout phase itself , and the possibility and limitations of quantitative comparison to experiment . the gyrofluid model used herein is given in ref . it is based on the original derivations given in refs.@xcite , with corrections given as motivated by free energy conservation . an alternative derivation using the conservation laws of the underlying gyrokinetic model is given in ref.@xcite . the equations are normalised to spatial scale @xmath8 and time scale @xmath9 where @xmath8 is the minor radius and @xmath10 is the reference sound speed given by @xmath11 , that is , in terms of the electron temperature ( in energy units ) and ion mass . the main parameters are the drift parameter , electron dynamical beta , and electron collisionality , respectively given by @xmath12 where @xmath13 and @xmath14 are the speed of light and fundamental charge , @xmath15 and @xmath16 are the pressure and braginskii collision time @xcite of the electrons , and @xmath17 is the drift scale . note that if @xmath4 then @xmath18 . the reference for @xmath19 is the deuterium mass @xmath20 . for the electrons the physical value of the mass @xmath21 is always used . with global geometry there is no single magnetic field line connection length , but one with a profile with @xmath22 a function of the minor radius coordinate , where @xmath23 is a constant giving the reference value of the major radius . each species @xmath24 has its own set of gyrofluid moment variable equations , one each for density @xmath25 , parallel velocity @xmath26 , parallel and perpendicular temperature @xmath27 and @xmath28 , and parallel and perpendicular energy components of the parallel heat flux @xmath29 and @xmath30 . they are coupled to the electrostatic and parallel magnetic potentials @xmath31 and @xmath32 through self consistent field equations . the tilde denotes a dependent variable to distinguish from constant parameters . the correspondence of these equations to the braginskii fluid model is given in ref . each species is characterised by a background charge density , temperature / charge ratio , and mass / charge ratio , given by the normalised parameters @xmath33 for electrons @xmath34 and @xmath35 . for the main ions @xmath36 while @xmath37 is kept as a parameter . for the trace ions @xmath38 is always taken . for hot trace ions @xmath39 is large while @xmath40 is still moderate . the species gyroradius @xmath41 is given by @xmath42 , which is always small for electrons , close to unity for main ions , and can be moderate or large for trace ions . for this study a simplified geometry is taken with @xmath43 except in curvature terms , so magnetic pumping and shaping effects @xcite are neglected . since dynamics at the electron gyroradius scale are neglected the gyroaveraging effects on @xmath32 are not treated . the moment equations are given by @xmath44 @xmath45 @xmath46 @xmath47 @xmath48 @xmath49 the linearised pressures are given by @xmath50 collisional dissipation is controlled by a collision parameter @xmath51 analogous to @xmath52 and a set of numerical constants for each species , @xmath53 which are the thermal force , thermal conductivity , and viscosity coefficients involved in parallel dissipation , with values given by the braginskii model @xcite , though that regime is never reached in core turbulence . for electrons these coefficients are @xmath54 and @xmath55 and @xmath56 , and for ions they are @xmath57 and @xmath58 and @xmath59 , respectively . the resistive dissipation is given by @xmath60 with parallel current given by @xmath61 noting that it is the electron heat fluxes that enter @xmath62 and as with @xmath63 the @xmath62 term enters the parallel motion in the same way for every species . the anisotropy dissipation is given by @xmath64 the thermal conduction components are given by @xmath65 @xmath66 @xmath67 the landau damping effects are modeled by @xmath68 with thermal velocity @xmath69 given by @xmath70 . the gyroaveraging is done through pad approximants to operators which would act in fourier space , @xmath71 and then the gyroaveraged potentials are @xmath72 noting that they are species dependent . the main exb advection and the flr correction pieces are given by @xmath73 - \nu_\perp\dpp^4 f + \nu_\parallel\ddpl f\ ] ] @xmath74\ ] ] where the @xmath75 $ ] symbols denote the nonlinear bracket defined below , and @xmath76 denote the artificial dissipation coefficients . the parallel derivative is given by @xmath77\ ] ] the curvature operator is given by @xmath78\ ] ]      the species are coupled through two equations which give the self consistent response of the field potentials to the evolution of the gyrofluid moment variables . the electrostatic potential is governed by quasineutrality , which sets the total space charge density to zero . the density for each species is given by a part due to the gyrocenters and another due to polarisation , which is what sets @xmath79 . this polarisation equation is given by @xmath80 where @xmath81 gives the gyroscreening of the main ions . in this version of the model the electron and trace ion contributions to gyroscreening are neglected , due to small @xmath82 and zero @xmath83 , respectively . on the right hand side the electrons and main ions enter with oppositely signed @xmath83 values and again here the trace ions do not enter . this establishes them as a trace species .    in a similar way the trace ions are left out of the induction equation due to the zero @xmath83 , so that @xmath84 determines @xmath32 noting the way that the normalisation scales it with respect to @xmath85 . more detail , including energy conservation and the relationship of these equations to it , is given in refs . @xcite .      in these expressions the operators @xmath86 and @xmath87 and the nonlinear brackets are determined by the representation of the geometry that is used .    the magnetic geometry is a simplified representation of an axisymmetric magnetic field using field aligned hamada coordinates . the field representation is @xmath88 where @xmath89 is the physical toroidal angle about the symmetry axis , @xmath90 is the magnetic flux function , and @xmath91 is given by the constraint of mhd equilibrium . without loss of generality one may define magnetic flux coordinates @xmath92 with the following properties , @xmath93 where @xmath94 is the volume enclosed by the surface with flux @xmath90 , then @xmath95 and @xmath96 are two further flux functions found , and the subscript denotes the partial derivative with respect to @xmath97 . the ratio @xmath98 is another flux function which gives the pitch of the magnetic field . the angles are cyclic on @xmath99 $ ] which determines the coordinate jacobian to be unity . hence the metric tensor @xmath100 satisfies @xmath101 . the function @xmath102 is found by the constraint that @xmath103 is cyclic on the unit interval @xmath104 where @xmath105 is a simple cyclic coordinate on the flux surface in the poloidal plane ( e.g. , path length on the curve , or the physical angle position about the magnetic axis with respect to any branch cut ) . then @xmath106 is found @xmath107 where the angle brackets give the flux surface average defined as @xmath108 now , @xmath109 is also defined and the toroidal angle coordinate is given by @xmath110 which completes the prescription .    the field aligning is a one - to - one and onto coordinate transform from @xmath92 to @xmath111 given by @xmath112 where @xmath8 is a reference minor radius and @xmath113 is an arbitrary function of @xmath114 which is chosen to make the off diagonal perpendicular metric element @xmath115 vanish at a particular location . this is called the shifted metric procedure @xcite . the point is that at any position in @xmath116 where perpendicular drifts or @xmath86 is evaluated , the coordinate elements are rectangular , and the information contained in magnetic shear enters only in derivatives in the parallel coordinate . the magnetic field components now satisfy @xmath117 at all locations in space for any choice of @xmath113 . hence although the perpendicular coordinate elements are rectangular only for one particular location in @xmath116 the field aligning is exactly satisfied everywhere . the coordinate metric elements @xmath118 are found in the @xmath92 representation . then the grid locations @xmath119 are all given their own members of the family of these field aligned coordinates via @xmath120 then we have @xmath121 at @xmath119 . the drift tensor elements are given by @xmath122 that is , @xmath123 is the only independent , nonvanishing component . we have @xmath124 which defines the bracket @xmath125 $ ] for any scalar field quantities @xmath79 and @xmath126 . with @xmath79 the electrostatic potential this gives the exb advection term , for electrons for which flr corrections are neglected by taking @xmath127 . for ions this is generalised by the flr corrections as given above .    for local equations as those used here ( nonlinearities kept only in advection and magnetic flutter ) the geometry model must be modified to retain free energy conservation by the equations . in particular , any function of @xmath114 such as a radially varying normalised parameter ( such as the temperature / charge ration @xmath39 ) multiplying any of the curvature terms will cause them to fail to conserve free energy . the compressibility is already split away from advection and placed into the curvature terms , so the brackets must retain the properties of incompressible advection . this means that the model for @xmath128 which is used must satisfy @xmath129 and to avoid any confusion the curvature terms are written with @xmath130 , which can be given any spatial dependence . correspondingly , the mhd version of the continuity equation ( neglecting diamagnetic fluxes ) has already been manipulated according to @xmath131 @xmath132 and then terms such as @xmath133 are neglected due to the ordering . the result is the same as the local continuity equation @xmath134 where @xmath135 is a normalising constant , if we identify the dependent variable @xmath136 with @xmath137 . it is important to note that with the self consistent background being evolved this also includes the profile : @xmath138 . the simplified geometry is now defined by its operators . the laplacian and gyroaverging operators neglect @xmath87 . the form in the ampere s law ( which does not involve the gyroradius ) is given by @xmath139 the form in the gyroaverging operations and polarisation is given by @xmath140 the brackets are defined by @xmath141=q[f , g]_{xy}+[f , g]_{xs}-q_{,x}(s - s_k)[f , g]_{ys}\ ] ] for any two scalar fields @xmath142 with each bracket piece defined by @xmath141_{ij } = \lp{\partial f\over\partial x^i}{\partial g\over\partial x^j }          - { \partial g\over\partial x^i}{\partial f\over\partial x^j}\rp\ ] ] with the third piece not used since these are always evaluated at @xmath119 . nonlinear computations on ideal or peeling ballooning mode elm scenarios in tokamaks have previously been based on magnetohydrodynamic ( mhd ) models @xcite , and the respective codes had largely been tailored towards mode structures expected by linear analysis . restriction to single - fluid mhd equations including two - fluid correction terms allows computations of large to mesoscale dynamics in realistically shaped 2d tokamak equilibria including x - point geometry crossing the separatrix , due to the lack of a need to use field - aligned coordinates . recent mhd simulations are able to reproduce spatio - temporal mode structures , divertor footprints and other characteristics with good agreement to experimental diagnostics , in particular for the growth and immediate blowout phase of the finger - like instability @xcite . however , as soon as the dynamics becomes nonlinear , the spectrum broadens to include scales normally associated with drift wave dynamics . this turbulence is basically of a drift - alfvn type with strong edge itg character ( see ref . @xcite for the signatures of the various mode structure types ) , made more powerful than otherwise by the energetic access to the long - wavelength mhd component . the itg character results from the steeper logarithmic gradients in both temperature profiles , and the lack of involvement of ion temperature fluctuations in the parallel nonadiabatic response of the electrons . although edge turbulence of this type ( cf . @xcite ) is strongly suppressed in the h - mode phase and is not initially involved in the instability phase , the experience of edge turbulence becomes relevant as these physical components become involved in all the nonlinear phases of the blowout . the saturation and aftermath of the blowout , which in an experiment carries a large part of the actual transport losses that lead to a degradation of the pedestal , should be expected to involve physics not contained in the mhd model . ion temperature gradient ( itg ) driven modes and ideal ballooning modes are both essentially caused by the gradient and curvature driven interchange instability and show similar character in their initial linear growth phases . instability is achieved above their respective critical gradients , determined by the ratio @xmath143 between density gradient to ion temperature gradient scale lengths for itg , and by the ideal mhd ballooning parameter @xmath144 for ibm , where @xmath145 is the magnetic shear parameter and @xmath146 is a critical threshold which depends on beta , collisionality , and toroidicity . the principal difference between these modes is the relative role of the parallel alfvn responses which tend towards enforcement of an adiabatic response in the electrons ( parallel force balance , with the electron pressure gradient ) . this is nonexistent in an mhd model and subdominant generally for an mhd instabililty . but it constrains the electrons if the @xmath147 parameter is below threshold , still allowing the itg instability because the adiabatic response does not involve the ions . the itg instability exists at all @xmath147 values , and furthermore for the edge situation the itg instability and mode structure transitions smoothly to a drift wave one in a nonlinear setting for @xmath148 . hence there is actually no finite threshold in values of @xmath149 or @xmath147 for the nonlinear edge situation . therefore , the existence of a threshold in the experiment is the same as the mechanism which maintains the h - mode , which is not yet well understood . it follows that a gyrofluid model computation of an elm scenario with enough resolution to allow for ion - gyroradius based dynamics has to face the lack of a simple linear threshold in the parameter space . the transition from initial ( micro-)instabilities to generic edge turbulence was studied in detail in ref . as the most unstable linear modes crystallise out of an initial random bath of small - amplitude perturbations , the linear growth rate rises and becomes steady . the maximum value of the instantaneous growth rate of total fluctuation free energy @xmath150 given by @xmath151 may be taken as the maximal linear growth rate . the curve of @xmath152 then falls very sharply to zero ( over about @xmath153 ) as saturation occurs . there is some structural adjustment over the next few @xmath154 as the spectrum fills out , and then the turbulence is fully developed . but over the adjustment phase the value of @xmath155 is well below its previous maximum . we will use the same diagnostic herein , except that the total ion exb heat flux @xmath156 averaged over the computational domain is used instead of @xmath150 because most of @xmath150 is represented by the self consistent profiles while @xmath156 is entirely due to fluctuations . visualisation of the tokamak elm burst : the edge region of a circular poloidal tokamak cross - section is shown on the left with computational perpendicular @xmath116-domains ( @xmath157 2 , 5 , 8 , 11 , 14 of @xmath158 total ) mapped onto the circle during the maximum growth phase ( @xmath159 in units of @xmath9 ) of an ideal ballooning mode . on the right , the outer midplane area ( @xmath160 with size @xmath161 and @xmath162 ) is magnified and the radial domain stretched by a factor of 2 for better visibility at times @xmath159 ( max . flux ) and @xmath163 ( turbulent aftermath).,width=188 ]    the ibm elm blowout scenario is similar to this , initially , except that the instability is not a microinstability . nevertheless , the scale differs by less than an order of magnitude : the toroidal mode numbers for linear itg instabilities are in the range of @xmath164 while the main ideal ballooning mode is near mode number @xmath165 , on the entire flux surface , for these typical @xmath166 values ( recall @xmath167 for @xmath168 , and the ratio @xmath169 is not far from unity in experimental cases ) . the itg mode numbers are determined by the ion gyroradius and the dimensionless scale ratio @xmath170 , while the mhd values are determined by the width of the pedestal . since @xmath5 is not larger than about @xmath171 in an h - mode pedestal , these scales are not disparate . due to the closeness of the native scales , the ibm instability very quickly transfers energy to the ion gyroradius scale , which is only about two or three steps away in a turbulence cascade which acts at a factor of two at each step . this brings the itg and drift alfvn physics into play , and in any nonlinear stages the latter involves stable shear alfvn component with substantial energy content . the moment of saturation is defined as the time at which the ion exb heat flux averaged over the computational domain is maximum ; we refer to this as the `` peak - flux '' time . the results shown below indicate that coincides with the establishment of the fully turbulent regime . representation of this phase requires resolution of all scales and inclusion of the appropriate drift wave physics in the model . how much resolution is actually needed depends on the strength of the blowout , which is situation dependent . we therefore include a resolution scan on the toroidal mode number spectrum . the result is that convergence is reached in the temporal phases just after peak flux only if the spectrum reaches down to @xmath172 . the self - consistent equilibrium coupling in our gyrofluid model only allows treatment of shifted circular @xmath173 geometry ( the difference to mhd codes is that the high - resolution flr - gyrofluid dynamics necessitates field - aligned coordinates ) . effects of flux surface shaping are postponed to later studies . our focus here is on basic physics issues regarding the ability of the model to capture mhd phenomena ( see also ref . @xcite concerning global alfvn oscillations with self consistent profiles ) , resolution , turbulent character , scalings and tests which are also accessible by @xmath173 geometry . in particular , the theory of explosive instability has been formulated for circular geometry @xcite and is therefore accessible by the present model . first - principles based local drift wave edge turbulence simulations are not able to obtain a realistic h - mode edge state with the known experimental characteristics : correct density and temperature pedestal profiles shapes or strength of flow shear are not obtained by self - consistent evolution by specifying core sources only , nor has a threshold transition character been found in any verified edge turbulence simulation @xcite . therefore some kind of `` modelling '' has to take place when the ibm instability ( as an h - mode phenomenon ) and its subsequent nonlinear evolution is simulated with a nonlinear gyrofluid turbulence code : although the realistic development of an edge transport barrier ( and thus a full elm cycle ) can not be directly obtained , one still may prescribe the h - mode pedestal profile before the onset of an elm , known from experimental data , as an initial state for the simulation . as a base case for the prescribed pedestal profiles the well diagnosed edge characteristics of asdex upgrade h - mode shot # 17151 is used here @xcite . the local parameters , taken as mid pedestal values , correspond to electron and ion temperatures @xmath174  ev , @xmath175  ev , densities @xmath176  m@xmath177 , magnetic field strength @xmath178  t , major torus radius @xmath179  m , aspect ratio @xmath180 , perpendicular temperature gradient length @xmath181  cm , density gradient length @xmath182  cm . the profile of the safety factor @xmath183 is parabolic yielding local values at the last closed flux surface ( lcfs , @xmath184 ) of @xmath185 and @xmath186 . the radial domain of the simulations covers a range of @xmath187 on either side of the lcfs . the nominal pressure values are @xmath188 and @xmath189 . all times are given in units of @xmath9 . the computational grid is given in terms of the @xmath111 coordinate domains . the spacing is equidistant . for each grid point in the third coordinate @xmath119 the @xmath190-coordinate is @xmath191 as defined in eqs.([eqcoords],[eqalphak ] ) so that for each grid point @xmath192 the coordinate system @xmath193 is the one which has @xmath194 at @xmath119 . for the nominal case @xmath162 perpendicular and @xmath195 parallel mesh points are used . the radial domain with @xmath196 spans the plasma edge region between the h - mode pedestal top , with plasma core parameters as inner boundary values , and the outer bounded scrape - off layer region ( @xmath197 ) . this represents a spatial range from the global scale to smaller than the ion gyroradius scale , with a ratio @xmath198 between ion gyroradius @xmath199 and minor torus radius @xmath8 . in the sol region ( @xmath200 ) the parallel boundary condition is replaced by a debye sheath model whose treatment is given in ref . @xcite . radial tokamak edge pedestal profiles of the plasma pressure before and after an elm blowout event across the separatrix into the scrape - off layer ( sol ) . the area shaded in grey shows the present radial numerical simulation domain.,width=283 ] the initial conditions are thus prescribed and are based on experimentally diagnosed radial temperature and density pedestal profiles @xmath201 and @xmath202 for each species ( electrons and ions ) . a consistent electrostatic potential @xmath203 is derived by numerically solving the neoclassical balances ( parallel dynamics , toroidal drifts , and collisional dissipation , but not transport ) in a pre - processing step with a modified ( zonally frozen ) gemr setup , resulting in a time - steady 2d dissipative solution . the @xmath173 mhd equilibrium in turn is set internally in gemr by computing the shafranov shift , and in the present setup would not be consistently described by prescribed in terms of external or coupled shaped equilibrium solvers @xcite . the parallel and perpendicular electron and ion temperatures , @xmath204 , @xmath205 , @xmath206 and @xmath207 , are directly adopted and fixed from experimentally derived values by filtering the zonal component out of the total time derivative @xmath208 . their values are zonally frozen by setting @xmath209 where @xmath210 is the zonal ( parallel and perpendicular ) average of all right - hand - side terms @xmath211 in the temperature evolution ( eqs . 5,6 ) including numerical dissipation .    on the other hand , gyrocenter densities @xmath212 and @xmath213 have to be set to obey relaxation relations that allow the vorticity to freely evolve into equilibrium . this is achieved by freezing the zonal component of the sum @xmath214 , where @xmath215 , as the density part of the pressure during the equilibration phase , but allowing the difference ( i.e. , vorticity ) to evolve freely . this ensures that the contribution of densities to the total pressure is zonally frozen through the relation @xmath216 while the densities relax regarding to @xmath217 and @xmath218    the numerical solution of the equilibration phase , starting directly from realistically steep pedestal profiles @xmath219 , into steady state is delayed by long , weakly damped global geodesic alfvn oscillations . convergence is expedited by ramping up all of the gradients gradually from zero to prescribed value over the first @xmath220 of the run by @xmath221    this pre - processing equilibration phase is run until convergence with reduced perpendicular resolution ( @xmath222)=(@xmath223 ) , and without the exb and magnetic flutter nonlinearities , which allows establishment of the 2d structure in a smooth manner . then , the resolution is increased to the nominal values ( @xmath224 ) , and a random turbulent bath with relative amplitude @xmath225 is added to the background pedestal profiles inside the closed flux surface region . this procedure reduces transient alfvnic and geodesic acoustic ringing and prepares a reproducible initial state . however , the following sudden release of the nonlinearities also leads to transient oscillations . depending on parameters , these may still be present at the onset of the instabilities , and can obscure a clean view on the nonlinear growth rates , which will be relevant for the discussion below on diagnosing linear or explosive instability .    a schematic sketch of the equilibrium pedestal profiles , representing an idealised asdex upgrade h - mode scenario @xcite , is shown in fig . [ f : profiles ] together with the final relaxed state after the elm blowout phase ( which is discussed in the next section ) . when this initialised pedestal pressure profile is ideal ballooning unstable , the ibm instability in gemr simulations is observed to be linearly growing in the pedestal region and at the onset of nonlinearity further overshoots and saturates , representing the beginning of the turbulent blowout phase during which a substantial fraction of the pedestal energy is thrown onto the sol . subsequent evolution involves turbulence in both the pedestal and sol regions as the original pedestal energy is dissipated . previous nonlinear approaches on elm ideal ballooning mode burst computations have treated only the initial growth and nonlinear phases , focused on low - wavelength modes and resolved only the mhd - relevant scales , excluding treatment or discussion of the ion gyroradius scales . this however precludes the development of fully developed microturbulence which is caused by the onset of nonlinearity ( robust transfer of free energy to smaller scales ) . ultimately , this microturbulence decides both the transition to nonlinearity and the eventual saturation ; that is , most of the post - peak transport curve . this affects the mhd scenario of nonlinear explosive growth which most previous work advances . however , we find that the exclusion of the ion gyroradius scales produces an under - resolved situation . we measure the quantitative character of the growth in both linear and nonlinear phases with a growth curve @xmath226 defined in terms of the heat flux @xmath156 instead of the fluctuation free energy as explained above . the heat flux is computed as a zonal ( flux surface ) average , @xmath227 @xmath228 and then the time trace @xmath229 , normalised in terms of @xmath230 at nominal parameters , is computed as a volume average over the central half of the radial domain .    the ibm instability , which follows the random seeding of the pre - processed pedestal profile in our computations , is very violent , growing for the nominal case at a rate @xmath231 , just below the ideal interchange rate . on the other hand the growth curve @xmath152 appears qualitatively the same as in basic turbulence cases . at all time points in the nonlinear phase @xmath152 is well below its previous maximum . just after initial saturation there is some nonlinear evolution in which @xmath152 crosses zero before settling down into long - term statistical saturation . at late times the initial blowout no longer imprints the results : with a fixed source one merely finds bursty turbulence thereafter . hence , there is no evidence for explosive instability in this nominal case , which , if present , would be visible in a finite time singularity in the fluctuating free energy @xmath232 rising with a power @xmath233 .    the burst and the resulting decaying turbulence act to transport the pedestal plasma across the lcfs , where energy and particles are lost within short times by parallel boundary outflow to the scrape - off layer ( sol ) limiter . without maintenance of a heat and particle source at the core boundary , this leads to decay of the initial gradient with a rapid onset at the initial transition to nonlinearity at about @xmath234 . this time point also corresponds to the peak - flux time . total energy @xmath235 ( top figure ) and ion heat flux @xmath229 ( bottom figure ) during the elm blowout : without source maintenance ( dashed curve ) ; with the source term switched off at the time of maximum instability ( here : at @xmath236 ) ; and with the source maintained for the whole simulation time.,title="fig:",width=283 ]   total energy @xmath235 ( top figure ) and ion heat flux @xmath229 ( bottom figure ) during the elm blowout : without source maintenance ( dashed curve ) ; with the source term switched off at the time of maximum instability ( here : at @xmath236 ) ; and with the source maintained for the whole simulation time.,title="fig:",width=283 ]    however , the ion temperature profile already starts decaying during the low - amplitude linear growth phase , due to neoclassical ( finite larmor radius ) transport ( for detail on how this works in a fluid model see ref.@xcite ) . this situation is shown in fig . [ f : energy ] ( left frame , dashed line ) , where the total pedestal energy is reduced by a quarter of its initial content before the peak flux phase . therefore we test sensitivity against sources . maintenance of the ion temperature profile by a fixed ion heat source localised to the inner boundary ( following the method for driven cases in ref . @xcite ) ensures nearly unchanged gradients during the linear growth phase . at the end of the simulation time ( @xmath237 ) after the elm burst , the energy content is reduced to around 10% of the initial content when the source is absent or is switched off at the peak flux time ( @xmath234 ) , and to around 20% when the source is maintained for the whole simulation time . the influence of either switching off the source or maintaining the zonal profiles on the nonlinear ion heat flux @xmath229 is shown in fig . [ f : energy ] ( right ) : the linear growth phase and the peak flux is not significantly changed by the source . as expected , the heat flux saturates for late times with maintained core boundary inflow and is otherwise decaying .     consistency check for different perpendicular resolutions by varying the grid point number @xmath238 between 64 and 512 . convergence is seen for @xmath239 and higher.,width=283 ]    consistency checks on the perpendicular resolution are made by varying @xmath238 between 64 , 128 , 256 and 512 . the evolution of @xmath229 for these cases is shown in fig . [ f : nyres ] . the @xmath240 case is clearly under - resolved and overestimates the total energy decay by a factor 2 . for @xmath241 the initial part of the post - peak decay phase disagrees with the @xmath162 case . convergence is found for @xmath239 and higher . the nonlinear saturation phase thus requires resolving the ion gyroradius scale . this is the clearest indication that turbulence for toroidal mode numbers beyond about @xmath242 extending down to the @xmath243 range is involved in the blowout saturation process . perpendicular mode number @xmath244 spectra for @xmath159 ( top ) and @xmath163 ( bottom ) around the peak flux phase . the vorticity spectrum ( bold lines ) is already flattening down to the ion gyroradius scale for @xmath163 only a few @xmath245 after the maximum linear growth phase.,title="fig:",width=283 ]   perpendicular mode number @xmath244 spectra for @xmath159 ( top ) and @xmath163 ( bottom ) around the peak flux phase . the vorticity spectrum ( bold lines ) is already flattening down to the ion gyroradius scale for @xmath163 only a few @xmath245 after the maximum linear growth phase.,title="fig:",width=283 ]    the transition from linear instability to turbulence is studied during the peak flux phase around @xmath234 . [ f : kyspecs ] shows perpendicular mode number @xmath244 spectra of the squared amplitude of various fluctuating plasma quantities ( density @xmath246 , ion temperature @xmath247 , electrostatic potential @xmath79 and vorticity @xmath248 ) for @xmath159 and @xmath163 . this time difference corresponds to @xmath249 in physical units and around @xmath250 in local drift units , which is only slightly faster than the overshoot and saturation times known from edge microinstability cases @xcite . initially , the ion temperature gradient ( itg ) driven microinstability and the ideal ballooning mode ( ibm ) compete in growth out of the random low - amplitude bath . the itg mode grows strongest near the separatrix due to radially local steepening by parallel sol diffusion , which in our simulations may be overestimated by using the standard fluid bohm outflow boundary conditions . for @xmath251 the linear ibm is clearly dominant near a toroidal mode number of 9 - 10 for the nominal aug parameters , consistent with experimental observations @xcite . around @xmath252 rapid formation of a turbulent cascade range in the spectra is observed and the vorticity spectrum is already spread out to the ion gyroradius scale ( @xmath253 ) . this is a manifestation of the role of self generated drift wave turbulence in the saturation process . because of the way both the diffusive mixing and vorticity scattering nonlinearities enter the drift wave physics @xcite the spectrum is held together as a unit , all scales down to @xmath0 are involved in the saturation phase of the overall elm blowout transport . the involvement of @xmath0 at unit order , in turn , necessitates a nonlinear gyrofluid model ( and at some later date , a nonlinear gyrokinetic model ) . the characteristics of the blowout at various plasma beta is studied for @xmath254 , with the absolute @xmath255 a factor of @xmath256 higher due to the @xmath257 value . a competition between unstable growth of itg high - n modes and ibm low - n modes is observed : for low @xmath255 an early onset of itg turbulence degenerates the profile and influences the strength of the ibm . for large @xmath255 the nonlinear saturation is acting faster and stronger . the two effects lead to a rather nondistinctive result concerning the peak elm flux as shown in fig .  [ f : beta ] , with a continuous transition between itg and ibm triggered onset of the turbulent aftermath . competing growth of itg , microtearing and ibm instabilities determine onset times and intensities of the blowout time traces @xmath229 when scaling with the local plasma beta : for @xmath258 only an itg micro - instability is growing around @xmath259 and saturating on a low h - mode transport level . for @xmath260 the ibm instability dominates and leads to a clean elm signal around @xmath261 . in the nominal case ( @xmath262 ) itg and ibm mhd growth compete initially with similar growth rates but the ideal ballooning mode takes over around @xmath263 and further determines the blowout . for @xmath264 a microtearing instability is saturating on a gyro - bohm transport level around @xmath261 and is transformed into more violent itg mhd turbulence after deterioration of the initial electron temperature gradient ( @xmath265 ) . , width=283 ]    the presence or absence of background current gradient terms , which in this model are set by replacing @xmath266 everywhere the electron parallel velocity @xmath267 appears in the equations , with @xmath268 given by the @xmath106 profile , was found to have no discernible effect on the result . indeed self - consistent inclusion of @xmath269 effects with possible role in the l - to - h transition was among the motivations of the gemr model in the first place @xcite . however , the magnitude of even an impressive pedestal current peak is only one to three times the nominal saturation current given by @xmath270 , while values closer to @xmath271 are required to enter the energetics effectively . since @xmath272 for the nominal case , current gradient effects are very weak . the main conclusion of this study is that the qualitative nature of the saturation and aftermath of the initial ibm blowout is the same as for generic edge turbulence given a small - amplitude start . the transition from linear mode structure and energetics to turbulence found for these blowout cases is the same as in ref . only the nature of the linear mode itself differs . the blowout saturates upon its own self generated drift - alfvn turbulence , with a strong ion temperature component given the gradients . the vorticity spectrum reaches quickly to the ion gyroradius ( @xmath0 ) scale , requiring the gyrofluid model and explaining the numerical difficulties seen with braginskii models  it can be argued that the cases given in ref . @xcite crash on entry to the nonlinear stage . convergence in the aftermath requires resolving at least @xmath0 . unfortunately , due to the lack of a self consistent h - mode state in a well resolved computation , no threshold is found . at lower beta values one simply finds generic edge turbulence driven by the temperature gradients . it is not clear that this scenario really describes an actual elm , although the energetic growth and decay curves are not unrealistic . however , with a well resolved transition to nonlinearity in both the energetic peak and aftermath phases , our studies find no evidence of an explosive mhd phenomenology . indeed it can be argued that nonlinear mhd processes are pre - empted by the efficient transfer to smaller scales through the two - fluid drift wave physics more commonly associated with microturbulence . this work was supported by the austrian science fund fwf under contract y398 , by a junior research grant ( `` nachwuchsfrderung '' ) from university of innsbruck , and by the european communities under the contract of associations between euratom and the austrian academy of sciences , carried out within the framework of the european fusion development agreement . the views and opinions herein do not necessarily reflect those of the european commission . h. strauss , l. sugiyama l. , c.s . chang , et al . 21st iaea fusion energy conf . ( chengdu , china , 2006 ) ( vienna : iaea ) cd - rom file th / p8 - 6 and http://www-naweb.iaea.org/napc/physics/fec/fec2006/html/index.htm
Abstract:

Answer: three dimensional electromagnetic gyrofluid simulations of the ideal ballooning mode blowout scenario for tokamak edge localized modes ( elms ) are presented . special emphasis is placed on energetic diagnosis , examining changes in the growth rate in the linear , overshoot , and decay phases . the saturation process is energy transfer to self generated edge turbulence which exhibits an ion temperature gradient ( itg ) mode structure . convergence in the decay phase is found only if the spectrum reaches the ion gyroradius . the equilibrium is a self consistent background whose evolution is taken into account . approximately two thirds of the total energy in the edge layer is liberated in the blowout . parameter dependence with respect to plasma pressure and the ion gyroradius is studied . despite the violent nature of the short - lived process , the transition to nonlinearity is very similar to that found in generic tokamak edge turbulence .


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 simultaneous treatment of disorder and interactions in strongly correlated electron systems has always formed a knotty challenge ; this is because of the dearth of manageable analytical techniques that can deal with disorder and interactions at the same time.@xcite the strongly correlated system of interest in this work is the two - dimensional electron gas ( 2deg ) in a perpendicular magnetic field at landau level filling factor @xmath0 , whose ground state is commonly known as the _ quantum hall ferromagnet_.    the quantum hall ferromagnet is the spin - polarized ground state of the 2deg at @xmath0 in which all electrons completely fill the lowest landau level with spin up polarization . such configuration minimizes the coulomb energy for fermionic systems . in general , it is a competition between kinetic and coulomb energies , which determines the ground state . in the case of the quantum hall ferromagnet having @xmath0 the kinetic energy is frozen and does not change with spin flip , thus , the ground state is ferromagnetic , even with zero zeeman splitting . typically , the zeeman splitting in the gaas heterojunctions turns out to be roughly 70 times smaller than the spacing between landau levels and an order of magnitude smaller than the coulomb energy per particle . the neutral elementary excitations are spin wave excitations , also called magnons . the spin waves can be described by the action of the spin lowering operator @xmath1 , projected to the lowest landau level , on the ferromagnetic ground state . it turns out , that the projected operator creates an exact excited eigenstate of the hamiltonian . in the regime of low momenta , the magnon s dispersion is quadratic and the coefficient of the quadratic term represents a phenomenological constant known as the _ spin stiffness_. the spin stiffness provides a measure of the free - energy increment associated with twisting the direction of the spins . a significant spin stiffness indicates that the system lies in the ferromagnetic phase , while a paramagnetic state corresponds to a vanishing spin stiffness . the spin wave dispersion at very large momenta saturates at a constant value given by the sum of the coulomb and zeeman energies . thus , at large momenta , the value corresponds to the energy of separate quasiparticle and quasihole excitations . one approach that has successfully dealt with strongly correlated electron systems is the so - called _ bosonization _ procedure . bosonization is a non - perturbative approximation scheme that essentially treats the electron - hole excitation , known as _ exciton _ , as a bosonic single - particle ; consequently , a fermionic hamiltonian can be recast into a bosonic one . in 1950 , tomonaga revealed , in a ground - breaking paper,@xcite that the application of the bosonization formalism to a one - dimensional electron gas ( 1deg ) yielded an exactly - solvable hamiltonian . the reason is that the electron and the hole propagate with nearly the same group velocity in the low - energy region . however , that is not the case in two dimensions . at a given momentum @xmath2 , the particle - hole pair excitation holds a continuous range of energies . therefore , it is less straightforward to construct a coherently propagating bosonic entity in two dimensions . the first attempt to extend the bosonization procedure for higher dimensions was done by luther @xcite and then revised by haldane.@xcite castro neto and fradkin,@xcite as well as houghton and marston,@xcite developed a bosonization technique for a fermi liquid in any number of dimensions . as regards the interacting 2deg subject to an external perpendicular magnetic field , westfahl jr . _ _ et al.__@xcite constructed a formalism that treated the elementary neutral excitations of the system , the _ magnons _ , in a bosonic framework such that the fermionic hamiltonian of the system was transmuted into a quadratic bosonic hamiltonian . the drawback is that this method is valid in the limit of weak magnetic fields , which amounts to large landau level filling factors @xmath3 . doretto _ _ et al.__@xcite extended the methodology of westfahl jr . _ _ et al.__@xcite to the case of the 2deg at @xmath0 ( i.e. for a very strong magnetic field ) . given that the system is now restricted to one landau level , the task greatly simplifies , since the landau level quantum degree of freedom can then be disregarded . projecting the original fermionic interacting hamiltonian of the system into the lowest landau level , which is completely filled ( @xmath0 ) , allows one to expand it in magnon states.@xcite it then turns out remarkably that the dispersion relation of the free magnons coincides with the result derived by kallin and halperin@xcite within the fermionic description at the _ random phase approximation _ ( rpa ) level and the quartic interacting part of the magnon hamiltonian might be related to the skyrmion - antiskyrmion neutral excitations of the hall ferromagnet.@xcite moreover , in the vicinity of the ground state , without magnon - magnon interactions , magnons behave like bosons . this allows to treat magnons approximately as bosons in the so - called single - mode approximation . @xcite here , we intend to calculate a quantum phase transition in the quantum hall ferromagnet driven by disorder , accounting for the coulomb interactions between electrons . we will use the bosonization technique allied to the usual self - consistent born - approximation for the disorder averaging procedure .    before presenting the results obtained in this paper , it is worth getting acquainted with the current status of research related to the field . to begin with , green @xcite propounded that the vanishing of the renormalized spin stiffness at a threshold value of the disorder strength signifies the occurrence of a depolarization transition from the ferromagnetic phase to a paramagnetic one . his finding is based upon a previous result established by fogler and shklovskii , @xcite who proffered the same idea in the case of higher landau levels . green established this proposition in the framework of non - linear sigma model , used as an effective low - energy theory in the regime of weak disorder . the other quantity that green computed is the disorder contribution to the optical conductivity , which he found to be unmeasurably small . finally , green established that the quantization of the hall conductivity is not affected by the presence of weak disorder in the system.@xcite    another work was carried out by sinova , macdonald and girvin,@xcite who established the occurrence of a phase transition from the paramagnetic state to the partially - polarized ferromagnetic one and then finally to the fully - polarized ferromagnetic one as the interaction strength increases relative to the disorder strength . they determined this result by computing the average value of the spin polarization as a function of the interaction strength relative to the disorder strength . et al.__@xcite did consider coulomb interactions within the framework of the hartree - fock approximation . moreover , the transition from the paramagnetic phase to the ferromagnetic one was found to take place when the coulomb energy scale is about twice as large as the landau - level - broadening disorder energy scale . as a final point , the authors inferred that no phase transition can take place in the strong disorder limit . the last germane paper was published by rapsch , lee and chalker.@xcite they established the occurrence of a phase transition from the ferromagnetic state to the so - called spin glass phase . this result was obtained by calculating the magnetization , the magnetic susceptibility and the spin stiffness as functions of the disorder strength . they assumed the disorder potential to be gaussian distributed and described the system in terms of a _ semiclassical _ spin model . in their model , they took into account coulomb interactions within the hartree - fock approximation but modelled them as being short - ranged . like green,@xcite rapsch _ _ et al.__@xcite computed the disorder contribution to the optical conductivity and found as well that it is undetectable . finally , they calculated the dielectric susceptibility of both the partially - polarized ferromagnetic phase and the spin glass one and they concluded that both regimes display an insulating behavior at low momenta and a metallic behavior at large momenta . let us now put our work in perspective . our objective is to establish the behavior of the renormalized spin stiffness as a function of the disorder strength in order to ascertain a potential quantum phase transition driven by disorder to a non - ferromagnetic state . indeed , if the spin stiffness vanishes for a critical value of the disorder strength , then this signals an instability in the ferromagnetic phase.@xcite on the other hand , the appearance of an imaginary component of the spin stiffness , which might be interpreted as a spin wave damping,@xcite at a certain disorder strength , might indicate the appearance of localized spin waves and a spin - glass phase transition . another important characteristic is the pauli susceptibility , which diverges at the point of the phase transition from a non - ferromagnetic to a ferromagnetic state , indicating spontaneous magnetization . we consider a fully _ quantum _ model , include a short - range weak disorder potential up to the @xmath4 order born approximation and treat the true _ long - range coulomb _ interactions up to the rpa level . the method that we employ consists of five steps . first of all , a bosonized expression of the total hamiltonian , which includes a contribution from disorder , is sought for . the dispersion relation of the free bosons corresponds to the one computed by macdonald _ _ @xcite and more explicitly by doretto _ et al._,@xcite which entails interactions between electrons up to the rpa level . the second step consists in obtaining the full green s function , and precisely its disorder self - average . in our case , because the impurities are randomly distributed throughout the system , the disorder self - average can also be taken by averaging over the impurity positions . the third stage is then to determine the self - energy of that disorder self - averaged green s function through the use of the dyson s equation . the self - energy is determined in the low - impurity density and weak disorder scattering approximations . as a result , the self - energy corresponds to a single diagram with one propagator line and two disorder potential lines . the propagator line is evaluated within two further possible approximations : the bare approximation , which consists in using the bare bosonic propagator , and the self - consistent approximation , which uses instead the full disorder self - averaged green s function . one must bear in mind that both propagators take into account interactions between electrons up to the rpa level . furthermore , the bare approximation is first taken in the long wavelength limit , which keeps the lowest order terms in momenta , and then in the general case , where all the momenta terms are taken into account . the fourth step consists in obtaining the renormalized dispersion in these approximations : bare and self - consistent approximations . the final stage is then to determine the spin stiffness in the approximations by taking the coefficient of the quadratic term in the renormalized dispersion . it is found that a naive extrapolation of the bare approximation to the regime of finite disorder strength predicts vanishing of the renormalized spin stiffness at a certain disorder strength @xmath5 , indicating a paramagnetic phase transition . a more realistic self - consistent approximation , however , predicts even faster decrease of the renormalized spin stiffness with growing disorder strength up to a certain critical value @xmath6 of the disorder . at this point , the renormalized spin stiffness drastically changes its behavior : it becomes nonanalytic , acquires an imaginary part , and the real part saturates at a certain positive value without reaching zero . such nonanalytic behavior can not be accessed by any finite number of perturbative corrections . in addition , our calculations show a strong indication that the pauli susceptibility also diverges at the same critical point @xmath6 , suggesting a phase transition , presumably to a spin glass phase . the outline of this paper is the following : in section ii we present the model and in section iii we derive the expression for the self - energy . then , we first solve the problem using the bare green s function in section iv . we present our numerical and analytical results for the self - consistent solution of the dyson equations in section v and draw our conclusions in section vi . the 2deg in the presence of both a perpendicular magnetic field ( @xmath7 ) at @xmath0 and disorder is described by the fermionic hamiltonian @xmath8 , with @xmath9 and @xmath10 here , @xmath11 and @xmath12 are , respectively , the fermionic creation and annihilation operators in coordinate space , @xmath13 denotes the effective mass of the electron , @xmath14 is the vector potential , @xmath15 stands for the effective land @xmath16-factor and @xmath17 is the bohr magneton . in addition , @xmath18 denotes the coulomb potential , with @xmath19 being the dielectric constant of the host semiconductor , and @xmath20 stands for the impurity potential , with @xmath21 being the random position of an impurity . the first step consists in obtaining a @xmath4 quantized version of the magnon hamiltonian of the system . in our model we consider only single magnon processes , which allow us to use a bosonic description . it was shown in ref . [ ] that the bosonized hamiltonian of the system in the absence of disorder is ( neglecting a constant term ) @xmath22 where @xmath23 and @xmath24 are , respectively , the bosonic creation and annihilation operators in @xmath25 space and the bosonic dispersion relation is given by @xmath26.\ ] ] here , @xmath27 stands for the coulomb energy scale ( @xmath28 being the magnetic length ) , @xmath29 denotes the modified bessel function of the first kind , and @xmath30 . it must be stressed that although the interaction between magnons is omitted from the discussion , the coulomb interaction between electrons up to rpa level is taken into account by the bosonic dispersion relation @xmath31.@xcite    we now focus on the impurity part of the hamiltonian . we begin with the fermionic expression of the @xmath4 quantized impurity hamiltonian , @xmath32 here , @xmath33 denotes the fourier transformed density function @xmath34 for the impurities and @xmath35 and @xmath36 are , respectively , the fermionic creation and annihilation operators in @xmath25 space . in order to obtain the bosonic form of the above , the fourier - transformed electronic density operator must be used . it is given by @xmath37 the electronic field operators are related to the single - electron operators by @xmath38 where @xmath39 is the area of the system . substituting the above back into eq.([arshavin1 ] ) gives @xmath40 then , substituting eq . ( [ eq7 ] ) back into eq.([arshavin2 ] ) yields @xmath41 the bosonized version of the electron density operator reads @xcite @xmath42 where @xmath43 is the landau level degeneracy and @xmath44 . the disorder hamiltonian then becomes @xmath45.\end{aligned}\ ] ] the constant term @xmath46 is now omitted since the quantity of interest is the green s function . the bosonized impurity hamiltonian is then finally written as @xmath47 where @xmath48    labelling @xmath49 the full bosonized hamiltonian of the quantum hall ferromagnet in the presence of impurities is then expressed as @xmath50    let us now say a few words on the dimensions of the disorder potential . there are two sources of disorder present in the system : impurities positioned at a certain distance away from the 2deg and impurities present in the 2deg . in the case of gaas heterostructures , @xcite most of the disorder potential is spawned by the coulomb interaction between the electrons and the impurities located away from the 2deg . these impurities correspond to ionized donor atoms situated in the @xmath51-type region , which itself is detached from the 2deg by an insulating layer of thickness @xmath52 . in the present calculations , the disorder potential will be taken as an effective two - dimensional potential . having obtained the bosonized hamiltonian in the presence of impurities , one is now able to determine the expression for the self - energy . in the same spirit as ref [ ] , one first looks for the green s function , @xmath53 | 0 \rangle.\ ] ] here , @xmath54 stands for the bosonic vacuum state , which is none other than the quantum hall ferromagnet : i.e. @xmath55 . thus , one has @xmath56b^{\dagger}_{\textbf{p } ' } ( 0 ) \right ] | 0 \rangle , \ ] ] where @xmath57 is the time ordering operator . now , using eq.([arda1 ] ) , one easily finds that @xmath58 = \omega_{\textbf{p } } b_{\textbf{p } } ( t ) + \sum_{\textbf{q } } u^{e } ( \textbf{q } , \textbf{p } - \textbf{q } ) g_{\textbf{q } } b_{\textbf{p } - \textbf{q } } ( t),\ ] ] such that one obtains for the second term in eq.([arda2 ] ) @xmath59b^{\dagger}_{\textbf{p } ' } ( 0 ) \right ] | 0 \rangle = \omega_{\textbf{p } } g ( \textbf{p } ' , \textbf{p};t ) \\ & + \sum_{\textbf{q } } u^{e}(\textbf{q } , \textbf{p } - \textbf{q } ) g_{\textbf{q } } g ( \textbf{p } ' , \textbf{p } - \textbf{q};t).\end{aligned}\ ] ] hence , the equation of motion of @xmath60 is written as @xmath61    the zero - order approximation to the solution of eq.([arda2bis ] ) yields @xmath62 where @xmath63 stands for the bare bosonic green s function . we now look for the expression for @xmath63 . firstly , one needs to find the heisenberg bosonic operator in the absence of the disorder potential . starting with @xmath64 = \omega_{\textbf{p } } b_{\textbf{p } } ( t ) $ ] , one then obtains @xmath65 . therefore , for the case @xmath66 , the free green s function is @xmath67 whereas for @xmath68 , it turns out to be @xmath69 this solution is indeed identical to the electronic one . now , the cynosure is on the generic solution of the differential equation ( [ arda2bis ] ) . by coupling the latter with the boundary equation ( [ arda3 ] ) yields the integral equation @xmath70 by fourier - transforming the time in eq . to frequency and shifting @xmath71 one finds @xmath72 here , the bare green s function reads @xmath73 where @xmath74 and @xmath75 is given by eq .. in the same way as for fermions , the solution of eq.([arda5 ] ) is obtained by iteration . one gets the so - called _ born series _ : @xmath76 where @xmath77 and for @xmath78 , @xmath79 expansion of eq.([saqibismyfriend1 ] ) then yields    @xmath80    due to disorder self - averaging in the limit of very large number of impurities @xmath81 , with constant density @xmath82 , the full bosonic one - particle green s function approaches its average value @xmath83 which is @xmath84 in the thermodynamic limit @xmath85 : @xmath86 moreover , one has @xmath87 substituting eqs . and eq.([firstorder ] ) into the expression for @xmath88 shows that the translational invariance is recovered after the averaging @xmath89 , where @xmath90 therefore , there is no @xmath91 order born scattering contribution to the bosonic self - energy . moreover , it is possible to show that all odd order contributions to the self - energy vanish ( see appendix a ) . this result is expressed diagrammatically in fig.[diaag9 ] . it was shown@xcite that the disorder averaged green s function can also be expressed as @xmath92    hence , the self - energy must now be computed . the low - density weak scattering approximation will be used throughout the calculations . low density means that the number of disorder atoms present in the system is taken to be much lower than the number of electrons , while the weak - scattering approximation signifies that the scattering potential induced by a given impurity atom is weak , such that only the first and second - order born scatterings are accounted for . the problem then reduces to solving the diagrammatic expression shown in fig.[diaag10 ] . the self - energy can be evaluated in two different manners : 1 ) the bare approximation that uses the bare propagator @xmath93 and 2 ) the self - consistent approximation that uses the full disorder self - averaged green s function @xmath94 . therefore , for generality we will use the propagator @xmath95 , which is going to be specified further for each particular case . this yields algebraically    @xmath96    where eq.([qzdr1 ] ) was substituted in the second line .    in this work , the impurity potential is assumed to be short - range , i.e. @xmath97 . an uniform potential in momentum space is attained from a delta function interaction in real space , @xmath98 ( such that @xmath99 ) . thus , this model assumes that the bosons ( and therefore the electrons ) collide directly with the impurity atoms that constitute the effective disorder potential ; in reality , most of the impurities are located away from the 2deg . thus , one has @xmath100 g ( \textbf{q } , \omega ) \left [ 2 i e^{- |\ell ( \textbf{p } - \textbf{q})|^2 / 4 } \sin \left ( \frac{(\textbf{p } - \textbf{q } ) \wedge \textbf{q}}{2 } \right ) \right ] \\ & = 4 n_{\rm imp } \left ( \frac{u}{a } \right)^2   \sum_{\textbf{q } } e^{- |\ell ( \textbf{q } - \textbf{p})|^2 / 2 } \sin ^2 \left ( \frac{\textbf{q } \wedge \textbf{p}}{2 } \right ) g ( \textbf{q } , \omega ) . \end{split}\label{fghw1}\ ] ] one then expands the argument of the exponential : @xmath101 here , @xmath102 denotes the angle between vectors @xmath25 and @xmath103 . furthermore , the summation is transmuted into an integration through the use of the formula , @xmath104 the angle @xmath105 is taken arbitrarily on the plane containing the vector @xmath25 , therefore , one is free to set @xmath106 . the sine squared term in eq.([fghw1 ] ) can be re - written as @xmath107 / 2 = [ 1 - \cos ( \ell^2 q p \sin \phi ) ] / 2 $ ] . we also assume rotation invariance of the green s function @xmath108 . hence , the self - energy is also rotation invariant and can be expressed as @xmath109,\ ] ] where @xmath110 stands for the impurity density . after a straightforward calculation ( see appendix b ) , we find @xmath111.\ ] ] rescaling the momenta by @xmath112 simplifies the self - energy to @xmath113 g ( q , \omega ) , \ ] ] where the various pre - factors , including the disorder potential strength and the impurity density , can be re - grouped into a single convenient parameter : @xmath114 which will be dubbed the _ disorder strength_. thus , @xmath115 is a dimensionless parameter that measures the disorder interaction strength relative to the coulomb interaction , @xmath116 . the above self - energy expression will be evaluated in two different ways : i ) first order corrections in @xmath115 and ii ) self- consistently .      in the bare approximation the self - energy ( [ fghw4 ] ) becomes @xmath117 g^0 ( q , \omega ) .\ ] ] after substituting eq.([berbatovatmanutd1 ] ) into the above , we obtain @xmath118 making use of the identity ( for @xmath74 ) , @xmath119 we find the real and imaginary parts of the self - energy @xmath120 \pi \delta ( \omega - \omega_{q } ) . \label{saaidismybro2bisbis}\end{aligned}\ ] ]    the above equations can be evaluated analytically in the long wavelength approximation , which is done in appendix c. here , one uses the complete bosonic dispersion relation given by eq .. as a result , one can only solve the imaginary self - energy numerically ; that task is not performed here . we concentrate , instead , on the real part . the renormalized energy of the bosons ( including the disorder contribution ) is obtained by looking at the poles of the full disorder self - averaged green s function in eq . , @xmath121 , such that the renormalized dispersion relation is determined from eq.([saaidismybro2 ] ) : @xmath122   + \frac{u}{4 } \epsilon_b^2 e^{-p^2 / 2 } \mathcal{p } \int^{\infty}_0 { \mathrm{d}}q \ q e^{-q^2 / 2 } \frac{i_0 ( q p ) -1}{\omega - \left [ g + \epsilon_b \left ( 1 - e^{- q^2 / 4 } i_0 \left ( \frac{q^2}{4 } \right ) \right ) \right ] } .\ ] ] the corresponding plot is illustrated on fig.[scottmurray3](a ) . one can notice that at not too large momenta ( i.e. near @xmath123 ) there exists already a substantial difference between the bare ( long wavelength ) and bare ( full @xmath2 ) approximations . approximation ( blue ) , in units of the coulomb energy @xmath124 , contrasted with the one in the bare long wavelength approximation ( red ) , both as functions of the momentum @xmath125 and at @xmath126 . ( b ) renormalized spin stiffness in the bare full @xmath2 approximation ( blue ) and in the bare long wavelength approximation ( red ) . notice that using the bare green s function @xmath127 we find a transition from a ferromagnetic to a paramagnetic phase , whereas using @xmath128 in the self - consistent approximation we find a transition into a spin glass phase ( see next section).,title="fig : " ]   approximation ( blue ) , in units of the coulomb energy @xmath124 , contrasted with the one in the bare long wavelength approximation ( red ) , both as functions of the momentum @xmath125 and at @xmath126 . ( b ) renormalized spin stiffness in the bare full @xmath2 approximation ( blue ) and in the bare long wavelength approximation ( red ) . notice that using the bare green s function @xmath127 we find a transition from a ferromagnetic to a paramagnetic phase , whereas using @xmath128 in the self - consistent approximation we find a transition into a spin glass phase ( see next section).,title="fig : " ]    now , the renormalized spin stiffness is sought for . for the sake of convenience , one begins by introducing the variables @xmath129 and re - writing eq.([kirilishere1 ] ) as @xmath130 } .\ ] ] one then expands the above in powers of p , @xmath131 } , \end{split}\ ] ] and one takes only the @xmath132 terms , @xmath133 }   + \ldots \\ & = \frac{p^2}{4 } \left[1 - \frac{u}{4 } ( 5.72 ) \right ] + \dots \end{split}\ ] ] the renormalized spin stiffness then reads @xmath134 eq.([mainbare ] ) is the main result of this section . the above expression was derived in the bare approximation , which takes into account only the lowest order corrections in @xmath115 . such assumption is only true in the realm of weak - disorder scattering.@xcite it can be seen that the renormalized spin stiffness decreases linearly in this approximation . a naive extrapolation of this dependence to the region of finite and strong disorder strength shows that there is a certain value @xmath135 , for which the renormalized spin stiffness vanishes ( in the long wavelength approximation @xmath136 ) , see fig.[scottmurray3](b ) . green @xcite explains that a vanishing renormalized spin stiffness at a threshold disorder strength means that the 2deg at @xmath0 undergoes a quantum phase transition from a ferromagnetic state to a paramagnetic one . thus , one can infer that the quantum hall ferromagnet undergoes a disorder - driven quantum phase transition to a paramagnetic state at critical disorder strength @xmath137 . it is also interesting to remark that green established this general finding in the domain of the weak disorder limit ( though in the context of a different model ) . the results obtained in this section can not be directly compared quantitatively with those of green @xcite , sinova _ et al . _ @xcite and rapsch _ et al . _ @xcite in addition to the fact that the model used in the studies of green is different , he does not complement his proposition on the vanishing of the renormalized spin stiffness with some quantitative results . sinova _ _ et al.__@xcite use a disparate variable in the ratio of the interaction strength to the landau - level broadening disorder energy scale . finally , rapsch _ _ et al.__@xcite perform their numerical calculations on a semiclassical spin model .    in the next section we evaluate the self - energy using the so - called self - consistent approximation and show that the renormalized spin stiffness drastically changes its behavior , which leads to completely different conclusions about the phase transition . the self - consistent approximation means that the self - energy is evaluated with the total disorder averaged green s function ( [ berabatamanutd4 ] ) instead of the bare one . therefore , one has ( see eq.([fghw4 ] ) ) @xmath138 \langle g_u ( q , \omega ) \rangle_{\rm imp},\ ] ] now , by referring to the computations carried out in the previous section and substituting eq.([berabatamanutd4 ] ) , one gets @xmath139 using that @xmath140 one has @xmath141 thus , one can write @xmath142 with @xmath143 using such expansion allows one to promptly get a numerical solution by iterations ( see fig.[diaag11 ] ) .    the convergence of the iterative solution is rather good up to some value of the disorder concentration @xmath6 . however , when @xmath144 , we find that @xmath145 . therefore , it would be desirable to derive an analytical solution in the neighborhood of @xmath6 . for convenience , we omit the arguments of @xmath146 in our notation in the next part . in general , @xmath147 or equivalently @xmath148    introducing for simplicity @xmath149 ^ 2}\ ] ] one finds @xmath150 introducing a matrix notation @xmath151 eq . reads @xmath152 its solution is found by computing the inverse matrix to eq . and has the form @xmath153 substituting this result into @xmath154 yields @xmath155 where @xmath156 with @xmath157 and @xmath158 where @xmath159 ^ 3}.\ ] ] suppose that @xmath160 when @xmath161 . in this case @xmath162 remains finite , as well as @xmath163 . this suggests that @xmath164 when @xmath165 , since the other terms are finite . moreover , if @xmath163 is a smooth function around @xmath6 , such that @xmath166 for some @xmath167 from the neighborhood of @xmath6 , then according to eq . there holds @xmath168=\sqrt{\kappa(u_c)(u_c - u)}+{\cal o}(u_c - u).\ ] ] it follows then from eq . that @xmath169 ^ 2/\kappa(u_0)$ ] as long as @xmath170 . however , the analysis of the infinite dimensional matrix @xmath171 and its determinant is quite complicated , which forces us to use an approximate solution , where we keep only the first @xmath172 terms in the expansion , thus reducing the dimension of the matrices to @xmath173 . in the absence of zeeman splitting ( @xmath174 ) , for @xmath175 , and @xmath176 one finds , setting @xmath177 , that @xmath178=0.0551776 $ ] and @xmath179 , which yields @xmath180 in excellent agreement with the numerical solution . the approximation also allows to check the validity of eq . , which yields @xmath181 at the point @xmath176 ( here the prime stands for the partial derivative with respect to @xmath115 ) . on the other hand , the numerical solution for the two points @xmath176 and @xmath182 yields @xmath183 , which agrees reasonably well with the previous result . the main difference stems from the fact that @xmath176 is rather close to the critical point @xmath6 , where the derivative diverges , so the value @xmath184 is still rather large and , of course , computational errors and approximation with finite number of terms make the result not very precise . furthermore , it follows from eq . that @xmath185 remains finite with @xmath165 . thus , @xmath186\sigma'_n(u_0)}{\sqrt{\kappa(u_0)}}\int^{u}_{u_0 } { \mathrm{d}}v(u_c - v)^{-1/2},\end{aligned}\ ] ] which leads to @xmath187\frac{\det[b(u)]}{\det[b(u_0)]}\ ] ] after performing the integration , where @xmath188 . from this analytic solution one may observe that @xmath146 and , consequently , @xmath189 acquires an imaginary part when @xmath190 . in particular , considering @xmath191 , for the case at hand @xmath192 and @xmath193 . defining @xmath194}{\sqrt{\kappa(u)}},\ ] ] the value of @xmath195 can be evaluated without any fitting parameters directly from eqs . and , which yields @xmath196 . it follows directly from the above that the renormalized spin stiffness now obeys @xmath197 where @xmath198 . both numerical and analytic results for @xmath199 are plotted in fig.[scottmurray5 ] , which shows that the analytic solution remains in excellent agreement with the numerical one even for those values of @xmath115 , which are far from the critical point @xmath6 . the whole behavior of the renormalized spin stiffness is very similar to the one obtained by chalker _ et al._,@xcite describing a spin glass phase transition . moreover , such dependence of the renormalized spin stiffness as a square root function of a control parameter was already observed previously by shender,@xcite as well as by avgin _ _ et al.__@xcite they considered the two- and tree- dimensional @xmath200 heisenberg spin glass model in a ferromagnetic ground state due to a strong external magnetic field . they found that for a certain value of the control parameter , @xmath199 acquires an imaginary part . the real part of @xmath199 is proportional to the spin wave stiffness , whereas the imaginary part is proportional to the damping of the spin wave excitations , thus signalling localization . it was argued that when the frequency of the spin - wave excitation @xmath201 multiplied by its lifetime @xmath202 is @xmath203/{\rm im}[\rho_s(u)]<1 $ ] , then the spin waves are completely localized . as we can see from the fig.[scottmurray5 ] , the condition of localization is already satisfied for the values of the disorder strength starting from @xmath204 . the calculations presented in the appendix d contain a strong indication that the pauli susceptibility diverge at the point @xmath205 , suggesting a phase transition from a ferromagnetic ground state to a spin glass state,@xcite since the spin waves become localized . our discussion was mainly concerned with the static case @xmath175 . however , our approach allows to find @xmath206 for any given @xmath201 . the dispersion spectrum in the self - consistent approximation then satisfies @xmath207 .     as functions of the disorder strength @xmath115 . ( square - analytic solution , star - numerical ) ] this paper accounts for the presence of both disorder and interactions in a 2deg at landau level filling factor @xmath0 , whose ground state constitutes the well - known quantum hall ferromagnet . the bosonization technique developed by doretto @xcite was employed in order to facilitate the treatment of both disorder and interactions in this strongly correlated system . the bosonization procedure consists in treating the spin wave ( magnon ) excitation as a boson such that the fermionic hamiltonian of the system can be approximately recast into a hamiltonian expressed in terms of bosonic operators . as a consequence , the interaction between electrons up to rpa level was incorporated within the bare propagator that represents the free boson . the intent was then to identify a disorder - driven quantum phase transition to a non - ferromagnetic state by analyzing the behavior of the renormalized spin stiffness as a function of the disorder strength , which itself corresponds to the ratio squared of the disorder energy scale to the coulomb energy one . to achieve this aim , firstly , we derived the bosonic expression for the hamiltonian of the system . in the second stage , the focus was on seeking out the disorder self - averaged green s function , which is the full bosonic green s function averaged over the impurity positions . then , by using the dyson s equation , we obtained a diagrammatic representation of the self - energy . the latter was subsequently computed within the framework of the low - density weak - scattering approximation . low density means that the number of disorder atoms present in the system is taken to be much lower than the number of electrons , while the weak - scattering approximation signifies that the scattering potential induced by a given impurity atom is weak , such that only the first and second - order born scatterings are accounted for . as a result , the self - energy corresponded to a single diagram . furthermore , the self - energy was evaluated in three different approximations : 1 ) the bare ( long wavelength ) approximation , which consists in using the bare bosonic propagator and keeping the lowest order terms in momenta , 2 ) the bare ( full @xmath2 ) approximation , which uses as well the bare bosonic propagator but with all the momenta terms kept in the calculation and , finally , 3 ) the self - consistent approximation , which uses the full disorder averaged green s function instead of the bare one in the self - energy diagram . then , the renormalized spin stiffness was determined by extracting the coefficient of the quadratic term in the dispersion relation together with the contribution from the self - energy . in the case of the bare ( long wavelength ) approximation , the spin stiffness was found to vanish linearly at the disorder strength @xmath136 . for the bare ( full @xmath2 ) scheme , the spin stiffness also vanished linearly , but at the disorder strength @xmath137 . these results suggest the occurrence of a disorder - driven quantum phase transition from the ferromagnetic phase to a paramagnetic one at the critical value @xmath137 . lastly , the self - consistent calculation revealed a completely different behavior : the real part of the renormalized spin stiffness also initially decreases with increasing the disorder strength @xmath115 , but then it saturates without reaching zero beyond a critical value @xmath6 , at which it ( and the self - energy ) acquires an imaginary component . according to the shender criterium,@xcite the spin waves become completely localized when the imaginary part of the renormalized spin stiffness becomes larger than the real part , which occurs in our system for @xmath208 ( see fig.[scottmurray5 ] ) . the physical mechanism behind a phase transition from the ferromagnetic ground state can be understood by considering electrons completely filling the lowest landau level ( @xmath0 ) in the presence of some inhomogeneous electrostatic background ( disorder ) . then , for sufficiently strong impurity potential , by adjusting the electron density to the electrostatic background , the system would gain more energy than is needed to rearrange the spin configuration . in this case the ferromagnetic state does not minimize the total energy of the system and a phase transition should take place . this quantum phase transition could be detected by calculating the behavior of the magnetic susceptibility as a function of the disorder strength . a sharp peak is anticipated at the transition point . in particular , if the energy cost for exciting a spin wave is less than the gain in the electrostatic energy , then the renormalized spin stiffness becomes negative and the system undergoes a phase transition to a paramagnetic state with zero local magnetization . on the other hand , as it was argued by rapsch _ et al . _ , @xcite in the case of a smoothly varying impurity potential , keeping nonzero local magnetization is still energetically favorable and the electrostatic energy is lowered by screening the impurity potential due to the formation of spin textures . at strong disorder such phase would correspond to a spin glass and the spin textures might be considered as the localized spin waves . thus , the character of the phase transition might depend on the nature of the disorder . the calculations performed within our model indicate that the pauli susceptibility diverges at the same critical point of the disorder strength @xmath6 , where an imaginary part of the renormalized spin stiffness appears , thus suggesting a phase transition to a spin glass phase . our approach can be extended for the case of bilayer systems in the presence of disorder . in fact , fertig and murthy @xcite have already considered such systems . thus , it would be interesting to apply our formalism to the case of a bilayer system with the total filling factor @xmath209 and compare the results . we acknowledge insightful discussions with t.  giamarchi and l.cugliandolo . we also would like to thank r.  l.  doretto for proposing us this interesting problem . this work was partially supported by the netherlands organization for scientific research ( nwo ) . let us now evaluate the @xmath210 order diagram contribution to the self - energy . its diagrammatic representation is shown in fig . 6.algebraically , we have @xmath211    [ diaag12 ]    here again , the impurity potential is short - range @xmath212 . by replacing all the functions defined previously , we find @xmath213 g ( \textbf{q } , \omega ) \\ & \qquad \times \left [ 2 i e^{- |\ell ( \textbf{q } ' - \textbf{q})|^2 / 4 } \sin \left ( \frac{(\textbf{q } ' - \textbf{q } ) \wedge \textbf{q}}{2 } \right ) \right ] g ( \textbf{q } ' , \omega ) \left [ 2 i e^{- |\ell ( \textbf{p } - \textbf{q}')|^2 / 4 } \sin \left ( \frac{(\textbf{p } - \textbf{q } ' ) \wedge \textbf{q}'}{2 } \right ) \right ] \\ & = n_{\rm imp } \sum_{\textbf{q } , \textbf{q } ' } \left ( \frac{u}{a } \right)^3   ( 2 i)^3 \underset{\textbf{term a}}{\underbrace { e^{-    \textbf{q}')|^2 / 4 } e^{- |\ell ( \textbf{q } ' - \textbf{q})|^2 / 4 } } } \underset{\textbf{term b}}{\underbrace{\sin \left ( \frac{\textbf{q } \wedge \textbf{p}}{2 } \right )   \sin \left ( \frac{\textbf{p } \wedge \textbf{q}'}{2 } \right ) } } \\ & \qquad \times \underset{\textbf{term c}}{\underbrace { \sin \left ( \frac{\textbf{q } ' \wedge \textbf{q}}{2 } \right ) } } g ( \textbf{q } , \omega ) g ( \textbf{q } ' , \omega ) . \end{split}\ ] ] it is clear that terms a and b are symmetric under the interchange @xmath214 while term c is antisymmetric . thus , one has @xmath215 this result holds true in both the full bare and self - consistent approximations .    as a matter of fact , due to the antisymmetric property of the wedge product within the sine term , it turns out that all odd order terms vanish . we prove here the expressions and . we begin with the expression of the self - energy given by eq.([fghw2 ] ) @xmath216 . \label{a1}\ ] ] one first deals with the polar integral , @xmath217 = \frac{1}{2 } \int^{2 \pi}_0 \frac{{\mathrm{d}}\phi } { 2 \pi } \ e^{\ell^2 q p \cos \phi } - \frac{1}{2 }   \int^{2 \pi}_0 \frac{{\mathrm{d}}\phi } { 2 \pi } \ e^{\ell^2 q p \cos \phi } \cos ( \ell^2 q p \sin \phi).\ ] ] the two terms are evaluated separately . for the first term , one must note that @xcite @xmath218 such that @xmath219 the second term @xmath220= \frac{1}{2}\ ] ] substituting eqs.([ghigbigb1 ] ) and ( [ ghigbigb2 ] ) back into eq.([ghigbigb3 ] ) then yields the simpler expression , @xmath221 = \frac{1}{2 } [ i_0 ( \ell^2 q p)-1 ] .\ ] ] now , substituting the above back into eq.([a1 ] ) we find @xmath222.\ ] ] to evaluate the self - energy within the long - wavelength approximation , we must return to eq.([fghw1 ] ) . firstly , one remarks that the sine squared term in eq.([fghw1 ] ) greatly simplifies , @xmath223 ^ 2 = \frac{1}{4 } ( \ell^2 |\textbf{q } \times \textbf{p}|)^2 = \frac{1}{4 } \ell^4 |\textbf{q}|^2 |\textbf{p}|^2 \sin^2 \phi.\ ] ] then , substituting eq.([longsine ] ) into eq.([fghw1 ] ) yields : @xmath224    the polar integral then turns out to match @xcite @xmath225.\ ] ] the series expansion for eq.([intbessel ] ) gives ( @xmath226 ) @xmath227 - \left [ \frac{x^2 } { 8 } + \ldots \right ] \\ & = 1 + \frac{x^2}{8 } + \ldots \end{split}\ ] ] eq.([stefanos1 ] ) already holds a @xmath228 term and therefore a @xmath229 term is not needed in the long wavelength approximation . thus , one assumes that @xmath230 moreover , the momenta are rescaled as @xmath231 . as a result , eq.([stefanos1 ] ) simplifies to @xmath232 where we replaced @xmath233 by its definition ( see eq . ) . one can then make use of the identity @xcite @xmath119 where @xmath234 symbolizes the cauchy principal value of the integral . consequently , one has @xmath235 let us first examine the real part of the self - energy , which actually denotes the physical self - energy . it has been shown@xcite that in the long wavelength approximation the bosonic dispersion relation for @xmath25 can be written as @xmath236 the physical self - energy then becomes @xmath237 let us then work temporarily with the new quantities @xmath238 such that the self - energy is re - written as @xmath239 now , one performs a change of variable in the @xmath240 momentum : @xmath241 . one must note that @xmath242 and that the integration limits are not altered . consequently , one gets @xmath243 a further change of the integration variable is performed @xmath244 , leading to    @xmath245 . \end{aligned}\ ] ]    term a corresponds to the definition of the exponential integral function;@xcite @xmath246 whereas term b can be straightforwardly integrated , @xmath247 where the shift of variable @xmath248 was used in the first step . thus , the physical self - energy becomes @xmath249 . \end{split}\ ] ]    the renormalized energy of the bosons is obtained by looking at the poles of the full disorder self - averaged green s function , @xmath250 consequently , in the long wavelength approximation , the renormalized dispersion relation takes the form @xmath251 . \end{split}\ ] ] it is straightforward to notice that the renormalized spin stiffness , which corresponds to the coefficient of the @xmath132 term , is given by @xmath252    we now turn to the imaginary part of the self - energy given by eq .. in the long wavelength approximation , the dirac delta function becomes @xmath253 by performing a change of variable in the @xmath240 momentum , @xmath241 and replacing eq . into eq . one gets @xmath254    finally , the scattering time , which amounts to the lifetime of the bosonic excitation , is given by @xmath255 it is clear that @xmath256 when @xmath257 , i.e. low energy quasiparticles are long - lived , with finite lifetime induced by disorder . the pauli susceptibility in case of linear response is given by the kubo formula : @xmath258 using the fourier transformation @xmath259 the susceptibility can be written as @xmath260 on the other hand , the operators @xmath261 can be written in the bosonized form@xcite @xmath262 where @xmath263 . thus , after substitution @xmath264 evaluation of the expectation value yields @xmath265 using the notation defined earlier in eq.([greendoniach ] ) . expanding the green s function @xmath266 into the born series and performing the disorder averaging one recovers the translational invariance @xmath267 . moreover , performing the fourier transformation in the time variable @xmath268 and introducing @xmath269 the susceptibility is      in the self consistent approximation the function @xmath271 obeys@xcite @xmath272.\ ] ] we are interested mostly in the static susceptibility @xmath273 . thus , in particular @xmath274.\ ] ] a spherically symmetric solution satisfies @xmath275.\ ] ] let us introduce a new function @xmath276 then eq.([pdef ] ) can be rewritten as @xmath277 or explicitly @xmath278 ^ 2 } h(q,\omega).\ ] ] notice that eq.([intgequ ] ) has the same form as eq.([selfderiv ] ) but with @xmath279 instead of @xmath280 , which is known to diverge @xmath281 when @xmath161 . in the next part we will demonstrate that @xmath282 also diverges , @xmath283 when @xmath165 . we are looking for a solution in the form @xmath284 substitution of eq.([solform ] ) into eq.([intgequ ] ) yields an expression , which looks similar to the equation previously obtained ( see eq.([matreq ] ) ) , @xmath285 where the function @xmath286 was defined earlier by eq.([fndef ] ) and @xmath287 ^ 2}.\ ] ] notice that @xmath288 and @xmath289 equivalently @xmath290 where @xmath291 was defined in eq.([matr ] ) . the solution is found by computing the inverse matrix to eq . and has the form @xmath292    therefore , @xmath293 and @xmath294 if the integral is convergent . otherwise , it has to be regularized , which we wo nt consider here . this leads to @xmath295 which can be further simplified by means of some algebraic transformations , @xmath296    despite the simplifications , the above expression is difficult to evaluate analytically , as well as numerically . however , since most of the terms there involve the inverse matrix , it is reasonable to suppose that if @xmath175 it diverges with @xmath297 as @xmath298^{-1}.\ ] ] on the other hand the susceptibility is given by @xmath299 thus , considering @xmath300 , we see that the integrand is divergent at @xmath175 with @xmath165 , which is definitely not enough to infer the divergence of the integral itself , but can be considered as an indication to such possibility .
Abstract:

Answer: this work treats the effects of disorder and interactions in a quantum hall ferromagnet , which is realized in a two - dimensional electron gas ( 2deg ) in a perpendicular magnetic field at landau level filling factor @xmath0 . we study the problem by projecting the original fermionic hamiltonian into magnon states , which behave as bosons in the vicinity of the ferromagnetic ground state . the approach permits the reformulation of a strongly interacting model into a non - interacting one . the latter is a non - perturbative scheme that consists in treating the two - particle neutral excitations of the electron system as a bosonic single - particle . indeed , the employment of bosonization facilitates the inclusion of disorder in the study of the system . it has been shown previously that disorder may drive a quantum phase transition in the hall ferromagnet . however , such studies have been either carried out in the framework of nonlinear sigma model , as an effective low - energy theory , or included the long - range coulomb interaction in a quantum description only up to the hartree - fock level . here , we establish the occurrence of a disorder - driven quantum phase transition from a ferromagnetic 2deg to a spin glass phase by taking into account interactions between electrons up to the random phase approximation level in a fully quantum description .