Abstract
Neutrinos emitted in the carbon, nitrogen, oxygen (CNO) fusion cycle in the Sun are a subdominant, yet crucial component of solar neutrinos whose flux has not been measured yet. The Borexino experiment at the Laboratori Nazionali del Gran Sasso (Italy) has a unique opportunity to detect them directly thanks to the detector’s radiopurity and the precise understanding of the detector backgrounds. We discuss the sensitivity of Borexino to CNO neutrinos, which is based on the strategies we adopted to constrain the rates of the two most relevant background sources, \(pep\) neutrinos from the solar ppchain and \(^{210}\)Bi beta decays originating in the intrinsic contamination of the liquid scintillator with \(^{210}\)Pb. Assuming the CNO flux predicted by the highmetallicity Standard Solar Model and an exposure of 1000 days \(\times \) 71.3 t, Borexino has a median sensitivity to CNO neutrino higher than 3 \(\sigma \). With the same hypothesis the expected experimental uncertainty on the CNO neutrino flux is 23%, provided the uncertainty on the independent estimate of the \(^{210}\text {Bi}\) interaction rate is 1.5 \(\hbox {cpd}/100~\hbox {ton}\) . Finally, we evaluated the expected uncertainty of the C and N abundances and the expected discrimination significance between the high and low metallicity Standard Solar Models (HZ and LZ) with future more precise measurement of the CNO solar neutrino flux.
Introduction
The Sun releases energy mainly through a nuclear fusion process known as the proton–proton chain (pp chain). Another process, called the carbon–nitrogen–oxygen (CNO) cycle (see Fig. 1), is expected to contribute about one percent of the total energy and neutrino production [1, 2]. The CNO cycle emits neutrinos with energies up to around 1.2 MeV and 1.7 MeV for its two main components, see Fig. 2. This process, thought to be the dominant energy production process for stars heavier than 1.3 solar masses [3] as well as solarlike stars in advanced evolutionary stages [4], has many implications for astrophysical problems. For example, the measurement of CNO neutrinos would allow the evaluation of the efficiency of the CNO cycle, helping with the determination of the age of globular clusters [5]. It would also provide a direct reading of the metal abundance in the solar core, which would in turn allow the study of the chemical evolution paradigm assumed by the standard solar model (SSM) [6].
Currently the CNO neutrino rate is only measured to be less than 8.1 counts per day per 100 ton (hereinafter as \(\hbox {cpd}/100~\hbox {ton}\) ) at 95% confidence level by Borexino [7]. This paper presents a detailed study of the sensitivity of the Borexino experiment to CNO neutrinos. The study relies on a precise and independent determination of the two main residual background components: \(^{210}\)Bi contamination of the liquid scintillator, which can be estimated using measurements of \(^{210}\)Po decays (as suggested in [8]), and pep solar neutrinos, which are a part of the pp chain.
After a brief overview of the Borexino detector in Sects. 2, 3 discusses Borexino’s sensitivity to CNO neutrinos. Section 4 present the strategy to constrain the backgrounds and the influence of additional purification of the liquid scintillator. In Sect. 5 we comment on the relevance of the measurement of the flux of CNO neutrinos in the context of solar physics, with an emphasis on the “solar metallicity (or abundance) problem”. This scientific puzzle originated when a redetermination of the surface metallicity of the Sun [9,10,11,12,13] indicated a lower value than previously assumed [14]. However, solar models incorporating these lower metal abundances [10] meet difficulties in reproducing the results from helioseismology, which support models with a higher metal content. Section 5.1 shows that, by following the approach proposed in [15, 16], it is possible to infer the carbon and nitrogen contents of the solar core independently of the assumed opacity of solar plasma by combining a CNO neutrino flux measurement with the very precise measurement of the \(^{8}\hbox {B}\) neutrino flux by the SuperKamiokande collaboration [17] (about 2% precision). Finally, Sect. 5.2 discusses the possibility of using a CNO neutrino measurement to discriminate among SSMs with different hypotheses about the Sun’s surface metallicity by combining a measurement of the CNO neutrino flux with existing \(^7\)Be and \(^8\)B data measured by Borexino [7].
The Borexino experiment
The Borexino detector [20] has been taking data since 2007 in the Laboratori Nazionali del Gran Sasso (Italy) at a depth of 3800 m.w.e. Borexino is an unsegmented calorimeter using about 280 ton of ultrapure organic liquid scintillator. The scintillator consists of PC (pseudocumene, 1,2,4trimethylbenzene) as a solvent with 1.5 g/l of fluor PPO (2,5diphenyloxazole) as a solute. The electron density is \((3.307\pm 0.015)\times 10^{31}~\hbox {e}^\) per 100 ton. This scintillator is contained in a \(125~{\upmu \hbox {m}}\)thick nylon vessel with a nominal radius of 4.25 m. It is surrounded by nonscintillating buffer fluid contained in a 6.85 m radius stainless steel sphere, which supports 2212 8in. ETL 9351 photomultipliers (PMTs). The stainless steel sphere is submerged in ultrapure water serving as an active muon veto and a passive shield against external radiation.
Neutrinos are detected via the elastic scattering reaction on electrons in the liquid scintillator. The electrons deposit their energy in the scintillator, which results in scintillation light. To measure this deposited energy, several energy estimators are developed [18, 19] and they give consistent results in solar neutrino analyses. In this work, we use the variable \(N_h\), the overall number of hits on the PMTs, as a baseline. The interaction vertex is reconstructed via the arrival times of the scintillation photons on the PMTs. As a reference, the detected photoelectron yield is \(\sim \)500 photoelectrons/MeV (normalized to 2000 PMTs). The energy and position reconstruction uncertainties at 1 MeV are \(\sim \,50\,\hbox {keV}\) and \(\sim \,10\,\hbox {cm}\), respectively.
The time distribution of the detected hits on the PMTs allows pulseshape discrimination (PSD). With PSD, \(\alpha \) particles, a product of \(^{210}\text {Po}\) decay in the scintillator, can be discriminated from \(\beta \) particles on an eventbyevent basis [21, 22]. It can also be used to constrain the ratio of the total counts of \(\beta ^+\) events, such as \({}^{11}\hbox {C}\) decays, to \(\beta ^\) events in specific energy ranges [18, 19].
Borexino sensitivity to CNO neutrinos
In the Borexino PhaseII analysis [7], the signal and background rates were determined using the multivariate fit method. In detail, as discussed in [18, 19], we started from a background model consisting of a list of unstable isotopes contaminating the scintillator (namely \({}^{14}\hbox {C}\), \(^{85}\)Kr, \(^{210}\)Po, \(^{210}\)Bi), isotopes of cosmogenic origin (\(^{11}\)C), and \(\gamma \)rays emitted in the regions outside of the sensitive volume and referred to as external backgrounds (originating from decays from \(^{208}\)Tl, \(^{214}\)Bi, and \(^{40}\)K).
We then performed a multivariate fit procedure of events selected in an optimized wallless fiducial volume (FV) of about 70 ton in the central part of the detector that included two observables: the energy spectra of the events in the region from 0.19 to 2.93 MeV and their reconstructed radial coordinate distribution [19]. This fitting procedure made it possible to disentangle the rates of events induced by neutrinos from those induced by backgrounds. At the same time, the fair goodnessoffit (\(p\) value of 0.7) verified the consistency of the background model.
In this study, we make assumptions of background rates and signal rates based on the multivariate fit analyses of PhaseIII data, which starts from 2016 June. The expected energy distribution of the signals and of the most relevant background components can be seen in Fig. 3, while the expected and measured interaction rates of solar neutrinos (PhaseII results) are listed in Table 1, and the expected signal and background rates considered in the present study are listed in Table 2.
There is an intrinsic difficulty in the measurement of solar neutrinos due to the similarity of the energy spectra of \(^{210}\text {Bi}\) electrons and those of electron recoils induced by pep and CNO neutrinos. This feature, which was discussed in [18, 19], is clearly visible in Fig. 3. The partial degeneracy of these spectral shapes induces significant correlations among these three components. Only by constraining the rates of two out of three components can one measure the rate of the third precisely.
However, the measurement of \({R_\mathrm{CNO}}\) (i.e., the interaction rate of CNO solar neutrinos) is now possible in Borexino PhaseIII as we have identified a strategy (discussed in Sect. 4) to measure or constrain the two main components, pep neutrinos and \(^{210}\text {Bi}\) in the scintillator. The thermal stabilization of the detector, achieved during the PhaseIII period, is fundamental to measure the rate of \(^{210}\text {Bi}\) . The low internal background levels, great depth of LNGS, passive detector shielding, and active removal (based on the threefoldcoincidence method [18], i.e. TFC cut) of the cosmogenic \({}^{11}\hbox {C}\) background are also crucial.
Further, we observed an energy region from about 0.8 MeV to 1 MeV, hereinafter referred to as region of interest (ROI), that is dominated by \(R_\mathrm{Bi}\), \(R_{pep}\), and \({R_\mathrm{CNO}}\) (see Fig. 4). More precisely, in this ROI the background contributions are of the same order as the statistical fluctuations of the total expected CNO rate. Thus, it is also possible to extract the sensitivity of Borexino to CNO neutrinos through a simple procedure consisting only of determining the background rates and subsequently counting all the events in this energy region. Section 3.1 first discusses this counting analysis. Then, Sect. 3.2 presents a sensitivity study performed with the multivariate spectral fit. Finally, Sect. 3.3 and 3.4 compare the resulting expected sensitivities and discovery significance from both methods.
In the following sections, our best estimates of the conditions of the Borexino PhaseIII are applied. The \(^{210}\text {Bi}\) rate \(R_\text {Bi}\) is assumed to be 10 \(\hbox {cpd}/100~\hbox {ton}\) . The exposure of \({1000}~\hbox {day}\times 70~\hbox {ton}\) considers the same fiducial volume as the one applied in the previous solar analysis [7] and amounts to \(192\,\text {year}\,\times \,\text {ton}\) in total. In the counting analysis, only events after the TFC cut are used, and the exposure loss is 36%. In the multivariate fit analysis, all events are used. In order to facilitate the comparison of the performances of the two methods, we use the same values of the constraint of the \(pep\) neutrino interaction rates as well as the \(^{210}\text {Bi}\) contamination of the scintillator. The pep constraint (Sect. 4.1) is based on the predictions of the Standard Solar Model, while \(^{210}\text {Bi}\) constraint (Sect. 4.2) is centered in the range considered achievable in Borexino.
We note, that with a high exposure, the intrinsic small difference between the spectral shapes of CNO and \(^{210}\text {Bi}\) makes it possible, in principle, to measure the interaction rate of CNO neutrinos without any constraint on \(R_\mathrm{Bi}\). With a spectral analysis similar to that discussed in this paper, we estimate that the statistical significance of about \(3\sigma \) can be reached for an exposure of \(450\,\text {year}\,\times \,\text {ton}\) assuming CNO neutrinos from the HZ SSM. The statistical sensitivity of such analysis improves as the inverse square root of exposure. However, the systematic uncertainties due to the detector response and energy scale modeling must be treated carefully.
Borexino as a counting experiment
This approach consists of counting the total number of events in the ROI to obtain the interaction rate of CNO neutrinos as the difference between the total number of events detected in ROI and that of the backgrounds, as evaluated in an independent way.
Our key assumptions lie in the availability of a realistic background model and the ability to measure the rate of the residual background components independently (or to label them as negligible). Another assumption is the availability of a sufficiently accurate model of the detector response function, needed to determine the fraction of each spectral component in the ROI.
As shown in Figs. 3 and 4, there are six contributions to the event rate in the ROI: CNO neutrinos, \(^{210}\text {Bi}\) electrons, pep neutrinos, residual \(^{11}\text {C}\) positrons, \(^{7}\)Be neutrinos, and external gamma backgrounds. The impact of the relevant parameters like the exposure, the choice of the ROI, and the precision of the estimation of pep neutrino and \(^{210}\text {Bi}\) rates can be understood by writing \(R_\text {CNO}\) and its statistical uncertainty \( \sigma ^{count}_\text {CNO}\) as:
where \(N^\text {ROI}_x\) (\(x=\)tot(al) or o(thers)) is the number of events in ROI of \(x\), \(\mathcal {E}\) is the exposure, \(\varepsilon \) is the fraction of events for each component falling inside the ROI, R and \(\sigma \) are the rate and the uncertainty, respectively, of components in the full energy range, and \(r_\text {o}\) is the total event rate of all other components in the ROI. The contributions to \(r_\text {o}\) are either negligible and/or independently determined. The \(^7\hbox {Be}\) solar neutrino rate is constrained to PhaseII results [7]. \(^{85}\hbox {Kr}\) is determined with the fast coincidence tagging of its minor branch [7]. \({}^{11}\hbox {C}\), external \(^{214}\hbox {Bi}\), \(^{208}\hbox {Tl}\), and \(^{40}\hbox {K}\) are determined together with events outside the ROI via the spectral fitting method. \(^{210}\text {Po}\) is negligible in the ROI. Variables with a tilde mark the values from independent determinations and the subscripts indicate the spectral component. The symbol \(\oplus \) is defined as \(a\oplus b = \sqrt{a^2+b^2}\).
The values of the assumed rates, uncertainties, and efficiencies are listed in Table 2. The ratios of the four terms in Eq. (1), from left to right, are \(10{:}3{:}3{:}2\). When \(\tilde{\sigma }_\text {Bi}\) is 1.5 \(\hbox {cpd}/100~\hbox {ton}\) and \(\tilde{\sigma }_{pep}\) is \(1.4\%\cdot \tilde{R}_{pep}\), or 0.04 \(\hbox {cpd}/100~\hbox {ton}\) , the ratios of the four terms in Eq. (2), from left to right, are \(9{:}11{:}1{:}3\). The largest contribution to \( \sigma ^{count}_\text {CNO}\) in Eq. (2) is the \(^{210}\text {Bi}\) term unless the uncertainty on the \(^{210}\text {Bi}\) constraint is lower than 1.3 \(\hbox {cpd}/100~\hbox {ton}\) when \(\sigma _\text {tot}=\varepsilon _\text {Bi} \tilde{R}_\text {Bi}\), in which case \( \sigma ^{count}_\text {CNO}\) is limited by a statistical error. The resulting uncertainty on the CNO neutrino interaction rate as a function \(\tilde{\sigma }_\text {Bi}\) and \(\tilde{\sigma }_{pep}\) can be seen in Fig. 5.
Multivariate fit analysis
We performed a study of the sensitivity of Borexino to CNO neutrinos by simulating tens of thousands of pseudoexperiments and fitting the pseudodata with the full multivariate procedure as the one used for real data [7]. This approach requires maximizing a two dimensional likelihood function built using two observables (energy estimator and radial position). We assumed, as throughout Sect. 3, that \({R_{pep}}\) and \(R_\mathrm{Bi}\) can be independently estimated. The constraints on \({R_{pep}}\) and \(R_\mathrm{Bi}\) were implemented by introducing two Gaussian pull terms in the likelihood function. The centroid of the pull terms are varied according to the uncertainty of the pull. The uncertainties of the two constraints are the most important parameters when determining the sensitivity.
Each simulated dataset was fitted assuming the same detector response function as that used in the simulation. The distribution of the bestfit results for the CNO interaction rate was used to assess the expected statistical uncertainty and discovery significance. The results were obtained using the multivariate fitting tools GooStats [25] and mstats [26].
Similarly to the data analysis of [7], the study was conducted using both Monte Carlo and analytical methods that are complementary to each other. In the Monte Carlo approach the detector response was modeled using an ab initio simulation software [24] that takes into account realistic and microscopic descriptions of the energy deposition, scintillation light generation and propagation, electronics response, and energy and position reconstruction. In the analytical approach, the detector response was modeled using analytical functions [19]. The analytical procedure provides an independent way of evaluating the impact of the systematic uncertainty related to detector response parameters (e.g., the light yield). The two approaches gave consistent results for the expected uncertainty and discovery significance.
Uncertainty of the CNO rate
In this section the uncertainty of the CNO rate resulting from both the counting and multivariate fit analyses is discussed. In the framework of the multivariate analysis, the distribution of the values of \( R_\mathrm{{CNO}} \) resulting from fits performed with fixed external constraints is built. The width of this distribution is used to estimate \(\sigma ^{fit}_\text {CNO}\). Figure 6 shows \(\sigma ^{fit}_\text{ CNO }\) as a function of the \(^{210}\text {Bi}\) constraint uncertainty for two different values of \(\tilde{\sigma }_{pep}\) as well as the uncertainty \(\sigma ^{count}_\text {CNO}\) resulting from the counting analysis presented in Sect. 3.1.
In the counting analysis, the rates of \(^{210}\text {Bi}\) and pep are fixed by the external constraints. Thus, according to Eq. (1), a bias on the constraints will transfer linearly to the reconstructed CNO rate. The situation becomes slightly more complicated when a full spectral analysis, combined with constraints on the background, is performed. Indeed, in the case of a bias in the external constraints, the impact on the value of the extracted CNO rate is mitigated by the tension between the shape information and the biased values.
The systematic uncertainty of the CNO rate is dominated by the bias of the \(^{210}\text {Bi}\) accuracy (see Sect. 4.2), followed by the accuracy on the light yield.
Expected discovery significance to CNO neutrinos
A frequentist hypothesis test was performed to assess the expected discovery significance to CNO neutrinos. In the search for the CNO signal, two hypotheses are considered: the null hypothesis, \(H_0\), where no CNO signal is present; and the alternative hypothesis, \(H_1\), that includes the presence of a CNO signal in addition to the background. We used the profile likelihood ratio \(q_0\) as the test statistic [27]:
where \(L(H_0)\) and \(L(H_1)\) are the the maximum values of the likelihood against \(H_0\) and \(H_1\). The significance of the signal of a measurement is quantified by evaluating the compatibility of the observed data with the null hypothesis \(H_0\) and represented by the \(p\) value. For a specific measurement with a test statistic value \(q_{0,\text {obs}}\), its \(p\) value is the probability that the test statistic is above \(q_{0,\text {obs}}\) under \(H_0\). In quantifying the expected sensitivity of \(H_1\) to discovery of a signal, the median discovery significance is used. The median discovery significance of \(H_1\) is the \(p\) value of the median value of the test statistic under \(H_1\). This requires the probability density functions (hereinafter as PDFs) of \(q_0\) for both the null hypothesis \(f(q_0H_0)\) and the alternative hypotheses \(f(q_0H_1=\text {HZ or LZ})\).
The PDFs of the test statistics \(q_0\) were obtained by analyzing the two sets of the simulated datasets (with and without the CNO injected), with both the simplified counting analysis and full multivariate fit. We generate 20,000 pseudoexperiments to construct the PDFs. When considering the null hypothesis, \(q_0\) is distributed according to \(1/2\delta (q_0)+\chi ^2_1(q_0)\), where \(\delta (q_0)\) is the DiracDelta function, and \(\chi ^2_1\) is the PDF of the chisquare distribution with one degree of freedom. For \(H_1\) and in the case of uncertainties of 0.04 \(\hbox {cpd}/100~\hbox {ton}\) and 1.5 \(\hbox {cpd}/100~\hbox {ton}\) on pep neutrinos and \(^{210}\text {Bi}\) , respectively, the median value of the test statistic is 8.0 for HZ SSMs and 4.1 for LZ SSMs.
The expected discovery significance depends on the strength of the signal and on the precision of the external constraints: a higher rate of CNO neutrinos, as well as stronger constraints, results in a higher sensitivity. For the CNO rate predicted by the HZ SSM and in the case of uncertainties of 0.04 \(\hbox {cpd}/100~\hbox {ton}\) and 1.5 \(\hbox {cpd}/100~\hbox {ton}\) on pep neutrinos and \(^{210}\text {Bi}\) , respectively, the median discovery significance to CNO neutrinos is \(4.5\sigma \). In case of the LZ SSM, it is \(3.2\sigma \). Results of the expected discovery significance under other conditions are shown in Fig. 7.
The difference between the median p value obtained from the simple counting analysis and the one from the multivariate analysis gets smaller as the precision of the constraints increase. This effect is coherent with the results on the statistical uncertainty of CNO presented in Sect. 3.3, indicating that the impact of the spectral shape information is larger when the constraints are relatively weak (high uncertainty), while it becomes negligible when the background constraints get more stringent (low uncertainty).
Impact of an upper limit on \(^{210}\text {Bi}\) rate
The method discussed in Sect. 4.2, used to obtain an independent measurement of the \(^{210}\text {Bi}\) rate, can only provide an upper limit on the \(^{210}\text {Bi}\) background under less stringent assumptions, that an additional contribution from migrated (diffusion and/or convection) \(^{210}\text {Po}\) (see Sect. 4.2) is present. The presence of migrated \(^{210}\text {Po}\) leads to a positive bias of the estimation of the \(^{210}\text {Bi}\) rate.
In this case, the constraint is implemented in the fitting procedure as a onesided Gaussian penalty:
where \(R_\text {mig Po}\) is the rate of the migrated \(^{210}\text {Po}\) as discussed in Sect. 4.2, and \(\mu \) and \(\sigma \) are the centroid and uncertainty of the \(^{210}\text {Bi}\) pull term. In so doing, the upper limit on the \(^{210}\text {Bi}\) rate is equivalent to a lower limit on the CNO rate. Therefore, the expected discovery significance for CNO is the same as would be obtained using a twosided Gaussian constraint on the \(^{210}\text {Bi}\) rate. This was confirmed by performing a hypothesis test analogous to that described in Sect. 3.4. The results for varying widths (\(\tilde{\sigma }_{\text {Bi}}\)) on the \(^{210}\text {Bi}\) penalty are shown in Fig. 8. The expected discovery significance obtained from the counting analysis and from applying a symmetric constraint are also displayed for comparison. As one can see, allowing an additional contribution from the migrated \(^{210}\text {Po}\) barely decreases the expected discovery significance for CNO.
Therefore, an upper limit on the \(^{210}\text {Bi}\) rate may be sufficient to claim detection of CNO neutrinos, but it would not allow a precise measurement of the CNO interaction rate required to solve the solar metallicity problem. Since there is no lower limit on the \(^{210}\text {Bi}\) rate, the upper limit of the CNO neutrino rate will not be stringent and the central CNO value might have a negative bias. Figure 9 shows the distribution of the best fit estimates, obtained from the fit of simulated datasets when applying a constraint to the pep neutrino rate and leaving the \(^{210}\text {Bi}\) rate free (Fig. 9a), constraining the \(^{210}\text {Bi}\) rate with a symmetric Gaussian penalty (Fig. 9b), and applying the upper limit described in this section (Fig. 9c).
The expected confidence intervals in Fig. 9c are asymmetric: the upper limit of the confidence interval is similar to the one obtained when leaving the \(^{210}\text {Bi}\) rate free, while its lower limit resembles the one resulting when the \(^{210}\text {Bi}\) rate is constrained.
Strategy for establishing background assessment
This section discusses the strategy for setting constraints on the pep neutrinos and \(^{210}\text {Bi}\) rates independently of both the counting analysis and the multivariate fit procedure.
pep constraint
The electron capture reaction \(p+e^{}+p \rightarrow d +\nu _{e}\), which generates pep neutrinos in the Sun, is linked to the \(\beta ^+\) decay process \(p+p\rightarrow d + e^{+} + \nu _{e}\), which is responsible for the pp neutrino production, by wellknown nuclear physics. Since the two processes depend on the same allowed nuclear matrix element, the ratio between their rates is determined by the available reaction phase spaces and by the electron density \(n_{e}\) of the solar plasma only. This ratio was calculated in [28, 29], and the effect of radiative corrections subsequently discussed in [30] (see, for example, [31] for a review). It can be determined with \(\sim 1\%\) precision [31] for the conditions of the solar interior and is mildly dependent on the properties of the solar plasma, roughly proportional to \(T_c^{1/2} n_{e}\) (where \(T_c\) is the temperature of the core of the Sun). As a consequence, the ratio \( \varPhi _\mathrm{pep}/\varPhi _\mathrm{pp} \) between the pep and pp neutrino fluxes is a robust prediction of SSMs, and it can be used to improve the sensitivity to the CNO neutrino signal. With this approach, Borexino’s direct observation of the pp neutrino flux [7] can be translated into a \(\sim 10\%\) determination of the pep neutrino component, motivating the analysis performed in Sect. 3, which assumes \(\tilde{\sigma }_{{ pep}\,} = 0.28\) \(\hbox {cpd}/100~\hbox {ton}\) .
The precision of the pp and pep neutrino flux determination can be further improved by performing a global analysis [32, 33] on all neutrino experiment results applying the socalled solar luminosity constraint [34,35,36]. In this way, it was shown in [33] that the pep neutrino flux is constrained with \(\sim 1\%\) precision of by solar neutrino data. By taking into account an additional \(\sim 1\%\) uncertainty in neutrino survival probability considering the matter effects in the Sun and the uncertainty of oscillation parameters, the event rate \(R_{pep}\) can be constrained with a precision of \(\sim 1.4\%\), motivating the assumption \(\tilde{\sigma }_\mathrm{pep} = 0.04\) \(\hbox {cpd}/100~\hbox {ton}\) considered in Sect. 3. Because the contribution of CNO neutrinos to the solar luminosity is around 1%, neglecting the dependence of this pep constraint on the assumed \(R_\mathrm{CNO}\) has almost no effects on the significance to CNO neutrinos, as shown in Fig. 5.
\(^{210}\text {Bi}\) constraint
\(^{210}\text {Bi}\) is a \(\beta \)emitting daughter of \(^{210}\text {Pb}\) with a mean lifetime of 7.23 days and a \(Q\) value of 1160 keV. Since the lifetime of \(^{210}\text {Bi}\) is small, it must be supported by its longlasting parent nucleus \(^{210}\text {Pb}\) (\(\tau = 32.2\) years) to maintain a constant decay rate. \(^{210}\text {Pb}\) is a part of the \(^{238}\hbox {U}\) chain, and due to its low \(Q\) value of 64 keV it does not represent a background in this analysis. The \(^{210}\)Pb and \(^{210}\)Bi are part of the of the \(^{238}\hbox {U}\) chain. The concentration of \(^{238}\hbox {U}\) is estimated to be less than \(9.4\times 10^{20}\) g/g (95% C.L.) by the absence of \(^{226}\)Ra and \(^{222}\)Rn, easily detectable through the fast \(^{214}\)Bi\(^{214}\)Po coincidence following their decays in the part of the chain above \(^{210}\)Pb. We note, that the secular equilibrium of the chain is broken at the level of \(^{210}\)Pb. This means, that while the background from the \(^{238}\)U chain down to \(^{210}\)Pb is negligible, the background from the chain below \(^{210}\)Pb is important. The purification campaign in 2011 [37], that marked the start of the Borexino PhaseII, has significantly reduced in short order the initial \(^{210}\text {Pb}\) contamination, leading to a residual rate in the range of a few tens of \(\hbox {cpd}/100~\hbox {ton}\) .
The strategy to independently determine the \(^{210}\text {Bi}\) rate using the decay rate of its decay product, \(^{210}\text {Po}\) , was suggested in [8]. The isotope \(^{210}\text {Po}\) , with a mean lifetime of 199.6 days, ends the chain by decaying into stable \(^{206}\)Pb thereby emitting a 5305 keV \(\alpha \) particle (visible energy around 400 keV electron equivalent):
When the above chain is in equilibrium, the \(^{210}\text {Bi}\) rate is equal to that of \(^{210}\text {Po}\) . The \(^{210}\text {Po}\) rate can be precisely determined in Borexino because \(^{210}\text {Po}\) emits \(\alpha \) particles that can be identified eventbyevent via pulseshape discrimination. Using a pulseshape analysis based on a multilayer perceptron discriminator [21, 38], \(^{210}\text {Po}\) events in Borexino are identified with an efficiency very close to 1.
Unfortunately, the measured \(^{210}\text {Po}\) rate is not only due to \(^{210}\text {Bi}\) decays, but consists of three components:

1.
Unsupported \(^{\mathbf{210}}{} \mathbf{Po}\) The residual \(^{210}\text {Po}\) left over by the water extraction phase of the scintillator purification campaign and not linked to local \(^{210}\text {Bi}\) .

2.
Migrated \(^{\mathbf{210}}{} \mathbf{Po}\) Originally produced from the \(^{210}\text {Pb}\) lead decays on the inner surface of the nylon vessel holding the liquid scintillator and brought into the fiducial volume by convective and diffusive motions of the liquid scintillator.

3.
Supported \(^{\mathbf{210}}{} \mathbf{Po}\) In secular equilibrium with local \(^{210}\text {Bi}\) present in the liquid scintillator.
The \(^{210}\text {Bi}\) rate is equal to the supported \(^{210}\text {Po}\) rate and thus less than the total \(^{210}\text {Po}\) rate in presence of the unsupported and migration terms. The rate of unsupported \(^{210}\text {Po}\) , following the law of radioactive decay, gets asymptotically closer to zero over time with a mean lifetime of 199.6 days. The migrated \(^{210}\text {Po}\) is the most intricate contribution to handle. In Borexino, diffusion of \(^{210}\text {Po}\) into the fiducial volume is a completely negligible process, because on average \(^{210}\text {Po}\) can only travel 15 cm via pure diffusion before it decays, while presence of migrated \(^{210}\text {Po}\) is found in the center region where \(^{210}\text {Po}\) need to travel by around 4.5 m. Therefore the movement of \(^{210}\text {Po}\) is mainly driven by convection motion in the liquid scintillator. Because this convection motion is generated by the inhomogeneity and instability of the detector temperature, the migrated \(^{210}\text {Po}\) is time and spacedependent and hard to model.
These considerations are summarized through the following equation:
where \(R_\text {Po}^\text {tot}\) is the observed \(^{210}\text {Po}\) rate, \(R_\text {Po}^\text {u}\) is the initial rate of unsupported \(^{210}\text {Po}\) decay, \(R_\text {Po}^\text {m}\) is the migrated \(^{210}\text {Po}\) decay rate, and \(R_\text {Po}^\text {s}\) is the supported \(^{210}\text {Po}\) rate and is equal to the \(^{210}\text {Bi}\) rate.
The convective motion of the liquid scintillator was suppressed in a major effort over several years since 2015 [39]. This included the installation of a passive thermal insulation system around the water tank and an active temperature control system to mitigate seasonal modulations. Our observations showed that the migrated \(^{210}\text {Po}\) was significantly reduced in the Borexino PhaseIII which starts from 2016 June.
In PhaseIII, the rate of unsupported \(^{210}\text {Po}\) has decayed and is less than 0.1 \(\hbox {cpd}/100~\hbox {ton}\) . It is also possible to get a more precise estimate of the rate of supported \(^{210}\text {Po}\) by identifying the region inside the fiducial volume, that is least affected by convective motions and thus with minimal migrated \(^{210}\text {Po}\) .
The residual migrated \(^{210}\text {Po}\) rate in this region can then be treated as a source of systematic uncertainty. Neglecting the migrated \(^{210}\text {Po}\) rate, the statistical uncertainty of the supported term depends on the magnitude of the initial outofequilibrium contamination and on the exposure. As an example, with an unsupported rate of 50 \(\hbox {cpd}/100~\hbox {ton}\) and a supported rate of 10 \(\hbox {cpd}/100~\hbox {ton}\) , an precision of the order of \(10\%\) on the supported rate can be achieved in about 8 months of datataking using the same fiducial volume of 71.3 ton as adopted in [7].
Furthermore, in order to extrapolate the value of \(^{210}\text {Bi}\) obtained in the region with minimum \(^{210}\text {Po}\) migration rate to the entire fiducial volume, two hypotheses must be verified: the \(^{210}\text {Bi}\) spatial distribution must be uniform within the entire fiducial volume and its rate must follow the slow decay rate of \(^{210}\text {Pb}\) , since there are no sources of either \(^{210}\text {Bi}\) or \(^{210}\text {Pb}\) in the liquid scintillator. The residual nonuniformity and instability of the \(^{210}\text {Bi}\) rate can also be treated as systematic uncertainties of the \(^{210}\text {Bi}\) constraint.
Impact of further scintillator purification
The purification campaign of the Borexino scintillator performed in the year 2011 [37], with 6 cycles of closedloop water extraction, significantly reduced the concentrations of several radioactive contaminants of the scintillator. The \(^{210}\text {Bi}\) rate was reduced by a factor \(\sim \) 2.3 [7, 19]. In principle, another cycle of purification could decrease the amount of \(^{210}\text {Bi}\) even further without introducing more \(^{210}\text {Po}\) in the liquid scintillator. However, in order to keep the amount of migrated \(^{210}\text {Po}\) small, thus for the sensitivity to CNO neutrinos to benefit from further purification, the thermal stability of the detector must be maintained.
In order to quantify the effect of a reduced rate of \(^{210}\text {Bi}\) , we rewrite Eq. 2 as
where \(a_1=1.07\) \(\hbox {cpd}/100~\hbox {ton}\) is the total rate of components other than \(^{210}\text {Bi}\) in the ROI, \(a_2=0.06\) \(\hbox {cpd}/100~\hbox {ton}\) is the total uncertainty of the last three terms in Eq. 2 and is dominated by \(\tilde{\sigma }_\text {Bi}\), and \(\mathcal {E}\) is the exposure in day\(\times \) 100 ton. From Eq. 7 we see that the importance of \(\mathcal {E}\) and \(R_\mathrm{Bi}\) is greater when \(a_2\) is smaller than the first term. Thus, in these conditions, reducing the rate of \(^{210}\text {Bi}\) is equivalent to increasing the exposure. This can also be concluded from Fig. 10.
CNO neutrinos in the standard solar model
Despite the small contribution of the CNO cycle to the Sun’s luminosity, a measurement of the CNO neutrino flux would be extremely valuable to expand our knowledge of the Sun.
Experimental measurements of solar neutrino fluxes are a fundamental test of the SSM, the benchmark of every stellar evolution model. Possible disagreements between predicted and measured neutrino fluxes may highlight issues in the model assumptions. Additionally, measurements of solar neutrino fluxes can be used to infer some properties of the Sun that are used as inputs for the SSM.
This is particularly interesting in the context of the solar metallicity problem. As mentioned in the introduction, SSMs implementing the most recent determinations of the solar surface composition (e.g., the low metallicity (LZ) admixture described in [10]), fail to reproduce helioseismic observations, which agree better with models that assume a higher metal content, as prescribed by the older determinations from [14]. So far, all attempts to understand the origin of this discrepancy remain inconclusive. New experimental results on solar neutrinos could provide key pieces of information to solve this problem. This applies in particular to the flux of CNO neutrinos, whose value differs by almost 30% between the predictions of the SSM HZ and LZ models (see Table 1).
However, the dependence of predicted solar neutrino fluxes on the metallicity and other SSM parameters is nontrivial. There is indeed a strong degeneracy between the impact of metallicity and radiative opacity through their effect on the core temperature of the Sun \(T_c\) (\(T_c\) is determined by the opacity profile, which in turn depends on the chemical composition of the solar interior). Ultimately, \(T_c\) acts as a proxy for the influence of those parameters (i.e., they impact neutrino fluxes via their effect on \(T_c\)). This degeneracy makes it difficult to disentangle metallicity from opacity with solar neutrino measurements, in particular after recent laboratory measurements [40] and theoretical advances [41,42,43,44,45] suggested that the uncertainty of the radiative opacity might be severely underestimated. Deviations of the radiative opacity from its nominal values can mimic the effect of a higher metal content in the Sun’s core. As it was discussed in, for example, [46,47,48], the agreement of SSM calculations implementing the recent LZ composition and helioseismic data could be restored by considering suitable modifications to the radiative opacity. Similarly, neutrino fluxes expected by the HZ SSM can be obtained by assuming the LZ composition and a larger opacity of the solar plasma.
However, the case of the CNO neutrinos is somewhat different. In fact, the dependence of the CNO neutrino flux on metallicity is twofold: as for the other ppchain neutrinos the CNO flux is indirectly linked to metallicity via \(T_c\), but it also depends directly on the amount of carbon and nitrogen in the core of the Sun. Considering the metallicity problem, other tests of the SSM that consider measurements of \(^{8}\hbox {B}\) and \(^{7}\hbox {Be}\) neutrino fluxes include a larger amount of information but are greatly affected by uncertainties on the radiative opacity.
Section 5.1 shows that by exploiting this fact, a measurement of the CNO neutrino flux with Borexino can be used to infer the content of carbon and nitrogen in the Sun’s core almost independently of the effect of radiative opacity, thus directly probing the solar metallicity. Section 5.2 discusses the impact of a future measurement of CNO neutrinos on the discrimination power of a hypothesis test between the HZ and LZ SSM.
CNO neutrinos as CN abundance messengers
As explained in [16], the temperature dependence for both the \(^{8}\hbox {B}\) and CNO neutrino fluxes can be described by a powerlaw \(\varPhi _i \propto T_c^{\gamma _i} (i \in \{^{8}\text {B}; \text {CNO}\})\). Therefore, one can build a weighted ratio \(\displaystyle \varPhi _{\text {CNO}}/(\varPhi _{ ^{8}\text {B}})^k\), with \(\displaystyle k = \gamma _{\text {CNO}}/\gamma _{^{8}\text {B}}\), that is very nearly (considering different uncertainties) independent of \(T_c\). This breaks the opacitycomposition degeneracy, making it possible to infer the abundance of carbon and nitrogen in the Sun in a way that is almost independent from variations in the SSM parameters that affect the temperature profile.
To elaborate, we followed the approach of [15, 16], which adopts the practice of factorizing the dependence of the neutrino flux \(\varPhi _i\) on the model’s input parameters \(\{\beta _j\}\):
where the \(x_j\) terms are the model inputs normalised to their respective nominal SSM value (i.e., \(\beta _j/\beta _j^\text {SSM}\)). These parameters are typically grouped in three categories: (i) the solar parameters, related to the Sun’s astrophysical (age \(A_{\odot }\), luminosity \(L_{\odot }\)) and nonnuclear physical properties (diffusion D, radiative opacity \(\kappa \)), (ii) the metallicity parameters (i.e., the abundances of C, N, O, F, Ne, Mg, Si, S, Ar, and Fe), and (iii) the nuclear cross sections of the relevant processes, described by the astrophysical Sfactors. The \(\alpha ( i,j )\) coefficients in the exponents are normalised logarithmic partial derivatives of the fluxes with respect to the input parameters \(\beta _j\)
which are typically provided along with SSM predictions.
Using Eq. 8 to express solar neutrino fluxes, the aforementioned weighted ratio between the \({}^{15}\hbox {O}\) and \(^{8}\hbox {B}\) fluxes reads:
To reduce the dependence on \(T_c\), the weight k is chosen such that it minimizes the uncertainty of the quantity defined in Eq. 10 due to the socalled environmental parameters (i.e., those SSM inputs which affect the Sun’s temperature profile the most [15, 16]). These parameters include the solar and metallicity parameters discussed above, but without the C and N abundances.
When applying this strategy to Borexino we refer, for simplicity, to the counting analysis described in Sect. 3.1. In this narrow energy window, neglecting insignificant \(^{17}\hbox {F}\) neutrinos, CNO neutrino events originate from contributions of \({}^{15}\hbox {O}\) and \({}^{13}\hbox {N}\) neutrinos. We thus define the flux of CNO neutrinos measured by Borexino \(\varPhi _{\mathrm{CNO}}^\text {BX}\) as a combination of \(\varPhi _{{}^{15}{\mathrm{O}}}\) and \(\varPhi _{{}^{13}{\mathrm{N}}}\). This can be written as:
where \(\xi \) is the ratio between the event rate of \({}^{15}\hbox {O}\) neutrinos and all CNO neutrinos in the ROI, as expected from the SSM:
Using the logarithmic derivatives of \({}^{15}\hbox {O}\) and \({}^{13}\hbox {N}\) neutrinos from the B16SSM [6], we can describe the dependence of the measured CNO neutrino flux on the inputs of the SSM as
Having computed those derivatives, we can express Eq. 10 using a CNO measurement by Borexino and the \(^{8}\hbox {B}\) neutrino flux measured by SuperKamiokande [17]. In the case of Borexino, the value of k, obtained from a minimization of the uncertainty from the environmental parameters, is found to be 0.716. To explicitly show the dependence on all SSM inputs:
where the contributions of the environmental parameters^{Footnote 1} are grouped in the second and third rows on the righthand side, while the contributions of nuclear cross sections are in the fourth row. Following [15, 16], the exponents in Eq. 14 were obtained using the values of the logarithmic derivatives of the solar neutrino fluxes with respect to the SSM inputs used by the HZ SSM [6].
Assuming, for simplicity, that the C and N abundances are modified by the same factor (i.e., \(x_{\mathrm{C}} = x_{\mathrm{N}} \equiv (N_{\mathrm{C}}+N_{\mathrm{N}})/(N_{\mathrm{C}}^\text {SSM} + N_{\mathrm{N}}^\text {SSM})\), where \(N_{\mathrm{C}}\) and \(N_{\mathrm{N}}\) indicate the number density of C and N with respect to hydrogen) and noticing in Eq. 14 that the sum of the exponents of \(x_{\mathrm{C}}\) and \(x_{\mathrm{N}}\) is \(0.814+0.191\simeq 1\), we can invert the formula, obtaining:
where \(R_{\mathrm{CNO}}^\text {BX}\) is the rate measured by Borexino and proportional to the flux \(\varPhi _{\mathrm{CNO}}^\text {BX}\). The quoted uncertainties were obtained by propagating the uncertainties of the environmental, nuclear, and diffusion parameters used in the SSM. The current uncertainty on the measurement of the boron neutrino flux is about 2% [33]. The dominant contributions to the uncertainty budget besides the CNO neutrino rate are from nuclear reactions, with the largest coming from \(S_{114}\) (\(7.3\%\)), \(S_{34}\) (\(3.4\%\)), and \(S_{17}\) (\(3.5\%\)).
The low level of the environmental contribution to the total error budget (which takes into account all SSM parameters except nuclear reactions and diffusion coefficients) demonstrates that a future CN flux measurement will be converted into an almost direct determination of the C+N content of the solar core, as discussed in [15, 16].
The impact of the experimental and model uncertainties in the determination of the C+N abundance inferred from a CNO neutrino measurement, assuming the rate predicted by the HZ SSM, is represented in Fig. 11. While a 1.5 \(\hbox {cpd}/100~\hbox {ton}\) precision is translated, according to Eq. 15, to a precision that is about three to four times larger than the precision of the GS98 and AGSS09 catalogues, a future measurement of \(R_{\mathrm{CNO}}\) with a precision \(\simeq \) 0.5 \(\hbox {cpd}/100~\hbox {ton}\) (achievable by next generation experiments [49,50,51]) will be able to constrain the C+N fraction of the Sun with an uncertainty \(\simeq 15\%\), which is comparable to the precision of the current estimations obtained from measurements of the photosphere.
High versus low metallicity standard solar model
If one considers the current treatment of the opacity and of its corresponding uncertainties in the SSM as reliable, then the measurements of neutrino fluxes due to the pp chain can also be used to test models of the composition of the Sun.
In the results of Borexino PhaseII [7], the measurement of the fluxes of \(^{7}\hbox {Be}\) and \(^{8}\hbox {B}\) neutrinos were used to this end because of their strong dependence on the Sun’s core temperature, which is affected by the metal content as described above. Including CNO neutrinos in a similar test is particularly interesting since their flux depends on the carbon and nitrogen abundances in the Sun. This strong dependence on the SSM parameters leads to a \(30\%\) difference between the predictions of the CNO neutrino flux of the HZ and LZ SSMs (See Table 1).
In order to assess the relevance of a future CNO flux determination in this context, we performed a frequentist hypothesis test analogue to the one described in [7] and considered the HZ and LZ SSMs as the alternative hypothesis (\(H_1\)) and the null hypothesis (\(H_0\)), respectively. Namely, we assumed a given model and determined the p value against the alternative hypothesis when the \({^{8}\mathrm{B}}\), \({^{7}\mathrm{Be}}\), and CNO neutrino rates were measured with prescribed accuracies. For \({^{8}\mathrm{B}}\) and \({^{7}\mathrm{Be}}\) neutrinos, we adopted the rate uncertainties from [7]. This makes it possible to directly compare our results with the one from [7] that did not incorporate CNO neutrinos in the analysis. To simulate future CNO neutrino measurements, we used uncertainties of 1.5 and 0.5 \(\hbox {cpd}/100~\hbox {ton}\) , as discussed in the previous sections.
Our results are shown in Fig. 12, where we plot the distributions of the test statistic defined as the difference between the \(\chi ^2\) computed for HZ and LZ predictions including both model and experimental uncertainties. The more the distributions are separated, the higher the probability of discriminating among the different hypotheses is.
When taking into account values for the \(^{8}\hbox {B}\) and \(^{7}\hbox {Be}\) neutrino fluxes as predicted by the HZ SSM and measured with the same precision as in [7], the median p value for the LZ predictions is 0.057, corresponding to a \(1.6\sigma \) significance in the exclusion of the wrong hypothesis. When including a measurement of CNO neutrinos with an uncertainty of 1.5 \(\hbox {cpd}/100~\hbox {ton}\) , the median discriminatory power for the LZ SSM hypothesis does not change significantly (\(p\text { value} = 0.047\), \(1.7\sigma \) significance), as the experimental uncertainty is 2–3 times larger than the model precision. A larger increase in significance can be obtained if a 0.5 \(\hbox {cpd}/100~\hbox {ton}\) precision in the determination of the CNO neutrino flux is achieved. In this case, assuming the HZ SSM predictions, the median p value for the LZ SSM is 0.016 (\(2.1\sigma \)). Similar significance levels were found when using the values of \(\varPhi _{{}^{8}\mathrm{B}}\) and \(\varPhi _{{}^{7}\mathrm{Be}}\) measured by Borexino in [7].
This test relies on stronger assumptions as the determination of C+N abundance presented in the previous section. It also considers a simplified “binary” hypothesis system where either the HZ and LZ SSM is assumed to be correct, with no possible alternative explanation. However, its discriminatory power is better than the one achievable from the inferred C+N abundance for the simple reason that while in the latter the impact of metallicity on the core temperature is cancelled out, in this study it is the dominant source of information, driving the sensitivity of the test through \(\varPhi _{{}^{8}\mathrm{B}}\) and \(\varPhi _{{}^{7}\mathrm{Be}}\). In this context, a measurement of CNO neutrinos is expected to have a sensible impact on the test only if the experimental precision approaches model uncertainty (i.e., \(\sigma _\mathrm{CNO}\simeq \)0.5 \(\hbox {cpd}/100~\hbox {ton}\) ). The precision of the solar model predictions indeed poses a strong limit to the overall sensitivity of such a test. Currently, the largest sources of uncertainty – besides the plasma opacity and the carbon abundance – are the \(S_{17}\) (\(^{8}\hbox {B}\)), \(S_{34}\) (\(^{7}\hbox {Be}\)), and \(S_{114}\) (CNO) astrophysical Sfactors [6]. Therefore, a reduction of the nuclear cross section uncertainties is crucial to allow for a more significant test of the SSM based on solar neutrino measurements.
Conclusions
Borexino has a a strong potential to detect CNO neutrinos due to its unique radiopurity, intrinsically low cosmogenic \({}^{11}\hbox {C}\) rate and successful TFC method for its further reduction, reduction of external gammas by passive detector shielding, as well as comprehensive detector and background modeling.
A detailed sensitivity study was performed which showed that if the backgrounds are constrained, the bulk of the sensitivity to the CNO neutrino signal can be established with a simple counting analysis. The multivariate analysis only mildly improves this sensitivity when the uncertainty on the background rates is relatively loose. The statistical uncertainty on the CNO rate depends on the uncertainty in the determination of the pep neutrinos and \(^{210}\text {Bi}\) backgrounds. If \(R_\mathrm{Bi}\) and \(R_{{ pep}\,}\) can be constrained with a 1.5 \(\hbox {cpd}/100~\hbox {ton}\) and 0.04 \(\hbox {cpd}/100~\hbox {ton}\) precision, respectively, the expected uncertainty of CNO neutrino rate is 1.2 \(\hbox {cpd}/100~\hbox {ton}\) , and the expected median discovery significance to CNO neutrinos is \(4.1\sigma \) (\(3.0\sigma \)) for the CNO flux assuming the HZ (LZ) SSM.
The sensitivity does not improve significantly if a new purification campaign were to further reduce the \(^{210}\text {Bi}\) content, unless the absolute precision of the \(^{210}\text {Bi}\) constraint is improved. This improvement is not guaranteed by the reduction of the \(^{210}\text {Bi}\) content alone and is mainly related to the capability of stabilizing the detector temperature and suppressing convective motions. However, for a given uncertainty on the \(^{210}\text {Bi}\) constraint and if the uncertainty on the rate of CNO neutrinos is dominated by statistical fluctuations, a reduction of \(^{210}\text {Bi}\) would shorten the time needed to achieve a given precision on the CNO neutrino interaction rate.
In order to break the correlations, the pep and \(^{210}\text {Bi}\) background rates must be independently constrained. The rate of pep neutrinos can be constrained by exploiting its relation with the pp neutrino rate as already done in the analysis presented in [7]. The uncertainty on the pep rate can be reduced to the level of around one percent by imposing an additional constraint based on solar luminosity measurements.
The \(^{210}\text {Bi}\) isotope is supported by their longlived parent \(^{210}\text {Pb}\) (\(\tau _\text {Pb}=32.2\) years). The most promising method to constrain the \(^{210}\text {Bi}\) background relies on a measurement of the daughter isotope \(^{210}\text {Po}\) , which is expected to be in equilibrium with the \(^{210}\text {Pb}\) when the liquid scintillator is thermally stable. \(^{210}\text {Po}\) events can efficiently be tagged using pulseshape discrimination. Since mid2015, the Borexino collaboration has been focusing its effort on the thermally stabilizing the detector.
By combing the measurements of \(\varPhi _{{}^{8}\mathrm{B}}\) and \(\varPhi _\mathrm{CNO}\), the carbon and nitrogen contents of the solar core can be determined independently of the opacity of solar plasma. Assuming \(R_{\mathrm{CNO}}\) is measured with an uncertainty \(\simeq \) 0.5 \(\hbox {cpd}/100~\hbox {ton}\) , the C+N fraction of the Sun can be determined within \(\simeq 15\%\), which is comparable to the precision of the current estimations obtained from measurements of the photosphere, and limited by the nuclear cross section precision.
The neutrino fluxes depend on the chemical composition of the Sun, and can be used to discriminate between SSM HZ and LZ models. Combining the measurements of \(\varPhi _{^{7}\text {Be}\,},\,\varPhi _{{}^{8}\mathrm{B}}\), and \(\varPhi _\mathrm{CNO}\), assuming the CNO rate to be determined at a precision of 0.5 \(\hbox {cpd}/100~\hbox {ton}\) , the median discrimination sensitivity is 2.1\(\sigma \) and limited by the uncertainties of SSM.
Data Availability Statement
This manuscript has associated data in a data repository. [Authors’ comment: The datasets generated during the current study are freely available in the repository https://bxopen.lngs.infn.it/. Additional information is available from the Borexino Collaboration spokesperson (spokespersonborex@lngs.infn.it) upon reasonable request.]
Notes
 1.
The radiative opacity is represented by two parameters, namely \(\kappa _{a}\equiv 1+a\) and \(\kappa _b \equiv 1+b\), which describe the variations of the solar opacity profile in terms of the parameters a and b defined in [6].
References
 1.
D. Clayton, Principles of stellar evolution and nucleosynthesis (University of Chicago Press, Chicago, 1983)
 2.
J. Bahcall, Neutrino astrophysics (Cambridge University Press, Cambridge, 1989)
 3.
M. Salaris, S. Cassisi, Evolution of stars and stellar populations (Wiley, New York, 2005)
 4.
R. Kippenhahn, A. Weigert, A. Weiss, Stellar structure and evolution (Springer, New York, 2012)
 5.
S. Degl’Innocenti, G. Fiorentini, B. Ricci, F. Villante, Phys. Lett. B 590(1), 13–20 (2004)
 6.
N. Vinyoles et al., Astrophys. J. 835(2), 202 (2017)
 7.
M. Agostini et al. (Borexino Coll.), Nature 562, 505–510 (2018)
 8.
F. Villante et al., Phys. Lett. B 701(3), 336–341 (2011)
 9.
M. Asplund, N. Grevesse, A. Sauval, Nucl. Phys. A 777, 1–4 (2006) (special issue on Nuclear Astrophysics)
 10.
M. Asplund, A. Grevesse, N. Sauval, P. Scott, Annu. Rev. Astron. Astrophys. 47(1), 481–522 (2009)
 11.
P. Scott et al., Astron. Astrophys. 573(A25), 1–19 (2015)
 12.
P. Scott, M. Asplund, N. Grevesse, M. Bergemann, A. Sauval, Astron. Astrophys. 573(A26), 1–33 (2015)
 13.
E. Caffau, H.G. Ludwig, M. Steffen, B. Freytag, P. Bonifacio, Sol. Phys. 268(2), 255–269 (2011)
 14.
N. Grevesse, A. Sauval, Space Sci. Rev. 85(1), 161–174 (1998)
 15.
A. Serenelli, C. Peña Garay, W.C. Haxton, Phys. Rev. D 87, 1–9 (2013)
 16.
W.C. Haxton, A.M. Serenelli, Astrophys. J. 687(1), 678–691 (2008)
 17.
K. Abe et al. (SuperKamiokande Coll.), Phys. Rev. D 94(5), 052010 (2016)
 18.
G. Bellini et al. (Borexino Coll.), Phys. Rev. D 89(11), 1–68 (2014)
 19.
M. Agostini, Zenodo (2020). https://doi.org/10.5281/zenodo.4209228
 20.
G. Alimonti et al. (Borexino Coll.), Nucl. Instrum. Methods A 600, 3, 568–593 (2009)
 21.
M. Agostini et al. (Borexino Coll.), Astropar. Phys. 92, 21–29 (2017)
 22.
M. Agostini et al. (Borexino Coll.), Phys. Rev. D 101(1), 12009 (2020)
 23.
I. Esteban et al., J. High Energy Phys. 2017(1), 87 (2017)
 24.
M. Agostini et al. (Borexino Coll.), Astropar. Phys. 97, 136–159 (2018)
 25.
X. Ding, J. Instrum. 13(12), P12018–P12018 (2018)
 26.
M. Agostini, MSTATS: framework for frequentist statistical analysis (2016). https://github.com/mmatteo/mstats
 27.
G. Cowan, K. Cranmer, E. Gross, O. Vitells, Eur. Phys. J. C 71(2), 1–19 (2011)
 28.
J. Bahcall, R. May, Astrophys. J. 155, 501 (1969)
 29.
J. Bahcall, Phys. Rev. D 41, 2964–2966 (1990)
 30.
A. Kurylov, M. RamseyMusolf, P. Vogel, Phys. Rev. C 67, 12 (2003)
 31.
E.G. Adelberger et al., Rev. Mod. Phys. 83, 195–245 (2011)
 32.
F. Capozzi, E. Lisi, A. Marrone, A. Palazzo, J. Phys. Conf. Ser. 1312, 012005 (2019). https://doi.org/10.1088/17426596/1312/1/012005
 33.
J. Bergström et al., J. High Energy Phys. 2016, 132 (2016)
 34.
J. Bahcall, Phys. Rev. C 65, 5 (2002)
 35.
V. Castellani, S. Degl’Innocenti, G. Fiorentini, M. Lissia, B. Ricci, Phys. Rept. 281, 309–398 (1997)
 36.
F. Vissani, in ed. by M. Meyer, K. Zuber (eds.) Solar Neutrinos (2019), pp. 121–141. arXiv:1808.01495
 37.
A. Ianni, Acta Phys. Pol., A 121(1), 101–103 (2012)
 38.
M. Agostini et al. (Borexino Coll.) (2020)
 39.
W.A. BravoBerguño et al., Nucl. Instrum. Methods Phys. Res. Sect. A 885, 38–53 (2018). https://doi.org/10.1088/17426596/1312/1/012005
 40.
J. E. Bailey et al., Nature 517(7532), 56–59 (2015)
 41.
M. Krief, Y. Kurzweil, A. Feigel, D. Gazit, Astrophys. J 856(2), 135 (2018)
 42.
J. Colgan et al., Astrophys. J 817(2), 116 (2016)
 43.
C.A. Iglesias, Mon. Not. R. Astron. Soc. 450(1), 2–9 (2015)
 44.
M. Krief, A. Feigel, D. Gazit, Astrophys. J 821(1), 45 (2016)
 45.
C. Blancard, P. Cossé, G. Faussurier, Astrophys. J 745(1), 1 (2012)
 46.
J. Bahcall, S. Basu, M. Pinsonneault, A. Serenelli, Astrophys. J. 618, 623–626 (2005)
 47.
J. ChristensenDalsgaard, M. Di Mauro, G. Houdek, F. Pijpers, Astron. Astrophys. 494(1), 205–208 (2009)
 48.
F. Villante, Astrophys. J. 724, 623–626 (2010)
 49.
D. Franco et al., J. Cosmol. Astropart. P. 2016, 8 (2016). arXiv:1510.04196
 50.
M. Askins, Z. Bagdasarian, N. Barros et al., THEIA: an advanced optical neutrino detector. Eur. Phys. J. C 80, 416 (2020). https://doi.org/10.1140/epjc/s1005202079778
 51.
J.F. Beacom et al., Chin. Phys. C 41, 2 (2017)
 52.
Jülich Supercomputing Centre, J. Large Scale Res. Facilities 4, A132 (2018)
Acknowledgements
The Borexino program is made possible by funding from Istituto Nazionale di Fisica Nucleare (INFN) (Italy), National Science Foundation (NSF) (USA), Deutsche Forschungsgemeinschaft (DFG) and HelmholtzGemeinschaft (HGF) (Germany), Russian Foundation for Basic Research (RFBR) (Grants No. 162913014ofim, No. 170200305A, and No. 190200097A) and Russian Science Foundation (RSF) (Grant No. 171201009) (Russia), and Narodowe Centrum Nauki (NCN) (Grant No. UMO 2017/26/M/ST2/00915) (Poland). We gratefully acknowledge the computing services of Bologna INFNCNAF data centre and ULite Computing Center and Network Service at LNGS (Italy), and the computing time granted through JARA on the supercomputer JURECA [52] at Forschungszentrum Jülich (Germany). This research was supported in part by PLGrid Infrastructure. We acknowledge also the generous hospitality and support of the Laboratori Nazionali del Gran Sasso (Italy).
Author information
Affiliations
Consortia
Corresponding authors
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP^{3}
About this article
Cite this article
Agostini, M., Altenmüller, K., Appel, S. et al. Sensitivity to neutrinos from the solar CNO cycle in Borexino. Eur. Phys. J. C 80, 1091 (2020). https://doi.org/10.1140/epjc/s10052020085342
Received:
Accepted:
Published: