Dark matter imprint on $^8$B neutrino spectrum

The next generation of solar neutrino detectors will provide a precision measure of the $^8$B electron-neutrino spectrum in the energy range from 1-15 MeV. Although the neutrino spectrum emitted by $^8$B $\beta$-decay reactions in the Sun's core is identical to the neutrino spectrum measured in the laboratory, due to vacuum and matter flavor oscillations, this spectrum will be very different from that measured on Earth by the different solar neutrino experiments. We study how the presence of dark matter (DM) in the Sun's core changes the shape of the $^8$B electron-neutrino spectrum. These modifications are caused by local variations of the electronic density and the $^8$B neutrino source, induced by local changes of the temperature, density and chemical composition. Particularly relevant are the shape changes at low and medium energy range $(E_\nu\le 10 {\; \rm MeV})$, for which the experimental noise level is expected to be quite small. If such a distortion in the $^8$B$\nu_e$ spectrum were to be observed, this would strongly hint in favor of the existence of DM in the Sun's core. The $^8$B electron-neutrino spectrum provides a complementary method to helioseismology and total neutrino fluxes for constraining the DM properties. In particular, we study the impact of light asymmetric DM on solar neutrino spectra. Accurate neutrino spectra measurements could help to determine whether light asymmetric DM exists in the Sun's core, since it has been recently advocated that this type of DM might resolve the solar abundance problem.


I. INTRODUCTION
Solar neutrino detectors have been one of the beacons of particle physics, both by leading the way in discovering the basic properties of particles, including the nature of neutrino flavor oscillations, and by being responsible for developing pioneering techniques in experimental neutrino detection [e.g., [1][2][3]. The next generation of detectors like the DUNE Experiment [4], the CJPL Laboratory [5], the JUNO Observatory [6], and the LENA detector [3], will measure with high precision the neutrino fluxes and neutrino spectra of a few key solar nuclear reactions, such as the electron-neutrino ( 8 Bν e ) spectrum produced by the β-decay of 8 B [7,8]. This will allow us to probe in detail the Sun's core, including the search for unknown physics processes. Moreover, the high quality of the data will enable the development of inversion techniques for determining basic properties of the solar plasma [e.g. , 9]. Specific examples can be found in Balantekin et al. [10], Davis [11] and Lopes [12]. Equally, solar neutrino data can be used to find specific features associated with new physical processes [e.g. 13,14], such as the possibility of an isothermal solar core associated with the presence of DM [15].
The 8 Bν e spectrum emitted by the nuclear reactions in the Sun's core is identical to that determined by current laboratory experiments [e.g. [16][17][18][19][20]. Bahcall [21] has shown that the corrections on the shape of neutrino energy spectra caused by the surrounding plasma in the Sun's core are negligible. For example, the corrections related with the thermal motions of the colliding ions are negligible, as the thermal velocity of ions is much smaller than the velocity of light. Given that the 8 Bν e spectrum shape is well-known and we know that the solar plasma does not influence significantly the nuclear reactions occurring in the Sun's core, the changes detected in the 8 Bν e spectrum will be mostly due to neutrino flavor oscillations induced by matter [a process also known as the Mikheyev-Smirnov-Wolfenstein: MSW effect, [22][23][24]. These modifications will change not only the overall neutrino flux but, more significantly, modify the neutrino spectrum by affecting in a differential manner the survival probability of electronneutrinos -depending upon the energy of the emitted neutrino.
The impact of neutrino flavor oscillations on the total neutrino fluxes is extensively documented in the literature [see 2, and references therein], but the impact on solar neutrino spectra has been discussed only briefly. The reason is that neutrino flavor oscillations are expected to be unimportant, because the temperature of the Sun's core is strongly constrained by the total neutrino fluxes. Nevertheless, as we discuss in this article, this is not necessarily the case, mostly because neutrino flavor oscillations due to the MSW effect in the Sun's core depend strongly of the local properties of the solar plasma. Unlike total neutrino fluxes, these give differential information about the physics of the Sun's core. In particular, if light DM is present in the solar core, the amount of electron-neutrinos converted to other flavors will be different from the value found in the standard solar model [SSM, e.g., 25,26], and consequently their total neutrino fluxes and neutrino spectra will also be different from the SSM. In recent years, significant improvements in the measurement accuracy of solar neutrino fluxes have been instrumental in allowing the use of the Sun to set constraints on the properties of dark matter, including the neutralino [e.g., 27,28]. and impose limits to the expected neutrino fluxes coming from the Sun due to DM annihilation [e.g., 29]. Moreover, a large number of different types of asymmetric DM have been discussed in the literature [30][31][32][33][34]. The presence of dark matter in the Sun's core could help solve the long-running solar composition problem [35], a discrepancy between the solar structure inferred from helioseismology and the one computed from a SSM by inputting the most up-to-date photospheric abundances [36,37]. This type of diagnostic has also been successfully extended to other stars, including other sun-like stars [e.g., 38,39] and neutron stars [40][41][42][43]. In addition, such types of studies have also been extended to the first generation of stars [44][45][46][47][48].
In this paper, we show that by measuring the 8 Bν e solar spectrum, it is possible to constrain the DM content in the Sun's core. This diagnostic complements the total neutrino flux analysis. This is a robust result, as the shape variation of the 8 Bν e spectrum is uniquely related to the radial variation of the plasma properties in the Sun's core, where the maximum accumulation of DM is expected to occur. This type of diagnostic is particularly useful for testing new types of DM models [e.g., 49,50], which have a more pronounced impact in the core of the Sun.

II. CURRENT STATUS OF DARK MATTER RESEARCH
In the last few years, several types of light DM particles have been suggested as an ideal DM candidate for a elementary DM particle, motivated by fundamental theoretical arguments in cosmology and particle physics, and by a few positive hints from some direct DM search experiments. Nevertheless, these results are controversial since other experimental detectors have excluded the same DM parameter space.
In favor of the theoretical argument, several DM models succeed in explaining the observed DM relic density [51,52], as a new type of light DM, usually referred to as asymmetric DM [43,53,54]. Unlike symmetric DM, this new type of DM is believed to be produced in the primordial universe by physical mechanisms identical to the production of baryons, known as darkogenesis [e.g., [55][56][57][58], and likewise composed of an unbalanced mixture of particles and antiparticles. The proportionality of DM particles relatively to DM antiparticles is measured by the asymmetry parameter η DM , which is identical to η B for baryons. In the case that the DM is symmetric, i.e., the DM particle is its own antiparticle, η DM = 0.
The production of asymmetric DM in the early universe is computed by a similar procedure to baryogenesis [e.g. 59,60]. The origin of such a DM asymmetry is not known; however, as suggested by some extensions to the standard model of particle physics, this could be related to the existence of electric and magnetic dipole moments of some standard particles and possibly new particles [e.g., 61]. For reference, the current upper limit on the electric dipole moment of the electron is set to 8.7 × 10 −29 e cm [62]. Similarly Lopes et al. [49] suggest that the magnetic dipole moment of light DM particle could not be larger than 1.6 × 10 −17 e cm.
On the experimental side, the findings of DM searches are a puzzle that is difficult to resolve [63]: several experimental collaborations in direct DM searches have found experimental hints that could be related with dark matter detection: DAMA/LIBRA [64,65] and possibly Co-GeNT [66,67] observed an annual modulation, CRESST [68] and CDMSSi [69,70] show hints of an excess of events. These experiments seem to indicate the existence of a DM candidate with a mass of ∼ 10 GeV and a scattering cross section on hydrogen and other chemical elements varying between 10 −41 and 10 −36 cm 2 . The specific value of the scattering cross section is strongly dependent on the DM model used to interpret the data [e.g., 67,71]. Presently, the null results constraints are from CDMSGe [72], XENON [73,74] and LUX [75]. These experiments found no evidence for an interaction of DM with baryons for the cited mass and scattering cross section range, at least in the case of a contact type of the DM-nucleus interaction models. Nonetheless, there are new theoretical proposals that resolve the differences between the different experimental results, the most successfully being the long-range DM-nucleus interactions. In these type of DM models, the interaction between DM and baryons is not contactlike, but occurs through a light particle mediator [e.g., [76][77][78][79]. The impact of such a DM particle in the Sun's interior can modify significantly its core structure [36].

III. DARK MATTER AND THE SUN
As is usually done in these studies, we consider that the Sun's evolution in a DM halo is identical to the SSM. Likewise, these solar models are required to reproduce the current Sun observables such as radius and luminosity. Therefore, the models to compute the impact of DM in the evolution of the Sun were obtained as follows: for each set of DM parameters, we compute a solar-calibrated model following the same procedure used to compute a SSM [80], i.e., by automatically adjusting the helium abundance and the convection mixing length parameter until the total luminosity and the solar radius are within 10 −5 of the present solar values. Typically, a calibrated DM solar model is obtained after a sequence of 10 to 20 intermediate models.
As the accretion of DM by the star produces minor differences in the Sun's core structure and almost no effect in the stellar envelope, these solar models follow the same path as the SSM in the Hertzsprung-Russell diagram. For the solar model of reference, we choose a SSM with a low-metallicity composition [81], usually referred to as low-Z metallicity SSM. This SSM was computed using an updated version of the stellar evolution code [82]. The code has up-to-date microscopic physics, and in particular uses the nuclear reaction rates from the NACRE Compilation [83]. This SSM predicts solar neutrino fluxes and helioseismic data that are consistent with other SSM models found in the literature [e.g., 26,84]. In relation to the properties of our standard solar model, this can be found in Lopes & Silk [85].
In a DM halo, a star captures DM from the beginning of the premain sequence until the present age (4.6 Gyr). The efficiency of the star in accumulating DM in its core is regulated by three leading processes: capture, annihilation and evaporation of DM particles.
The total number of particles N χ (t) that accumulates inside the Sun at a certain epoch is computed by solving the following differential equation where Γ c , Γ a and Γ e are the capture, annihilation and evaporation rates. A detailed account about these quantities can be found in Jungman et al. [86] and Bertone et al. [87]: -Γ c determines the amount of DM particles captured by the star. This quantity, among others, depends on the radius and escape velocity at each step of the star's evolution. Nevertheless, it is the scattering of the DM particles with baryons which is the leading process in the capture rate. The scattering cross section depends on the mass and spin of the baryon nuclei. As usual, the scattering cross sections of DM particles with nuclei σ χ can be either a spin-dependent or spin-independent cross section, that are represented by σ χ,SD and σ χ,SI . For all of the chemical elements excluding hydrogen, the interaction with a DM particle is of spin-independent type (coherent scattering), for which the scattering cross section scales as the fourth power of the baryon nucleus mass number [e.g., 87]. For hydrogen, the spin-dependent interaction (incoherent scattering) is also taken into account. In our code the Γ c expression is computed following the original expression of Gould [88], Gondolo et al. [89] as described in Lopes et al. [90].
-Γ a depends on the annihilation cross section σv χ of particles and antiparticles. In the current sets of DM models, we are uniquely concerned with σv χ ≈ 0 as discussed in the previous section. A detailed account about the differences between the s-wave and p-wave DM annihilation channels can be found in Lopes & Silk [33].
-Γ e determines the amount of particles that evaporates from the Sun. In our study we use an approximate expression computed for sun-like stars by Busoni et al. [91] from the original work of Griest & Seckel [92]. Nevertheless, this should not much affect our result as we restrict our analysis to DM particles with a mass above 4 GeV, for which evaporation is not significant [93].
In this study, the focus is on the interaction of DM with chemical elements heavier than hydrogen. The impact related to the capture of DM by the scattering off hydrogen was previously studied by Frandsen & Sarkar [31], among others.
Our DM models, if not stated otherwise, have the following properties: the DM particles in the halo follow a Maxwell-Boltzmann velocity distribution, with a thermal velocity v th = 270 km/s; the density of the DM halo is equal to 0.38 GeVcm −3 [e.g., 94]; and the stellar velocity of the Sun is v ⋆ = 220 km/s. The mass of the DM particle m χ , and the spin-independent and dependent scattering cross sections with baryons σ χ,SD and σ χ,SI were chosen to be in agreement with the current experimental bounds for light DM particles. In particular, the spin-dependent and scattering cross section σ χ,SD is equal to 10 −46 cm 2 . It is worth noting that, unlike in previous studies, we solve numerically equation (1), the equation that regulates the accumulation of DM inside the star [49].
The DM impact on the star at each stage of evolution is determined mostly by N χ (t), the number of DM particles accumulated by the star. Once the DM particles are captured by the Sun, these drift towards the Sun's central region, providing the star with a new energy transport mechanism, which then removes energy from the core towards the more external layers of the star. The efficiency of this transport of energy depends mainly on the ratio between the mean free path of the DM particles through the solar plasma l χ , and the characteristic radius of the DM particles distribution in the core of the star r χ [e.g., 48,95]. For most of the DM-nuclei scattering cross sections σ χ (σ SD, χ or σ SI, χ ) considered here, in which l χ ≥ r χ , the energy transport by DM is nonlocal. On the other hand, for large values of DM-nuclei scattering cross sections, in which l χ ≤ r χ , the DM particles are in local thermal equilibrium with the baryons. This latter regime applies only to values of σ χ which are not considered in this work (σ χ ≥ 10 −33 cm 2 ). However, we follow the prescription described in Gould & Raffelt [96] that extends the formalism developed for the local thermal equilibrium to other regimes by the use of tabulated suppression factors. Moreover, once the characteristic radius of the DM core decreases with the mass of the DM particle, such as r χ ∝ m −1/2 χ [e.g., 95]: stellar models computed for DM particles with different masses will produce the 8 B solar neutrino spectra with different shapes.
The main effect of this additional transport of energy is a decrease of temperature in the core of the Sun in relation to the standard solar model. This temperature variation is followed by an increase in the radial profiles of the density ρ(r), and the mean molecular weight per electron µ e (r). But as the increase of the density dominates over the increase of mean molecular weight per electron, and the electron density n e (r) is proportional to the ratio ρ(r)/µ e (r), this leads to an overall increase of n e (r) at core of the star [e.g., 97]. Moreover, as the MSW effect (i.e., the conversion of electron-neutrinos to other neutrino flavors) increases with n e (r), this process leads to a decrease of survival probability of electron-neutrino, as it will be discussed in Sec. V. Furthermore, as the proton-proton chain and carbon-nitrogen-oxygen cycle of nuclear reactions are much more sensitive to the local variations of temperature than density, for a nuclear reaction such as the 8 B β-decay process in the 8 B nuclear reaction rate, this temperature reduction necessarily leads to smaller 8 B solar neutrino flux.

IV. 8 B SOLAR ELECTRON-NEUTRINO SPECTRUM AND FLAVOR OSCILLATIONS
The 8 Bν e spectrum emitted by the 8 B reaction in the Sun's core has been shown to be equivalent to several experimental determinations of the 8 . This calculation used an up-to-date SSM (see text).
shown that the 8 Bν e neutrino spectrum emitted in the Sun's core is equal to the spectrum measured in the laboratory, as the surrounding solar plasma does not affect this type of nuclear reaction. Moreover, the 8 Bν e experimental spectrum agrees remarkably well with the theoretical prediction for neutrinos with an energy below 12 MeV. In particular, the 8 Bν e neutrino spectra deduced from four laboratory experiments [16,[18][19][20]100] agrees within about 1% at high neutrino energies, whereas before they differed by 4% [18]. Figure 1 shows Ψ e ⊙ (E ν ), the 8 Bν e spectrum emitted by the 8 B solar reaction in the Sun's core, with Ψ e ⊕ (E ν ) and Ψ µτ ⊕ (E ν ), the two components of the 8 Bν e neutrino spectrum measured on Earth. We note that the Ψ e ⊙ (E ν ) spectrum (cf. Fig. 1) is identical to the 8 Bν e spectrum measured in the laboratory. Therefore, the only variation expected in the electron-neutrino spectrum measured by solar neutrino detectors, i.e., Ψ e ⊕ (E ν ), is uniquely related to the neutrino flavor oscillations.
The fraction of electron-neutrinos that changes flavor depends on the parameters associated with vacuum and matter oscillations, and this latter process depends also on the local properties of the solar plasma [22]. This is the reason why Ψ e ⊕ (E ν ) is significantly different from Ψ e ⊙ (E ν ). These quantities are related as follows: where P ν e (E ν ) is the electron-neutrino survival probability of a neutrino of energy E ν . P ν e (E ν ) reads  15 Oν e (magenta curve), 17 Fν e (cyan curve). The reference curve (red dashed curve) defines the survival probability of electron-neutrinos in the center of the Sun. The generic properties of such curves can be found in Lopes [22].
where Φ ν (r) is the 8 B electron-neutrino emission source. As usual, r is the solar radius, ρ(r) is the density and A is a normalization constant. In the absence of matter-induced oscillations due to the Earth's atmosphere, P ν e (E ν , r) corresponds to the electron-neutrino survival probabilities on Earth during the day. It follows that P ν e (E ν , r) = cos 4 θ 13 P 2ν e (E ν , r) + sin 4 θ 13 , where P 2ν e (E ν , r) is the probability of a two-flavor neutrino oscillation model [e.g., [101][102][103] and θ 13 a neutrino mixing angle in vacuum. P 2ν e (E ν , r) is given by where ∆m 12 is the mass difference between two flavors, θ 21 is a flavor mixing angle in vacuum and θ m is the matter mixing angle inside the Sun. θ m reads where V m (E ν , r) = 2 √ 2G f n e (r) E ν cos 2 (θ 13 )/∆m 21 , G f is the Fermi constant and n e (r) is the electron density of the solar plasma. Equations (3)(4)(5) determine the probability of electronneutrinos to be converted to other flavors when propagating in matter. This process can affect all solar neutrino sources, but it is more pronounced on the 8 Bν e spectrum. Figure 2 illustrates this specific point. In the figure it is shown the "theoretical" (dashed black curves) and the "observable" (colored curves) survival probability P ν e (E ν ) of electron-neutrinos 1 as a function of neutrino energy for the SSM. P ν e (E ν ) was computed for 8 Bν e , as well as for other neutrino source reactions of the proton-proton chain and carbon-nitrogen-oxygen cycle. Although all neutrino's nuclear reaction sources occur in the Sun's core, the only ones that can be affected the structure changes due to accretion DM in the Sun's core are the ones that produce the neutrinos with the higher energy. This corresponds to the 8 B neutrinos (blue curve) and hep neutrinos (yellow curve), as shown in Figure 2. Nevertheless, the former occur near the center of the Sun and are measured with much better precision that the hep neutrino spectrum. Therefore, the 8 Bν e spectrum will be the most affected by the presence of the DM in the Sun's core.  ) with the fractional solar radius r. The cyan and blue curves correspond to Φ ν (r) and E c ν (r) for the SSM, and the magenta and red curves are the equivalent ones for a DM solar model with the m χ = 5 GeV and σ χ,SI = 10 −36 cm 2 . Both sets of curves have relatively similar shapes. Cyan and magenta curves: Φ ν 's is drawn as a function of r such that Φ ν (r) = F −1 df (r)/dr for which f (r) is the 8 B neutrino flux in s −1 and F is the total neutrino flux for 8 B nuclear reaction rate. Blue and red curves: E c ν 's represented as a function of r corresponds to the minimum neutrino energy E ν that a neutrino must have in order to experience a resonance (see text). Accordingly, electron neutrinos such as E ν ≥ E c ν (r) will experience matter flavor oscillations in the solar core, otherwise this effect is negligible.
The impact of DM on the P ν e (E ν ) or 8 Bν e spectrum can be described as follows: In the Sun's interior, a neutrino of energy E ν can be converted to other flavors if E ν ≥ E c ν (r). The quantity E c ν (r) defines the minimum (critical) energy that a neutrino must have to be strongly affected by flavor oscillations. E c ν (r) is determined by the condition V m (E c ν , r) = cos (2θ 12 ) (from Eq. 5), it follows that E c ν (r) = ∆m 21 /(2 √ 2G f ) cos (2θ 12 )/cos 2 (θ 13 ) n −1 e (r).
1 The "theoretical" P νe (E ν ) although not directly related with the real solar spectrum unlike the "observable" P νe (E ν ) . This quantity illustrates well the effect that the energy dependence of neutrino matter oscillations have on the flux of electron-neutrinos. The survival probabilities of electron-neutrinos and E c ν (r) were computed by using the fundamental parameters of solar neutrino oscillations in the vacuum: ∆m 12 and θ 12 as determined by the KamLAND experiment [102]. Although the contribution related to θ 13 is very small, we take its contribution into account by choosing θ 13 = 9 deg, a value that is in agreement with current experimental measurements [104,105]. Figure 3 shows the critical value E c ν (r) for current SSM and other solar models: neutrinos experiment MSW flavor oscillations in regions of the Sun's core where the neutrino energy is such that E ν ≥ E c ν (yellow region in Fig. 3), otherwise the effect is insignificant. The magnitude of flavor oscillations caused by matter depends on the local value of n e (r), namely the values of density and metallicity. These oscillations are only significant in the Sun's core and negligible in most of the radiative region and solar convection zone. The fraction of electron-neutrinos converted to other flavors depends also on the location of the neutrino source, as well as the local temperature as shown in Fig. 3.

V. DARK MATTER SIGNATURE ON 8 B NEUTRINO SPECTRUM
The presence of DM in the Sun's core changes its thermodynamic structure, modifying the temperature, density and chemical composition, as well as n e (r). Although the effect is relatively small, as neutrinos are very sensitive to the temperature of the Sun's core, minor variations in temperature produce variations in the 8 Bν e spectrum. Consequently, the 8  source (cf. Fig. 3). In addition, the variation of n e (r) distorts the 8 Bν e spectrum, due to an alteration of the survival probability of electron-neutrinos which determines the fraction of electron-neutrinos converted to other flavors. Different DM models have different critical neutrino energies E c ν , leading to distinct P ν e (E ν ) for 8 Bν e and other neutrino sources (cf. Figs 2, 3 and 4). The combination of these different physical processes modifies the shape of Ψ e ⊕ (E ν ), i.e., the 8 Bν e spectrum measured in terrestrial detectors.
The 8 Bν e spectrum is strongly dependent on the temperature, but also on the density and chemical composition. Actually, Ψ e ⊕ (E ν ) the 8 Bν e spectrum shape of the electron neutrino, is related to the variation of the density by three possibilities: the production rate of electron-neutrinos leading to the neutrino function Φ ν (r), the location of the maximum of Φ ν (r) and the survival probability of electron neutrinos (i.e., the conversion of electron-neutrinos to other flavors, Eq. 3): -The first two effects result from the fact that variations in total neutrino flux φ (or equally on the production rate of 8 B neutrinos) depend on the temperature T and density ρ as ∆φ/φ ≈ ∆ρ/ρ + α∆T /T where α = 24.5 is obtained from Turck-Chieze & Lopes [80]. Accordingly, a 10% variation on ∆φ/φ is either attributed to a variation in 10% of density, 0.4% in temperature, or a combination of both. Moreover, a similar variation on the molecular weight is also expected. Nevertheless, as mentioned by several authors [e.,g. 27,30,31] as solar models are calibrated to have the observed solar radius and luminosity, the effective variation of ρ is smaller than the previous estimate, the mitigation coming from the temperature and chemical composition readjustment. The variation of Φ ν (r) leads to a slight change in the location of the maximum of Φ ν (r). This variation also influences the amount of electron-neutrinos converted to other flavors as described by equation (3).
-Equally from Eq. (4), the variation of ∆P e /P e is propor-tional to the variation of electronic density ∆n e /n e (or density and molecular weight). The impact of DM on the electronicneutrino survival probability function P ν e (E ν ) is shown in Fig. 4. As pointed out by previous authors [e.g., 22, and references therein] the effect on P ν e (E ν ) at first order is relatively small, as at low energies the neutrino oscillations are vacuumrelated and therefore insensitive to the Sun's structure; for the higher energy neutrinos, the flavor oscillations are vacuum and density-related (see Fig. 3). The effect of the Sun's structure on P ν e (E ν ) is more pronounced for neutrinos with intermediate energies (from 0.1 to 1.0 MeV). As shown in Fig. 4, the variation of electronic density (density and molecular weight) with solar radius slightly changes the profile of P ν e (E ν ) , leading to small changes in the shape of the 8 B ν e spectrum (see Sec. III). Moreover, Fig. 4 shows the part of the 8 Bν e spectrum that is more affected. This corresponds to neutrinos with an energy in the range: 1 to 10 MeV. This variation is more pronounced for light DM particles with the largest scattering cross sections. This effect reduces the 8 Bν e electron survival probability curve by as much as 6% in relation to the standard case. It is important to observe that such effect on the electron survival probability will distort the electron-neutrino 8 B spectrum in the same neutrino energy range. Although this shape deformation is small, once future measurements of 8 B electron neutrinos will be able to possibly detect such types of effect, if observed it could provide a hint of the existence of dark matter. It is worth remembering that the shape of the 8 B neutrino spectrum is very well measured by current laboratory experiments (see introduction and references therein). For DM solar models discussed in this paper, the maximum effect observed in Ψ e ⊕ (E ν ) uniquely related with P ν e (E ν ) is of the order of 6.5% and occurs near 6 MeV.
The identification by a future solar neutrino detector of a strong distortion in the shape of Ψ e ⊕ (E ν ) compared to that predicted by the SSM, would constitute a strong hint for the presence of DM in the Sun's core. The magnitude of the distortion should give some indication about the amount of DM and the extension of the DM core. Figure 5 shows the difference for the 8 Bν e spectrum for a DM solar model. This is due to the fact that ν e neutrinos of different energy have a different sensitivity to the local distribution of electron density of the Sun's core, specifically, only the more energetic neutrinos are affected by matter flavor oscillations.
In this study we have explored how the presence of DM in the Sun's core changes the shape of solar neutrino spectra, for instance the 8 Bν e neutrino spectrum. In many cases, the impact of DM in the Sun's core can be determined by variations on the total neutrino fluxes due to local temperature changes. Nonetheless, there is an important point to make: even for an identical percentage variation Ψ e ⊕ (E ν ) and φ, there is a fundamental different between both quantities, as the former gives the location where the DM effect occurs (cf. Fig. 5). Indeed, the presence of DM in the solar core will distort Ψ e ⊕ (E ν ) very likely around E ν ∼ 6 MeV (although depending on the DM models, Ψ e ⊕ (E ν ) could be quite singular) and not uniformly distributed, information that is not possible to obtain from φ. As φ is an integrated quantity Φ ν (r), this only give us very limited information about the radial distribution of the DM in the core.
Although this work is focused on studying the 8 Bν e spectrum, our study is easily extended to other neutrino sources as shown in Fig. 2. However, there are currently no neutrino experiments planned to measure the spectrum of other solar neutrino sources in the near future.

VI. SUMMARY
In this paper, we have shown that a detailed measurement of the 8 Bν e spectrum in the range from 1-15 MeV by future solar neutrino experiments will permit us to probe in great detail the core of the Sun (below 0.1R ⊙ ) in a search for traces of DM. We have also shown that this type of DM diagnostic can be extended to other solar neutrino spectra, once the experimental data becomes available (cf. Fig. 2).
The SSM predicts that the 8 Bν e earth spectrum, expected to be measured by solar neutrino detectors, is very different from the 8 Bν e solar emitted spectrum, due to vacuum and matter oscillations which neutrinos experience when travelling to Earth. The presence of DM in the Sun's core will change the magnitude and shape of the 8 Bν e spectrum for terrestrial observers in a very distinct manner. Since there are many astrophysics processes that are not yet included in the standard solar model, that could also affect the physics of Sun's core and the solar neutrino fluxes [e.g., 106], the distortion of the 8 Bν spectrum is an additional important signature that could play a determinant role in disentangling the impact of DM from other possible physical processes. The next generation of solar neutrino detectors like JUNO [6] and LENA [3] should be able to achieve the required precision to test such solar DM models. This will be achieved by simultaneously increasing the preci-sion on the measurements of the solar neutrino spectrum (or the survival electron-neutrino probability) and also by increasing the energy resolution, without which it is not possible to precisely measure the shape distortion of the 8 B ν e spectrum. Moreover, it is expected that the LENA detector after only 5 years of measurements will be able to obtain a probability survival for electron-neutrino (or the equivalent 8 Bν flux) with an experimental error smaller than 0.025 [107]. This precision is sufficient to put constrains in some solar DM models [108], since for some of them the survival electron-neutrino probability variation is of the order of 0.06 (cf. Fig. 4). For instance, if we assume that such experimental accuracy is attained on LENA measurements, solar DM models with m χ ≤ 5GeV and σ χ,SI ≥ 10 −35 cm 2 as the ones shown in Figure 4 can be excluded using the putative LENA data set. This type of diagnostic could help to determine if light asymmetric DM is indeed present in the Sun's core, as this type of DM has been suggested as a nonstandard solution to resolve the solar abundance problem.