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.

Nuclear astrophysics, the union of nuclear physics and astronomy, went through an impressive expansion during the last twenty years. This could be achieved thanks to milestone improvements in astronomical observations, cross section measurements, powerful computer simulations and much refined stellar models. Italian groups are giving quite important contributions to every domain of nuclear astrophysics, sometimes being the leaders of worldwide unique experiments. In this paper we will discuss the astrophysical scenarios where nuclear astrophysics plays a key role and we will provide detailed descriptions of the present and future of the experiments on nuclear astrophysics which belong to the scientific programme of INFN (the National Institute for Nuclear Physics in Italy).

Time modulations at per mil level have been reported to take place in the decay constant of several nuclei with period of one year (most cases) but also of about one month or one day. On the other hand, experiments with similar or better sensitivity have been unable to detect any modulation. In this letter we give the results of the activity study of two different sources: 40K and 226Ra. The two gamma spectrometry experiments have been performed underground at the Gran Sasso Laboratory, this way suppressing the time dependent cosmic ray background. Briefly, our measurements reached the sensitivity of 3.4 and 3.5 parts over 10^6 for 40K and 226Ra, respectively (1 sigma) and they do not show any statistically significant evidence of time dependence in the decay constant. We also give the results of the activity measurement at the time of the two strong X-class solar flares which took place in September 2017. Our data do not show any unexpected time dependence in the decay rate of 40K in correspondence with the two flares. To the best of our knowledge, these are the most precise and accurate results on the stability of the decay constant as function of time.

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.

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.

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.

We give the results of a study on the 222Rn decay we performed in the Gran Sasso Laboratory (LNGS) by detecting the gamma rays from the radon progeny. The motivation was to monitor the stability of radioactivity measuring several times per year the half-life of a short lifetime (days) source instead of measuring over a long period the activity of a long lifetime (tens or hundreds of years) source. In particular, we give a possible reason of the large periodical fluctuations in the count rate of the gamma rays due to radon inside a closed canister which has been described in literature and which has been attributed to a possible influence of a component in the solar irradiation affecting the nuclear decay rates. We then provide the result of four half-life measurements we performed underground at LNGS in the period from May 2014 to January 2015 with radon diffused into olive oil. Briefly, we did not measure any change of the 222Rn half-life with a 8*10^-5 precision. Finally, we provide the most precise value for the 222Rn half-life: 3.82146(16){stat}(4){syst} days.

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 observation of $^{26}$Al gives us the proof of active nucleosynthesis in the Milky Way. However the identification of the main producers of $^{26}$Al is still a matter of debate. Many sites have been proposed, but our poor knowledge of the nuclear processes involved introduces high uncertainties. In particular, the limited accuracy on the $^{25}$Mg($\alpha$,n)$^{28}$Si reaction cross section has been identified as the main source of nuclear uncertainty in the production of $^{26}$Al in C/Ne explosive burning in massive stars, which has been suggested to be the main source of $^{26}$Al in the Galaxy. We studied this reaction through neutron spectroscopy at the CN Van de Graaff accelerator of the Legnaro National Laboratories. Thanks to this technique we are able to discriminate the ($\alpha$,n) events from possible contamination arising from parasitic reactions. In particular, we measured the neutron angular distributions at 5 different beam energies (between 3 and 5 MeV) in the \ang{17.5}-\ang{106} laboratory system angular range. The presented results disagree with the assumptions introduced in the analysis of a previous experiment.

Time modulations at per mil level have been reported to take place in the decay constant of about 15 nuclei with period of one year (most cases) but also of about one month or one day. In this paper we give the results of the activity measurement of a 40K source and a 232Th one. The two experiments have been done at the Gran Sasso Laboratory during a period of about 500 days, above ground (40K) and underground (232Th) with a target sensitivity of a few parts over 10^5. We also give the results of the activity measurement at the time of the X-class solar flares which took place in May 2013. Briefly, our measurements do not show any evidence of unexpected time dependence in the decay rate of 40K and 232Th.

The deacay rate of three different radioactive sources 40K, 137Cs and natTh has been measured with NaI and Ge detectors. Data have been analyzed to search for possible variations in coincidence with the two strongest solar flares of the years 2011 and 2012. No significant deviations from standard expectation have been observed, with a few 10-4 sensitivity. As a consequence, we could not find any effect like that recently reported by Jenkins and Fischbach: a few per mil decrease in the decay rate of 54Mn during solar flares in December 2006.

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.

In this review we discuss the main theoretical aspects and experimental effects of neutrino electromagnetic properties. We start with a general description of the electromagnetic form factors of Dirac and Majorana neutrinos. Then, we discuss the theory and phenomenology of the magnetic and electric dipole moments, summarizing the experimental results and the theoretical predictions. We discuss also the phenomenology of a neutrino charge radius and radiative decay. Finally, we describe the theory of neutrino spin and spin-flavor precession in a transverse magnetic field and we summarize its phenomenological applications.

The primordial abundance of 7Li as predicted by Big Bang Nucleosynthesis (BBN) is more than a factor 2 larger than what has been observed in metal-poor halo stars. Herein, we analyze the possibility that this discrepancy originates from incorrect assumptions about the nuclear reaction cross sections relevant for BBN. To do this, we introduce an efficient method to calculate the changes in the 7Li abundance produced by arbitrary (temperature dependent) modifications of the nuclear reaction rates. Then, considering that 7Li is mainly produced from 7Be via the electron capture process 7Be + e -> 7Li + nu_e, we assess the impact of the various channels of 7Be destruction. Differently from previous analysis, we consider the role of unknown resonances by using a complete formalism which takes into account the effect of Coulomb and centrifugal barrier penetration and that does not rely on the use of the narrow-resonance approximation. As a result of this, the possibility of a nuclear physics solution to the 7Li problem is significantly suppressed. Given the present experimental and theoretical constraints, it is unlikely that the 7Be + n destruction rate is underestimated by the 2.5 factor required to solve the problem. We exclude, moreover, that resonant destruction in the channels 7Be + t and 7Be + 3He can explain the 7Li puzzle. New unknown resonances in 7Be + d and 7Be + alpha could potentially produce significant effects. Recent experimental results have ruled out such a possibility for 7Be+d. On the other hand, for the 7Be + alpha channel very favorable conditions are required. The possible existence of a partially suitable resonant level in 11C is studied in the framework of a coupled-channel model and the possibility of a direct measurement is considered.

Starting from June 2011, the activity of a 137Cs source has been measured by means of a HPGe detector installed deep underground in the Gran Sasso Laboratory. In total about 5100 energy spectra, one hour measuring time each, have been collected. These data allowed the search for time variations of the decay constant with periods from a few hours to 1 year. No signal with amplitude larger than 9.6{\cdot}10-5 at 95% C.L. has been detected. These limits are more than one order of magnitude lower than the values on the oscillation amplitude reported in literature. In particular, for 1 year period an oscillation amplitude larger than 8.5{\cdot}10-5 has been excluded at 95% C.L., independently of the phase. The same data give a value of 29.96{\pm}0.08 years for the 137Cs half life, in good agreement with the world mean value of 30.05{\pm}0.08 years.

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.

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.

LUNA, Laboratory for Underground Nuclear Astrophysics at Gran Sasso, is measuring fusion cross sections down to the energy of the nucleosynthesis inside stars. Outstanding results obtained up to now are the cross-section measurements within the Gamow peak of the Sun of $^{3}He(^{3}He,2p)^{4}He$ and the $D(p,\gamma)^{3}He$. The former plays a big role in the proton-proton chain, largely affecting the calculated solar neutrino luminosity, whereas the latter is the reaction that rules the proto-star life during the pre-main sequence phase. The implications of such measurements will be discussed. Preliminary results obtained last year on the study of $^{14}N(p,\gamma)^{15}O$, the slowest reaction of the CNO cycle, will also be shown.

We report on the results of the \hethet\ experiment at the underground accelerator facility LUNA (Gran Sasso). For the first time the lowest projectile energies utilized for the cross section measurement correspond to energies below the center of the solar Gamow peak ($E_{\rm 0}$=22 keV). The data provide no evidence for the existence of a hypothetical resonance in the energy range investigated. Although no extrapolation is needed anymore (except for energies at the low-energy tail of the Gamow peak), the data must be corrected for the effects of electron screening, clearly observed the first time for the \hethet\ reaction. The effects are however larger than expected and not understood, leading presently to the largest uncertainty on the quoted $S_{\rm b}(E_{\rm 0})$ value for bare nuclides (=5.40 MeV b).