The relevance of nuclear reactions for Standard Solar Models construction

The fundamental processes by which nuclear energy is generated in the Sun have been known for many years. However, continuous progress in areas such as neutrino experiments, stellar spectroscopy and helioseismic data and techniques requires ever more accurate and precise determination of nuclear reaction cross sections, a fundamental physical input for solar models. In this work, we review the current status of (standard) solar models and present a detailed discussion on the relevance of nuclear reactions for detailed predictions of solar properties. In addition, we also provide an analytical model that helps understanding the relation between nuclear cross sections, neutrino fluxes and the possibility they offer for determining physical characteristics of the solar interior. The latter is of particular relevance in the context of the conundrum posed by the solar composition, the solar abundance problem, and in the light of the first ever direct detection of solar CN neutrinos recently obtained by the Borexino collaboration. Finally, we present a short list of wishes about the precision with which nuclear reaction rates should be determined to allow for further progress in our understanding of the Sun.


INTRODUCTION
The history of solar models, or standard solar models (SSMs) to be more precise, is formed by three large chapters related to the type of observational and experimental data about the solar interior that existed at any given time. The first part of this history comprises the period over which only neutrino data were available, and it spans about 20 years, from the mid 60s to the early 80s of the past century. During that period, the solar neutrino problem was seen by many as having an origin in the complexities involved in building accurate and precise SSMs, a fundamental part of which is determined by the nuclear reaction rates involved in the generation of the solar nuclear energy. Around the end of that era, the precision of nuclear reaction rates involved in the chains of reactions leading to the production of the different solar neutrino fluxes were on the order to 20 to 30%. These uncertainties may seem large for present day standards. However, if some faith was put in their accuracy, these uncertainties were small enough that associating the solar neutrino problem to nuclear cross section measurements was highly unlikely (Bahcall et al., 1982).
In the mid 80s helioseismology, the study of solar oscillations, evolved into a precision branch of solar physics. The sensitivity of the frequency spectrum of these global pressure waves to the details of the interior solar structure allowed their reconstruction by means of inversion methods, see e.g. Deubner and Gough (1984); Christensen-Dalsgaard et al. (1985), in particular of the solar interior sound speed. This (r)evolution peaked during the second half of the 1990s with the establishment of the Global Oscillation Network Group (GONG), a network of six instruments established around the world that carried out resolved radial velocity measurements of the solar surface (Harvey et al., 1996) and with the launch of the SoHO satellite, both of which provided rich helioseismic datasets. In turn, this led to determination of the solar interior properties with precision of better than 1% (and in some cases even an order of magnitude better) (Gough et al., 1996). These results led to the appearance of a new generation of SSMs (Bahcall et al., 1995;Christensen-Dalsgaard et al., 1996), which were successful in satisfying the tight observational constraints imposed by helioseismology, leaving little room for an astrophysical solution to the solar neutrino problem, as had originally been suggested a few years earlier (Elsworth et al., 1990). Simultaneously, Super-Kamiokande (Fukuda et al., 1998(Fukuda et al., , 2001 led to the precise measurement of 8 B neutrino flux which, in combination with the results of radiochemical experiments Homestake (Cleveland et al., 1998), Gallex (Hampel et al., 1999) and SAGE (Abdurashitov et al., 1999) strongly hinted at the existence of solar neutrino oscillations, result confirmed just a few years later by SNO results (Ahmad et al., 2001(Ahmad et al., , 2002. The needs of refined nuclear reaction rates imposed by the type and quality of the new observational and experimental data led to famous revisions of nuclear reaction rates such as NACRE (Angulo et al., 1999) and in particular that of Solar Fusion I (Adelberger et al., 1998). In the latter, a critical analysis of the accumulated experimental and theoretical data was performed and consensus values were provided for all relevant nuclear reactions affecting energy generation and neutrino production in the Sun. The improvement in the uncertainties, in particular, was about a factor of to 2, leading to typical errors around 10%. Simultaneously, several authors used helioseismic inversion of the solar sound speed to determine, or at least set constraints, on the proton-proton reaction rate, showing that its value had to be within about 15% of its theoretically determined value (Degl'Innocenti et al., 1998;Schlattl et al., 1999;Antia and Chitre, 1999;Turck-chièze et al., 2001;Antia and Chitre, 2002).
The combination of helioseismic constraints and the discovery of neutrino oscillations changed the focus of interest of SSMs. In particular, the accurate and precise determination of neutrino fluxes from individual reactions started playing a fundamental role in the determination of the neutrino oscillation parameters. SSMs became a fundamental source of information, a reference, not just for astrophysics, but for particle physics as well. In 2007, the final and present chapter in this history started when Borexino presented the first measurement of the 7 Be neutrinos (Arpesella et al., 2008), originating from a subdominant branch of reactions, the so-called pp-II branch of the pp-chain that accounts for about 10% of the energy generation of the Sun. Further work by Borexino led to an almost complete characterization of the spectrum of neutrinos from the pp-chain (Agostini et al., 2018). Together with the very precise measurement of the 8 B flux from SNO (Aharmim et al., 2013) and Super-Kamiokande (Abe et al., 2016), we have come full circle and results from solar neutrino experiments can now be used to learn about the properties of the Sun. This is timely. There is a lingering dispute about which is the detailed chemical composition of the Sun, the solar abundance problem (Sect. 2.1), that is intimately linked to the uncertainties in our knowledge of radiative opacities in the solar interior. Solar neutrino data can in principle be used to disentangle this problem (Haxton and Serenelli, 2008;Serenelli et al., 2013;Villante et al., 2014), in particular if the promising results by Borexino on solar CN neutrinos (Agostini et al., 2020a) can be further improved. But progress along this line depends crucially on the accuracy and precision with which nuclear reaction rates are known. The latest compilation, Solar Fusion II (Adelberger et al., 2011), and subsequent work on specific reactions (Section 2), show on average a factor of two improvement with respect to the status 10 to 15 years ago, and 5% uncertainties are nowadays typical. But further work is still needed; uncertainties from nuclear reactions still have a non negligible role in the overall SSMs error budget.
In Section 2 we summarize the current status of SSMs, review the solar abundance problem, the SSM predictions on the solar neutrino spectrum and the status of nuclear reaction rates affecting model predictions. Section 3 presents an analytical formulation of the relation between nuclear reaction rates and solar model properties both for reactions from the pp-chains and CNO-cycles. Section 4 reviews results from numerical SSM calculations, including a detailed assessment of uncertainties and highlighting where progress is most needed, and revises the possibility of using future CN neutrino measurements to determine the solar core C+N abundance.

STANDARD SOLAR MODELS
SSMs are a snapshot in the evolution of a 1 M star, calibrated to match present-day surface properties of the Sun. Two basic assumptions in SSM calculations are: 1) after the phase of star formation the Sun was chemically homogenized as a result of the fully convective phase during its contraction along the Hayashi track and before nuclear reactions start altering its initial composition and, 2) at all moments during its evolution up to the present solar age τ = 4.57 Gyr mass loss is negligible. The calibration is done by adjusting the mixing length parameter (α MLT ) and the initial helium and metal mass fractions (Y ini and Z ini respectively) in order to satisfy the constraints imposed by the present-day solar luminosity L = 3.8418 × 10 33 erg s −1 , radius R = 6.9598 × 10 10 cm (Bahcall et al., 2006), and surface metal to hydrogen abundance ratio (Z/X) , see sect.2.1. As a result of this procedure, SSM has no free parameters and completely determines the physical properties of the Sun. It can be then validated (or falsified) by other observational constraints, in particular by those provided by solar neutrino fluxes measurements and helioseismic frequencies determinations.
The physics input in the SSM is rather simple and it accounts for: convective and radiative transport of energy, chemical evolution driven by nuclear reactions, microscopic diffusion of elements which comprises different processes but among which gravitational settling dominates. Over more than 25 years, since the modern version of the SSM was established with the inclusion of microscopic diffusion (Bahcall and Pinsonneault, 1992;Christensen-Dalsgaard et al., 1993), the continuous improvement of the constitutive physics has brought about the changes and the evolution of SSMs. In particular, a lot of effort has gone into experimental and theoretical work on nuclear reaction rates. But changes in radiative opacities and the equation of state were also relevant. We take here as a reference the results of recent SSM calculations by Vinyoles et al. (2017), the so-called Barcelona 2016 (B16, for short) SSMs, which are based on the following state of the art ingredients. The equation of state is calculated consistently for each of the compositions used in the solar calibrations by using FreeEOS (Cassisi et al., 2003). Atomic radiative opacities are from the Opacity Project (OP) (Badnell et al., 2005), complemented at low temperatures with molecular opacities from Ferguson et al. (2005). Nuclear reaction rates for the pp-chain and CNO-bicycle, which are described in more details in the following section, are from the Solar Fusion II compilation (Adelberger et al., 2011) with important updates for the rates of p(p, e + ν e )d (Marcucci et al., 2013;Tognelli et al., 2015;Acharya et al., 2016), 7 Be(p, γ) 8 B (Zhang et al., 2015) and 14 N(p, γ) 15 O (Marta et al., 2011) reactions. Microscopic diffusion coefficients are computed as described in Thoul et al. (1994). Convection is treated according to the mixing length theory (Kippenhahn and Weigert, 1990). The atmosphere is grey and modeled according to a Krishna-Swamy T − τ relationship (Krishna Swamy, 1966).

The solar composition problem
The solar surface composition, determined with spectroscopic techniques, is a fundamental input in the construction of SSMs. The development of three dimensional hydrodynamic models of the solar atmosphere, of techniques to study line formation under non-local thermodynamic conditions and the improvement in atomic properties (e.g. transition strengths) have led since 2001 to a complete revision of solar abundances. Table 1 lists the abundances determined by different authors for the most relevant metals in solar modeling: GN93 (Grevesse and Noels, 1993), GS98 (Grevesse and Sauval, 1998), AGSS09 (Asplund et al., 2009), C11 (Caffau et al., 2011) and AGSS15 (Scott et al., 2015b,a;Grevesse et al., 2015). Abundances are given in the standard astronomical scale i = log 10 (n i /n H ) + 12, where n i is the number density of a given atomic species.
Note that only abundances relative to hydrogen can be obtained from spectroscopy because the intensity of spectroscopic lines is measured relative to a continuum that is determined by the hydrogen abundance in the solar atmosphere. The last row in the table gives the total photospheric present-day metal-to-hydrogen ratio (Z/X) and it is the quantity used as observational constraint to construct a solar model. In fact, the solar composition set used in solar models determines not only (Z/X) but also the relative abundances of metals in the models. In this sense, Z ini acts as a normalization factor that, together with Y ini and the relation X ini + Y ini + Z ini = 1, determines completely the initial composition of the model.
There is no complete agreement among authors, and some controversy still remains as to what the best values for the new spectroscopic abundances are. However, there is consensus in that all determinations of the solar metallicity based on the new generation of spectroscopic studies yield a solar metallicity lower than older spectroscopic results (Grevesse and Noels, 1993;Grevesse and Sauval, 1998), in particular for the volatile and most abundant C, N, and O. For refractories elements, like Fe, Si, Mg and S that have important role in solar modeling being important contributors to the radiative opacity, meteorites offer a very valuable alternative method (see e.g. Lodders et al. (2009)) and, in fact, elemental abundances determined from meteorites have been historically more robust than spectroscopic ones.
Considering that uncertainties in element abundances are difficult to quantify, it has become customary to consider two canonical sets of abundances to which we refer to as high metallicity (HZ) and low metallicity (LZ) solar admixtures, see e.g. Serenelli et al. (2011);Vinyoles et al. (2017) as reference assumptions for SSM calculations. These are obtained by using the photospheric (volatiles) + meteoritic (refractories) abundances from GS98 and AGSS09 respectively, and are reported in Table 2. In the last column, we give the fractional differences δz i ≡ z HZ i /z LZ i − 1 where z i ≡ Z i /X is the ratio of the i−element abundance with that of hydrogen, to facilitate comparison among the two admixtures. Even if GS98 abundances are presumably surpassed by the more recent determinations, they are still considered as a valid option to construct solar models because they lead to a temperature stratification that well reproduces the helioseismic constraints.  Table 2. The two canonical HZ and LZ solar mixtures given as i = log 10 (n i /n H ) + 12. The two compilations are obtained by using the photospheric (volatiles) + meteoritic (refractories) abundances from GS98 and AGSS09 respectively, and correpond to the admixture labelled as GS98 and AGSS09met in Vinyoles et al. (2017).
This can be better appreciated by considering Tab.3 and Fig.1 where we compare theoretical predictions of SSMs implementing HZ and LZ surface composition with helioseismic determinations of the surface helium abundance Y s , of the convective envelope depth R CZ and the solar sound speed c (r). We see that solar models implementing the LZ abundances fail to reproduce all helioseismic probes of solar properties. This disagreement constitutes the so-called solar abundance problem (Basu and Antia, 2004;Bahcall et al., 2005a;Delahaye and Pinsonneault, 2006) that has defied a complete solution. All proposed modifications to physical processes in SSMs offer, at best, only partial improvements in some helioseismic probes (e.g. Guzik et al. (2005); Castro et al. (2007); Basu and Antia (2008); Guzik and Mussack (2010); Serenelli et al. (2011)). An alternative possibility is to consider modifications to the physical inputs of SSMs at the level of the constitutive physics, radiative opacities in particular. The effective opacity profile in the solar interior results from the combination of the reigning thermodynamic conditions, including composition, and the atomic opacity calculations at hand. Early works (Bahcall et al., 2005b;Montalban et al., 2004) already suggested that a localized increase in opacities could solve or, at least, alleviate the disagreement of low-Z solar models with helioseismology. Christensen-Dalsgaard et al. (2009) andVillante (2010) have concluded that a tilted increase in radiative opacities, with a few percent increase in the solar core and a larger (15-20%) increase at the base of the convective envelope could lead to low-Z SSMs that would satisfy helioseismic probes equally as well as SSMs based on the older, higher, metallicities.
Recent years have seen a surge of activity in theoretical calculations of atomic radiative opacities. Updated calculations (Badnell et al., 2005) by the Opacity Project have led the way, followed by OPAS (Blancard et al., 2012;Mondet et al., 2015), STAR  and a new version of OPLIB, the opacities from Los Alamos (Colgan et al., 2016). For conditions in solar interiors, all theoretical opacities agree with each other within few %. Interestingly, Bailey et al. (2015) have presented the first ever measurement of opacity under conditions very close to those at the bottom of the solar convective envelope. While the experiment has been carried out only for iron, their conclusion is that all theoretical calculations predict a too low Rosseland mean opacity, at a level of 7 ± 4%, for the temperature and density combinations realized in the experiment. Further experimental work on chromium and nickel opacities (Nagayama et al., 2019) helps to evaluate discrepancies between experimental and theoretical results on iron opacity. Results point towards shortcomings that affect models, particularly in the case of open electronic L-shell configurations such   Table 3. Main characteristics of SSMs with different surface composition (Vinyoles et al., 2017). The observational values for Y s and R CZ are taken from Basu and Antia (2004) and Basu and Antia (1997), respectively. The quantity δc/c = (c − c mod )/c mod is the fractional difference between sound speed helioseismic determination and model prediction.
as is present in iron at the base of the convective envelope. Also, the disagreement between theoretical and measured line shapes for the three elements indicates shortcomings in the theoretical understanding of atomic interaction with the plasma. On the other hand, the results also indicate that the quasicontinuum opacity determined experimentally agrees well with the chromium and nickel experiments, contrary to results from the iron experiment. However, the chromium and nickel experiments were carried out at lower temperatures than those used in the original iron experiment, which suggests that the problem of missing quasicontinuum opacity might have an unknown temperature dependence, or that a systematic error affected the high temperature iron measurements. Moreover,  in a recent theoretical analysis of line broadening modeling in opacity calculations, have found that uncertainties linked to this are larger at the base of the convective envelope than in the core. These arguments suggest that opacity calculations are more accurate in the solar core than in the region around the base of the convective envelope.
To take this into account, opacity uncertainty was modelled in B16-SSM calculations in terms of two parameters, κ a and κ b , that can change both the scale and the temperature dependence of opacity according to δk(T ) = κ a + (κ b /∆) log(T /T c ), where δκ is the fractional opacity variation, ∆ = log(T c /T CZ ), T c = 15.6 × 10 6 K and T CZ = 2.3 × 10 6 K are the temperatures at the solar center and at the bottom of the convective zone, respectively. The parameters κ a and κ b have been treated as independent random variables with mean equal to zero and dispersions σ a = 2% and σ b = 6.7%, corresponding to opacity uncertainty σ in = σ a = 2% at the solar center and σ out = (σ 2 a + σ 2 b ) 1/2 = 7% at the base of the convective region.

Nuclear reactions in the Sun
The overall effect of nuclear reactions in the Sun, as in any other star in hydrogen burning stage, is the conversion: with the production of a fixed amount of energy Q = 4 m p + 2m e − m4 He = 26.7 MeV per synthesized 4 He nucleus. Most of this energy is released in the solar plasma and slowly diffuses toward the solar surface supporting the radiative luminosity of the Sun. A small fraction of it, that depends on the specific channel by which hydrogen burning proceeds, is emitted in neutrinos. According to SSM calculations, the two neutrinos carry away about 0.6 MeV on the average.  The SSM predicts that most of the solar energy (> 99%) is produced by the pp-chain, i.e. the hydrogen fusion reaction chain displayed in the left panel of Fig.2. The pp-chain is is mostly initiated by p(p, e + ν e )d reaction and, to a minor extent, by electron capture reaction p(pe − , ν e )d and has several possible terminations that depend on the specific mechanism by which helium-3 nuclei, which are produced by d(p, γ) 3 He reaction, are converted to heavier elements. In the Sun, the dominant mechanism is 3 He( 3 He, 2p) 4 He that corresponds to the so-called pp-I termination of the pp-chain. Alternatively, helium-3 can undergo 3 He( 4 He, γ) 7 Be reaction with the effect of producing beryllium-7. Depending on the destiny of 7 Be, that can be processed either by the electron capture 7 Be(e − , ν e ) 7 Li or by the (largely sub-dominant) proton capture reaction 7 Be(p, γ) 8 B, one obtains the pp-II or the pp-III terminations of the chain. Finally, a very small amount of helium-4 nuclei is produced by 3 He(p, e + ν e ) 4 He reaction. The relative importance of the different branches of the pp-chain depends primarily on the core temperature of the Sun and on the cross section of specific reactions, as will be discussed in next section. The numbers given in Fig.2 show the branching ratios in the present Sun. According to SSM calculations, the central temperature and density of the present Sun are T c 15.6 × 10 6 K and ρ c 150 g cm −3 and they decrease as a function of the solar radius as it is shown in Fig.3. Most of the solar luminosity is produced in the region r < ∼ 0.2 R that contains about 30% of the total mass of the Sun. In this region we observe a relevant increase (decrease) of the helium-4 (hydrogen) mass fraction Y (X), as a result of hydrogen burning during the Sun lifetime. The helium-3 mass fraction (X 3 ) has a non monotonic behaviour, explained by the fact that 3 He burning time is larger than the age of the Sun for r > ∼ 0.3R and thus helium-3 accumulates proportionally to the efficiency of d(p, γ) 3 He reaction. In the energy-producing core, however, 3 He nuclei are efficiently converted to heavier elements by nuclear processes (mainly by 3 He( 3 He, 2p) 4 He), and the abundance X 3 is equal to the equilibrium value.
An alternative hydrogen burning mechanism is provided by the CNO-bicycle that is displayed in the right panel of Fig.2. The CNO-bicycle uses carbon, nitrogen and oxygen nuclei that are present in the core of the Sun as catalysts for hydrogen fusion. It is composed by two different branches, i.e. the CN-cycle and the NO-cycle, whose relative importance depends on the outcome of proton capture reaction on nitrogen-15. In the Sun, the 15 N(p, α) 12 C channel is largely dominant and so, in practice, the CNO-bicycle is reduced to the CN-cycle with a marginal contribution by the NO-cycle. Note that the CN-cycle conserves the total number of 12 C and 14 N nuclei in the core of the Sun, but alters their distribution as it burns into equilibrium, eventually achieving equilibrium abundances proportional to the inverse of the respective rates, see right panel of Fig.3. The reactions controlling conversion of 12 C and 14 N in the solar core and the approach to equilibrium are 12 C(p, γ) 13 N and 14 N(p, γ) 15 O: these are the next-to-slowest and slowest rates in the CN-cycle, respectively. The temperature above which the 12 C burning time through 12 C(p, γ) 13 N is smaller than the Sun's lifetime is T ∼ 10 7 K. In the SSM, the entire energy-producing core, r < ∼ 0.2R and m < ∼ 0.3M is at temperature larger than this value, so that nearly all of the core's carbon-12 is converted to nitrogen-14. The slower 14 N(p, γ) 15 O reaction determines whether equilibrium is achieved. The 14 N burning time is shorter than the age of the Sun for T > ∼ 1.3 × 10 7 K. Therefore equilibrium for the CN cycle is reached only for R < ∼ 0.1R , corresponding to the central 7% of the Sun by mass. Consequently, over a significant portion of the outer core, 12 C is converted to 14 N, but further reactions are inhibited by the Table 4. Solar neutrino fluxes predicted by SSMs with different surface composition (Vinyoles et al., 2017). Units are: 10 10 (pp), 10 9 ( 7 Be), 10 8 (pep, 13 N, 15 O), 10 6 ( 8 B, 17 F), 10 5 (eN, eO) and 10 3 (hep, eF) cm −2 s −1 .
A very effective tool to investigate nuclear energy generation in the Sun is provided by neutrinos which are necessarily produced along with 4 He nuclei during hydrogen burning, in order to satisfy lepton number conservation. Neutrinos free stream in the solar plasma and reach the Earth in about 8 minutes where they can be detected by solar neutrino experiments. While the total amount of neutrinos produced in the Sun can be easily estimated from the solar luminosity constraint, i.e. the assumption that the luminosity radiated from the surface of the Sun is exactly counterbalanced by the amount of energy produced by hydrogen fusion reactions in the solar core (see e.g. Bahcall (2002);Degl'Innocenti et al. (1997); Vissani (2019) for a detailed discussion), the evaluation of their spectrum requires the knowledge of the individual rates of neutrino producing reactions and thus the construction of a complete solar model. We report in Fig.4 and Tab.4, the SSM predictions for the different components of the solar neutrino flux, named according to the specific reaction by which they are produced (Vinyoles et al., 2017). We also include, for completeness, ecCNO neutrinos, i.e. neutrinos produced by electron capture reaction in the CNO-bicycle (in addition to the "standard" CNO neutrinos produced by β decays of 13 N, 15 O and 17 F) that were originally calculated in Bahcall (1990) and Stonehill et al. (2004) and recently reevaluated in Villante (2015) 1 . The two columns "B16-HZ" and "B16-LZ" reported in Tab.4 are obtained by considering two different options for the solar surface composition, as it discussed in Sect.2.1. During the last few decades, solar neutrino experiments have allowed us to determine with great accuracy most of the components of the solar flux. As an example, 7 Be and 8 B neutrino fluxes are measured with accuracy better than ∼ 3% by Borexino (Agostini et al., 2018), Super-Kamiokande (Abe et al., 2016) and SNO (Aharmim et al., 2013). The pp and pep-neutrino flux can be determined with < ∼ 1% accuracy by assuming the solar luminosity constraint, see e.g. Bergstrom et al. (2016). These fluxes, however, have been also directly measured by Borexino (Bellini et al., 2012(Bellini et al., , 2014Agostini et al., 2018) with ∼ 10% and ∼ 17% accuracy, respectively. Finally, Borexino has recently obtained the experimental identification of CNO neutrinos (Agostini et al., 2020a), providing the first direct evidence that CNO-bicycle is active in the Sun.

Nuclear reaction rates
The cross sections of nuclear reaction in pp-chain and in CNO-bicycle are fundamental inputs for SSM calculations. Even if the focus of this work is on the role of nuclear rates for solar modeling (more than on reviewing the present situation for cross section measurements and calculations), we believe that it is useful to briefly discuss the adopted assumptions for the B16-SSM (Vinyoles et al., 2017). whose results have been previously discussed. The nuclear rates adopted for these models are from the Solar Fusion II compilation (Adelberger et al., 2011) with few relevant changes summarized in the following.
• p(p, e + ν e )d: The astrophysical factor S 11 (E) has been recalculated in Marcucci et al. (2013) by using chiral effective field theory framework, including the P-wave contribution that had been previously neglected. For the leading order they obtain S 11 (0) = (4.03 ± 0.006) · 10 −25 MeV b. More recently, and also using chiral effective field theory, S 11 (E) was calculated by Acharya et al. (2016), resulting in S 11 (0) = 4.047 +0.024 −0.032 · 10 −25 MeV b. This is in very good agreement with result from Marcucci et al. (2013). Acharya et al. (2016) have performed a more thorough assessment of uncertainty sources leading to an estimated error of 0.7%, much closer to the 1% uncertainty which was obtained by Adelberger et al. (2011). In B16-SSM calculations, the astrophysical factor S 11 (E) is taken from Marcucci et al. (2013) with a conservative 1% error estimate (Vinyoles et al., 2017).
• 7 Be(p, γ) 8 B: Solar Fusion II recommended value is S 17 (0) = (2.08 ± 0.07 ± 0.14) · 10 −5 MeV b (Adelberger et al., 2011), where the first error term comes from uncertainties in the different experimental results and the second one from considering different theoretical models employed for the low-energy extrapolation of the rate. Zhang et al. (2015) presented a new low-energy extrapolation S 17 (0) = (2.13 ± 0.07) · 10 −5 MeV b, based on Halo Effective Field Theory, which allows for a continuous parametric evaluation of all low-energy models. Marginalization over the family of continuous parameters then amounts to marginalizing the results over the different lowenergy models. In B16-SSM calculations, it was conservatively adopted an intermediate error between those from Zhang et al. (2015) and Adelberger et al. (2011). The adopted value is S 17 (0) = (2.13 ± 0.1) · 10 −5 MeV b. The derivatives of the astrophysical factor were updated by using the recommended values in Zhang et al. (2015).  Table XI in that work) this leads to S 114 (0) = (1.59 · 10 −3 ) MeV b, about 4% lower than the previous recommended value in Adelberger et al. (2011). The derivatives and the errors remain unchanged.
• 3 He( 4 He, γ) 7 Be: Two recent analyses (deBoer et al., 2014;Iliadis et al., 2016) have provided determinations of the astrophysical factor that differs by about 6% (to be compared with a claimed accuracy equal to 4% and 2% for deBoer et al. (2014) and Iliadis et al. (2016), respectively). Considering that the results from deBoer et al. (2014) and Iliadis et al. (2016) bracket the previously adopted value from Adelberger et al. (2011), the latter was considered as preferred choice in Vinyoles et al. (2017).
Finally, Salpeter's formulation of weak screening (Salpeter, 1954) is adopted. The validity of this formulation for solar conditions, where electrons are only weakly degenerate, has been discussed in detail in Gruzinov and Bahcall (1998), where a more sophisticated approach was shown to lead, to within differences of about 1%, to Salpeter's result. Other proposed deviations from this formulation have been discussed at length in Bahcall et al. (2002), including different approaches to dynamic screening, and shown to be flawed or not well physically motivated. More recent calculations of dynamic screening (Mao et al., 2009;Mussack and Dappen, 2011) still leave, however, some room for discussion on this topic. In the weak screening limit, and in conditions under which screening is not numerically large, the dominant scaling is with the product of the charge of the two reacting nuclei. In the solar core, screening enhancement is about 5% for p(p, e + ν e )d, 20% for 3 He( 4 He, γ) 7 Be and 7 Be(p, γ) 8 B, and 40% for 14 N(p, γ) 15 O.

THE ROLE OF NUCLEAR REACTIONS
In the following, we discuss the role of nuclear reactions in SSM construction. Among nuclear processes, the p(p, e + ν e )d reaction is the only one that can affect the temperature stratification of the Sun. Indeed, this process determines the global efficiency of hydrogen burning in the Sun. The other reactions in the pp-chain and in the CNO-cycle have a minor importance in this respect. However, they have a crucial role in determining the relative rates of the different pp-chain terminations and the efficiency of the CNO-cycle, thus affecting the predictions for the different components of the solar neutrino spectrum.

The pp-reaction rate and the central temperature of the Sun
In SSM calculations, where the Sun is assumed to be in thermal equilibrium, the rate of the pp-reaction is basically determined by the solar luminosity. Indeed, by considering that helium-4 is mainly produced by 3 He( 3 He, 2p) 4 He, we arrive at the conclusion that the integrated pp-rate in the Sun is λ 11 ∼ 2L /Q I , where Q I = Q − 2 E ν pp ∼ 26.2 MeV is the energy released in the solar plasma when 4 He is synthesized through pp-I termination. In the previous expression, we considered that the average energy of neutrinos produced by p(p, e + ν e )d is E ν pp = 0.265 MeV and we took into account that, at equilibrium, the pp-I termination involves twice the pp-reaction in order to feed the process 3 He( 3 He, 2p) 4 He.
Being the reaction rate fixed by the observed luminosity, the cross section of p(p, e + ν e )d determines the central temperature of the Sun, as it is explained in the following. The rate λ 11 can be expressed as: where ρ is the density, m u is the atomic mass unit, X is the hydrogen mass fraction and σv 11 is the reaction rate per particle pair of the p(p, e + ν e )d reaction. The above integral involves, in principle, the entire solar structure but it gets a non-vanishing contribution only from the inner core of the Sun at r ≤ 0.3 R . This can be appreciated by looking at Fig. 5 where we show the differential rates (1/λ ij ) (dλ ij /dr) for the p(p, e + ν e )d (black), p(pe − , ν e )d (blue), 3 He( 4 He, γ) 7 Be (red), 7 Be(p, γ) 8 B (green) and 3 He(p, e + ν e ) 4 He (purple) reactions as a function of the solar radius. The different curves are all normalized to one in order to facilitate comparison among them. These curves also corresponds to the normalized production rates of pp, pep, 7 Be, 8 B and hep neutrinos, respectively 2 .
Taking into account that p(p, e + ν e )d reaction is active in a narrow region of the Sun at r 0 0.1R whose physical conditions are similar to those at the solar center, we write the approximate scaling law: where the notation Q c indicates that the generic quantity Q is evaluated at the center of the Sun, S 11 is the astrophysical factor of the pp-reaction and we considered that σv 11 ∝ S 11 T γ 11 c with γ 11 4. Eq.(3) implies the following linearised relationship: δλ 11 2δρ c + 2δX c + γ 11 δT c + δS 11 , where δQ indicates the fractional variation of the quantity Q with respect to the reference SSM value. The above expression contains input parameters for solar model construction, i.e. the astrophysical factor S 11 , and structural parameters, like e.g. the temperature, density and hydrogen abundance in the core of the Sun which are the result of solar model self-calibrated calculations. In principle, a modification of S 11 induces a change of the solar structure and, thus, the different terms in the r.h.s of eq.(4) are correlated. In order to keep δλ 11 0, an increase of the astrophysical factor δS 11 ≥ 0 has to be counterbalanced by an opposite contribution 2δρ c + 2δX c + γ 11 δT c ≤ 0. This is achieved by varying the initial helium and metal abundance of the Sun according to δY ini = 0.6 δS 11 and δZ ini −0.10 δS 11 with the effect of obtaining a (slightly) colder solar core. We obtain numerically: δT c ∼ −0.13 δS 11 (4) that will be useful in the following to understand the effects of S 11 variations on the various components of the solar neutrino spectrum. In Fig.6, we show the effect of a 10% increase of S 11 on the temperature profile of SSMs and on the helioseismic observable quantities δc(r) and δρ(r).

The dependence of neutrino fluxes on the central temperature of the Sun and on nuclear reaction cross sections
Even a small modification of the central temperature of the Sun reflects into large variations of solar neutrino fluxes. By considering the arguments discussed in Bahcall and Ulmer (1996) and Degl'Innocenti et al. (1997), we discuss the dependence of solar neutrino fluxes on the core temperature of the Sun, highlighting the role of nuclear reactions for determining the branching ratios of the different pp-chain terminations and the efficiency of the CNO-bicycle.

PP-chain neutrino fluxes
The pp-neutrino flux: The vast majority of the solar neutrino emission is due to pp-neutrinos whose flux Φ(pp) is directly linked to λ 11 being Φ(pp) = λ 11 /(4πD 2 ) where D = 1 A.U. is the Sun-Earth distance. According to discussion in the previous Section, the rate λ 11 is directly fixed by solar luminosity and thus Φ(pp) is expected to be independent from the central temperature of the Sun and nuclear reaction cross sections. This result is obtained by assuming that pp-I is the only mechanism for helium-4 production by nuclear reaction in the Sun. A more accurate description can be obtained by taking into account the contribution the secondary branches of the pp-chain (namely, the pp-II termination) initiated by the 3 He( 4 He, γ) 7 Be which provides an alternative 3 He burning mechanism to the most common 3 He( 3 He, 2p) 4 He. In this assumption, we have: where λ 33 and λ 34 are the integrated rate of the 3 He( 3 He, 2p) 4 He and 3 He( 4 He, γ) 7 Be reactions, while Q I = 26.20 MeV and Q II = 25.65 MeV give the amount of energy, corrected for neutrino emission, delivered in the plasma when 4 He is produced through pp-I and pp-II terminations, respectively. By considering that λ 11 = 2λ 33 + λ 34 at equilibrium, we arrive at the conclusion that (Bahcall and Ulmer, 1996): where we considered that Q I Q II . While the first term in the r.h.s. of the above equation is constant, the rate λ 34 depends on the temperature of the plasma and on nuclear reaction cross sections. If we take into account that λ 34 ∝ S 34 · (S 11 /S 33 ) 1/2 · T β Be c with β Be ∼ 11, as motivated later in this section, we obtain the following relationship δΦ(pp) = −η δS 34 − η 2 (δS 11 − δS 33 ) + β pp δT c (7) that gives the fractional variation of the flux δΦ(pp) as a function of fractional variation of the core temperature δT c and of the astrophysical factors δS ij . The coefficients in the above equation correspond to the logarithmic derivatives of Φ(pp) with respect to these quantities and are given by η = λ 34 /λ 11 Φ( 7 Be)/Φ(pp) 0.08 and β pp = −ηβ Be −0.9, showing that the pp-neutrino flux is a decreasing function of the central temperature of the Sun.
The pep-neutrino flux: The pep-neutrinos are produced by electron capture reaction p(pe − , ν e )d which is linked to the β-decay process p(p, e + ν e )d by well-known 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. It can be determined with ∼ 1% precision for the conditions of the solar interior and is mildly dependent on the properties of the solar plasma, being roughly proportional to T −1/2 c n e (see e.g. Adelberger et al. (2011) for a review). We can thus assume Φ(pep) ∝ T 1/2 c Φ(pp), allowing us to conclude: where β pep = β pp − 1/2 −1.4, and we neglected effects related to possible density and chemical composition variations in the solar core.
The 7 Be-neutrino flux: The formation of beryllium-7 through 3 He( 4 He, γ) 7 Be leads to neutrino production through the electron capture reaction 7 Be(e − , ν e ) 7 Li. This process largely dominates over the competing proton-capture reaction whose effects are discussed in the following paragraph. Taking this into account, the Be-neutrino flux can be directly estimated from the rate of the 3 He( 4 He, γ) 7 Be reaction by using Φ( 7 Be) = λ 34 /(4πD 2 ). The rate λ 34 is given by: where Y (X 3 ) is the helium-4 (helium-3) mass fraction and σv 34 is the reaction rate per particle pair of 3 He( 4 He, γ) 7 Be. The amount of helium-4 nuclei in the present Sun is determined by the assumed initial abundance Y ini and by nuclear processes that have converted hydrogen into helium during the Sun evolution. We may thus expect that Y depends on nuclear cross sections, in particular on S 11 that determines the global efficiency of hydrogen burning. This dependence is however marginal because the product L τ essentially provides an observational determination of the integrated solar luminosity (and thus of the total amount of helium synthesized by nuclear reactions during the Sun lifetime). The helium-3 abundance in the solar core depends instead on the temperature T c and on the cross sections of the p(p, e + ν e )d and 3 He( 3 He, 2p) 4 He reactions. It can be indeed calculated by using the equilibrium condition X 3 X 3,eq = 3 X σv 11 2 σv 33 where X is the hydrogen mass fraction. Considering that σv ij ∝ S ij T γ ij c , this can be rewitten as X 3,c ∝ (S 11 /S 33 ) 1/2 · T (γ 11 −γ 33 )/2 c where we neglected effects related to possible hydrogen abundance variations 3 . This expression, combined with Eq. (9), allows us to conclude that: Φ( 7 Be) ∝ S 34 · (S 11 /S 33 ) 1/2 · T β Be c (11) or, equivalently, δΦ( 7 Be) = δS 34 + 1 2 (δS 11 − δS 33 ) + β Be δT c where β Be = γ 34 + (γ 11 − γ 33 )/2 ∼ 11. Note that the 7 Be-neutrino flux does not depends on the cross section of 7 Be(e − , ν e ) 7 Li, due to the fact that (almost) the totality of beryllium-7 nuclei produced by 3 He( 4 He, γ) 7 Be are expected to decay through this reaction.
The 8 B-neutrino flux: The 8 B neutrinos constitute a largely subdominant component of the solar flux which is produced when 7 Be nuclei capture a proton (instead of an electron) producing 8 B (instead of 7 Li). The 8 B-neutrino flux is thus given by Φ( 8 B) = r Φ( 7 Be) where r ≡ λ 17 /λ e7 is the ratio between proton and electron capture rates on beryllium-7. The parameter r scales as r ∝ (S 17 /S e7 ) · T α c where α = γ 17 + (1/2) and we have considered that σv e7 ∝ S e7 T −1/2 c for electron capture reaction. Taking this into account, we obtain the following scaling law: that also corresponds to: with β B = β Be + γ 17 + 1/2 24. The large value of β B indicates that 8 B neutrinos are a very sensitive probe of the core temperature of the Sun.

The CNO neutrino fluxes
The neutrino fluxes produced in the CN-cycle by β-decay (and electron capture reactions) of 13 N and 15 O nuclei, besides depending on the solar central temperature, are approximately proportional to the stellar-core number abundance of CN elements. This dependence is relevant to understand the role of cross section for CNO-neutrino production. Moreover, as it is discussed in Haxton and Serenelli (2008); Haxton et al. (2013); Serenelli et al. (2013), it permits us to use CNO neutrinos, in combination with other neutrino fluxes, to directly probe the chemical composition of the Sun.
The 15 O-neutrino flux: This component of the solar neutrino spectrum is determined by the production rate of oxygen-15 by 14 N(p, γ) 15 O reaction in the core of the Sun. It can be calculated as Φ( 15 O) = λ 114 /(4πD 2 ) where the rate λ 114 , given by: 3 We evaluate the exponents γ ij by using γ ij (E 0 ) ij − 2/3 where (E 0 ) ij is the Gamov peak energy of the considered reaction, see e.g. Bahcall and Ulmer (1996) is proportional to the nitrogen-14 mass fraction X 14 in the solar core (see Fig.3) and to the reaction rate per particle pair σv 114 of the 14 N(p,γ) 15 O reaction. The above integral get a non vanishing contribution from a narrow region at r < ∼ 0.1R whose conditions are similar to that at the solar center, see Fig.5. We thus write the approximate scaling law: where S 114 is the astrophysical factor of the 14 N(p,γ) 15 O reaction, we considered that σv 114 ∝ S 114 T γ 114 c and we defined β O = γ 114 20. Eq. (16) implies the following linearised relationship: In the above expressions, we neglected effect related to possible variations of the density and of the hydrogen abundance in the solar core, since these are expected to be small. We instead explicitly considered the dependence of Φ( 15 O) on the central abundance of nitrogen-14 which is essentially determined, as it is explained in the following, by the total abundances of CN elements in the solar core. It is useful to remark that, being the CNO cycle sub-dominant, a modification of its efficiency does not alter the solar luminosity and does not require a readjustement of the central temperature. Moreover, carbon and nitrogen give a marginal contribution to the opacity of the solar plasma and thus a variation of their abundances do not alter the temperature stratification. As a result of this, we can consider the different terms in eq.(17) as independent.
The 13 N-neutrino flux: The flux of 13 N-neutrinos can be calculated Φ( 13 N) = λ 112 /(4πD 2 ) where λ 112 is the total rate of the 12 C(p, γ) 13 N reaction in the Sun. This is given by: where X 12 is the carbon-12 mass fraction and σv 112 is the reaction rate per particle pair of 12 C(p,γ) 13 N. We can write: where the quantity: gives the contribution to the total rate produced in the region of the Sun where the CN-cycle is incomplete. The above integral vanishes indeed for r ≤ 0.13R where the equilibrium condition for the CN-cycle ensures that (X 12 /12) σv 112 − (X 14 /14) σv 114 = 0. This can be appreciated in the right panel of Fig.5 where we show the differential rate dλ 114 /dr and dλ 112 /dr of 14 N(p, γ) 15 O (black) and 12 C(p, γ) 13 N (red) reactions as a function of the solar radius r.
Eq.(19) implies that Φ( 13 N) can be decomposed as the sum: where the quantity Φ( 13 N) (ne) ≡ λ (ne) 112 /(4πD 2 ) represents the neutrino flux produced in the region 0.13 < ∼ r/R < ∼ 0.25, where 14 N(p, γ) 15 O reaction is not effective. This component of the flux scales as: where we considered that σv 112 ∝ S 112 T γ 112 c with γ 112 18 and we neglected effects related to possible variations of density and hydrogen abundance. Note that the carbon-12 mass fraction in Eq. (22) is evaluated at r ne 0.16 R where the out-of-equilibrium 13 N-neutrino production rate is maximal, see Figs.3 and 5.
In principle, the temperature should be also evaluated at this position. However, we can take the central value T c as representative for the entire energy producing region, motivated by the fact that T (r) (differently from X 12 (r)) is slowly varying in the solar core. Eq.(22) implies the following relationship: that combined with eq.(17) gives: where f = Φ( 15 O)/Φ( 13 N) = 0.74 is the ratio between 15 O and 13 N neutrino fluxes in SSMs (Vinyoles et al., 2017).
The abundance of carbon and nitrogen in the core of the Sun Eqs. (17) and (24) describe the dependence of the CN-neutrino fluxes from the abundances of nitrogen X 14,c and carbon X 12 (r ne ) at the center of the Sun and close to r ne = 0.16R , respectively. These abundances are determined by the formation and chemical evolution history of the Sun, i.e by the initial solar composition and by the subsequent action of nuclear reactions and elemental diffusion, as it is described in the following. Let us first consider that the CN-cycle conserves the total number of CN-nuclei in the core of the Sun. This is shown in Fig.3 by the behaviour of the quantity: N ≡ X 12 /12 + X 13 /13 + X 14 /14 which is proportional to the total carbon+nitrogen number density (X 13 represents the carbon-13 mass abundance) and it is nearly constant in the solar core despite the action of nuclear reactions. In the SSM paradigm, the radial dependence of N is only due to elemental diffusion so that we can write: where N ini is the initial carbon+nitrogen abundance that is assumed to be uniform in the solar structure while the function ∆(r) describes the effects of gravitational settling. It takes the value ∆ c = 0.06 at the center of the Sun that can be considered also representative for r ne = 0.16R , and ∆ s = −0.09 in external convective envelope according to SSM calculations (Vinyoles et al., 2017). It is useful to connect the core composition to photospheric abundances since these are observationally constrained by spectroscopic measurements. We thus write: where N s (N c ) is the carbon+nitrogen abundance in the external convective envelope (at the center) of the Sun while ∆ (cs) = (∆ c − ∆ s )/(1 + ∆ s ) = 0.16 represents the fractional difference between core and surface abundances.
The abundance X 14,c that controls the equilibrium production of CN-neutrinos is directly related to total abundance of carbon and nitrogen in the core of the Sun. Indeed, for r < ∼ 0.1R the CN-cycle is complete and all available carbon is essentially transformed into nitrogen, giving X 14,c 14 N c (see Fig.3). We thus obtain the relation δX 14,c = δN c that, by taking advantage of eqs. (25,27), can be rewritten as: where b = 1/(1 + 0.16) = 0.86, a = 6ξ/(6ξ + 7) 0.20 and ξ = (X N, s /X C, s ) 0.30 is the surface nitrogen-to-carbon ratio in SSM. The first two terms of the r.h.s. in the above equation describe the effects produced by a variation of the surface composition. A modification of the chemical composition profile that is instead produced either "primordially" (e.g. by assuming that the Sun was not born chemical homogenous) or during the evolution (e.g. by anomalous diffusion) on time scales longer than carbon and nitrogen burning time at the solar center, is instead described in terms of a variation of ∆ (cs) from the SSM value, i.e. by assuming ∆ (cs) − 0.016 = 0.
A slightly more involved expression is obtained for the abundance X 12 (r ne ) that controls the nonequilibrium production of 13 N-neutrinos. In the relevant region 0.13 < ∼ r/R < ∼ 0.25, the carbon-12 abundance differs from the surface value X 12, s due to the action of elemental diffusion and 12 C(p, γ) 13 N reaction only, since further reactions are inhibited by the 14 N(p, γ) 15 O bottleneck. It can be approximately described as where the quantity D 112 represents the carbon-12 burning rate averaged over the Sun lifetime, see appendix for details. The maximal neutrino production is achieved at r ne 0.16 R where the integrated burning rate is D 112 (r ne ) t 1. Indeed, in the inner core where D 112 t 1, carbon-12 abundance is too low to efficiently feed 12 C(p, γ) 13 N reaction. On the other hand, the carbon-12 burning time is much larger than solar age (and thus 12 C(p, γ) 13 N reaction is not effective) in more external regions where D 112 (1/t ), as can be understood by considering that D 112 D 112 . Taking this into account, we obtain the following relation: where we considered that D 112 (r ne ) ∝ S 112 T γ 112 c .
The final expressions the CN neutrino fluxes By using the above equations, we are able to calculate the dependence of neutrino fluxes produced in the CN-cycle on the properties of the Sun. By using eqs.  15 and a ≡ f a = 0.15. Note that, in the derivation of the second equation, we took into account that the third and the fourth terms in the r.h.s of eq.(31) cancels the dependence of Φ( 13 N) (ne) on S 112 and T c expressed in eq.(23). This is due to the fact that, as far as the 13 N-neutrino (non equilibrium) production rate is concerned, the effect of 12 C(p, γ) 13 N cross section enhancement is compensated by the reduction of residual carbon-12 abundance due the more efficient carbon burning.

NUMERICAL RESULTS AND NUCLEAR UNCERTAINTIES
The expressions obtained for the neutrino fluxes can be compared with the results of SSMs calculations. In particular, the numerical coefficients in eqs. (7,8,12,14,32) should reproduce the logarithmic derivatives of the neutrino fluxes with respect to the astrophysical factors of the relevant nuclear cross sections reported in Tab. 5. We see that a good agreement exists, indicating that all the major physical effects are included in our discussion and correctly described. In the case of S 11 , we have to take into account that the role of this parameter is twofold; indeed, besides altering the efficiency of pp-reaction (at fixed temperature), this parameter also induces a variation of the central temperature of the Sun as described by Eq.(4). This effect, combined with the strong temperature dependence of the fluxes, allow us to understand the large values for logarithmic derivatives reported in the first column of Tab.5.
For completeness, we also discuss in the last two rows of Tab. 5 the dependence of the helioseismic observable quantities Y s (surface helium abundance) and R CZ (depth of the convective envelope) on nuclear reactions cross sections. We see that S 11 is the only nuclear parameter that affects the predictions for these quantities. The effects of S 11 modifications on sound speed and density profiles are shown in the right panel of Fig.6. Finally, Tab. 6 gives the logarithmic derivatives of neutrino fluxes and helioseismic quantities on other input parameters (beside nuclear cross sections) which are necessary to construct SSMs. These are: the solar age (age), luminosity (lumi) and the diffusion coefficients (diffu); the opacity of solar plasma whose uncertainty is described in terms of two parameters κ a and κ b defined in Sect. 2.1; the surface abundances of key elements (C, N, O, Ne, Mg, Si, S, Ar, Fe) which are determined through spectroscopic measurements as discussed in Sect. 2.1. We can see that the logarithmic derivatives of the CN-neutrino fluxes with respect to the surface carbon and nitrogen abundances are correctly predicted by Eqs. (32).
The uncertainties in solar properties leading to environmental effects and chemical composition parameters, together with uncertainties in nuclear reaction cross sections propagate to SSM predictions which are affected by a theoretical (or model) error that can be estimated by Monte-Carlo techniques and/or linear propagation. By using this approach, the fractional error σ Q on a generic SSM prediction Q can be obtained as the sum (in quadrature) of different contributions, according to: where I = age, lumi, . . . indicates a specific input, σ I represents its fractional uncertainty and α(Q, I) ≡ d ln Q/d ln I is the logarithmic derivative of Q with respect to I. Tab.7 contains the uncertainties σ I that have been considered for the construction of B16-SSMs (the surface composition errors are reported in Tab.2), see Vinyoles et al. (2017) for details. By using these values, one is able to estimate the contribution δQ I ≡ α(Q, I)σ I of each input parameter to the total error budget of Q. The dominant error sources for solar neutrino fluxes and helioseismic quantities are given in Tab. 8. 4 Focusing on nuclear reactions, we note that, despite the progress in the field, they are still an important uncertainty source for neutrino fluxes. In particular, the error contributions from S 34 and S 17 are comparable to or larger than the uncertainties in the experimental determinations of Φ( 8 B) and Φ( 7 Be). As discussed in Vinyoles et al. (2017), the ability of solar neutrinos produced in the pp-chain to play a significant role in constraining physical conditions in the solar interior depends, although it is not the only factor, on pinning down errors of nuclear reaction rates to just ∼ 2%. For CN fluxes, we see that S 114 is the dominant error source if composition is left aside. This is particularly relevant, especially in consideration of the fact that Borexino has just opened the era of CNO neutrino detection, obtaining for the first time ∼ 5σ direct experimental evidence for a non vanishing flux from the Sun (Agostini et al., 2020a).
For a correct evaluation of the importance of nuclear cross section, it should be remarked that, while neutrino fluxes generally change with variation in any of the input parameters, SSM predictions are strongly correlated with a single output parameter, the core temperature T c (Bahcall and Ulmer, 1996;Degl'Innocenti et al., 1997;Haxton and Serenelli, 2008;Serenelli et al., 2013). As a consequence, a multi-dimensional set of variations of enviromental and chemical composition parameters {δI} often collapses to a one-dimensional dependence on δT c , where δT c is an implicit function of the variations {δI}. The dominance of T c as the controlling parameter for neutrino fluxes can be exploited to cancel out uncertainties in the analysis of solar neutrino data. One can indeed form weighted ratios Φ(ν 1 )/Φ(ν 2 ) x 12 , or equivalently weighted fractional differences δΦ(ν 1 ) − x 12 δΦ(ν 1 ) with respect to SSM predictions that are nearly independent of T c and thus marginally affected by environmental effects and chemical composition, using the residual dependence on selected parameters to learn something about them.
In Haxton and Serenelli (2008) and Serenelli et al. (2013), it was suggested to combine the CN-neutrino fluxes with the boron neutrino flux that, due to the exquisite precision of current experimental results and the large temperature sensitivity can be efficiently used as solar thermometer. As can be understood by considering Eqs. (14) and (32) where x = β O /β B 0.8 and x = f x 0.6, that are independent from δT c . This possibility is extremely important because it allows us to cancels out the dependence on the radiative opacity (implicit in δT c ). The uncertainty of available opacity calculations is indeed not easily quantified and may be potentially underestimated. Moreover, it breaks the degeneracy between composition and opacity effects on solar observable properties. Indeed, the considered flux combinations only depend on the carbon and nitrogen abundance in the solar core allowing us to test the chemical composition and evolution of the Sun. The first two terms in the r.h.s. of Eqs. (34,35) quantify the effects of a variation of the surface C and N abundances. A change of the diffusion efficiency is instead described in terms of a variation of ∆ (cs) from the SSM value, i.e. by assuming ∆ (cs) − 0.016 = 0. It should be remarked that the ability to probe solar composition by using this approach is only limited by experimental accuracy of flux determinations and by nuclear cross section uncertainties.
While the above relationships are based on the simplified arguments discussed in the previous section, the optimal combinations δΦ(ν 1 ) − x 12 δΦ(ν 2 ), or equivalently weighted ratios Φ(ν 1 )/Φ(ν 2 ) x 12 , can be determined by using the power-law coefficients from Vinyoles et al. (2017) given in Tab.6. The parameter x 12 is obtained by minimizing the residual where the sum extends to the N input parameters whose dependence we want to cancel out and σ I are the corresponding uncertainties. The minimal value for ρ gives the intrinsic error in the considered approach. This method, originally proposed by Haxton and Serenelli (2008) and Serenelli et al. (2013), has been recently adapted to Borexino in Agostini et al. (2020b). By taking into account that the measured CNO neutrino signal in Borexino is basically probing δφ BX CNO ≡ ξ δΦ( 15 O) + (1 − ξ) δΦ( 13 N) with ξ = 0.764 , it was concluded that the surface composition of the Sun can be probed by the combination: δR BX CNO − 0.716 δΦ( 8 B) = 0.814 δX 12, s + 0.191 δX 14, s ± 0.5% (env) ± 9.1% (nucl) ± 2.8% (diff) where δR BX CNO is the fractional difference of the observed CNO signal with respect to SSM expectations and the quoted uncertainties are obtained by propagating errors of SSM input parameters. The error budget is presently dominated by the uncertainty of the CNO signal Borexino measurement. However, a relevant error (∼ 10%) is also provided by nuclear reactions, with the largest contributions coming from S 114 (7.6%), S 34 (3.4%), and S 17 (3.5%). In the perspective of future improvements of the CNO signal determination, it is evidently important to have reliable and accurate determinations of these cross sections.  Table 5. The logarithmic derivatives α(Q, I) of the solar neutrino fluxes with respect to nuclear input parameters calculated in B16-HZ SSMs.

CONCLUDING REMARKS
A fundamental part in solar model calculations is the knowledge of the rates of nuclear reactions involved in the generation of solar nuclear energy. During the last decades, we experienced a substantial progress in the accuracy of SSM calculations that was made possible, among the other ingredients, by the continuous improvements of nuclear cross sections that are now typically determined with ∼ 5% accuracy. However, SSMs have now to challenge new puzzles, like e.g. the solar composition problems. Moreover, SSM neutrino flux predictions, which are directly affected by nuclear cross sections uncertainties, have to be F.L. Villante et al.   compared against very accurate observational determinations, having errors at few% level or better e.g. for Φ( 7 Be) and Φ( 8 B).
As a consequence, further work is needed on the side of nuclear reactions. Indeed, nuclear uncertainties have a non negligible role in SSMs error budget. As an example, the error contributions from S 34 and S 17 are about a factor 2 larger than the uncertainties in the experimental determinations of Φ( 7 Be) and Φ( 8 B). As it is discussed in Vinyoles et al. (2017), the few percent systematics in the determination of these reaction rates is still a relevant source of difficulty in using neutrino fluxes as constraints to solar model properties. The astrophysical factor S 114 is morevoer a relevant error source for CN neutrino fluxes. This last point is particularly important after Borexino opened the era of CNO neutrino detection, obtaining the first ever direct evidence of a non vanishing CN neutrino signal from the Sun. In the perspective of future and more accurate measurements, nuclear uncertainties can become a limiting factor in the possibility to use the CN-neutrinos, in combination with 8 B neutrinos, to directly probe the solar composition, thus addressing the solar composition problem. At the moment, the nuclear error contribution to CN-core abundance uncertainty is ∼ 10%, see eq. (37). This is comparable to the error in CN-surface abundance determinations (0.05 dex in LZ composition) and only a factor ∼ 2 smaller than the difference between HZ and LZ results, which can be regarded as an estimate of the systematic shift in the surface abundances produced by advances in stellar spectroscopy during the last 20 years. We remark that a high accuracy determination of the solar core composition could be used not only to discriminate among different solar surface admixtures but also to test the chemical evolution scheme employed by SSMs, e.g. by verifying the effect of elemental diffusion according to which core abundances are expected to be ∼ 15% larger than surface values.
In conclusion, it would be desirable to further improve our knowledge of nuclear cross sections, in particular for 3 He( 4 He, γ) 7 Be, 7 Be(p, γ) 8 B and 14 N(p, γ) 15 O reactions. As we discussed in the introduction, the history of SSMs appears to be formed by three large chapters, during which the knowledge of nuclear rates improved at each stage by about a factor two with respect to the previous period, up to the present situation in which the leading cross section in pp-chain and CN-cycle are typically determined with ∼ 5% accuracy. The ambitious goal for the next stage could be a further factor ∼ 2 reduction, in such a way that nuclear reactions uncertainties will not represent a limiting factor in constraining the physical conditions of solar interior.