Sensitivity of a low threshold directional detector to CNO-cycle solar neutrinos

A first measurement of neutrinos from the CNO fusion cycle in the Sun would allow a resolution to the current solar metallicity problem. Detection of these low-energy neutrinos requires a low-threshold detector, while discrimination from radioactive backgrounds in the region of interest is significantly enhanced via directional sensitivity. This combination can be achieved in a water-based liquid scintillator target, which offers enhanced energy resolution beyond a standard water Cherenkov detector. We study the sensitivity of such a detector to CNO neutrinos under various detector and background scenarios, and draw conclusions about the requirements for such a detector to successfully measure the CNO neutrino flux. A detector designed to measure CNO neutrinos could also achieve a few-percent measurement of pep neutrinos.


Solar neutrinos
Solar neutrino experiments were pivotal in the groundbreaking discovery of neutrino oscillation and, hence, massive neutrinos, while at the same time confirming our understanding of fusion processes in the Sun. The Sudbury Neutrino Observatory (SNO) experiment resolved the so-called Solar Neutrino Problem by detecting the Sun's "missing" neutrinos, confirming the theory of neutrino flavor change. The combination of a charged-current (CC) measurement from SNO with SuperKamiokande's high-precision elastic scattering (ES) measurement demonstrated that the electron neutrinos produced in the Sun were transitioning to other flavors prior to detection [1], a result later confirmed at 5σ by SNO's measurement of the flavor-independent 8 B flux using the neutral current (NC) interaction [2]. The KamLAND a e-mail: rbonventre@lbl.gov reactor experiment confirmed this flavour change as being due to oscillation [3]. This opened the door to a precision regime, allowing neutrinos to be used to probe the structure of the Sun, as well as the Earth and far-distant stars. Solar neutrinos remain the only sector of neutrinos with a confirmed observation of the effect of matter on neutrino oscillation at high significance, providing a unique opportunity to further probe this interaction to search for non-standard interactions and other effects.
The Borexino experiment made the first direct measurements of the 7 Be, pp and pep fluxes [4][5][6], as a result of which the pp fusion cycle in our Sun has been well studied, with measurements of all neutrino sources bar the highenergy hep neutrinos. The subdominant CNO cycle is less well understood and yet has the potential to shed light on remaining mysteries within the Sun. One of the critical factors that engendered confidence in the Standard Solar Model [7] (SSM) was the excellent agreement (∼ 0.1%) of SSM predictions for the speed of sound with helioseismological measurements. The speed of sound predicted by the SSM is highly dependent on solar dynamics and opacity, which are affected by the Sun's composition [8]. In recent years the theoretical prediction for the abundance of metals (elements heavier than H or He) in the photosphere has fallen due to improvements in the modeling of the solar atmosphere, including replacing previous one-dimensional models with fully three-dimensional modeling, and inclusion of effects such as stratification and inhomogeneities [9]. The new results are more consistent with neighboring stars of similar type, and yield improved agreement with absorptionline shapes [10], but at the same time reduce the prediction for the metal abundance by ∼30%. When these new values are input to the SSM, the result is a discrepancy in the speed of sound with helioseismological observations. This new disagreement has become known as the "Solar Metallicity Problem" [11,12]. A measurement of the CNO neutrino flux may help in resolving this problem [13,14]. The impact of metal-licity on pp-chain neutrinos is small relative to theoretical uncertainties, but the neutrino flux from the sub-dominant CNO cycle depends linearly on the metallicity of the solar core, and the predictions for the two models differ by greater than 30% [15]. The theoretical uncertainty on these predictions is roughly 14-18% in the so-called AGS05-Opt model, although greater in other models [15]. However, these uncertainties can be reduced to < 10% using correlations in the theoretical uncertainties between the CNO and 8 B neutrino fluxes: the two have similar dependence on environmental factors, thus a precision measurement of the 8 B neutrino flux can be used to "calibrate" the core temperature of the Sun and, thus, constrain the CNO neutrino flux prediction [13]. In [13], the final uncertainty is dominated by the nuclear physics. A precision measurement of the CNO flux then has the potential to resolve the current uncertainty in heavy element abundance.
Borexino has placed the most stringent limits on the CNO neutrino flux to date [6], and continues to pursue a first observation. However, extraction of this flux is extremely challenging due to the similarity of the spectrum of ES recoil electrons with background 210 Bi decays in the target. Borexino propose to use the time evolution of the α decay of the daughter, 210 Po, to constrain the level of 210 Bi. This method requires both a stable α-particle detection efficiency and a lack of external sources of the 210 Po daughter, which can be challenging to achieve [16]. A recent paper discusses the sensitivity of several current and future experiments to the CNO flux [17]. A detector with directional sensitivity could discriminate between the directional solar neutrino signal and the isotropic background, without the need for a time-series analysis.

Low-threshold directional detection
Water Cherenkov detectors (WCD) are limited in energy threshold by the relatively low Cherenkov photon yield. Scintillator-based detectors can achieve the thresholds required to observe CNO neutrinos, but lose the advantage of the directional information provided by Cherenkov light. The novel water-based liquid scintillator (WbLS) target medium [18] offers the potential to benefit from both the abundant scintillation and directional Cherenkov signals, thus achieving a massive, low-threshold directional detector.
The sensitivity of a 50 kT pure LS detector has been studied by the LENA collaboration [19]. While the threshold of a WbLS detector will not be as low as for a pure LS target (LENA studies assumed a 250 keV threshold), the additional information provided by the Cherenkov component provides a strong benefit in signal/background separation. WbLS offers a uniquely broad, multi-parameter phase space that can be optimized to maximize sensitivity to a particular physics goal. The WbLS "cocktail" can range from a high-LS fraction, oil-like mixture, with > 90% LS, to a water-like mixture with anything from 10% to sub-percent levels of LS. The choice of fluor affects the scintillation yield, timing, and emission spectrum, and metallic isotopes can be deployed to provide additional targets for neutrino interaction [20].
In this article we study the sensitivity of a large WbLS detector to CNO neutrinos under a range of detector scenarios, including target size, LS fraction, photo-coverage, angular resolution, and the level of intrinsic radioactive contaminants. These studies can inform the design of a future experiment targeting a CNO flux observation, such as Theia [21,22].
There is much interest in the community in developing low-threshold directional detectors. Monte Carlo studies in [23,24] discuss how the potential for separation of a Cherenkov signal in a scintillating target could be used to extract particle direction. The CHESS experiment has recently demonstrated first detection of a Cherenkov signal in pure LS (both LAB and LAB/PPO) [25,26]. While highenergy muons were used for this demonstration, they were in the MIP regime and thus the energy deposited along the fewcm track in the CHESS target was only a few MeV, within the regime relevant for this work. Studies based on data from the KamLAND detector show the potential for directional reconstruction using time-of-flight of the isotropic scintillation light [27,28].
In Sect. 2 we describe the analysis method for evaluating the uncertainty on the CNO and pep neutrino fluxes, and the simulation of each signal and background source. In Sect. 3 we describe various scenarios for both detector configuration and background assumptions under which the CNO flux is evaluated. Section 4 presents the results, and Sect. 5 describes the conclusions.

Analysis methods
Neutrino flux sensitivities are determined using a binned maximum likelihood fit over two-dimensional PDFs in energy and direction relative to the Sun, cos θ . The energy dimension allows separation of neutrino fluxes from each other, as well as discrimination from certain background events. The direction dimension is critical for a full separation of the CNO flux from radioactive background.

Simulation of expected signals
Simulations for this paper were produced using the RAT-PAC software (https://github.com/rat-pac/rat-pac), which is based on Geant4. Optical properties for the WbLS were constructed from weighted combinations of water and scintillator optics as determined by the SNO+ collaboration. The optical simulation was tuned to allow multicomponent absorption and reemission with separate absorption lengths and reemis- sion probabilities for each component of the WbLS cocktail. Radioactive decays were simulated using the decaychain generator developed by Joe Formaggio and Jason Detwiler (Private communication), and solar neutrino interactions were simulated using an elastic scattering generator also developed by Joe Formaggio. Solar signals were simulated assuming fluxes from the BS05OP solar model [7], which assumes the higher solar metallicity, and using LMA-MSW survival probabilities from the three flavor best fit oscillation values from [29]. The reconstructed effective electron energy spectrum for each signal was determined semi-analytically. First, the distribution of the number of PMT hits (NHit) per event for each signal was found, handling the scintillation and Cherenkov components separately. Given a WbLS cocktail and detector configuration, the scintillation contribution to the NHit for an event at a specific position scales linearly with the quenched energy deposition. This scaling was determined by simulating electrons at each position with the Cherenkov light production disabled. The Cherenkov contribution to the NHit was determined by simulating each signal with the scintillation light yield set to zero, but with absorption and reemission of the Cherenkov light in the scintillator and wavelength shifter enabled. The expected NHit for each event was taken to be the sum of the Cherenkov Nhit plus a number of scintillation hits drawn from a Poisson distribution with a mean given by the scintillator energy deposition in that event times the scaling factor for the relevant event position. This method made efficient use of computing resources to simulate a full set of background event types. The result of this procedure was compared with the NHit distribution from a full simulation for 1 MeV electrons, and both the mean and width were found to agree to within 0.5%, as shown in The conversion from NHit to reconstructed effective electron energy was determined using a position-and directiondependent lookup table. The lookup table was generated by simulating electrons at various positions, directions, and energies using the above procedure. The resultant energy PDFs for 5% scintillator are shown in Fig. 2.
For each signal the cos θ distribution was determined fully analytically. All non-neutrino signals were assumed to be flat. For the solar signals the ES generator determined the electron direction relative to the Sun based on the differential cross sections. This was then convolved with a chosen angular resolution (Sect. 3).

Baseline fit configuration
The normalizations for 11 signals are floated in the fit. 238 U and 234 Th chain backgrounds are assumed to be in equilibrium except for 210 Bi, 210 Po, and 210 Pb, and the backgrounds in each chain are floated together as a single parameter. The various flavor components of the 8 B, 7 Be, pep, and CNO signals are combined into one parameter per flux, assuming survival probabilities from [29]. The CNO signal contains the sum of the 17 F, 15 O, and 13 N solar neutrino signals. The pp solar signal is not included as it falls below the energy threshold, and the hep solar signal is fixed in magnitude as it is too small to be reconstructed. 210 Bi, 40 K, 85 Kr, and cosmogenically activated 11 C are each included as a separate parameter. The 39 Ar and 210 Po backgrounds are floated together as their energy spectra above 600 keV are similar.
The baseline fit uses 40 bins in cos θ and 20 keV bins in energy from 600 keV to 6.5 MeV. A 5 year livetime is assumed. Projection of two dimensional fit in cos θ for energies between 1 and 1.5 MeV for a randomly generated fake dataset with the baseline detector configuration and background assumptions A 50% fiducial volume cut is applied in order to reduce the impact of background contaminants in external regions of the detector, such as γ s from 208 Tl in the PMTs, to negligible.
Full bias and pull studies were performed, and the fit was observed to be unbiased, with the expected pull distribution. Figure 3 shows the pull distribution for fits with the baseline configuration. Figure 4 shows the projection of one two-dimensional fit in the cos θ dimension. This figure illustrates the importance of angular resolution in extracting the solar neutrino signal even with backgrounds many order of magnitudes larger, using the high statistics achievable in a large detector.

Detector configuration and background assumptions
The sensitivity to CNO neutrinos is studied under a range of detector scenarios and background assumptions: For the purposes of comparison we define the baseline detector configuration to be a 50-kT detector with 90% PMT coverage, 5% WbLS, and 25 • angular resolution, with baseline background levels as given in Table 2 in Sect. 3. All results assume a five year livetime.
Target volume We consider both a 25 and 50 kT total detector volume, corresponding to a 31.7 or 40-m sized rightcylindrical vessel. The PMTs are positioned at the edge of this volume, and a fiducial volume is selected for analysis (Sect. 2). The 50% fiducial volume corresponds to a 4.15-m buffer between the PMTs and the target volume for the 50 kT detector, and a 3.27-m buffer for the 25 kT detector.
Increasing the target volume scales the exposure accordingly, but at the same time reduces the overall light collection of the detector due to absorption, thus impacting the achievable energy resolution. This additional absorption would also negatively impact the angular resolution. Since this work does not perform a full directional reconstruction, this correlation is not explicitly included. However, a range of possible angular resolutions are considered for each target volume.
Angular resolution A critical factor in this work is the assumed angular resolution. The angle between the incoming particle direction and the direction to the Sun, cos θ , can be used to differentiate signal from background. Due to the kinematics of the ES interaction, solar neutrinos are predominantly directed away from the Sun. This provides a key handle to discriminate solar neutrino events from an isotropic radioactive background.
All radioactive backgrounds are assumed to be isotropically distributed and thus have a flat distribution in cos θ . For solar neutrinos the direction of the electron relative to the Sun was determined as a function of energy via simulations that take as input the full differential cross sections. The resulting electron direction was then convolved with a detector angular resolution of the form with a width given by σ . The angular resolution was assumed to be constant with energy, with the value defined at threshold. Any improvement in sensitivity as the angular resolution improves at energies significantly above the energy threshold was observed to be a second order effect. This work does not attempt a full reconstruction of event direction. Instead, we consider a range of possible angular resolutions in order to determine the impact on the final neutrino flux sensitivity. In the best case, we consider a resolution of 25 • -similar to that achieved by SNO. In SNO, a 1-kT heavy-water detector with 55% PMT coverage, an angular resolution of 26.7 • was achieved for 16 N events (approximately 5 MeV), which had an average of 36 PMT hits [30]. In Super Kamiokande, a 50-kT light-water detector with 40% coverage, an angular resolution of approximately 35 • was achieved at 6 MeV, with 41 PMT hits [31]. Table 1 shows the expected number of hits from Cherenkov photons in the proposed 0.5% WbLS detector for a range of energy thresholds. As shown, by 1 MeV we expect as many hits from Cherenkov photons as in SNO and Super Kamiokande at 5 or 6 MeV. While many of these photons may be scattered, or absorbed and reemitted by the scintillator, absorption lengths in WbLS are significantly longer than a pure LS detector. Additionally, some directional information is retained by considering the offset of the reemitted photon from the original production point. The proposed detector has several further advantages that could allow for greatly improved angular resolution at lower energies. Increased coverage and use of high quantum efficiency PMTs allows detection of many more direct Cherenkov photons than in SNO and Super Kamiokande at equivalent energies. The inclusion of wavelength shifter in the WbLS absorbs and reemits Cherenkov photons that would otherwise be at too small a wavelength to be detected by the PMTs and, depending on the absorption length, may also retain some directional information. As demonstrated in [27,28], even a pure LS detector can provide some directional information by considering photon time of flight. With these advantages, a resolution of 25 • may be achievable.
As the detector size increases, coverage is reduced, or the LS fraction is increased, the resolution will naturally be degraded by increased scattering and absorption, and reduced light collection. We find that a large fraction of wavelength shifted photons are absorbed within a short distance, and so it is possible a lower concentration of PPO would be desirable to better retain directional information. To estimate the impact of these effects we consider degraded resolutions of 35 • , 45 • , and 55 • for each detector scenario.

WbLS cocktail
The scintillator component of the WbLS cocktail is taken to be LAB with 2g/L of PPO as a fluor (hereafter referred to as LAB/PPO). LAB/PPO properties have been determined by the SNO+ collaboration [32]. We consider fractions of LS from 0.5 to 5%, as well as a pure LS  The Cherenkov light yield is unaffected at first order by this change to the target cocktail, although there is a non-zero impact through absorption and reemission in the LS, which is fully modeled in the simulation. While this change can therefore be expected to affect the angular resolution, this effect is not included in the PDFs since this work does not perform a full directional reconstruction. Instead, we consider the impact of a range of angular resolutions for each target cocktail.
Photocathode coverage We assume instrumentation with Hamamatsu R11780 12 inch high quantum efficiency (HQE) PMTs, which have a peak efficiency of 32% at 390 nm [33]. We study photocathode coverages of up to 90%, which would require approximately 100k PMTs for the 50 kT detector. Changing the photocathode coverage effectively scales the  Figure 6 shows the PMT QE, and the emission spectra for both Cherenkov and scintillation light.

Energy scale and resolution
The energy resolution is determined by the overall light yield, which is dependent on the WbLS cocktail and the photocathode coverage, the impact of each of which is studied separately (Sect. 3). Here, we consider the impact of systematic uncertainties in the energy scale and resolution on the analysis. We investigate the effect of uncertainties in energy by modifying the PDFs relative to the spectra from which fake datasets are drawn. For the energy scale, a linear shift in energy is applied to the energy spectra. For the resolution a Gaussian smearing is applied to the energy spectra, with the width of the Gaussian set to a constant value.
Energy threshold The ability to reach a lower energy threshold than a water detector is one of the main advantages of a WbLS detector besides increased energy resolution. We consider energy thresholds from 600 keV up to 1 MeV. The possibility of reaching thresholds this low depend on the background rate, the PMT dark noise rate, the trigger setup, and the sustainable data rate. At 0.5% scintillator, we expect a total of 19.3 PMT hits at 0.6 MeV in a 50 kT detector, and by 5% scintillator we expect 93.0, compared to 5.4 PMT hits in a water only detector. Assuming a trigger window of 200 ns as used in Super Kamiokande, with ∼ 94,000 PMTs, we expect 18.8 noise hits per trigger window per kHz PMT dark rate, and so at 0.5% scintillator a dark rate lower than 1 kHz would be required. In a 25 kT detector we expect 21.1 PMT hits at 0.6 MeV, and 11.7 noise hits per trigger window per kHz PMT dark rate.  [4] a The 40 K level in water is taken to be 0.1× the Borexino measurement [34] b The 85 Kr, 39 Ar, and 210 Bi levels in water are taken to be the Borexino measured level in scintillator [36], although levels increased by several orders of magnitude are explored Background assumptions Radioactive contaminants in the target material are calculated as the sum of contamination from the LS and the water components, weighted by the corresponding mass fractions of the WbLS cocktail. Table 2 details the numbers assumed in the baseline analysis, taken from measurements by SNO and Borexino [34][35][36][37]. In the analysis we consider a range of levels for the intrinsic contamination in the WbLS target, as well as the degree of α-β separation and Bi-Po pile-up rejection achievable.
-U-and Th-chain: The LS components of uranium-and thorium-chain background levels are assumed to be at Borexino levels and the water components are assumed to be at the level of the heavy water in SNO [35,37]. The baseline alpha rejection is assumed to be 95%. For 212 Bi-212 Po and 214 Bi-214 Po events, an event window of 400 ns is assumed, with a 95% rejection for in-window coincident events and 100% rejection for tagging out of window events. While this level of discrimination has not yet been demonstrated in a WbLS target, this is a future goal of the CHESS experiment [25]. The discrimination achieved will depend on both the α quenching and overall light yield of the target, as well as specific timing properties. These microphysical properties must be fully understood in order to quantify the level of α and Bi-Po coincidence rejection that can be achieved. The impact of the efficiency of both α rejection and Bi-Po tagging is studied. -40 K: The level of 40 K in LS is a conservative estimate from the upper limit of Borexino's initial measurement. The level in water is taken from an upper limit measured in the Borexino Counting Test Facility (CTF) [34]. SNO measured a level of 2e−9 gK/gH 2 O in the light water, although this background was below threshold in SNO and thus little effort was made to reduce it [38]. These measurements are therefore taken as conservative upper bounds on the level. We use 0.1× the Borexino level in water as the baseline level for this study, and investigate the impact of a contamination an order of magnitude higher. -210 Bi: 210 Pb, 210 Bi, and 210 Po cannot be assumed to be in equilibrium with the rest of the uranium chain due to the long half life of 210 Pb and the possibility of Rn contamination. The baseline level of 210 Bi in LS is taken from Borexino. The level of out-of-equilibrium background achievable in ultra-pure water has not been measured to the precision needed for this kind of experiment. The uranium-chain contribution from the water component is orders of magnitude larger than the 210 Bi level measured in Borexino, thus the contribution from any out-ofequilibrium component in water must be many orders of magnitude larger than in LS in order to impact the sensitivity. Values of 10×, 100× and 1000× the Borexinomeasured value in scintillator are explored for the contamination in water. to shield from cosmogenic backgrounds. The 11 C level in Borexino, with 3800 m.w.e. overburden, is approximately 1e5 events per kiloton year. This is used as a conservative initial estimate for the rate, adjusted for the carbon content of the different target materials. Rates of an order of magnitude higher and lower are studied, to simulate the effect of different possible detector sites. Production of cosmogenic backgrounds on the water component of the target were considered according to [39], which provides a complete list of potential spallation products. The dominant sources in the energy range considered in this work are 11 C and 15 O. Inclusion of these backgrounds was observed to change the CNO sensitivity by 0.03% for the baseline configuration, and the pep sensitivity by an unobservable amount. Thus, these backgrounds were omitted for the remainder of this work. -Externals: Background contributions from radioactive contamination external to the target region (for example 208 Tl γ s from the PMTs and any support structures) are assumed to be negligible inside a chosen 50% fiducial volume. In future studies, vertex reconstruction could be used to constrain such sources and thus potentially expand the fiducial volume.

Results
The fit uncertainty for each signal with the baseline detector configuration and background assumptions is shown in Table 3.
Detector size, target, and angular resolution The CNO solar neutrino sensitivity as a function of detector size, LS fraction, and angular resolution for the baseline background assumptions is shown in Table 4 and Fig. 7. As a comparison we look at simulated spectra from a 50 kT pure LS detector and study the results of a one dimensional fit in energy (under the assumption that there would be no directional resolution in pure LS). Here we find a CNO sensitivity of 3.5%, but pull distributions show that the fit is not able to converge on the full errors in the CNO and 210 Bi signals (Fig. 8). The correlation between the two values in the fit is −0.84. This suggests than energy alone is not sufficient to distinguish these signals even at very high statistics. PMT coverage The above results show that the impact of the improved energy resolution due to a larger scintillator fraction (0.5-5%) is marginal. However, changes in PMT coverage can have a potentially greater effect. At low energies a significant fraction of the resolution comes from Cherenkov photons, which does not scale with scintillator fraction. Reducing the PMT coverage will reduce the angular resolution as well as the overall light collection. The fit uncertainty for CNO for 5 years of data using a detector with 60% PMT coverage instead of 90% is shown in Table 5.
We can see as predicted that at low scintillator fractions the change in PMT coverage has a larger effect than an equivalent fractional change in scintillator fraction. At higher scintillator fractions the effect is smaller, and suggests that the main consideration for the PMT coverage requirement will be the achievable angular resolution.  (Figs. 9, 10). For the 50 kT detector at 5% WbLS and 25 • angular resolution the energy scale uncertainty can be constrained to less than 0.006%, where the change in the fitted CNO normalization is 1.7% (as shown in Fig. 9).
By smearing the PDFs with a Gaussian of pre-defined width, we investigate the impact of an uncertainty on the energy resolution. A scan of the likelihood space for the 50 kT detector at 5% WbLS and 25 • angular resolution demonstrates the capability to constrain this uncertainty with the data itself down to 5.5 keV, at which point the systematic change in the fitted CNO normalization is 3.6%, as shown in Fig. 10.

Energy threshold
The CNO solar neutrino sensitivity as a function of the energy threshold for the likelihood fit is shown in Table 6. These results show that a precise measurement of the CNO flux requires maintaining sensitivity below 1 MeV. It can be seen in Fig. 2 that this energy region is where the CNO and pep solar neutrino signals can be distinguished from each other, and at higher energy thresholds these signals become highly correlated.
Background assumptions As seen in Fig. 2, 40 K is a dominant background at low energies, even assuming an order of magnitude improvement over the level measured in water by Borexino and SNO. This is due to the much higher contamination in the water component of the WbLS compared to the   relatively cleaner scintillator. The 40 K background in water was not critical for these previous measurements, and so it may be possible to further reduce the level with additional effort. The SNO water processing plant could be improved by increasing the frequency of replacing ion exchange columns or by distilling the water. Table 7 shows that at the Borexino measured level this background greatly decreases the sensitivity in all configurations compared to the baseline background assumptions which includes the order of magnitude reduction. Table 8 shows the sensitivity to each floated signal in the fit for the baseline detector configuration.
Due to the small fraction of scintillator in the WbLS, the relative 11 C background is already reduced by an order of magnitude or more compared to pure scintillator experiments. Combined with directional reconstruction, the solar sensitivity becomes mostly insensitive to the 11 C event rate at around the Borexino level, as shown in Table 9. This suggests that the overburden is not critical for a solar measurement with an experiment such as Theia. Table 9 Fit uncertainty on the CNO neutrino flux for 5 years of data with the baseline detector configuration as a function of isotope contamination level relative to baseline background level. For 210 Bi this is the level of the out-of-equilibrium component in water relative to the out-of-equilibrium level in scintillator. For 85 Kr and 39 Ar this is the level of contamination in water relative to baseline Fraction of nominal Isotope contamination 11  No BiPo in-window rejection 5.4 No α, no BiPo in-window rejection 5.4 No α, no BiPo rejection 6.4 210 Pb, 210 Bi, and 210 Po are not necessarily in equilibrium with the rest of the uranium chain due to the long halflife of 210 Pb and the possibility of Rn contamination. The level of out-of-equilibrium background achievable in ultrapure water has not been measured to the precision needed for this kind of experiment. The uranium chain contribution from the water component is already orders of magnitude larger than the non-equilibrium level measured in Borexino, so the nonequilibrium component in water must be many orders of magnitude larger than in scintillator as well to have any effect on the sensitivity. Table 9 shows that there begins to be a moderate impact when the out-of-equilibrium level in water is 1000× the out-of-equilibrium level in scintillator. 85 Kr and 39 Ar are also found to have little impact on the sensitivity at up to 10,000× the expected level in scintillator, as shown in Table 9.
We find that alpha and Bi-Po rejection is not critical for this measurement. The largest effect comes from the ability to reject Bi-Po coincidences that occur in separate trigger windows, but there is little impact even assuming zero rejection, as shown in Table 10.

Conclusions
The feasibility of a low energy solar neutrino measurement with a large WbLS detector depends on the backgrounds and angular resolution achievable. At currently measured levels, the limiting background appears to be 40 K in water. As this has not been a critical background in previous measurements, it is possible that lower levels may be achievable with more effort. The remaining unknown backgrounds in water have been shown to have little effect if they are kept to 1000-10,000 times the level achieved in scintillator. The 11 C background is shown to be unimportant, which suggests that the overburden does not matter for this measurement. Finally, changes in scintillator and PMT coverage have been shown to have relatively small effects, which suggest that even by a 0.5% scintillator fraction, energy resolution is no longer the critical parameter to optimize. Instead, the impact of these changes on threshold, angular resolution, and energy systematic uncertainties will be the deciding factor. With a baseline detector of 50-kT total volume (50% fiducial), 90% PMT coverage and a 5% WbLS target, assuming a 25 • angular resolution, a precision of several percent is possible for both CNO and pep neutrino fluxes.
Further studies hinge on additional R&D, including the Cherenkov detection efficiency enhancement provided by deployment of fast photon sensors [40][41][42][43], and demonstration of quenching and particle ID capabilities in WbLS. Development of a directional reconstruction algorithm would allow a direct demonstration of the required angular resolution discussed in this article. An exciting avenue for further exploration would be isotope loading of the WbLS target. Loading the WbLS target with an isotope such as 7 Li for CC detection would provide an additional handle for signal/background separation via improved spectral information. Neutrinos interact in a pure (Wb)LS detector via ES. While the ES differential cross section is almost maximally broad, thus providing little handle on the underlying neutrino spectrum, that for the CC interaction on 7 Li is extremely sharply peaked, resulting in the potential for high-precision measurement of the underlying neutrino energy spectrum. 7 Li was proposed as an additive to a WCD in [44] for this reason; addition of this isotope to a WbLS detector would further improve the discrimination power. Initial studies are presented in [22]. A more quantitative study could yield significant improvements to the results presented here.
to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .