the lma msw solution  @xcite has been identified  @xcite@xmath6@xcite as the correct solution of the solar neutrino problem .
the @xmath7 conversion probability of this solution gives a very good description of all available data : no statistically significant deviation has been found so far .
new physics effects beyond lma , if exist , are below few per cent .
the program of future solar neutrino studies includes    1 ) .
further tests of the lma solution , in particular , searches for signatures of this solution such as the day - night asymmetry and the distortion (  upturn " ) of the boron @xmath0 spectrum at low energies ;    2 ) .
precise determination of the oscillation parameters , especially the 1 - 2 mixing angle ;    3 ) .
searches for the sub - leading effects which originate from    - 1 - 3 mixing ,    - sterile neutrino mixing ,    - non - standard neutrino interactions ,    - spin - flavor flip in the magnetic fields of the sun ,    - violation of the fundamental symmetries ( cpt , equivalence principle , _
etc_. ) . +
already the present solar neutrino measurements have sensitivity at the level of few per cent .
for instance , the predicted day - night asymmetry of the superkamiokande signal is about @xmath8 which is comparable with the existing @xmath9 experimental error  @xcite . at sno one expects the @xmath10 asymmetry , consistent with the experimental result  @xcite at the @xmath9 level .
future experiments will have substantially higher sensitivity  @xcite .
the solar neutrino studies enter a phase of precision measurements .
+ in this connection it is important    * to give precise description of the lma conversion , both in the sun and in the earth , taking into account various corrections ; * to estimate accuracy of the approximations made ; * to find the precise _ analytic _ expressions for probabilities and observables as functions of the oscillation parameters ( @xmath11 , @xmath12 ) .
this will help to test the lma solution and to search for physics beyond lma .
we address these issues in the present paper . in section [ sec2 ]
we consider the non - adiabatic corrections to the lma conversion probability .
we calculate these corrections for propagation inside the sun and the earth . in section [ sec3 ]
we obtain the analytical formula for the probability averaged over the distribution of neutrino sources . in section [ sec4 ]
we derive the analytic formula for the @xmath0 regeneration effect in the earth .
we present our conclusions in section [ sec5 ] . in the appendices a and b
, alternative derivations of formulas for the regeneration factor are given .
according to the lma msw solution of the solar neutrino problem , a conversion of the solar electron neutrinos is driven by mixing of the two active neutrinos , @xmath13 : _ f = u ( ) _ mass , where , in general , the mixing matrix is determined as u ( ) = , [ matr ] and @xmath14 is the vector of mass states .
the main feature of the lma solution is the adiabaticity of conversion . according to lma the averaged @xmath15 survival probability of the electron neutrinos is given by the adiabatic formula  @xcite : p_ee = ( 1 + 2 ^ 0_m 2 ) . [ adprob ] here @xmath16 is the vacuum mixing angle , @xmath17 is the mixing angle in matter in the neutrino production point , @xmath18 , and the mixing angle in matter is determined by 2_m(v ) = ,   v= g_f n_e(x ) .
[ cost ] here @xmath11 is the mass squared difference , @xmath19 is the neutrino energy , @xmath20 is the potential , @xmath21 is the fermi coupling constant and @xmath22 is the number density of electrons in the point @xmath23 .
how precise the expression ( [ adprob ] ) is , and what are the non - adiabatic corrections ? to answer these questions , we will elaborate on the adiabatic perturbation theory .
dynamics of neutrino conversion is described in terms of the instantaneous eigenstates of the hamiltonian in matter , @xmath24 , @xmath25 . representing an arbitrary neutrino state as @xmath26 ,
we can write the evolution equation in the base @xmath27 , @xmath28 as  @xcite i = , [ evoleigen ] where ( x ) m^2 [ massdiff ] and _ m(x ) = . [ tdot ]
the adiabatic approximation corresponds to a situation when 1 , [ adcond ] and the off - diagonal terms in the hamiltonian ( [ evoleigen ] ) can be neglected . in this case
there are no transitions between the eigenstates , and therefore the eigenstates propagate independently .
the solution of ( [ evoleigen ] ) is straightforward : _
m^ad(x ) = s^ad ( x , x_0 ) _ m ( x_0 ) , where _ m^ad(x ) = ,   _ m ( x_0 ) = ,   [ defpsi ] and the adiabatic evolution matrix is s^ad(x , x_0 ) = s^ad ( ) = .
[ evad ] the adiabatic phase @xmath29 equals ( x ) = _ x_0^x dx ( x ) .    a state initially produced as the electron neutrino , @xmath30 , @xmath31 , evolves as ( x ) = _
m^0 e^i ( x)_1 m + _ m^0 e^-i(x ) _ 2 m . the incoherent survival probability ( [ adprob ] ) can be immediately obtained by averaging the matrix element squared , @xmath32 , over the oscillations .
+      the non - adiabatic corrections correspond to transitions between the instantaneous eigenstates .
we calculate these corrections by solving the equation ( [ evoleigen ] ) .
we will implement a perturbation theory using the fact that for lma the adiabaticity parameter ( [ adcond ] ) is very small .
let us search for the solution of the equation ( [ evoleigen ] ) in the form & & = + & & = , [ evolnew ] where @xmath33 is supposed to hold everywhere along the neutrino trajectory .
( we check this _ a posteriori _ ) .
the expression ( [ evolnew ] ) can be rewritten as _ m ( x ) = s(x , x_0 ) _ m ( x_0 ) , where the evolution matrix equals s(x , x_0 ) c s^ad .
[ ev ]    inserting ( [ evolnew ] ) into ( [ evoleigen ] ) we find the differential equation for @xmath34 from the condition that the off - diagonal elements of the evolution equation for @xmath35 are zero : i c(x)= - c(x)-i _ m(x ) .
[ evolc ] here the first order terms in @xmath34 and @xmath36 are kept only . in this approximation
the energy gap between the states coincides with the adiabatic split @xmath37 given in eq .
( [ massdiff ] ) .    the solution of equation ( [ evolc ] ) can be written in the following form c(x ) = - ^x_x_0 dx .
[ newsolutb ] the integration constant is fixed by the condition : @xmath38 as @xmath39 , so that @xmath40 .    since for the lma solution the phase @xmath29 is a fast oscillating function , the integral in ( [ newsolutb ] ) can be calculated using the following formula ( essentially the integration by parts ) ^b_a f(x ) e^i g(x ) dx = e^i g(x)|^b_a + o(1/g^3 ) [ approxm ] which is valid for smooth functions of @xmath41 and @xmath42 . here
@xmath43 and @xmath44 .
the formula gives very good approximation if @xmath45 . in the case of integral ( [ newsolutb ] ) this condition
coincides with the adiabaticity condition ( [ adcond ] ) which is well satisfied . + according to ( [ approxm ] ) and ( [ newsolutb ] ) we find c(x_f ) & & = -i    & & = -i sign(_m ) |^x = x_f_x = x_0 , [ csol1 ] or explicitly c(x_f)= - i    sign@xmath46
.    we will apply this formula for neutrino propagation inside the sun in section [ sec2.3 ] and inside the earth in section [ sec4 ] .
the survival probability with the adiabaticity violation effect included can be written as p_ee = , [ surprob ] where @xmath47 is the jump probability @xmath6 the probability of transition @xmath48 on the way from @xmath18 to @xmath49 .
let us calculate @xmath50 .
notice that for the lma solution , one can not use the landau - zener formula  @xcite for @xmath50 for the following reasons .
the mixing angle in the final point of evolution , @xmath49 , is large .
the adiabaticity parameter @xmath51 is of the same order for all points inside @xmath52 of the solar radius  @xcite .
the point of maximal adiabaticity violation is not the resonant point , though not far from it .
moreover , the resonance layer defined as @xmath53 , is broad since the vacuum mixing angle is large .
futhermore , significant part of the neutrino flux is produced inside the resonance region or does not cross the resonance region at all .
the double exponential formula  @xcite is not valid too .
it requires production of neutrinos far above the resonance region in the density scale .
this formula is not applied in the range 2~ ( 1.6 - 8.0)10 ^ -6  , [ range ] for which the density in the production point turns out to be close to the resonance density . for the best fit values of the lma oscillation parameters the range ( [ range ] ) corresponds to @xmath54 mev , that is , to the region of interest .
let us apply the results of section 2.2 for calculation of the non - adiabatic corrections .
notice that at the surface of the sun the effective potential @xmath20 is negligible and @xmath55 can be taken zero . using eq .
( [ csol2 ] ) we find the transition amplitude @xmath56 on the way from the production point to the surface of the sun in the leading order approximation as c(x_f)= i |_x = x_0 .
[ nonadinsunb ] then the probability of non - adiabatic transition is given by p_c = |c(x_f)|^2 = = ^2 , [ nonadinsun ] where h(x ) ^-1 ,   l_osc(x ) are the density height and the oscillation length in matter .
the transition probability @xmath50 ( [ nonadinsun ] ) depends only on parameters of the initial ( production ) point .
one can understand this by noting that @xmath57 .
therefore many oscillation lengths are obtained on the distance in which the potential changes sizably .
the non - adiabatic corrections are averaged out being negligible along the trajectory of the neutrino except for the boundaries of trajectory , _ i.e. _ , around the production point or the point at the surface of the sun . at the surface of the sun
the contribution can be neglected because the potential is zero .
the probability is determined basically by square of the ratio of the oscillation length and the density height .
second factor in ( [ nonadinsun ] ) is of the order one .
so , essentially @xmath58 ^ 2 $ ] . using the best fit values of the lma oscillation parameters we find from ( [ nonadinsun ] ) p_c = ( 10 ^ -9 - 10 ^ -7)()^2 .
[ pc ] here the numerical prefactor depends on the production point . as a function of @xmath18 , the probability @xmath50 reaches maximum at around @xmath59 , where the potential does nt drop down substantially and @xmath60 reaches its almost minimal value due to increase of the density gradient .
the corrections are negligible in the whole relevant range of neutrino energies and production points .
the probability ( [ pc ] ) strongly differs from what one would get using double - exponential formula  @xcite : @xmath61 .    notice that the jump probability ( [ nonadinsun ] ) equals ( up to factor 4 ) the adiabaticity parameter in the production point squared , as is expected in the adiabatic perturbation theory
this contrasts the landau - zener probability , @xmath62 , which is essentially non - perturbative effect .
.__the average value of potential @xmath63 and the corresponding value of @xmath64 for different component of the solar neutrino spectrum . _ _ [ cols="^,^,^,^,^,^,^,^,^",options="header " , ]     in the adiabatic approximation the survival probability depends on the potential ( density ) in neutrino production point @xmath65 : p_ee = p_ee(v_0 ) ,   v_0=v(x_0=r ) .
observables at the earth are determined by the survival probability averaged over the neutrino production region : p_k = ,   k = pp , pep , be , n , o , f , b , hep , [ aveprob1 ] where @xmath66 is the distribution of sources of the @xmath67 component of neutrino spectrum .
the distributions are different for different components .
let us introduce the average value of the potential in the production region for the type @xmath67 neutrinos : |v_k .
[ avepot ] we will use the fact that in the effective production region , @xmath68 deviates weakly from @xmath63 .
therefore the survival probability can be expanded in series around @xmath63 : p_ee = p_ee(|v_k ) + ( ) _ v = |v_k(v - |v_k ) + ( ) _ v = |v_k(v - |v_k)^2 + inserting this expression into ( [ aveprob1 ] ) and using the definition ( [ avepot ] ) we find p_k = p_ee(|v_k ) - v^2_k , [ aveprobb ] where v^2_k .
[ delpot2 ] notice that the correction to @xmath69 appears in the second order of the deviation of potential from the average value .
the expression for probability ( [ aveprobb ] ) can be rewritten as p_k = + ( 1-_k ) 2_m(|v_k ) 2 , [ aveprob ] where the correction @xmath70 equals _ k = .
[ deltan ] in the table [ tab : para ] we present the average values of potentials and the corresponding second order deviations from the average values for eight types of solar neutrinos .
we use the distributions of neutrino sources from the bp2000 model  @xcite .
the expansion parameters @xmath71 are all small , especially for the boron neutrinos which have the narrowest distribution of sources .    in fig .
[ figdiff ] we compare the probability @xmath72 obtained from the approximate analytic formula ( [ aveprob ] ) with results of numerical calculations , @xmath73 .
@xmath74 is shown as a function of @xmath75 .
the lines have been cut for @xmath76 because of their lower energies in comparison with the energies of @xmath77 and boron neutrinos .
the plot shows that the analytic formula is rather precise .
in particular , the deviations are extremely small ( @xmath78 ) for small and large values of @xmath79 .
relatively large deviations can be seen in the intermediate region of @xmath79 .
for example , for @xmath80 the magnitude of @xmath81 reaches maximum @xmath82 at around @xmath83  ev@xmath84 .
the corrections @xmath70 are important : _
e.g. _ , for the @xmath77 neutrinos the deviation would be up to @xmath85 without @xmath70 .
the solar neutrinos arrive at the surface of the earth as incoherent fluxes of the mass states .
the mass states oscillate in the matter of the earth producing partial regeneration of the electron neutrino flux  @xcite .
previously the effect has been described in one or two layers approximation . in the later case , interference effects of contributions from the core and the mantle
have been discussed  @xcite . in this section
we will study effects for the realistic density profile of the earth .
the probability of @xmath86 transition can be written as p(_2 _ e ) ^2 + f_reg , [ freg ] where @xmath87 is the regeneration factor which describes the earth matter effect . in the absence of matter ( _ i.e. _ , during the day ) @xmath88 . using the definition ( [ freg ] ) we find the @xmath0 survival probability with the regeneration effect included as p_ee = ( 1 + 2_m^0 2 ) - 2_m^0 f_reg .
notice that the mixing angle @xmath89 in the neutrino production point in the sun determines the mass ( @xmath90 , @xmath91 ) composition of the neutrino flux which arrives at the earth .
the essential feature of the lma solution is that the earth matter effect is small .
this smallness is characterized by the ratio = 0.024 ( ) ( ) ,  v_a= g_f n_a , where @xmath92 is the avogadro number .
we will use @xmath93 as the expansion parameter .
the @xmath86 transition probability can be written as p(_2 _ e ) & & = |_e|u(_mr ) s(x_f , x_0 ) u^(_mr)u()|_2|^2 + & & = |_e|u(_mr ) s(x_f , x_0 ) u^(_mr - ) |_2|^2 , [ p2e ] where @xmath94 is the mixing angle in matter at the surface of the earth and the matrix @xmath95 given in ( [ ev ] ) describes evolution of the neutrino eigenstates in matter . noting that in matter of the earth ( _ mr - ) 21 , we find from ( [ evolnew ] ) , ( [ freg ] ) and ( [ p2e ] ) an expression for the regeneration factor in the lowest order in @xmath34 and @xmath96 as f_reg = ^2 2 ^ 2 ( x_f ) + 2 re\{c(x_f ) } , [ regcor ] where @xmath97 is the potential at the surface of the earth and @xmath98 is the total phase acquired along the trajectory in the earth .    for the profile with slowly changing density ( lowest adiabatic approximation ) , @xmath99
, we obtain f^ad_reg= ^2 2 ^ 2(x_f ) .
[ regold ] the prefactor ( the depth of oscillations ) is determined by the potential at the surface of the earth , whereas the phase is given by the integral along whole trajectory . for one layer with constant potential ( density ) , and
therefore @xmath100 , the regeneration factor ( [ regcor ] ) or ( [ regold ] ) is reduced to the well known expression : f_reg= ^2 2 ^2 .
[ regold1 ] here @xmath101 is the distance traveled by neutrino in the earth and @xmath102 is the oscillation length in matter .
let us consider a neutrino propagation in realistic density profile of the earth .
the profile can be described by @xmath103 nearly spherical shells of matter with sharp ( step - like ) density changes at the borders of shells and slow variation of density in layers between the borders . according to the prem model @xmath104  @xcite .
so , in @xmath105th shell ( @xmath106 ) , the potential @xmath107 is a smooth function of the radial distance @xmath65 . crossing @xmath108 shells
corresponds to crossing @xmath109 layers ( see fig .
[ figearth ] ) .
we denote by @xmath110 the outer radius of @xmath105th shell , so that @xmath111 corresponds to the radius of the earth : @xmath112 .
the trajectory of the neutrino is characterized by the zenith angle @xmath113 .
we determine a position of neutrino along trajectory by the coordinate @xmath23 with origin in the center of trajectory , so that x ,  x^2=r^2-r_e^2 ^ 2_z .
[ trajcoor ] here @xmath101 is the total length of the trajectory in the earth .
the length of the part of trajectory inside border @xmath114 is given by l_i = , [ length ] and @xmath115 by definition .
we introduce the adiabatic phase @xmath116 acquired by neutrinos in the interval @xmath117 along the trajectory , that is , inside the outer border of the @xmath118th shell : _ i & & = ^l_i/2_-l_i/2 dx + & & ^l_i/2_-l_i/2 dx ,  i=0,,n-1 .
[ phasei ] @xmath37 is given in eq .
( [ massdiff ] ) . here
we keep the order @xmath119 term since due to integration its contribution to the phase is not negligible .    at the borders of shells
there are jumps of the potential and hence the discontinuities of the mixing angle in matter .
we denote them as v_i & & v_i+1(r_i)-v_i(r_i ) ,  i=0,,n-1 , [ jump1 ] + _ mi & & _ m(v_i+1(r_i))- _ m(v_i(r_i ) ) ,   i=0 , , n-1 .
[ jump ] at the surface of the earth , we obtain v_0= v_r ,  _ m0 = _ mr- .
[ jump2 ] furthermore , we find _ mi 2 ,  _ mi 1 ,  i=0,,n-1 .
[ approxth ] corrections to ( [ approxth ] ) are of the order @xmath120 , and hence negligible .
smooth variation of the potential in each shell of the earth can be approximated by the analytic formula  @xcite : v_i = v_a ,  i=1,,n . [
poteni ]    let us find an analytic expression for the regeneration factor using the density profile described above .
the problem can be solved in two steps : ( 1 ) computation of the non - adiabatic corrections to propagation within a given layer , @xmath121 ; ( 2 ) computation of effect of the borders between layers , @xmath122 . so that f_reg = f^ad_reg + _ i f_i+_i f_i^jump .
the virtue of the lma solution is that it enables us to study both effects using the same formalism of the adiabatic perturbation theory .
let us compute the non - adiabatic corrections for one layer .
suppose a neutrino trajectory crosses the @xmath105th layer with the borders at @xmath123 and @xmath124 . using ( [ trajcoor ] ) and ( [ poteni ] ) , the potential in this layer can be expressed in terms of the trajectory coordinate as v_i = v_a .
[ potenib ] here _i=_i+ _
i ^2_z+_i ^4_z ,   _i=_i+2 _ i ^2_z ,   _i=_i .
[ parametb ]    according to ( [ length ] ) and ( [ parametb ] ) the gradients of potential at the borders equal |_x = l_i 2 & & = ( _ i+ 2 _ i y_i^2 ) , [ gradi ] +    @xmath125 .
then for this layer eq .
( [ csol2 ] ) gives the amplitude of non - adiabatic transition c_i & & = -2i e^2 i _ 0 e^-2 i ( x ) |^x = l_i-1 2_x = l_i 2 + & & = - v_a e^i _ 0 , [ ampli ] where phases @xmath116 are defined in ( [ phasei ] ) . inserting this expression into ( [ regcor ] ) we obtain the non - adiabatic correction from this layer , @xmath126 to the regeneration factor as f_i & & = v_a .
the ratio of the absolute value of correction and the adiabatic term equals ~ ~ , [ ratio ] where @xmath127 plays the role of typical scale of the density change . as an example , let us consider the layer between @xmath128 and @xmath129 . in this layer @xmath130 , @xmath131 and @xmath132  @xcite . from ( [ ratio ] )
we obtain that the non - adiabatic correction to the regeneration factor is about @xmath133 at @xmath134 mev .
notice that for some particular values of energies and @xmath113 , the contributions @xmath121 from different layers @xmath105 may sum up  constructively " producing larger effect . in this connection
let us notice the following .
\1 ) the enhancement effect may occur for exceptional values of @xmath19 and @xmath113 and therefore any realistic averaging over @xmath19 and integration over @xmath113 will wash it out .
\2 ) the enhancement can not be large ( proportional to the number of layers , @xmath103 ) since ( i ) only few layers give significant contribution and for the rest , the effect is below @xmath135 ; ( ii ) there is a systematic cancellation of contributions from the upper limit of integration in @xmath121 and the lower limit of integration in @xmath136 ( the adiabatic phases are the same for both ) ; ( iii ) typically , contributions from two layers of the same shell have opposite signs .
so , we conclude that the non - adiabaticity within layers of the earth can be safely neglected
.      the jumps of potential between the layers strongly violate the adiabaticity and on the first glance , the adiabatic perturbation theory can not be applied .
we show , however , that the results for non - adiabatic case obtained in section [ sec2.2 ] can be also used here .
the key point is that for the lma parameters the earth matter effects are small , whatever the density profile in the earth is . as a consequence , variations of the mixing angle in matter
are small : @xmath137 , and essentially the expansion parameter here is @xmath93 .    consider a neutrino trajectory which crosses @xmath138 layers ( @xmath103 shells ) . in the points of the trajectory
@xmath139 neutrinos cross the borders of shell with @xmath140 , as illustrated in fig .
[ figearth ] .
the corresponding potential jumps equal @xmath141 for @xmath139 . using the expression ( [ tdot ] )
we obtain in the lowest approximation _
m(x ) = ^n-1_i=1 v_i .
[ theder ] as it has been shown in section [ sec4.2 ] to a good approximation one can take @xmath142 everywhere outside the borders .
the evolution equation ( [ evoleigen ] ) can be averaged in small intervals @xmath143 to eliminate @xmath144-functions which originate from @xmath145 . however , this is not necessary since in the expression for @xmath34 in ( [ evolc ] ) @xmath145 is integrated anyway .    plugging expression ( [ theder ] ) into eq .
( [ newsolutb ] ) we obtain c(l/2 ) = e^i _ 0 ^n-1_i=1 v_i ( e^-i _ i - e^i _ i ) , where @xmath116 are defined in ( [ phasei ] ) .
inserting the real part of @xmath146 into ( [ regcor ] ) we find the regeneration factor for the case of @xmath103 shells crossing : f_reg = ^2 2 ^2 _ 0 + ^n-1_i=1 ^2 2_i _ 0 .
[ regnew ] in the appendix a we present the rigorous derivation of this factor considering evolution in the sequent layers explicitly .
the results of two approaches coincide exactly in the first order in @xmath147 .    noting that @xmath148 is the jump of potential at the surface of the earth
, we can rewrite the expression ( [ regnew ] ) in the following compact form : f_reg = _ 0 ^n-1_i=0 v_i _ i. [ regnewb ] so , @xmath87 is proportional to the sum of similar terms which correspond to the borders of the shells .
each term is the product of the potential jump at a given border and sine of the total adiabatic phase acquired on the part of trajectory inside a given border ( that is , from @xmath149 to @xmath150 for the border @xmath105 ) .
the sum runs over all borders including the surface of the earth . the expression ( [ regnewb ] )
corresponds to the symmetric density profile .
the zenith angle dependence of the regeneration factor appears via the phases : @xmath151 .
the formula ( [ regnewb ] ) ( which is the main result of our study ) allows us to get complete understanding of the earth matter effects including effects of complicated shell structure . apparently , this is not possible using the one layer approximation ( [ regold ] ) , where the interference terms induced by different shells are absent .    in fig .
[ figreg ] we compare the zenith angle dependence of the regeneration factor computed using the analytic formula ( [ regnew ] ) with the one obtained by the exact numerical integration for the prem profile .
two results coincide extremely well .
one can see that the analytic formula reproduces quite precisely the magnitude and the phase structure of the regeneration factor .
let us mark some features .
the change of the oscillatory behaviors for @xmath152 is induced by the sharp density jump at the border between the mantle and the core of the earth , at @xmath153 .
notice that at @xmath152 the amplitude of oscillations for some periods increases , however , the frequency of large peaks becomes lower . so that the average value of @xmath87 does not increase in comparison with the value for @xmath154 .    for small @xmath155
the dependence of @xmath87 on @xmath113 is a result of interference of terms in ( [ regnewb ] ) which correspond to the outer shells of the earth . to understand this dependence , it is convenient to introduce the phase @xmath156 : _ i _ 0 - _ i , [ phased ] so that @xmath157 is the phase acquired by neutrino on the way from the surface of the earth to @xmath140 . using @xmath156 we can rewrite the expression for regeneration factor ( [ regnewb ] ) as f_reg = ^n-1_i=0 v_i .
[ newreg1 ] apparently , if averaging over @xmath156 occurs only the term with @xmath158 survives ( @xmath159 ) in @xmath87 which is reduced to the adiabatic expression for one layer .
the increase of regeneration factor with @xmath155 in the range @xmath160 is related to the effect of three density jumps near the surface of the earth : according to the prem profile they are situated at depths 10 km , 22 km and 31 km correspondingly .
the distance , @xmath161 , on which @xmath156 is acquired depends on the zenith angle as l - l_i 2 ,   _ z > ( @xmath162 is the depth from the surface of the earth to the borders of the @xmath105th shell ) . in the case of small @xmath155 ( @xmath163 ) , the distance @xmath161 can be of several hundreds kilometers which is comparable to or larger than the oscillation length . furthermore ,
@xmath161 and @xmath156 are fast changing functions of @xmath113 .
so , @xmath156 are large and different for different @xmath105 .
therefore , the terms @xmath164 for different @xmath105 ( @xmath165 ) are quite different and therefore partially cancel each other ( `` interfere destructively '' ) .
in constrast , for @xmath166 , the distances @xmath161 ( @xmath165 ) for the outer shells become much smaller than the oscillation length and they slowly change with @xmath155 . in this case the phases @xmath156 ( @xmath165 )
are all small and @xmath167 .
so , for @xmath168 the effects of outer shells `` interfere constructively '' producing larger regeneration factor .
it is then possible to account these close layers effectively as a single layer , as it was done in ref .
increase of the regeneration factor in the transition region @xmath169 corresponds to convergence of the term @xmath170 ( @xmath165 ) to @xmath171 .      for the lma solution
the oscillation length in the earth is small : @xmath172 .
since the time of the neutrino detection is well known , averaging over the zenith angle can be avoided , and in fact , in the unbinned analysis of the data developed recently  @xcite one needs to know the zenith angle dependence without averaging . at the same time since the recoil electron ( and not neutrino ) energy
is measured and a detector has finite energy resolution , averaging over the neutrino energy occurs .    in the leading approximation
the phase equals @xmath173 .
therefore the energy resolution @xmath174 corresponds to averaging over the interval of phases : _ i _ i . if @xmath175 , so that @xmath176 , the interval @xmath177 can be large , thus leading to strong averaging of terms @xmath170 in ( [ regnewb ] ) .
this happens to the contributions from structures situated far from the surface of the earth . in fig .
[ figave ] we show the result of averaging of the regeneration factor folded with the cross section of the neutrino - electron elastic scattering over two different energy intervals . comparing with fig .
[ figreg ] , one sees that the complicated oscillatory pattern produced by the density jumps in the central regions of the earth is strongly averaged when @xmath155 is large ( for general analysis of this effect see  @xcite ) . according to fig .
[ figave ] for @xmath178 the regeneration factor oscillates with small depth around @xmath179 .
this happens because the main term @xmath180 is strongly averaged too .
in contrast , for @xmath163 only the outer structures of the earth can contribute and the averaging is not as efficient as for large @xmath155 .
indeed , for the borders of outer shells ( @xmath165 ) the distance where @xmath156 is acquired , @xmath181 , can be about several hundreds kilometers , that is , comparable with the oscillation length in matter . in this case
the phase interval @xmath182 is still not large enough to give sufficient averaging . for @xmath183 ,
@xmath184 is also small and averaging is weak .
furthermore , in the interval @xmath185 the regeneration factor increases with @xmath155 .
the reason is that in this interval the main term , @xmath171 , starts to `` interfere constructively '' with the terms produced by the outer shells , @xmath170 ( @xmath165 ) , as it has been discussed in section [ sec4.3 ] .    in fig .
[ snocc1 ] and [ snocc2 ] we show dependence of the regeneration factor for the charged current events at sno . here
we have taken into account the energy resolution of the sno detector and also performed integration over various energy bins of the observed kinetic energy .
on the basis of our analytic formulas and discussion , the interpretation of results of fig .
[ snocc1 ] and [ snocc2 ] is straightforward .
there are small scale structures in the outer mantle of the earth of depth @xmath186 km in which matter has quite different densities ( _ e.g. _ , ocean , rock and soil ) .
in contrast to the ideal prem model , these structures are not isotropically distributed and can be quite complicated .    in section [ sec4.3 ]
we have shown for the ideal prem model that contributions produced by structures close the surface of the earth interfere destructively for small @xmath155 . furthermore , averaging over the energy does nt smooth the dependence of these contributions on @xmath155 completely .
this produces an uncertainty for the future high statistics solar neutrino experiments ( see also comments in  @xcite ) unless the local density distribution is well known  @xcite . in view of this
we will consider general ( not spherically symmetric ) density profile .
suppose neutrinos cross @xmath187 layers of matter .
the density jumps occur in the points @xmath188 ( @xmath189 ) ; @xmath190 and @xmath191 are the points where neutrinos enter and leave the matter correspondingly . similarly to ( [ theder ] )
we parameterize @xmath55 as _
m= ^k-1_i=1 v_i  ( x - x_i ) , [ thederasy ] where the jump of potential at the @xmath105th border between layers equals v_i = v(x_i+)-v(x_i- ) ,  i=1 , , k-1 .
[ jumpnewden ] @xmath192 is the infinitesimally small distance . noting that the potential is zero for neutrinos before entering the earth and after leaving the earth , we define also v_0 = v(x_0 ) ,  v_k = -v(x_k ) .
[ jumpnewden1 ] plugging the potential jumps into ( [ newsolutb ] ) gives c(x_k)= - ^k-1_i=1 v_i  e^2 i _
i , where _ i = ^x_k_x_i dx ,  i=0 , , k [ phasex ] is the phase acquired from a given border @xmath105 to the final point of the trajectory ( detector ) .
now it is straightforward to compute the regeneration factor , and in the leading order in @xmath147 we obtain f_reg & & = | _ e| u(_m(x_k ) ) s(x_k , x_0 ) u^(_m(x_0))u ( ) | _ 2 |^2 -^2 + & & = - ^k_i=0 v_i 2_i .
[ regasy ] in the appendix b , a direct computation of @xmath87 for this case is given .
its result coincides with ( [ regasy ] ) in the leading order in @xmath147 . using this formula one can easily reproduce ( [ regnewb ] ) assuming a symmetric density profile and taking into account that @xmath193 for @xmath194 and @xmath195 .    using ( [ regasy ] )
it is easy to study averaging effects following the discussion in section [ sec4.4 ] .
structures situated far from the detector have @xmath196 .
so that , after averaging , remote structures do not produce significant effect .
however , if the energy resolution is improved , effects of these remote structures can be large .
this agrees with general consideration in  @xcite .
thus , we arrive at the following conclusion . if @xmath155 is large , small scale structures near the entering point can be taken effectively as a single layer . furthermore , averaging over the energy makes contributions of these small structures to be unimportant .
if @xmath155 is small , we can not consider small scale structures near the entering point as a single layer .
however uncertainties produced by these structures can be significantly reduced if averaging is performed over broad energy interval , _
i.e. _ @xmath197 . after averaging , the regeneration factor still shows an oscillatory behavior in the region of small @xmath155 , and this effect is produced by contributions of the shells close to the detector  @xcite .
we have performed detailed analytic study of the lma msw conversion of the solar neutrinos .
our main result is the precise analytic formula for the survival probability which includes non - adiabatic corrections , averaging over the neutrino production region and the earth regeneration effect . for the @xmath67 component of the solar neutrino spectrum ( @xmath198 )
it can be written as p_k = + ( 1-_k ) 2_m(|v_k ) 2-(1-_k ) 2_m(|v_k ) f_reg . here
the correction due to averaging effect , @xmath70 , is given in eq .
( [ deltan ] ) ; the average values of matter potential in the production regions of @xmath67 components , @xmath63 , are defined in ( [ avepot ] ) and their numerical values are presented in the table 1 .
the regeneration factor @xmath87 is given in ( [ regnew ] ) for the symmetric density profile and in ( [ regasy ] ) for general asymmetric density profile .
effect of averaging over the neutrino production region in the sun is reduced to specific value of the initial mixing angle in matter which should be taken for the average value of the potential , @xmath199 , and to the appearance of the correction @xmath70 .
we have compared the analytic results with the results of numerical computation and found that maximal deviation @xmath200 happens for the @xmath77 neutrinos .
for the boron neutrinos the precision is better than 0.2% .
we have obtained precise analytic formula for the regeneration effect in the earth using the realistic density profile .
we present simple derivation of this formula which uses the adiabatic perturbation theory . performing also explicit calculations of the evolution in sequent layers
we show that this derivation is correct .
the analytic formula reproduces results of numerical computations with accuracy determined by @xmath201 .
essentially the regeneration effect is the sum of contributions from different shells which are determined by jumps of the potential at the borders and by the adiabatic phase acquired inside the outer borders of the corresponding shells .
the dependence of regeneration factor on the zenith angle can be understood in terms of interference of contributions from different borders .
the derived analytical formula allows us to understand the effect of averaging over the neutrino energy . using the analytical formula we have considered effects of small scale structures ( @xmath202 km ) of the earth profile .
these effects can be important for small values of @xmath203 .
we stress that local `` perturbations '' of the density profile can produce sizable uncertainties in @xmath87 .
+ * acknowledgment *    one of the authors ( p.c.h . ) would like to thank fapesp for financial support .
a.y.s . thanks tokyo metropolitan university where this work has been accomplished for hospitality .
let us derive the expression for the regeneration factor by considering neutrino evolution in sequent layers of the earth explicitly .
first , we find the complete evolution matrix , @xmath204 , in the basis of the mass eigenstates @xmath205 : ( x_f)=(x_f , x_0 ) ( x_0 ) . as discussed in section [ sec4.2 ] , the adiabaticity violation effect within layers is suppressed by @xmath206 in comparison with the leading order earth matter effect ( @xmath207 ) .
therefore , we neglect the adiabaticity violation within layers . 1 ) . in the case of neutrino propagation in one shell ( one layer )
we can simply project the adiabatic evolution matrix ( [ evad ] ) obtained for the matter eigenstates on to the basis of the mass states . in the leading order in @xmath147 we find & & = _ 1 ( , - ) + & & = u^ ( ) u(_mr ) s^ad(,- ) u^(_mr ) u ( ) = u(_m0 ) s^ad(_0 ) u^(_m0 ) + & & = s^ad(,- ) + _ m0 ( e^-i _ 0-e^i _ 0 ) + & & = + ( e^-i _ 0-e^i _ 0 ) .
[ matr1 ] here @xmath208 is the jump of the mixing angle at the surface of the earth . we have used eq .
( [ approxth ] ) , and @xmath209 is defined in eq .
( [ phasei ] ) .    2 ) .
in the case of two shells crossing , the neutrino encounters three layers ( the outer shell is crossed twice ) .
the evolution matrix can be similarly obtained by using the adiabatic evolution matrix ( [ evad ] ) in each layer and by rotation from the matter eigenstates basis in the layer before the border to the basis after the border . as a result ,
we find = u(_m0 ) s^ad ( , ) s_1 ( , - ) s^ad(-,- ) u^(_m0 ) , where @xmath210 is the evolution matrix in the inner shell which has a form similar to eq . ( [ matr1 ] ) and
it can be written as s_1 & & = u(_m1 ) s^ad(,- ) u^(_m1 ) + & & = s^ad(,- ) + 2(e^-i _ 1-e^i _ 1 ) .
[ matr2b ] here @xmath211 is the jump of mixing angle on the border between the first and the second shells , and @xmath212 is defined in eq .
( [ phasei ] ) .    combining the last two formulas we find to the order @xmath147 & & = s_2(,- ) + & & = s^ad(,- ) + _
i=0 ^ 1 ( e^-i _ i - e^i _ i ) .
[ matr2 ]    3 ) .
suppose the evolution matrix for @xmath108 shells crossing ( @xmath109 layers ) equals & & = s_j ( , - ) + & & = s^ad(,- ) + _ i=0^j-1 ( e^-i _ i - e^i _ i ) .
[ matrj ] consider now the trajectory with @xmath213 shells crossings .
the evolution matrix is & & = s_j+1 + & & = ^j-1_i=0 s_j ( , - ) ^0_i = j-1 , [ sss ] where @xmath214 defined in ( [ jump ] ) , is the jump of the mixing angle in matter at the border @xmath114 .
@xmath215 is the evolution matrix in the central shell which can be written similarly to eq .
( [ matr2b ] ) as s_j = s^ad(,- ) + 2(e^-i _ j - e^i _ j ) [ matrjb ] .
@xmath216 is given in eq .
( [ phasei ] ) .
after insertion into ( [ sss ] ) , the first term of ( [ matrjb ] ) leads to @xmath217 .
the second term in ( [ matrjb ] ) is already of the order @xmath147 .
note that @xmath214 is small , as is shown in ( [ approxth ] ) .
so , we can approximate @xmath218 by the unit matrix when plugging the second term in ( [ matrjb ] ) into ( [ sss ] ) . as a result , we find = s_j + ( e^-i _ j - e^i _ j ) .
using then expression ( [ matrj ] ) for @xmath217 , the formula ( [ matrj ] ) is immediately extended to the case of crossing @xmath213 shells , thus accomplishing the proof .
the result for the case of @xmath103 shells crossing is = s^ad(,- ) + _ i=0^n-1 ( e^-i _ i - e^i _ i ) .
[ matrn ]    using ( [ matrn ] ) we obtain the expression for the regeneration factor in the leading order in @xmath147 as f_reg & & =| e^-i _ 0+_i=0^n-1 ( e^-i _ i - e^i _ i ) |^2 -^2 + & & = _ 0 _ i=0^n-1v_i _ i , [ regen ] where @xmath116 is given in ( [ phasei ] ) and @xmath219 is defined in ( [ jump1 ] ) .
this expression coincides with ( [ regnew ] ) or ( [ regnewb ] ) which have been obtained in section [ sec4.3 ] using the adiabatic perturbation theory .
as in the section [ sec4.5 ] , we define _ i _
m(x_i+)-_m(x_i- ) ,  i=0 , , k , where @xmath220 and @xmath221 , and @xmath188 are the points of density jumps . @xmath18 and @xmath222 are the initial and final points of neutrino trajectory in matter . we will use the following expression : _ i 2 ,  _ i 1 , [ approxb ] which is a good approximation in the leading order in @xmath223 .
@xmath219 is given in ( [ jumpnewden ] ) .
neglecting the adiabaticity violation within each layer , we obtain the evolution matrix @xmath224 as = ( ^1_i = k u^(_i ) s^ad(x_i , x_i-1 ) ) u^(_0 ) .
[ matrasy ] @xmath225 is the adiabatic evolution matrix given in ( [ evad ] ) .
we approximate @xmath226 as u^(_i ) = + q(_i ) , where q(_i ) .
straightforward computation gives the following result in the leading order in @xmath147 & & = s^ad(x_k , x_0)+q(_k ) s^ad(x_k , x_0 ) + s^ad(x_k , x_0 ) q(_0 ) + & & + ^k-1_i=1 s^ad(x_k , x_i ) q(_i ) s^ad(x_i , x_0 ) + & & = s^ad(x_k , x_0)+ ^k_i=0 .
@xmath227 is defined in ( [ phasex ] ) .
then , using ( [ approxb ] ) , the regeneration factor can be directly computed in the first order in @xmath147 as f_reg & & = | e^-i _ 0-^k_i=0_i e^i(2 _ i-_0 ) |^2 -^2 + & & = - ^k_i=0 v_i 2 _ i. it coincides with ( [ regasy ] ) . +
* note added * + this note has been added on request of the referee
.    1 ) . after the present paper had appeared in the hep - ph archive [ hep - ph/0404042 ] , the preprint by akhmedov _ et .
, [ hep - ph/0404083 ] , has been published in which the analytic integral formula is given for the the regeneration effect in the earth in the three neutrino framework . in the first version of [ hep - ph/0404083 ] the correct oscillation phase in this integral formula has been introduced on the `` heuristic '' basis : it does not follow from their perturbation theory .
correct integral formula ( with the correct phase ) has been derived for the first time in the paper by ioannisian and smirnov , [ hep - ph/0404060 ] . later in the journal version jhep 0405 ( 2004 ) 057 , akhmedov _ et .
_ , have also presented derivation of correct phase .      in ioannisian and smirnov paper hep - ph/0404060 and akhmedov _
_ paper jhep 0405 ( 2004 ) 057 the integral formula has been obtained using the improved perturbation theory in the small parameter @xmath228 . in the present paper
we use the adiabatic perturbation theory .
it can be shown that in the lowest order in @xmath93 both approaches coincide . indeed , inserting expression for @xmath55 from ( [ tdot ] ) into ( [ newsolutb ] ) and performing integration by parts in eq .
( [ newsolutb ] ) of the present paper one can derive the integral formula .
this formula is valid for arbitrary density profile provided that the condition @xmath229 is satisfied . inserting the potential of prem model ( eqs .
( [ jump1 ] , [ potenib ] ) ) into the integral formula one can reproduce the result ( [ regnewb ] ) . however technically
the use of formula ( [ newsolutb ] ) of the present paper is more convenient for the derivation of ( [ regnewb ] ) since the derivative @xmath230 gives @xmath144- functions at the borders of layers and further integration becomes trivial .    2 ) .
the kamland collaboration has published recently results of improved measurements of oscillations on the basis of 766.3 ton - year exposure [ t. araki _ et .
al . _ , hep - ex/0406035 ] . in assumption of the cpt conservation , the global analysis of the solar neutrino data and kamland gives slightly ( @xmath231 ) higher best fit value @xmath232 ev@xmath233 than it was before .
the increase of @xmath11 leads to the corresponding small decrease of the adiabaticity parameter @xmath234 and the expansion parameter @xmath93 for a given energy of neutrinos .
therefore increase in @xmath11 ( i ) further improves the adiabatic perturbation theory and implies that the non - adiabatic correction for probability in the sun is smaller ; ( ii ) diminishes the earth matter regeneration effect ; ( iii ) makes our analytic formula for @xmath87 preciser .
notice that the analytic study of this paper has a general character and does not rely on particular values of @xmath11 .
we use specific value of @xmath11 for illustration only .
s. p. mikheyev and a. yu .
smirnov , yad .
42 * , 1441 ( 1985 ) [ sov .
. phys . * 42 * , 913 ( 1985 ) ] ; nuovo cim . * c9 * , 17 ( 1986 ) ; s. p. mikheyev and a. yu .
smirnov , zhetf , * 91 * , ( 1986 ) , [ sov .
jetp , * 64 * , 4 ( 1986 ) ] ( reprinted in `` solar neutrinos : the first thirty years '' , eds .
j.n.bahcall _ et .
al . _ ) .
next - generation cherenkov detector hyper - kamiokande , k. nakamura , talk at neutrinos and implications for physics beyond the standard model 2002 , http://insti.physics.sunysb.edu/itp/conf/neutrino/talks/nakamura.pdf .
s. p. mikheyev and a. yu .
smirnov , _
86 massive neutrinos in astrophysics and in particle physics _ , proceedings of the sixth moriond workshop , edited by o. fackler and j. tr@xmath236n thanh v@xmath236n ( editions frontires , gif - sur - yvette , 1986 ) , pp . 355 .
baltz and j. weneser , phys . rev .
* d35 * , 528 ( 1987 ) ; a. dar _ et .
_ , phys . rev . *
d 35 * ( 1987 ) 3607 ; s. p. mikheyev and a. yu .
smirnov , sov .
30 ( 1987 ) 759 - 790 ; l. cherry and k. lande , phys .
rev d * 36 * 3571 ( 1987 ) ; s. hiroi , h. sakuma , t. yanagida , m. yoshimura , phys .
lett . b*198
* 403 , ( 1987 ) and prog . theor .
phys . * 78 * 1428 , ( 1987 ) ; a. j. baltz and j. weneser , phys . rev . * d37 * , 3364 ( 1988 ) .
m. spiro and d. vignaud , phys .
b * 242 * 297 ( 1990 ) .
a. j. baltz and j. weneser , phys . rev .
* d50 * , 5971 ( 1994 ) ; _ ibid _ * d51 * , 3960 ( 1994 ) ; m. maris and s.t .
petcov , phys .
d*62*(2000)093006 ; n. hata and p. langacker , phys .
d * 48*(1993)2937 and phys .
d * 50*(1994)632 ; e. lisi , d. montanino , phys . rev .
d*56*(1997)1792 ; j. m. gelb , wai - kwok kwong , s. p. rosen , phys .
lett . * 78*(1997)2296 .
q. y. liu , m. maris and s. t. petcov phys .
d * 56 * 5991 ( 1997 ) ; m. maris , s.t .
petcov , phys .
d * 56 * 7444 , 1997 . m. narayan , g. rajasekaran , r. sinha ; mod .
a13 : 1915 ( 1998 ) ; a. de gouvea , a. friedland , h. murayama , hep - ph/9910286 ; g.l .
fogli , e. lisi , d. montanino , phys . rev .
d * 61 * 073009 ( 2000 ) ; g.l .
fogli , e. lisi , d. montanino , a. palazzo , hep - ph/0008012 .
j. n. bahcall , p. i. krastev , a. yu .
smirnov , phys .
* d60*(1999 ) 093001;p . c. de holanda , c. pena - garay , m.c .
gonzalez - garcia , j. w. f. valle , phys .
d*60*(1999)093010 ; a. h. guth , l. randall , m. serna , jhep * 9908*(1999)018 .
p. i. krastev and a. yu .
smirnov , phys .
b * 226*(1989)341 ; j.m .
losecco , phys .
d*47*(1993)2032 ; q. y. liu and a. yu .
smirnov , nucl .
b * 524*(1998)505 ; q. y. liu , s. p. mikheyev , and a. yu .
smirnov , phys .
b*440*(1998)319 ; s.t .
petcov , phys .
b * 434*(1998)321 ; kh .
akhmedov , a. dighe , p. lipari , and a. yu .
smirnov , nucl .
b*542*(1999)3 ; e. kh .
akhmedov , hep - ph/9903302 ; e.k .
akhmedov , nucl .
b*538*(1999)25 ; m.v .
chizhov , s.t .
petcov , phys .
rev . lett . * 83*(1999)1096 ; e.k .
akhmedov , a.yu .
smirnov , phys .
* 85*(2000)3978 .
e. lisi and d. montanino , phys .
a. ioannissian , a. smirnov , hep - ph/0404060 .
g. l. fogli , e. lisi , a. marrone , d. montanino , a. palazzo , phys .
a. ioannissian , a. smirnov , in preparation .