A bubble chamber has been developed to be used as an active target system for low energy nuclear astrophysics experiments. Adopting ideas from dark matter detection with superheated liquids, a detector system compatible with gamma-ray beams has been developed. This detector alleviates some of the limitations encountered in standard measurements of the minute cross sections of interest to stellar environments. While the astrophysically relevant nuclear reaction processes at hydrostatic burning temperatures are dominated by radiative captures, in this experimental scheme we measure the time-reversed processes. Such photodisintegrations allow us to compute the radiative capture cross sections when transitions to excited states of the reaction products are negligible. Due to the transformation of phase space, the photodisintegration cross sections are up to two orders of magnitude higher. The main advantage of the new target-detector system is a density several orders of magnitude higher than conventional gas targets. Also, the detector is virtually insensitive to the gamma-ray beam itself, thus allowing us to detect only the products of the nuclear reaction of interest. The development and the operation as well as the advantages and disadvantages of the bubble chamber are discussed.

We have devised a technique for measuring some of the most important nuclear reactions in stars which we expect to provide considerable improvement over previous experiments. Adapting ideas from dark matter search experiments with bubble chambers, we have found that a superheated liquid is sensitive to recoils produced from \gamma-rays photodisintegrating the nuclei of the liquid. The main advantage of the new target-detector system is a gain in yield of six orders of magnitude over conventional gas targets due to the higher mass density of liquids. Also, the detector is practically insensitive to the \gamma-ray beam itself, thus allowing it to detect only the products of the nuclear reaction of interest. The first set of tests of a superheated target with a narrow bandwidth \gamma-ray beam was completed and the results demonstrate the feasibility of the scheme. The new data are successfully described by an R-matrix model using published resonance parameters. With the increase in luminosity of the next generation \gamma-ray beam facilities, the measurement of thermonuclear rates in the stellar Gamow window would become possible.

A method based on Monte Carlo techniques is presented for evaluating thermonuclear reaction rates. We begin by reviewing commonly applied procedures and point out that reaction rates that have been reported up to now in the literature have no rigorous statistical meaning. Subsequently, we associate each nuclear physics quantity entering in the calculation of reaction rates with a specific probability density function, including Gaussian, lognormal and chi-squared distributions. Based on these probability density functions the total reaction rate is randomly sampled many times until the required statistical precision is achieved. This procedure results in a median (Monte Carlo) rate which agrees under certain conditions with the commonly reported recommended "classical" rate. In addition, we present at each temperature a low rate and a high rate, corresponding to the 0.16 and 0.84 quantiles of the cumulative reaction rate distribution. These quantities are in general different from the statistically meaningless "minimum" (or "lower limit") and "maximum" (or "upper limit") reaction rates which are commonly reported. Furthermore, we approximate the output reaction rate probability density function by a lognormal distribution and present, at each temperature, the lognormal parameters miu and sigma. The values of these quantities will be crucial for future Monte Carlo nucleosynthesis studies. Our new reaction rates, appropriate for bare nuclei in the laboratory, are tabulated in the second paper of this series (Paper II). The nuclear physics input used to derive our reaction rates is presented in the third paper of this series (Paper III). In the fourth paper of this series (Paper IV) we compare our new reaction rates to previous results.

We present new rates for the 22Ne(alpha, n)25Mg and 22Ne(alpha,gamma)26Mg reactions, with uncertainties that have been considerably reduced compared to previous estimates, and we study how these new rates affect the production of the heavy magnesium isotopes in models of intermediate mass Asymptotic Giant Branch (AGB) stars of different initial compositions. All the models have deep third dredge-up, hot bottom burning and mass loss. Calculations have been performed using the two most commonly used estimates of the 22Ne + alpha rates as well as the new recommended rates, and with combinations of their upper and lower limits. The main result of the present study is that with the new rates, uncertainties on the production of isotopes from Mg to P coming from the 22Ne + alpha-capture rates have been considerably reduced. We have therefore removed one of the important sources of uncertainty to effect models of AGB stars. We have studied the effects of varying the mass-loss rate on nucleosynthesis and discuss other uncertainties related to the physics employed in the computation of stellar structure, such as the modeling of convection, the inclusion of a partial mixing zone and the definition of convective borders. These uncertainties are found to be much larger than those coming from 22Ne + alpha-capture rates, when using our new estimates. Much effort is needed to improve the situation for AGB models.

In the light of recent recalculations of the \el{19}{F}(alpha,p)\el{22}{Ne} reaction rate we present results of the expected yield of \el{19}{F} from Wolf-Rayet (WR) stars. We have computed models using the upper and lower limits for the rate in addition to the recommended rate and hence we constrain the uncertainty in the yield with respect to this reaction. We find a yield of 3.1x10^{-4}M_sol of \el{19}{F} with our recommended rate and a difference of a factor of two between the yields computed with the upper and lower limits. In comparison with previous work we find a difference in the yield of approximately a factor of 4, connected to a different choice of mass loss. Model uncertainties must be carefully evaluated in order to obtain a reliable estimate of the yield of fluorine from WR stars together with its uncertainties.

We present nucleosynthesis calculations and the resulting 19F stellar yields for a large set of models with different masses and metallicity. We find that the production of fluorine depends on the temperature of the convective pulses, the amount of primary 12C mixed into the envelope by third dredge up and the extent of the partial mixing zone. Then we perform a detailed analysis of the reaction rates involved in the production of 19F and the effects of their uncertainties. We find that the major uncertainties are associated with the 14C(alpha,gamma)18O and the 19F(alpha,p)22Ne reaction rates. For these two reactions we present new estimates of the rates and their uncertainties. The importance of the partial mixing zone is reduced when using our estimate for the 14C(alpha,gamma)18O rate. Taking into account both the uncertainties related to the partial mixing zone and those related to nuclear reactions, the highest values of 19F enhancements observed in AGB stars are not matched by the models. This is a problem that will have to be revised by providing a better understanding of the formation and nucleosynthesis in the partial mixing zone, also in relation to reducing the uncertainties of the 14C(alpha,gamma)18O reaction rate. At the same time the possible effect of Cool Bottom Processing at the base of the convective envelope should be included in the computation of AGB nucleosynthesis. This process could in principle help matching the highest 19F abundances observed by decreasing the C/O ratio at the surface of the star, while leaving the 19F abundance unchanged.