Sub-Chandrasekhar CO white dwarfs accreting helium have been considered as candidates for Type Ia supernova(SNIa) progenitors since the early 1980s (helium shell mass $> 0.1 M_\odot $). These models, once detonated did not fit the observed spectra and light curve of typical SNIa observations. New theoretical work examined detonations on much less massive ($< 0.05 M_\odot $) envelopes. They find stable detonations that lead to light curves, spectra and abundances that compare relatively well with the observational data. The exact mechanism leading to the ignition of helium detonation is a key issue, since it is a mandatory first step for the whole scenario. As the flow of the accreted envelope is unstable to convection long before any hydrodynamic phenomena develops, a multidimensional approach is needed in order to study the ignition process. The complex convective reactive flow is challenging to any hydrodynamical solver. To the best of our knowledge, all previous 2D studies ignited the detonation artificially. We present here, for the first time, fully consistent results from two hydrodynamical 2D solvers that adopt two independent accurate schemes. For both solvers an effort was made to overcome the problematics raised by the finite resolution and numerical diffusion by the advective terms. Our best models lead to the ignition of a detonation in a convective cell. Our results are robust and the agreement between the two different numerical approaches is very good.

When carbon is ignited off-centre in a CO core of a super-AGB star, its burning in a convective shell tends to propagate to the centre. Whether the C flame will actually be able to reach the centre depends on the efficiency of extra mixing beneath the C convective shell. Whereas thermohaline mixing is too inefficient to interfere with the C-flame propagation, convective boundary mixing can prevent the C burning from reaching the centre. As a result, a C-O-Ne white dwarf (WD) is formed, after the star has lost its envelope. Such a "hybrid" WD has a small CO core surrounded by a thick ONe zone. In our 1D stellar evolution computations, the hybrid WD is allowed to accrete C-rich material, as if it were in a close binary system and accreted H-rich material from its companion with a sufficiently high rate at which the accreted H would be processed into He under stationary conditions, assuming that He could then be transformed into C. When the mass of the accreting WD approaches the Chandrasekhar limit, we find a series of convective Urca shell flashes associated with high abundances of 23Na and 25Mg. They are followed by off-centre C ignition leading to convection that occupies almost the entire star. To model the Urca processes, we use the most recent well-resolved data for their reaction and neutrino-energy loss rates. Because of the emphasized uncertainty of the convective Urca process in our hybrid WD models of SN Ia progenitors, we consider a number of their potentially possible alternative instances for different mixing assumptions, all of which reach a phase of explosive C ignition, either off or in the centre. Our hybrid SN Ia progenitor models have much lower C to O abundance ratios at the moment of the explosive C ignition than their pure CO counterparts, which may explain the observed diversity of the SNe Ia.

When carbon is ignited off-center in a CO core of a super-AGB star, its burning in a convective shell tends to propagate to the center. Whether the C flame will actually be able to reach the center depends on the efficiency of extra mixing beneath the C convective shell. Whereas thermohaline mixing is too inefficient to interfere with the C-flame propagation, convective boundary mixing can prevent the C burning from reaching the center. As a result, a C-O-Ne white dwarf (WD) is formed, after the star has lost its envelope. Such a "hybrid" WD has a small CO core surrounded by a thick ONe zone. In our 1D stellar evolution computations the hybrid WD is allowed to accrete C-rich material, as if it were in a close binary system and accreted H-rich material from its companion with a sufficiently high rate at which the accreted H would be processed into He under stationary conditions, assuming that He could then be transformed into C. When the mass of the accreting WD approaches the Chandrasekhar limit, we find a series of convective Urca shell flashes associated with high abundances of 23Na and 25Mg. They are followed by off-center C ignition leading to convection that occupies almost the entire star. To model the Urca processes, we use the most recent well-resolved data for their reaction and neutrino-energy loss rates. Because of the emphasized uncertainty of the convective Urca process in our hybrid WD models of SN Ia progenitors, we consider a number of their potentially possible alternative instances for different mixing assumptions, all of which reach a phase of explosive C ignition, either off or in the center. Our hybrid SN Ia progenitor models have much lower C to O abundance ratios at the moment of the explosive C ignition than their pure CO counterparts, which may explain the observed diversity of the SNe Ia.

Classical novae are the result of thermonuclear flashes of hydrogen accreted by CO or ONe white dwarfs, leading eventually to the dynamic ejection of the surface layers. These are observationally known to be enriched in heavy elements, such as C, O and Ne that must originate in layers below the H-flash convection zone. Building on our previous work, we now present stellar evolution simulations of ONe novae and provide a comprehensive comparison of our models with published ones. Some of our models include exponential convective boundary mixing to account for the observed enrichment of the nova ejecta even when accreted material has a solar abundance distribution. Our models produce maximum temperature evolution profiles and nucleosynthesis yields in good agreement with models that generate enriched ejecta by assuming that the accreted material was pre-mixed. We confirm for ONe novae the result we reported previously, i.e.\ we found that $^3$He could be produced {\it in situ} in solar-composition envelopes accreted with slow rates ($\dot{M} < 10^{-10}\,M_\odot/\mbox{yr}$) by cold ($T_{\rm WD} < 10^7$ K) CO WDs, and that convection was triggered by $^3$He burning before the nova outburst in that case. In addition, we now find that the interplay between the $^3$He production and destruction in the solar-composition envelope accreted with an intermediate rate, e.g.\ $\dot{M} = 10^{-10}\,M_\odot/\mbox{yr}$, by the $1.15\,M_\odot$ ONe WD with a relatively high initial central temperature, e.g.\ $T_{\rm WD} = 15\times 10^6$ K, leads to the formation of a thick radiative buffer zone that separates the bottom of the convective envelope from the WD surface. (Abridged)

A common model of the explosion mechanism of Type Ia supernovae is based on a delayed detonation of a white dwarf. A variety of models differ primarily in the method by which the deflagration leads to a detonation. A common feature of the models, however, is that all of them involve the propagation of the detonation through a white dwarf that is either expanding or contracting, where the stellar internal velocity profile depends on both time and space. In this work, we investigate the effects of the pre-detonation stellar internal velocity profile and the post-detonation velocity of expansion on the production of alpha-particle nuclei, including Ni56, which are the primary nuclei produced by the detonation wave. We perform one-dimensional hydrodynamic simulations of the explosion phase of the white dwarf for center and off-center detonations with five different stellar velocity profiles at the onset of the detonation. We observe two distinct post-detonation expansion phases: rarefaction and bulk expansion. Almost all the burning to Ni56 occurs only in the rarefaction phase, and its expansion time scale is influenced by pre-existing flow structure in the star, in particular by the pre-detonation stellar velocity profile. We find that the mass fractions of the alpha-particle nuclei, including Ni56, are tight functions of the empirical physical parameter rho_up/v_down, where rho_up is the mass density immediately upstream of the detonation wave front and v_down is the velocity of the flow immediately downstream of the detonation wave front. We also find that v_down depends on the pre-detonation flow velocity. We conclude that the properties of the pre-existing flow, in particular the internal stellar velocity profile, influence the final isotopic composition of burned matter produced by the detonation.

After off-center C ignition in the cores of super-AGB stars the C flame propagates all the way down to the center, trailing behind it the C-shell convective zone, and thus building a degenerate ONe core. This standard picture is obtained in stellar evolution simulations if the bottom C-shell convection boundary is assumed to be a discontinuity associated with a strict interpretation of the Schwarzschild condition for convective instability. However, this boundary is prone to additional mixing processes, such as thermohaline convection and convective boundary mixing. Using hydrodynamic simulations we show that, contrary to previous results, thermohaline mixing is too inefficient to interfere with the C-flame propagation. However, even a small amount of convective boundary mixing removes the physical conditions required for the C-flame propagation all the way to the center. This result holds even if we allow for some turbulent heat transport in the CBM region. As a result, super AGB stars build in their interiors hybrid C-O-Ne degenerate cores composed of a relatively large CO core (M_CO ~ 0.2 M_sun) surrounded by a thick ONe zone (M_ONe ~ 0.85 M_sun) with another thin CO layer above. If exposed by mass loss, these cores will become hybrid C-O-Ne white dwarfs. Otherwise, the ignition of C-rich material in the central core, surrounded by the thick ONe zone, may trigger a thermonuclear supernova explosion. The quenching of the C-flame may have implications for the ignition mechanism of SN Ia in the double-degenerate merger scenario.

Classical novae are the results of surface thermonuclear explosions of hydrogen accreted by white dwarfs (WDs) from their low-mass main-sequence or red-giant binary companions. Chemical composition analysis of their ejecta shows that nova outbursts occur on both carbon-oxygen (CO) and more massive oxygen-neon (ONe) WDs, and that there is cross-boundary mixing between the accreted envelope and underlying WD. We demonstrate that the state-of-the-art stellar evolution code MESA and post-processing nucleosynthesis tools of NuGrid can successfully be used for modeling of CO and ONe nova outbursts and nucleosynthesis. The convective boundary mixing (CBM) in our 1D numerical simulations is implemented using a diffusion coefficient that is exponentially decreasing with a distance below the bottom of the convective envelope. We show that this prescription produces maximum temperature evolution profiles and nucleosynthesis yields in good agreement with those obtained using the commonly adopted 1D nova model in which the CBM is mimicked by assuming that the accreted envelope has been pre-mixed with WD's material. In a previous paper, we have found that 3He can be produced in situ in solar-composition envelopes accreted with slow rates (dM/dt < 1e-10 M_sun/yr) by cold (T_WD < 1d7 K) CO WDs, and that convection is triggered by 3He burning before the nova outburst in this case. Here, we confirm this result for ONe novae. Additionally, we find that the interplay between the 3He production and destruction in the solar-composition envelope accreted with an intermediate rate, e.g. dM/dt = 1e-10 M_sun/yr, by the 1.15 M_sun ONe WD with a relatively high initial central temperature, e.g. T_WD = 15e6 K, leads to the formation of a thick radiative buffer zone that separates the bottom of the convective envelope from the WD surface.

Classical novae are the results of surface thermonuclear explosions of H-rich material accreted by white dwarfs (WDs) from their low-mass main-sequence or red-giant binary companions. Chemical composition analysis of their ejecta shows that nova outbursts occur on both carbon-oxygen (CO) and oxygen-neon (ONe) WDs, and that there is cross-boundary mixing between the accreted envelope and underlying WD. We have combined the stellar evolution code MESA and post-processing nucleosynthesis tools of NuGrid into a framework that allows to produce up-to-date models of nova outbursts and compute detailed nucleosynthesis in novae occurring on CO and ONe WDs. The convective boundary mixing (CBM) in our 1D numerical simulations is implemented using a diffusion coefficient that is exponentially decreasing with a distance below the bottom of the convective envelope. This MESA CBM prescription is based on the findings in 3D hydrodynamic simulations that the velocity field, and along with it the mixing expressed in terms of a diffusion coefficient, decays exponentially in the stable layer adjacent to a convective boundary. The framework can also use the commonly adopted 1D nova model in which the CBM is mimicked by assuming that the accreted envelope has been pre-enriched with WD's material. In this preliminary report, we present the most interesting new results related to CO and ONe nova outbursts that have been obtained with the Nova Framework.

We present here, for the first time, a 2D study of the overshoot convective mechanism in nova outbursts for a wide range of possible compositions of the layer underlying the accreted envelope. Previous surveys studied this mechanism only for solar composition matter accreted on top of carbon oxygen (C-O) white dwarfs. Since, during the runaway, mixing with carbon enhances the hydrogen burning rates dramatically, one should question whether significant enrichment of the ejecta is possible also for other underlying compositions (He, O, Ne, Mg), predicted by stellar evolution models. We simulated several non-carbon cases and found significant amounts of those underlying materials in the ejected hydrogen layer. Despite large differences in rates, time scales and energetics between the cases, our results show that the convective dredge up mechanism predicts significant enrichment in all our non-carbon cases, including helium enrichment in recurrent novae. The results are consistent with observations.

We present ultraviolet (UV) spectroscopy and photometry of four Type Ia supernovae (SNe 2004dt, 2004ef, 2005M, and 2005cf) obtained with the UV prism of the Advanced Camera for Surveys on the Hubble Space Telescope. This dataset provides unique spectral time series down to 2000 Angstrom. Significant diversity is seen in the near maximum-light spectra (~ 2000--3500 Angstrom) for this small sample. The corresponding photometric data, together with archival data from Swift Ultraviolet/Optical Telescope observations, provide further evidence of increased dispersion in the UV emission with respect to the optical. The peak luminosities measured in uvw1/F250W are found to correlate with the B-band light-curve shape parameter dm15(B), but with much larger scatter relative to the correlation in the broad-band B band (e.g., ~0.4 mag versus ~0.2 mag for those with 0.8 < dm15 < 1.7 mag). SN 2004dt is found as an outlier of this correlation (at > 3 sigma), being brighter than normal SNe Ia such as SN 2005cf by ~0.9 mag and ~2.0 mag in the uvw1/F250W and uvm2/F220W filters, respectively. We show that different progenitor metallicity or line-expansion velocities alone cannot explain such a large discrepancy. Viewing-angle effects, such as due to an asymmetric explosion, may have a significant influence on the flux emitted in the UV region. Detailed modeling is needed to disentangle and quantify the above effects.

Most of the leading explosion scenarios for Type Ia supernovae involve the nuclear incineration of a white dwarf star through a detonation wave. Several scenarios have been proposed as to how this detonation may actually occur, but the exact mechanism and environment in which it takes place remain unknown. We explore the effects of an off-center initiated detonation on the spatial distribution of the nucleosynthetic yield products in a toy model -- a pre-expanded near Chandrasekhar-mass white dwarf. We find that a single-point near edge-lit detonation results in asymmetries in the density and thermal profiles, notably the expansion timescale, throughout the supernova ejecta. We demonstrate that this asymmetry of the thermodynamic trajectories should be common to off-center detonations where a small amount of the star is burned prior to detonation. The sensitivity of the yields on the expansion timescale results in an asymmetric distribution of the elements synthesized as reaction products. We tabulate the shift in the center of mass of the various elements produced in our model supernova and find an odd-even pattern for elements past silicon. Our calculations show that off-center single-point detonations in carbon-oxygen white dwarfs are marked by significant composition asymmetries in their remnants which bear potentially observable signatures in both velocity and coordinate space, including an elemental nickel mass fraction which varies by a factor of two to three from one side of the remnant to the other.

There is a wide consensus in the astrophysics community that the mechanism underlying the observed Classical Nova eruptions is a surface thermonuclear runaway. We start this short review with the main observational facts that lead to the theoretical model of a thermonuclear runaway that takes place in an accreted hydrogen rich envelope placed on top of a cool degenerate core of a white dwarf. According to the theory, the accreted envelope becomes unstable to convection days to weeks prior to the runaway. During the extreme stages of the runaway itself, when the burning is most efficient, the envelope is fully convective. Therefore, the elements processed under such extreme conditions are lifted to the outermost regions of the star. A significant fraction of the envelope is ejected during the outburst. The complicated combination of hydrodynamic instabilities and explosive hydrogen burning, close to the surface of the star, gives us a unique opportunity to study this complex reactive flow. The range of core masses, core temperatures and accretion rates introduce a whole range of burning temperatures and densities. Following the description of the "standard" cases, we then focus on rare, but still possible, portions of the relevant parameter space, in which "breakout" of the traditional CNO cycle can occur and lead to heavy element enrichment patterns caused only by breakout burning. We conclude our review with the main challenges that nova theorists face today, with special emphasis on problems related to the nucleosynthesis issues.

The origin of the iron-group elements titanium to zinc in nature is understood to occur under explosive burning conditions in both Type Ia (thermonuclear) and Type II (core collapse) supernovae. In these dynamic environments, the most abundant product is found to be 56Ni ({\tau} = 8.5 days) that decays through 56Co ({\tau} = 111.5 days) to 56Fe. For the case of SNe Ia, the peak luminosities are proportional to the mass ejected in the form of 56Ni. It follows that the diversity of SNe Ia reflected in the range of peak luminosity provides a direct measure of the mass of 56Ni ejected. In this contribution, we identify and briefly discuss the factors that can influence the 56Ni mass and use both observations and theory to quantify their impact. We address specifically the variations in different stellar populations and possible distinctions with respect to SNe Ia progenitors.

Recently several type Ib supernovae (SNe; with the prototypical SN 2005E) have been shown to have atypical properties. These SNe are faint (absolute peak magnitude of ~ -15) and fast SNe that show unique composition. They are inferred to have low ejecta mass (a few tenths of a solar mass) and to be highly enriched in calcium, but poor in silicon elements and nickel. These SNe were therefore suggested to belong to a new class of calcium-rich faint SNe explosions. Their properties were proposed to be the result of helium detonations that may occur on helium accreting white dwarfs. In this paper we theoretically study the scenario of helium detonations, and focus on the results of detonations in accreted helium layers on low mass carbon-oxygen (CO) cores. We present new results from one dimensional simulations of such explosions, including their light curves and spectra. We find that when the density of the helium layer is low enough the helium detonation produces large amounts of intermediate elements, such as calcium and titanium, together with a large amount of unburnt helium. Our results suggest that the properties of calcium-rich faint SNe could indeed be consistent with the helium-detonation scenario on small CO cores. Above a certain density (larger CO cores) the detonation leaves mainly 56Ni and unburnt helium, and the predicted spectrum will unlikely fit the unique features of this class of SNe. Finally, none of our studied models reproduces the bright, fast evolving light curves of another type of peculiar SNe suggested to originate in helium detonations (SNe 1885A, 1939B and 2002bj).

We present Spitzer Space Telescope observations of the "peculiar variable" DZ Cru, identified by Rushton et al. (2008, MNRAS, 386, 289) as a classical nova. A dust shell, on which are superimposed a number of features, is prominent in the 5-35micron range some 4 years after eruption. We suggest that the dust in DZ Cru is primarily hydrogenated amorphous carbon in which aliphatic bands currently predominate, and which may either become predominantly aromatic as the dust is photo-processed by ultraviolet radiation from the stellar remnant, or more likely completely destroyed.

The astrophysical site of the r-process is still uncertain, and a full exploration of the systematics of this process in terms of its dependence on nuclear properties from stability to the neutron drip-line within realistic stellar environments has still to be undertaken. Sufficiently high neutron to seed ratios can only be obtained either in very neutron-rich low-entropy environments or moderately neutron-rich high-entropy environments, related to neutron star mergers (or jets of neutron star matter) and the high-entropy wind of core-collapse supernova explosions. As chemical evolution models seem to disfavor neutron star mergers, we focus here on high-entropy environments characterized by entropy $S$, electron abundance $Y_e$ and expansion velocity $V_{exp}$. We investigate the termination point of charged-particle reactions, and we define a maximum entropy $S_{final}$ for a given $V_{exp}$ and $Y_e$, beyond which the seed production of heavy elements fails due to the very small matter density. We then investigate whether an r-process subsequent to the charged-particle freeze-out can in principle be understood on the basis of the classical approach, which assumes a chemical equilibrium between neutron captures and photodisintegrations, possibly followed by a $\beta$-flow equilibrium. In particular, we illustrate how long such a chemical equilibrium approximation holds, how the freeze-out from such conditions affects the abundance pattern, and which role the late capture of neutrons originating from $\beta$-delayed neutron emission can play.

We study the initiation of the detonation in the gravitationally confined detonation (GCD) model of Type Ia supernovae (SNe Ia). Initiation of the detonation occurs spontaneously in a region where the length scale of the temperature gradient extending from a flow (in which carbon burning is already occurring) into unburned fuel is commensurate to the range of critical length scales which have been derived from 1D simulations that resolve the initiation of a detonation. By increasing the maximum resolution in a truncated cone that encompasses this region, beginning somewhat before initiation of the detonation occurs, we successfully simulate in situ the first gradient-initiated detonation in a whole-star simulation. The detonation emerges when a compression wave overruns a pocket of fuel situated in a Kelvin-Helmholtz cusp at the leading edge of the inwardly directed jet of burning carbon. The compression wave pre-conditions the temperature in the fuel in such a way that the Zel'dovich gradient mechanism can operate and a detonation ensues. We explore the dependence of the length scale of the temperature gradient on spatial resolution and discuss the implications for the robustness of this detonation mechanism. We find that the time and the location at which initiation of the detonation occurs varies with resolution. In particular, initiation of a detonation had not yet occurred in our highest resolution simulation by the time we ended the simulation because of the computational demand it required. We suggest that the turbulent shear layer surrounding the inwardly directed jet provides the most favorable physical conditions, and therefore the most likely location, for initiation of a detonation in the GCD model.

This white paper, directed to the Stars and Stellar Evolution panel, has three objectives: 1) to provide the Astro2010 Decadal Survey with a vista into the goals of the nuclear physics and nuclear astrophysics community; 2) to alert the astronomical community of joint opportunities for discoveries at the interface between nuclear physics and astronomy; and 3) to delineate efforts in nuclear physics and describe the observational and theoretical advances in astrophysics necessary to make progress towards answering the following questions in the Nuclear Science 2007 Long Range Plan: 1) What is the origin and distribution of the elements? 2) What are the nuclear reactions that power stars and stellar explosions? 3) What is the nature of dense matter? The scope of this white paper concerns the specific area of "low energy" nuclear astrophysics. We define this as the area of overlap between astrophysics and the study of nuclear structure and reactions. Of the questions listed above, two -- What is the origin of the elements? and What is the nature of dense matter? -- were specifically listed in the National Academies Study, "Connecting Quarks with the Cosmos".

Some explosion models for Type Ia supernovae (SN Ia), such as the gravitationally confined detonation (GCD) or the double detonation sub-Chandrasekhar (DDSC) models, rely on the spontaneous initiation of a detonation in the degenerate C/O material of a white dwarf. The length scales pertinent to the initiation of the detonation are notoriously unresolved in multi-dimensional stellar simulations, prompting the use of results of 1D simulations at higher resolution, such as the ones performed for this work, as guidelines for deciding whether or not conditions reached in the higher dimensional full star simulations successfully would lead to the onset of a detonation. Spontaneous initiation relies on the existence of a suitable gradient in self-ignition (induction) times of the fuel, which we set up with a spatially localized non-uniformity of temperature -- a hot spot. We determine the critical (smallest) sizes of such hot spots that still marginally result in a detonation in white dwarf matter by integrating the reactive Euler equations with the hydrodynamics code FLASH. We quantify the dependences of the critical sizes of such hot spots on composition, background temperature, peak temperature, geometry, and functional form of the temperature disturbance, many of which were hitherto largely unexplored in the literature. We discuss the implications of our results in the context of modeling of SNe Ia.

While the high-entropy wind (HEW) of Type II supernovae remains one of the more promising sites for the rapid neutron-capture (r-) process, hydrodynamic simulations have yet to reproduce the astrophysical conditions under which the latter occurs. We have performed large-scale network calculations within an extended parameter range of the HEW, seeking to identify or to constrain the necessary conditions for a full reproduction of all r-process residuals N_{r,\odot}=N_{\odot}-N_{s,\odot} by comparing the results with recent astronomical observations. A superposition of weighted entropy trajectories results in an excellent reproduction of the overall N_{r,\odot}-pattern beyond Sn. For the lighter elements, from the Fe-group via Sr-Y-Zr to Ag, our HEW calculations indicate a transition from the need for clearly different sources (conditions/sites) to a possible co-production with r-process elements, provided that a range of entropies are contributing. This explains recent halo-star observations of a clear non-correlation of Zn and Ge and a weak correlation of Sr - Zr with heavier r-process elements. Moreover, new observational data on Ru and Pd seem to confirm also a partial correlation with Sr as well as the main r-process elements (e.g. Eu).

For very slow white dwarf accretors in CV's Townsley and Bildsten (2004) found a relation between the accretion rate and the central temperature of the white dwarf Tc. According to this relation for accretion rates less than 10^-10 solar masses per year Tc is much lower than 10^7 K. Motivated by this study we follow the thermonuclear runaway on massive white dwarfs (M_WD=1.25 - 1.40 solar masses) with Tc lower than 10^7 K, accreting matter of solar composition. We demonstrate that in that range of the relevant parameter space (Tc,M_WD and accretion rate) the slope of the relation between the peak temperatures achieved during the runaway and Tc becomes much steeper than its value for Tc above 10^7 K. The peak temperatures we derive can lead to nuclear breakout from the conventional "hot carbon-nitrogen-oxygen" cycle. When breakout conditions are achieved the heavy element abundances can show a much wider variety than what is possible with the common enrichment mechanisms.

Proton-rich material in a state of nuclear statistical equilibrium (NSE) is one of the least studied regimes of nucleosynthesis. One reason for this is that after hydrogen burning, stellar evolution proceeds at conditions of equal number of neutrons and protons or at a slight degree of neutron-richness. Proton-rich nucleosynthesis in stars tends to occur only when hydrogen-rich material that accretes onto a white dwarf or neutron star explodes, or when neutrino interactions in the winds from a nascent proto-neutron star or collapsar-disk drive the matter proton-rich prior to or during the nucleosynthesis. In this paper we solve the NSE equations for a range of proton-rich thermodynamic conditions. We show that cold proton-rich NSE is qualitatively different from neutron-rich NSE. Instead of being dominated by the Fe-peak nuclei with the largest binding energy per nucleon that have a proton to nucleon ratio close to the prescribed electron fraction, NSE for proton-rich material near freeze-out temperature is mainly composed of Ni56 and free protons. Previous results of nuclear reaction network calculations rely on this non-intuitive high proton abundance, which this paper will explain. We show how the differences and especially the large fraction of free protons arises from the minimization of the free energy as a result of a delicate competition between the entropy and the nuclear binding energy.

We study the gravitationally confined detonation (GCD) model of Type Ia supernovae through the detonation phase and into homologous expansion. In the GCD model, a detonation is triggered by the surface flow due to single point, off-center flame ignition in carbon-oxygen white dwarfs. The simulations are unique in terms of the degree to which non-idealized physics is used to treat the reactive flow, including weak reaction rates and a time dependent treatment of material in nuclear statistical equilibrium (NSE). Careful attention is paid to accurately calculating the final composition of material which is burned to NSE and frozen out in the rapid expansion following the passage of a detonation wave over the high density core of the white dwarf; and an efficient method for nucleosynthesis post-processing is developed which obviates the need for costly network calculations along tracer particle thermodynamic trajectories. Observational diagnostics are presented for the explosion models, including abundance stratifications and integrated yields. We find that for all of the ignition conditions studied here, a self regulating process comprised of neutronization and stellar expansion results in final \iso{Ni}{56} masses of $\sim$1.1\msun. But, more energetic models result in larger total NSE and stable Fe peak yields. The total yield of intermediate mass elements is $\sim0.1$\msun and the explosion energies are all around 1.5$\times10^{51}$ ergs. The explosion models are briefly compared to the inferred properties of recent Type Ia supernova observations. The potential for surface detonation models to produce lower luminosity (lower \iso{Ni}{56} mass) supernovae is discussed.

We report the results of a series of three-dimensional (3-D) simulations of the deflagration phase of the gravitationally confined detonation mechanism for Type Ia supernovae. In this mechanism, ignition occurs at one or several off-center points, resulting in a burning bubble of hot ash that rises rapidly, breaks through the surface of the star, and collides at a point opposite breakout on the stellar surface. We find that detonation conditions are robustly reached in our 3-D simulations for a range of initial conditions and resolutions. Detonation conditions are achieved as the result of an inwardly-directed jet that is produced by the compression of unburnt surface material when the surface flow collides with itself. A high-velocity outwardly-directed jet is also produced. The initial conditions explored in this paper lead to conditions at detonation that can be expected to produce large amounts of $^{56}$Ni and small amounts of intermediate mass elements. These particular simulations are therefore relevant only to high luminosity Type Ia supernovae. Recent observations of Type Ia supernovae imply a compositional structure that is qualitatively consistent with that expected from these simulations.

We present ground-based SpectroCam-10 mid-infrared, MMT optical, and Spitzer Space Telescope IRS mid-infrared spectra taken 7.62, 18.75, and 19.38 years respectively after the outburst of the old classical nova QU Vulpeculae (Nova Vul 1984 #2). The spectra of the ejecta are dominated by forbidden line emission from neon and oxygen. Our analysis shows that neon was, at the first and last epochs respectively, more than 76 and 168 times overabundant by number with respect to hydrogen compared to the solar value. These high lower limits to the neon abundance confirm that QU Vul involved a thermonuclear runaway on an ONeMg white dwarf and approach the yields predicted by models of the nucleosynthesis in such events.

We develop an improved method for tracking the nuclear flame during the deflagration phase of a Type Ia supernova, and apply it to study the variation in outcomes expected from the gravitationally confined detonation (GCD) paradigm. A simplified 3-stage burning model and a non-static ash state are integrated with an artificially thickened advection-diffusion-reaction (ADR) flame front in order to provide an accurate but highly efficient representation of the energy release and electron capture in and after the unresolvable flame. We demonstrate that both our ADR and energy release methods do not generate significant acoustic noise, as has been a problem with previous ADR-based schemes. We proceed to model aspects of the deflagration, particularly the role of buoyancy of the hot ash, and find that our methods are reasonably well-behaved with respect to numerical resolution. We show that if a detonation occurs in material swept up by the material ejected by the first rising bubble but gravitationally confined to the white dwarf (WD) surface (the GCD paradigm), the density structure of the WD at detonation is systematically correlated with the distance of the deflagration ignition point from the center of the star. Coupled to a suitably stochastic ignition process, this correlation may provide a plausible explanation for the variety of nickel masses seen in Type Ia Supernovae.

We report optical photometry and optical through mid-infrared spectroscopy of the classical nova V1186 Sco. This slowly developing nova had an complex light curve with multiple secondary peaks similar to those seen in PW Vul. The time to decline 2 magnitudes, t$_2$, was 20 days but the erratic nature of the light curve makes determination of intrinsic properties based on the decline time (e.g., luminosity) problematic, and the often cited MMRD relationship of Della Valle and Livio (1995) fails to yield a plausible distance. Spectra covering 0.35 to 35 $\mu$m were obtained in two separate epochs during the first year of outburst. The first set of spectra, taken about 2 months after visible maximum, are typical of a CO-type nova with narrow line emission from \ion{H}{1}, \ion{Fe}{2}, \ion{O}{1} and \ion{He}{1}. Later data, obtained between 260 and 380 days after maximum, reveal an emerging nebular spectrum. \textit{Spitzer} spectra show weakening hydrogen recombination emission with the emergence of [\ion{Ne}{2}] (12.81 $\mu$m) as the strongest line. Strong emission from [\ion{Ne}{3}] (15.56 $\mu$m) is also detected. Photoionization models with low effective temperature sources and only marginal neon enhancement (Ne $\sim$ 1.3 Ne$_{\odot}$) are consistent with these IR fine-structure neon lines indicating that V1186 Sco did not occur on a ONeMg white dwarf. In contrast, the slow and erratic light curve evolution, spectral development, and photoionization analysis of the ejecta imply the outburst occurred on a low mass CO white dwarf. We note that this is the first time strong [\ion{Ne}{2}] lines have been detected so early in the outburst of a CO nova and suggests that the presence of mid-infrared neon lines is not directly indicative of a ONeMg nova event.

Abundances of heavier elements (barium and beyond) in many neutron-capture-element-rich halo stars accurately replicate the solar system r-process pattern. However, abundances of lighter neutron-capture elements in these stars are not consistent with the solar system pattern. These comparisons suggest contributions from two distinct types of r-process synthesis events -- a so called main r-process for the elements above the 2nd r-process peak and a weak r-process for the lighter neutron-capture elements. We have performed r-process theoretical predictions to further explore the implications of the solar and stellar observations. We find that the isotopic composition of barium and the elemental Ba/Eu abundance ratios in r-process-rich low metallicity stars can only be matched by computations in which the neutron densities are in the range 23< log n_n < 28, values typical of the main r-process. For r-process conditions that successfully generate the heavy element pattern extending down to A=135, the relative abundance of I129 produced in this mass region appears to be at least 90% of the observed solar value. Finally, in the neutron number density ranges required for production of the observed solar/stellar 3rd r-process-peak (A~200), the predicted abundances of inter-peak element hafnium (Z=72, A~180) follow closely those of 3rd-peak elements and lead. Hf, observable from the ground and close in mass number to the 3rd r-process-peak elements, might also be utilized as part of a new nuclear chronometer pair, Th/Hf, for stellar age determinations.

Al Cameron, who died recently (October 3, 2005) at 80, was one of the giants in astrophysics. His insights were profound and his interests were wide-ranging. Originally trained as a nuclear physicist, he made major contributions in a number of fields, including nuclear reactions in stars, nucleosynthesis, the abundances of the elements in the Solar System, and the origin of the Solar System and the Moon. In 1957, Cameron and, independently, Burbidge, Burbidge, Fowler and Hoyle, wrote seminal papers on nuclear astrophysics. Most of our current ideas concerning element formation in stars have followed from those two pioneering and historical works. Al also made many contributions in the field of Solar System physics. Particularly noteworthy in this regard was Cameron's work on the formation of the Moon. Al was also a good friend and mentor of young people. Al Cameron will be missed by many in the community both for his scientific contributions and for his friendship.

Neutron stars, with their strong surface gravity, have interestingly short timescales for the sedimentation of heavy elements. Recent observations of unstable thermonuclear burning (observed as X-ray bursts) on the surfaces of slowly accreting neutron stars ($< 0.01$ of the Eddington rate) motivate us to examine how sedimentation of CNO isotopes affects the ignition of these bursts. We further estimate the burst development using a simple one-zone model with a full reaction network. We report a region of mass accretion rates for weak H flashes. Such flashes can lead to a large reservoir of He, the unstable burning of which may explain some observed long bursts (duration $\sim 1000$ s).

We develop and calibrate a realistic model flame for hydrodynamical simulations of deflagrations in white dwarf (Type Ia) supernovae. Our flame model builds on the advection-diffusion-reaction model of Khokhlov and includes electron screening and Coulomb corrections to the equation of state in a self-consistent way. We calibrate this model flame--its energetics and timescales for energy release and neutronization--with self-heating reaction network calculations that include both these Coulomb effects and up-to-date weak interactions. The burned material evolves post-flame due to both weak interactions and hydrodynamic changes in density and temperature. We develop a scheme to follow the evolution, including neutronization, of the NSE state subsequent to the passage of the flame front. As a result, our model flame is suitable for deflagration simulations over a wide range of initial central densities and can track the temperature and electron fraction of the burned material through the explosion and into the expansion of the ejecta.

Neutron stars, with their strong surface gravity, have interestingly short timescales for the sedimentation of heavy elements. Motivated by observations of Type I X-ray bursts from sources with extremely low persistent accretion luminosities, $L_X < 10^{36}\usp\ergspersecond (\simeq 0.01\ensuremath{L_{\mathrm{Edd}}}$), we study how sedimentation affects the distribution of isotopes and the ignition of H and He in the envelope of an accreting neutron star. For local mass accretion rates $\mdot \lesssim 10^{-2}\medd$ (for which the ignition of H is unstable), where $\medd = 8.8\times 10^{4}\nsp\gpscps$, the helium and CNO elements sediment out of the accreted fuel before reaching a temperature where H would ignite. Using one-zone calculations of the thermonuclear burning, we find a range of accretion rates for which the unstable H ignition does not trigger unstable He burning. This range depends on the emergent flux from reactions in the deep neutron star crust; for $F = 0.1\nsp\MeV(\dot{m}/\mb)$, the range is $3\times 10^{-3}\medd\lesssim\mdot\lesssim 10^{-2}\medd$. We speculate that sources accreting in this range will build up a massive He layer that later produces an energetic and long X-ray burst. At mass accretion rates lower than this range, we find that the H flash leads to a strong mixed H/He flash. Surprisingly, even at accretion rates $\mdot \gtrsim 0.1\medd$, although the H and He do not completely segregate, the H abundance at the base of the accumulated layer is still reduced. While following the evolution of the X-ray burst is beyond the scope of this introductory paper, we note that the reduced proton-to-seed ratio favors the production of \iso{12}{C}--an important ingredient for subsequent superbursts.

Since the pioneering observations of Spite & Spite in 1982, the constant lithium abundance of metal-poor ([Fe/H]<-1.3) halo stars near the turnoff has been attributed a cosmological origin. Closer analysis revealed that the observed abundance lies at $\Delta$ Li~0.4 dex below the predictions of Big Bang Nucleosynthesis. The measurements of deuterium abundances on the lines of sight toward quasars and the recent data from the Wilkinson Microwave Anisotropy Probe independently confirmed this gap. We suggest here that part of the discrepancy is explained by the first generation of stars that efficiently depleted lithium. Assuming that the models for lithium evolution in the halo turnoff stars and $\Delta$ Li estimates are correct, we infer that between 1/3 and 1/2 of the baryonic matter of the early halo (~10^9 Mo) was processed through Population III stars. This new paradigm proposes a very economical solution to the lingering difficulty of understanding the properties of the Spite Plateau and its lack of star-to-star scatter down to [Fe/H]=-2.5. It is moreover in agreement both with the absence of lithium detection in the most metal-poor star presently known (HE 1327-2326), and also with new trends of the Plateau suggesting its low metallicity edge may be reached around [Fe/H]=-2.5. We discuss the turbulent mixing associated with enhanced supernovae explosions in the early interstellar medium in this picture. Finally we show that other chemical properties of the extremely metal-poor stars are in agreement with a significant Population III processing in the halo, provided these models include mass-loss and rotationally-induced mixing.

The damped Lyman-alpha (DLA) and sub-DLA quasar absorption lines provide powerful probes of the evolution of metals, gas, and stars in galaxies. One major obstacle in trying to understand the evolution of DLAs and sub-DLAs has been the small number of metallicity measurements at z < 1.5, an epoch spanning \~70 % of the cosmic history. In recent surveys with the Hubble Space Telescope and Multiple Mirror Telescope, we have doubled the DLA Zn sample at z < 1.5. Combining our results with those at higher redshifts from the literature, we find that the global mean metallicity of DLAs does not rise to the solar value at low redshifts. These surprising results appear to contradict the near-solar mean metallicity observed for nearby (z ~ 0) galaxies and the predictions of cosmic chemical evolution models based on the global star formation history. Finally, we discuss direct constraints on the star formation rates (SFRs) in the absorber galaxies from our deep Fabry-Perot Ly-alpha imaging study and other emission-line studies in the literature. A large fraction of the observed heavy-element quasar absorbers at 0 < z < 3.4 appear to have SFRs substantially below the global mean SFR, consistent with the low metallicities observed in the spectroscopic studies.

Multidimensional reactive flow models of accreted hydrogen rich envelopes on top of degenerate cold white dwarfs are very effective tools for the study of critical, non spherically symmetric, behaviors during the early stages of nova outbursts. Such models can shed light both on the mechanism responsible for the heavy element enrichment observed to characterize nova envelope matter and on the role of perturbations during the early stages of ignition of the runaway. The complexity of convective reactive flow in multi-dimensions makes the computational model itself complex and sensitive to the details of the numerics. In this study, we demonstrate that the imposed outer boundary condition can have a dramatic effect on the solution. Several commonly used choices for the outer boundary conditions are examined. It is shown that the solutions obtained from Lagrangian simulations, where the envelope is allowed to expand and mass is being conserved, are consistent with spherically symmetric solutions. In Eulerian schemes which utilize an outer boundary condition of free outflow, the outburst can be artificially quenched.

We present new abundance determinations of neutron-capture elements Ge, Zr, Os, Ir, and Pt in a sample of 11 metal-poor (-3.1 <= [Fe/H] <= -1.6) Galactic halo giant stars, based on Hubble Space Telescope UV and Keck I optical high-resolution spectroscopy. The stellar sample is dominated by r-process-rich stars such as the well-studied CS 22892-052 and bd+173248, but also includes the r-process-poor, bright giant HD 122563. Our results demonstrate that abundances of the 3rd r-process peak elements Os, Ir and Pt in these metal-poor halo stars are very well-correlated among themselves, and with the abundances of the canonical r-process element Eu (determined in other studies), thus arguing for a common origin or site for r-process nucleosynthesis of heavier (Z>56) elements. However, the large (and correlated) scatters of [Eu,Os,Ir,Pt/Fe] suggests that the heaviest neutron-capture r-process elements are not formed in all supernovae. In contrast, the Ge abundances of all program stars track their Fe abundances, very well. An explosive process on iron-peak nuclei (e.g., the alpha-rich freeze-out in supernovae), rather than neutron capture, appears to have been the dominant synthesis mechanism for this element at low metallicities -- Ge abundances seem completely uncorrelated with Eu.

Most models of cosmic chemical evolution predict that the mass-weighted mean interstellar metallicity of galaxies should rise with time from a low value $\sim 1/30$ solar at $z \sim 3$ to a nearly solar value at $z = 0$. In the absence of any selection effects, the damped Lyman-alpha absorbers (DLAs) in quasar spectra are expected to show such a rise in global metallicity. However, it has been difficult to determine whether or not DLAs show this effect, primarily because of the very small number of DLA metallicity measurements at low redshifts. In an attempt to put tighter constraints on the low-redshift end of the DLA metallicity-redshift relation, we have observed Zn II and Cr II lines in four DLAs at $0.09 < z < 0.52$, using the Space Telescope Imaging Spectrograph (STIS) onboard the Hubble Space Telescope (HST). These observations have provided the first constraints on Zn abundances in DLAs with $z < 0.4$. In all the three DLAs for which our observations offer meaningful constraints on the metallicity, the data suggest that the metallicities are much lower than the solar value. These results are consistent with recent imaging studies indicating that these DLAs may be associated with dwarf or low surface brightness galaxies. We combine our results with higher redshift data from the literature to estimate the global mean metallicity-redshift relation for DLAs. We find that the global mean metallicity shows at most a slow increase with decreasing redshift. ...(Please see the paper for the complete abstract).

We present the first high-resolution three-dimensional simulations of the deflagration phase of Type Ia supernovae that treat the entire massive white dwarf. We report the results of simulations in which ignition of the nuclear burning occurs slightly off-center. The subsequent evolution of the nuclear burning is surprisingly asymmetric with a growing bubble of hot ash rapidly rising to the stellar surface. Upon reaching the surface, the mass of burned material is $\approx 0.075 M_\sun$ and the kinetic energy is $4.3 \times 10^{49}$ ergs. The velocity of the top of the rising bubble approaches 8000 km s$^{-1}$. The amount of the asymmetry found in the model offers a natural explanation for the observed diversity in Type Ia supernovae. Our study strongly disfavors the classic central-ignition pure deflagration scenario by showing that the result is highly sensitive to details of the initial conditions.

Many classical nova ejecta are enriched in CNO and Ne. Rosner et al. recently suggested that the enrichment might originate in the resonant interaction between large-scale shear flows in the accreted H/He envelope and gravity waves at the interface between the envelope and the underlying C/O white dwarf. The shear flow amplifies the waves, which eventually form cusps and break. This wave breaking injects a spray of C/O into the superincumbent H/He. Using two-dimensional simulations, we formulate a quantitative expression for the amount of C/O per unit area that can be entrained into the H/He at saturation. The fraction of the envelope that is enriched depends on the horizontal distribution of shear velocity and the density contrast between the C/O white dwarf and the H/He layer but is roughly independent of the vertical shape of the shear profile. Using this parameterization for the mixed mass, we then perform several one-dimensional Lagrangian calculations of an accreting white dwarf envelope and consider two scenarios: that the wave breaking and mixing is driven by the convective flows; and that the mixing occurs prior to the onset of convection. In the absence of enrichment prior to ignition, the base of the convective zone, as calculated from mixing-length theory with the Ledoux instability criterion, does not reach the C/O interface. As a result, there is no additional mixing, and the runaway is slow. In contrast, the formation of a mixed layer during the accretion of H/He, prior to ignition, causes a more violent runaway. The envelope can be enriched by <25% of C/O by mass (consistent with that observed in some ejecta) for shear velocities, over the surface, with Mach numbers <0.4.

Recent Chandra and XMM-Newton observations of galaxy cluster cooling flows have revealed X-ray emission voids of up to 30 kpc in size that have been identified with buoyant, magnetized bubbles. Motivated by these observations, we have investigated the behavior of rising bubbles in stratified atmospheres using the FLASH adaptive-mesh simulation code. We present results from 2-D simulations with and without the effects of magnetic fields, and with varying bubble sizes and background stratifications. We find purely hydrodynamic bubbles to be unstable; a dynamically important magnetic field is required to maintain a bubble's integrity. This suggests that, even absent thermal conduction, for bubbles to be persistent enough to be regularly observed, they must be supported in large part by magnetic fields. Thermal conduction unmitigated by magnetic fields can dissipate the bubbles even faster. We also observe that the bubbles leave a tail as they rise; the structure of these tails can indicate the history of the dynamics of the rising bubble.

We investigate whether a single population of first stars could have influenced both the metal enrichment and reionization of the high-redshift intergalactic medium (IGM), by calculating the generated ionizing radiation per unit metal yield as a function of the metallicity of stellar populations. We examine the relation between the ionizing radiation and carbon created by the first stars, since the evidence for the widespread enrichment of the IGM at redshifts $z$ about 3-4 comes from the detection of C IV absorption. We find that the number of ionizing photons per baryon generated in association with the detected IGM metallicity may well exceed that required for a late hydrogen reionization at $z$ of about 6, by up to a factor of 10-20 for metal-free stars in a present-day initial mass function (IMF). This would be in agreement with similar indications from recent observations of the microwave background and the high-$z$ IGM. In addition, the contribution from intermediate-mass stars to the total metal yield, neglected in past works, substantially impacts such calculations. Lastly, a top-heavy IMF is not necessarily preferred as a more efficient high-$z$ source of ionizing radiation, based on nucleosynthetic arguments in association with a given level of IGM enrichment.

Critically understanding the `standard candle'-like behavior of Type Ia supernovae requires understanding their explosion mechanism. One family of models for Type Ia Supernovae begins with a deflagration in a Carbon-Oxygen white dwarf which greatly accelerates through wrinkling and flame instabilities. While the planar speed and behavior of astrophysically-relevant flames is increasingly well understood, more complex behavior, such as the flame's response to stretch and curvature, has not been extensively explored in the astrophysical literature; this behavior can greatly enhance or suppress instabilities and local flame-wrinkling, which in turn can increase or decrease the bulk burning rate. In this paper, we explore the effects of curvature on both nuclear flames and simpler model flames to understand the effect of curvature on the flame structure and speed.

In this paper, we revisit the problem of the determination of the frequency of occurrence of galactic nova outbursts which involve an oxygen-neon (ONe) white dwarf. The improvement with respect to previous work on the subject derives from the fact that we use the results that our evolutionary calculations provide for the final mass and for the chemical profiles of intermediate-to-massive primary components of close binary systems. In particular, the final evolutionary stages, such as the carbon burning phase, have been carefully followed for the whole range of masses of interest. The chemical profiles obtained with our evolutionary code are of interest in determining the chemical composition of the ejecta after being processed through the thermonuclear runaway, although such other factors as the efficiency of the mixing between the accreted material and that of the underlying white dwarf must also be considered. In our calculations of the frequency of occurrence of nova outbursts involving an ONe white dwarf, we also take into account the observational selection effects introduced by the different recurrence times of the outbursts and by the spatial distribution of novae. In spite of the very different evolutionary sequences, we find that approximately 1/3 of the novae observed in outburst should involve an oxygen-neon white dwarf, in agreement with previous theoretical estimates.

We explore the idea that the observed variations in the peak luminosities of Type Ia supernovae originate in part from a scatter in metallicity of the main-sequence stars that become white dwarfs. Previous, numerical, studies have not self-consistently explored metallicities greater than solar. One-dimensional Chandrasekhar mass models of SNe Ia produce most of their 56Ni in a burn to nuclear statistical equilibrium between the mass shells 0.2 and 0.8 solar masses, for which the electron to nucleon ratio is constant during the burn. We show analytically that, under these conditions, charge and mass conservation constrain the mass of 56Ni produced to depend linearly on the original metallicity of the white dwarf progenitor. Detailed post-processing of W7-like models confirms this linear dependence. The effect that we identify is most evident at metallicities larger than solar, and is in agreement with previous self-consistent calculations over the metallicity range common to both calculations. The observed scatter in the metallicity (1/3--3 times solar) of the solar neighborhood is enough to induce a 25% variation in the mass of 56Ni ejected by Type Ia supernovae. This is sufficient to vary the peak V-band brightness by approximately 0.2. This scatter in metallicity is present out to the limiting redshifts of current observations (z < 1). Sedimentation of 22Ne can possibly amplify the variation in 56Ni mass up to 50%. Further numerical studies can determine if other metallicity-induced effects, such as a change in the mass of the 56Ni-producing region, offset or enhance this variation.

High-resolution spectra obtained with three ground-based facilities and the Hubble Space Telescope (HST) have been combined to produce a new abundance analysis of CS 22892-052, an extremely metal-poor giant with large relative enhancements of neutron-capture elements. A revised model stellar atmosphere has been derived with the aid of a large number of Fe-peak transitions, including both neutral and ionized species of six elements.Several elements, including Mo, Lu, Au, Pt and Pb, have been detected for the first time in CS 22892-052, and significant upper limits have been placed on the abundances of Ga, Ge, Cd, Sn, and U in this star. In total, abundance measurements or upper limits have been determined for 57 elements, far more than previously possible. New Be and Li detections in CS 22892-052 indicate that the abundances of both these elements are significantly depleted compared to unevolved main-sequence turnoff stars of similar metallicity. Abundance comparisons show an excellent agreement between the heaviest n-capture elements (Z >= 56) and scaled solar system r-process abundances, confirming earlier results for CS 22892-052 and other metal-poor stars. New theoretical r-process calculations also show good agreement with CS 22892-052 abundances as well as the solar r-process abundance components.The abundances of lighter elements (40<= Z <= 50), however, deviate from the same scaled abundance curves that match the heavier elements, suggesting different synthesis conditions or sites for the low-mass and high-mass ends of the abundance distribution. The detection of Th and the upper limit on the U abundance together imply a lower limit of 10.4 Gyr on the age of CS 22892-052, quite consistent with the Th/Eu age estimate of 12.8 +/- ~= 3 Gyr. An average of several chronometric ratios yields an age 14.2 +/- ~= 3 Gyr.

V1494 Aql (Nova Aql 1999 No. 2) was discovered on 2 December 1999. We obtained Chandra ACIS-I spectra on 15 April and 7 June 2000 which appear to show only emission lines. Our third observation, on 6 August, showed that its spectrum had evolved to that characteristic of a Super Soft X-ray Source. We then obtained Chandra LETG+HRC-S spectra on 28 September (8 ksec) and 1 October (17 ksec). We analyzed the X-ray light curve of our grating observations and found both a short time scale ``burst'' and oscillations. Neither of these phenomena have previously been seen in the light curve of a nova in outburst. The ``burst'' was a factor of 10 rise in X-ray counts near the middle of the second observation, and which lasted about 1000 sec; it exhibited at least two peaks, in addition to other structure. Our time series analysis of the combined 25 ksec observation shows a peak at 2500 s which is present in independent analyses of both the zeroth order image and the dispersed spectrum and is not present in similar analyses of grating data for HZ 43 and Sirius B. Further analyses of the V1494 Aql data find other periods present which implies that we are observing non-radial g+ modes from the pulsating, rekindled white dwarf.

The heavy elements formed by neutron capture processes have an interesting history from which we can extract useful clues to and constraints upon both the characteristics of the processes themselves and the star formation and nucleosynthesis history of Galactic matter. Of particular interest in this regard are the heavy element compositions of extremely metal-deficient stars. At metallicities [Fe/H] <= -2.5, the elements in the mass region past barium (A >= 130-140 have been found (in non carbon-rich stars) to be pure r-process products. The identification of an environment provided by massive stars and associated Type II supernovae as an r-process site seems compelling. Increasing levels of heavy s-process (e.g., barium) enrichment with increasing metallicity, evident in the abundances of more metal-rich halo stars and disk stars, reflect the delayed contributions from the low- and intermediate-mass (M \~ 1-3 Msol) stars that provide the site for the main s-process nucleosynthesis component during the AGB phase of their evolution. New abundance data in the mass region 60 <~ A <~ 130 is providing insight into the identity of possible alternative r-process sites. We review recent observational studies of heavy element abundances both in low metallicity halo stars and in disk stars, discuss the observed trends in light of nucleosynthesis theory, and explore some implications of these results for Galactic chemical evolution, nucleosynthesis, and nucleocosmochronology.

We look in detail at the process of mapping an astrophysical initial model from a stellar evolution code onto the computational grid of an explicit, Godunov type code while maintaining hydrostatic equilibrium. This mapping process is common in astrophysical simulations, when it is necessary to follow short-timescale dynamics after a period of long timescale buildup. We look at the effects of spatial resolution, boundary conditions, the treatment of the gravitational source terms in the hydrodynamics solver, and the initialization process itself. We conclude with a summary detailing the mapping process that yields the lowest ambient velocities in the mapped model.

Observed novae abundances and explosion energies estimated from observations indicate that there must be significant mixing of the heavier material of the white dwarf (C+O) into the lighter accreted material (H+He). Accordingly, nova models must incorporate a mechanism that will dredge up the heavier white dwarf material, and fluid motions from an early convection phase is one proposed mechanism. We present results from two-dimensional simulations of classical nova precursor models that demonstrate the beginning of a convective phase during the `simmering' of a Nova precursor. We use a new hydrostatic equilibrium hydrodynamics module recently developed for the adaptive-mesh code FLASH. The two-dimensional models are based on the one-dimensional models of Ami Glasner, and were evolved with FLASH from a pre-convective state to the onset of convection.

We present the results of a simulation of a wind-driven non-linear gravity wave breaking on the surface of a white dwarf. The ``wind'' consists of H/He from an accreted envelope, and the simulation demonstrates that this breaking wave mechanism can produce a well-mixed layer of H/He with C/O from the white dwarf above the surface. Material from this mixed layer may then be transported throughout the accreted envelope by convection, which would enrich the C/O abundance of the envelope as is expected from observations of novae.

We present detailed abundance observations of the metal-poor ([Fe/H] = --2.0) Galactic halo star BD+17 3248. We have detected gold in this star, for the first time in any halo star. We also detected a weak line at 3859.60 A in the spectrum of BD+17 3248 that we tentatively identify with uranium. This would be only the second such detection of this radioactive element in any metal-poor halo stars. Similarly to the case of CS 22892--052, the heavier n-capture elements (Z >= 56) fall on the scaled solar system r-process abundance curve, while the elements with Z = 40--50 in BD+17 3248 fall below that same solar curve. Employing the newly detected Th, U and third r-process peak element (e.g., Pt, Ir) abundances, chronometric age estimates have been made for BD+17 3248. The average value of the various chronometric pairs suggests an age of 13.8 +/- 4 Gyr for this star, consistent, within error limits, with other chronometric age determinations for metal-poor Galactic halo stars.

We have combined new high-resolution spectra obtained with the Hubble Space Telescope (HST) and ground-based facilities to make a comprehensive new abundance analysis of the metal-poor, halo star BD +17^\circ 3248. We have detected the third r-process peak elements osmium, platinum, and (for the first time in a metal-poor star) gold, elements whose abundances can only be reliably determined using HST. Our observations illustrate a pattern seen in other similar halo stars with the abundances of the heavier neutron-capture elements, including the third r-process peak elements, consistent with a scaled solar system r-process distribution. The abundances of the lighter neutron-capture elements, including germanium and silver, fall below that same scaled solar r-process curve, a result similar to that seen in the ultra-metal-poor star CS 22892--052. A single site with two regimes or sets of conditions, or perhaps two different sites for the lighter and heavier neutron-capture elements, might explain the abundance pattern seen in this star. In addition we have derived a reliable abundance for the radioactive element thorium. We tentatively identify U II at 3859 A in the spectrum of BD +17^\circ 3248, which makes this the second detection of uranium in a very metal-poor halo star. Our combined observations cover the widest range in proton number (from germanium to uranium) thus far of neutron-capture elements in metal-poor Galactic halo stars. Employing the thorium and uranium abundances in comparison with each other and with several stable elements, we determine an average cosmochronological age for BD +17^\circ 3248 of 13.8 +/- 4 Gyr, consistent with that found for other similar metal-poor halo stars.

Using the results of recent work in shear instabilities in stratified fluids, we show that the resonant interaction between large-scale flows in the accreted H/He envelope of white dwarf stars and interfacial gravity waves can mix with the star's envelope with the white dwarf's surface material, leading to the enhancement of the envelope's C/O abundance to levels required by extant models for nova outbursts.

Nucleocosmochronology involves the use of the abundances of radioactive nuclear species and their radiogenic decay daughters to establish the finite age of the elements and the time scale for their formation. While there exist radioactive products of several specific nucleosynthesis mechanisms that can reveal the histories of these mechanisms, it is the long lived actinide isotopes {232}Th, {235}U, and {238}U, formed in the r-process, that currently play the major role in setting the time scale of Galactic nucleosynthesis. Age determinations in the context of Galactic chemical evolution studies are constrained by intrinsic model uncertainties. Recent studies have rather taken the alternative approach of dating individual stars. Thorium/europium dating of field halo stars and globular cluster stars yields ages on the order of 15 +/- 4 Gyr. A solid lower limit on stellar ages is available as well, for stars for which one both knows the thorium abundance and has an upper limit on the uranium abundance. For the cases of the two extremely metal deficient field halo stars CS 22892-052 and HD 115444, lower limits on the nuclear age of Galactic matter lie in the range 10-11 Gyr. Th/U dating of the star CS 31082-00 gives an age of 12.5 +/- 3 Gyr. Observations of thorium and uranium abundances in globular cluster stars should make possible nuclear determinations of the ages of clusters that can be compared with ages from conventional stellar evolution considerations.

Abundance observations indicate the presence of rapid-neutron capture (i.e., r-process) elements in old Galactic halo and globular cluster stars. These observations demonstrate that the earliest generations of stars in the Galaxy, responsible for neutron-capture synthesis and the progenitors of the halo stars, were rapidly evolving. Abundance comparisons among several halo stars show that the heaviest neutron-capture elements (including Ba and heavier) are consistent with a scaled solar system r-process abundance distribution, while the lighter such elements do not conform to the solar pattern. These comparisons suggest two r-process sites or at least two different sets of astrophysical conditions. The large star-to-star scatter observed in the neutron-capture/iron ratios at low metallicities -- which disappears with increasing [Fe/H] -- suggests an early, chemically unmixed and inhomogeneous Galaxy. The stellar abundances indicate a change from the r-process to the slow neutron capture (i.e., s-) process at higher metallicities in the Galaxy. The detection of thorium in halo and globular cluster stars offers a promising, independent age-dating technique that can put lower limits on the age of the Galaxy.

Recent observations of heavy elements produced by rapid neutron capture (r-process) in the halo have shown a striking and unexpected behavior: within a single star, the relative abundances of r-process elements heavier than Eu are the same as the same as those of solar system matter, while across stars with similar metallicity Fe/H, the r/Fe ratio varies over two orders of magnitude. In this paper we present a simple analytic model which describes a star's abundances in terms of its ``ancestry,'' i.e., the number of nucleosynthesis events (e.g., supernova explosions) which contributed to the star's composition. This model leads to a very simple analytic expression for the abundance scatter versus Fe/H, which is in good agreement with the data and with more sophisticated numerical models. We investigate two classes of scenarios for r-process nucleosynthesis, one in which r-process synthesis events occur in only \sim 4% of supernovae but iron synthesis is ubiquitous, and one in which iron nucleosynthesis occurs in only about 9% of supernovae. (the Wasserburg- Qian model). We find that the predictions in these scenarios are similar for [Fe/H] \ga -2.5, but that these models can be readily distinguished observationally by measuring the dispersion in r/Fe at [Fe/H] \la -3.

In a brief review of abundances neutron-capture elements (Z > ~30) in metal-poor halo stars, attention is called to their star-to-star scatter, the dominance of r-process synthesis at lowest metallicities, the puzzle of the lighter members of this element group, and the possibility of a better r-/s-process discriminant.

We review recent observational studies of heavy element abundances in low metallicity stars and explore some implications of these results for nucleosynthesis and early Galactic chemical evolution.

Observations of metal-poor Galactic halo stars indicate that the abundance pattern of the (heaviest) neutron-capture elements is consistent with the scaled solar system r-process abundances. Utilizing the radioactive (r-process) element thorium, age determinations have been made for several of these same stars, placing constraints on both Galactic and cosmological age estimates.

At present, there are several feasible observational probes of the first stars in the universe. Here, we examine the constraints on early stellar activity from the metallicity of the high-redshift Ly-$\alpha$ clouds, from the effects of stellar ionizing photons on reionization and the cosmic microwave background (CMB), and from the implications of gravitational microlensing results for the presence of stellar remnants in galactic halos. We also discuss whether the above signatures are consistent with each other, i.e., if they reflect the same population of stars.

We discuss some of the computational challenges encountered in simulating the evolution of clusters of galaxies. Eulerian adaptive mesh refinement (AMR) techniques can successfully address these challenges but are currently being used by only a few groups. We describe our publicly available AMR code, FLASH, which uses an object-oriented framework to manage its AMR library, physics modules, and automated verification. We outline the development of the FLASH framework to include collisionless particles, permitting it to be used for cluster simulation.

We consider the contribution of neutron stars and black holes to the dynamical mass of galactic halos. In particular, we show that if these compact objects were produced by an early generation of stars with initial metallicity less than about $10^{-4} Z_\odot$, they can contribute at most 30--40% of the Galactic halo mass without creating supersolar levels of enrichment. We show that the case for halo neutron stars and black holes cannot be rejected on metal overproduction arguments alone, due to the critical factor of the choice of progenitor metallicity in determining the yields. We show that this scenario satisfies observational constraints, similar to but no more severe than those faced by halo white dwarfs. We also discuss the recent results on halo microlensing, the presence of enriched hot gas in clusters and groups of galaxies, and other observations. If there are halo neutron stars and black holes, they will be detected by microlensing experiments in the future as longer-timescale events.

Motivated by the recent measurement of the primeval abundance of deuterium, we re-examine the nuclear inputs to big-bang nucleosynthesis (BBN). Using Monte-Carlo realization of the nuclear cross-section data to directly estimate the theoretical uncertainties for the yields of D, 3-He and 7-Li, we show that previous estimates were a factor of 2 too large. We sharpen the BBN determination of the baryon density based upon deuterium, rho_B = (3.6 +/- 0.4) * 10^{-31} g/cm^3 (Omega_B h^2 = 0.019 +/- 0.0024), which leads to a predicted 4-He abundance, Y_P = 0.246 +/- 0.0014 and a stringent limit to the equivalent number of light neutrino species: N_nu < 3.20 (all at 95% cl). The predicted 7-Li abundance, 7-Li/H = (3.5 + 1.1 - 0.9) * 10^{-10}, is higher than that observed in pop II stars, (1.7 +/- 0.3) * 10^{-10} (both, 95% cl). We identify key reactions and the energies where further work is needed.

We present first results of three-dimensional (3D-) calculations of turbulent and degenerate hydrogen-burning on top of a C+O white dwarf of one solar mass. The simulations are carried out by means of a code which solves Euler's equation for an arbitrary equation of state together with a nuclear reaction network and the energy input from nuclear reactions on a Cartesian grid covering a fraction of the white dwarf's surface and accreted atmosphere. The flow patterns we obtain are very different from those of earlier 2D simulations using the same initial conditions and the same numerical resolution. The possibility of self-enrichment of the accreted hydrogen-rich atmosphere with carbon and oxygen from the surface layers of the white dwarf during the violent phase of the burning is investigated, and it is demonstrated that self-enrichment proceeds too slowly if the accreted gas has near-solar CNO-abundances at the onset of the thermonuclear runaway. As a result, we do not find a fast nova outburst. This conclusion remains valid if the initial metallicity of the accreted gas is raised by a factor of five. Therefore we conclude that fast nova outbursts indeed require huge enrichments of C and O, as postulated from spherically symmetric models, and that the mechanism which leads to such enhancements must operate prior to the outburst.

Ground- and space-based observations of stellar heavy element abundances are providing a clearer picture of the chemical evolution of the Galaxy. A large number of (r)apid and (s)low neutron capture process elements, including the first Hubble Space Telescope observations of Pt, Os, Pb and Ge, have been identified in metal-poor, galactic halo stars. In the very low metallicity (i.e. [Fe/H] $<$ --2.0) stars, the abundance pattern of the elements from Ba through the third neutron-capture peak (Os-Pt) is consistent with a scaled solar {\it r}-process distribution. These results support previous observations that demonstrate the operation of the {\it r}-process, including the synthesis of the heaviest such elements, early in the history of the Galaxy. New ground-based observations further confirm that the {\it s}-process element Ba and the {\it r}-process element Eu were both synthesized solely by the r-process at low metallicities, and indicate the onset of the {\it s}-process occurred near [Fe/H] = --2. Over a range of metallicity (--2.90 $<$ [Fe/H] $<$ --0.86) the data indicate that there exist real star-to-star differences in the ratios of the [n-capture/Fe] abundances as well as in the actual spectra of the stars.

We present the results of two-dimensional calculations of turbulent nuclear burning of hydrogen-rich material accreted onto a white dwarf of 1 solar mass. The main aim of the present paper is to investigate the question as to whether and how the general properties of the burning are affected by numerical resolution effects. In particular, we want to see whether or not convective overshooting into the surface layers of the C+O white dwarf can lead to self-enrichment of the initially solar composition of the hydrogen-rich envelope with carbon and oxygen from the underlying white dwarf core. Our explicit hydrodynamic code is based on the PPM-method and computes the onset of the thermonuclear runaway on a Cartesian grid. In contrast to previous works we do not observe fast mixing of carbon and oxygen from the white dwarf's surface into the envelope by violent overshooting of large eddies. The main features of the flow fields in our simulations are the appearance of small persistent coherent structures of very high vorticity (and velocity) compared to the background flow. Their typical linear scales are about 10 to 20 grid zones and thus their physical size depends on the numerical resolution, i.e, their size decreases with increasing resolution. The two simulations (low and high resolution) which are presented here show only moderate differences in spatially integrated quantities such as laterally averaged temperature, energy generation rate, and chemical composition. We have not expanded both simulations equally long, but for the physical time under consideration the major difference seems to be that the highly resolved simulation is a bit less violent. In conclusion, we do find some self-enrichment, but on time-scales much longer than in previous calculations.

We present a quantitative analysis of the abundances of heavy elements in damped Ly-alpha galaxies in the sample of Lu et al. (1996). In particular, we compare the observed gas-phase abundances with those expected when the intrinsic (i.e., nucleosynthetic) pattern is the same as that in either the Sun or in Galactic halo stars and when the depletion pattern is the same as that in the warm Galactic interstellar medium, but with various values of the dust-to-metals ratio. We find that the observations are equally consistent with the solar and halo-star intrinsic patterns and that they favor some depletion, the typical dust-to-metals ratio being 40%-90% of that in the Milky Way today. However, neither intrinsic pattern matches the observations perfectly. For the solar pattern, the discrepancy is mainly with [Mn/Fe], while for the halo-star pattern, the discrepancy is with [Zn/Fe], [Ni/Fe], and possibly [Al/Fe]. Our analysis does not support the claim by Lu et al. that the damped Ly-alpha galaxies have halo-star abundance patterns and no dust depletion.

Plausible ratios of deuterium to hydrogen D/H as a function of metallicity, time, and redshift are investigated. Guided by the heavy element abundance patterns observed locally in Galactic dwarf stars and at large redshift in quasi-stellar object absorption line systems, empirical evolution of the relative abundance ratios Li/D, B/D, N/D, O/D, and F/D for QSO absorption line systems are given for the possible evolutionary patterns in D/H. Shortened abstract.

We explore possible depletion factors for Li6 in the three stars where it has been detected. To this end, we assume that 6Li/H scales as 16O/H throughout the galactic evolution, an assumption motivated by the observations of a similar scaling for 9Be. We examine the uncertainties attached to this modeling of 6Li galactic evolution; we also incorporate a recent evaluation of the primordial production of 6Li. The depletion factor D6, for 6Li, in HD84937, is constrained to be less than 4 at 95% c.l.. This depletion factor is an extreme upper limit on the depletion factor for 7Li. We conclude that the primordial abundance of 7Li cannot be higher that twice the observed value, since no current stellar model is able to reproduce the observed trends of 7Li with temperature and metallicity, and yet deplete 7Li by more than a factor 2, while satisfying D6<4. All three 6Li observations are in excellent agreement with all standard expectations: Big-Bang nucleosynthesis with 2<$\eta_{10}$<6.5, and standard isochrones of 6Li survival. We discuss possible deviations from our assumption of 6Li/H scaling with 16O/H, due to alpha-alpha creation of 6Li in various sites. We notably reject stellar flares as a significant source of 6Li in HD84937.

The prediction of standard chemical evolution models of higher abundances of He3 at the solar and present-day epochs than are observed indicates a possible problem with the yield of He3 for stars in the range of 1-3 solar masses. Because He3 is one of the nuclei produced in Big Bang Nucleosynthesis (BBN), it is noted that galactic and stellar evolution uncertainties necessarily relax constraints based on He3. We incorporate into chemical evolution models which include outflow, the new yields for He3 of Boothroyd & Malaney (1995) which predict that low mass stars are net destroyers of He3. Since these yields do not account for the high \he3/H ratio observed in some planetary nebulae, we also consider the possibility that some fraction of stars in the 1 - 3 solar mass range do not destroy their He3 in theirpost main-sequence phase. We also consider the possibility that the gas expelled by stars in these mass ranges does not mix with the ISM instantaneously thus delaying the He3 produced in these stars, according to standard yields, from reaching the ISM. In general, we find that the Galactic D and He3 abundances can be fit regardless of whether the primordial D/H value is high (2 x 10^{-4}) or low (2.5 x 10^{-5}).

In this paper, we revisit in detail the effects of primordial magnetic fields on big bang nucleosynthesis (BBN) including a discussion of the magnetic field geometry and the anomalous magnetic moment. The presence of magnetic fields affects BBN by (1) increasing the weak reaction rates; (2) increasing the electron density due to changes to the electron phase space; and (3) by increasing the expansion rate of the universe, due both to the magnetic field energy density and to the modified electron energy density. Of the effects considered, the increase in the expansion rate due to the magnetic field energy is the most significant for the interests of BBN. The allowed magnetic field intensity at the end of nucleosynthesis (0.01 MeV) is about $2 \times 10^{9}$G and corresponds to an upper limit on the magnetic field energy density of about 28\% of the neutrino energy density ($\rho_B \le 0.28 \rho_\nu$).

Deuterium plays a crucial role in testing big-bang nucleosynthesis. Its chemical evolution, while simple (it is burned to $^3$He), is intertwined with the more complicated evolution of $^3$He. Gloeckler \& Geiss' new measurement of the $^3$He abundance and the HST measurement of D, both in the local ISM today, can be compared to the pre-solar nebula abundances of D and $^3$He. Within the uncertainties, the sum of D + $^3$He relative to hydrogen is unchanged. This provides some validation of the cosmological utility of D + $^3$He, first suggested by Yang et al (1984), and further, indicates that over the past 4.5 Gyr there has been at most modest stellar production of $^3$He, in contradiction with stellar modeling, or modest stellar destruction of $^3$He, in contradiction with some ``solar spoons.'' While the earlier Galactic evolution of D + $^3$He cannot be constrained directly, it is expected to be dominated by massive stars, which deplete their $^3$He and produce metals. Based on the Galactic metallicity and the constancy of D + $^3$He over the past 4.5 Gyr, we derive a more empirically based lower bound to the cosmological baryon density; while not dramatically different from the original bound of Yang et al (1984) based on D + $^3$He, it alleviates some of the cosmic tension between the big-bang $^4$He abundance and those of D and $^3$He.

Determinations of the \he3 concentrations in Galactic matter serve to impose interesting and important constraints both on cosmological models and on models of Galactic chemical evolution. At present, observations of \he3 in the solar system and in the interstellar medium today suggest that the \he3 abundance has not increased significantly over the history of the Galaxy, while theoretical models of Galactic chemical evolution (utilizing current nucleosynthesis yields from stellar evolution and supernova models) predict a rather substantial increase in \he3. We consider the possibility that the solar \he3 abundance may have been affected by stellar processing in the solar neighborhood prior to the formation of the solar system. Such a discrepancy between solar abundances and average galactic abundances by as much as a factor of two, may be evidenced by several isotopic anomalies. Local destruction of \he3 by a similar amount could serve to help to reconcile the expected increase in the \he3 abundance predicted by models of galactic chemical evolution. We find however, that the production of heavier elements, such as oxygen, places a strong constraint on the degree of \he3 destruction. We also explore the implications of both alternative models of Galactic chemical evolution and the stellar yields for \he3 in low mass stars, which can explain the history of the \he3 concentration in the Galaxy.

Abundance histories for QSO absorption line systems as a function of redshift are presented for all metals lighter than gallium. Coupling various conventional cosmological models with a simple, first--order model for the chemical evolution of the QSO absorption line systems allows transformation of the observed abundance trends in the Galaxy to abundance histories for the gas in QSO absorption line systems. Comparison of the transformed abundance patterns with the zinc to hydrogen [Zn/H] and silicon to hydrogen [Si/H] ratios observed in damped Lyman-$\alpha$ systems finds agreement over more than two orders of magnitude in abundance with a $\Lambda$=0, \ $\Omega$=0.2 -- 1.0, \ and $\tau_{\rm {delay}}$=3 Gyr cosmology, where $\tau_{\rm {delay}}$ is the time between the start of the Big Bang and the beginning of galaxy formation. Alternative meanings for the required time delay are explored, along with extensions to the simple first--order model. Predictions are made for the abundance histories of elements not yet observed in QSO absorption line systems.

We consider the galactic evolutionary history of \he3 in models which deplete deuterium by as much as a factor of 2 to $\sim$ 15 from its primordial value to its present day observed value in the ISM. We show that when \he3 production in low mass stars (1 -- 3 $M_\odot$) is included over the history of the galaxy, \he3 is greatly over-produced and exceeds the inferred solar values and the abundances determined in galactic \hii regions. Furthermore, the ISM abundances show a disturbing dispersion which is difficult to understand from the point of view of standard chemical evolution models. In principle, resolution of the problem may lie in either 1) the calculated \he3 production in low mass stars; 2) the observations of the \he3 abundance; or 3) an observational bias towards regions of depleted \he3. Since \he3 observations in planetary nebula support the calculated \he3 production in low mass stars, option (1) is unlikely. We will argue for option (3) since the \he3 interstellar observations are indeed made in regions dominated by massive stars in which \he3 is destroyed. In conclusion, we note that the problem with \he3 seems to be galactic and not cosmological.

The ages of globular clusters inferred from observations depends sensitively on assumptions like the initial helium abundance and the mass loss rate. A high helium abundance (e.g., $Y\approx$0.28) or a mass loss rate of $\sim$10$^{-11}M_\odot$ yr$^{-1}$ near the main sequence turn-off region lowers the current age estimate from 14 Gyr to about 10--12 Gyr, significantly relaxing the constraints on the Hubble constant, allowing values as high as 60km/sec/Mpc for a universe with the critical density and 90km/sec/Mpc for a baryon-only universe. Possible mechanisms for the helium enhancement in globular clusters are discussed, as are arguments for an instability strip induced mass loss near the turn-off. Ages lower than 10 Gyr are not possible even with the operation of both of these mechanisms unless the initial helium abundance in globular clusters is $>0.30$, which would conflict with indirect measurements of helium abundances in globular clusters.

The effects of magnetic fields on Big Bang Nucleosynthesis(BBN) have been calculated, and the impact on the abundances of the light elements have been investigated numerically. An upper limit on the strength of primordial magnetic fields compatible with observations of light element abundances has been thus obtained. In the framework of standard BBN theory, the maximum strength of the primordial magnetic fields, on scales greater than $10^4$cm but smaller than the event horizon at the BBN epoch ($\sim 1$ min., $\sim 2\times 10^{12}$cm), is $\le 10^{11}$ Gauss. This limit is shown to allow magnetic fields at the time of recombination no stronger than $\sim 0.1$ gauss on scales $\ge 10^{11}{\rm cm}$. Our results also strongly indicate that, at the BBN epoch, and for field strengths $B\le10^{13}$gauss, the effects of magnetic fields on the primordial abundances of light elements are dominated by effects from reaction rates in the presense of primeval magnetic fields rather than by magnetic density effects on the expansion rate.

In this paper, we have derived the the effects of strong magnetic fields ${\bf \vec B}$ on nucleon and particle reaction rates of astrophysical significance. We have explored the sensitivity to the presence of arbitrary degeneracy and polarization. The possible astrophysical applications of our results are discussed.