during cell life , membranes are submitted to an inhomogeneous and variable environment .
local inhomogeneities can be strongly related to biological processes , which has led to experiments investigating the effect of local modifications on biomimetic membranes  @xcite .
for instance , in the inner membrane of mitochondria , the enzymes that use the local ph difference across the membrane to synthesize adenosine triphosphate , the cell s fuel , are located in membrane invaginations called cristae  @xcite .
experiments on model lipid membranes have shown that a local ph change can induce a local dynamical membrane deformation  @xcite , and in particular the formation of cristae - like invaginations  @xcite : membrane shape is tightly coupled to local ph inhomogeneities .
other concentration inhomogeneities in the environment of a cell have a crucial biological role , for instance in chemotaxis or in paracrine signaling .
it is therefore of great interest to study the response of a biological membrane to a local modification of its environment .    motivated by experiments conducted on biomimetic membranes ,
we have developed a theoretical description of the dynamics of a lipid bilayer membrane submitted to a local concentration increase of a substance that reacts reversibly and instantaneously with the membrane lipid headgroups .
we focus on the regime of small deformations , and we treat linear membrane dynamics in the spirit of ref .  @xcite .
while our first works focused on the simple case of a constant modification of the membrane involving only one wavelength  @xcite , we have recently extended our theoretical description in order to take into account the spatiotemporal profile of the fraction of chemically modified lipids resulting from the local reagent concentration increase  @xcite .
this profile is determined by the diffusion of the reagent in the solution that surrounds the membrane . in ref .
@xcite , we compared the predictions of this theoretical description to experimental measurements of the deformation of the membrane of giant unilamellar vesicles caused by local microinjection of a basic solution , and we obtained good agreement between theory and experiment .    in the present article , we pursue the theoretical investigation of the effect of a local chemical modification on a lipid membrane . in general , the dynamics that results from a local reagent concentration increase is quite complex , as it involves the evolution of the reagent concentration profile simultaneously as the response of the membrane .
this is the case in the experimental data analyzed in ref .
@xcite , which corresponds to microinjection steps lasting a few seconds . here
, we show that the effect of the evolution of the reagent concentration profile on the membrane dynamics becomes negligible some time after the beginning of the reagent concentration increase , after what the dynamics corresponds to the response of the membrane to a chemical modification imposed instantaneously .
we find that studying this regime enables to extract interesting properties of the membrane response .
the article is organized as follows .
first , in sec
.  [ mbr_dyn ] , we review the linear dynamics of a membrane submitted to a local chemical modification .
then , in sec
.  [ res ] , we study separately the dynamics associated with each of the two effects that can arise from a chemical modification , namely a spontaneous curvature change and an equilibrium density change of the external monolayer .
we find that a local asymmetric density perturbation between the two monolayers of the membrane relaxes by spreading diffusively in the whole membrane .
intermonolayer friction plays a crucial part in this behavior .
subsequently , in sec .  [ res_b ]
, we treat the general case where both effects are present , and we show how the ratio of the spontaneous curvature change to the equilibrium density change induced by the local chemical modification can be extracted from the dynamics .
this ratio can not be deduced from the study of global modifications of vesicle equilibrium shapes in light of the area - difference elasticity model  @xcite .
finally , sec .
[ ccl ] is a conclusion .
for the article to be self - contained , the present section reviews the linear dynamics of a membrane submitted to a local chemical modification , starting from first principles .
the main points of this description were presented in ref .
@xcite . in that article , we compared theoretical predictions to experimental measurements of the deformation of a membrane submitted to a local and brief ph increase . here , our aim will be to go further in the analysis of our theoretical description in order to understand the fundamental properties of the response of a membrane to a local chemical modification .
our description of the bilayer membrane is based on a local version of the area - difference elasticity membrane model  @xcite .
we focus on small deformations of an infinite flat membrane , and we denote the upper monolayer by  @xmath0 and the lower one by  @xmath1 .    in the absence of a chemical modification , the local state of each monolayer is described by two variables : the local total curvature @xmath2 defined on the membrane midlayer , which is common to both monolayers , and the local scaled density @xmath3 , defined on the midlayer of the membrane , @xmath4 being a reference density . the sign convention for the curvature
is chosen in such a way that a spherical vesicle has @xmath5 .
the free energy @xmath6 per unit area in monolayer @xmath7 reads  @xcite : @xmath8 where @xmath9 represents the tension of the bilayer and @xmath10 its bending modulus , while @xmath11 is the stretching modulus of a monolayer , and @xmath12 denotes the distance between the neutral surface  @xcite of a monolayer and the midsurface of the bilayer . as we assume that the two monolayers of the membrane are identical before the chemical modification , these constants are the same for both monolayers .
the spontaneous curvatures of the two monolayers have the same absolute value @xmath13 and opposite signs , since their lipids are oriented in opposite directions .
the expression for @xmath6 in eq .
( [ fpm0 ] ) corresponds to a general second - order expansion in the small dimensionless local variables @xmath14 and @xmath15 , around the reference state which corresponds to a flat membrane with uniform density @xmath16 .
it is valid for small deformations around this reference state : @xmath17 and @xmath18 , where @xmath19 is a small dimensionless parameter used for bookkeeping purposes , which characterizes the amplitude of the small deformations of the membrane around the reference state .
mathematically , @xmath19 is considered infinitesimal .
note that in general , both @xmath2 and @xmath14 , which describe local small deformations around the reference state , are functions of time and of position on the membrane .
let us now focus on the way the membrane free energy is affected by the local chemical modification .
we consider that the reagent source , which corresponds to the micropipette tip in an experiment , is localized in the water above the membrane .
besides , membrane permeation and flip - flop are neglected given their long timescales .
hence , the chemical modification only affects the upper monolayer , i.e. , monolayer @xmath0 , and not the lower one .
let us denote by @xmath20 the local mass fraction of the lipids of the upper monolayer that are chemically modified : @xmath20 depends on time and position since it arises from the local chemical modification .
we assume that the reagent concentration is small enough for @xmath20 to remain small at every time and position on the membrane , and we characterize this smallness through @xmath21 . in order to describe the chemically modified membrane
, we have to include the third small variable @xmath20 in our second - order expansion of @xmath22 .
we obtain  @xcite : @xmath23 where the constants @xmath24 , @xmath25 , and @xmath26 describe the response of the membrane to the chemical modification .
these constants depend on the reagent that is injected .
their physical meaning will be explained in the next paragraph . besides
, the non - analytical mixing entropy term @xmath27 has been added to our second - order expansion  @xcite .
note that we assume that the three small dimensionless local variables @xmath20 , @xmath14 and @xmath28 are of the same order .
in fact , in the present work , the deformation of the membrane and the density variation are caused by the local chemical modification , i.e. , they are a response to @xmath20 , which justifies that @xmath28 and @xmath14 are of the same order as @xmath20 .
we refer the reader to ref .
@xcite for more details on the derivation of eqs .
( [ fpm0 ] ) and  ( [ fmod ] ) .
the effect of the chemical modification ( i.e. , of @xmath20 ) on the upper monolayer is twofold .
first , the scaled equilibrium density on the neutral surface of the upper monolayer is changed by the amount @xmath29 to first order .
second , the spontaneous curvature of the upper monolayer is changed by the amount @xmath30 to first order , with @xmath31 .
these results are obtained by minimization of the free energy of a homogeneous monolayer with constant mass ( see  [ ap_sc_dens ] ) .
hence , the constants @xmath24 and @xmath32 describe the linear response of the monolayer equilibrium density and of its spontaneous curvature , respectively , to the chemical modification .
this explains the physical meaning of the constants @xmath24 and @xmath26 in eq .
( [ fmod ] ) .
note that @xmath25 corresponds to the quadratic response of the membrane to the chemical modification , but it will not have any relevant effect in the following .      the elastic force densities in a monolayer described by the free - energy densities in eqs .
( [ fpm0][fmod ] ) have been derived in ref .
@xcite to first order in @xmath19 , using the principle of virtual work . as we focus on small deformations of an infinite flat membrane , it is convenient to describe it in the monge gauge by the height @xmath33 , @xmath34 , of its midlayer with respect to the reference plane @xmath35
. then , @xmath36 to second order .
such a description is adapted to practical cases where the distance between the reagent source and the membrane is much smaller than the vesicle radius .
the force per unit area of the reference plane , which we call `` force density '' , then reads to first order in @xmath19 @xmath37 where @xmath38 is the tangential component of the force density in monolayer
@xmath7 " , while @xmath39 is the total normal force density in the membrane . in these formulas , we have introduced the antisymmetric scaled density @xmath40 , and the constant @xmath41 .
these expressions show that both the equilibrium density change and the spontaneous curvature change ( i.e. , both @xmath24 and @xmath32 ) can yield a normal force density , and thus a deformation of the membrane , while only the equilibrium density change can yield a tangential force density and induce tangential lipid flow .
an illustration of the role of these force densities in the membrane response to a local chemical modification is provided in fig .
[ mem_fig ] .    ) .
( a ) : membrane at equilibrium in the reference state ( flat shape , uniform lipid densities in both monolayers ) .
( b ) : membrane just after the instantaneous local chemical modification : the membrane has not deformed yet at this stage , and the densities are still uniform .
some lipids of monolayer @xmath0 are modified : they are represented in dark blue ( @xmath20 represents their local mass fraction ; for clarity @xmath20 is locally large in the figure , but in our work , @xmath21 ) . here
we assume that the result of the chemical modification is to increase the effective size of their headgroups .
hence , the equilibrium density in monolayer @xmath0 is locally decreased ( @xmath42 ) , while the spontaneous curvature is increased in absolute value ( @xmath43 ) .
this results into a normal force density @xmath44 at the center of the modified zone of the membrane ( see eq .
( [ pn_b ] ) ) , and into a tangential force density in monolayer @xmath0 , oriented toward the exterior of the modified zone ( see eq .
( [ pip_b ] ) ) .
both are indicated on the figure by arrows .
( c ) due to the normal force density , the membrane starts deforming upwards .
the maximum height @xmath45 and the full width at half maximum @xmath46 of the deformation are indicated .
note that tangential flows also appear due to tangential force densities.[mem_fig ] ]    using the elastic force densities in eqs .
( [ pip_b][pn_b ] ) , it is possible to describe the dynamics of the membrane to first order in the spirit of ref .
let us outline the derivation of these equations  @xcite , which is described in more detail in  [ apeq ] .
we use a normal force balance for the bilayer ( see  [ ap_mbdyn ] ) , which involves @xmath47 and the normal viscous stress exerted by the surrounding fluid , which is derived in  [ ap_hyd ] .
besides , as each monolayer is a two - dimensional fluid , we write down generalized stokes equations ( see  [ ap_mbdyn ] ) , which involve @xmath38 , the two - dimensional viscous stress associated with the lipid flow , the tangential viscous stress exerted by the surrounding fluid ( see  [ ap_hyd ] ) , and the intermonolayer friction @xmath48 , where @xmath49 is the velocity in monolayer @xmath7 and @xmath50 is the intermonolayer friction coefficient  @xcite .
this intermonolayer friction is a tangential force that comes into play when one monolayer slides relatively to the other one .
finally , we use mass conservation in each monolayer . these dynamical equations are best expressed using two - dimensional fourier transforms of the various fields involved , denoted with hats : for any field @xmath51 which depends on @xmath52 and on time @xmath53 , @xmath54 is such that @xmath55 combining all the above - mentioned equations yields a system of first - order linear differential equations on the two - dimensional variable @xmath56 : @xmath57 where we have introduced the matrix which describes the dynamical response of the membrane  @xcite : @xmath58 here , we have assumed that @xmath59 and @xmath60 .
this is true for all the wavevectors with significant weight in @xmath61 , if the modified lipid mass fraction @xmath20 has a smooth profile with a characteristic width larger than 1 @xmath62 m .
indeed , @xmath63 for water , and typically @xmath64 and @xmath65  @xcite . besides , the forcing term in eq .
( [ ed ] ) reads  @xcite : @xmath66    eqs .
( [ ed]-[ed_defs ] ) show that the membrane deformation is coupled to the antisymmetric density : the symmetry breaking between the monolayers causes the deformation of the membrane . here
, the symmetry breaking is caused by the chemical modification of certain membrane lipids in the external monolayer , i.e. , to the presence of @xmath20 . and
indeed , eq .  ( [ ed_defs_2 ] ) shows that the forcing term in eq .
( [ ed ] ) is proportional to @xmath67 .      before solving eq .
( [ ed ] ) , we need to determine @xmath68 , which is involved in the forcing term @xmath69 . the profile @xmath70 of the mass fraction of chemically modified lipids in the external monolayer arises from the local reagent concentration increase . we focus on reagents that react reversibly with the membrane lipid headgroups . besides , we assume that the reaction between the lipids and the reagent is diffusion - controlled ( see , e.g. , ref .
in other words , the molecular reaction timescales are very small compared to the diffusion timescales .
for such a reversible diffusion - controlled chemical reaction , @xmath70 is instantaneously determined by the local reagent concentration on the membrane , which results from the diffusion of the reagent in the fluid above the membrane .
note that these hypotheses are verified in the experiments analyzed in ref .
@xcite , where the reagent is sodium hydroxide .
we consider that the reagent source is localized in @xmath71 , which would represent the position of the micropipette tip in an experiment ( see ref .
the cylindrical symmetry of the problem then implies that the fields involved in our description only depend on @xmath72 .
we focus on the regime of small deformations @xmath73 , and we work at first order in @xmath74 . besides , we study the linear regime where @xmath75 is proportional to the reagent concentration on the membrane : denoting by @xmath76 the reagent concentration field , we have @xmath77 . to first order in the membrane deformation , this can be simplified into @xmath78 .
the field @xmath79 is determined by the diffusion of the reagent from the local source in the fluid above the membrane . in microinjection experiments ,
there is also a convective transport of the reagent due to the injection , but the pclet number remains so small that diffusion dominates .
besides , since the membrane is a surface and @xmath80 , we can neglect the number of reagent molecules that react with the membrane when calculating @xmath79 . hence , @xmath79 can be obtained by solving the diffusion equation @xmath81 where @xmath82 is the diffusion coefficient of the reagent in the fluid above the membrane , and the source term reads @xmath83 where @xmath84 is heaviside s function .
this corresponds to a constant injection flow from the source starting at @xmath85 .
in addition , the membrane imposes a neumann boundary condition @xmath86 .
this relation , which corresponds to a vanishing flux across the membrane , can be simplified to first order into @xmath87    the solution to this diffusion problem reads @xmath88 where the causal green function @xmath89 of our diffusion problem can be expressed using the method of images  @xcite : @xmath90 where we have introduced the infinite - volume causal green function of the diffusion equation @xmath91 combining eqs .
( [ cint ] , [ g1 ] , [ g2 ] ) provides an analytical expression for @xmath76 , and for @xmath92 . since @xmath93
, we thus obtain an analytical expression for @xmath61 , which reads @xmath94 where erf denotes the error function .    as @xmath95
, @xmath96 converges towards @xmath97      the membrane deformation @xmath98 is given by the solution to eq .
( [ ed ] ) with @xmath61 expressed in eq .
( [ phiq2 ] ) and with the initial condition @xmath99 , corresponding to a non - perturbed membrane .
we have @xmath100 . \label{solution}\ ] ] + here , we have introduced the eigenvalues @xmath101 and @xmath102 of the matrix @xmath103 ( see eq .
( [ ed_defs ] ) ) , the associated eigenvectors @xmath104 and @xmath105 , and the solutions @xmath106 and @xmath107 of the linear system @xmath108 the integral in eq .
( [ solution ] ) and the inverse fourier transform of @xmath98 can be calculated numerically .
this gives the spatiotemporal evolution of the membrane deformation @xmath109 .      both for physical understanding and for computational efficiency , it is interesting to put our description in dimensionless form . using the buckingham pi theorem and choosing @xmath110 as the distance unit , @xmath111 as the time unit and @xmath112 as the mass unit , the ten parameters in our equations yield seven dimensionless numbers : @xmath113    the parameters @xmath114 and @xmath115 compare the characteristic lengthscales @xmath116 and @xmath12 of the membrane to @xmath110 , while @xmath117 is a dimensionless version of the membrane tension @xmath9 .
besides , the parameter @xmath118 quantifies the importance of the reagent diffusion on the membrane dynamics .
we will discuss briefly the effect of varying @xmath118 in the following .    in the case of a lipid membrane in water ,
the only dimensional parameters that can span various orders of magnitudes are @xmath9 and @xmath110 . for @xmath119 , i.e. , for realistic membrane tensions  @xcite , and for @xmath120 m , we will see in the following that the only relevant dimensionless parameters in the dynamics of the membrane are @xmath62 and @xmath121 .
the first one , @xmath62 , then corresponds to the ratio of the two eigenvalues of @xmath103 for @xmath122 ( see sec .  [ anin ] ) : hence , it is a crucial element of the membrane dynamical response .
the second one , @xmath123 , quantifies the relative weight of the spontaneous curvature change and of the equilibrium density change of the external monolayer due to the chemical modification ( see sec .
[ res_b ] ) .
note indeed that @xmath124 , and that @xmath89 and @xmath125 describe the two effects of the reagent on the membrane , through @xmath32 and @xmath24 .
the effect of varying @xmath62 and @xmath123 will thus be discussed in the following .    for our numerical calculations , we initially take @xmath126 these values correspond to the injection of naoh in water ( @xmath127  @xcite ) from a source at @xmath128 m above a floppy membrane with typical constants @xmath129 ( see ref .
@xcite and references therein ) , @xmath130 , @xmath131 and @xmath132  nm  @xcite .
the values of the intermonolayer friction coefficient @xmath50 and of the water viscosity @xmath133 are those given above .
note that the time unit @xmath111 is then equal to one second .
in this section , we will study the two extreme cases @xmath134 and @xmath135 .
in other words , we will study separately the dynamics associated with an equilibrium density change and with a spontaneous curvature change of the external monolayer .
we will see that these two manifestations of the chemical modification of the membrane due to the reagent concentration increase yield different spatiotemporal evolutions of the membrane deformation . the general case where both effects are present
will then be discussed in sec .
[ res_b ] .    _ a priori _ ,
the dynamics of the membrane deformation is quite complex , as it involves the evolution of the reagent concentration profile , due to diffusion , simultaneously as the response of the membrane . in order to investigate the effect of the reagent diffusion on the dynamics of the deformation , and to see when the effect of the evolution of the reagent concentration profile becomes negligible
, we will compare the two following cases : + ( i ) the realistic case where @xmath96 is given by eq .
( [ phiq2 ] ) , + ( ii ) the theoretical case where the stationary modification @xmath136 in eq .
( [ phis ] ) is imposed instantaneously : @xmath137 .
note that this case corresponds to the limit of very large @xmath118 .
let us first focus on the case where @xmath134 , in which only the equilibrium density of the upper monolayer is affected by the chemical modification .
[ hw_density ] shows the evolution with dimensionless time @xmath138 of the membrane deformation height @xmath139 in front of the source , and of the full width at half - maximum @xmath140 of the deformation ( see fig .
[ mem_fig](c ) for an illustration of the definitions of @xmath45 and @xmath46 ) .
the cases ( i ) and ( ii ) introduced above are presented . on these graphs ,
we also show @xmath141 and the full width at half - maximum @xmath142 of @xmath20 in case ( i ) .
[ hw_density ] shows that the membrane undergoes a transient deformation that relaxes to zero while getting broader and broader .
the local asymmetric density perturbation relaxes by spreading in the whole membrane .
however , this process is slowed down by intermonolayer friction , so the density asymmetry is transiently solved by a deformation of the membrane  @xcite .
this deformation is downwards if @xmath143 , and upwards if @xmath144 .      in the realistic case ( i ) , at short times , the membrane dynamics is governed by the evolution of @xmath20 due to the reagent diffusion .
conversely , at long times , once @xmath20 is close enough to its steady - state profile @xmath136 , the dynamics of the membrane deformation is similar in case ( i ) and in case ( ii ) .
this can be seen in fig .
[ hw_density ] : first , @xmath45(i ) and @xmath46(i ) follow @xmath145 and @xmath142 , and then , they have an evolution very similar to those of @xmath45(ii ) and @xmath46(ii ) .
thus , after some time , the effect of reagent diffusion on the dynamics of the membrane deformation becomes negligible , and the dynamics of the realistic case ( i ) can be well approximated by that of the simpler case ( ii ) , which corresponds to the response of the membrane to an instantaneously imposed modification .
, where only the equilibrium density is changed .
the values taken for the other dimensionless numbers are those in eq .
( [ values ] ) . both the realistic case ( i ) where reagent diffusion is accounted for , and the simpler case ( ii ) where the chemical modification is imposed instantaneously , are considered .
( a ) logarithmic plot of the height of the membrane deformation @xmath45 and of the fraction @xmath145 of modified lipids in front of the reagent source versus dimensionless time @xmath146 .
both in case ( i ) and in case ( ii ) , @xmath45 is plotted in units of the extremal value it attains in case ( ii ) .
( b ) logarithmic plot of the width @xmath46 of the membrane deformation and of the width @xmath142 of the fraction of modified lipids versus @xmath146 . both @xmath46 and @xmath142 are plotted in units of @xmath110 .
it can be seen on graphs ( a ) and ( b ) that cases ( i ) and ( ii ) yield similar dynamics for @xmath147 .
the thin red ( gray ) lines correspond to the analytical laws mentioned in the text .
[ hw_density ] ]    the transition time @xmath148 between the diffusion - dominated regime and the membrane - response  dominated regime is determined by the convergence of @xmath20 to @xmath136 .
as the reagent takes a dimensionless time @xmath149 to diffuse from the source to the membrane , we expect @xmath150 .
we studied the dynamics of the deformation height for @xmath151 $ ] by integrating eq .
( [ ed ] ) numerically , and we found that this law is very well verified . for our standard value of @xmath118 ( see eq .
[ values ] ) , used in fig .
[ hw_density ] , we have @xmath152 .
let us discuss analytically the simple case ( ii ) .
the long - time behaviors obtained in this case are especially interesting , since they also apply to the realistic case ( i ) .
let us focus on @xmath119 , i.e. , on realistic membrane tensions , and let us keep the standard values of the other parameters involved in @xmath103 ( see below eq .
( [ values ] ) ) , as these parameters can not vary significantly for a membrane in water .
the eigenvalues of @xmath103 , which correspond to the two independent relaxation rates of a perturbation of the membrane with wavelength @xmath153 , can then be approximated by @xmath154 for all wavevectors with significant weight in @xmath155 if @xmath120 m . within this approximation , for @xmath122
, we have @xmath156 .
this leads to a simple expression of the solution of the dynamical equation eq .
( [ ed ] ) : @xmath157 with @xmath158 .
this expression shows that the only parameter that is relevant in the dynamics is @xmath62 .
note that , in the realistic case ( i ) , the value of @xmath118 is relevant too .
in the long - time limit , as the deformation spreads , we have @xmath159 for all the wavevectors with significant weight in @xmath160 .
in this case , calculating the inverse fourier transform of eq .
( [ lim1 ] ) for @xmath161 yields @xmath162 and the numerical results for @xmath46 are in excellent agreement with @xmath163 as can be seen in fig .
[ hw_density](b ) .
this law can also be obtained analytically if @xmath136 is replaced by a gaussian .    in the short - time limit ,
( [ lim1 ] ) yields @xmath164 and @xmath165 .
these asymptotic behaviors , both for the short - time limit and for the much more interesting long - time limit , are plotted in red ( gray ) lines on fig .
[ hw_density ] .
the transition between the two asymptotic regimes is determined by the value of @xmath62 , which is the only parameter that controls the dynamics in case ( ii ) .
in particular , it is possible to show that the long - term scaling laws are valid for @xmath166 .
we studied the dynamics of the deformation height for @xmath167 $ ] by integrating eq .
( [ ed ] ) numerically , and we found that this law is very well verified .
the long - term scaling @xmath168 shows that the local antisymmetric density perturbation spreads diffusively .
let us emphasize that this diffusive behavior is not related to the diffusion of the reagent in the solution above the membrane , as in the long - term limit , @xmath20 has reached its steady - state profile @xmath136 . the fact that the long - term scaling @xmath168 holds in case ( ii ) , where the profile @xmath136 is established instantaneously , as well as in case ( i ) ( see fig .
[ hw_density](b ) ) , illustrates that this scaling law is not related to the reagent diffusion .
the long - term diffusive spreading of the antisymmetric density perturbation can be understood as follows . at long times , the difference between the stokes equations for each monolayer eq .
( [ baltg+ ] ) and eq .
( [ baltg- ] ) can be approximated by @xmath169 in real space . to obtain this equation from the difference between eq .
( [ baltg+ ] ) and eq .
( [ baltg- ] ) , we have used @xmath59 and @xmath60 ( see below eq .
( [ ed_defs ] ) ) , and also @xmath170 , which holds for @xmath171 .
the last relation can be shown using eq .
( [ lim1 ] ) in the long - time limit .
using eq .
( [ balapp ] ) , the antisymmetric convective mass current @xmath172 can be expressed to first order as @xmath173 this current has a diffusive form . combining this with the mass conservation relations in eq .
( [ massc ] ) finally yields the diffusion equation @xmath174 hence , the local asymmetry in density between the two monolayers finally relaxes by spreading diffusively , with an effective diffusion coefficient @xmath175 .
this slow relaxation is due to intermonolayer friction .
note that we can now interpret the dimensionless number @xmath118 defined in eq .
( [ nbadim ] ) as ( half ) the ratio of the reagent diffusion coefficient @xmath82 to this effective diffusion coefficient .
the fact that a local density asymmetry in the two monolayers of a membrane relaxes by spreading diffusively is generic . here
, the asymmetry is due to a local chemical modification of one monolayer , but it can also be caused , e.g. , by a sudden flip of some lipids from one monolayer to the other  @xcite .
our description of this diffusive behavior generalizes that of ref .
@xcite , which focused on a perturbation with a spherical cap shape and a uniform density asymmetry .
note that this diffusive behavior was first mentioned in ref .
@xcite , within a different theoretical framework , and restricting to axially - symmetric deformations .
let us now focus on the second extreme case , where @xmath135 , i.e. , where only the spontaneous curvature of the upper monolayer is affected by the chemical modification .
[ hw_curv ] shows that the deformation converges to a deformed profile , in contrast with the previous case : while a local density asymmetry spreads on the whole membrane , a local spontaneous curvature modification leads to a locally curved equilibrium shape .
the deformation is upwards if @xmath176 , and downwards if @xmath177 .
, where only the spontaneous curvature is changed .
the values taken for the other dimensionless numbers are those in eq .
( [ values ] ) . both the realistic case ( i ) where reagent diffusion is accounted for , and the simpler case ( ii ) where the chemical modification is imposed instantaneously , are considered .
( a ) logarithmic plot of the height of the membrane deformation @xmath45 and of the fraction @xmath145 of modified lipids in front of the reagent source versus dimensionless time @xmath146 . both in case ( i ) and in case ( ii ) , @xmath45 is plotted in units of the extremal value it attains in case ( ii ) .
( b ) logarithmic plot of the width @xmath46 of the membrane deformation and of the width @xmath142 of the fraction of modified lipids versus @xmath146 . both @xmath46 and @xmath142
are plotted in units of @xmath110 . it can be seen on graphs ( a ) and ( b ) that cases ( i ) and ( ii ) yield similar dynamics for @xmath147 .
the thin red ( gray ) lines correspond to the analytical laws mentioned in the text .
[ hw_curv ] ]    as in sec .
[ eqdensc ] , in the realistic case ( i ) , the membrane dynamics is governed by the reagent diffusion at short times @xmath178 , while at long times @xmath179 , it is similar in case ( i ) and in case ( ii ) .
the long - time dynamics thus corresponds to the response of the membrane to an instantaneously imposed chemical modification .    in the simple case ( ii ) ,
the approximations on @xmath101 and @xmath102 introduced in sec .  [ anin ] yield @xmath180 where we have used the notation @xmath158 as above . calculating the inverse fourier transform of this function for @xmath161 yields @xmath181 for all @xmath182 .
this analytical law for the deformation height @xmath45 is plotted in red ( gray ) lines in fig .
[ hw_curv](a ) .
let us focus on the long - time limit , which is especially interesting , since the results found in case ( ii ) also apply to the realistic case ( i ) .
( [ ancourb ] ) shows that in this limit , @xmath183 , so that the long - time profile of the membrane deformation is fully determined by that of @xmath20 .
in particular , for @xmath184 , @xmath185 and @xmath186 .
for @xmath187 , we find @xmath165 again .
these asymptotic behaviors regarding the deformation width @xmath46 are plotted in red ( gray ) lines in fig .
[ hw_curv](b ) , both for the short - time limit and for the much more interesting long - time limit .    as in sec .
[ eqdensc ] , in case ( ii ) , @xmath62 is the only relevant parameter in the dynamics of the membrane deformation ( see eq .
( [ ancourb ] ) ) , and it determines the transition time between the two asymptotic regimes . in the realistic case ( i ) , the value of @xmath118 is relevant too .
generically , a chemical modification will affect both the equilibrium density and the spontaneous curvature of the upper monolayer .
the general solution of eq .
( [ ed ] ) is a linear combination of the two solutions obtained for the two extreme cases : @xmath134 ( see sec .
[ eqdensc ] ) and @xmath135 ( see sec .  [ scc ] ) . in this section
, we will only discuss the realistic case ( i ) where the reagent diffusion is taken into account .
as in the two extreme cases , the effect on the membrane dynamics of the evolution of @xmath20 due to reagent diffusion is crucial for @xmath178 , while it becomes negligible for @xmath179 .    for a lipid membrane in water such
that @xmath119 and for @xmath188 m , we have seen that in both extreme cases , the dynamics is influenced only by @xmath62 and @xmath118 ( see secs .
[ eqdensc][scc ] ) .
hence , in the general case , the dynamics is influenced by @xmath62 , @xmath118 , and by the parameter @xmath189 which quantifies the relative importance of the spontaneous curvature change and of the equilibrium density change of the external monolayer due to the chemical modification .
indeed , eq .  ( [ pn_b ] ) shows that the ratio of the normal force density arising from the spontaneous curvature change to the normal force density arising from the equilibrium density change is equal to @xmath190 .
the effect of varying @xmath118 was discussed in secs .
[ eqdensc][scc ] , and in addition , @xmath118 can not vary much for a reagent injected in water above a lipid membrane .
hence , we will focus on the influence of @xmath62 and @xmath123 on the membrane dynamics .
the value of the ratio @xmath123 of the spontaneous curvature change to the equilibrium density change is _ a priori _ unknown , as @xmath24 and @xmath32 , which are involved in @xmath125 and @xmath89 , respectively ( see eq .
[ nbadim ] ) , are unknown .
the actual values of @xmath24 and @xmath32 depend on the reagent as well as on the membrane itself , as these two parameters describe the linear response of the membrane to a reagent ( see eq .
[ fmod ] ) .
let us assume that @xmath191 , i.e. , that the equilibrium density change and the spontaneous curvature change induce deformations in the same direction .
this is true , e.g. , for a chemical modification that affects the lipid headgroups in such a way that it yields an effective change of the preferred area per headgroup  @xcite . a rough microscopic model , where lipids are modeled as cones that favor close - packing , then yields @xmath192 , which corresponds to the case where the two effects yield destabilizing normal force densities of similar magnitudes  @xcite .
thus , we expect @xmath193 .
[ hw_mixte](a ) shows the evolution of @xmath45 in case ( i ) for three different values of @xmath123 .
the deformation height features an extremum @xmath194 , and then a relaxation .
this behavior is due to the change of the equilibrium density ( see sec .
[ eqdensc ] ) .
conversely , the nonzero asymptotic deformation @xmath195 arises from the change of the spontaneous curvature ( see sec .
[ scc ] ) .
the relative importance of @xmath194 to @xmath195 thus depends on @xmath123 .
this means that studying the dynamics of the membrane deformation in response to a local chemical modification provides information on the ratio of the spontaneous curvature change to the equilibrium density change induced by this chemical modification .    in fig .
[ hw_mixte](b)(c ) , the ratio @xmath196 is plotted versus @xmath123 for different values of @xmath62 .
indeed , as mentioned above , the membrane dynamics is determined both by @xmath123 and by @xmath62 ( at constant @xmath118 ) .
the order of magnitude of @xmath194 can be estimated assuming an equilibrium density change and/or a spontaneous curvature change of a few percent : @xmath194 is of order @xmath197 m for flaccid membranes , and smaller than @xmath198 m if @xmath199 .
hence , we choose @xmath200\,\mathrm{n / m}$ ] , which yields @xmath201 $ ] for @xmath128 m : such values are taken in fig .
[ hw_mixte](b)(c ) .     in front of the reagent source versus dimensionless time @xmath146 in the realistic case ( i ) , for different values of the parameter @xmath123 , which quantifies the relative importance of the spontaneous curvature change to the equilibrium density change .
the values taken for the other dimensionless numbers are those in eq .
( [ values ] ) . in each case , @xmath45 is plotted in units of its asymptotic value in the case @xmath202 .
the lines show the asymptotic values @xmath195 of @xmath45 for @xmath184 in each case , and @xmath194 denotes the extremal value of @xmath45 .
( b ) and ( c ) : logarithmic plot of @xmath196 versus @xmath123 in case ( i ) , for different values of the dimensionless parameter @xmath62 , which affects the membrane dynamics as well as @xmath123 . from down to up : @xmath203 .
[ hw_mixte ] ]    our study shows that one can deduce @xmath123 , and thus the relative importance of the spontaneous curvature change to the equilibrium density change due to a chemical modification , from the measurement of @xmath196 .
this is very interesting , given that such information can not be deduced from the study of static and global membrane modifications  @xcite .
indeed , the equilibrium vesicle shapes in the area - difference elasticity model are determined by the combined quantity @xmath204 which involves both the equilibrium density and the spontaneous curvature . in this expression , @xmath205 is the dimensionless preferred area difference between the two monolayers , which is related to the asymmetric equilibrium density change , while @xmath206 denotes the dimensionless spontaneous curvature of the bilayer , which is related to the asymmetric spontaneous curvature change  @xcite .
the parameter @xmath207 in eq .
( [ combqty ] ) is a dimensionless number involving the elastic constants of the membrane .
note that the vesicle shape variations due to global modifications of the vesicle environment are usually interpreted as coming only from a change of the spontaneous curvature , under the assumption that the preferred area per lipid is not modified  @xcite .    in order to determine the ratio @xmath123 of the spontaneous curvature change to the equilibrium density change in a practical case , it is necessary to know the value of @xmath62 ( see fig .  [ hw_mixte](b)(c ) ) .
however , as @xmath62 does not involve any parameter that depends on the reagent ( see eq .
( [ nbadim ] ) ) , it is possible to compare the effects of different reagents on the same membrane , i.e. , their values of @xmath123 , even without knowing the precise value of @xmath62 .
we have analyzed theoretically the spatiotemporal response of a membrane submitted to a local heterogeneity , in the realistic case of a local concentration increase of a substance that reacts reversibly and instantaneously with the membrane lipid headgroups . in general , the dynamics of the membrane deformation is quite complex , as it involves the evolution of the reagent concentration profile due to diffusion in the solution above the membrane , simultaneously as the response of the membrane .
we have shown that , some time after the beginning of the reagent concentration increase , the effect of the evolution of the reagent concentration becomes negligible .
studying this regime enables to extract interesting properties of the membrane response .
we have shown that a local density asymmetry between the two monolayers relaxes by spreading diffusively in the whole membrane .
intermonolayer friction plays a crucial part in this behavior .
in addition , we have shown how the relative importance of the spontaneous curvature change to the equilibrium density change can be extracted from the dynamics of the membrane response to the local chemical modification .
this is a significant result since such information can not be deduced from the study of a static and global modification using the area - difference elasticity model .
our description provides a theoretical framework for experiments involving the microinjection of a reagent close to biomimetic membranes . in ref .
@xcite , we used the theoretical model presented here to analyze experimental results corresponding to brief microinjections of a basic solution in the regime of small deformations , and we obtained good agreement between theory and experiment .
the results of the present theoretical work show that it would be interesting to conduct experiments with a continuous injection phase , since the ratio @xmath123 of the spontaneous curvature change to the equilibrium density change could then be determined for various membrane compositions and reagents .    in biomimetic membranes as well as in cells ,
remarkable phenomena occur in the regime of larger deformations : cristae - like invaginations  @xcite , tubulation  @xcite , pearling  @xcite , budding , exo- or endocytosis  @xcite , etc . to study such phenomena
, it would be useful to pursue our study in the nonlinear regime .
we thank miglena i. angelova for motivating us to study this subject .
we thank miglena i. angelova and nicolas puff for interesting discussions .
the spontaneous curvature and the equilibrium density in monolayer @xmath0 can be obtained by minimizing @xmath208 with respect to @xmath209 and @xmath2 for a homogeneous monolayer with constant mass .
first , using the expression of @xmath22 in eq .
( [ fmod ] ) , the minimization with respect to @xmath209 gives , to first order in @xmath19 : @xmath210 then , the minimization with respect to @xmath2 yields to first order , using eq .
( [ min1 ] ) : @xmath211 with @xmath31 ( as defined in the main text ) .
note that , since we assume that @xmath212 and @xmath18 , we must have @xmath213 and @xmath214 for our description to be valid for the values of @xmath209 and @xmath2 that minimize @xmath208 . this property has been used to simplify the results of the minimization .
the scaled lipid density @xmath215 on the neutral surface of the monolayer is related to @xmath209 through @xmath216 to first order .
this relation arises from the geometry of parallel surfaces  @xcite , given that the membrane midlayer and the monolayer neutral surface are parallel surfaces separated by a distance @xmath12 .
hence , eq .
( [ min1 ] ) can be rewritten as @xmath217 this result is independent of the curvature @xmath2 , contrary to that in eq .
( [ min1 ] ) . indeed , by definition , on the neutral surface
, curvature and density are decoupled  @xcite , while these two variables are coupled on other surfaces .
( [ min1b ] ) shows that , due to the chemical modification , the scaled equilibrium density on the neutral surface of monolayer @xmath0 is changed by the amount @xmath218 to first order . besides , eq .
( [ min2 ] ) indicates that the spontaneous curvature of monolayer @xmath0 is changed by the amount @xmath219 to first order .
one might wonder why the spontaneous curvature @xmath220 found by minimization is not simply @xmath221 for @xmath222 ( see eq .
( [ min2 ] ) ) .
this is due to the fact that we work on the membrane midsurface , which is more convenient to study the membrane dynamics .
if the monolayer free energy had been originally written using variables and constants defined on the neutral surface , then the spontaneous curvature found by minimization would correspond exactly to the constant @xmath223 which plays the part of @xmath13 when everything is defined on the neutral surface .
we wish to determine the velocity field @xmath224 in the fluid above ( @xmath0 ) and below ( @xmath1 ) the membrane .
this flow is caused by the deformation of the membrane and by the lateral flow in the membrane : mathematically , it is determined by the boundary conditions corresponding to the continuity of velocity at the interface between the fluid and the membrane .    given the short lengthscales considered
, the dynamics of @xmath225 can be described using stokes equation . adding the incompressibility condition
, we have : @xmath226 where @xmath227 is the excess pressure field in the fluid : the total pressure is @xmath228 , where @xmath229 is a constant and @xmath227 goes to zero far from the membrane . as in the main text , hats indicate two - dimensional fourier transforms . taking the scalar product of @xmath230 and eq .
( [ stokes ] ) , and using eq .
( [ inc ] ) yields @xmath231 taking the scalar product of @xmath232 and eq .
( [ stokes ] ) , and using eq .
( [ interm ] ) yields @xmath233 solving these equations , with the boundary conditions at infinity : @xmath234 for @xmath235 , and at the water
membrane interface : @xmath236 and @xmath237 , finally yields @xmath238\,e^{\mp q z}\,,\label{3d_1}\\ \bm{\hat w}^\pm&=\left\{\bm{\hat v}^\pm-\left[i\bm{q } \,(\partial_t \hat h ) \pm q\ , \bm{\hat v}^\pm\right ] z\,\right\}e^{\mp q z}\,,\label{3d_2}\\
\hat w_z^\pm&=\left\{\partial_t \hat h\pm\left[q \,(\partial_t \hat h ) \mp i \bm{q}\cdot\bm{\hat v}^\pm\right ] z\,\right\}e^{\mp q z}\,.\label{3d_3}\end{aligned}\ ] ]    for our study of the dynamics of the membrane , what is needed is the stress exerted by the fluid on the membrane , i.e. , since we are working at first order in the membrane deformation , the stress in @xmath35 .
the viscous stress tensor of the fluid is defined by @xmath239 where @xmath240 and @xmath241 .
hence , using eqs .
( [ 3d_1 ] , [ 3d_2 ] , [ 3d_3 ] ) , we obtain @xmath242 where we have introduced the tangential part @xmath243 of the stress tensor of the fluid above and below the membrane .      the first dynamical equation we use is a balance of forces per unit area acting normally to the membrane .
it involves the normal elastic force density in the membrane given by eq .
( [ pn_b ] ) and the normal viscous stresses exerted by the fluid above and below the membrane , which are derived in  [ ap_hyd ] ( see eq .
( [ st1 ] ) ) .
it reads : @xmath244 where @xmath133 denotes the viscosity of the fluid above and below the membrane .
besides , as each monolayer is a two - dimensional fluid , we write down generalized stokes equations describing the tangential force balance in each monolayer .
the first force involved is the density of elastic forces given by eqs .
( [ pip_b ] ) and ( [ pim_b ] ) .
the second one arises from the viscous stress in the two - dimensional flow of lipids .
the third one comes from the viscous stress exerted by the water , which is derived in  [ ap_hyd ] ( see eq .
( [ st2 ] ) ) .
the last force that has to be included is the intermonolayer friction  @xcite .
we thus obtain : @xmath245 where @xmath49 denotes the velocity in monolayer @xmath7 , while @xmath246 is the two - dimensional viscosity of the lipids and @xmath50 is the intermonolayer friction coefficient  @xcite .
25 natexlab#1#1[2]#2 , , , , ( ) .
, , , , , , , , , ( ) . , , , , , , ( ) . , , , , , ( ) . , , , , , ( ) .
, , , , ( ) . , , , , , , , ( ) .
, , , ( ) .
, , , , ( ) .
, , , , ( ) . , , , , ( ) .
, , , ( ) . , , , ( ) .
, , , , ( ) . , , ( ) .
, , , , , .
, , , ( ) .
, , , , , , ( ) . , , , . , , ( ) . , , , , , ( ) . , , , , ( ) . , , , , ( ) . , , , , , , ( ) .


in time - reversal experiments a signal emitted by a localized source is recorded at an array of transducers .
it is then re - emitted into the medium reversed in time , that is , the part of the signal that is recorded first is sent back last . because of the time - reversibility of the wave equation the back - propagated signal refocuses approximately at the location of the original source because the array is limited in size .
a striking experimental observation is that the presence of inhomogeneities in the medium improves the refocusing resolution .
the explanation for this super - resolution is multipathing : waves in complex media that are captured by the recording array have undergone multiple scattering making it effectively larger than its physical size .
another important feature of super - resolution in time reversal is that the refocused signal does not depend on the realization of the random medium .
that is , the refocused signal is deterministic .
super - resolution and self - averaging of refocused signals in complicated media has been observed both in laboratory experiments ( see reviews @xcite and references therein ) and in underwater acoustic wave propagation over long distances ( tens of kilometers ) @xcite .
time - reversal techniques have numerous applications ranging from medicine to communications and , more recently , imaging in random media @xcite .
the first mathematical analysis of time - reversal in random media was given by clouet and fouque @xcite , who analyzed refocusing and self - averaging of time - reversed pulses in a one - dimensional layered random medium .
their result was extended to a three - dimensional layered medium in @xcite .
super - resolution in spatial refocusing and its statistical stability for multi - dimensional waves in random media was analyzed in @xcite , in a remote - sensing regime where the paraxial or parabolic wave equation can be used .
the refocusing of the average signal in a full three - dimensional medium , in the regimes of random geometrical optics and radiative transfer ( transport ) , was studied in @xcite .
we also mention that another source of multipathing is the mixing of waves by the boundaries in an ergodic cavity .
this has been studied experimentally in @xcite and mathematically in @xcite .
the purpose of this paper is to analyze time reversal in the radiative transfer regime using the parabolic wave equation , when the waves interact fully with the random inhomogeneities .
we prove mathematically that the refocused signal is self - averaging , which means that it does not depend on the realization of the random medium .
the mathematical quantitiy that we analyze is the wigner measure of a pair of oscillatory solutions of the random schrdinger equation . in the present setting
, the random potential depends in a markovian way on the variable @xmath1 , the main direction of propagation of the waves .
this allowed us to use in @xcite a martingale method to prove that the average of the wigner distribution converges to a solution of the radiative transfer equation . in this paper
we use additional regularity of the wigner measure , available in time - reversal when there is some blurring at the recording array , to show that the whole wigner distribution , and not only its average , converges weakly , as a schwartz distribution and in probability , to the deterministic solution of the transport equation .
the blurring at the recording array provides a priori bounds for the wigner transform in @xmath0 .
these bounds and the markovianity of the random potential in the direction of propagation make the time - reversal problem more tractable mathematically and allow us to prove in a fairly simple and straightforward manner self - averaging of the time - reversed signal .
we recall that the wigner transform is a convenient tool to analyze high frequency wave propagation in deterministic @xcite and random media @xcite
. introduced by wigner in 1932 @xcite , it has been used extensively in the mathematical literature recently . convergence of the average wigner distribution to the solution of the radiative transfer equation was first proved by h. spohn in @xcite for time - independent potentials on small time intervals .
this result was extended to global in time convergence by l. erds and h .-
yau @xcite .
these proofs involve infinite neumann ( diagrammatic ) expansions for the solution of the schrdinger equation and are quite involved technically .
the corresponding problem with time - dependent potentials is much simpler mathematically .
it was treated by us in @xcite in the markovian case , and by f. poupaud and a. vasseur @xcite in the case of finite - range time correlations . in this paper
we use the fact that the wigner family arising in time - reversal problems is more regular than the usual one because blurring is added at the recording array .
this provides some additional regularity usually obtained by considering mixtures of states as , for instance , in @xcite .
the paper is organized as follows : we describe the scaling and obtain an expression for the back - propagated signal in terms of the wigner transform in section [ sec : signal ] . the main result and assumptions on the random medium
are formulated in section [ sec : mainresult ] .
the proofs are presented in section [ sec : wigner ] .
this work was partially supported by onr grants n00014 - 02 - 1 - 0088 and n00014 - 02 - 1 - 0089 .
g.  bal was supported in part by nsf grant dms-0072008 , g. papanicolaou by grants afosr f49620 - 01 - 1 - 0465 and nsf dms-9971972 , and l. ryzhik in part by nsf grant dms-9971742 and the alfred p. sloan foundation .
the pressure field @xmath2 satisfies the scalar wave equation @xmath3 here @xmath4 is the local wave speed that we will assume to be random , and the laplacian operator includes both direction of propagation , @xmath1 , and the transverse variable @xmath5 . in the physical setting , we have @xmath6 .
we consider dimensions @xmath7 to stress that our analysis is independent of the dimension .
if we assume that at time @xmath8 , the wave field has a `` beam - like '' structure in the @xmath1 direction , and if back - scattering may be neglected , we can replace the wave equation by its parabolic approximation @xcite .
more precisely , the pressure @xmath9 may be approximated as @xmath10 where @xmath11 satisfies the schrdinger equation @xmath12 with @xmath13 the transverse laplacian in the variable @xmath14 .
the refraction index @xmath15 , and @xmath16 in ( [ eq : ansatz ] ) is a reference speed .
the rigorous passage to the parabolic approximation from the wave equation has been analyzed in @xcite in the deterministic case , and @xcite in a one dimensional random medium .
a formal derivation of the paraxial approximation that leads to the radiative transfer regime is given below in section [ sec : scaling ] .
we assume that the original source is located in the plane @xmath17 , and the time - reversal mirror is located in the plane @xmath18 as depicted in figure [ fig1 ] .        during the first stage of a time reversal experiment the signal propagates for a time @xmath19 , or equivalently , over a distance @xmath20 , so that the signal arriving at the time - reversal mirror is given by @xmath21 where the green s function solves @xmath22 then the signal is time reversed . for three - dimensional acoustic pulses
, this means that the pressure field is kept unchanged and that the sign of its time derivative is reversed . in the parabolic approximation , this is equivalent to phase conjugation , where @xmath11 is replaced by its complex conjugate @xmath23 .
we assume that the recording array occupies a compact subset of the plane @xmath18 , and introduce a real - valued aperture function @xmath24 .
it represents the restriction of the signal onto the array , and possible amplification by the array that may vary from one receiver to another . in the absence of amplification it is given by the characteristic function of an array set @xmath25 .
we also allow for some blurring of the recorded signal , modeled by a convolution with kernel @xmath26 .
the signal at @xmath27 after time reversal takes then the form @xmath28 the last step consists in letting the signal propagate back to the origin @xmath17 and time reversing it one more time @xmath29 the last step , phase - conjugating at the source , is not performed in real physical experiments but is convenient if we need to compare the back - propagated signal to the original signal .
it affects neither the degree of refocusing nor the self - averaging effect .
the back - propagated signal in time is given at the plane @xmath17 by @xmath30 we will be interested in the sequel in the refocusing and self - averaging properties of the back - propagated signal @xmath31 for each fixed frequency .
self - averaging of the time signal @xmath32 then follows from ( [ eq : psib - time ] ) when @xmath33 has compact support in @xmath34 .
an interesting conclusion of this paper is that the additional blurring by convolution with the kernel @xmath35 makes the signal @xmath31 self - averaging for every frequency @xmath34 .
this should be contrasted with the situation studied in @xcite , where no blurring was introduced and self - averaging was observed only for the full time signal .
we recall the formal passage from the reduced wave equation to the parabolic equation ( [ eq : sch ] ) described in the appendix to @xcite and explain how the radiative transfer scaling arises in this context .
we start with the reduced wave equation @xmath36 and look for solutions of ( [ red - wave ] ) in the form @xmath37 then we obtain @xmath38 the refraction index @xmath39 is weakly fluctuating so that it has the form @xmath40 where @xmath41 models random fluctuations that will be described in detail in section [ sec : random - refr - ind ] . here , @xmath42 and @xmath43 are the correlation lengths of @xmath41 in the transverse and longitudinal directions , respectively , and the small parameter @xmath44 measures the strength of the fluctuation .
the waves propagate over a distance @xmath45 in the @xmath14-plane and @xmath46 in the @xmath1-direction , and we rescale @xmath14 and @xmath1 accordingly .
we also introduce a carrier wave number @xmath47 and replace @xmath48 , @xmath49 being the non - dimensional wave number ( we drop the prime below ) .
the physical parameters determined by the medium are the length scales @xmath42 , @xmath43 and the non - dimensional parameter @xmath50 .
we now explain the scaling of the parameters @xmath45 , @xmath46 and @xmath47 in the radiative transfer regime .
equation ( [ red - unscaled ] ) in the non - dimensional variables @xmath51 , @xmath52 becomes after we drop the primes @xmath53 we introduce the following small parameters : @xmath54 and rewrite ( [ red - scaled ] ) as @xmath55 let us assume now that @xmath56 and denote @xmath57 .
then ( [ red - param ] ) becomes after multiplication by @xmath58 @xmath59 observe that when @xmath60 and @xmath61 , the first term in ( [ red - param1 ] ) is small and may be neglected in the leading order since @xmath62
. then ( [ red - param1 ] ) becomes @xmath63 which is the parabolic wave equation ( [ eq : sch ] ) in the radiative transfer scaling .
the second condition in ( [ parab - assump ] ) implies that the carrier wave number should be chosen as @xmath64 .
then @xmath65 implies that @xmath66 . therefore the correlation length in the longitudinal direction @xmath1
should be much longer than in the transverse plane @xmath14 .
furthermore , we should have @xmath67 which in turn implies that the fluctuation strength @xmath68 .
given that these constraints are satisfied the last condition in ( [ parab - assump ] ) implies that the spatial scales @xmath45 and @xmath46 should be chosen according to @xmath69 note that the constraint @xmath66 implies that @xmath70 , which is the usual constraint for the validity of the parabolic approximation .
the connection between the relations of the physical parameters and the resulting scaling of the parabolic approximation is further discussed in @xcite . the rigorous passage to the parabolic approximation in one dimension in a random medium in a similar scaling
is discussed in @xcite .    to define what we mean
by the quality of the refocused signal , we assume that the initial pulse is centered around a point @xmath71 with support comparable to the transverse correlation length .
then the initial data for the schrdinger equation ( [ eq : sch - rescaled ] ) becomes in the rescaled variables @xmath72 the transducers should obviously be capable of capturing signals of frequency @xmath73 so that the blurring occurs on the scale of the source .
we therefore replace @xmath26 by @xmath74 . finally , we are interested in the refocusing properties of @xmath75 in the vicinity of @xmath71 and consequently introduce the scaling @xmath76 .
we then recast ( [ eq : psib ] ) as @xmath77 where @xmath78 we observe that @xmath79 , so that @xmath80    let us now introduce the auxiliary function @xmath81 as in @xcite , @xmath82 which satisfies the initial value problem @xmath83 it physically describes the component with wave vector @xmath84 that is sent back to the plane @xmath17 by the array .
we then define the wigner transform as @xmath85 where @xmath86 the wigner transform @xmath87 is more regular than @xmath88 because of additional averaging in the wave vector @xmath84 .
more , precisely , while the family @xmath88 is not uniformly bounded in @xmath0 , the family @xmath87 has a uniform bound in @xmath0 .
this regularizing effect is essentially the same as the one obtained by considering mixtures of states for the schrdinger equation as in @xcite .
we can then recast ( [ eq : psibeps2 ] ) as @xmath89 the above formula shows that the asymptotic behavior of @xmath90 as @xmath91 is characterized by that of the wigner transform @xmath92 .
this section presents our two main results .
the first one describes the self - averaging of the back - propagated signal @xmath93 in ( [ eq : psibeps3 ] ) .
its proof is based on the second theorem , of independent interest , which shows that the wigner transform converges in probability to a unique deterministic limit , solution of a transport equation .
[ thm1 ] let the array function @xmath94 and the radially symmetric filter @xmath95 be in @xmath96 , while @xmath97 for a given @xmath98
. assume also that the refraction index @xmath39 satisfies assumptions in section [ sec : random - refr - ind ] below .
then for each @xmath99 the back - propagated signal @xmath100 given by ( [ eq : psibeps3 ] ) converges in probability and weakly in @xmath101 as @xmath102 to the deterministic function @xmath103 the function @xmath104 satisfies the transport equation @xmath105 with initial data @xmath106 and where the operator @xmath107 is defined by @xmath108 here @xmath109 is the fourier transform of the correlation function of @xmath41 , defined by ( [ eq : tilder ] ) below .
the convergence of the wigner transform is described by the following theorem .
[ thm2 ] under the assumptions of theorem [ thm1 ] the wigner distribution @xmath87 converges in probability and weakly in @xmath110 to the solution @xmath104 of the transport equation ( [ eq : treqformal ] ) .
more precisely , for any test function @xmath111 the process @xmath112 converges to @xmath113 in probability as @xmath102 , uniformly on finite intervals @xmath114 .    here
, @xmath115 is the usual scalar product in @xmath116 .
theorem [ thm1 ] is then a corollary of theorem [ thm2 ] and lemma [ lemma2 ] below .
theorem [ thm1 ] provides a mathematical explanation for the important properties of the time - reversal experiment .
first , it ensures that the back - propagated signal is essentially independent of the realization of the random medium in the high - frequency limit since @xmath117 is deterministic .
second , it provides a quantitative description of the back - propagated field , which may be written as @xmath118 the kernel @xmath119 is the fourier transform in @xmath120 of the solution of ( [ eq : treqformal ] ) .
we deduce from this formula that the back - propagated signal has much better refocusing properties in random media than in homogeneous media .
as is explained in @xcite , the reason is that the limiting wigner measure @xmath121 is more regular in its @xmath122 variable when the right - hand side in ( [ eq : treqformal ] ) is present .
this leads to better spatial decay of @xmath123 and hence @xmath124 is localized tighter in a random medium than in a homogeneous medium .
the improvement of the refocusing in random media has been carefully analyzed theoretically and numerically in @xcite , in the fokker - plank approximation of ( [ eq : treqformal ] ) .
we summarize here the main steps of the derivation of the results stated in theorems [ thm1 ] and [ thm2 ] .
we first obtain in lemma [ lemma1 ] uniform bounds in @xmath0 for the wigner transform @xmath87 independently of the realization of the random medium .
lemma [ lemma2 ] then shows that theorem [ thm1 ] is a consequence of theorem [ thm2 ] .
the main mathematical restriction is to assume that the random potential @xmath125 is markovian in its first variable .
we refer to section [ sec : random - refr - ind ] for the details of its construction .
following our previous work in @xcite , we show the convergence of the expectation @xmath126 for every test function @xmath127 .
its proof is simpler ( and less general ) than in @xcite because of the available a priori @xmath0 bounds on the wigner transform , which satisfies the following cauchy problem @xmath128 with @xmath129.\ ] ] since @xmath34 is a fixed parameter , we set @xmath130 without loss of generality .
the solution of ( [ eq : weps ] ) with initial data in @xmath0 is understood in the sense that for every smooth test function @xmath131 , we have @xmath132 here , we have used that @xmath133 is a self - adjoint operator for @xmath115 .
therefore for a smooth function @xmath134 we obtain @xmath135 , where @xmath136 is the solution of the backward problem @xmath137 with the terminal condition @xmath138 .
this defines the process @xmath139 in @xmath140 and generates a corresponding measure @xmath141 on the space @xmath142;l^2({\mathbb r}^{2d}))$ ] of continuous functions in time with values in @xmath0 .
the measure @xmath141 is actually supported on paths inside a ball @xmath143 with the constant @xmath144 as in lemma [ lemma1 ] . the set @xmath145 is the state space for the random process
we have proved in @xcite and will use in the sequel that the family @xmath141 is tight :    [ lemma - tight ] the family of measures @xmath141 is weakly compact .
the proof of convergence of @xmath87 to its deterministic limit is obtained in two steps .
let us fix a deterministic test function @xmath131 .
we use the markovian property of the random field @xmath125 in @xmath1 to construct a first functional @xmath146;x)\!\to c[0,l]$ ] by @xmath147(z)=\langle w,\lambda\rangle(z)-\int_0^z\langle w,{\dfrac{\partial{\lambda}}{\partial{z}}}+{\mathbf k}\cdot\nabla_{\mathbf x}\lambda+ { \cal l}\lambda\rangle(\zeta)d\zeta\ ] ] and show that it is an approximate @xmath141-martingale .
more precisely , we show that @xmath148(z)|{\cal f}_s\right\ } -g_\lambda[w](s)\right|\le c_{\lambda , l}\sqrt{{\varepsilon}}\ ] ] uniformly for all @xmath149;x)$ ] and @xmath150 .
lemma [ lemma - tight ] implies that there exists a subsequence @xmath151 so that @xmath152 converges weakly to a measure @xmath153 supported on @xmath142;x)$ ] .
weak convergence of @xmath141 and the strong convergence ( [ eq : approxmart ] ) together imply that @xmath154(z)$ ] is a @xmath153-martingale so that @xmath155(z)|{\cal f}_s\right\}-g_\lambda[w](s)=0.\ ] ] taking @xmath156 above we obtain as in @xcite the transport equation ( [ eq : treqformal ] ) for @xmath157 in its weak formulation .
the second step is to show that for every test function @xmath131 the new functional @xmath158(z)=\langle w,\lambda\rangle^2(z)- 2\int_0^z\langle w,\lambda\rangle(\zeta ) \langle w,{\dfrac{\partial{\lambda}}{\partial{z}}}+{\mathbf k}\cdot\nabla_{\mathbf x}\lambda+ { \cal l}\lambda\rangle(\zeta)d\zeta\end{aligned}\ ] ] is also an approximate @xmath141-martingale .
we then obtain that @xmath159 , which implies convergence in probability .
it follows that the limit measure @xmath153 is unique and deterministic , and that the whole sequence @xmath141 converges .    that @xmath160(z)$ ] is an approximate @xmath141-martingale uses very explicitly the uniform a priori @xmath0 bound on the wigner distribution @xmath87 .
when the a priori bound is available only in a much larger space ( the space @xmath161 of @xcite ) , we are not able to prove that the functional @xmath162 is an approximate @xmath141-martingale , and actually suspect that the result is not true .
we expect to observe that some long - range correlations survive for sufficiently singular initial data , which would imply that the limiting measure @xmath153 is no longer deterministic .
we describe here the construction of the random potential @xmath125 .
the refraction index @xmath39 in non - dimensional variables is a random function of the form @xmath163 where the non - dimensional parameter @xmath44 measures the strength of fluctuations .
the mathematical analysis of the wigner distribution in random media can be obtained by the method of diagrammatic expansions of the solution of the schrdinger equation as in @xcite . however , here we assume that the random field @xmath125 is a markov process in @xmath1 , as in @xcite , and we are able to analyze the wigner transform @xmath87 directly .
the markovian hypothesis is crucial to simplify the mathematical analysis because it allows us to treat the process @xmath164 as jointly markov .
in addition , @xmath125 is assumed to be stationary in @xmath14 and @xmath1 , mean zero , and is constructed in the fourier space as follows .
let @xmath165 be the set of measures of bounded total variation with support inside a ball @xmath166 @xmath167 and let @xmath168 be a mean - zero markov process on @xmath169 with generator @xmath170 .
the random potential @xmath125 is given by @xmath171 and is real and uniformly bounded : @xmath172 we assume that the process @xmath125 is stationary in @xmath1 and @xmath14 with correlation function @xmath173 @xmath174 in terms of the process @xmath175 this means that given any two bounded continuous functions @xmath176 and @xmath177 we have @xmath178 here @xmath179 is the usual duality product on @xmath180 , and the power spectrum @xmath181 is the fourier transform of @xmath173 in @xmath14 : @xmath182 we assume that @xmath183 for simplicity and define @xmath109 as @xmath184 which is the space - time fourier transform of @xmath185 .
we assume that the generator @xmath170 is a bounded operator on @xmath186 with a unique invariant measure @xmath187 @xmath188 and that there exists @xmath189 such that if @xmath190 then @xmath191 the simplest example of a generator with gap in the spectrum and invariant measure @xmath192 is a jump process on @xmath165 where @xmath193 given ( [ eq : expdecaynew ] ) , the fredholm alternative holds for the poisson equation @xmath194 provided that @xmath195 satisfies @xmath196 .
it has a unique solution @xmath35 with @xmath197 and @xmath198 .
the solution @xmath35 is given explicitly by @xmath199 and the integral converges absolutely because of ( [ eq : expdecaynew ] ) .
the particular markovian model adopted in this paper is somewhat restrictive .
however , the only information about the process that enters into the main result , the transport equation for the wigner distribution in the limit @xmath102 , is the two - point correlation function of @xmath41 .
therefore , we expect that self - averaging of the wigner transform and convergence to the solution of the transport equation hold for much more general classes of random perturbations of the refraction index .
* proof . * it is shown in @xcite that equation ( [ eq : weps ] ) preserves the @xmath0-norm of @xmath87 : @xmath203 .
the initial conditions for @xmath87 are given by @xmath204 so that @xmath205 and ( [ weps - bd ] ) follows .
the above calculation also shows that the initial wigner distribution @xmath206 converges to @xmath207 notice however that one can not expect strong convergence in @xmath0 of @xmath139 to @xmath208 , the solution of the transport equation ( [ eq : treqformal ] ) , because the latter has an @xmath0-norm that decreases as @xmath1 increases while the @xmath0-norm of @xmath87 is preserved .              to obtain the approximate martingale property ( [ eq : approxmart ] )
, one has to consider the conditional expectation of functionals @xmath215 with respect to the probability measure @xmath216 on the space @xmath142;{\cal v}\times x)$ ] generated by @xmath217 and the cauchy problem ( [ eq : weps ] ) .
the only functions we need to consider are actually of the form @xmath218 with @xmath219;{\cal s}({\mathbb r}^{2d})))$ ] . given a function @xmath215
let us define the conditional expectation @xmath220 the weak form of the infinitesimal generator of the markov process generated by @xmath216 is given by @xmath221\right)\lambda\right\rangle,\ ] ] hence @xmath222\right)\lambda\right\rangle(s)ds\ ] ] is a @xmath216-martingale . the operator @xmath223 is defined by @xmath224\psi({\mathbf x},{{\boldsymbol{\eta}}},{\mathbf k},\hat v)= \frac 1i\int_{\rm^d }
\frac{d\hat v({\mathbf p})}{(2\pi)^d } e^{i{\mathbf p}\cdot { { \boldsymbol{\eta } } } } \left[\psi({\mathbf x},{{\boldsymbol{\eta}}},{\mathbf k}-\frac
{ \mathbf p}2)-\psi({\mathbf x},{{\boldsymbol{\eta}}},{\mathbf k}+\frac { \mathbf p}2)\right].\ ] ] the generator ( [ generator ] ) comes from equation ( [ eq : weps ] ) written in the form @xmath225 w_{\varepsilon}.\ ] ]    given a test function @xmath226;{\cal s})$ ] we construct a function @xmath227 with @xmath228 bounded in @xmath229 uniformly in @xmath230 $ ] . the functions @xmath231 will be chosen so that @xmath232 for all @xmath230 $ ] . here
@xmath233 is defined by ( [ g - eps ] ) with @xmath127 replaced by @xmath234 , and @xmath235 is defined by ( [ g - lambda ] )
. the approximate martingale property ( [ eq : approxmart ] ) follows from this .
the functions @xmath236 and @xmath237 are as follows .
let @xmath238 be the mean - zero solution of the poisson equation @xmath239 it is given explicitly by @xmath240.\ ] ] then we let @xmath241 .
furthermore , the second order corrector is given by @xmath242 where @xmath243 is the mean - zero solution of @xmath244 which exists because @xmath245 , and is given by @xmath246(z,{\mathbf x},{{\boldsymbol{\eta}}}+r{\mathbf k},{\mathbf k},\hat v)\right].\end{aligned}\ ] ] using ( [ lambda1eq ] ) and ( [ lambda2eq ] ) we have @xmath247 + \frac 1{\varepsilon}q\right ) \left(\lambda+\sqrt{{\varepsilon}}\lambda_1^{\varepsilon}+{\varepsilon}\lambda_2^{\varepsilon}\right)\right\rangle \\ & & = \left\langle w ,
\left({\dfrac{\partial{}}{\partial{z}}}+{\mathbf k}\cdot\nabla_{\mathbf x}\right )    \lambda+{\cal l}\lambda\right\rangle + \left\langle w ,
\left({\dfrac{\partial{}}{\partial{z}}}+{\mathbf k}\cdot\nabla_{\mathbf x}\right)\left(\sqrt{{\varepsilon}}\lambda_1^{\varepsilon}+ { \varepsilon}\lambda_2^{\varepsilon}\right ) + \sqrt{{\varepsilon}}{\cal k}[\hat v,\frac { \mathbf x}{\varepsilon}]\lambda_2^{\varepsilon}\right\rangle\\ & & = \left\langle w,\left({\dfrac{\partial{}}{\partial{z}}}+{\mathbf k}\cdot\nabla_{\mathbf x}\right)\lambda+{\cal l}\lambda\right\rangle+\sqrt{{\varepsilon}}\langle w,\zeta_{\varepsilon}^\lambda\rangle\end{aligned}\ ] ] with @xmath248\lambda_2^{\varepsilon}.\ ] ] the terms @xmath249 above are understood as differentiation with respect to the slow variable @xmath14 only , and not with respect to @xmath250 .
it follows that @xmath233 is given by @xmath251 and is a martingale with respect to the measure @xmath216 defined on @xmath252;x\times{\cal v})$ ] , the space of right - continuous paths with left - side limits @xcite .
the estimate ( [ eq : approxmart ] ) follows from the following two lemmas .
[ lemma1new ] let @xmath253;{\cal s}({\mathbb r}^{2d}))$ ] .
then there exists a constant @xmath254 independent of @xmath230 $ ] so that the correctors @xmath255 and @xmath256 satisfy the uniform bounds @xmath257 and @xmath258      indeed , ( [ eq : lambda1bdnew ] ) implies that @xmath263 for all @xmath264 and @xmath261 , while ( [ eq : lambda2bdnew ] ) and lemma [ lemma3 ] imply that for all @xmath230 $ ] @xmath265 for all @xmath261 so that ( [ eq : approxmart ] ) follows .        * proof of lemma [ lemma1new ] .
* we only prove ( [ eq : lambda1bdnew ] ) .
since @xmath266 , there exists a constant @xmath259 so that @xmath267 then we obtain using ( [ defcalv ] ) and ( [ eq : expdecaynew ] ) @xmath268\right| \\ & & \le { c}\int_0^\infty dr e^{-\alpha r}\sup_{\hat v }    \int_{\rm^d } |d\hat v({\mathbf p})| \left[|\lambda(z,{\mathbf x},{\mathbf k}-\frac{{\mathbf p}}{2})|+      we show next that @xmath237 is uniformly bounded .
we have @xmath270\right].\end{aligned}\ ] ] the second term above may be written as @xmath271\\ & & = -\int_{\rm^d } \frac{d\hat v({\mathbf p})}{(2\pi)^{d}}e^{i{\mathbf p}\cdot({\mathbf x}/{\varepsilon}+r{\mathbf k } ) } \int_0^\infty ds e^{sq}\int_{\rm^d } \frac{d\hat v({\mathbf q})}{(2\pi)^d } e^{is({\mathbf k}-{\mathbf p}/2)\cdot { \mathbf q}+i({\mathbf x}/{\varepsilon}+r{\mathbf k})\cdot { \mathbf q}}\\ & & \times \left[\lambda({\mathbf x},{\mathbf k}-\frac { \mathbf p}2-\frac { \mathbf q}2)- \lambda({\mathbf x},{\mathbf
k}-\frac { \mathbf p}2+\frac { \mathbf q}2)\right ] \\ & & + \int_{\rm^d } \frac{d\hat v({\mathbf p})}{(2\pi)^d } e^{i{\mathbf p}\cdot({\mathbf x}/{\varepsilon}+r{\mathbf k } ) } \int_0^\infty ds e^{sq}\int_{\rm^d } \frac{d\hat v({\mathbf q})}{(2\pi)^d } e^{is({\mathbf k}+{\mathbf p}/2)\cdot { \mathbf q}+i({\mathbf x}/{\varepsilon}+r{\mathbf k})\cdot { \mathbf q}}\\ & & \times\left [ \lambda({\mathbf x},{\mathbf k}+\frac { \mathbf p}2-\frac { \mathbf q}2)- \lambda({\mathbf x},{\mathbf k}+\frac
{ \mathbf p}2+\frac { \mathbf q}2)\right].\end{aligned}\ ] ] therefore we obtain @xmath272\\ & & \le c\left[|{\cal l}\lambda({\mathbf x},{\mathbf k})|+ \frac{1}{(1+|{\mathbf x}|^{5d})(1+(|{\mathbf k}|-l)^{5d}\chi_{|{\mathbf k}|\ge 5l}({\mathbf k}))}\right]\end{aligned}\ ] ] and the @xmath0-bound on @xmath237 in ( [ eq : lambda1bdnew ] ) follows because the operator @xmath273 is bounded .
the proof of ( [ eq : lambda2bdnew ] ) is very similar and is omitted .
lemma [ lemma1new ] and lemma [ lemma3 ] together with ( [ geps2 ] ) imply the bound ( [ glambda - glambdaeps ] ) .
the tightness of measures @xmath141 given by lemma [ lemma - tight ] implies then that the expectation @xmath274 converges weakly in @xmath140 to the solution @xmath275 of the transport equation for each @xmath230 $ ] .
we now prove that for any test function @xmath127 the second moment @xmath276 converges to @xmath277 .
this will imply the convergence in probability claimed in theorem [ thm2 ] .
the proof is similar to that for @xmath278 and is based on constructing an appropriate approximate martingale for the functional @xmath279 , where @xmath280 is a test function , and @xmath281 .
we need to consider the action of the infinitesimal generator on functions of @xmath282 and @xmath283 of the form @xmath284 where @xmath285 is a given function .
the infinitesimal generator acts on such functions as @xmath286 where @xmath287\mu+ { \mathbf k}^j\cdot\nabla_{{\mathbf x}^j}\mu,\ ] ] with @xmath288\mu=\frac 1i   \int_{\rm^d } d\hat v({\mathbf p})e^{i({\mathbf p}\cdot { { \boldsymbol{\eta}}}_1 ) } \left[\mu({\mathbf k}_1-\frac { \mathbf p}2,{\mathbf k}_2)-\mu({\mathbf k}_1+\frac { \mathbf p}2,{\mathbf k}_2)\right]\ ] ] and @xmath289\mu=\frac 1i   \int_{\rm^d } d\hat v({\mathbf p})e^{i({\mathbf p}\cdot { { \boldsymbol{\eta}}}_2 ) } \left[\mu({\mathbf k}_1,{\mathbf k}_2-\frac { \mathbf p}2)-\mu({\mathbf k}_1,{\mathbf k}_2+\frac { \mathbf p}2)\right].\ ] ] therefore the functional @xmath290     + { \cal k}_2[\hat v,\dfrac{{\mathbf x}_2}{{\varepsilon } } ] )
\big ) \mu\right\rangle(s)ds   \nonumber\end{aligned}\ ] ] is a @xmath291 martingale .
we let @xmath292 be a test function independent of @xmath283 , where @xmath293 , and @xmath294 .
we define an approximation @xmath295 we will use the notation @xmath296 and @xmath297 .
the functions @xmath298 and @xmath299 are to be determined .
we now use ( [ generator2 ] ) to get @xmath300\mu \right\rangle \nonumber\\ & & +   \left\langle   w\otimes w , \left(q+{\displaystyle\sum}_{j=1}^2   { \mathbf k}^j\cdot\nabla_{{{\boldsymbol{\eta}}}^j } \right ) \mu_2 + { \displaystyle\sum}_{j=1}^2 { \cal k}_j\left[\hat v,{{\boldsymbol{\eta}}}^j\right ] \mu_1 + { \dfrac{\partial{\mu}}{\partial{z}}}+ { \displaystyle\sum}_{j=1}^2{\mathbf k}^j\cdot\nabla_{{\mathbf x}^j } \mu   \right\rangle \nonumber\\ & & + \sqrt{\varepsilon}\left\langle w\otimes w , { \displaystyle\sum}_{j=1}^2 { \cal k}_j\left[\hat v,{{\boldsymbol{\eta}}}^j\right ] \mu_2 + \left({\dfrac{\partial{}}{\partial{z}}}+{\displaystyle\sum}_{j=1}^2{\mathbf k}^j\cdot\nabla_{{\mathbf x}^j}\right ) ( \mu_1+\sqrt{\varepsilon}\mu_2)\right\rangle\nonumber . \end{aligned}\ ] ] the above expression is evaluated at @xmath301 .
the term of order @xmath73 in @xmath302 vanishes since @xmath285 is independent of @xmath41 and the fast variable @xmath303 .
we cancel the term of order @xmath304 in the same way as before by defining @xmath298 as the unique mean - zero ( in the variables @xmath283 and @xmath305 ) solution of @xmath306\mu=0.\ ] ] it is given explicitly by @xmath307\\ & & ~~~~~~~~~~~~~~~~~~+ \frac 1i\int_0^\infty dr e^{rq}\int_{\rm^d } d\hat v({\mathbf p } ) e^{ir({\mathbf k}_2\cdot { \mathbf p})+i({{\boldsymbol{\eta}}}_2\cdot { \mathbf p})}\left[\mu({\mathbf k}_1,{\mathbf k}_2-\frac { \mathbf p}2)- \mu({\mathbf k}_1,{\mathbf k}_2+\frac { \mathbf p}2)\right].\end{aligned}\ ] ] when @xmath285 has the form @xmath308 , then @xmath298 has the form @xmath309 with the corrector @xmath269 given by ( [ lambda1-eq ] ) .
let us also define @xmath299 as the mean zero with respect to @xmath310 solution of @xmath311 \mu_1          = \overline{{\displaystyle\sum}_{j=1}^2 { \cal k}_j\big[\hat v,{{\boldsymbol{\eta}}}^j\big ] \mu_1},\ ] ] where @xmath312 .
the function @xmath299 is given by @xmath313 \mu_1({\mathbf x},{{\boldsymbol{\eta}}}+r{\mathbf k},{\mathbf k})}\\ & & ~~~~~~~~~~~~~~~~~~~~-[{\cal k}_1[\hat v,{{\boldsymbol{\eta}}}_1+r{\mathbf k}_1 ] \mu_1]({\mathbf x},{{\boldsymbol{\eta}}}+r{\mathbf k},{\mathbf k},\hat v)\big]\nonumber\\ & & ~~~~~~~~~~~~~~~~~~~~- \int_0^\infty dr e^{rq } \big[\overline{{\cal k}_2[\hat v,{\mathbf k}_2+r{{\boldsymbol{\eta}}}_2 ] \mu_1({\mathbf x},{{\boldsymbol{\eta}}}+r{\mathbf k},{\mathbf k } ) } \nonumber \\ & & ~~~~~~~~~~~~~~~~~~~~ -[{\cal k}_2[\hat v,{{\boldsymbol{\eta}}}_2+r{\mathbf k}_2]\mu_1 ]    ( { \mathbf x},{{\boldsymbol{\eta}}}+r{\mathbf k},{\mathbf k},\hat v ) \big].\nonumber\end{aligned}\ ] ] unlike the first corrector @xmath298 , the second corrector @xmath299 may not be written as an explicit sum of tensor products even if @xmath285 has the form @xmath308 because @xmath298 depends on @xmath283 .
the @xmath291-martingale @xmath314 is given by @xmath315 where @xmath316 is given by @xmath317\mu_2^{\varepsilon}+ \left({\dfrac{\partial{}}{\partial{z}}}+{\displaystyle\sum}_{j=1}^2{\mathbf k}^j\cdot\nabla_{{\mathbf x}^j}\right ) ( \mu_1^{\varepsilon}+\sqrt{\varepsilon}\mu_2^{\varepsilon})\end{aligned}\ ] ] and the operator @xmath318 is defined by @xmath319 \\   + \big[e^{i{\mathbf p}\cdot\frac{{\mathbf x}_2-{\mathbf x}_1}{{\varepsilon}}+ir{\mathbf k}_2\cdot { \mathbf p } }    ( \mu(z,{\mathbf x}_1,{\mathbf k}_1+\frac{{\mathbf p}}{2},{\mathbf x}_2,{\mathbf k}_2-\frac{{\mathbf p}}{2})-     \mu(z,{\mathbf x}_1,{\mathbf k}_1+\frac{{\mathbf p}}{2},{\mathbf x}_2,{\mathbf k}_2+\frac{{\mathbf p}}{2}))\\         -e^{i{\mathbf p}\cdot \frac{{\mathbf x}_2-{\mathbf x}_1}{{\varepsilon}}+ir{\mathbf k}_2\cdot { \mathbf p } }    ( \mu(z,{\mathbf x}_1,{\mathbf k}_1-\frac{{\mathbf p}}{2},{\mathbf x}_2,{\mathbf k}_2-\frac{{\mathbf p}}{2})-     \mu(z,{\mathbf x}_1,{\mathbf k}_1-\frac{{\mathbf p}}{2},{\mathbf x}_2,{\mathbf k}_2+\frac{{\mathbf p}}{2}))\big]\\   + \big[e^{ir{\mathbf k}_1\cdot { \mathbf p}+
i\frac{{\mathbf x}_1-{\mathbf x}_2}{{\varepsilon}}\cdot { \mathbf p } }    ( \mu(z,{\mathbf x}_1,{\mathbf k}_1-\frac{{\mathbf p}}{2},{\mathbf x}_2,{\mathbf k}_2+\frac{{\mathbf p}}{2})-     \mu(z,{\mathbf x}_1,{\mathbf k}_1-\frac{{\mathbf p}}{2},{\mathbf x}_2,{\mathbf k}_2-\frac{{\mathbf p}}{2}))\\         -e^{ir{\mathbf k}_1\cdot { \mathbf p}+ i\frac{{\mathbf x}_1-{\mathbf x}_2}{{\varepsilon}}\cdot { \mathbf p } }    ( \mu(z,{\mathbf x}_1,{\mathbf k}_1+\frac{{\mathbf p}}{2},{\mathbf x}_2,{\mathbf k}_2+\frac{{\mathbf p}}{2})-     \mu(z,{\mathbf x}_1,{\mathbf k}_1+\frac{{\mathbf p}}{2},{\mathbf x}_2,{\mathbf k}_2-\frac{{\mathbf p}}{2}))\big ]
\\    + \big[e^{ir({\mathbf k}_2+\frac{{\mathbf p}}{2})\cdot { \mathbf p } }    ( \mu(z,{\mathbf x}_1,{\mathbf k}_1,{\mathbf x}_2,{\mathbf k}_2)-\mu(z,{\mathbf x}_1,{\mathbf k}_1,{\mathbf x}_2,{\mathbf k}_2+{\mathbf p}))\\    -e^{ir({\mathbf k}_2-\frac{{\mathbf p}}{2})\cdot { \mathbf p } }    ( \mu(z,{\mathbf x}_1,{\mathbf k}_1,{\mathbf x}_2,{\mathbf k}_2-{\mathbf p})-\mu(z,{\mathbf x}_1,{\mathbf k}_1,{\mathbf x}_2,{\mathbf k}_2))\big ] .
\end{array}\ ] ] we have used in the calculation of @xmath318 that for a sufficiently regular function @xmath35 , we have @xmath320     = { \displaystyle\int}_0^\infty dr { \displaystyle\int}_{\rm^d }         \tilde r(r,{\mathbf p } ) f(r,{\mathbf p},-{\mathbf p } ) d{\mathbf p}.\ ] ] the bound on @xmath316 is similar to that on @xmath321 obtained previously as the correctors @xmath322 satisfy the same kind of estimates as the correctors @xmath323 :        unlike the first moment case , the averaged operator @xmath318 still depends on @xmath200 .
we therefore do not have strong convergence of the @xmath291-martingale @xmath314 to its limit yet .
however , the a priori bound on @xmath87 in @xmath0 allows us to characterize the limit of @xmath314 and show strong convergence .
this is shown as follows .
the first and last terms in ( [ eq : l-2eps ] ) that are independent of @xmath200 give the contribution : @xmath329\\ & &   = { \displaystyle\int}_{\rm^d }     \dfrac{d{\mathbf p}}{(2\pi)^{d } }   \hat r(\frac{p^2-k_1 ^ 2}{2},{\mathbf p}-{\mathbf k}_1 ) ( \mu(z,{\mathbf x}_1,{\mathbf p},{\mathbf x}_2,{\mathbf k}_2)-\mu(z,{\mathbf x}_1,{\mathbf k}_1,{\mathbf x}_2,{\mathbf k}_2))\\ & & + \hat r(\dfrac{p^2-k_2 ^ 2}{2},{\mathbf p}-{\mathbf k}_2 ) ( \mu(z,{\mathbf x}_1,{\mathbf k}_1,{\mathbf x}_2,{\mathbf p})-\mu(z,{\mathbf x}_1,{\mathbf k}_1,{\mathbf x}_2,{\mathbf k}_2 ) ) .
\end{aligned}\ ] ] the two remaining terms give a contribution that tends to @xmath330 as @xmath91 for sufficiently smooth test functions .
they are given by @xmath331
we have @xmath332 by bochner s theorem . since @xmath333 and @xmath127
are real quantities , we can take the real part of the above term and , after the change of variables @xmath334 and @xmath335 , obtain @xmath336 we have ( since @xmath285 is real - valued ) @xmath337 using density arguments we may assume that @xmath285 has the form @xmath338 then we have @xmath339 where @xmath340 .
we introduce @xmath341 and use the cauchy - schwartz inequality in @xmath342 and @xmath343 : @xmath344 we use again the cauchy - schwartz inequality , now in @xmath345 , to get @xmath346 this proves that @xmath347 as @xmath102 .
notice that oscillatory integrals of the form @xmath348 are not small in the bigger space @xmath161 , which is natural in the context of wigner transforms and was used in @xcite . in this bigger space
, we can not control @xmath349 and actually suspect that the limit measure @xmath153 may no longer be deterministic .
we therefore deduce that @xmath350 is an approximate @xmath216 martingale .
the limit of the second moment @xmath351 thus satisfies ( weakly ) the transport equation @xmath352 with initial data @xmath353 .
moreover , the operator @xmath354 acting on a tensor product @xmath355 has the form @xmath356={\cal l}\lambda\otimes\lambda+ \lambda\otimes{\cal l}\lambda.\ ] ] this implies that @xmath357 by uniqueness of the solution to the above transport equation with initial conditions given by @xmath358 .
this proves that the limiting measure @xmath359 is deterministic and unique ( because characterized by the transport equation ) and that the sequence @xmath360 converges in probability to @xmath361 .
w. kuperman , w. hodgkiss , h. song , t. akal , c. ferla and d. jackson , phase conjugation in the ocean : experimental demonstration of an acoustic time reversal mirror , j. acoust .
* 103 * , 1998 , 25 - 40 .