The formation and evolution of stars depends on various physical aspects of stellar matter, including the equation of state (EOS) and transport properties. Although often dismissed as `ideal gas-like' and therefore simple, states occurring in stellar matter are dense plasmas, and the EOS has not been established precisely. EOS constructed using multi-physics approaches found necessary for laboratory studies of warm dense matter give significant variations in stellar regimes, and vary from the EOS commonly used in simulations of the formation and evolution of stars. We have investigated the sensitivity of such simulations to variations in the EOS, for sun-like and low-mass stars, and found a strong sensitivity of the lifetime of the Sun and of the lower luminosity limit for red dwarfs. We also found a significant sensitivity in the lower mass limit for red dwarfs. Simulations of this type are also used for other purposes in astrophysics, including the interpretation of absolute magnitude as mass, the conversion of inferred mass distribution to the initial mass function using predicted lifetimes, simulations of star formation from nebulae, simulations of galactic evolution, and the baryon census used to bound the exotic contribution to dark matter. Although many of the sensitivities of stellar physics to the EOS are large, some of the inferred astrophysical quantities are also constrained by independent measurements, although the constraints may be indirect and non-trivial. However, it may be possible to use such measurements to constrain the EOS more than presently possible by established plasma theory.

We have constructed a series of non-rotating quasi-hydrostatic evolutionary models for the M2 Iab supergiant Betelgeuse ($\alpha~Orionis$). Our models are constrained by multiple observed values for the temperature, luminosity, surface composition and mass loss for this star, along with the parallax distance and high resolution imagery that determines its radius. We have then applied our best-fit models to analyze the observed variations in surface luminosity and the size of detected surface bright spots as the result of up-flowing convective material from regions of high temperature in the surface convective zone. We also attempt to explain the intermittently observed periodic variability in a simple radial linear adiabatic pulsation model. Based upon the best fit to all observed data, we suggest a best progenitor mass estimate of $ 20 ^{+5}_{-3} M_\odot$ and a current age from the start of the zero-age main sequence of $8.0 - 8.5$ Myr based upon the observed ejected mass while on the giant branch.

We use the 3D stellar structure code DJEHUTY to model the ingestion of protons into the intershell convection zone of a 1 solar mass asymptotic giant branch star of metallicity Z=10^-4. We have run two simulations: a low resolution one of around 300,000 zones, and a high resolution one consisting of 2,000,000 zones. Both simulations have been evolved for about 4 hours of stellar time. We observe the existence of fast, downward flowing plumes that are able to transport hydrogen into close proximity to the helium burning shell before burning takes place. The intershell in the 3D model is richer in protons than the 1D model by several orders of magnitude and so we obtain substantially higher hydrogen-burning luminosities -- over 10^8 solar luminosities in the high resolution simulation -- than are found in the 1D model. Convective velocities in these simulations are over 10 times greater than the predictions of mixing length theory, though the 3D simulations have greater energy generation due to the enhanced hydrogen burning. We find no evidence of the convective zone splitting into two, though this could be as a result of insufficient spatial resolution or because the models have not been evolved for long enough. We suggest that the 1D mixing length theory and particularly the use of a diffusion algorithm for mixing do not give an accurate picture of these events. An advective mixing scheme may give a better representation of the transport processes seen in the 3D models.

Three-dimensional stellar modeling has enabled us to identify a deep-mixing mechanism that must operate in all low mass giants. This mixing process is not optional, and is driven by a molecular weight inversion created by the 3He(3He,2p)4He reaction. In this paper we characterize the behavior of this mixing, and study its impact on the envelope abundances. It not only eliminates the problem of 3He overproduction, reconciling stellar and big bang nucleosynthesis with observations, but solves the discrepancy between observed and calculated CNO isotope ratios in low mass giants, a problem of more than 3 decades' standing. This mixing mechanism, which we call `$\delta\mu$-mixing', operates rapidly (relative to the nuclear timescale of overall evolution, ~ 10^8 yrs) once the hydrogen burning shell approaches the material homogenized by the surface convection zone. In agreement with observations, Pop I stars between 0.8 and 2.0$\Msun$ develop 12C/13C ratios of 14.5 +/- 1.5, while Pop II stars process the carbon to ratios of 4.0 +/- 0.5. In stars less than 1.25$\Msun$, this mechanism also destroys 90% to 95% of the 3He produced on the main sequence.

We continue our study of the core helium flash using the three dimensional hydrodynamics code Djehuty. Continuing from earlier calculations, we now take relaxed 3D configurations and add various amounts of rotation. We find that rotation periods consistent with those observed in white dwarfs produce negligible changes in the structure and evolution of the core flash, at least for the very small timescales we have yet been able to investigate. There is no sign of any extra mixing due to the rotation. There is some inconclusive evidence for a slight change in the luminosity, at the 1% level.

Low-mass stars, ~1-2 solar masses, near the Main Sequence are efficient at producing 3He, which they mix into the convective envelope on the giant branch and distribute into the Galaxy by way of envelope loss. This process is so efficient that it is difficult to reconcile the observed cosmic abundance of 3He with the predictions of Big Bang nucleosynthesis. In this paper we find, by modeling a red giant with a fully three-dimensional hydrodynamic code and a full nucleosynthetic network, that mixing arises in the supposedly stable and radiative zone between the hydrogen-burning shell and the base of the convective envelope. This mixing is due to Rayleigh-Taylor instability within a zone just above the hydrogen-burning shell. In this zone the burning of the 3He left behind by the retreating convective envelope is predominantly by the reaction 3He + 3He -> 4He + 2p, a reaction which, untypically for stellar nuclear reactions, {\it lowers} the mean molecular weight, leading to a local minimum. This local minimum leads to Rayleigh-Taylor instability, and turbulent motion is generated which will continue ultimately up into the normal convective envelope. Consequently material from the envelope is dragged down sufficiently close to the burning shell that the 3He in it is progressively destroyed. Thus we are able to remove the threat that 3He production in low-mass stars poses to the Big Bang nucleosynthesis of 3He. Some slow mixing mechanism has long been suspected, that connects the convective envelope of a red giant to the burning shell. It appears to be necessary to account for progressive changes in the 12C/13C and 14N/12C ratios on the First Giant Branch. We suggest that these phenomena are also due to the Rayleigh-Taylor-unstable character of the 3He-burning region.

Low-mass stars, ~1-2 solar masses, near the Main Sequence are efficient at producing He-3, which they mix into the convective envelope on the giant branch and should distribute into the Galaxy by way of envelope loss. This process is so efficient that it is difficult to reconcile the low observed cosmic abundance of He-3 with the predictions of both stellar and Big Bang nucleosynthesis. In this paper we find, by modeling a red giant with a fully three-dimensional hydrodynamic code and a full nucleosynthetic network, that mixing arises in the supposedly stable and radiative zone between the hydrogen-burning shell and the base of the convective envelope. This mixing is due to Rayleigh-Taylor instability within a zone just above the hydrogen-burning shell, where a nuclear reaction lowers the mean molecular weight slightly. Thus we are able to remove the threat that He-3 production in low-mass stars poses to the Big Bang nucleosynthesis of He-3.

We model the core helium flash in a low-mass red giant using Djehuty, a fully three-dimensional (3D) code. The 3D structures were generated from converged models obtained during the 1D evolutionary calculation of a 1$\Msun$ star. Independently of which starting point we adopted, we found that after some transient relaxation the 3D model settled down with a briskly convecting He-burning shell that was not very different from what the 1D model predicted.

Current practice in stellar evolution is to employ one-dimensional calculations that quantitatively apply only to a minority of the observed stars (single non-rotating stars, or well detached binaries). Even in these systems, astrophysicists are dependent on approximations to handle complex three-dimensional processes like convection. Understanding the structure of binary stars, like those that lead to the Type Ia supernovae used to measure the expansion of the universe, are grossly non-spherical and await a 3D treatment. To approach very large problems like multi-dimensional modeling of stars, the Lawrence Livermore National Laboratory has invested in massively parallel computers and invested even more in developing the algorithms to utilize them on complex physics problems. We have leveraged skills from across the lab to develop a 3D stellar evolution code, Djehuty (after the Egyptian god for writing and calculation) that operates efficiently on platforms with thousands of nodes, with the best available physical data (opacities, EOS, etc.). Djehuty has incorporated all basic physics for modeling stars including an accurate equation of state, radiation transport by diffusion, thermonuclear reaction rates, and hydrodynamics, and we have begun testing it in a number of applications.

The DJEHUTY project is an intensive effort at the Lawrence Livermore National Laboratory (LLNL) to produce a general purpose 3-D stellar structure and evolution code to study dynamic processes in whole stars.

New stellar models which track the production and destruction of $^3$He (and D) have been evolved for a range of stellar masses $(0.65\leq M/M_{\odot}\leq 100)$, metallicities $(0.01 \leq Z/Z_{\odot} \leq 1)$ and initial (main sequence) $^3$He mass fractions $(10^{-5} \leq X_{3,MS} \leq 10^{-3})$. Armed with the $^3$He yields from these stellar models we have followed the evolution of D and $^3$He using a variety of chemical evolution models with and without infall of primordial or processed material. Production of new $^3$He by the lower mass stars overwhelms any reasonable primordial contributions and leads to predicted abundances in the presolar nebula and/or the present interstellar medium in excess of the observationally inferred values. This result, which obtains even for zero primordial D and $^3$He, and was anticipated by Rood, Steigman \& Tinsley (1976), is insensitive to the choice of chemical evolution model; it is driven by the large $^3$He yields from low mass stars. In an attempt to ameliorate this problem we have considered a number of non-standard models in which the yields from low mass stars have been modified. Although several of these non-standard models may be consistent with the $^3$He data, they may be inconsistent with observations of $^{12}$C/$^{13}$C, $^{18}$O and, most seriously, the super-$^3$He rich planetary nebulae (Rood, Bania \& Wilson 1992). Even using the most extreme of these non-standard models (Hogan 1995), we obtain a generous upper bound to pre-galactic $^3$He: X$_{3P} \leq 3.2 \times10^{-5}$ which, nonetheless, leads to a stringent lower bound to the universal density of nucleons.

Without assuming any solar models and neutrino flavor conversions, $^8$B neutrinos seen by the Kamiokande experiment should contribute 2.6$\pm 0.45$ SNU to the chlorine experiment. When this rate is compared with the total event rate of 2.3$\pm 0.2$ SNU observed by the Homestake experiment which should inlude a 0.2 SNU contribution from uncertainty-free pep neutrinos, there may still be a possible evidence that $^7$Be neutrinos are more severely suppressed than the $^8$B neutrinos with respect to the predictions of standard solar models, which cannot be explained by any known astrophysics solution. Given a Kamiokande event rate of larger than 36$\%$ (2$\sigma$) of the prediciton of Bahcall and Pinsonneault's standard solar model, variations of standard solar models yield minimal rates of 3.6 SNU for the Homestake experiment and 114 SNU for GALLEX and SAGE. Therefore, variations of standard solar models as solutions to the solar neutrino problem are so far inconsistent with the Homestake experiment and only marginally allowed by the gallium experiments. If the gallium experiments confirm a flux significantly below 114 SNU, it would imply unconventional physics.

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.