Abstract

We discuss constraints on the coefficient which is introduced to simulate the effect of weaker or stronger matter potential for electron neutrinos with the current and future solar neutrino data. The currently available solar neutrino data leads to a bound at 1 (3) CL, which is consistent with the Standard Model prediction . For weaker matter potential (), the constraint which comes from the flat 8B neutrino spectrum is already very tight, indicating the evidence for matter effects. However for stronger matter potential (), the bound is milder and is dominated by the day-night asymmetry of 8B neutrino flux recently observed by Super-Kamiokande. Among the list of observables of ongoing and future solar neutrino experiments, we find that (1) an improved precision of the day-night asymmetry of 8B neutrinos, (2) precision measurements of the low-energy quasi-monoenergetic neutrinos, and (3) the detection of the upturn of the 8B neutrino spectrum at low energies are the best choices to improve the bound on .

1. Introduction

Neutrino propagation in matter is described by the Mikheyev-Smirnov-Wolfenstein (MSW) theory [14]. It was successfully applied to solve the solar neutrino problem [5], the discrepancy between the data [615], and the theoretical prediction of solar neutrino flux [16], which blossomed into the solution of the puzzle, the large-mixing-angle (LMA) MSW solution. The solution is in perfect agreement with the result obtained by KamLAND [17] detector which measured antineutrinos from nuclear reactors, where the flavor conversion corresponds to vacuum oscillations with subpercent corrections due to matter effects.

The MSW theory relies on neutrino interaction with matter dictated by the standard electroweak theory and the standard treatment of refraction which is well-founded in the theory of refraction of light. Therefore, it is believed to be on a firm basis. On the observational side it predicts a severer reduction of the solar neutrino flux at high energies due to the adiabatic flavor transition in matter than at low energies where the vacuum oscillation effect dominates. Globally, the behavior is indeed seen in the experiments observing 8B solar neutrinos at high energies [1114] and in radiochemical experiments detecting low-energy pp and 7Be neutrinos [610] and more recently by the direct measurement of 7Be neutrinos by Borexino [15]. For a summary plot of the current status of high- and low-energy solar neutrinos, see the review of solar neutrinos in this series. Therefore, one can say that the MSW theory is successfully confronted with the available experimental data.

Nevertheless, we believe that further test of the MSW theory is worth pursuing. First of all, it is testing the charged current (CC) contribution to the index of refraction of neutrinos of the Standard Model (SM), which could not be tested anywhere else. Furthermore, in analyses of future experiments to determine and the mass hierarchy, the MSW theory is usually assumed to disentangle the genuine effect of CP phase from the matter effect. Therefore, to prove it to the accuracy required by measurement of is highly desirable to make discovery of CP violation robust in such experiments that could have matter effect contamination. This reasoning was spelled out in [18]. Since the survival probability does not depend on solar neutrinos provide us with a clean environment for testing the theory of neutrino propagation in matter.

We notice that in solar neutrinos, the transition from low- to high-energy behaviors mentioned above has not been clearly seen in a single experiment in a solar-model-independent manner. The Borexino and KamLAND experiments tried to fill the gap by observing 8B neutrinos at relatively low energies [19, 20]. SNO published the results of analyses with lower threshold energy of 3.5 MeV [21, 22], and the similar challenge is being undertaken by the Super-Kamiokande (SK) group [23]. In addition to 7Be, a new low-energy neutrino line, pep neutrinos, was observed by Borexino [24]. Recently, the SK group announced their first detection of the day-night asymmetry of 8B neutrinos [23]. As we will see in Section 3 it gives a significant impact on our discussions. With these new experimental inputs, as well as all the aforementioned ones, it is now quite timely to revisit the question of how large deviation from the MSW theory is allowed by data.

In this paper, we perform such a test of the theory of neutrino propagation in the environments of solar and earth matter. For this purpose, we need to specify the framework of how deviation from the MSW theory is parametrized. We introduce, following [25, 26], the parameter defined as the ratio of the effective coupling of weak interactions measured with coherent neutrino matter interactions in the forward direction to the Fermi coupling constant . We first analyze the currently available solar neutrino data to obtain the constraints on and find that the features of the constraints differ depending upon or . We will discuss interpretations of the obtained constraints including this feature and provide a simple qualitative model to explain the bound at , more nontrivial one. We then discuss the question of to what extent the constraints on can be made stringent by various future solar neutrino observables.

Our framework of testing the theory of neutrino propagation in matter requires comments. It actually involves the three different ingredients: (1) non-SM weak interactions in the forward direction parametrized as , (2) refraction theory of neutrino propagation in matter which includes the resonant enhancement of neutrino flavor conversion [14] and (3) electron number densities in the sun and in the earth. However, on ground of well-founded refraction theory and because no problem can be arguably raised in the formulation of the MSW mechanism we do not question the validity of (2). We also note that the electron number density in the sun is reliably calculated by the standard solar model (SSM), and the result is cross-checked by helioseismology to an accuracy much better than the one discussed here. We can also take the earth matter density and chemical compositions calculated by the preliminary reference earth model (PREM) [27] as granted. It is the case because the earth matter-dependent observable, the day-night variation of solar neutrino flux, is insensitive to the precise profile of the earth matter density. Therefore, we assume that our test primarily examines the aspect (1), namely, whether neutrino matter coupling in the forward direction receives additional contribution beyond those of SM.

Are non-SM weak interactions parametrized as general enough? Most probably not because in many models with new non-SM interactions they have flavor structure. Flavor-dependent new neutral current interactions have been discussed in the framework of nonstandard interactions (NSI) of neutrinos [2833], and constraints on effective neutrino matter coupling were obtained with this setting, for example, in [34, 35]. With solar neutrinos see [36] for discussion of NSI. If we denote the elements of NSI as (), our may be interpreted as , assuming that for , . To deal with the fully generic case, however, we probably have to enlarge the framework of constraining the NSI parameters by including other neutrino sources, in particular, the accelerator and atmospheric neutrinos. See Section 4 for more comments.

2. Simple Analytic Treatment of Matter Effect Dependences

In this section, we give a simple analytic description of how various solar neutrino observables depend upon the matter effect. It should serve for intuitive understanding of the characteristic features which we will see in the later sections. The reader will find a physics discussion in the flavor conversion review of this series. In the following, we denote the matter densities inside the sun and in the earth as and , respectively. Solar neutrino survival or appearance probabilities depend on three oscillation parameters: the solar oscillation parameters (, ) and . Smallness of the recently measured value of [3741] and its small error greatly restrict the uncertainty introduced by this parameter on the determination of matter effects.

To quantify possible deviation from the MSW theory, we introduce the parameter by replacing the Fermi coupling constant by [25, 26]. The underlying assumption behind such simplified framework is that the deviation from the Fermi coupling constant is universal over fermions, in particular up and down quarks.

The survival probability in the absence of the earth matter effect, that is, during the day, is well described by [4244]: Here is the mixing angle at the production point inside the sun: where is the mixing angle in matter of density , In (2.4), is defined as the ratio of the neutrino oscillation length in vacuum, , to the refraction length in matter, : where In (2.4) and (2.5), is the matter density, is the number of electrons per nucleon, and is the nucleon mass. In the last term we have used the best fit of the global analysis  eV. The average electron number densities at the production point of various solar neutrino fluxes are tabulated in Table 1. These numbers serve to show the differences in solar densities probed by the different sources of neutrinos, but the precise calculations are correctly done by averaging the survival probability with the production point distribution of the corresponding source [5, 16, 45].

We observe that in (2.1) depends on neutrino energy and in the particular combination . The property may have the following implications to constraints on : (1) since shifting is equivalent to shifting our analysis which calculates as a function of is inevitably affected by the whole spectrum. (2) Nonetheless, we generically expect that the constraint at () principally comes from neutrino spectrum at high (low) energies. It appears that the apparently contradictory remarks are both true in view of the results in Section 3.

2.1. Energy Spectrum

Solar neutrino observables taken in a single experiment have not shown an energy dependence yet. The neutrino oscillation parameters are such that we cannot expect strong energy dependences. At low neutrino energies, small , (2.1) can be approximated by whereas at high energies, small , the oscillation probability (2.1) can be approximated, keeping only the first energy-dependent term as Notice that the correction to the asymptotic behavior is linear in at low energies while it is quadratic in at high energies. It may mean that the energy spectrum at low energies could be more advantageous in tightening up the constraint on provided that these formulas with leading order corrections are valid.

It is well known that, in the LMA MSW mechanism, 8B neutrino spectrum must show an upturn from the asymptotic high energy (MeV) to lower energies. The behavior is described by the correction term in (2.7) but only at a qualitative level. It indicates that the upturn component in the spectrum is a decreasing function of . On the other hand, at low energies populated by pp, 7Be, and pep neutrinos, the solar neutrino energy spectrum displays vacuum-averaged oscillations or decoherence, (2.6). The deviation from this asymptotic low energy limit can be described by the correction term in (2.6) again at the (better) qualitative level. The term depends upon linearly so that the correction term is an increasing function of . Because of the negative sign in the correction term in (2.6), larger values of lead to smaller absolute values of in both low- and high-energy regions. (The simpler way to reach the same conclusion is to use the property mentioned earlier. Then, for larger () corresponds to the one at higher energy. Since is a monotonically decreasing function of , larger the , smaller the .)

To see how accurate is the behavior predicted by the above approximate analytic expressions, we have computed numerically (using the PREM profile) the average as a function of electron energy. Here, means taking average of over neutrino energies with neutrino fluxes times the differential cross-sections integrated over the true electron energy with response function. In the above expression, with and being the cross-sections of and scattering, respectively. The computed results confirm qualitatively the behavior discussed above based on our analytic approximations. Thus, the energy spectrum of solar neutrinos at low and high energies can constrain in this way, as will be shown quantitatively in Section 3.

2.2. Day-Night Variation

The survival probability at night during which solar neutrinos pass through the earth can be written, assuming adiabaticity, as [46, 47] where is the one given in (2.1). denotes the regeneration effect in the earth and is given as , where is the transition probability of second mass eigenstate to . Under the constant density approximation in the earth, is given by [46, 47] for passage of distance , where we have introduced . In (2.9), and stand for the mixing angle and the parameter (see (2.4)) with matter density in the earth. Within the range of neutrino parameters allowed by the solar neutrino data, the oscillatory term averages to in good approximation when integrated over zenith angle. Then, the equation simplifies to At MeV, which is a typical energy for 8B neutrinos, and for the average density g/cm3 and the electron fraction in the earth. Then, is given as for and . This result is in reasonable agreement with more detailed estimate using the PREM profile [27] for the earth matter density.

We now give a simple estimate of the day-night asymmetry assuming constant matter density approximation in the earth and its dependence. Under the approximation of small regeneration effect , the day-night asymmetry for the CC number of counts NCC measurement is approximately given by where in the right-hand side we have approximated by the asymmetry of survival probabilities in day and in night at an appropriate neutrino energy and ignored the terms of order . Notice that the effects of the solar and the earth matter densities are contained only in and , respectively.

At MeV, , , , and hence , about 3% day-night asymmetry for . Note that for , so that the impact of on gives only a minor modification. Though based on crude approximations, the value of at obtained above is in excellent agreement with the one evaluated numerically for SNO CC measurement.

SNO and SK observe the day-night asymmetry by measurement of CC reactions and elastic scattering (CC + NC), respectively. We have computed as a function of numerically (with PREM profile) without using analytic approximation. The result of scales linearly with in a good approximation, . Similarly, the day-night asymmetry for elastic scattering measurement can be easily computed. Its relationship to the can be estimated in the similar manner as in (2.11), taking into account the modification due to NC scattering. Using approximate values, and , the factor in the square bracket can be estimated to be , giving a reasonable approximation for the ratio of to . A better approximation to the computed results of the dependence of the asymmetry is given by .

3. Constraints on by Solar Neutrino Observables

In this section we investigate quantitatively to what extent can be constrained by the current and the future solar neutrino data. The results of our calculations are presented in Figure 1, supplemented with the relevant numbers in Table 2. We will discuss the results and their implications to some details in a step-by-step manner. We first discuss the constraints by the data currently available (Section 3.1). Then, we address the question of how the constraint on can be tightened with the future solar neutrino data, the spectral upturn of 8B neutrinos (Section 3.2), the low-energy 7Be and pep neutrinos (Section 3.3), and finally the day-night asymmetry of the solar neutrino flux (Section 3.4). We pay special attention to the question of how the constraints on depend upon the significance of these measurements.

3.1. Current Constraint on

We include in our global analyses the KamLAND and all the available solar neutrino data [615, 17, 1923]. To obtain all the results quoted in this paper we marginalize over the mixing angles and , the small mass squared difference , and the solar neutrino fluxes [16, 48] imposing the luminosity constraint [49, 50]. We include in the analysis the dependence derived from the analysis of the atmospheric, accelerator, and reactor data included in [51] as well as the recent measurement of by [3741]. The used is defined by where Marg implies to marginalize over the parameters shown but not over . Further details of the analysis methods can be found in [48].

The currently available neutrino data (blue solid line), which include SNO lower-energy threshold data [21, 22] and SK IV [23], do not allow a very precise determination of the parameter. A distinctive feature of the parabola shown in Figure 1 is the asymmetry between the small and large regions. At the parabola is already fairly steep, and the “wall” is so stiff that can barely be changed by including the future data. While at the slope is relatively gentle. More quantitatively, at 1 (3) CL. The best fit point with the present data is significantly larger than unity, . It was 1.32 before and has driven to the larger value mostly by the new SK data which indicates a stronger matter effect than those expected by the MSW LMA region preferred by the KamLAND data. The larger best fit value could also partly be due to an artifact of the weakness of the constraint in region. Notice that the Standard Model MSW theory value is off from the 1 region but only by a tiny amount, as seen in Table 2. Let us understand these characteristics.

The lower bound on mostly comes from the SK and the SNO data which shows that 8B neutrino spectrum at high energies is well described by the adiabatic LMA MSW solution (). The energy spectrum is very close to a flat one with which can be approximated by with corrections due to the contribution of the energy-dependent term (see (2.7)). The value is inconsistent with the vacuum oscillation, and hence the point is highly disfavored, showing the evidence for the matter effect.

One would think that the upper bound on should come from either the low energy solar neutrino data or the deviation from the flat spectra at high energies. But, we still lack precise information on low-energy solar neutrinos, and the spectral upturn of 8B neutrinos has not been observed beyond the level in [19, 20]. Then, what is the origin of the upper bound at about CL? We argue that it mainly comes from the day-night asymmetry of 8B neutrino flux which is contained in the binned data of SK and SNO. Recently, the SK collaboration reported a positive indication of the day-night asymmetry though the data is still consistent with no asymmetry at 2.3 CL [23].

To show the point, we construct a very simple model for for the day-night asymmetry . It is made possible by the approximate linearity of to . Let us start from the data of day-night asymmetry at SK I–IV obtained with the D/N amplitude method [23]: %, giving the total error 1.2% if added in quadrature. The expectation of by the LMA solution is % for  eV2. Then, one can create an approximate model as .

Despite admittedly crude nature it seems to capture the qualitative features of with the current data (blue solid line) in Figure 1 in region . It is true that it predicts a little too steep rise of and leads to at , whereas in Figure 1. In the actual numerical analysis for Figure 1, however, parabola can naturally become less steep because various other parameters are varied to accommodate such a large values of . Therefore, we find that about 2 evidence of in the SK data is most likely the main cause of the sensitivity to in the region . The simple model cannot explain the behavior of in region in Figure 1, because the other more powerful mechanism is at work to lead to stronger bound on , as discussed above.

To what extent an improved knowledge of affects ? It was suggested that a dedicated reactor neutrino experiment can measure to 2% accuracy [52, 53]. It is also expected that precision measurement of pp spectrum could improve the accuracy of determination to a similar extent [48]. Therefore, it is interesting to examine to what extent an improved knowledge of affects the constraint on . Therefore, we recompute the curves presented in Figure 1 by adding the artificial term in the assuming 2% accuracy in determination. The result of this computation is given in Table 2 in parentheses. As we see, size of the effect of improved knowledge is not very significant.

3.2. Spectrum of Solar Neutrinos at High Energies

Evidence for the upturn of 8B neutrino spectrum must contribute to constrain the larger values of because could be very large without upturn, if day-night asymmetry is ignored. We discuss the impact on of seeing the upturn in recoil electron energy spectrum with 3 significance, which we assume to be in the region MeV. To calculate we assume the errors estimated by the SK collaboration [23]. Adding the simulated data to the currently available dataset produces the magenta dashed line in Figure 1. We find a 25% reduction of the 3 allowed range, at 1 (3) CL. We can see that it does improve the upper bound on , for which the current constraint (blue solid line) is rather weak, but the improvement in the precision of is still moderate.

Some remarks are in order about the minimum point of . The best fit point with the present data is at as we saw above. For the analysis with future data discussed in this and the subsequent subsections, we assume that the minimum is always at for simulated data. Therefore, the analysis with the present plus simulated data tends to pull the minimum toward smaller values of , and at the same time make the parabola narrower around the minimum. By conspiracy between these two features the current constraint (blue solid line) is almost degenerate to the other lines at , the ones with spectral upturn (magenta dashed line) and low-energy neutrinos (red dash-dotted line). These features can be observed in Figure 1 and in Table 2.

3.3. Spectrum of Solar Neutrinos at Low Energies

Now, let us turn to the low-energy solar neutrinos, 7Be and pep lines. The Borexino collaboration has already measured the 7Be neutrino-electron scattering rate to an accuracy of ±5% [15], which we assume throughout this section. For pep neutrinos we assume measurement with 3% precision in the future. See [24] for the first observation of pep neutrinos and its current status of the uncertainties.

The measurement of the pep flux has two important advantages, when compared to the 7Be flux, in determining : (a) the neutrino energy is higher, 1.44 MeV, so the importance of the solar matter effects is larger and (b) the uncertainty in the theoretical estimate is much smaller. Firstly, the ratio of the pep to the pp neutrino flux is robustly determined by the SSM calculations, so it can be determined more accurately than the individual fluxes because the ratio depends only weakly on the solar astrophysical inputs. Secondly, a very precise measurement of the 7Be flux, with all the other solar data and assuming energy conservation (luminosity constraint) leads to a very precise determination of the pp and pep flux, at the level of 1% accuracy [48]. On the other hand, to determine 7Be flux experimentally, we have to use the SSM flux to determine the neutrino survival probability, and therefore, the uncertainties in the theoretical estimate [16] limit the precision of the 7Be flux measurement.

The red dash-dotted line in Figure 1 shows the result of the combined analysis of future low energy data, an improved 7Be measurement with 5% precision, and a future pep measurement with 3% precision, added to the current data. The obtained constraint on is at 1 (3) CL. The resultant constraint on from above is much more powerful than the one obtained with spectrum upturn of high-energy 8B neutrinos at 3.

By having solar neutrino spectrum information both at high and low energies it is tempting to ask how tight the constraint becomes if we combine them. The result of this exercise is plotted by the red dashed line in Figure 1 and is also given in Table 2. The resultant constraint on is at 1 (3) CL.

3.4. Day-Night Asymmetry

To have a feeling on to what extent constraint on can be tightened by possible future measurement, we extend the simple-minded model discussed in Section 3.1, but with further simplification of assuming as the best fit. Let us assume that the day-night asymmetry can be determined with % accuracy, an evidence for the day-night asymmetry at CL. Then, the appropriate model is given under the same approximations as in Section 3.1 as . We boldly assume that the day-night asymmetry at CL would be a practical goal in SK. It predicts , which means that can be constrained to the accuracy of 33% uncertainty at 1 CL.

Now, we give the result based on the real simulation of data. The black dotted line in Figure 1 shows the constraint on obtained by future CL measurement of the day-night asymmetry, which is added to the present solar neutrino data. As we see, the day-night asymmetry is very sensitive to the matter potential despite our modest assumption of CL measurement of . (Given the powerfulness of the day-night asymmetry for constraining , it is highly desirable to measure it at higher CL in the future. Of course, it would be a challenging task and probably requires a megaton class water Cherenkov or large volume liquid scintillator detectors with solar neutrino detection capability. They include, e.g., Hyper-Kamiokande [54], UNO [55], or the ones described in [56]). The obtained constraint on is at 1 (3) CL. The obtained upper bound on is actually stronger than the one expected by our simple-minded model . Apart from the shift of the bast fit to a larger value of , the behavior of is more like in the region . It can also be seen in Figure 1 that the upper bound on due to the day-night asymmetry at 3 CL (black dotted line) is stronger than the one from combined analysis of all the expected measurements of the shape of the spectrum (red dashed line) discussed at the end of Section 3.3.

3.5. Global Analysis

We now discuss to what extent the constraint on can become stringent when all the data of various observable are combined. The solid green line in Figure 1 shows the constraint on obtained by the global analysis combining all the datasets considered in our analysis. The obtained sensitivity reads at 1 (3) CL. Therefore, the present and the future solar neutrino data, under the assumptions of the accuracies of measurement stated before, can constrain to 15% (40%) at 1 (3) CL from below and to 20% (60%) at 1 (3) CL from above. If we compare this to the current constraint the improvement of the errors for over the current precision is, very roughly speaking, a factor of 1.5–2 in region , and it is a factor of 2 at . Noticing that the efficiency of adding more data to have tighter constraint at is weakened by shift of the minimum of , improvement of the constraint on is more significant at .

4. Summary

In this paper, we have discussed the question of to what extent tests of the MSW theory can be made stringent by various solar neutrino observables. First, we have updated the constraint on , the ratio of the effective coupling constant of neutrinos to , the Fermi coupling constant with the new data including SNO 8B spectrum and SK day-night asymmetry. Then, we have discussed in detail how and to what extent the solar neutrino observable in the future tightens the constraint on .

The features of the obtained constraints can be summarized as follows.(i) Interpretation of solar neutrino data at high energies by the vacuum oscillation is severely excluded by the SNO and SK experiments, which leads to a strong and robust lower bound on . On the other hand, the day-night asymmetry at 2σ level observed by SK dominates the bound at high side. We find that present data lead to at 1 (3) CL. The Standard Model prediction is outside the 1 CL range but only by tiny amount.(ii) We have explored the improvements that could be achieved by solar neutrinos experiments, ongoing and in construction. We discussed three observables that are sensitive enough to significantly improve the limits on , particularly in the region : (a) upturn of the 8B solar neutrino spectra at low energies at 3 CL, (b) high precision measurement of monoenergetic low-energy solar neutrinos, 7Be (5% precision), and pep (3% precision) neutrinos, and (c) day-night asymmetry of the 8B solar neutrino flux at 3 CL. They lead to the improvement of the bound as follows: (a) at 1 (3) CL,(b) at 1 (3) CL,(c) at 1 (3) CL.It could be expected that future measurement by SNO+ [57, 58] and KamLAND [59] may detect spectrum modulation of B neutrinos at low energies at CL higher than .Finally, by combining all the datasets we have considered we obtain at 1 (3) CL.(iii)As mentioned in Section 1, the issue of effective neutrino matter coupling in a wider context may be better treated in the framework of NSI. If we think about the extended setting together with accelerator and atmospheric neutrino measurement to look for effects of NSI, the off-diagonal elements () can be better constrained by long-baseline experiments. In fact, in a perturbative treatment with small parameter with the assumption , the terms with and are of second order in , while comes only at third order in [60]. The analyses show that the sensitivity to is indeed lower at least by an order of magnitude compared to the ones to or . See the analysis in [61] and the references cited therein. It is also known that can be severely constrained by atmospheric neutrinos [62]. Hence, we feel that the solar neutrinos are still a powerful and complementary probe for in such extended setting.

In conclusion, testing the theory of neutrino propagation in matter deserves further endeavor. The lack of an accurate measurement of the matter potential felt by solar neutrinos reflects the fact that solar neutrino data only do not precisely determine the mass square splitting. The good match of the independently determined mass square splitting by solar neutrino data and by reactor antineutrino data will confirm the Standard Model prediction of the relative index of refraction of electron neutrinos to the other flavor neutrinos. The lack of match of both measurements would point to new physics like the one tested here.

Acknowledgments

We thank the Galileo Galilei Institute for Theoretical Physics and the organizers of the workshop “What is ?” for warm hospitality. C. Peña-Garay is supported in part by the Spanish MICINN Grants FPA-2007-60323 and FPA2011-29678, the Generalitat Valenciana Grant PROMETEO/2009/116, and the ITN INVISIBLES (Marie Curie Actions, PITN-GA-2011-289442). H. Minakata is supported in part by KAKENHI, Grant-in-Aid for Scientific Research No. 23540315, Japan Society for the Promotion of Science.