Among the reactions involved in the production and destruction of deuterium during Big Bang Nucleosynthesis, the deuterium-burning D(p,gamma)3He reaction has the largest uncertainty and limits the precision of theoretical estimates of primordial deuterium abundance. Here we report the results of a careful commissioning of the experimental setup used to measure the cross-section of the D(p,gamma)3He reaction at the Laboratory for Underground Nuclear Astrophysics of the Gran Sasso Laboratory (Italy). The commissioning was aimed at minimising all sources of systematic uncertainty in the measured cross sections. The overall systematic error achieved (< 3 %) will enable improved predictions of BBN deuterium abundance.

The experimental study of nuclear reactions of astrophysical interest is greatly facilitated by a low-background, high-luminosity setup. The Laboratory for Underground Nuclear Astrophysics (LUNA) 400 kV accelerator offers ultra-low cosmic-ray induced background due to its location deep underground in the Gran Sasso National Laboratory (INFN-LNGS), Italy, and high intensity, 250-500 $\mu$A, proton and $\alpha$ ion beams. In order to fully exploit these features, a high-purity, recirculating gas target system for isotopically enriched gases is coupled to a high-efficiency, six-fold optically segmented bismuth germanate (BGO) $\gamma$-ray detector. The beam intensity is measured with a beam calorimeter with constant temperature gradient. Pressure and temperature measurements have been carried out at several positions along the beam path, and the resultant gas density profile has been determined. Calibrated $\gamma$-intensity standards and the well-known $E_p$ = 278 keV $\mathrm{^{14}N(p,\gamma)^{15}O}$ resonance were used to determine the $\gamma$-ray detection efficiency and to validate the simulation of the target and detector setup. As an example, the recently measured resonance at $E_p$ = 189.5 keV in the $^{22}$Ne(p,$\gamma$)$^{23}$Na reaction has been investigated with high statistics, and the $\gamma$-decay branching ratios of the resonance have been determined.

The essential ingredients of nuclear astrophysics are the thermonuclear reactions which shape the life and death of stars and which are responsible for the synthesis of the chemical elements in the Universe. Deep underground in the Gran Sasso Laboratory the cross sections of the key reactions responsible for the hydrogen burning in stars have been measured with two accelerators of 50 and 400 kV voltage right down to the energies of astrophysical interest. As a matter of fact, the main advantage of the underground laboratory is the reduction of the background. Such a reduction has allowed, for the first time, to measure relevant cross sections at the Gamow energy. The qualifying features of underground nuclear astrophysics are exhaustively reviewed before discussing the current LUNA program which is mainly devoted to the study of the Big-Bang nucleosynthesis and of the synthesis of the light elements in AGB stars and classical novae. The main results obtained during the study of reactions relevant to the Sun are also reviewed and their influence on our understanding of the properties of the neutrino, of the Sun and of the Universe itself is discussed. Finally, the future of LUNA during the next decade is outlined. It will be mainly focused on the study of the nuclear burning stages after hydrogen burning: helium and carbon burning. All this will be accomplished thanks to a new 3.5 MV accelerator able to deliver high current beams of proton, helium and carbon which will start running under Gran Sasso in 2019. In particular, we will discuss the first phase of the scientific case of the 3.5 MV accelerator focused on the study of $^{12}$C+$^{12}$C and of the two reactions which generate free neutrons inside stars: $^{13}$C($\alpha$,n)$^{16}$O and $^{22}$Ne($\alpha$,n)$^{25}$Mg.

The $^{14}\textrm{N(p,}\gamma\textrm{)}^{15}\textrm{O}$ reaction is the slowest reaction of the carbon-nitrogen cycle of hydrogen burning and thus determines its rate. The precise knowledge of its rate is required to correctly model hydrogen burning in asymptotic giant branch stars. In addition, it is a necessary ingredient for a possible solution of the solar abundance problem by using the solar $^{13}$N and $^{15}$O neutrino fluxes as probes of the carbon and nitrogen abundances in the solar core. After the downward revision of its cross section due to a much lower contribution by one particular transition, capture to the ground state in $^{15}$O, the evaluated total uncertainty is still 8\%, in part due to an unsatisfactory knowledge of the excitation function over a wide energy range. The present work reports precise S-factor data at twelve energies between 0.357-1.292~MeV for the strongest transition, capture to the 6.79~MeV excited state in $^{15}$O, and at ten energies between 0.479-1.202~MeV for the second strongest transition, capture to the ground state in $^{15}$O. An R-matrix fit is performed to estimate the impact of the new data on astrophysical energies. The recently suggested slight enhancement of the 6.79~MeV transition at low energy could not be confirmed. The present extrapolated zero-energy S-factors are $S_{6.79}(0)$~=~1.24$\pm$0.11~keV~barn and $S_{\rm GS}(0)$~=~0.19$\pm$0.05~keV~barn.

Among the solar fusion reactions, the rate of the $^7$Be($p,\gamma$)$^8$B reaction is one of the most difficult to determine rates. In a number of previous experiments, its astrophysical $S$-factor has been measured at $E$ = 0.1-2.5 MeV center-of-mass energy. However, no experimental data is available below 0.1 MeV. Thus, an extrapolation to solar energies is necessary, resulting in significant uncertainty for the extrapolated $S$-factor. On the other hand, the measured solar neutrino fluxes are now very precise. Therefore, the problem of the $S$-factor determination is turned around here: Using the measured $^7$Be and $^8$B neutrino fluxes and the Standard Solar Model, the $^7$Be($p,\gamma$)$^8$B astrophysical $S$-factor is determined at the solar Gamow peak. In addition, the $^3$He($\alpha$,$\gamma$)$^7$Be $S$-factor is redetermined with a similar method.

Stardust grains recovered from meteorites provide high-precision snapshots of the isotopic composition of the stellar environment in which they formed. Attributing their origin to specific types of stars, however, often proves difficult. Intermediate-mass stars of 4-8 solar masses are expected to contribute a large fraction of meteoritic stardust. However, no grains have been found with characteristic isotopic compositions expected from such stars. This is a long-standing puzzle, which points to serious gaps in our understanding of the lifecycle of stars and dust in our Galaxy. Here we show that the increased proton-capture rate of $^{17}$O reported by a recent underground experiment leads to $^{17}$O/$^{16}$O isotopic ratios that match those observed in a population of stardust grains, for proton-burning temperatures of 60-80 million K. These temperatures are indeed achieved at the base of the convective envelope during the late evolution of intermediate-mass stars of 4-8 solar masses, which reveals them as the most likely site of origin of the grains. This result provides the first direct evidence that these stars contributed to the dust inventory from which the Solar System formed.

We investigate the impact of the new LUNA rate for the nuclear reaction $^{22}$Ne$(p,\gamma)^{23}$Na on the chemical ejecta of intermediate-mass stars, with particular focus on the thermally-pulsing asymptotic giant branch (TP-AGB) stars that experience hot-bottom burning. To this aim we use the PARSEC and COLIBRI codes to compute the complete evolution, from the pre-main sequence up to the termination of the TP-AGB phase, of a set of stellar models with initial masses in the range $3.0\,M_{\odot} - 6.0\,M_{\odot}$, and metallicities $Z_{\rm i}=0.0005$, $Z_{\rm i}=0.006$, and $Z_{\rm i} = 0.014$. We find that the new LUNA measures have much reduced the nuclear uncertainties of the $^{22}$Ne and $^{23}$Na AGB ejecta, which drop from factors of $\simeq 10$ to only a factor of few for the lowest metallicity models. Relying on the most recent estimations for the destruction rate of $^{23}$Na, the uncertainties that still affect the $^{22}$Ne and $^{23}$Na AGB ejecta are mainly dominated by evolutionary aspects (efficiency of mass-loss, third dredge-up, convection). Finally, we discuss how the LUNA results impact on the hypothesis that invokes massive AGB stars as the main agents of the observed O-Na anti-correlation in Galactic globular clusters. We derive quantitative indications on the efficiencies of key physical processes (mass loss, third dredge-up, sodium destruction) in order to simultaneously reproduce both the Na-rich, O-poor extreme of the anti-correlation, and the observational constraints on the CNO abundance. Results for the corresponding chemical ejecta are made publicly available.

Context. Material processed by the CNO cycle in stellar interiors is enriched in 17O. When mixing processes from the stellar surface reach these layers, as occurs when stars become red giants and undergo the first dredge up, the abundance of 17O increases. Such an occurrence explains the drop of the 16O/17O observed in RGB stars with mass larger than 1.5 M_\solar. As a consequence, the interstellar medium is continuously polluted by the wind of evolved stars enriched in 17O . Aims. Recently, the Laboratory for Underground Nuclear Astrophysics (LUNA) collaboration released an improved rate of the 17O(p,alpha)14N reaction. In this paper we discuss the impact that the revised rate has on the 16O/17O ratio at the stellar surface and on 17O stellar yields. Methods. We computed stellar models of initial mass between 1 and 20 M_\solar and compared the results obtained by adopting the revised rate of the 17O(p,alpha)14N to those obtained using previous rates. Results. The post-first dredge up 16O/17O ratios are about 20% larger than previously obtained. Negligible variations are found in the case of the second and the third dredge up. In spite of the larger 17O(p,alpha)14N rate, we confirm previous claims that an extra-mixing process on the red giant branch, commonly invoked to explain the low carbon isotopic ratio observed in bright low-mass giant stars, marginally affects the 16O/17O ratio. Possible effects on AGB extra-mixing episodes are also discussed. As a whole, a substantial reduction of 17O stellar yields is found. In particular, the net yield of stars with mass ranging between 2 and 20 M_\solar is 15 to 40% smaller than previously estimated. Conclusions. The revision of the 17O(p,alpha)14N rate has a major impact on the interpretation of the 16O/17O observed in evolved giants, in stardust grains and on the 17O stellar yields.

The $^{22}$Ne(p,$\gamma$)$^{23}$Na reaction is the most uncertain process in the neon-sodium cycle of hydrogen burning. At temperatures relevant for nucleosynthesis in asymptotic giant branch stars and classical novae, its uncertainty is mainly due to a large number of predicted but hitherto unobserved resonances at low energy. Purpose: A new direct study of low energy $^{22}$Ne(p,$\gamma$)$^{23}$Na resonances has been performed at the Laboratory for Underground Nuclear Astrophysics (LUNA), in the Gran Sasso National Laboratory, Italy. Method: The proton capture on $^{22}$Ne was investigated in direct kinematics, delivering an intense proton beam to a $^{22}$Ne gas target. $\gamma$ rays were detected with two high-purity germanium detectors enclosed in a copper and lead shielding suppressing environmental radioactivity. Results: Three resonances at 156.2 keV ($\omega\gamma$ = (1.48\,$\pm$\,0.10)\,$\cdot$\,10$^{-7}$ eV), 189.5 keV ($\omega\gamma$ = (1.87\,$\pm$\,0.06)\,$\cdot$\,10$^{-6}$ eV) and 259.7 keV ($\omega\gamma$ = (6.89\,$\pm$\,0.16)\,$\cdot$\,10$^{-6}$ eV) proton beam energy, respectively, have been observed for the first time. For the levels at 8943.5, 8975.3, and 9042.4 keV excitation energy corresponding to the new resonances, the $\gamma$-decay branching ratios have been precisely measured. Three additional, tentative resonances at 71, 105 and 215 keV proton beam energy, respectively, were not observed here. For the strengths of these resonances, experimental upper limits have been derived that are significantly more stringent than the upper limits reported in the literature. Conclusions: Based on the present experimental data and also previous literature data, an updated thermonuclear reaction rate is provided in tabular and parametric form. The new reaction rate is significantly higher than previous evaluations at temperatures of 0.08-0.3 GK.

With the R$^{3}$B-LAND setup at GSI we have measured exclusive relative-energy spectra of the Coulomb dissociation of $^{18}$C at a projectile energy around 425~AMeV on a lead target, which are needed to determine the radiative neutron-capture cross sections of $^{17}$C into the ground state of $^{18}$C. Those data have been used to constrain theoretical calculations for transitions populating excited states in $^{18}$C. This allowed to derive the astrophysical cross section $\sigma^{*}_{\mathrm{n}\gamma}$ accounting for the thermal population of $^{17}$C target states in astrophysical scenarios. The experimentally verified capture rate is significantly lower than those of previously obtained Hauser-Feshbach estimations at temperatures $T_{9}\leq{}1$~GK. Network simulations with updated neutron-capture rates and hydrodynamics according to the neutrino-driven wind model as well as the neutron-star merger scenario reveal no pronounced influence of neutron capture of $^{17}$C on the production of second- and third-peak elements in contrast to earlier sensitivity studies.

The $^{22}$Ne(p,$\gamma$)$^{23}$Na reaction takes part in the neon-sodium cycle of hydrogen burning. This cycle affects the synthesis of the elements between $^{20}$Ne and $^{27}$Al in asymptotic giant branch stars and novae. The $^{22}$Ne(p,$\gamma$)$^{23}$Na reaction rate is very uncertain because of a large number of unobserved resonances lying in the Gamow window. At proton energies below 400\,keV, only upper limits exist in the literature for the resonance strengths. Previous reaction rate evaluations differ by large factors. In the present work, the first direct observations of the $^{22}$Ne(p,$\gamma$)$^{23}$Na resonances at 156.2, 189.5, and 259.7\,keV are reported. Their resonance strengths have been derived with 2-7\% uncertainty. In addition, upper limits for three other resonances have been greatly reduced. Data were taken using a windowless $^{22}$Ne gas target and high-purity germanium detectors at the Laboratory for Underground Nuclear Astrophysics in the Gran Sasso laboratory of the National Institute for Nuclear Physics, Italy, taking advantage of the ultra-low background observed deep underground. The new reaction rate is a factor of 5 higher than the recent evaluation at temperatures relevant to novae and asymptotic giant branch stars nucleosynthesis.

The 3He({\alpha},{\gamma})7Be reaction affects not only the production of 7Li in Big Bang nucleosynthesis, but also the fluxes of 7Be and 8B neutrinos from the Sun. This double role is exploited here to constrain the former by the latter. A number of recent experiments on 3He({\alpha},{\gamma})7Be provide precise cross section data at E = 0.5-1.0 MeV center-of-mass energy. However, there is a scarcity of precise data at Big Bang energies, 0.1-0.5 MeV, and below. This problem can be alleviated, based on precisely calibrated 7Be and 8B neutrino fluxes from the Sun that are now available, assuming the neutrino flavour oscillation framework to be correct. These fluxes and the standard solar model are used here to determine the 3He(alpha,gamma)7Be astrophysical S-factor at the solar Gamow peak, S(23+6-5 keV) = 0.548+/-0.054 keVb. This new data point is then included in a re-evaluation of the 3He({\alpha},{\gamma})7Be S-factor at Big Bang energies, following an approach recently developed for this reaction in the context of solar fusion studies. The re-evaluated S-factor curve is then used to re-determine the 3He({\alpha},{\gamma})7Be thermonuclear reaction rate at Big Bang energies. The predicted primordial lithium abundance is 7Li/H = 5.0e-10, far higher than the Spite plateau.

The main background above 3\,MeV for in-beam nuclear astrophysics studies with $\gamma$-ray detectors is caused by cosmic-ray induced secondaries. The two commonly used suppression methods, active and passive shielding, against this kind of background were formerly considered only as alternatives in nuclear astrophysics experiments. In this work the study of the effects of active shielding against cosmic-ray induced events at a medium deep location is performed. Background spectra were recorded with two actively shielded HPGe detectors. The experiment was located at 148\,m below the surface of the Earth in the Reiche Zeche mine in Freiberg, Germany. The results are compared to data with the same detectors at the Earth's surface, and at depths of 45\,m and 1400\,m, respectively.

The $^{22}$Ne(p,$\gamma$)$^{23}$Na reaction takes part in the neon-sodium cycle of hydrogen burning. This cycle is active in asymptotic giant branch stars as well as in novae and contributes to the nucleosythesis of neon and sodium isotopes. In order to reduce the uncertainties in the predicted nucleosynthesis yields, new experimental efforts to measure the $^{22}$Ne(p,$\gamma$)$^{23}$Na cross section directly at the astrophysically relevant energies are needed. In the present work, a feasibility study for a $^{22}$Ne(p,$\gamma$)$^{23}$Na experiment at the Laboratory for Underground Nuclear Astrophysics (LUNA) 400\,kV accelerator deep underground in the Gran Sasso laboratory, Italy, is reported. The ion beam induced $\gamma$-ray background has been studied. The feasibility study led to the first observation of the $E_{\rm p}$ = 186\,keV resonance in a direct experiment. An experimental lower limit of 0.12\,$\times$\,10$^{-6}$\,eV has been obtained for the resonance strength. Informed by the feasibility study, a dedicated experimental setup for the $^{22}$Ne(p,$\gamma$)$^{23}$Na experiment has been developed. The new setup has been characterized by a study of the temperature and pressure profiles. The beam heating effect that reduces the effective neon gas density due to the heating by the incident proton beam has been studied using the resonance scan technique, and the size of this effect has been determined for a neon gas target.

The nucleosynthesis of elements beyond iron is dominated by neutron captures in the s and r processes. However, 32 stable, proton-rich isotopes cannot be formed during those processes, because they are shielded from the s-process flow and r-process beta-decay chains. These nuclei are attributed to the p and rp process. For all those processes, current research in nuclear astrophysics addresses the need for more precise reaction data involving radioactive isotopes. Depending on the particular reaction, direct or inverse kinematics, forward or time-reversed direction are investigated to determine or at least to constrain the desired reaction cross sections. The Facility for Antiproton and Ion Research (FAIR) will offer unique, unprecedented opportunities to investigate many of the important reactions. The high yield of radioactive isotopes, even far away from the valley of stability, allows the investigation of isotopes involved in processes as exotic as the r or rp processes.

The 40Ca(alpha,gamma)44Ti reaction is believed to be the main production channel for the radioactive nuclide 44Ti in core-collapse supernovae. Radiation from decaying 44Ti has been observed so far for two supernova remnants, and a precise knowledge of the 44Ti production rate may help improve supernova models. The 40Ca(alpha,gamma)44Ti astrophysical reaction rate is determined by a number of narrow resonances. Here, the resonance triplet at E_alpha = 4497, 4510, and 4523 keV is studied both by activation, using an underground laboratory for the gamma counting, and by in-beam gamma spectrometry. The target properties are determined by elastic recoil detection analysis and by nuclear reactions. The strengths of the three resonances are determined to omega gamma = (0.92+-0.20), (6.2+-0.5), and (1.32+-0.24) eV, respectively, a factor of two more precise than before. The strengths of this resonance triplet may be used in future works as a point of reference. In addition, the present new data directly affect the astrophysical reaction rate at relatively high temperatures, above 3.5 GK.

The production of the stable isotope Li-6 in standard Big Bang nucleosynthesis has recently attracted much interest. Recent observations in metal-poor stars suggest that a cosmological Li-6 plateau may exist. If true, this plateau would come in addition to the well-known Spite plateau of Li-7 abundances and would point to a predominantly primordial origin of Li-6, contrary to the results of standard Big Bang nucleosynthesis calculations. Therefore, the nuclear physics underlying Big Bang Li-6 production must be revisited. The main production channel for Li-6 in the Big Bang is the 2H(alpha,gamma)6Li reaction. The present work reports on neutron-induced effects in a high-purity germanium detector that were encountered in a new study of this reaction. In the experiment, an {\alpha}-beam from the underground accelerator LUNA in Gran Sasso, Italy, and a windowless deuterium gas target are used. A low neutron flux is induced by energetic deuterons from elastic scattering and, subsequently, the 2H(d,n)3He reaction. Due to the ultra-low laboratory neutron background at LUNA, the effect of this weak flux of 2-3 MeV neutrons on well-shielded high-purity germanium detectors has been studied in detail. Data have been taken at 280 and 400 keV alpha-beam energy and for comparison also using an americium-beryllium neutron source.

Proton captures on Mg isotopes play an important role in the Mg-Al cycle active in stellar H-burning regions. In particular, low-energy nuclear resonances in the $^{25}$Mg(p,$\gamma$)$^{26}$Al reaction affect the production of radioactive $^{26}$Al$^{gs}$ as well as the resulting Mg/Al abundance ratio. Reliable estimations of these quantities require precise measurements of the strengths of low-energy resonances. Based on a new experimental study performed at LUNA, we provide revised rates of the $^{25}$Mg(p,$\gamma$)$^{26}$Al$^{gs}$ and the $^{25}$Mg(p,$\gamma$)$^{26}$Al$^{m}$ reactions with corresponding uncertainties. In the temperature range 50 to 150 MK, the new recommended rate of the $^{26}$Al$^{m}$ production is up to 5 times higher than previously assumed. In addition, at T$=100$ MK, the revised total reaction rate is a factor of 2 higher. Note that this is the range of temperature at which the Mg-Al cycle operates in an H-burning zone. The effects of this revision are discussed. Due to the significantly larger $^{25}$Mg(p,$\gamma$)$^{26}$Al$^{m}$ rate, the estimated production of $^{26}$Al$^{gs}$ in H-burning regions is less efficient than previously obtained. As a result, the new rates should imply a smaller contribution from Wolf-Rayet stars to the galactic $^{26}$Al budget. Similarly, we show that the AGB extra-mixing scenario does not appear able to explain the most extreme values of $^{26}$Al/$^{27}$Al, i.e. $>10^{-2}$, found in some O-rich presolar grains. Finally, the substantial increase of the total reaction rate makes the hypothesis of a self-pollution by massive AGBs a more robust explanation for the Mg-Al anticorrelation observed in Globular-Cluster stars.

Classical novae are important contributors to the abundances of key isotopes, such as the radioactive ^{18}F, whose observation by satellite missions could provide constraints on nucleosynthesis models in novae. The ^{17}O(p,\gamma)^{18}F reaction plays a critical role in the synthesis of both oxygen and fluorine isotopes but its reaction rate is not well determined because of the lack of experimental data at energies relevant to novae explosions. In this study, the reaction cross section has been measured directly for the first time in a wide energy range Ecm = 200 - 370 keV appropriate to hydrogen burning in classical novae. In addition, the E=183 keV resonance strength, \omega \gamma=1.67\pm0.12 \mueV, has been measured with the highest precision to date. The uncertainty on the ^{17}O(p,\gamma)^{18}F reaction rate has been reduced by a factor of 4, thus leading to firmer constraints on accurate models of novae nucleosynthesis.

The NO cycle takes place in the deepest layer of a H-burning core or shell, when the temperature exceeds T {\simeq} 30 {\cdot} 106 K. The O depletion observed in some globular cluster giant stars, always associated with a Na enhancement, may be due to either a deep mixing during the RGB (red giant branch) phase of the star or to the pollution of the primordial gas by an early population of massive AGB (asymptotic giant branch) stars, whose chemical composition was modified by the hot bottom burning. In both cases, the NO cycle is responsible for the O depletion. The activation of this cycle depends on the rate of the 15N(p,{\gamma})16O reaction. A precise evaluation of this reaction rate at temperatures as low as experienced in H-burning zones in stellar interiors is mandatory to understand the observed O abundances. We present a new measurement of the 15N(p,{\gamma})16O reaction performed at LUNA covering for the first time the center of mass energy range 70-370 keV, which corresponds to stellar temperatures between 65 {\cdot} 106 K and 780 {\cdot}106 K. This range includes the 15N(p,{\gamma})16O Gamow-peak energy of explosive H-burning taking place in the external layer of a nova and the one of the hot bottom burning (HBB) nucleosynthesis occurring in massive AGB stars. With the present data, we are also able to confirm the result of the previous R-matrix extrapolation. In particular, in the temperature range of astrophysical interest, the new rate is about a factor of 2 smaller than reported in the widely adopted compilation of reaction rates (NACRE or CF88) and the uncertainty is now reduced down to the 10% level.

The rate of the carbon-nitrogen-oxygen (CNO) cycle of hydrogen burning is controlled by the 14N(p,gamma)15O reaction. The reaction proceeds by capture to the ground states and several excited states in O-15. In order to obtain a reliable extrapolation of the excitation curve to astrophysical energy, fits in the R-matrix framework are needed. In an energy range that sensitively tests such fits, new cross section data are reported here for the four major transitions in the 14N(p,gamma)15O reaction. The experiment has been performed at the Laboratory for Underground Nuclear Astrophysics (LUNA) 400 kV accelerator placed deep underground in the Gran Sasso facility in Italy. Using a composite germanium detector, summing corrections have been considerably reduced with respect to previous studies. The cross sections for capture to the ground state and to the 5181, 6172, and 6792 keV excited states in O-15 have been determined at 359, 380, and 399 keV beam energy. In addition, the branching ratios for the decay of the 278 keV resonance have been remeasured.

Nuclear astrophysics strives for a comprehensive picture of the nuclear reactions responsible for synthesizing the chemical elements and for powering the stellar evolution engine. Deep underground in the Gran Sasso laboratory the cross sections of the key reactions of the proton-proton chain and of the Carbon-Nitrogen-Oxygen (CNO) cycle have been measured right down to the energies of astrophysical interest. The salient features of underground nuclear astrophysics are summarized here. The main results obtained by LUNA in the last twenty years are reviewed, and their influence on the comprehension of the properties of the neutrino, of the Sun and of the Universe itself are discussed. Future directions of underground nuclear astrophysics towards the study of helium and carbon burning and of stellar neutron sources in stars are pointed out.

We summarize and critically evaluate the available data on nuclear fusion cross sections important to energy generation in the Sun and other hydrogen-burning stars and to solar neutrino production. Recommended values and uncertainties are provided for key cross sections, and a recommended spectrum is given for 8B solar neutrinos. We also discuss opportunities for further increasing the precision of key rates, including new facilities, new experimental techniques, and improvements in theory. This review, which summarizes the conclusions of a workshop held at the Institute for Nuclear Theory, Seattle, in January 2009, is intended as a 10-year update and supplement to Reviews of Modern Physics 70 (1998) 1265.

The 15N(p,gamma)16O reaction controls the passage of nucleosynthetic material from the first to the second carbon-nitrogen-oxygen (CNO) cycle. A direct measurement of the total 15N(p,gamma)16O cross section at energies corresponding to hydrogen burning in novae is presented here. Data have been taken at 90-230 keV center-of-mass energy using a windowless gas target filled with nitrogen of natural isotopic composition and a bismuth germanate summing detector. The cross section is found to be a factor two lower than previously believed.

The rate of the hydrogen-burning carbon-nitrogen-oxygen (CNO) cycle is controlled by the slowest process, 14N(p,gamma)15O, which proceeds by capture to the ground and several excited states in 15O. Previous extrapolations for the ground state contribution disagreed by a factor 2, corresponding to 15% uncertainty in the total astrophysical S-factor. At the Laboratory for Underground Nuclear Astrophysics (LUNA) 400 kV accelerator placed deep underground in the Gran Sasso facility in Italy, a new experiment on ground state capture has been carried out at 317.8, 334.4, and 353.3 keV center-of-mass energy. Systematic corrections have been reduced considerably with respect to previous studies by using a Clover detector and by adopting a relative analysis. The previous discrepancy has been resolved, and ground state capture no longer dominates the uncertainty of the total S-factor.

Solar neutrino fluxes depend both on astrophysical and on nuclear physics inputs, namely on the cross sections of the reactions responsible for neutrino production inside the Solar core. While the flux of solar 8B neutrinos has been recently measured at Superkamiokande with a 3.5% uncertainty and a precise measurement of 7Be neutrino flux is foreseen in the next future, the predicted fluxes are still affected by larger errors. The largest nuclear physics uncertainty to determine the fluxes of 8B and 7Be neutrinos comes from the 3He(alpha,gamma)7Be reaction. The uncertainty on its S-factor is due to an average discrepancy in results obtained using two different experimental approaches: the detection of the delayed gamma rays from 7Be decay and the measurement of the prompt gamma emission. Here we report on a new high precision experiment performed with both techniques at the same time. Thanks to the low background conditions of the Gran Sasso LUNA accelerator facility, the cross section has been measured at Ecm = 170, 106 and 93 keV, the latter being the lowest interaction energy ever reached. The S-factors from the two methods do not show any discrepancy within the experimental errors. An extrapolated S(0)= 0.560+/-0.017 keV barn is obtained. Moreover, branching ratios between the two prompt gamma-transitions have been measured with 5-8% accuracy.

The flux of 7Be and 8B neutrinos from the Sun and the production of 7Li via primordial nucleosynthesis depend on the rate of the 3He(alpha,gamma)7Be reaction. In extension of a previous study showing cross section data at 127 - 167 keV center of mass energy, the present work reports on a measurement of the 3He(alpha,gamma)7Be cross section at 106 keV performed at Italy's Gran Sasso underground laboratory by the activation method. This energy is closer to the solar Gamow energy than ever reached before. The result is sigma = 0.567 +- 0.029(stat) +- 0.016(syst) nbarn. The data are compared with previous activation studies at high energy, and a recommended S(0) value for all 3He(alpha,gamma)7Be activation studies, including the present work, is given.

The nuclear physics input from the 3He(alpha,gamma)7Be cross section is a major uncertainty in the fluxes of 7Be and 8B neutrinos from the Sun predicted by solar models and in the 7Li abundance obtained in big-bang nucleosynthesis calculations. The present work reports on a new precision experiment using the activation technique at energies directly relevant to big-bang nucleosynthesis. Previously such low energies had been reached experimentally only by the prompt-gamma technique and with inferior precision. Using a windowless gas target, high beam intensity and low background gamma-counting facilities, the 3He(alpha,gamma)7Be cross section has been determined at 127, 148 and 169 keV center-of-mass energy with a total uncertainty of 4%. The sources of systematic uncertainty are discussed in detail. The present data can be used in big-bang nucleosynthesis calculations and to constrain the extrapolation of the 3He(alpha,gamma)7Be astrophysical S-factor to solar energies.

In stars with temperatures above 20*10^6 K, hydrogen burning is dominated by the CNO cycle. Its rate is determined by the slowest process, the 14N(p,gamma)15O reaction. Deep underground in Italy's Gran Sasso laboratory, at the LUNA 400 kV accelerator, the cross section of this reaction has been measured at energies much lower than ever achieved before. Using a windowless gas target and a 4pi BGO summing detector, direct cross section data has been obtained down to 70 keV, reaching a value of 0.24 picobarn. The Gamow peak has been covered by experimental data for several scenarios of stable and explosive hydrogen burning. In addition, the strength of the 259 keV resonance has been remeasured. The thermonuclear reaction rate has been calculated for temperatures 90 - 300 *10^6 K, for the first time with negligible impact from extrapolations.

The transition between the Main Sequence and the Red Giant Branch in low mass stars is powered by the onset of the CNO burning, whose bottleneck is the $^{14}$N(p,$\gamma)^{15}$O. The LUNA collaboration has recently improved the low energy measurements of the cross section of this key reaction. We analyse the impact of the revised reaction rate on the estimate of the Globular Clusters ages, as derived from the turnoff luminosity. We found that the age of the oldest Globulars should be increased by about 0.7-1 Gyr with respect to the current estimates.