we describe in this paper analysis procedures that can answer two of the most important questions of neutrino research .
how can one determine the total solar neutrino fluxes ( @xmath12 , and @xmath10 ) for comparison with solar model predictions ?
how can one determine the sterile contribution to the total solar neutrino fluxes ?
our answers allow for the possibility of an arbitrary mixture in solar neutrino oscillations of active and sterile neutrinos , but require the correctness of the lma solution of the solar neutrino problems and careful attention to all the sources of error ( theoretical as well as experimental).be and @xmath0b neutrinos are unchanged and the numerical results have not been affected significantly , but the present version is more up - to - date with respect to the input data .
we have also added , inspired by the sage discussion , a detailed analysis of what one can learn about @xmath10 neutrinos before there is a dedicated experiment to measure just the @xmath10 neutrino flux . ]
we focus first on determining total fluxes by comparing charged current ( cc ) observables measured in different experiments ; this method yields results as independent as possible of uncertainties due to the presence of sterile neutrinos .
we then describe how similar techniques can be applied to determine total solar neutrino fluxes using a cc experiment plus a neutrino - electron scattering experiment [ or a neutral current ( nc ) measurement ] , which yields results that depend more on the sterile neutrino mix but which can nevertheless be relatively accurate .
the numerical values we estimate for the expected precision with which different quantities can be measured rely upon simulations of the performance of the relevant experiments .
therefore the accuracies that we quote are illustrative ; the actual accuracies that are obtainable can only be determined once the experimental uncertainties are known .
neutrinos change flavors as they travel to the earth from the center of the sun .
this flavor change was seen directly by the comparison of the sudbury neutrino observatory ( sno ) measurement  @xcite of the charged current reaction for @xmath0b solar neutrinos with the super - kamiokande measurement  @xcite of the neutrino - electron scattering rate ( charged plus neutral current ) .
even more clearly , flavor change has been demonstrated by comparing the neutral current measurement by sno with the sno cc measurement  @xcite .
the conclusion that flavor changes occur among solar neutrinos , if based solely upon the comparison of the sno and super - kamiokande event rates , is valid statistically at about the @xmath13 confidence level  @xcite .
the neutral current measurement of sno increases the significance level for flavor changes among solar neutrinos to the @xmath14 confidence level .
the sno and super - kamiokande results demonstrate simply that new physics is required to resolve the long - standing solar neutrino problem  @xcite , i.e. , to understand the origin of the discrepancy between the predictions of the standard solar model  @xcite and the observed solar neutrino event rates @xcite .
if one includes the results of the chlorine @xcite , kamiokande @xcite , sage @xcite , gallex @xcite , and gno @xcite experiments together with the sno ( cc ) and super - kamiokande results , then the combined measurements require  @xcite new physics at @xmath15 and , if the relative temperature scaling of the @xmath8be and @xmath0b neutrino production reactions is taken into account , at @xmath16 .
helioseismological measurements confirm the predicted sound speeds of the standard solar model to better than @xmath17% and show that stellar physics can not account for the discrepancies between standard predictions and the observed solar neutrino rates  @xcite .      the combination of the charged current ( cc ) and the charged plus neutral - current measurement with super - kamiokande has been used by several groups  @xcite to determine the flux of active @xmath0b neutrinos independent of the solar model .
these model - independent determinations of the active flux exploit the similarity between the response functions in the sno and super - kamiokande detectors  @xcite .
in addition , if one includes all the experimental data in a global oscillation solution in which the @xmath0b flux is a free parameter , one obtains a similar ( but slightly smaller ) allowed range for the @xmath0b neutrino flux  @xcite . the measurement of the nc rate by sno  @xcite provides an independent determination of the active @xmath0b neutrino flux .
all of the analyses yield the same result : if electron neutrinos oscillate into only active neutrinos , then the total @xmath0b neutrino flux is in excellent agreement with the flux predicted by the standard solar model .
this close agreement of the active @xmath0b neutrino flux with the total flux predicted by the standard solar model ( ssm )  @xcite is , if the flux of sterile neutrinos is small , an important confirmation of the quantitative theory of stellar evolution .
we summarize below the current best - estimates and the associated @xmath2 uncertainties for the active @xmath0b neutrino flux , @xmath18 .
* standard solar model ( bp00 ) prediction  @xcite : @xmath19 * comparison of sno and super - kamiokande event rates ( cf . esp . ref .
@xcite ) : @xmath20 * sno neutral current measurement ( assuming undistorted spectrum )  @xcite .
@xmath21 * global neutrino oscillation analysis  ( see table  [ tab : fbranges ] this paper and ref .
@xcite ) : @xmath22    the agreement , summarized in eqs .  ( [ eq : itemfb8ssm])([eq : itemfb8global ] ) , between the ssm calculated flux and the measured active flux suggests that the sterile neutrino contribution to the @xmath0b neutrino flux may be small . in this paper
, we ignore this tempting suggestion and instead concentrate on developing methods to determine experimentally the total @xmath0b and @xmath8be neutrino fluxes emitted by the sun , independent of the active - sterile mixture  ( for a discussion of earlier investigations of sterile neutrinos see refs .
@xcite ) .
if we want to understand the particle physics implications of solar neutrino research , we must determine if sterile neutrinos are present in the solar neutrino flux .
moreover , the original  and still valid  goal of solar neutrino research was  @xcite to compare solar model predicted and experimentally measured ( total ) solar neutrino fluxes .
what is the situation if sterile neutrinos exist ?
the total flux of @xmath0b neutrinos could in this case be much larger than the standard solar model prediction ; a major fraction of the total flux that reaches the earth could arrive in a form that is not detected in solar neutrino experiments .
the existing data disfavor ( at @xmath23 cl ) oscillation into purely sterile neutrinos .
nevertheless , a large sterile component is allowed  @xcite if oscillations occur into a combination of active and sterile neutrino states ( see refs .
@xcite for a description of the formalism adopted here ) .
the flux of sterile neutrinos could in principle be large enough to destroy the apparently excellent agreement between the flux predicted by the ssm and the true flux of @xmath0b neutrinos [ which is assumed to be pure active neutrinos in the comparison shown above in eqs .
( [ eq : itemfb8ssm])-([eq : itemfb8global ] ) ] .
a measurement of the total @xmath0b solar neutrino flux , including the sterile component ( if any ) , will provide information that is important for astrophysics and for particle physics .
the motivation for investigating sterile neutrinos is not dependent upon the lsnd  @xcite results that might suggest the existence of sterile neutrinos .
of course , the lsnd results are not supported by other experimental results ( cf .
@xcite ) nor by theoretical predictions that have been precisely confirmed ( like the helioseismological verifications of the standard solar model ) , as is the case for the inference of flavor change based upon the sno - super - kamiokande comparison .    fortunately ,
the kamland reactor neutrino experiment  @xcite , when combined with the sno measurement of the cc flux , is capable of providing a precise determination of the total , i.e. the active plus the sterile , @xmath0b neutrino flux .
we assume throughout this paper the correctness of the currently favored large mixing angle ( lma ) solution to the solar neutrino problem . if the lma solution is not correct , then all global analyses of the available solar and reactor data indicate that either the mass difference , @xmath24 , or the vacuum mixing angle , @xmath25 , will be too small to produce a measurable effect in the kamland experiment ( see e.g. refs .
@xcite ) . in this case
, kamland will be unable to provide the information required to determine the total @xmath0b neutrino flux .
here is the basic physical idea of the method we propose for measuring the total @xmath0b solar neutrino flux . for the kamland  @xcite experiment
, one will know accurately the flux of anti - neutrinos from the @xmath26 reactors that contribute significantly to the measured anti - neutrino events . from measurements of the total event rate and the energy spectrum induced by the surviving @xmath27
, the kamland experimentalists can determine with precision  @xcite the anti - neutrino propagation parameters , @xmath28 and @xmath29 .
both the kamland measurement and the cc sno measurements are disappearance experiments for neutrinos ( or anti - neutrinos ) of similar energies . for the cc measurement made with sno
, one does not know the total @xmath0b neutrino flux created in the sun .
but , assuming conservation of cpt , one can use the propagation parameters @xmath30 and @xmath31 determined by kamland and the measured ( by sno ) cc rate to solve for the flux that gives the observed result .
summarizing , for the kamland experiment one knows the total flux but not the propagation parameters , which are measured . for the sno cc experiment
, one will know ( from kamland ) the propagation parameters and therefore can measure the total flux .
figure  [ fig : solarvskamlandreal ] shows the results of a refined global solution for the solar neutrino oscillation parameters that was made ( see sec .
[ sec : presentknowledge ] and the appendix ) using all the available solar and reactor data .
the figure displays the allowed solar neutrino oscillation contours at @xmath2 , @xmath32 , and @xmath33 .
the results are obtained by the procedures described most recently in ref .
@xcite , where we have used in the present paper the analysis strategy ( a ) ( of ref .
@xcite ) including the 1496 day super - kamiokande data sample  @xcite as well as the sno cc , nc and day - night observations @xcite .
we have also made some improvements ( see the appendix ) in the treatment of the uncertainties in the neutrino cross sections and in the correlation of errors .
for comparison , we also show in fig .
[ fig : solarvskamlandreal ] the small size of the expected allowed region for kamland if this reactor anti - neutrino experiment observes a signal corresponding to the current best - fit point of the solar neutrino analysis . in calculating the kamland allowed region ,
we have made a conservative estimate ( see sec .  [ subsec : precisionb8 ] ) , following the principles discussed in refs .  @xcite .
figure  [ fig : solarvskamlandreal ] shows clearly that kamland has the potential for making a precise measurement of the solar neutrino oscillation parameters , provided that the lma is the correct oscillation solution .    the only complication involved in determining the total @xmath0b flux from a comparison of the kamland and sno cc measurements results from the fact that for the favored large mixing angle ( lma ) solar neutrino oscillation solution matter effects in the sun and the earth can be significant .
matter effects are unimportant for the kamland reactor experiment .
the role of matter effects in solar neutrino experiments depends somewhat upon the _ a priori _ unknown active - sterile mixture , which introduces a calculable uncertainty in the inferred total @xmath0b neutrino flux .    in summary
, the combination of the sno cc neutrino measurement and the kamland anti - neutrino measurement will determine the total ( active plus sterile ) flux , @xmath34 , of @xmath0b solar neutrinos . by subtracting the previously determined flux ( see above ) , @xmath18 , from @xmath35 , one can determine the flux of sterile solar neutrinos .    using similar reasoning , we shall also show that the combined super - kamiokande and kamland measurements can be analyzed to yield an accurate value for the total @xmath0b flux , although in this case the results are somewhat more sensitive to the active - sterile admixture .
the flux of @xmath8be solar neutrinos , which in the ssm is predicted@xcite to be @xmath36 , can be determined in a model independent way from a global analysis of the solar neutrino data assuming only active neutrino oscillations . for example , the latest analysis by garzelli and giunti  @xcite yields @xmath37 at 99% cl .
we shall also show in this paper that one can extract the value of the @xmath8be neutrino flux from measurements of the gallium solar neutrino experiments , gallex , sage , and gno and the results of the sno and kamland measurements .
the value of the @xmath8be flux that will be derived in this way is relatively insensitive to the assumed neutrino oscillation parameters , although it does depend somewhat on the assumed contributions of the cno , @xmath38 , and @xmath10 neutrino fluxes which we adopt from the standard solar model .
the constraint provided by the chlorine experiment is not very useful for determining the total @xmath8be neutrino flux .
one can obtain an independent measurement of the total @xmath8be solar neutrino flux by comparing data from the kamland experiment with data from the borexino  @xcite solar neutrino experiment . if both the kamland and the borexino detectors work as expected , then this method will be more accurate than the methods involving the gallium and chlorine radiochemical detectors .
the determination of the total solar neutrino fluxes , and even more so the determination of the sterile components of these neutrino fluxes , requires precision in both the experimental measurements and the theoretical calculations and analyses .
we present in the appendix a refined discussion of the theoretical errors , and their correlations , for the absorption cross for the gallium and chlorine solar neutrino experiments .
the outline of this paper is as follows . in sec .
[ sec : presentknowledge ] , we describe the current experimental knowledge of the @xmath0b solar neutrino flux . our results are summarized in table  [ tab : fbranges ] both for the special case of oscillations between purely active neutrinos and for the general case of oscillations between electron neutrinos and an active - sterile neutrino admixture .
we limit our analysis to the allowed lma region of solar neutrino oscillations .
we show in sec .
[ sec : determiningb8total ] how one can use the cc measurements with sno and kamland to determine an accurate total @xmath0b solar neutrino flux including experimental and theoretical uncertainties and the possibility of an appreciable active - sterile admixture .
we switch to the @xmath8be flux in sec .
[ sec : be7 ] and evaluate how well one can determine the total @xmath8be solar neutrino flux by also using the results of the gallium ( gallex , sage , and gno ) solar neutrino experiments or the chlorine experiment . in sec .
[ sec : klandes ] , we investigate how well the total @xmath0b and @xmath8be neutrino fluxes can be determined using the combined measurements of kamland and @xmath39 scattering observed in the super - kamiokande ( sec .
[ subsec : kamlandsuperkb8 ] ) and borexino ( sec .
[ subsec : kamlandborexino ] ) detectors .
we show that even in the presence of active - sterile admixtures the total @xmath8be solar neutrino flux may be measured with relatively high accuracy by comparing results from the kamland and the borexino experiments .
we describe and analyze in sec .
[ sec : determiningpp ] three strategies for determining the total @xmath10 solar neutrino flux in the absence of a dedicated experiment that measures separately the @xmath10 neutrinos .
we summarize and discuss our principal conclusions in sec .
[ sec : discuss ] .
we urge the reader to turn first to sec .
[ sec : discuss ] and read there the summary and discussion of our main results and their implications .
the rest of the paper can then be understood more easily .
we generalize in this section the determination of the @xmath0b neutrino flux to the case in which @xmath40 oscillates into a state that is a linear combination of active ( @xmath41 ) and sterile ( @xmath42 ) neutrino states ,    @xmath43    where @xmath44 is the parameter that describes the active - sterile admixture .
this admixture arises naturally in the framework of 4-@xmath45 mixing  @xcite .
the total @xmath0b neutrino flux can be written @xmath46 where @xmath47 .
clearly , the larger the sterile component , the larger the value of @xmath48 that is inferred from the experimental data .    we have performed a global analysis of the solar neutrino data treating the total @xmath0b neutrino flux as a free parameter .
the details of the analysis procedure are the same as those used in ref .
@xcite except where we explicitly state otherwise .
we concentrate here on the lma region , @xmath49 ,   @xmath50 .    to take account of the possibility of oscillations into sterile neutrinos , we determine the allowed regions in the parameter space defined by @xmath24 , @xmath51 , and a third parameter , @xmath52 , that is defined by eq .
( [ eq : linearcombination ] ) .
it is convenient to introduce the dimensionless parameter , @xmath53 where @xmath54 is the total @xmath0b neutrino flux in units of the predicted standard solar model flux .
.the allowed range of the total @xmath0b neutrino flux .
the table presents the allowed range of @xmath55 [ defined by eq .
( [ eq : defnfb ] ) ] that was found in a global solution of all the currently available solar and reactor neutrino data .
the second column gives the allowed range of @xmath55 for an arbitrary mixture of active and sterile neutrinos and the third column gives the range for the case where only active neutrinos are considered .
the results shown were obtained using eq .
( [ eq : fb8determination3 ] ) .
[ tab : fbranges ] [ cols="^,^,^",options="header " , ]     the largest uncertainties in determining @xmath56 using the chlorine data are caused by the experimental errors in the chlorine absorption rate ( @xmath57 ) and the sno cc absorption rate ( @xmath58 ) .
the uncertainty in the chlorine neutrino absorption cross sections are also significant ( @xmath59 to @xmath60 , depending upon the neutrino oscillation parameters ) . to be explicit , the uncertainties in @xmath56 due to the theoretical uncertainties in calculating the chlorine neutrino absorption cross sections that were used in constructing table  [ tab : fbecl ] were evaluated from the following equation , @xmath61 ^ 2 + \left[\sum_{{\rm j=   ^8b},hep}(f_{\rm be}({\rm c.s.}_j\pm 1 \sigma_{\rm c.s._j } ) - f_{\rm be}(i)\right]^2 \label{appex : chlorineferrors}~.\ ] ] ( see the appendix for more details regarding the treatment of the cross section errors . ) even the cno fluxes ( @xmath62 ) and the sno cc absorption cross section ( @xmath63 ) contribute non - negligible errors .
the uncertainty from sterile neutrinos is small ( @xmath64 ) and does not significantly affect the total uncertainty .
remarkably , the individual uncertainties in @xmath56 which arise from several different sources are larger than the current best - fit values of @xmath56 ( see table  [ tab : fbecl ] ) .    for the chlorine based determination ,
we have computed a representative ( average ) shift in @xmath56 by a monte carlo sampling of the allowed kamland oscillation region shown in fig .
[ fig : solarvskamlandreal ] .
we find @xmath65 this uncertainty is so large as to render not very useful the determination of @xmath56 with the aid of chlorine data .
in this section , we show how data from kamland can be combined with @xmath39 scattering data obtained with the super - kamiokande and borexino experiments in order to determine , respectively , the total @xmath0b ( sec .  [ subsec : kamlandsuperkb8 ] ) and @xmath8be ( sec .  [ subsec : kamlandborexino ] ) solar neutrino fluxes . as discussed before ,
the advantage of using purely cc measurements to determine the fluxes is that the answers depend only mildly on the active - sterile admixture .
on the other hand , @xmath39 scattering measurements have the advantage of smaller uncertainties in the interaction cross sections .
moreover , for the determination of the @xmath8be flux , borexino depends less strongly on other neutrino fluxes predicted by the standard solar model than do the radiochemical experiments , chlorine , gallex , sage , and gno .      in this section ,
we show how data from the super - kamiokande and kamland experiments can be combined to measure the total @xmath0b neutrino flux .
we shall see that the determination using super - kamiokande and kamland yields , on average , a value for the total flux than is comparable in precision to what is expected to be obtained using sno and kamland .
the systematic uncertainties are different in the two experiments , sno and super - kamiokande .
therefore , it will be important to compare the total @xmath0b neutrino flux that is inferred using super - kamiokande and kamland with the value that is obtained using sno and kamland .
if the uncertainties are correctly estimated , then the two methods should agree within their quoted errors .
the advantage of using purely charged current measurements with sno and kamland to determine the total @xmath0b flux is that the answer depends only very slightly upon the unknown active - sterile mixture ( as discussed in sec .
[ sec : determiningb8total ] ) .
the principal advantages of the super - kamiokande experiment in the present context is that the neutrino interaction cross section is accurately known  @xcite and the statistical and systematic errors have already been extensively investigated  @xcite .
however , as we shall see below , when using super - kamiokande there is a significant uncertainty ( @xmath67% and @xmath68% ) in the total @xmath0b flux due to the active - sterile mixture , at least with our present knowledge of @xmath44 .
the average survival probability for @xmath0b solar neutrinos can be written in the form    @xmath69    where @xmath70 is the probability of oscillating into active neutrinos and @xmath71 is the ratio of the the @xmath72 and @xmath73 elastic scattering cross - sections .
the expected sensitivity of @xmath54 to the principal sources of errors can be calculated from the following equation :    @xmath74    we suppose that kamland will observe ( with associated uncertainties that we simulate with a monte carlo ) the rate predicted for the global best - fit point shown in fig .  [ fig : solarvskamlandreal ] . for purely active oscillations
, we find that @xmath75= 1\left[1 \pm 0.054\right ] , ~~\cos^2 \eta = 1.0 , \label{eq : skfboactive}\ ] ] where the first error corresponds to the super  kamiokande experimental uncertainty and the second error is caused by the finite size of the allowed kamland region in oscillation parameter space .
the combined kamland and super - kamiokande measurement will , for purely active neutrinos , yield a determination for the total @xmath0b flux that is more accurate than can be obtained with the sno cc measurement and kamland . within the 1 @xmath76 lma region ,
the average uncertainty in the combined kamland and super - kamiokande measurement for the total active @xmath0b neutrino flux is expected to be @xmath77    for purely active neutrinos , super - kamiokande and kamland may provide us with a more accurate determination of the total @xmath0b neutrino flux than sno cc and kamland [ cf .
( [ eq : averagefbsigma ] ) and eq .
( [ eq : activesigmab8kamlandsk ] ) ] .    allowing for the currently allowed @xmath2 ( @xmath78% ) sterile admixture
, we find for the best fit point @xmath79,~~\cos^2
\eta \geq 0.75 \label{eq : borexinofbosterile}\ ] ]    the last error in eq .
( [ eq : borexinofbosterile ] ) corresponds to the present 1@xmath76 uncertainty from the active - sterile admixture @xmath80 . the average uncertainty in the combined kamland and super - kamiokande measurement for @xmath81 is expected to be @xmath82% , within the 1 @xmath76 lma region and if the sterile admixture is as large as currently allowed .
in this section , we show how data from the kamland and borexino experiments can be combined to measure the total @xmath8be neutrino flux .
the principal advantage of using the borexino experiment for this purpose is that the signal in the borexino experiment  @xcite is predicted to be dominated by @xmath8be neutrinos , whereas @xmath8be solar neutrinos are expected to contribute only a relatively small ( and unlabeled ) fraction to the observed event rate in the gallium and chlorine radiochemical measurements ( see table  [ tab : garates ] ) scattering .
the radiochemical detectors do not have energy resolution , only an energy threshold . ] .
however , unlike the cases involving the gallium and chlorine experiments that were discussed in sec .
[ sec : be7 ] , which include only @xmath40 ( cc ) absorption , the borexino experiment detects both @xmath83 scattering and @xmath84 and @xmath85 scattering .
one must consider in the present case the extent to which the active - sterile neutrino admixture influences the detected event rate for each set of oscillation parameters , @xmath86 , determined by kamland . as we shall see quantitatively in the following discussion , this uncertainty regarding
the sterile admixture does not decrease significantly the overall accuracy of the inferred total @xmath8be neutrino flux .
the average survival probability for active neutrinos can be written conveniently in the form    @xmath87    where @xmath88 is the oscillation probability of neutrino fluxes of source @xmath89 into active neutrinos and @xmath90 is the ratio of the the @xmath91 and @xmath73 elastic scattering cross - sections .
the expression for @xmath56 has the same form for the case in which the kamland and borexino experiments are considered together as for the previously discussed cases involving the gallium experiments [ see eq .
( [ eq : fbe ] ) ] and the chlorine experiment [ see eq .
( [ eq : fbechlorine ] ) ] . using the average survival probability defined in eq .
( [ eq : survivalscattering ] ) , we can write    @xmath92 \label{eq : fborexino}\end{aligned}\ ] ]    in the last term of eq .
( [ eq : fborexino ] ) , we include the contributions from @xmath10 , and cno neutrinos .
the contribution from @xmath0b neutrinos is negligible because the observed @xmath0b neutrino flux is about a factor of @xmath93 smaller than the predicted ssm @xmath8be neutrino flux and because @xmath0b neutrinos primarily produce high energy recoils electrons that the borexino detector can distinguish from the low energy recoil electrons produced by @xmath8be neutrinos .
the expected sensitivity of @xmath56 to different sources of errors is given by @xmath94    in order to estimate the accuracy of this method for determining @xmath56 , we suppose that kamland will observe ( with associated uncertainties that we simulate with a monte carlo ) the rate predicted for the global best - fit point , see ref .
@xcite . ] shown in fig .
[ fig : solarvskamlandreal ] and that borexino also will observe a signal corresponding to this best - fit point ( with associated uncertainties ) . in the absence of any published data , we estimate that - based upon experience in previous solar neutrino experiments  @xcite - that borexino will achieve a systematic uncertainty of between @xmath95% and @xmath4% .
we estimate in this way that the combined experiments will yield a determination of @xmath56 that is , for purely active oscillations , @xmath96= 1\left[1 \pm { 0.057 } ( \pm { 0.103)}\right ] , ~~\cos^2 \eta = 1.0 , \label{eq : borexinofbeactive}\ ] ] where the first error corresponds to the borexino experimental uncertainty , the second to the uncertainty in the reconstructed kamland region , and the last error is due to the theoretical uncertainty in the predicted cno fluxes . within the 1 @xmath76 lma region ,
the average uncertainty in the combined kamland and borexino measurement for the total active @xmath8be neutrino flux is expected to be @xmath97 if we allow for a 25% sterile admixture , we find @xmath98 \label{eq : borexinofbesterile}~.\ ] ] the last error in eq .
( [ eq : borexinofbesterile ] ) corresponds to the present 1@xmath76 uncertainty from the active - sterile admixture @xmath80 .
the numbers in parentheses in eq .
( [ eq : borexinofbeactive ] ) and eq .
( [ eq : borexinofbesterile ] ) correspond to assuming the larger systematic uncertainty , @xmath99% , for the borexino measured rate .    for @xmath100
, the average uncertainty in the combined kamland and borexino measurement is expected to be @xmath101% , within the 1 @xmath76 lma region and if the sterile admixture is as large as currently allowed ,
in this section , we analyze three strategies for determining the total @xmath10 solar neutrino flux without requiring a dedicated experiment that measures only the @xmath10 neutrinos . we first describe in sec .
[ subsec : existingdata ] how one can make a crude determination of the @xmath10 neutrino flux using the measured gallium , chlorine , and sno ( cc ) event rates  @xcite and the standard solar model predictions for all but the @xmath0b , @xmath8be , and @xmath10 solar neutrino fluxes .
this part of the discussion is similar to the analysis described in ref .
@xcite , although we evaluate explicitly the uncertainty caused by the finite size of the allowed region in oscillation parameter space .
we then determine in sec .
[ subsec : galliumsno ] how well one can infer the @xmath10 flux using just the gallium and the sno measurements and the bp00 predictions for the other neutrino fluxes , especially the @xmath8be neutrino flux .
the unknown sterile - active mixture contributes only a negligible error using the strategies described in sec .
[ subsec : existingdata ] and sec .
[ subsec : galliumsno ] .
next , we show in sec .
[ subsec : ppkamlandborexino ] how the precision of the determination of the @xmath10 flux can be improved by using , in the future , the results from kamland ( to constrain the neutrino oscillation parameters ) and from borexino ( to constrain the @xmath8be neutrino flux ) .
the principle of this strategy has also been discussed by the sage collaboration  @xcite .
the theoretical error on the @xmath10 neutrino flux is @xmath103% ( see ref .
@xcite ) , which is an order of magnitude smaller than the estimated errors that we find in this section on the empirically derived values of @xmath104 . to achieve a precision of @xmath4% or better in the determination of the solar neutrino @xmath10 flux will require a dedicated and accurate experiment devoted to measuring the @xmath10 flux .    throughout this section
, we treat the @xmath10 flux and the @xmath38 flux as a single variable because they are so closely linked physically  @xcite .
the two reactions are linked because the @xmath38 reaction is obtained from the @xmath10 reaction by exchanging a positron in the final state with an electron in the initial state .
thus the rates for the two reactions are proportional to each other to high accuracy  @xcite . for convenience ,
we denote the sum of the @xmath10 and @xmath38 fluxes as simply @xmath10 .      in this subsection ,
we show how data from the chlorine , gallium , and sno experiments can be combined with the bp00 predictions for the cno and @xmath105 fluxes to determine the total @xmath10 neutrino flux .
the reduced @xmath10 solar neutrino flux , @xmath104 defined ( by analogy with @xmath54 ) with respect to the predicted standard solar model sum of the @xmath10 and @xmath38 fluxes , can be written as : @xmath106 & & - f_{\rm be } r^{\rm ^7be , ssm}_{\rm ga }   \langle p_{ee}(\delta m^2,\tan^2\theta)\rangle^{\rm ^7be}_{\rm ga } -r^{\rm^8b , ssm}_{\rm ga } \frac{r^{\rm cc , exp}_{\rm sno } } { r^{\rm cc , ssm}_{\rm sno } } \frac{\langle p_{ee}(\delta m^2,\tan^2\theta)\rangle^{\rm ^8b}_{\rm ga } } { \langle p_{ee}(\delta m^2,\tan^2\theta)\rangle^{\rm ^8b}_{\rm sno } } \big ] .\label{eq : fpp}\end{aligned}\ ] ] here the sum over @xmath107 in eq .
( [ eq : fpp ] ) refers to the three cno neutrino sources and the @xmath105 neutrinos ( see table  [ tab : garates ] ) .
we insert in eq .
( [ eq : fpp ] ) the expression for @xmath56 in terms of the chlorine and sno rates and the standard solar model cno and @xmath105 neutrino fluxes .
explicitly , @xmath108   & & -r^{\rm ^8b , ssm}_{\rm cl } \frac{r^{\rm cc , exp}_{\rm sno } } { r^{\rm cc , ssm}_{\rm sno}}\frac { \langle p_{ee } ( \delta m^2,\tan^2\theta)\rangle^{\rm ^8b}_{\rm cl } } { \langle p_{ee } ( \delta m^2,\tan^2\theta)\rangle_{\rm sno } } \big ] . \label{eq : fppbechlorine}\end{aligned}\ ] ]    for the special case of the best fit point in the lma solution region , we find @xmath109 the first error in eq .
( [ eq : fppbestlma ] ) , ( @xmath110 ) , results from the weighted average experimental error for the gallium experiment  @xcite .
the second error , ( @xmath111 ) , reflects the uncertainties in the calculated neutrino absorption cross sections on gallium  @xcite .
the third error , ( @xmath112 ) , is caused by the uncertainties in the calculated standard solar model cno fluxes  @xcite .
the fourth error , ( @xmath113 ) , derives from the measurement errors in the sno cc experimental data  @xcite and neutrino absorption cross section  @xcite .
the fifth error , ( @xmath114 ) , results from the experimental error for the chlorine event rate  @xcite and the sixth error , ( @xmath115 ) , reflects the uncertainties in the calculated neutrino absorption cross sections for chlorine  @xcite .
the physical constraint , @xmath116 , was imposed in obtaining these errors .    within the @xmath2 lma allowed region ,
we find that the central value of @xmath104 can vary in the range @xmath117 and the best - fit value can be determined with an average uncertainty of @xmath118% .
the range of central values of @xmath104 always exceeds the bp00 value because the observed chlorine rate implies a low value of @xmath119 .
since the largest contributions to the gallium rate are from @xmath10 neutrino and @xmath8be neutrinos , a low @xmath8be flux implies a relatively high @xmath10 flux .
we can summarize as follows the results of the simulations carried out within the @xmath2 allowed lma region for the @xmath10 flux : @xmath120 the first error in eq .
( [ eq : fppsummaryknown ] ) results from the allowed range of neutrino oscillation parameters and the second error results from the uncertainties in all other recognized sources of errors , combined quadratically .
the result given in eq .
( [ eq : fppsummaryknown ] ) should be compared with the quoted estimate by the sage collaboration  @xcite , @xmath121 .
this agreement is very welcome since the sage paper points out that they made  several approximations ... whose nature can not be easily quantified . "      in this subsection , we determine the range of allowed values for @xmath104 using the bp00 predictions ( and uncertainties ) for the @xmath8be neutrinos as well as the cno and @xmath105 neutrinos .
we use in this subsection data only from the gallium and sno solar neutrino experiments .
following the same line of reasoning as in section  [ subsec : existingdata ] , we find for the best fit point in the lma allowed oscillation region , @xmath122 just as in eq .
( [ eq : fppbestlma ] ) , the first error in eq .
( [ eq : fbbestlmabe7 ] ) , ( @xmath123 ) , is the experimental error from the weighted average event rate in the gallium experiments , the second error,(@xmath124 ) , is due to the neutrino absorption cross section on gallium , the third error , ( @xmath125 ) , is due to the uncertainties in the bp00 predictions of the cno neutrino fluxes , the fourth error , ( @xmath126 ) , contains the uncertainty in the sno cc experimental rate and the calculated absorption cross sections on deuterium .
the uncertainty in the bp00 prediction for the @xmath8be neutrino flux contributes the last error , ( @xmath127 ) , in eq .
( [ eq : fbbestlmabe7 ] ) .    within
the @xmath2 allowed lma region , we find the central value of @xmath104 varies in the range @xmath128 ; the best - fit value of @xmath104 can be determined with an average uncertainty of @xmath129% .
we can summarize the determination of @xmath104 as follows : @xmath130 where the first error is the uncertainty due to the neutrino oscillation parameters and the second error contains the uncertainties from all other sources .      in this subsection
, we show how the uncertainty in @xmath104 could be improved by using the kamland data to determine the neutrino oscillation parameters and the borexino data , together with the kamland oscillation parameters , to constrain @xmath131 .
we start with eq .
( [ eq : fpp ] ) but now use @xmath56 as determined from borexino and kamland data [ eq .  ( [ eq : fborexino ] ) ] .
we find that for the best - fit lma solution ( assuming that borexino finds the rate expected for this best fit point ) , @xmath132 the first error shown in eq .
( [ eq : fppkamlandborexino ] ) , ( @xmath133 ) , is the experimental error on the measured gallium event rate , the second error , ( @xmath124 ) , represents the uncertainties from the calculated gallium cross sections .
the third error , @xmath134 , is from the predicted cno fluxes and the fourth error , ( @xmath126 ) , contains the uncertainties from the sno experimental data and calculated deuterium absorption cross sections .
the range of oscillation parameters within the kamland reconstructed region gives the fifth error , ( @xmath135 ) , and the sixth error , ( @xmath136 [ @xmath137 ) , results from the uncertainty in the measured borexino event rate [ which we take to be 5% ( 10% ) ] .
the unknown active - sterile admixture ( which goes in the direction of lowering @xmath104 ) contributes the last error shown in eq .
( [ eq : fppkamlandborexino ] ) .    on average , the
the precision expected using the gallium , kamland , and sno experiments is @xmath138 .
\label{eq : avfppborexinokamland}\ ] ] the dominant sources of error in eq .
( [ eq : avfppborexinokamland ] ) are the experimental error in the gallium rate and the uncertainties in the calculated gallium absorption cross section .
in this section , we review and discuss our principal conclusions .
we begin in sec .
[ subsec : discusscurrentlyallowed ] by summarizing the experimentally - allowed range that currently exists for the total @xmath0b solar neutrino flux , taking account of the possibility that sterile neutrinos exist .
we then summarize in sec .
[ subsec : discussb8 ] how well the total @xmath0b neutrino flux , and separately the sterile @xmath0b neutrino flux , can be obtained by combining kamland and sno measurements .
next we describe in sec .
[ subsec : discussbe7 ] how well the total @xmath8be flux can be determined using data from the kamland , gallium ( gallex , sage , and gno ) , and sno experiments , i.e. , using only cc disappearance experiments . in this same subsection , we summarize the measurement accuracy that can be obtained using data from just the kamland and the borexino ( @xmath39 scattering ) experiments .
we describe three strategies ( making use of existing and future solar neutrino measurements and guidance from the bp00 solar model ) for determining the @xmath10 solar neutrino flux without the help of an experiment that measures separately the @xmath10 neutrinos .
we point out in sec .
[ subsec : discusscrosssections ] ( and in the appendix ) that in order to determine the total @xmath0b or @xmath8be neutrino flux , and _ a fortiori _ to determine the sterile contributions to these fluxes , the correlations among theoretical errors for neutrino absorption cross sections and for solar neutrino fluxes must be treated more accurately than in previous discussions .
section  [ subsec : discussfocus ] summarizes the main focus of the present paper .
we concentrate on procedures to determine experimentally the total solar @xmath0b neutrino flux and the total solar @xmath8be neutrino flux in a universe in which sterile neutrinos might exist .
our methods can work only if the lma solution of the solar neutrino problem is correct .
the numerical values given here for the precision with which different quantities can be determined are obtained using simulations of how well different experiments may perform .
therefore , the numerical values are intended only as illustrative guides to what may be possible .      the combined sno cc data and the super - kamiokande @xmath83 scattering data together yield a widely acclaimed agreement between the @xmath0b solar neutrino flux that is predicted by the standard solar model and the empirically inferred @xmath0b neutrino flux .
however , this agreement between solar model prediction and solar neutrino measurement is not a unique interpretation of the existing measurements if one allows for the possibility that the incident solar neutrino flux could contain a significant component of sterile neutrinos .
we show in table  [ tab : fbranges ] and in sec .
[ sec : presentknowledge ] that if one takes account of the possibility that sterile neutrinos may exist then the total solar @xmath0b neutrino flux could be as large as @xmath139 times the flux predicted by the standard solar model . in principle , the existing solar neutrino data could be inconsistent with the standard solar model predictions .
the total @xmath0b solar neutrino flux , active plus sterile neutrinos , can be determined with a typical @xmath2 accuracy of about @xmath4% by comparing the charged current measurements from the kamland reactor experiment and the sno experiments ( see sec .
[ sec : determiningb8total ] ) .
the active @xmath0b neutrino flux has been measured by comparing the sno cc flux ( which measures @xmath40 ) with the super - kamiokande neutrino - electron scattering rate ( which measures @xmath40 plus , with less sensitivity , @xmath140 ) .
the sno neutral current measurement will provide an additional and , ultimately , more accurate measurement of the total active @xmath0b solar neutrino flux .    by subtracting the independently - measured active @xmath0b neutrino flux from the total ( active plus sterile ) @xmath0b neutrino flux , one can determine empirically the sterile component of the solar neutrino flux .
the active @xmath0b neutrino flux will eventually be determined accurately by the sno neutral current measurement  @xcite .
we estimate therefore that the procedure described here has the potential of measuring , or setting an upper limit on , the sterile component of the @xmath0b neutrino flux that is as small as @xmath141% of the total @xmath0b solar neutrino flux .
the measurement of the total @xmath0b neutrino flux , and the sterile component of this flux , are independent of solar model considerations . in order to establish the quantitative conclusions ,
we have performed detailed simulations of the accuracy of the kamland reactor experiment in determining neutrino oscillation parameters ( see fig .
[ fig : solarvskamlandreal ] ) and have evaluated the theoretical and experimental uncertainties that affect the different flux determinations ( see sec .
[ subsec : precisionb8 ] and the appendix ) .
the combined measurements of the super - kamiokande and kamland experiments can be used to determine independently a value for the total @xmath0b neutrino flux .
this determination may be as accurate as @xmath6% for purely active neutrinos . with the current limits on the active - sterile admixture
, the total @xmath0b neutrino flux could be inferred to an accuracy of @xmath7% or better , as described in sec .
[ subsec : kamlandsuperkb8 ] .
it will be important to compare the value of the total @xmath0b neutrino flux inferred by combining the kamland and sno charged current measurements with the value obtained using the kamland and super - kamiokande experiments
. this comparison will be an important test of whether the systematic uncertainties in the experiment and in the analyses are understood .
the total @xmath8be solar neutrino flux , active plus sterile neutrinos , can be determined to a @xmath2 accuracy of about @xmath142% by combining measurements from kamland , sno , and the gallium experiments ( see sec .
[ sec : be7 ] ) .
unlike the purely empirical determination that is possible for the @xmath0b flux , the measurement of the total @xmath8be solar neutrino flux requires some assumption regarding the cno solar neutrino fluxes . in our estimates
, we have assumed that the standard solar model predictions for the cno fluxes , and their uncertainties , are at least approximately valid .
table  [ tab : fbe ] shows that as long as the cno fluxes are not a factor of three or more larger than the standard solar model predictions , then they will not significantly limit the accuracy with which the total @xmath8be neutrino flux can be determined . the measured capture rate in the gallium experiments ( gallex , sage , and gno ) currently constitutes the largest recognized uncertainty in the determination of the total @xmath8be flux by the method described here ( see table  [ tab : fbe ] ) .
the constraints provided by the chlorine experiment are not very useful in providing an accurate determination of the @xmath8be neutrino flux [ see table  [ tab : fbecl ] and eq .
( [ eq : averagebe7shiftcl ] ) ] .
one can also determine the allowed range of the total @xmath8be solar neutrino flux using the data from the kamland reactor experiment and the borexino solar neutrino experiment .
we show in sec .
[ subsec : kamlandborexino ] that with this method one may hope to obtain a @xmath2 accuracy of @xmath1% or better for the total @xmath8be solar neutrino flux , which is more accurate than we estimate will be possible with the gallium and chlorine radiochemical experiments .
one can determine an upper limit to the sterile component of the @xmath8be solar neutrino flux by combining the measured rate in the neutrino - electron scattering experiment , borexino , with the @xmath8be total flux inferred from measurements of the kamland , sno , and gallium experiments .
if one assumes that the entire signal measured in borexino is due to @xmath40 , then one obtains a minimum value for the active component of the @xmath8be neutrino flux . subtracting this minimum value from the total @xmath8be flux ,
one will obtain an upper limit to the sterile component of the flux .
it seems unlikely that the procedure described here has the sensitivity to measure a value for the sterile component unless the sterile flux is larger than @xmath142% of the total @xmath8be neutrino flux .
however , limits on the sterile neutrino admixture can be obtained from the analysis of @xmath0b and kamland neutrino measurements described in sec .
[ subsec : discussb8 ] .    in order to make a direct and precise measurement of the sterile component of the @xmath8be solar neutrino flux , we need a charged current measurement of the @xmath8be flux .
a @xmath8be solar neutrino absorption experiment , e.g. , with a lithium target  @xcite or with lens  @xcite , would make possible an accurate determination of the sterile @xmath8be neutrino flux by providing a set of experimental constraints that is analogous to what will exist for the sno , super - kamiokande , and kamland experiments .      the theoretical uncertainty in the calculated @xmath10 solar neutrino flux is estimated to be only @xmath143  @xcite .
therefore , a precise determination of the @xmath10 solar neutrino flux will be of great interest as a crucial test of the theory of stellar evolution .
the measurement of the @xmath10 neutrino flux will also provide a critical test of whether the neutrino oscillation theory , which works well at energies above @xmath95 mev , also describes accurately the lower energy neutrino phenomena ( energies less than @xmath144 mev ) .
we show in sec .
[ sec : determiningpp ] that the @xmath10 flux can be determined with modest accuracy , of order @xmath145% to @xmath146% , using a combination of existing experimental data and some guidance from the bp00 standard solar model . in the future , it should be possible to determine the @xmath10 neutrino flux to an accuracy of @xmath11% using experimental data from the borexino , kamland , and sno solar neutrino experiments and predictions ( in a non - critical way ) from the standard solar model ( see eq .
[ eq : avfppborexinokamland ] ) .
the unknown flux of sterile neutrinos does not significantly affect the quoted estimates on the accuracy with which the @xmath10 flux can be determined . to measure the @xmath10 flux with an accuracy sufficient to test
stringently the standard solar model prediction will require a dedicated and accurate experiment that measures separately the @xmath10 neutrino flux .      in order to determine the total @xmath0b or @xmath8be solar neutrino flux , and their sterile components
, one must evaluate carefully all known sources of error .
in the course of this investigation , we realized that the evaluations of the neutrino absorption cross section uncertainties in previous neutrino oscillation studies , including our own , have not properly taken account of the correlations among the theoretical uncertainties in the cross section calculations ( for an insightful discussion of this point , see ref .
we discuss in the appendix how the cross section errors can be treated more correctly .
we also emphasize here that it is necessary to treat as fully correlated the uncertainties in the principal cno neutrino fluxes that are obtained from standard solar model predictions .
these effects are small compared to other uncertainties in determining solar neutrino oscillation parameters ( the principal goal of nearly all previous oscillation studies ) .
the correlations among cross section uncertainties become important only when one wants to make accurate inferences regarding the neutrino fluxes themselves .
the main focus in this paper is on determining experimentally the total @xmath0b , @xmath8be , and @xmath10 solar neutrino fluxes in order to make possible more precise tests of solar model predictions .
we have also shown that the contribution of sterile neutrinos to the total flux can be measured , or a useful upper limit can be set , for @xmath0b solar neutrinos and for @xmath8be solar neutrinos .
jnb acknowledges support from nsf grant no .
mcg - g is supported by the european union marie - curie fellowship hpmf - ct-2000 - 00516 .
this work was also supported by the spanish dgicyt under grants pb98 - 0693 and fpa2001 - 3031 .
we determine , as is conventional for many analyses of solar neutrino data , the allowed regions in the neutrino oscillation space using a @xmath147 function that includes all the relevant data . in the construction of the @xmath147 function
, we have followed closely the prescription of ref .
@xcite ( see also ref .
@xcite ) , but we have included some modifications to this prescription in order to account in more detail for the energy dependence and the correlation of the cross section errors for the chlorine and gallium solar neutrino experiments .    for a given experiment @xmath148 ( for example , gallium or chlorine ) , the expected number of events can be written as @xmath149 where @xmath150 labels the solar neutrino fluxes , @xmath89 , and @xmath151 labels the energy bins of energy @xmath152 .
the quantity @xmath153 is the cross section for the interaction of a neutrino of energy @xmath152 in the experiment @xmath148 ; @xmath154 is the survival probability of @xmath40 for a given energy @xmath152 at a specified point , @xmath155 , in neutrino oscillation space , because different neutrino species have different probability distributions for the location of their production within the sun . for matter oscillations ,
the survival probability obviously depends upon where in the sun the neutrino was produced . ] .
the error matrix for the neutrino absorption cross sections can be derived using eq .
( [ eqapp : rate ] ) .
since the errors are uncorrelated between different experiments , we need to evaluate the following expression for a specified experiment : @xmath156 where @xmath157 is best - estimate for the ( gallium or chlorine ) experimental capture rate and the brackets indicate an average over the probability distribution of uncertainties in the neutrino cross sections .
writing out the various terms in eq .
( [ appeq : defnsigma ] ) , we find @xmath158 in the standard ref .
@xcite , the cross section error matrix for each experiment is given as : @xmath159 where @xmath160 is the average uncertainty in the interaction cross section of neutrinos of flux type @xmath107 in experiment @xmath148 .
equation  ( [ sigma2 ] ) is obtained from eq .
( [ eqapp : rate ] ) by neglecting the correlations of the cross section errors among the different neutrino fluxes and by neglecting the energy dependence of the cross section errors . for the gallium detector ,
neither of these assumptions is correct and even for chlorine the cross section errors for different neutrino sources are strongly correlated .
for all but the @xmath0b and @xmath105 neutrino fluxes , either ground state to ground state transitions are the only energetically possible transitions ( which is the case for the chlorine detector ) or the ground state to ground state transitions dominate ( which is the case for the gallium detector ) . thus all the cross sections for the lower energy neutrinos move up or down together , proportional to the square of the dominant matrix element .
when considering the gallium experiments , we include the energy dependence of the cross section errors and assume full correlations among the cross section errors of the different neutrino sources that contribute at a specific energy .
specifically , we use the results from ref .
@xcite for the cross section errors . for continuum sources ,
we take the fractional error of the interaction cross section in the energy bin @xmath151 to be @xmath161,\ ] ] where @xmath162 are the 3@xmath76 upper and lower limit cross sections given in tables iii and iv of ref .
@xcite . for the line sources , @xmath8be and @xmath38
, we use the errors given in eqs .
( 41 ) and ( 42 ) of ref .
@xcite . to be complete ,
we have checked that the shapes of these neutrino lines ( see ref .
@xcite ) do not affect significantly the error estimates .
the gallium contribution to the cross section error matrix is therefore given by @xmath163 where @xmath164 is the correlation matrix for the cross section errors of the different energy bins . for low energy neutrino sources , the dominant process is the transition to the germanium ground state and for higher energy neutrino sources ( @xmath165 mev ) transitions to excited states dominate .
therefore , we assume the cross section errors to be fully correlated between energy bins either below or above @xmath166 mev ( @xmath167 for @xmath168 mev or @xmath169 mev ) .
we take the errors to be uncorrelated between one energy bin below @xmath170 mev and one energy bin above @xmath166 mev ( @xmath171 for @xmath172 mev and @xmath173 mev , or @xmath172 mev and @xmath173 mev ) .    for chlorine
, we take as fully correlated the cross section errors for the @xmath10 , @xmath38 , @xmath8be , and cno neutrinos .
the cross section errors for the @xmath0b and @xmath105 neutrinos are uncorrelated with the errors for the lower energy neutrinos but are fully correlated with each other .
we can neglect for chlorine the energy dependence of the chlorine to argon cross section errors because the uncertainty for the lower energy neutrinos is determined almost entirely from the ground - state to ground - state matrix element .
( forbidden corrections are unimportant for these low energy neutrinos . ) for the @xmath0b and @xmath105 neutrinos , the absorption cross sections are also dominated by a single ( but different ) transition , in this case the super - allowed transition  @xcite , and therefore we can also neglect the energy dependence for the higher energy neutrinos .
we adopt the values of the averaged chlorine cross section errors given in refs .
@xcite . to be explicit , the chlorine contribution to the cross section error matrix is given by @xmath174 where @xmath175 for @xmath176 or @xmath177 .
also , @xmath178 for @xmath179 and @xmath180 .
99 q.  r.  ahmad _ et al .
_ , sno collaboration , nucl - ex/0204008
. q.  r.  ahmad _ et al .
_ , sno collaboration , nucl - ex/0204009
. j. n. abdurashitov _
et al . _ ,
sage collaboration , astro - ph/0204245 .
q.  r.  ahmad _ et al .
_ , sno collaboration , phys .
lett .   * 87 * , 071301 ( 2001 ) .
s.  fukuda _ et al .
_ , super - kamiokande collaboration , phys.rev .
lett .   * 86 * , 5651 ( 2001 ) .
g.  l.  fogli , e.  lisi , d.  montanino , and a.  palazzo , phys .
d * 64 * , 093007 ( 2001 ) . c. giunti , phys .
d * 65 * , 033006 ( 2002 ) ; m.  v.  garzelli and c.  giunti , phys .
d * 65 * , 093005 ( 2002 ) .
g. fiorentini , f. l. villante , and b. ricci , hep - ph/0109275 .
j. n. bahcall , phys .
c. * 65 * , 015802 ( 2002 ) .
j.  n.  bahcall , phys .
* 12 * , 300 ( 1964 ) ; r.  davis jr .
* 12 * , 303 ( 1964 ) ; j.  n.  bahcall and r.  davis , science * 191 * , 264 ( 1976 ) . j.  n.  bahcall , m.  h.  pinsonneault , and s.  basu , astrophys .  j.   * 555 * , 990 ( 2001 ) .
b. t. cleveland , t. daily , r. davis , jr . , j. r. distel , k. lande , c. k. lee , and p. s. wildenhain , astrophys .
j. * 496 * , 505 ( 1998 ) .
y. fukuda _
_ , kamiokande collaboration , phys .
lett . * 77 * , 1683 ( 1996 ) .
j. n. abdurashitov _
et al . _ ,
sage collaboration , phys .
c * 60 * , 055801 ( 1999 ) .
w.  hampel _ et al .
_ , gallex collaboration , phys .
b * 447 * , 127 ( 1999 ) .
e. belloti , nucl .
b ( proc . suppl . )
* 91 * , 44 ( 2001 ) .
( ` http://nu2000.sno.laurentian.ca ` ) .
j.  n.  bahcall , s.  basu , and m.  pinsonneault , phys .
b * 433 * , 1 ( 1998 ) .
w. kwong and s. p. rosen , phys .
d * 54 * , 2043 ( 1996 ) .
f.  l.  villante , g.  fiorentini , and e.  lisi , phys .
d * 59 * , 013006 ( 1999 ) .
j.  n.  bahcall , m.  c.  gonzalez - garcia , and c.  pena - garay , j. high energy phys .
* 08 * , 014 ( 2001 ) .
p.  i.  krastev and a.  y.  smirnov , phys .
d * 65 * , 073022 ( 2002 ) .
j.  n.  bahcall , m.  c.  gonzalez - garcia , and c.  pena - garay , j. high energy phys .
* 04 * , 007 ( 2002 ) .
a. r. junghans , e. c. mohrmann , k. a. snover , t. d. steiger , e. g. adelberger , j. m. casandjian , h. e. swanson , l. buchmann , s. h. park , and a. zyuzin , phys .
* 88 * , 041101 ( 2002 ) .
s. m. bilenky , c. giunti , and w. grimus , prog .
* 43 * , 1 ( 1999 ) .
m. c. gonzalez - garcia and y. nir , hep - ph/0202058 .
v.  barger , d.  marfatia , and k.  whisnant , phys .
lett . * 88 * , 011302 ( 2002 ) ; v.  barger , d.  marfatia , and k.  whisnant , phys .
b * 509 * , 19 ( 2001 ) .
d. dooling , c. giunti , k. kang , and c. w. kim , phys .
d * 61 * , 073011 ( 2000 ) . c.  giunti , m.  c.  gonzalez - garcia , and c.  pena - garay , phys .
d * 62 * , 013005 ( 2000 ) .
j. aguilar _ et al .
_ , lsnd collaboration , phys . rev .
d * 64 * , 112007 ( 2001 ) .
k. eitel , karmen collaboration , nucl .
b ( proc . suppl . ) * 91 * , 191 ( 2001 ) .
a. piepke , kamland collaboration , nucl .
b ( proc . suppl . )
* 91 * , 99 ( 2001 ) ; j.  busenitz _ et al . , _ , kamland collaboration , `` proposal for us participation in kamland , '' march 1999 ( unpublished ) .
may be downloaded from:`http://kamland.lbl.gov/ ` ; k.  ishihara for the kamland collaboration , talk at the nufact01 workshop in tsukuba , japan ( may 2430 , 2001 ) . transparencies at:`http://psux1.kek.jp/~nufact01/index.html ` .
a.  bandyopadhyay , s.  choubey , s.  goswami , and k.  kar , phys .
b * 519 * , 83 ( 2001 ) .
v.  barger , d.  marfatia , and b. p.  wood , phys .
b * 498 * , 53 ( 2001 ) .
r.  barbieri and a.  strumia , j. high energy phys .
* 12 * , 016 ( 2000 ) . h.  murayama and a.  pierce , phys .
d * 65 * , 013012 ( 2002 ) .
j.  n.  bahcall , m.  c.  gonzalez - garcia , and c.  pena - garay , hep - ph/0204314 .
m. smy , proceedings of the noon2001 workshop , kashiwa , japan , december 5 - 8 , 2001 , hep - ex/0202020 a.  de gouvea and c.  pena - garay , phys .
d * 64 * , 113011 ( 2001 ) .
m.  c.  gonzalez - garcia and c.  pena - garay , phys .
b * 527 * , 199 ( 2002 ) .
g. alimonti _
et al . _ ,
borexino collaboration , astropart .
* 16 * , 205 ( 2002 ) ; a. piepke _
et al . _ ,
kamland collaboration , nucl .
suppl . ) * 91 * , 99 ( 2001 ) .
s. nakamura , t. sato , s. ando , t .- s .
park , f. myhrer , v. gudkov , and k. kubodera , nucl - th/0201062 .
see also , a.  kurylov , m.  j.  ramsey - musolf , and p.  vogel , phys .  rev .
c * 65 * , 055501 ( 2002 ) ; j.  f.  beacom and s.  j.  parke , phys .  rev .
d * 64 * , 091302 ( 2001 ) .
m. butler , j .- w .
chen , and x. kong , phys .
c * 63 * , 035501 ( 2001 ) .
r. davis jr .
, d. s. harmer , and k. c. hoffman , phys .
lett . * 20 * , 1205 ( 1968 ) .
j.  n.  bahcall , phys .
c * 56 * , 3391 ( 1997 ) .
j.  n.  bahcall phys .
lett . * 71 * , 2369 ( 1993 ) ; phys .
d * 49 * , 3923 ( 1994 ) .
j. n. bahcall and h. a. bethe , phys .
lett . * 65 * , 2233 ( 1990 ) .
j. n. bahcall , m. kamionkowski , and a. sirlin , phys .
d * 51 * , 6146 ( 1995 ) .
j. n. bahcall , _ neutrino astrophysics _
( cambridge university press , cambridge , 1989 ) .
j. n. bahcall and r. m. may , astrophys .
j.   * 155 * , 501 ( 1969 ) .
j. n. bahcall , e. lisi , d. e. alburger , l. de braeckeleer , s. j. freedman , and j. napolitano , phys .
c * 54 * , 411 ( 1996 ) .
j.  n.  bahcall , p.  i.  krastev , and a.  y.  smirnov , phys .
d * 62 * , 093004 ( 2000 ) .
m. galeazzi , g. gallinaro , f. gatti , p. meunier , s. vitale , a. v. kopylov , v. v. petukhov , e. a. yanovich , and g. t. zatsepin , phys .
b * 398 * , 187 ( 1997 ) ; a. kopylov and v. petukhov , hep - ph/0201042 .
r. s. raghavan , phys .
lett . * 78 * , 3618 ( 1997 ) ; r. s. raghavan , hep - ex/0106054 . m.  v.  garzelli and c.  giunti , phys .
b * 488 * , 339 ( 2000 ) .
g. l. fogli and e. lisi , astropart . phys . * 3 * , 185 ( 1995 ) with updated uncertainties as discussed in g. l. fogli , e. lisi , d. montanino , and a. palazzo , phys .
d * 62 * , 013002 ( 2000 ) .
j.  n.  bahcall , phys .
c * 56 * , 3391 ( 1997 ) .