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 absence of diffusion of waves in disordered media was predicted by anderson 50 years ago @xcite . originally proposed in the context of electronic transport in disordered crystals , it has since been observed for different types of waves , including light and sound . recently , direct observations of the anderson localization by disorder @xcite and of a localization transition by quasiperiodic potentials @xcite of quasi one - dimensional ( 1d ) matter waves of ultracold atoms were reported . these experiments pave the way to the observation of new phenomena and shed new light on long standing problems , amongst which the question of possible anderson localization in presence of interactions . in the present paper we consider the case of an atomic vapor described as a weakly interacting bose gas in the presence of a weak disorder ( what is meant by `` weak '' here will be made quantitative in sec . [ model ] ) . in this configuration it has been shown theoretically in refs . @xcite and supported by numerical simulations in ref . @xcite that a small amount of disorder does not drastically alter the _ equilibrium _ properties of the system , but merely decreases the condensate and superfluid fractions . furthermore , even in the 1d limit considered in the present work , it has been experimentally demonstrated in refs . @xcite that one can observe global phase coherence in the presence of disorder and remain far from , say , the bose glass phase originally proposed by giamarchi and schulz and fisher _ et al . _ @xcite . here we are interested in _ transport _ properties . specifically , we study a quasi 1d , weakly interacting bose einstein condensate ( bec ) , propagating through a disordered potential . in this context , localization has been theoretically studied mainly for effective _ attractive _ interactions ( see , e.g. , @xcite and references therein ) , with less attention on the _ repulsive _ case we consider here ( see , however , refs . @xcite ) . in the absence of an external potential , ( repulsive ) interactions make the system superfluid and introduce a new characteristic speed in the system , the speed of sound @xmath2 . as mentioned above , when the speed @xmath3 of the bec relative to the external potential tends to zero , the addition of a weak random potential preserves superfluidity , although with a reduced superfluid fraction . what happens as @xmath3 increases ? this question was investigated in a previous publication @xcite , where the disordered potential , of length @xmath0 , was modeled by a series of randomly located delta peaks . for small velocities @xmath4 , perturbation theory shows that the superfluidity is preserved , e.g. , the flow is dissipationless and with a perfect transmission . in contrast , in the high speed limit @xmath5 , where the kinetic energy dominates over the interaction energy , the transport properties of the bec are deeply altered , and tend to those of the non - interacting gas , displaying an exponential damping of the transmission with length @xmath0 , a behavior characteristic of the strong anderson localization . thus , two limiting cases of stationary flow have been identified @xcite , with contrasting transport properties : superfluidity in the deep subsonic regime , and anderson localization in the deep supersonic one . in between , in the region @xmath6 where both interaction and kinetic energies are important , it was shown that stationary scattering solutions do not exist : one reaches a regime of time - dependent flows with more or less ( depending on the speed ) complex density excitations . the range of speeds around @xmath2 where this phenomenon is observed increases as the length @xmath0 increases . the different types of existing flows are summarized in fig . [ fig1 ] .     through a disordered potential @xmath7 consisting in a series of uncorrelated delta peaks extending over a domain of size @xmath0 [ cf . ( [ e2a ] ) and the discussion in sec . [ sec - lstar ] ] . dark region : time dependent flow ; light gray ( light blue online ) regions : stationary flow . in the supersonic case , the yellow solid line corresponds to the threshold @xmath1 between these two domains as determined from eq . ( [ tstar - impl ] ) . the blue dashed line is the localization length @xmath8 ( [ exp1 ] ) . the supersonic region below @xmath8 denoted as `` ballistic '' corresponds to the region where the perturbation theory of sec . [ pts ] applies . note the enlarged scale for @xmath9 $ ] . [ fig1 ] ]    in the present study we concentrate on the supersonic stationary region of the phase diagram [ gray ( light blue online ) @xmath10 region in fig . [ fig1 ] ] . in this domain we provide analytical and numerical evidence of anderson localization in the presence of interaction for different types of disorder . we compute analytically the interaction - dependent localization length as well as the corresponding distribution of transmission coefficients . we also explain the disappearance of the supersonic stationary flow observed at a given velocity for increasing length of the disordered sample . this onset of time dependence is an important qualitative effect revealed by our study . we show that it is directly connected to interaction effects and provide an analytical estimate of the length @xmath1 of the disordered region above which most of the realizations of the random potential lead to time dependent flows ( see fig . [ fig1 ] ) . the paper is organized as follows . in sec . [ model ] we present the model and identify its range of validity . in sec . [ seciii ] we take some time to properly define the transmission coefficient of a bose - condensed beam over an obstacle . in [ seciv ] we introduce the different types of disordered potentials studied in the present work . in sec . [ ssr ] we present analytical and numerical results showing that anderson localization is indeed possible in the supersonic regime . we consider the three possible supersonic regimes : perturbative ( sec . [ pts ] ) , anderson localized ( sec . [ non - pert ] ) and onset of time dependence ( sec . [ sec - lstar ] ) . in sec . [ experiment ] we discuss experimental strategies and possible signatures for the observation of anderson localization in an interacting bose - einstein condensate . finally we present our conclusions in sec . [ conclusion ] . some technical points are given in the appendixes . in appendix [ app1 ] we derive the probability distribution of transmission in a special case ( perturbative regime and correlated gaussian potential ) . in appendix [ app2 ] we present the derivation of the fokker - plank equation ( [ dmpk ] ) for the distribution of the transmission coefficients . we study here the transport properties of a quasi - one - dimensional ( 1d ) bose - einstein condensate formed of particles of mass @xmath11 , experiencing a repulsive effective interaction ( characterized by the 3d s - wave scattering length @xmath12 ) , in the presence of an obstacle represented by the external potential @xmath13 . the potential is not necessarily disordered at this point , the only restriction we impose throughout the present work is that it should have a finite extent , i.e. , @xmath14 when @xmath15 . the configuration we consider corresponds to the `` 1d mean field regime '' @xcite ( see also the discussion in ref . @xcite ) , where the system is described by a 1d order parameter @xmath16 depending on a single spatial variable : the coordinate @xmath17 along the direction of propagation . @xmath16 obeys the nonlinear schrdinger equation @xmath18\psi\ ; .\ ] ] in all the present work we choose to work in the `` laboratory frame '' where the condensate is initially at rest . ( [ e1 ] ) describes its 1d dynamics in the presence of an obstacle moving at constant velocity @xmath3 in this frame , which corresponds to the experimental situation where an obstacle is swept through a condensate initially at rest , see e.g. refs . @xcite . on the theoretical side , one should imagine that , from an initial static configuration where the condensate is at rest with @xmath19 , the potential intensity and speed have been slowly ramped up to a point where the condensate dynamics is described by eq . ( [ e1 ] ) . we choose @xmath20 , this corresponds to a potential moving from left to right in the laboratory frame . the reduction of the motion of the condensate to a single spatial dimension is typically achieved through a transverse harmonic confining potential of pulsation @xmath21 . we choose a normalization such that @xmath22 is the linear density of the condensate . in this case , the interaction amongst particles results in eq . ( [ e1 ] ) in the nonlinear term @xmath23 , with @xmath24 @xcite .    in the stationary regime , where the flow is time  independent in the frame moving with the potential , @xmath25 depends on @xmath17 and @xmath26 only through the variable @xmath27 . the appropriate boundary condition is @xmath28 ( where @xmath29 is a constant ) ( see @xcite and the discussion in section [ stationary_regime ] below ) . the condensate is then characterized by a chemical potential @xmath30 , a speed of sound @xmath31 and a healing length @xmath32 . it is customary to characterize the transverse confinement via the `` harmonic oscillator length '' @xmath33 . with @xmath34 denoting a typical order of magnitude of @xmath35 , the 1d mean field regime in which eq . ( [ e1 ] ) is valid corresponds to a density range such that @xmath36 in this domain the wave function of the condensate can be factorized in a transverse and a longitudinal part @xcite . the transverse wave function is gaussian ( this is ensured by the condition @xmath37 ) , the longitudinal one is of the form @xmath38 and @xmath16 satisfies eq . ( [ e1 ] ) @xcite . the left - hand side ( l.h.s . ) inequality in eq . ( [ s1.0 ] ) prevents the system to enter in the tonks - girardeau regime . more precisely , a general analysis of 1d bose gas shows that at zero temperature no bec is possible @xcite . this results in a algebraic decrease of the one body density matrix monitored by phase fluctuations occurring over a phase - coherence length @xmath39 @xcite . hence the results obtained using eq . ( [ e1 ] ) are valid if they describe structures with a characteristic length scale smaller than @xmath40 . the l.h.s . inequality in eq . ( [ s1.0 ] ) ensures that @xmath40 is exponentially large compared to the healing length . if one considers , for instance , @xmath41rb or @xmath42na atoms in a guide with a transverse confinement characterized by @xmath43 hz , the ratio @xmath44 is roughly of order @xmath45 and restriction ( [ s1.0 ] ) still allows the density to vary over four orders of magnitude . even if the mean field approach is legitimate in 1d , the effects of disorder have to be taken into account with some care . it may well be that the introduction of a disordered potential @xmath13 in eq . ( [ e1 ] ) modifies the properties of the ground state . this is indeed the case as shown in refs . @xcite  : a disordered potential decreases the condensate and the superfluid fraction , but the effects are weak provided the intensity of the disorder remains weak ( see ref . @xcite for an extension to finite temperature ) . more precisely , in the case of a disorder formed by randomly spaced delta impurities with density @xmath46 ( see sec .  [ delta ] ) one can show @xcite that , in the dilute impurity limit , at @xmath47 the non - superfluid fraction ( normal part ) is proportional to @xmath48 [ the notations are those of eq . ( [ e2a ] ) ] and thus remains small provided the dimensionless coefficient @xmath49 is small ( weak disorder limit ) . at finite @xmath3 , the normal fraction is multiplied by a factor @xmath50^{-3/2}$ ] ( see ref . @xcite ) , which diverges when @xmath51 . one thus expects the mean field approach to fail near the region @xmath52 of fig . this is supported by the numerical results presented in @xcite . hence , in the center of the time - dependent region of fig . [ fig1 ] we can not trust the results obtained from eq . ( [ e1 ] ) . however , far from this region , the 1d mean field approach is expected to be valid even in presence of ( weak ) disorder , as experimentally demonstrated in refs .  @xcite . in the present work we characterize the localization properties of the condensate in the random potential by studying the transmission coefficient . ( [ e1 ] ) being non linear , the definition of transmission and reflection coefficients needs to be treated with special care . this is the purpose of the present section where we first define the stationary regime ( section [ stationary_regime ] ) and then the transmission coefficient within this regime ( section [ transmis ] ) . it is customary to perform a madelung transformation and to write @xmath53 where @xmath35 is the density and @xmath54 the local velocity . from ( [ e1 ] ) one can check that they verify the continuity equation @xmath55 the stationary regime is defined as the regime where the system is at rest in the frame moving with the obstacle . in this case , in the laboratory frame @xmath25 , @xmath56 , @xmath57 and @xmath58 are time dependent , but they depend on @xmath17 and @xmath26 only through the variable @xmath27 . it is then possible to get a first integral of ( [ e2 ] ) under the form @xmath59    in the case of subsonic ( @xmath60 ) and stationary motion , the flow is superfluid and the order parameter is only affected in the vicinity of the obstacle , with @xmath61 and @xmath62 @xcite .    for @xmath63 , a regime of stationary flow also exists but in this case the obstacle induces density oscillations with a pattern stationary in its rest frame @xcite . this means that in the laboratory frame the phase velocity of these waves is identical to the velocity @xmath3 of the obstacle . on the other hand , the energy transferred from the obstacle to the fluid propagates with the group velocity , which in the case of bogoliubov excitations is greater than the phase velocity , i.e.  as just argued  than @xmath3 . as a consequence , radiation conditions require that the wake is always located ahead of the obstacle , i.e. , upstream , with no long - range perturbation of the fluid on the downstream side @xcite . this means that in this case the flow far in the downstream region remains unperturbed , with @xmath64 and @xmath65 . the two possible stationary configurations ( subsonic and supersonic ) are represented in fig . [ fig2 ] . ( 9,5.5 ) ( 1,0 ) ( @xmath66 ) . in both plots the potential is represented by a thick solid line ( hatched down to zero ) and the density profile is represented by a thin solid line.[fig2],title="fig:",width=302 ] ( 0.5,3.)downstream ( 6.7,3.)upstream ( 1,2.7)(1,0)7.5 ( 6.6,2.3)@xmath27 ( 1.4,1.7)@xmath67 ( 1.4,5.3)@xmath67 ( 3.2,1.)@xmath68 ( 3.8,5.)@xmath68 ( 0,1.5)@xmath29 ( 0.4,1.55)(1,0)0.5 ( 0,5.03)@xmath29 ( 0.4,5.08)(1,0)0.5    hence , in any stationary configuration ( subsonic or supersonic ) , the above reasoning fixes the integration constant in the right hand side ( r.h.s . ) of eq . ( [ e3 ] ) to its value at @xmath69 , i.e. , @xmath70 . in the stationary regime one gets from eqs . ( [ e1 ] ) and ( [ e3 ] ) @xmath71 where @xmath72 and @xmath73 \ ; .\ ] ]      in this section we restrict the analysis to the stationary regime of section [ stationary_regime ] , and define the transmission of the condensate through the obstacle represented by a potential @xmath13 ( not necessarily disordered ) verifying @xmath74 . as the wave equation ( [ e1 ] ) is nonlinear , one can not , in general , properly define reflection and transmission coefficients , since it is generally not possible to disentangle incoming and reflected waves in the nonlinear flow upstream the obstacle . however , following a procedure devised in ref . @xcite ( see also @xcite ) , we will show that one can define a transmission and a reflection coefficient in the limit of small nonlinearity as well as in the limit of weak reflection and arbitrary nonlinearity . outside the scattering region , @xmath75 and one can get a first integral of eq . ( [ e4 ] ) under the form @xmath76 which defines the `` free '' asymptotic density profiles . @xmath77 in eq . ( [ t1 ] ) are integration constants . the boundary condition discussed in the previous section imposes @xmath78 and @xmath79 when @xmath80 . this fixes the value @xmath81 . the value of @xmath82 at @xmath83 has to be determined by the integration of eq . ( [ e1 ] ) ( cf . ref . ( [ t1 ] ) expresses the energy conservation for a fictitious classical particle with `` mass '' @xmath84 , `` position '' @xmath85 and `` time '' @xmath86 , evolving in a potential @xmath87 ( whose typical shape is displayed in fig . [ fig3 ] ) . this type of analysis is common in the study of nonlinear equations such as eq . ( [ e1 ] ) , see e.g. , the review @xcite ( the first time we found it used is in ref . it is employed here as a convenient tool for getting intuition about the behavior of the solution of the gross - pitaevskii equation ( see below ) .     as a function of @xmath88 ( drawn for @xmath89 ) . the fictitious particle has a classical energy @xmath82 when @xmath90 . the ( red ) dashed line corresponds to an approximation of @xmath91 by @xmath92 , obtained by keeping only the first term in expansion ( [ t3a]).[fig3 ] ]    from now on we restrict to the supersonic stationary regime where an imperfect transmission occurs ( in the subsonic stationary regime one has perfect transmission ) . in this case the fictitious particle is initially ( i.e. , when @xmath66 ) at rest at the bottom of the potential @xmath87 with @xmath93 . the behavior of @xmath85 for @xmath94 depends on the value of @xmath95 . a stationary solution exists only if @xmath96 remains bounded , i.e. , if @xmath97 ( where @xmath98 corresponds to the local maximum of @xmath87 , see fig . [ fig3 ] ) . in this case , the asymptotic behavior of @xmath99 corresponds to oscillations between the values @xmath100 and @xmath101 defined in fig . 3 .    for future references we note that @xmath91 is zero when @xmath102 and when @xmath103 , and that the derivative @xmath104 is zero when @xmath102 and when @xmath105 , with @xmath106 at large velocity , when @xmath107 , one has @xmath108 and @xmath109    writing @xmath110 we now argue that , following ref . @xcite , one can write a perturbative version of ( [ t1 ] ) in a limit where @xmath111 we emphasize that restriction ( [ t2 ] ) corresponds to @xmath112 ( i.e. , small oscillations ) only when @xmath3 is of order of @xmath2 or smaller . the approach developed below is however able to tackle large relative density oscillations ( @xmath113 ) at large velocities ( @xmath107 ) @xcite . in this sense it will allow us to penetrate in the non perturbative regime where the upstream density oscillations are large and the transmission is low . using the variable @xmath114 , we write eq . ( [ t1 ] ) in the upstream region ( @xmath90 ) as @xmath115 where @xmath116 . a simple limited expansion around @xmath117 yields @xmath118 where @xmath119 and @xmath120 the second term in the r.h.s . of expansion ( [ t3a ] ) is small compared to the first one precisely in the limit ( [ t2 ] ) . in the following we restrict to this regime and neglect the second term of the r.h.s . of ( [ t3a ] ) . this corresponds to approximating the exact @xmath91 by the ( red online ) dashed line in fig . [ fig3 ] and to write eq . ( [ t3 ] ) under the form @xmath121 where the dimensionless parameter @xmath122 is defined by @xmath123 the solution of eq . ( [ t4 ] ) is @xmath124 where @xmath125 and @xmath126 is a phase factor . we recall that eq . ( [ t6 ] ) describes the density oscillations in the upstream region @xmath90 . these oscillations can be described as the sum of incident and reflected waves ( @xmath127 and @xmath128 ) of the form @xmath129 this analysis allows one to determine the reflection and the transmission coefficients as @xmath130 of course the sum of the incident @xmath127 and the reflected @xmath131 waves ( [ t6c ] ) is an approximate solution of the nonlinear schdinger equation ( [ e1 ] ) which is only valid in regime ( [ t2 ] ) , i.e. , in the regime of arbitrary interaction and small transmission ( @xmath132 ) , or in the regime of arbitrary transmission and small interaction ( @xmath107 ) . up to this point we presented a theory valid for any potential of finite extent . from now on we concentrate on the particular case of random potentials . we denote @xmath133 an arbitrary random potential , and use a subscript when dealing with one of the particular cases defined below .      the first potential of interest , analyzed in ref . @xcite , is a series of @xmath135 randomly located identical delta peaks of the form  : @xmath136 the intensity of the peaks is measured by the ( non random ) positive quantity @xmath137 . the scatterers have random uncorrelated positions @xmath138 , with mean density @xmath139 and average separation @xmath140 . hence the potential extends over a mean length @xmath141 .    denoting henceforth the disorder average by @xmath142 , for @xmath17 and @xmath143 inside the disordered region one gets the mean value @xmath144 and the irreducible two - point correlation function @xmath145 where @xmath146      another commonly used model of disorder is provided by gaussian random processes with zero average . we consider here potentials which are non zero only over a region of finite extent ( with typical size @xmath0 ) and generate them in the following way ( see , e.g. , chap . 5 of ref . @xcite and references therein ) : let s consider a gaussian white noise @xmath147 extending over all the real axis , with zero mean and second moment @xmath148 . then for a given function @xmath149 one defines @xmath150 where @xmath151 is a parameter characterizing the disorder and whose meaning is explained below . if @xmath152 were a delta function , then @xmath153 would be a gaussian white noise over @xmath154 $ ] ( and zero everywhere else ) . the actual function @xmath149 has a finite extension , and this induces finite correlations in the disordered potential .    from ( [ e2d ] ) it is clear that @xmath155 . if the domain of integration in the r.h.s . of ( [ e2d ] ) were extended to all @xmath156 , @xmath153 would have a gaussian distribution @xmath157}{\sqrt{2\pi\sigma^2 } } \ ; , \ ] ] where @xmath158 defining the correlation function @xmath159 as @xmath160 one would get in this case @xmath161 with a fourier transform @xmath162 where @xmath163 is the fourier transform of @xmath152 . imposing here the normalization condition @xmath164 leads to a two - point correlation function ( [ e2f0 ] ) which is  as in eq . ( [ e2b ] )  of the form of @xmath165 multiplied by a function whose integral over @xmath17 equals unity [ @xmath166 in ( [ e2f0 ] ) instead of @xmath167 in ( [ e2b ] ) ] . thus , with definition ( [ e2d ] ) and normalization ( [ e2cbis ] ) , @xmath151 plays for disorder ( [ e2d ] ) the same role as @xmath168 [ eq . ( [ e2c ] ) ] for disorder ( [ e2a ] ) : it characterizes the amplitude of the fluctuations of the potential . the typical extent of @xmath149 will in turn characterize the range of the correlations .    since @xmath153 as given by eq . ( [ e2d ] ) is typically non zero only over a region of finite extent , eqs . ( [ e2f0 ] ) and ( [ e2f ] ) are only correct if @xmath17 and @xmath143 are inside this region . more precisely , they should be in this region , at a distance from 0 or @xmath0 larger than the typical extent @xmath169 of the function @xmath152 . in the following we always consider the case where @xmath0 is very large compared to @xmath169 ( otherwise one could simply not speak of a disordered region ) and it is clear that the characteristics of the disorder are properly defined only inside the disordered region . we consider two special cases of correlation corresponding to different forms of @xmath152 : a lorentzian @xmath170 and a gaussian @xmath171 we denote the corresponding potentials by @xmath172 and @xmath173 . for the correlation functions one gets respectively @xmath174 and @xmath175 in both cases @xmath169 is the typical correlation radius .    the choice of a lorentzian correlated disordered potential originates from experimental and theoretical results in the case of micro - fabricated guides . in this type of setting , the atoms are magnetically guided over a chip @xcite . unavoidable imperfections and irregularities in the design of the circuit induce fluctuations in the current which , in turn , result in a random contribution to the magnetic field used for guiding the atoms . thus the potential seen by the atoms has a random component which is typically lorentzian correlated , with a correlation length @xmath169 which decreases when the distance between the guide and the chip increases @xcite . the gaussian correlated potential @xmath173 is more academic but , by comparison with the results obtained with @xmath172 it allows one to check what is really specific to the lorentzian case , and what is a mere effect of finite correlation length . another experimentally relevant type of disorder is the so called speckle potential which is generated by an optical speckle field produced by a laser beam passing through a diffusing plate @xcite . the corresponding potential will be denoted by @xmath176 and may be mathematically generated as follows @xcite : @xmath177    { \rm d}y \right|^2 , \ ] ] where @xmath178 and @xmath179 are two independent gaussian white noise processes of zero mean with @xmath180 ( @xmath181 and @xmath182 or 2 ) .    here also we characterize the disorder by studying its statistical properties in the limit where the domain of integration in the r.h.s . of ( [ e2l ] ) is extended to all @xmath156 . in this case one gets @xmath183 where @xmath184 is given by formula ( [ s2 ] ) ( replacing @xmath152 by @xmath185 ) . this yields @xmath186 and the correlation function defined in eq . ( [ e2f0 ] ) reads here @xmath187 \\ = & 4 \left [ \int_{\mathbb{r}}\!\!{\rm d}y\ , w_{\sss s}(x - x'+y)w_{\sss s}(y ) \right]^2   \ ; . \end{split}\ ] ] contrarily to the choice ( [ e2cbis ] ) , @xmath185 should not be normalized to unity here because  from eq . ( [ e2l ] )  this is homogeneity - wise impossible . instead , the choice @xmath188 corresponds to the typical experimental situations @xcite and leads to a correlation function @xmath189 whose integral over @xmath17 equals unity and whose fourier transform is @xmath190 hence , definition ( [ e2l ] ) and choice ( [ ws ] ) correspond here also to characterizing the amplitude of the disorder s fluctuations by the parameter @xmath151 and the range of the correlations by @xmath169 . as explained in ref . @xcite , and recalled in the introduction , anderson localization in a weakly repulsive bose - einstein condensate is only possible in the supersonic regime ( cf . [ fig1 ] ) which we consider now . in the present section we first analyze the transmission across a short disordered sample , in which case perturbation theory is applicable ( sec . [ pts ] ) . we then turn to generic non - perturbative configurations ( sec . [ non - pert ] ) where anderson localization is expected . in this regime we obtain evidences of the occurrence of anderson localization in the presence of interaction . finally we discuss the upper limit of the localized regime and the onset of time dependent flows for long disordered samples ( sec . [ sec - lstar ] ) . in the supersonic stationary regime , simple perturbation theory yields @xmath191 where @xcite @xmath192 \ ; , \ ] ] and @xmath193 is given by eq . ( [ pt4 ] ) . perturbation theory always predicts a stationary density profile . this is certainly wrong when @xmath3 is close to @xmath2 ( cf . fig . [ fig1 ] ) , but in this case @xmath193 gets very small and one precisely goes out of the domain of validity of perturbation theory ( @xmath194 as given by ( [ spt1 ] ) is no longer small compared to @xmath29 ) .    far ahead of the obstacle ( in a region where @xmath195 is larger than @xmath169 and @xmath196 ) ( [ spt1 ] ) gives @xmath197 the perturbative regime in which eqs . ( [ spt1],[spt2 ] ) are valid is also the one where the constant @xmath122 in ( [ t5],[t6 ] ) is small compared to unity . this can be inferred from the comparison of ( [ spt2 ] ) and ( [ t6 ] ) which indeed shows that @xmath132 in the regime where ( [ spt2 ] ) holds and that , in this case , @xmath198 ) , yielding @xmath199 from ( [ spt3 ] ) it is clear that the average reflection coefficient is @xmath200 furthermore , one can show that the corresponding probability distribution of the reflection coefficient is poissonian with @xmath201 note that for properly normalizing this probability distribution for @xmath202 $ ] , one should include in the prefactor of the r.h.s . of ( [ spt4 ] ) a correcting term of order @xmath203 which can be safely neglected in the limit ( [ spt5 ] ) . we give in appendix [ app1 ] a demonstration of result ( [ spt4 ] ) for the special case of a correlated gaussian potential @xmath153 of type ( [ e2d ] ) . below , we show that the same result holds for a potential @xmath7 of type ( [ e2a ] ) [ see eq . ( [ dmpk0 ] ) in section [ non - pert ] ] , and we checked numerically that it is also the case for the speckle potential @xmath176 ( [ e2l ] ) ( cf fig . [ fig4 ] ) . in all these cases , the average reflection coefficient reads ( up to the above discussed exponentially small correction ) @xmath204 where @xmath205 we recall that the function @xmath206 depends on the type of disorder considered . for a potential of form ( [ e2a ] ) one has @xmath207 , for the other potentials considered in this work it is given by ( [ e2j ] ) , ( [ e2k ] ) and ( [ csq ] ) .    concomitantly to distribution ( [ spt4 ] ) of reflection coefficients one gets for the transmission @xmath208 from ( [ e7b ] ) [ or ( [ spt10 ] ) ] one gets @xmath209 the perturbative approach holds when @xmath210 i.e. , when @xmath211 . this corresponds to the region which is denoted as `` ballistic '' in fig . [ fig1 ] @xcite . its accuracy is shown for @xmath212 in fig . [ fig4 ] for a speckle potential @xmath176 of type ( [ e2l ] ) ( we also checked this prediction for the potentials @xmath213 and @xmath173 , with also excellent results ) .     for the transmission coefficient @xmath214 in a potential @xmath176 with @xmath215 , @xmath216 and @xmath217 moving at velocity @xmath218 in a condensate of initial constant density @xmath219 . the corresponding localization length is @xmath220 . the histogram corresponds to a statistical analysis of the results of the numerical solution of eq . ( [ e1 ] ) for 10000 different random potentials . the solid line is the perturbative result ( [ e7b]).[fig4 ] ]    at this stage , @xmath8 is simply a notation for expression ( [ spt11 ] ) , but it will be shown to be the actual localization length of the matter wave in a disordered potential ( in section [ non - pert ] ) . the results derived here also hold for a noninteracting gas , obtained by taking the limit @xmath221 in eq . ( [ e1 ] ) , in which case @xmath222 and @xmath223 . ( [ spt11 ] ) , with @xmath193 replaced by @xmath224 , then coincides with the antsygina - pastur - slyusarev formula for the localization length @xcite and the distribution of transmissions ( [ e7b ] ) holds , with @xmath225 .    in the present work , those formulas are modified to include interactions . the generalization simply consists in replacing the wave vector @xmath226 by @xmath227 . this replacement has , as an important physical consequence , the effect of diminishing , at a given speed @xmath3 , the localization length ( there is an effective reduction of the available kinetic energy by the repulsive interactions ) . for instance , in the case of a potential @xmath7 , since @xmath228 , there is a relative difference @xmath229 between @xmath230 and @xmath231 , that is 11 @xmath232 for @xmath233 . this is illustrated in fig . [ fig5 ] , which displays the average @xmath234 as a function of @xmath0 for a disorder @xmath7 of type ( [ e2a ] ) , with and without interactions .     for a potential @xmath7 ( characterized by @xmath235 and @xmath236 ) . in the interacting case @xmath237 and @xmath238 . the non - interacting case is drawn for the same velocity and corresponds to a value @xmath239 . in both cases the dashed line is the analytical result ( [ e7 ] ) and the solid line corresponds to a statistical analysis of the results of the numerical solution of eq . ( [ e1 ] ) for 15000 different random potentials . the departure of the numerical results from the dashed lines occurs when the systems leaves the perturbative regime.[fig5 ] ]      when the size @xmath0 of the sample is large compared to the value @xmath8 determined in section [ pts ] , the perturbative approach fails . we now propose a non perturbative method allowing to treat both the regimes @xmath240 and @xmath241 and showing that @xmath8 , as defined in eq . ( [ spt11 ] ) , is indeed the localization length in the presence of interactions .    within the framework of the non perturbative approach , we are able to provide approximate analytical results in the case of the model disorder potential ( [ e2a ] ) . this potential being zero between two impurities , one can write a series of first integrals of eq . ( [ e4 ] ) in each segment @xmath242x_n , x_{n+1}[$ ] as follows  : @xmath243 = e^{(n)}_{\rm cl } \ ; .\ ] ] in the region @xmath244 the integration constant @xmath245 of eq . ( [ t1 ] ) is denoted as @xmath246 in ( [ np1 ] ) ( taking @xmath247 ) and is zero , whereas in the region @xmath248 one has @xmath249 ( and @xmath250 ) .     as a function of @xmath85 ( drawn for @xmath89 ) . for @xmath244 , the fictitious particle is initially at rest at the bottom of potential @xmath87 with @xmath251 . the value of the classical energy changes from @xmath252 to @xmath253 at each impurity @xmath254 . the lower panel displays the corresponding oscillations of @xmath99 , with two impurities at @xmath255 and @xmath256 ( their position is indicated by vertical dashed lines).[fig6],title="fig : " ]   as a function of @xmath85 ( drawn for @xmath89 ) . for @xmath244 , the fictitious particle is initially at rest at the bottom of potential @xmath87 with @xmath251 . the value of the classical energy changes from @xmath252 to @xmath253 at each impurity @xmath254 . the lower panel displays the corresponding oscillations of @xmath99 , with two impurities at @xmath255 and @xmath256 ( their position is indicated by vertical dashed lines).[fig6],title="fig : " ]    from ( [ e4 ] ) it is a simple matter to show that the matching condition of the density at impurity position @xmath254 is @xmath257 where @xmath258 [ @xmath259 denotes the limit of the derivative @xmath260 at the left [ at the right ] of @xmath254 . relation ( [ np1a ] ) between the derivatives of the amplitude results [ from eq . ( [ np1 ] ) ] in a relation between the classical energies : @xmath261 \ ; .\ ] ]    hence , eq .  ( [ np1 ] ) allows to draw a classical analogous of the solution of the nonlinear schdinger equation in the presence of potential ( [ e2a ] ) formed by a series of delta peaks : the fictitious classical particle defined in sec . [ transmis ] evolves in the potential @xmath87 and experiences kicks at `` times '' @xmath254 . each kick changes the `` energy '' according to ( [ np2 ] ) , as illustrated in fig . the key point in the remaining of this section will be to derive the probability distribution of @xmath262 which then directly allows one to get the distribution of @xmath122 s and of the transmission coefficients [ through eqs . ( [ t5 ] ) and ( [ t7 ] ) ] . let us introduce the quantities @xmath263 the parameters @xmath122 and @xmath264 defined in ( [ t5 ] ) and ( [ t6b ] ) are related to the ones of eq . ( [ np2a ] ) by @xmath265 and @xmath266 ( i.e. , @xmath122 is the last of @xmath267 s ; the same holds for @xmath264 ) . denoting by @xmath268 the value of the phase @xmath126 [ appearing in([t6 ] ) ] for @xmath269x_{n-1},x_{n}[$ ] , one gets in this domain [ the derivation is exactly the same as for eq . ( [ t6 ] ) ] @xmath270 and one can rewrite eq . ( [ np2 ] ) as @xmath271 \nonumber \\   & + & \frac{4\hbar^2}{m b^2 } \lambda_{n-1 } \sqrt{\kappa^2 b^2 + 1 }    \ ; \zeta_{n-1 } \ ; , \end{aligned}\ ] ] where @xmath272 \ ; .\ ] ] using definition ( [ np2a ] ) one can rewrite eq . ( [ np3 ] ) in terms of the parameter @xmath267 as @xmath273    eqs . ( [ np3],[np5 ] ) are valid provided ( [ t2 ] ) holds , _ i.e. _ , provided @xmath274 , which reads @xmath275 in the following we also impose the condition @xmath276 precisely , we neglect all the quantities of order @xmath277 . this is an important technical point . it facilitates the analysis by allowing one to get simple formulas as we now illustrate in the perturbative case : eq . ( [ np5 ] ) allows for instance to compute the average value of the reflection coefficient in the perturbative regime ( as already done in section [ pts ] , eq . ( [ spt10 ] ) ) . in this regime , additionally to condition ( [ np5b ] ) one has @xmath278 . then eq . ( [ np5 ] ) implies at leading order @xmath279 which , together with the initial condition @xmath280 , leads immediately to @xmath281 which is identical to result ( [ spt10 ] ) in the case of a potential @xmath7 for which @xmath282 [ cf . ( [ spt11 ] ) ] . let us now proceed and consider the generic non perturbative regime where @xmath283 may become large compared to unity and where ( [ t2 ] ) and ( [ np5b ] ) are still valid . taking the condition ( [ np4a ] ) into account , eq . ( [ np5 ] ) reads @xmath284 it is natural to assume that the phase of the cosine in the r.h.s . of ( [ np4b ] ) is uniformly distributed in @xmath285 $ ] and independent of the phase at step @xmath286 . this could be called a `` phase randomization '' approximation . this relies on hypothesis ( [ np4a ] ) and on the assumption that there is a large number of density oscillations over the ( random ) length between @xmath287 and @xmath288 , i.e. , @xmath289 then , the argument of the cosine in definition ( [ np4b ] ) is uniformly distributed , @xmath290 s are uncorrelated random variables , with all the same law characterized by its average @xmath291 and variance @xmath292 note that the regimes ( [ np4a ] ) and ( [ gnp2bis ] ) imply that @xmath293 which in turn implies that the kinetic energy @xmath294 is much larger than @xmath295 ; _ i.e. _ , one is exactly in the anderson regime where the incident kinetic energy is much larger than the typical value of the ( disordered ) potential representing the obstacle . hence a classical particle would flow almost unperturbed over the potential but , as we shall see , a quantum particle experiences an exponentially small transmission . let @xmath296 be the probability that @xmath267 lies in the interval @xmath122 , @xmath297 . going to the continuous limit and defining the continuous variable @xmath298 [ where @xmath299 is the parameter ( [ spt11 ] ) in the case of a potential @xmath7 ] it is shown in appendix [ app2 ] that @xmath301 verifies the following fokker - planck equation @xmath302\ ; .\ ] ] equation ( [ dmpk ] ) follows directly from eq . ( [ np6 ] ) in the regime where conditions ( [ np4a ] ) and ( [ gnp2bis ] ) hold . it is precisely the dorokhov - mello - pereyra - kumar ( dmpk ) equation @xcite for the transmission in a single disordered channel [ with @xmath303 . equation ( [ dmpk ] ) is sometimes referred to as melnikov s equation ( after ref . @xcite ) but has a much longer history ( see the discussion in refs . @xcite ) .    since before entering the disordered region the particle has a classical energy @xmath251 corresponding to @xmath304 , eq . ( [ dmpk ] ) has to be solved for the initial condition @xmath305 where @xmath306 is the one - sided delta function : @xmath307 . in the limit of small @xmath26 ( i.e. , in the perturbative regime @xmath308 ) , @xmath122 remains small and one can approximate in the r.h.s . of ( [ dmpk ] ) the term @xmath309 by @xmath122 . it is then simple to verify that the solution of this approximate equation that satisfies ( [ np8 ] ) is @xmath310 this result for the small @xmath26 solution of the dmpk equation has been already obtained in ref . @xcite ( see also the discussion in ref . the distribution law ( [ dmpk0 ] ) is exactly equivalent to distribution ( [ spt4 ] ) of the reflection coefficient in the perturbative regime and this proves the validity of the poissonian distribution ( [ e7b ] ) for a potential @xmath7 of type ( [ e2a ] ) @xcite .    in the general case ( i.e. , for all @xmath311 ) the solution of ( [ dmpk ] ) with the initial condition ( [ np8 ] ) is ( see , _ e.g. _ , refs . @xcite and references therein ) @xmath312 where @xmath313 .    from distribution ( [ np9 ] ) , a lengthy computation or alternatively the direct use of the dmpk equation ( [ dmpk ] ) @xcite yields @xmath314    in the large @xmath26 limit , distribution ( [ np9 ] ) tends to a log - normal distribution , i.e. , the distribution of the variable @xmath315 is gaussian ( see ref . @xcite ) @xmath316 from this distribution one gets the correct average @xmath317 , eq .  ( [ np10 ] ) , and a standard deviation @xmath318^{1/2 } = \sqrt{2 t}$ ] , which is in agreement with the exact result in the limit @xmath319 @xcite . at the extremity of a sample of length @xmath0 one has @xmath320 , and the distribution ( [ np11 ] ) is the log - normal distribution of transmission typical for anderson localization in the regime @xmath321 ( see , e.g. , ref . @xcite ) . as a side product of this analysis , eqs . ( [ np10 ] ) and ( [ np11 ] ) confirm that @xmath8 is indeed the localization length as was anticipated in the notation . we have tested the validity of the dmpk approach for a bose - einstein beam of interacting particles propagating in a disordered potential @xmath7 of type ( [ e2a ] ) . the numerical results for the probability distribution @xmath322 are compared on fig . [ fig7 ] with the dmpk prediction ( [ np9 ] ) . the agreement is seen to be excellent . the distribution evolves from the poissonian result ( [ e7b ] ) ( for low values of @xmath323 ) towards a distribution peaked at low @xmath214-values for large @xmath323 . in this latter case one can check that the distribution tends to a log - normal by plotting @xmath324 .     of type ( [ e2a ] ) ( characterized by @xmath325 and @xmath326 ) plotted for different values of the ratio @xmath327 with @xmath328 . the black solid lines are the dmpk result ( [ np9 ] ) and the colored histograms correspond to the numerical simulations ( 50000 samples used for each value of @xmath26 ) . cases ( a ) , ( b ) , ( c ) and ( d ) correspond respectively to @xmath329 , 0.5 , 1 and 2 . the inset displays @xmath330 as a function of @xmath26 . the thick solid line is extracted from numerical simulations and the thin ( red ) solid line is the dmpk prediction ( [ np10 ] ) . they can be distinguished only around @xmath331 as shown in the blowup of the ( gray ) shaded region for @xmath332 . [ fig7 ] ]    we have also checked the validity of the present approach over a sizable range of lengths of disordered region and of intensities of disordered potential by plotting in the inset of fig . [ fig7 ] the average @xmath333 as a function of @xmath323 . the agreement of the numerical results with the dmpk prediction ( [ np10 ] ) is excellent . note however the beginning of a small departure around @xmath334 ; this effect will be studied more thoroughly in sec . [ sec - lstar ] ( cf . [ fig9 ] ) . finally , we discuss numerical results obtained for the disordered potentials introduced in secs . [ gaussian ] and [ speckle ] . although we do not have an analytical derivation of the dmpk equation for these potentials , the numerical results indicate a very good quantitative agreement for a disordered potential @xmath173 and for a speckle potential @xmath176 . we display the comparison of the numerical data with the dmpk predictions for a speckle potential in figure [ fig8 ] . the same agreement is obtained for a gaussian potential @xmath173 . hence , the behavior analytically predicted for the potential @xmath7 appears to be of general validity , meaning that the above defined regime of `` phase randomization '' can probably be extended to correlated potentials , leading to a regime of single parameter scaling . however , we have noticed that , although showing an overall good agreement with the dmpk prediction , the lorentzian correlated potential @xmath172 exhibits some deviations in the tail of the distribution , the details of which will be studied elsewhere . ( characterized by @xmath335 and @xmath336 @xmath337 ) plotted for different values of the ratio @xmath327 with @xmath338 . curves ( a ) , ( b ) and ( c ) corresponds respectively to @xmath339 , 0.52 and 0.68 . for each curve , the black solid line is the dmpk result ( [ np9 ] ) and the colored histogram is the result of 10000 numerical simulations.[fig8 ] ]      in the previous sections the main effect of interaction has been shown to be a renormalization of the localization length @xmath8 . interaction induces a modification of the wave vector : expression ( [ spt11 ] ) for the localization length coincides with the noninteracting one but computed for an effective interaction - dependent wave vector @xmath193 given by eq . ( [ pt4 ] ) , instead of @xmath226 . the repulsive interaction diminishes the available kinetic energy and therefore reduces the localization length with respect to the noninteracting case ( since @xmath340 ) . we now discuss another , more spectacular , effect of interactions on the localization properties of a propagating bec on a disordered potential .    in the previous sections [ pts ] and [ non - pert ] , we completely neglected the presence of an upper limit for the classical energy @xmath341 , which is given by the local maximum of the effective potential @xmath91 , namely @xmath342 ( see fig . [ fig3 ] ) . trajectories that pass beyond @xmath343 would become unstable and develop singularities with infinitely large density at @xmath344 . in practice this implies , on the level of eq . ( [ e1 ] ) , that a stationary flow can not be maintained in this case and that the disorder induces time - dependent dynamics of the condensate .    in the vicinity of @xmath343 , the density profile of the condensate in between two adjacent scatterers becomes quite different from the cosine shape ( [ t6 ] ) that was derived for weak nonlinearities and/or low density modulations , and resembles more to a periodic train of gray solitons @xcite . in a crude approximation , we neglect this complication and assume that the spatial evolution of the density is still given by eq . ( [ t6 ] ) for all classical energies until @xmath345 . trajectories that happen to pass beyond @xmath343 are considered to be `` lost '' , i.e. , they do no longer contribute to the probability distribution for the transmission . this formally amounts to introducing a `` sink '' in the stochastic equation ( [ np6 ] ) , namely at @xmath346 . in the corresponding fokker - planck equation ( [ dmpk ] ) , this sink is appropriately modeled by imposing the boundary condition @xmath347 as a consequence of this boundary condition , the integrated probability distribution @xmath348 is no longer conserved , but decreases with increasing @xmath26 , i.e. , increasing length @xmath0 of the disorder region .    in the following , we show how this affects the dmpk predictions of section [ non - pert ] and how the `` survival probability '' , i.e. , the fraction of trajectories that remain below this boundary at given length @xmath0 , can be analytically computed in the limit @xmath107 . in this limit , from eq . ( [ t1c ] ) one gets @xmath349 and thus @xmath350 . modifications of the probability density @xmath301 due to the presence of the sink appear only when the typical value of @xmath122 is of order @xmath351 which , as just remarked , is large compared to unity in the case @xmath107 . in this case @xmath301 is already negligibly small around @xmath352 . we therefore make the approximation @xmath353 in the fokker - planck equation ( [ dmpk ] ) , which then reads @xmath354\ ; .\ ] ] using , from now on , the probability distribution @xmath355 for finding a given value of @xmath315 at fixed @xmath356 , we obtain in this limit @xmath357 where we introduce @xmath358 . clearly , the log - normal distribution ( [ np11 ] ) corresponds to a solution of eq . ( [ lognorm - eq ] ) in the absence of any additional boundaries .    in the presence of the sink , which is imposed by the boundary condition @xmath359 with @xmath360 we can straightforwardly find the solution of eq . ( [ lognorm - eq ] ) by subtracting from the log - normal distribution ( [ np11 ] ) a `` mirror '' distribution centered at some @xmath361 ( namely @xmath362 ) with a suitable prefactor . this yields the distribution @xmath363\ ; ,   \label{lognorm - decay}\end{aligned}\ ] ] which is defined for @xmath364 . clearly , this distribution satisfies the evolution equation ( [ lognorm - eq ] ) as well as the boundary condition @xmath359 for all @xmath26 and the initial condition @xmath365 for @xmath364 . the presence of the sink at @xmath366 explains a phenomenon barely noticeable in fig . 7 , but exemplified in fig . [ fig9 ] , namely the departure of the observed average @xmath367 from the usual dmpk result @xmath368 . this departure is due to the fact that the numerically computed average only takes into account the stationary solutions which as will be seen from eq . ( [ psurv]) become less and less numerous when @xmath26 increases . hence what is computed numerically is the average of @xmath369 over the distribution ( [ lognorm - decay ] ) . this reads @xmath370 -   \nonumber \\ & & e^{z_{\rm max } }   \left(\frac{t}{2}+ z_{\rm max}\right)\ , \textrm{erfc } \left(\frac{t+z_{\rm max}}{2\sqrt{t}}\right ) \ ; .\label{zmoy_leak}\end{aligned}\ ] ] where the error function is defined by @xmath371 and @xmath372 . expression ( [ zmoy_leak ] ) is compared in fig . [ fig9 ] with the results of a numerical simulation performed in the case @xmath373 ( corresponding to @xmath374 ) for 10000 random potentials @xmath7 of type ( [ e2a ] ) characterized by @xmath326 and @xmath375 ( leading to @xmath376 ) . the agreement is seen to be very good . since the sink cuts the solutions which are strongly scattered by the random potential , the remaining stationary states have a higher transmission coefficient . this effect increases with the sample length @xmath0 , which explains the behavior of the curve in fig . [ fig9 ] .     plotted as a function of @xmath327 in the case of random potentials @xmath7 characterized by @xmath326 and @xmath375 . the curve is drawn in the case @xmath373 . the red solid line is the numerical result and the black dashed line is the analytical result ( [ zmoy_leak ] ) . the straight ( thin dashed ) line is the usual dmpk result @xmath377 [ eq . ( [ np10])].[fig9 ] ]    as an other test of the validity of our approach ( which amounts to model the upper boundary @xmath378 by a perfect sink and to neglect nonlinear deformations of the density pattern of the flow close to the threshold ) we now determine the probability for a trajectory to remain below the boundary . this survival probability reads @xmath379 \nonumber \\ & & - \frac{e^{z_{\rm max}}}{2 } \ , \mathrm{erfc } \left(\frac{z_{\rm max } + t}{\sqrt{4 t}}\right ) \ ; .\label{psurv}\end{aligned}\ ] ] as anticipated , @xmath380 clearly decreases from @xmath381 ( at @xmath382 ) to 0 ( for large @xmath26 ) . the knowledge of @xmath380 allows to determine the value @xmath1 of the length of the disordered region beyond which most of the random realizations lead to a non - stationary flow of the condensate . we can , most conveniently , define @xmath1 through the condition @xmath383 with @xmath384 . this leads to the implicit equation for the threshold value @xmath385 : @xmath386 this equation can be explicitly solved in the limiting case of large @xmath387 . as it is natural to assume that @xmath385 ought to be of the order of @xmath378 , which is the only relevant scale in this equation , we make the ansatz @xmath388 and assume ( which is to be verified _ a posteriori _ ) that @xmath389 is of the order of unity , whereas @xmath390 . this yields to lowest non vanishing order @xmath391 \ ; , \ ] ] for the left - hand side of eq . ( [ tstar - impl ] ) and @xmath392 \nonumber \\ & = & \frac{1}{\sqrt{\pi\ , z_{\rm max } } } \left [ 1 + \mathcal{o}\left(z_{\rm max}^{-1}\right ) \right ] \ ; , \label{tstar - right}\end{aligned}\ ] ] for the right - hand side of eq . ( [ tstar - impl ] ) . this finally results in @xmath393 neglecting terms of the order of @xmath394 , we therefore obtain for the threshold length @xmath395 \ .\ ] ] we emphasize that eq . ( [ lstar ] ) holds for @xmath396 , i.e. , for @xmath397 [ see eq . ( [ defzmax ] ) ] . this is much more restrictive than the condition @xmath107 which is assumed to hold true when deriving eqs . ( [ lognorm - decay ] ) and ( [ tstar - impl ] ) .     plotted as a function of the length @xmath0 of the disordered region . the condensate flow through disorder potentials @xmath7 of type ( [ e2a ] ) was numerically computed for this purpose ( red solid line ) , at parameters for which @xmath398 . the black dashed line shows the analytical prediction of this survival probability @xmath399 according to eq . ( [ psurv ] ) . the vertical dashed line marks prediction ( [ lstar ] ) for the threshold length @xmath1 at which @xmath400 ( horizontal dashed line ) , namely @xmath401  .[fig10 ] ]    figure [ fig10 ] shows a comparison of the analytical predictions ( [ psurv ] ) and ( [ lstar ] ) with numerical data obtained from the integration of the time - dependent gross - pitaevskii equation ( [ e1 ] ) . the condensate flows through a disorder potential @xmath7 of type ( [ e2a ] ) with @xmath398 . we see that the fraction of stationary trajectories @xmath380 is very well described by eq . ( [ psurv ] ) , and that the approximate expression ( [ lstar ] ) predicts very well the length @xmath1 at which the crossover length from stationary to time - dependent flow occurs .    for velocities @xmath3 not extremely large compared to the speed of sound , the condition @xmath396 will not be fulfilled and estimate ( [ lstar ] ) will not be valid , while @xmath402 might still hold and the average evolution of the system might still be fairly well described by the simplified fokker - planck equation ( [ dmpk - large ] ) . in that case , the implicit equation ( [ tstar - impl ] ) has to be solved numerically . in figure [ fig1 ] one can see that the numerical solution of eq . ( [ tstar - impl ] ) ( yellow solid line ) provides a very reasonable estimate of the boundary between the bright supersonic region ( stationary flows ) and the dark time - dependent region , in a regime of not extremely large @xmath403 , where eq . ( [ lstar ] ) fails to properly predict the threshold length @xmath1 . the simulations are performed by solving eq . ( [ e1 ] ) numerically using a potential of type @xmath7 ( characterized by @xmath404 and @xmath405 ) . for each @xmath3 and @xmath0 we consider 100 realizations of such a potential and statistically determine the quantity @xmath406 , i.e. , the fraction of stationary solutions . @xmath406 is plotted in fig . [ fig1 ] using a gray scale ( dark , @xmath407 ; light blue / gray , @xmath408 ) as a function of the normalized variables @xmath409 and @xmath403 ( this normalization rescales interaction effects ) . the qualitative agreement of figure [ fig1 ] is made quantitative in figure [ fig11 ] . in this figure the numerical solution of eq . ( [ tstar - impl ] ) is compared with its determination extracted from numerical simulations in the supersonic regime . more precisely , the solid ( red ) line in fig . [ fig11 ] is simply the contour @xmath410 in fig . this corresponds exactly to definition ( [ tstar ] ) of @xmath1 . the agreement between the numerical result and the theory of the present section [ dashed curve , solution of eq . ( [ tstar - impl ] ) ] is seen to be excellent @xcite .     as a function of @xmath3 in dimensionless units . the black dashed line corresponds to the solution of eq . ( [ tstar - impl ] ) and the red solid line corresponds to the value of @xmath1 extracted from numerical simulations for a potential @xmath7 with the same characteristics as in fig . [ fig1 ] ( see the text).[fig11 ] ]    we conclude this section by emphasizing that the existence of an upper threshold @xmath1 corresponding to lengths of the disordered region beyond which most of the flows are time - dependent is a genuine nonlinear effect [ absent if one sets @xmath411 in eq . ( [ e1 ] ) ] . actually , whereas interactions only weakly modify the precise value of the localization length , the existence of the threshold @xmath1 is a remarkable qualitative effect induced by nonlinearity . moreover , as illustrated in figures [ fig9 ] and [ fig10 ] , this effect persists even in the limit @xmath107 where naively one would expect no noticeable consequence of interaction . on the basis of the results obtained in the previous section we present here what are the more favorable experimental configurations for observing anderson localization in an interacting bose - einstein beam . we also discuss a possible experimental signature of localization .      in the non - interacting regime the only condition for observing anderson localization in 1d is that the size of the disordered region should be larger than the localization length . then one can observe an exponential decay of transmission with a log - normal distribution ( in the limit @xmath321 ) . the situation is more complex when interactions are turned on . what is particularly interesting is the interplay between localization and superfluidity . indeed these two phenomena are conflicting one with the other : superfluidity is the ( counterintuitive ) ability to pass over an obstacle without reflection whereas anderson localization corresponds to a large reflection in a domain where one would expect almost perfect transmission . as a result of the interplay between these two extreme phenomena , and depending on the fluid velocity and on the sample size , the flow may be stationary and superfluid , dissipative and time dependent or stationary supersonic ( and also dissipative ) @xcite . anderson localization does not occur in the superfluid region ( where the transmission is perfect ) and either does not exist or can not be clearly identified in the time dependent regime ( where interference effects are washed out @xcite ) ; but is truly observed in the supersonic stationary regime , as demonstrated in sec . [ ssr ] .    in that regime , a first experimentally relevant effect is the modification of the localization length with respect to its value in the absence of interactions . this effect is very well described by renormalizing the wave vector @xmath224 to @xmath193 [ eq .  ( [ spt11 ] ) ] , which means that part of the kinetic energy available to the flow is taken by interactions . however , as already discussed in ref . @xcite this effect is only sizable in a regime where @xmath3 is not too large compared to @xmath2 , and is thus relevant only in the perturbative regime ( cf . [ fig5 ] ) . a second experimentally observable effect is the modification of the localization length due to the correlations of the disordered potential . this is described by formula ( [ spt11 ] ) where @xmath206 is the fourier transform of the two - point correlation function of the disorder . for the different potentials considered here , @xmath412 for a potential @xmath7 or is alternatively given by eqs . ( [ e2j ] ) , ( [ e2k ] ) and ( [ csq ] ) for correlated potentials . explicitly this yields @xmath413 for a potential @xmath7 of type ( [ e2a ] ) ; @xmath414 for a potential @xmath173 of type ( [ e2d ] ) ; @xmath415 for a potential @xmath172 of type ( [ e2d ] ) ; and @xmath416 for a potential @xmath176 of type ( [ e2l ] ) ( when @xmath417 ) . the validity of these expressions has been tested in sec . [ non - pert ] . in the non - interacting case ( i.e. , @xmath223 ) , expressions ( [ exp1 ] ) to ( [ exp4 ] ) correspond to a high energy limit and can be obtained through a first order born expansion within the phase formalism of refs . @xcite . in all three cases , one sees that the localization length is drastically enhanced due to the non - zero correlation length with respect to the uncorrelated disorder , eq . ( [ exp1 ] ) . in the gaussian and the lorentzian cases the localization length scales exponentially with @xmath418 and @xmath419 , respectively [ see eqs . ( [ exp2],[exp3 ] ) ] . in the case of a speckle potential , the effect is even stronger : one sometimes speaks of an `` effective mobility edge '' @xcite , meaning that beyond a critical wave - vector ( or a critical velocity ) the localization length ( [ exp4 ] ) is infinite . this is an artifact of the antsygina - pastur - slyusarev formula ( [ spt11 ] ) which can be corrected by going to higher orders ( see refs . @xcite ) : the corrections to this result give a localization length which is finite , but typically larger than any other relevant scale in experimental systems .    hence , in all the cases the dependence of the localization length with respect @xmath193 ( i.e. , with velocity ) is amplified by correlations . mathematically this is due to the fact that the denominator in the antsygina - pastur - slyusarev formula ( [ spt11 ] ) for the localization length in presence of correlations tends to zero when @xmath420 . in order to minimize this effect one needs to impose the following condition : @xmath421 in the r.h.s . of eq . ( [ exp5 ] ) we replaced @xmath193 by @xmath422 because in practice condition ( [ exp5 ] ) is verified in regimes where @xmath423 , i.e. , when the approximation @xmath424 is sound . note that this condition is arbitrary and is only superficially analogous to the 3d ioffe - regel criterion @xcite . the latter defines a true mobility edge that separates a metallic from a localized phase whereas eq . ( [ exp5 ] ) only requires that the localization length does not get too large . understood in this sense , the criterion ( [ exp5 ] ) is exactly equivalent to the definition of an `` effective mobility edge '' sometimes used in the literature .    in the absence of interactions it is always possible ( at least theoretically ) to define a system with a length @xmath425 which verifies ( [ exp5 ] ) ; i.e. , a system where one can observe anderson localization . if we now turn on interactions , a major effect is the appearance of a length scale @xmath1 which signals the onset , for @xmath426 , of a regime of time - dependent flows ( cf . [ sec - lstar ] ) . in this regime , anderson localization disappears , and the time - averaged transmission coefficient scales as @xmath427 @xcite . this is the most spectacular effect of interactions in the transport properties of the system . in order to observe anderson localization , the system size should therefore satisfy @xmath428 . in practice , one should be in a regime of parameters such as illustrated in fig . [ fig12 ] : the crossing @xmath429 has to occur at a velocity lower than @xmath430 .     and the solid line is the threshold length @xmath1 . @xmath431 is the typical velocity beyond which it is almost impossible to observe anderson localization in a realistic system ( see text ) . the ( blue ) colored zone corresponds to the region where anderson localization can be experimentally observed in presence of interaction.[fig12 ] ]    based on the numerical solution of eq . ( [ tstar - impl ] ) one can show that the crossing @xmath432 occurs at a velocity @xmath433 @xcite ( see also fig . [ fig1 ] ) . this condition only allows the system to reach a ( stationary ) regime where @xmath434 . but if one wants to observe anderson localization one should be able to reach a regime where @xmath435 say , in order to get as close as possible to the domain of log - normal distribution of transmissions still remaining in the region of stationary flows . this imposes @xmath436 . this must be supplemented by condition ( [ exp5 ] ) , i.e. @xmath437 . hence the correlation length @xmath169 should be smaller or equal to @xmath438 . fig . [ fig13 ] shows the phase diagrams of a one dimensional interacting beam of condensed atoms moving through a speckle potential in this regime . for plotting this diagram one has generated 16 random potentials and studied in each case if a stationary solution exists or not . the dark blue region corresponds to a domain where no stationary solution exists while the light blue one corresponds to a domain where all the potentials admit a stationary solution ( the color code is the same as in fig . [ fig1 ] and is explained in sec . [ sec - lstar ] ) . the region between @xmath8 and @xmath1 in figure [ fig13 ] is the region where one can observe anderson localization .     for a beam with velocity @xmath3 moving in a speckle disorder of extent @xmath0 . the figure has been drawn for a potential @xmath176 of type ( [ e2l ] ) characterized by @xmath335 and @xmath439 . the light blue region corresponds to a domain of stationary flow ( @xmath408 : 100 @xmath232 of the solutions are stationary , see the explanation in the text ) ; the dark region corresponds to time dependent flow . the curves indicating the values of @xmath8 and @xmath1 correspond to the analytical results ( [ exp4 ] ) and to the numerical solution of eq . ( [ tstar - impl]).[fig13 ] ]    let us now evaluate the orders of magnitude of the different parameters allowing one to reach the appropriate regime . for concreteness we consider a beam of @xmath41rb atoms such as the one of the atom optics group at laboratoire charles fabry de linstitut doptique . for a correlation length of 0.26 @xmath440 m the velocity cut - off for observing anderson localization is roughly @xmath441 mm / s . note that in ref . @xcite the velocity of the expanding condensate is about 1.6 mm / s , i.e. , smaller that @xmath430 as it should . if we use the parameters of ref . @xcite ( @xmath442 , @xmath443 mm / s , @xmath444 nm , @xmath445 atoms/@xmath440 m , @xmath446 mm / s and @xmath447 @xmath440 m ) it is impossible to satisfy condition ( [ exp5 ] ) because @xmath443 mm / s @xmath448 and also because @xmath449 . however the atom optics group has recently improved the sensitivity of its detectors which is now close to being able to detect a density as small as 1 atom/@xmath440 m . this allows one to work with a smaller density and to improve the ratio @xmath450 which can be tuned down to the value of @xmath451 . then the localization length can be selected by tuning the speckle amplitude . for instance @xmath452 mm can be obtained by choosing @xmath453 hz at @xmath454 mm / s . these parameters are close to those used in fig . [ fig13 ] and are reachable experimentally . however , for observing anderson localization one needs to keep the beam stable for almost 1 s ( 0.31 s if we want @xmath455 ) whereas in the current experiment this is only possible during 0.1 s , hence it is still a matter of debate to decide if the observation of anderson localization of a bose condensed beam in the presence of interaction is within the reach of present time technology . it is also interesting to make a connection between the physics described here and the recent experiment observation of anderson localization of a condensate expanding in a disordered potential performed in the same group @xcite . contrarily to the propagation of a beam studied in the present work , ref . @xcite considers the spreading of a wave packet ( initially at rest ) in a speckle potential . after a first stage of expansion , mainly driven by interactions , the experimental cloud expands with a constant velocity @xmath456 mm / s but the particle density and the sound velocity are functions of the position . therefore it is not possible to place this experiment on a single point of the phase diagram displayed in fig . . however one can evaluate the ratio @xmath403 at different positions , that is at different @xmath0 . for instance , at the typical value @xmath434 , fig . 2 of reference @xcite allows to calculate the sound velocity as well as the healing length @xmath457 , yielding the typical experimental values @xmath458 and @xmath459 . moreover the ratio between the typical disorder amplitude @xmath460 and the chemical potential @xmath461 is the same in fig . [ fig13 ] and in experiment : @xmath462 ( note however that in the experimental case this is the value of the _ local _ chemical potential that matters ) . hence , although the experimental setup of ref . @xcite forbids a direct comparison with the results of the present work , the estimates of the typical values @xmath458 and @xmath459 indeed locate the experimental system within the regime of anderson localization of fig . [ fig13 ] . once the appropriate regime of parameters for observing anderson localization in a bose condensed beam has been determined , it is also important to identify possible experimental signatures . in our theoretical approach we use the transmission coefficient @xmath214 of the beam over the disordered region as the relevant parameter . however the measure of @xmath214 might be experimentally involved , and we propose here to use an other related quantity , namely the rate of energy dissipation @xcite : @xmath463 where @xmath464 is the drag force exerted by the beam on the obstacle . definition ( [ exp6 ] ) is quite natural : the force exerted on the obstacle is the mean value of the operator @xmath465 over the condensate wave function . it is rigorously justified by the analysis of ref . @xcite in terms of stress tensor . in the stationary case , changing integration from @xmath17 to @xmath86 in ( [ exp6 ] ) , a simple integration by parts yields @xmath466 in the r.h.s . of eq . ( [ exp10 ] ) we made use of relations ( [ e4 ] ) and ( [ t1 ] ) . it has been shown in section [ transmis ] that @xmath81 and that under assumption ( [ t2 ] ) ( small nonlinearity and arbitrary transmission or weak transmission and arbitrary nonlinearity ) one has ( see eqs . ( [ t7 ] ) and ( [ t5 ] ) ) @xmath467 which yields @xmath468 in the regime @xmath469 , using eq . ( [ spt3 ] ) one recovers from ( [ exp11 ] ) the perturbative result already obtained in ref . @xcite : @xmath470 ) is rather simple and it is worth spending some time to discuss it . consider an incident beam of particles with density @xmath471 and momentum @xmath472 moving from @xmath83 towards an obstacle at rest . part of the particles is transmitted ( a fraction @xmath214 ) and the other part is reflected ( a fraction @xmath473 ) . the collisions are elastic and each of the reflected particles experiences an exchange of momentum @xmath474 with the obstacle . during a time @xmath389 there are @xmath475 collisions and by the law of action and reaction the obstacle experiences a force @xmath476 in the r.h.s . of ( [ exp12 ] ) , one has written that @xmath477 is the flux of particles colliding with the obstacle i.e. , @xmath478 . ( [ exp11 ] ) and ( [ exp12 ] ) are identical because what we call @xmath29 is the downstream density of the beam ( cf . [ fig2 ] ) , i.e. , precisely @xmath479 . depending on which quantity is held constant ( @xmath29 as in the present paper , or @xmath471 ) eq . ( [ exp11 ] ) or ( [ exp12 ] ) is more appropriate ( cf . the discussion of the fixed input and fixed output problem in ref . this is somewhat reminiscent of the controversy on the correctness of the landauer formula ( see , e.g. , the discussion in ref . @xcite ) .    on the basis of ( [ exp11 ] ) one sees that the measure of @xmath480 gives direct informations on the transmission of the interacting beam through the disordered region , allowing one to reveal in which configuration is the system . for instance in the localized regime the energy dissipation is high ( @xmath481 ) and grows exponentially with the size @xmath0 of the disordered region , whereas in the perturbative regime @xmath480 is much lower and scales as @xmath0 . in the present work we have presented an analysis of the transmission of a weakly interacting bose gas incident on a disordered potential . we have shown on the basis of numerical and analytical results that there is a regime of anderson localization in this system and proposed experimental signature of this phenomenon . in order to properly identify a `` localized regime '' we have studied the transmission coefficient and its probability distribution . the transmission coefficient @xmath214 is well defined under assumption ( [ t2 ] ) , which holds in the following regimes : ( i ) small nonlinearity and arbitrary transmission or ( ii ) weak transmission and arbitrary nonlinearity . in other cases there is no obvious way to define the transmission of the non - linear beam because one can not separate in the up - stream region an incident flow from a reflected one . however , our analysis in terms of @xmath341 and @xmath122 ( defined in sec . [ transmis ] ) is always valid , even when condition ( [ t2 ] ) is not fulfilled . this just means that , out of regime ( [ t2 ] ) , the connection ( [ t7 ] ) between @xmath122 and @xmath214 is invalid . but for instance this does not invalidate at all the analysis leading to the dmpk equation ( [ dmpk ] ) , and the experimental signature proposed in sec . [ signature ] also remains valid even when it is not possible to properly define @xmath214 . we note that the validity of the dmpk approach for non - interacting particles is a well established fact in the theory of disordered systems . what is achieved in the present work is its extension to the case of interacting particles . other studies of anderson localization in the presence of interactions have concentrated on the long time behavior of the time evolution of initial wave packets @xcite . although those results are still a matter of active debate in the community , the results of the present work produce strong evidence of the existence of anderson localization for weakly interacting bose particles ( with effective repulsive interaction ) propagating through disordered samples of finite size @xmath482 . although the present study leads to the important conclusion that anderson localization in the presence of interaction is possible , it is rather disappointing to remark that it can be clearly identified only when @xmath483 , i.e. , in a regime where interactions do not play a major role ( see the discussion of sec . [ apconf ] and also ref . @xcite ) . in this respect , the more interesting and new effect of interactions is the existence of an upper threshold @xmath1 for the length of the disordered region  : for @xmath426 no stationary flow is possible . as shown in sec . [ sec - lstar ] , @xmath1 is directly connected to the probability distribution of the parameter @xmath122 and to the localization properties of the system . it would be very interesting to lead a systematic study of the transmission in the interaction - induced time dependent regime ( @xmath426 ) where the numerical results of ref . @xcite indicate a power law decay of the time - averaged transmission , a signature generally considered of loss of phase coherence and onset of ohmic behavior @xcite . work in this direction is in progress . it is a pleasure to thank b. altshuler , a. comtet , j .- l . pichard and c. texier for inspiring discussions . this work was supported by anr grants no . 05nano00802 , no . nt05242103 and no . 08blan016501 , by the ifraf institute and by the alexander von humboldt foundation . we gratefully acknowledge funding by the excellence initiative of the german research foundation ( dfg ) through the heidelberg graduate school of fundamental physics ( grant no gsc 129/1 ) and the global networks mobility measures the frontier innovation fund of the university of heidelberg . we furthermore acknowledge support through the dfg forschergruppe 760 `` scattering systems with complex dynamics '' . we give here a demonstration of the perturbative results ( [ spt4 ] ) , ( [ spt10 ] ) and ( [ spt11 ] ) in the special case of a gaussian disordered potential @xmath153 of type ( [ e2d ] ) . a simple way to obtain this results starts by noticing that any gaussian noise verifying @xmath484 and @xmath485 [ and here we are specifically interested in @xmath147 that appears in eq . ( [ e2d ] ) ] can be written as ( see , e.g. , ref . @xcite ) @xmath486 where @xmath487 s are random positions uniformly distributed on the real axis with density @xmath488 and mean spacing @xmath489 and @xmath490 is a random variable ( with @xmath491 and @xmath492 ) .    in order to calculate the probability distribution of the reflection coefficient @xmath473 [ whose value is given by ( [ spt3 ] ) ] one should first consider the distribution of @xmath493 the quantity @xmath494 as given by ( [ a2 ] ) is formally equivalent to the position @xmath495 of a particle performing a random walk in the complex plane after @xmath496 iterations . the particle is initially at the origin and performs jumps of constant amplitude @xmath497 with random direction . denoting by @xmath498 the probability to find the particle in the domain @xmath499 around @xmath495 after @xmath135 steps , if @xmath500 [ which is ensured by taking the limit @xmath501 in ( [ a2 ] ) ] , the central limit theorem yields @xmath502 it is then a simple exercise to get the distribution of @xmath503 . one obtains @xmath504 \ ; , \ ] ] where @xmath505 from relation ( [ spt3 ] ) , eq .  ( [ a5 ] ) immediately yields the announced probability distribution ( [ spt4 ] ) with @xmath506 for a potential @xmath153 of type ( [ e2d ] ) , @xmath507 [ see eq . ( [ e2 g ] ) ] and eq . ( [ a6 ] ) demonstrates in this case the validity of eqs . ( [ spt10 ] ) and ( [ spt11 ] ) . in this appendix , we explain how to obtain the dmpk equation ( [ dmpk ] ) starting from the discrete langevin equation ( [ np6 ] ) . let us consider a generic situation where @xmath267 obeys a stochastic recursion relation of the type @xmath508 with uncorrelated random variables @xmath290  : @xmath509 it is clear that under assumption ( [ gnp2bis ] ) , eq . ( [ np6 ] ) is of type ( [ b1 ] ) with all the odd @xmath135 averages in ( [ b2 ] ) being zero and @xmath510 [ cf . ( [ varxi ] ) ] . let @xmath511 be the probability that @xmath267 lies in the interval @xmath122 , @xmath297 . one can express @xmath512 as @xmath513 this yields @xmath514 using the fact that @xmath267 depends on the variables @xmath515 , @xmath516  @xmath517 but not on @xmath290 ( as can be seen directly from ( [ b1 ] ) ) one can write the last of eqs . ( [ b4 ] ) as @xmath518 in the case of eq . ( [ np6 ] ) one has @xmath519 and the successive moments of @xmath520 read @xmath521 @xmath522 with all the other moments being of order @xmath277 or more , i.e. , negligible in regime ( [ np4a ] ) . ( [ b5 ] ) thus reads @xmath523 =   \nonumber \\ & & \nonumber \\ & & \displaystyle -\frac{\partial}{\partial\lambda}\left[(1 + 2\lambda)p\right ] +   \frac{\partial^2}{\partial\lambda^2}\left[(\lambda^2+\lambda)p\right ] =   \nonumber \\ & & \nonumber \\   & & \displaystyle \frac{\partial}{\partial\lambda } \left[\lambda(\lambda+1)\frac{\partial p}{\partial\lambda}\right ]    \ ; .\end{aligned}\ ] ] in the continuous limit , defining @xmath524 , the l.h.s . of eq . ( [ b8 ] ) is simply the first derivative of @xmath525 with respect to @xmath26 and ( [ b8 ] ) reduces to eq . ( [ dmpk ] ) of the main text . 99 p. w. anderson , phys . rev . * 109 * , 1492 ( 1958 ) . et al . _ , nature * 453 * , 891 ( 2008 ) . et al . _ , nature * 453 * , 895 ( 2008 ) . k. huang and h. f. meng , phys . lett . * 69 * , 644 ( 1992 ) . s. giorgini , l. pitaevskii , and s. stringari , phys . b * 49 * , 12938 ( 1994 ) . g. e. astrakharchik , j. boronat , j. casulleras , and s. giorgini , phys . a * 66 * , 023603 ( 2002 ) . d. clment , ph . bouyer , a. aspect , and l. sanchez - palencia , phys . a * 77 * , 033631 ( 2008 ) . yong p. chen . _ et al . _ , a * 77 * , 033632 ( 2008 ) . t. giamarchi and h. j. schulz , phys . b * 37 * , 325 ( 1988 ) . m. p. a. fisher , p. b. weichman , g. grinstein , and d. s. fisher , phys . rev . b * 40 * , 546 ( 1989 ) . s. a. gredeskul and y. s. kivshar , phys . rep . * 216 * , 1 ( 1992 ) . n. bilas and n. pavloff , phys . lett . * 95 * , 130403 ( 2005 ) . t. paul , p. leboeuf , n. pavloff , k. richter , and p. schlagheck , phys . rev . a * 72 * , 063621 ( 2005 ) . t. paul , p. schlagheck , p. leboeuf and n. pavloff , phys . * 98 * , 210602 ( 2007 ) . c. menotti and s. stringari , phys . a * 66 * , 043610 ( 2002 ) . d. s. petrov , g. v. shlyapnikov , and j. t. m. walraven , phys . * 85 * , 3745 ( 2000 ) . c. raman _ et al . _ , phys . rev . lett . * 83 * , 2502 ( 1999 ) . _ et al._. phys . lett . * 85 * , 2228 ( 2000 ) . p. engels and c. atherton , phys . lett . * 99 * , 160405 ( 2007 ) . m. olshanii , phys . lett . * 81 * , 938 ( 1998 ) . a. d. jackson , g. m. kavoulakis , and c. j. pethick , phys . a * 58 * , 2417 ( 1998 ) . p. leboeuf and n. pavloff , phys . a * 64 * , 033602 ( 2001 ) . l. p. pitaevskii and s. stringari , j. low temp . physics * 85 * , 377 ( 1991 ) . m. schwartz , phys . b * 15 * , 1399 ( 1977 ) . d. s. petrov , ph.d . thesis , university of amsterdam , 2003 , available online at http://www-old.amolf.nl/publications/theses/petrow  . a. v. lopatin and v. m. vinokur , phys . lett . * 88 * , 235503 ( 2002 ) . g. e. astrakharchik and l. p. pitaevskii , phys . a * 70 * , 013608 ( 2004 ) . t. ernst , t. paul , and p. schlagheck , arxiv:0905.4750 ( submitted to phys . a ) . v. hakim , phys . e * 55 * , 2835 ( 1997 ) . h. lamb , _ hydrodynamics _ , ( cambridge university press , cambridge , 1997 ) . p. leboeuf , n. pavloff , and s. sinha , phys . a * 68 * , 063608 ( 2003 ) . t. paul , m. hartung , k. richter , and p. schlagheck , phys . a * 76 * , 063605 ( 2007 ) . b. i. ivlev and n. b. kopnin , adv 33 * , 47 ( 1984 ) . j. s. langer and v. ambegaokar , phys . rev . * 164 * , 498 ( 1967 ) . note in particular that condition ( [ t2 ] ) is always fulfilled in absence of interaction since in this case the speed of sound is @xmath526 . b. i. shklovskii and a. l. efros , _ electronic properties of doped semiconductors _ , ( springer - verlag , 1984 ) . j. fortgh and c. zimmermann , rev . phys . * 79 * , 235 ( 2007 ) . s. kraft , a. gnther , h. ott , d. wharam , c. zimmermann , and j. fortgh , j. phys . b * 35 * , l469 ( 2002 ) . j.  estve _ et al . a * 70 * , 043629 ( 2004 ) ; t.  schumm _ et al . _ j. d * 32 * , 171 ( 2005 ) . wang , m. d. lukin , and e. demler , phys . rev . lett . * 92 * , 076802 ( 2004 ) . j w goodman , _ speckle phenomena in optics : theory and applications _ , roberts and compagny publishers ( greenwood village , 2007 ) . d. clment _ et al . _ , new j. phys . * 8 * , 165 ( 2006 ) . l. fallani , c. fort , and m. inguscio , adv . at . . phys . * 56 * , 119 ( 2008 ) . ( [ e2l ] ) corresponds to a typical physical situation where the potential seen by the particles is proportional to the intensity of an electric field whose components are random gaussian variables . in ref . @xcite this region has been improperly called `` ohmic '' . t. n. antsygina , l. a. pastur and v. a. slyusarev , sov . j. low temp . phys . * 7 * , 1 ( 1981 ) [ fiz . nizk . temp . * 7 * , 5 ( 1981 ) ] . i. m. lifshits , s. a. gredeskul , and l. a. pastur , _ introduction to the theory of disordered systems _ , john wiley ( new - york 1988 ) . o. n. dorokhov , jetp lett . * 36 * , 318 ( 1982 ) [ pisma zh . fiz . * 36 * , 259 ( 1982 ) ] ; p. a. mello , p. pereyra , and n. kumar , ann . ( n.y . ) * 181 * , 290 ( 1988 ) . v. i. melnikov , sov . solid state * 23 * , 444 ( 1981 ) [ fis . tela ( leningrad ) * 23 * , 782 ( 1981 ) ] . c. w. j. beenakker , rev . * 69 * , 731 ( 1997 ) . p. a. mello , j. math . phys . * 27 * , 2876 ( 1986 ) . c. w. j. beenakker and j. a. melsen , phys . b * 50 * , 2450 ( 1994 ) . in order to complete the proof of this equivalence one has to show that the mean value @xmath527 obtained from ( [ dmpk0 ] ) for a sample of length @xmath528 is identical to the result ( [ spt5 ] ) . this amounts to show that for a potential of the form ( [ e2a ] ) one has @xmath529 , which is easily obtained . a. a. abrikosov , solid state comm . * 37 * , 997 ( 1981 ) . b. a. van tiggelen in _ diffusive waves in complex media _ , j. p. fouque ed . , p. 1 ( kluwer academic publishers , dordrecht , 1999 ) . note that in ref . @xcite the theoretical determination of @xmath1 was obtained using a criterion different from ( [ tstar ] ) resulting in a less precise estimate at high @xmath403 . f. m. izrailev and a. a. krokhin , phys . lett . * 82 * , 4062 ( 1999 ) . l. sanchez - palencia _ et al . _ , phys . rev . lett . * 98 * , 210401 ( 2007 ) . l. tessieri , j. phys . a * 35 * , 9585 ( 2002 ) . e. gurevich and o. kenneth , phys . a * 79 * , 063617 ( 2009 ) . a * 80 * , 023605 ( 2009 ) . a. f. ioffe and a. r. regel , prog . semicond . * 4 * , 237 ( 1960 ) . note that the value @xmath530 does not depend on the type of correlation of the disordered potential considered since it is simply the value of @xmath403 for which eq . ( [ tstar - impl ] ) admits the solution @xmath531 . w. gurin _ et al . _ , lett . * 97 * , 200402 ( 2006 ) . n. pavloff , phys . a * 66 * , 013610 ( 2002 ) . note that there is a factor of 2 difference between the present definition ( [ pt4 ] ) of @xmath193 and the one used in ref . y. imry , _ introduction to mesoscopic physics _ , oxford university press ( oxford , 1997 ) . d. l. shepelyansky , phys . lett . * 70 * , 1787 ( 1993 ) ; g. kopidakis , s. komineas , s. flach , and s. aubry , phys . lett . * 100 * 084103 ( 2008 ) ; a. s. pikovsky and d. l. shepelyansky , phys . lett . * 100 * , 094101 ( 2008 ) ; s. fishman , a. iomin , and k. mallick phys . e * 78 * , 066605 ( 2008 ) . s. datta , _ electronic transport in mesoscopic systems _ ( cambridge university press , 1995 ) . m. v. berry and s. klein , eur . * 18 * , 222 ( 1997 ) . c. texier , ph.d . thesis , universit paris 6 , 1999 , available online at http://www.lptms.u-psud.fr/membres/texier/research.html
Abstract:

Answer: we consider the phase coherent transport of a quasi one - dimensional beam of bose - einstein condensed particles through a disordered potential of length @xmath0 . among the possible different types of flow we identified [ t. paul , p. schlagheck , p. leboeuf and n. pavloff , phys . rev . lett . * 98 * , 210602 ( 2007 ) ] , we focus here on the supersonic stationary regime where anderson localization exists . we generalize the diffusion formalism of dorokhov - mello - pereyra - kumar to include interaction effects . it is shown that interactions modify the localization length and also introduce a length scale @xmath1 for the disordered region , above which most of the realizations of the random potential lead to time dependent flows . a fokker - planck equation for the probability density of the transmission coefficient that takes this new effect into account is introduced and solved . the theoretical predictions are verified numerically for different types of disordered potentials . experimental scenarios for observing our predictions are discussed .