Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

Convective boundary mixing is one of the major uncertainties in stellar evolution. In order to study its dependence on boundary properties and turbulence strength in a controlled way, we computed a series of 3D hydrodynamical simulations of stellar convection during carbon burning with a varying boosting factor of the driving luminosity. Our 3D implicit large eddy simulations were computed with the PROMPI code. We performed a mean field analysis of the simulations within the Reynolds-averaged Navier-Stokes framework. Both the vertical RMS velocity within the convective region and the bulk Richardson number of the boundaries are found to scale with the driving luminosity as expected from theory. The positions of the convective boundaries were estimated through the composition profiles across them, and the strength of convective boundary mixing was determined by analysing the boundaries within the framework of the entrainment law. We find that the entrainment is approximately inversely proportional to the bulk Richardson number. Although the entrainment law does not encompass all the processes occurring at boundaries, our results support the use of the entrainment law to describe convective boundary mixing in 1D models, at least for the advanced phases. The next steps and challenges ahead are also discussed.

The turbulent burning of nuclei is a common phenomenon in the evolution of stars. Here we examine a challenging case: the merging of the neon and oxygen burning shells in a 23 M$_{\odot}$ star. A previously unknown quasi-steady state is established by the interplay between mixing, turbulent transport, and nuclear burning. The resulting stellar structure has two burning shells within a single convection zone. We find that the new neon burning layer covers an extended region of the convection zone, with the burning peak occurring substantially below where the Damk\"ohler number first becomes equal to unity. These characteristics differ from those predicted by 1D stellar evolution models of similar ingestion events. We develop the mean-field turbulence equations that govern compositional evolution, and use them to interpret our data set. An important byproduct is a means to quantify sub-grid-scale effects intrinsic to the numerical hydrodynamic scheme. For implicit large eddy simulations, the analysis method is particularly powerful because it can reveal where and how simulated flows are modified by resolution, and provide straightforward physical interpretations of the effects of dissipation or induced transport. Focusing on the mean-field composition variance equations for our analysis, we recover a Kolmogorov rate of turbulent dissipation without it being imposed, in agreement with previous results which used the turbulent kinetic energy equation.

Massive stars have a strong impact on their surroundings, in particular when they produce a core-collapse supernova at the end of their evolution. In these proceedings, we review the general evolution of massive stars and their properties at collapse as well as the transition between massive and intermediate-mass stars. We also summarise the effects of metallicity and rotation. We then discuss some of the major uncertainties in the modelling of massive stars, with a particular emphasis on the treatment of convection in 1D stellar evolution codes. Finally, we present new 3D hydrodynamic simulations of convection in carbon burning and list key points to take from 3D hydrodynamic studies for the development of new prescriptions for convective boundary mixing in 1D stellar evolution codes.

We present the first detailed three-dimensional hydrodynamic implicit large eddy simulations of turbulent convection for carbon burning. The simulations start with an initial radial profile mapped from a carbon burning shell within a 15 solar mass stellar evolution model. We considered 4 resolutions from 128^3 to 1024^3 zones. These simulations confirm that convective boundary mixing (CBM) occurs via turbulent entrainment as in the case of oxygen burning. The expansion of the boundary into the surrounding stable region and the entrainment rate are smaller at the bottom boundary because it is stiffer than the upper boundary. The results of this and similar studies call for improved CBM prescriptions in 1D stellar evolution models.

Alastair Graham Walker Cameron was an astrophysicist and planetary scientist of broad interests and exceptional originality. A founder of the field of nuclear astrophysics, he developed the theoretical understanding of the chemical elements’ origins and made pioneering connections between the abundances of elements in meteorites to advance the theory that the Moon originated from a giant impact with the young Earth by an object at least the size of Mars. Cameron was an early and persistent exploiter of computer technology in the theoretical study of complex astronomical systems—including nuclear reactions in supernovae, the structure of neutron stars, and planetary collisions.

We present the first detailed three-dimensional (3D) hydrodynamic implicit large eddy simulations of turbulent convection of carbon burning in massive stars. Simulations begin with radial profiles mapped from a carbon burning shell within a 15$\,\textrm{M}_\odot$ one-dimensional stellar evolution model. We consider models with $128^3$, $256^3$, $512^3$ and $1024^3$ zones. The turbulent flow properties of these carbon burning simulations are very similar to the oxygen burning case. We performed a mean field analysis of the kinetic energy budgets within the Reynolds-averaged Navier-Stokes framework. For the upper convective boundary region, we find that the numerical dissipation is insensitive to resolution for linear mesh resolutions above 512 grid points. For the stiffer, more stratified lower boundary, our highest resolution model still shows signs of decreasing sub-grid dissipation suggesting it is not yet numerically converged. We find that the widths of the upper and lower boundaries are roughly 30% and 10% of the local pressure scale heights, respectively. The shape of the boundaries is significantly different from those used in stellar evolution models. As in past oxygen-shell burning simulations, we observe entrainment at both boundaries in our carbon-shell burning simulations. In the large P\'eclet number regime found in the advanced phases, the entrainment rate is roughly inversely proportional to the bulk Richardson number, Ri$_{\rm B}$ ($\propto $Ri${\rm_B}^{-\alpha}$, $0.5\lesssim \alpha \lesssim 1.0$). We thus suggest the use of Ri$_{\rm B}$ as a means to take into account the results of 3D hydrodynamics simulations in new 1D prescriptions of convective boundary mixing.

Stellar evolution models of massive stars are important for many areas of astrophysics, for example nucleosynthesis yields, supernova progenitor models and understanding physics under extreme conditions. Turbulence occurs in stars primarily due to nuclear burning at different mass coordinates within the star. The understanding and correct treatment of turbulence and turbulent mixing at convective boundaries in stellar models has been studied for decades but still lacks a definitive solution. This paper presents initial results of a study on convective boundary mixing (CBM) in massive stars. The 'stiffness' of a convective boundary can be quantified using the bulk Richardson number ($\textrm{Ri}_B$), the ratio of the potential energy for restoration of the boundary to the kinetic energy of turbulent eddies. A 'stiff' boundary ($\textrm{Ri}_B \sim 10^4$) will suppress CBM, whereas in the opposite case a 'soft' boundary ($\textrm{Ri}_B \sim 10$) will be more susceptible to CBM. One of the key results obtained so far is that lower convective boundaries (closer to the centre) of nuclear burning shells are 'stiffer' than the corresponding upper boundaries, implying limited CBM at lower shell boundaries. This is in agreement with 3D hydrodynamic simulations carried out by Meakin and Arnett [The Astrophysical Journal 667:448-475, 2007]. This result also has implications for new CBM prescriptions in massive stars as well as for nuclear burning flame front propagation in Super-Asymptotic Giant Branch stars and also the onset of novae.

Several parametrizations for overshooting in 1D stellar evolution calculations coexist in the literature. These parametrizations are used somewhat arbitrarily in stellar evolution codes, based on what works best for a given problem, or even for historical reasons related to the development of each code. We bring attention to the fact that these different parametrizations correspond to different physical regimes of overshooting, depending whether the effects of radiation are dominant, marginal, or negligible. Our analysis is based on previously published theoretical results, as well as multidimensional hydrodynamical simulations of stellar convection where the interaction between the convective region and a stably-stratified region is observed. Although the underlying hydrodynamical processes are the same, the outcome of the overshooting process is profoundly affected by radiative effects. Using a simple picture of the scales involved in the overshooting process, we show how three regimes are obtained, depending on the importance of radiative effects. These three regimes correspond to the different behaviors observed in hydrodynamical simulations so far, and to the three types of parametrizations used in 1D codes. We suggest that the existing parametrizations for overshooting should coexist in 1D stellar evolution codes, and should be applied consistently at convective boundaries depending on the local physical conditions.

In this contribution we present initial results of a study on convective boundary mixing (CBM) in massive stellar models using the GENEVA stellar evolution code. Before undertaking costly 3D hydrodynamic simulations, it is important to study the general properties of convective boundaries, such as the: composition jump; pressure gradient; and `stiffness'. Models for a 15Mo star were computed. We found that for convective shells above the core, the lower (in radius or mass) boundaries are `stiffer' according to the bulk Richardson number than the relative upper (Schwarzschild) boundaries. Thus, we expect reduced CBM at the lower boundaries in comparison to the upper. This has implications on flame front propagation and the onset of novae.

Both observations and direct numerical simulations are discordant with predictions of conventional stellar evolution codes for the latest stages of a massive star's life prior to core collapse. We suggest that the problem lies in the treatment of turbulent convection in these codes, which ignores finite amplitude fluctuations in velocity and temperature, and their nonlinear interaction with nuclear burning. The hydrodynamic instabilities that may arise prompt us to discuss a number of far-reaching implications for the fates of massive stars. In particular, we explore connections to enhanced presupernova mass loss, unsteady nuclear burning and consequent eruptions, swelling of the stellar radius that may trigger violent interactions with a companion star, and potential modifications to the core structure that could dramatically impact calculations of the core-collapse mechanism itself. These modifications may be of fundamental importance to the interpretation of measured isotopic anomalies in meteorites, changing the predictions of both mixing and detailed nucleosynthesis, and of young core-collapse supernova remnants. They may also make possible the development of an early warning system for the detection of impending core collapse.

We present a statistical analysis of turbulent convection in stars within our Reynolds-Averaged Navier Stokes (RANS) framework in spherical geometry which we derived from first principles. The primary results reported in this document include: (1) an extensive set of mean-field equations for compressible, multi-species hydrodynamics, and (2) corresponding mean-field data computed from various simulation models. Some supplementary scale analysis data is also presented. The simulation data which is presented includes: (1) shell convection during oxygen burning in a 23 solar mass supernova progenitor, (2) envelope convection in a 5 solar mass red giant, (3) shell convection during the helium flash, and (4) a hydrogen injection flash in a 1.25 solar mass star. These simulations have been partially described previously in Meakin [2006], Meakin and Arnett [2007a,b, 2010], Arnett et al. [2009, 2010], Viallet et al. [2011, 2013a,b] and Mocak et al. [2009, 2011]. New data is also included in this document with several new domain and resolution configurations as well as some variations in the physical model such as convection zone depth and driving source term. The long term goal of this work is to aid in the development of more sophisticated models for treating hydrodynamic phenomena (e.g., turbulent convection) in the field of stellar evolution by providing a direct link between 3D simulation data and the mean fields which are modeled by 1D stellar evolution codes. As such, this data can be used to test previously proposed turbulence models found in the literature and sometimes used in stellar modeling. This data can also serve to test basic physical principles for model building and inspire new prescriptions for use in 1D evolution codes.

The 2012 explosion of SN2009ip raises questions about our understanding of the late stages of massive star evolution. Here we present a comprehensive study of SN2009ip during its remarkable re-brightening(s). High-cadence photometric and spectroscopic observations from the GeV to the radio band obtained from a variety of ground-based and space facilities (including the VLA, Swift, Fermi, HST and XMM) constrain SN2009ip to be a low energy (E~ 10^50 erg for an ejecta mass ~ 0.5 Msun) and likely asymmetric explosion in a complex medium shaped by multiple eruptions of the restless progenitor star. Most of the energy is radiated as a result of the shock breaking out through a dense shell of material located at 5x10^14 cm with M~0.1 Msun, ejected by the precursor outburst ~40 days before the major explosion. We interpret the NIR excess of emission as signature of dust vaporization of material located further out (R>4x 10^15 cm), the origin of which has to be connected with documented mass loss episodes in the previous years. Our modeling predicts bright neutrino emission associated with the shock break-out if the cosmic ray energy is comparable to the radiated energy. We connect this phenomenology with the explosive ejection of the outer layers of the massive progenitor star, that later interacted with material deposited in the surroundings by previous eruptions. Future observations will reveal if the luminous blue variable (LBV) progenitor star survived. Irrespective of whether the explosion was terminal, SN2009ip brought to light the existence of new channels for sustained episodic mass-loss, the physical origin of which has yet to be identified.

The Procyon AB binary system (orbital period 40.838 years, a newly-refined determination), is near and bright enough that the component radii, effective temperatures, and luminosities are very well determined, although more than one possible solution to the masses has limited the claimed accuracy. Preliminary mass determinations for each component are available from HST imaging, supported by ground-based astrometry and an excellent Hipparcos parallax; we use these for our preferred solution for the binary system. Other values for the masses are also considered. We have employed the TYCHO stellar evolution code to match the radius and luminosity of the F5 IV-V primary star to determine the system's most likely age as 1.87 +/- 0.13 Gyr. Since prior studies of Procyon A found its abundance indistinguishable from solar, the solar composition of Asplund, Grevesse & Sauval (Z=0.014) is assumed for the HR Diagram flitting. An unsuccessful attempt to fit using the older solar abundance scale of Grevesse & Sauval (Z=0.019) is also reported. For Procyon B, eleven new sequences for the cooling of non-DA white dwarfs have been calculated, to investigate the dependence of the cooling age on (1) the mass, (2) the core composition, (3) the helium layer mass, and (4) heavy-element opacities in the helium envelope. Our calculations indicate a cooling age of 1.19+/-0.11 Gyr, which implies that the progenitor mass of Procyon B was 2.59(+0.44,-0.26) Msun. In a plot of initial vs final mass of white dwarfs in astrometric binaries or star clusters (all with age determinations), the Procyon B final mass lies several sigma below a straight line fit.

We present 3D implicit large eddy simulations (ILES) of the turbulent convection in the envelope of a 5 Msun red giant star and in the oxygen-burning shell of a 23 Msun supernova progenitor. The numerical models are analyzed in the framework of 1D Reynolds-Averaged Navier-Stokes (RANS) equations. The effects of pressure fluctuations are more important in the red giant model, owing to larger stratification of the convective zone. We show how this impacts different terms in the mean-field equations. We clarify the driving sources of kinetic energy, and show that the rate of turbulent dissipation is comparable to the convective luminosity. Although our flows have low Mach number and are nearly adiabatic, our analysis is general and can be applied to photospheric convection as well. The robustness of our analysis of turbulent convection is supported by the insensitivity of the mean-field balances to linear mesh resolution. We find robust results for the turbulent convection zone and the stable layers in the oxygen-burning shell model, and robust results everywhere in the red giant model, but the mean fields are not well converged in the narrow boundary regions (which contain steep gradients) in the oxygen-burning shell model. This last result illustrates the importance of unresolved physics at the convective boundary, which governs the mixing there.

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.

Issues concerning the structure and evolution of core collapse progenitor stars are discussed with an emphasis on interior evolution. We describe a program designed to investigate the transport and mixing processes associated with stellar turbulence, arguably the greatest source of uncertainty in progenitor structure, besides mass loss, at the time of core collapse. An effort to use precision observations of stellar parameters to constrain theoretical modeling is also described.

Supernovae (SNe) are thought to arise from two different physical processes. The cores of massive, short-lived stars undergo gravitational core collapse and typically eject a few solar masses during their explosion. These are thought to appear as as type Ib/c and II SNe, and are associated with young stellar populations. A type Ia SN is thought to arise from the thermonuclear detonation of a white dwarf star composed mainly of carbon and oxygen, whose mass approaches the Chandrasekhar limit. Such SNe are observed in both young and old stellar environments. Here we report our discovery of the faint type Ib SN 2005E in the halo of the nearby isolated galaxy, NGC 1032. The lack of any trace of recent star formation near the SN location (Fig. 1), and the very low derived ejected mass (~0.3 M_sun), argue strongly against a core-collapse origin for this event. Spectroscopic observations and the derived nucleosynthetic output show that the SN ejecta have high velocities and are dominated by helium-burning products, indicating that SN 2005E was neither a subluminous nor a regular SN Ia (Fig. 2). We have therefore found a new type of stellar explosion, arising from a low-mass, old stellar system, likely involving a binary with a primary white dwarf and a helium-rich secondary. The SN ejecta contain more calcium than observed in any known type of SN and likely additional large amounts of radioactive 44Ti. Such SNe may thus help resolve fundamental physical puzzles, extending from the composition of the primitive solar system and that of the oldest stars, to the Galactic production of positrons.

We investigate simulated turbulent flow within thermally driven stellar convection zones. Different driving sources are studied, including cooling at the top of the convectively unstable region, as occurs in surface convection zones; and heating at the base by nuclear burning. The transport of enthalpy and kinetic energy, and the distribution of turbulent kinetic energy dissipation are studied. We emphasize the importance of global constraints on shaping the quasi-steady flow characteristics, and present an analysis of turbulent convection which is posed as a boundary value problem that can be easily incorporated into standard stellar evolution codes for deep, efficient convection. Direct comparison is made between the theoretical analysis and the simulated flow and very good agreement is found. Some common assumptions traditionally used to treat quasi-steady turbulent flow in stellar models are briefly discussed. The importance and proper treatment of convective boundaries are indicated.

A survey that can cover the sky in optical bands over wide fields to faint magnitudes with a fast cadence will enable many of the exciting science opportunities of the next decade. The Large Synoptic Survey Telescope (LSST) will have an effective aperture of 6.7 meters and an imaging camera with field of view of 9.6 deg^2, and will be devoted to a ten-year imaging survey over 20,000 deg^2 south of +15 deg. Each pointing will be imaged 2000 times with fifteen second exposures in six broad bands from 0.35 to 1.1 microns, to a total point-source depth of r~27.5. The LSST Science Book describes the basic parameters of the LSST hardware, software, and observing plans. The book discusses educational and outreach opportunities, then goes on to describe a broad range of science that LSST will revolutionize: mapping the inner and outer Solar System, stellar populations in the Milky Way and nearby galaxies, the structure of the Milky Way disk and halo and other objects in the Local Volume, transient and variable objects both at low and high redshift, and the properties of normal and active galaxies at low and high redshift. It then turns to far-field cosmological topics, exploring properties of supernovae to z~1, strong and weak lensing, the large-scale distribution of galaxies and baryon oscillations, and how these different probes may be combined to constrain cosmological models and the physics of dark energy.

We present new reddening measurements towards the young supernova remnant Cassiopeia A, using two techniques not previously applied to this object. Our observations of the near-infrared [Fe II] 1.257 micron and 1.644 micron lines show the extinction to be highly variable across the remnant, increasing towards the west and the south, consistent with previous radio and X-ray observations. While the absolute value of AV as determined by the [Fe II] lines is uncertain due to conflicting calculations and observations of their intrinsic flux ratio, parts of the remnant without previous optical measurements show comparatively higher reddening. We find AV = 6.2 +/- 0.6 from the broadband shape of the infrared synchrotron emission of a knot within 13" of the expansion center. Given this reddening, the apparent faintness of the nascent supernova, and iron mass constraints from X-ray observations, we estimate an ejected mass of 56Ni of 0.058 - 0.16 MSun. Taken with gamma-ray observations of the 44Ti decay chain, this nickel mass is broadly consistent with the solar 44Ca/56Fe ratio.

We examine observed heavy element abundances in the Cassiopeia A supernova remnant as a constraint on the nature of the Cas A supernova. We compare bulk abundances from 1D and 3D explosion models and spatial distribution of elements in 3D models with those derived from X-ray observations. We also examine the cospatial production of 26Al with other species. We find that the most reliable indicator of the presence of 26Al in unmixed ejecta is a very low S/Si ratio (~0.05). Production of N in O/S/Si-rich regions is also indicative. The biologically important element P is produced at its highest abundance in the same regions. Proxies should be detectable in supernova ejecta with high spatial resolution multiwavelength observations.

We analyze stellar convection with the aid of 3D hydrodynamic simulations, introducing the turbulent cascade into our theoretical analysis. We devise closures of the Reynolds-decomposed mean field equations by simple physical modeling of the simulations (we relate temperature and density fluctuations via coefficients); the procedure (CABS, Convection Algorithms Based on Simulations) is terrestrially testable and is amenable to systematic improvement. We develop a turbulent kinetic energy equation which contains both nonlocal and time dependent terms, and is appropriate if the convective transit time is shorter than the evolutionary time scale. The interpretation of mixing-length theory (MLT) as generally used in astrophysics is incorrect; MLT forces the mixing length to be an imposed constant. Direct tests show that the damping associated with the flow is that suggested by Kolmogorov. The eddy size is approximately the depth of the convection zone, and this dissipation length corresponds to the "mixing length". New terms involving local heating by turbulent dissipation should appear in the stellar evolution equations. The enthalpy flux ("convective luminosity") is directly connected to the buoyant acceleration, and hence the velocity scale. MLT tends to systematically underestimate this velocity scale. Quantitative comparison with a variety of 3D simulations reveals a previously recognized consistency. Examples of application to stellar evolution will be presented in subsequent papers in this series.

This paper summarizes our theoretical understanding of supernova events in a "back of the envelope" way. It is intended to aid in the recognition and understanding of those events which are not "standard", and which may provide the most insight.

(Abridged) We describe the results of three-dimensional (3D) numerical simulations designed to study turbulent convection in the stellar interiors, and compare them to stellar mixing-length theory (MLT). Simulations in 2D are significantly different from 3D, both in terms of flow morphology and velocity amplitude. Convective mixing regions are better predicted using a [dynamic boundary condition] based on the bulk Richardson number than by purely local, static criteria like Schwarzschild or Ledoux. MLT gives a good description of the velocity scale and temperature gradient for a mixing length of $\sim 1.1 H_p$ for shell convection, however there are other important effects that it does not capture near boundaries. Convective "overshooting" is best described as an elastic response by the convective boundary, rather than ballistic penetration of the stable layers by turbulent eddies. We find that the rate at which material entrainment proceeds at the boundaries is consistent with analogous laboratory experiments as well as simulation and observation of terrestrial atmospheric mixing. In particular, the normalized entrainment rate E=$u_E/\sigma_H$, is well described by a power law dependence on the bulk Richardson number $Ri_B = \Delta b L/\sigma_H^2$ for the conditions studied, $20\lesssim Ri_B \lesssim 420$. We find $E = A Ri_B^{-n}$, with best fit values, $\log A = 0.027 \pm 0.38$, and $n = 1.05 \pm 0.21$. We discuss the applicability of these results to stellar evolution calculations.

In this paper we compare fully compressible (Meakin & Arnett 2006a,b) and anelastic (Kuhlen, Woosley, & Glatzmaier 2003) simulations of stellar oxygen shell burning. It is found that the two models are in agreement in terms of the velocity scale (v_c ~ 1e7 cm/s) and thermodynamic fluctuation amplitudes (e.g., rho'/<rho> ~ 2e-3) in the convective flow. Large fluctuations (~11%) arise in the compressible model, localized to the convective boundaries, and are due to internal waves excited in stable layers. Fluctuations on the several percent level are also present in the compressible model due to composition inhomogeneities from ongoing entrainment events at the convective boundaries. Comparable fluctuations (with amplitudes greater than ~1%) are absent in the anelastic simulation because they are due to physics not included in that model. We derive an analytic estimate for the expected density fluctuation amplitudes at convective boundaries by assuming that the pressure fluctuations due to internal waves at the boundary, p_w', balance the ram pressure of the convective motions, rho*v_c2. The predicted amplitudes agree well with the simulation data. The good agreement between the anelastic and the compressible solution within the convection zone and the agreement between the stable layer dynamics and analytic solutions to the non-radial wave equation indicate that the compressible hydrodynamic techniques used are robust for the simulated stellar convection model, even at the low Mach n umbers found M~0.01.

We have simulated 2.5$\times10^3$ s of the late evolution of a $23 \rm M_\odot$ star with full hydrodynamic behavior. We present the first simulations of a multiple-shell burning epoch, including the concurrent evolution and interaction of an oxygen and carbon burning shell. In addition, we have evolved a 3D model of the oxygen burning shell to sufficiently long times (300 s) to begin to assess the adequacy of the 2D approximation. We summarize striking new results: (1) strong interactions occur between active carbon and oxygen burning shells, (2) hydrodynamic wave motions in nonconvective regions, generated at the convective-radiative boundaries, are energetically important in both 2D and 3D with important consequences for compositional mixing, and (3) a spectrum of mixed p- and g-modes are unambiguously identified with corresponding adiabatic waves in these computational domains. We find that 2D convective motions are exaggerated relative to 3D because of vortex instability in 3D. We discuss the implications for supernova progenitor evolution and symmetry breaking in core collapse.

We compare a suite of 3D explosion calculations and stellar models incorporating advanced physics with observational constraints on the progenitor of Cassiopeia A. We consider binary and single stars from 16 to 40 solar masses with a range of explosion energies and geometries. The parameter space allowed by observations of nitrogen rich high velocity ejecta, ejecta mass, compact remnant mass, and 44Ti and 56Ni abundances individually and as an ensemble is considered. A progenitor of 15-25 solar masses which loses its hydrogen envelope to a binary interaction and undergoes an energetic explosion can match all the observational constraints.

The Sirius AB binary system has masses that are well determined from many decades of astrometric measurements. Because of the well-measured radius and luminosity of Sirius A, we employed the TYCHO stellar evolution code to determine the age of the Sirius A,B binary system accurately, at 225-250 Myr. Note that this fit requires the assumption of solar abundance, and the use of the new Asplund et al. primordial solar metallicity. No fit to Sirius A's position is possible using the old Grevesse & Sauval scale. Because the Sirius B white dwarf parameters have also been determined accurately from space observations, the cooling age could be determined from recent calculations by Fontaine et al. or Wood to be 124 +/- 10 Myr. The difference of the two ages yields the nuclear lifetime and mass of the original primary star, 5.056 +0.374/-0.276 solar masses. This result yields in principle the most accurate data point at relatively high masses for the initial-final mass relation. However, the analysis relies on the assumption that the primordial abundance of the Sirius stars was solar, based on membership in the Sirius supercluster. A recent study suggests that its membership in the group is by no means certain.

Recent multidimensional hydrodynamic simulations have demonstrated the importance of hydrodynamic motions in the convective boundary and radiative regions of stars to transport of energy, momentum, and composition. The impact of these processes increases with stellar mass. Stellar models which approximate this physics have been tested on several classes of observational problems. In this paper we examine the implications of the improved treatment on supernova progenitors. The improved models predict substantially different interior structures. We present pre-supernova conditions and simple explosion calculations from stellar models with and without the improved mixing treatment at 23 solar masses. The results differ substantially.

We examine contributions of second order physical processes to results of stellar evolution calculations amenable to direct observational testing. In the first paper in the series (Young et al. 2001) we established baseline results using only physics which are common to modern stellar evolution codes. In the current paper we establish how much of the discrepancy between observations and baseline models is due to particular elements of new physics. We then consider the impact of the observational uncertainties on the maximum predictive accuracy achievable by a stellar evolution code. The sun is an optimal case because of the precise and abundant observations and the relative simplicity of the underlying stellar physics. The Standard Model is capable of matching the structure of the sun as determined by helioseismology and gross surface observables to better than a percent. Given an initial mass and surface composition within the observational errors, and no additional constraints for which the models can be optimized, it is not possible to predict the sun's current state to better than ~7%. Convectively induced mixing in radiative regions, seen in multidimensional hydrodynamic simulations, dramatically improves the predictions for radii, luminosity, and apsidal motions of eclipsing binaries while simultaneously maintaining consistency with observed light element depletion and turnoff ages in young clusters (Young et al. 2003). Systematic errors in core size for models of massive binaries disappear with more complete mixing physics, and acceptable fits are achieved for all of the binaries without calibration of free parameters. The lack of accurate abundance determinations for binaries is now the main obstacle to improving stellar models using this type of test.

The detection of NaCl at large radii in the Egg Nebula, CRL 2688, requires densities of 10^7 - 10^8 cm^-3 in a thick shell of r ~ a few X 10^17 cm. To explain these results, a mechanism is needed for producing high densities at a considerable distance from the central star. In two dimensional simulations of the interaction of the fast wind with an ambient medium, the material becomes thermally unstable. The resulting clumps can achieve the requisite conditions for NaCl excitation. We present 2D models with simple physics as proof-of-principle calculations to show that the clumping behavior is robust. Clumping is a natural outcome of cooling in the colliding wind model and comparable to that inferred from observations.

We present an analysis of the response of a radiative region to waves generated by a convective region of the star; this wave treatment of the classical problem of ``overshooting'' gives extra mixing relative to the treatment traditionally used in stellar evolutionary codes. The interface between convectively stable and unstable regions is dynamic and nonspherical, so that the nonturbulent material is driven into motion, even in the absence of ``penetrative overshoot.'' These motions may be described by the theory of nonspherical stellar pulsations, and are related to motion measured by helioseismology. Multi-dimensional numerical simulations of convective flow show puzzling features which we explain by this simplified physical model. Gravity waves generated at the interface are dissipated, resulting in slow circulation and mixing seen outside the formal convection zone. The approach may be extended to deal with rotation and composition gradients. Tests of this description in the stellar evolution code TYCHO produce carbon stars on the asymptotic giant branch (AGB), an isochrone age for the Hyades and three young clusters with lithium depletion ages from brown dwarfs, and lithium and beryllium depletion consistent with observations of the Hyades and Pleiades, all without tuning parameters. The insight into the different contributions of rotational and hydrodynamic mixing processes could have important implications for realistic simulation of supernovae and other questions in stellar evolution.

We compare eighteen binary systems with precisely determined radii and masses from 23 to 1.1 M_sol, and stellar evolution models produced with our newly revised code TYCHO. ``Overshooting'' and rotational mixing were suppressed in order to establish a baseline for isolating these and other hydrodynamic effects. Acceptable coeval fits are found for sixteen pairs without optimizing for heavy element or helium abundance. The precision of these tests is limited by the accuracies of the observed effective temperatures. High dispersion spectra and detailed atmospheric modeling should give more accurate effective temperatures and heavy element abundances. PV Cas, a peculiar early A system, EK Cep B, a known post-T Tauri star, and RS Cha, a member of a young OB association, are matched by pre-main sequence models. Predicted mass loss agrees with upper limits from IUE for CW Cep A and B. Relatively poor fits are obtained for binaries having at least one component in the mass range 1.7 < M/M_sol <2.6, whose evolution is sensitive to mixing. These discrepancies are robust and consistent with additional mixing in real stars. The predicted apsidal motion implies that massive star models are systematically less centrally condensed than the real stars. If these effects are due to overshooting, then the overshooting parameter alpha_OV increases with stellar mass. The apsidal motion constants are controlled by radiative opacity under conditions close to those directly measured in laser experiments, making this test more stringent than possible before.

Two dimensional hydrodynamical simulations of convective oxygen burning shell in the presupernova evolution of a 20 solar-mass star are extended to later times. We used the VULCAN code to simulate longer evolution times than previously possible. Our results confirm the previous work of Bazan and Arnett (98) over their time span of 400s. However, at 1200s, we could identify a new steady state that is significantly different than the original one dimensional model. There is considerable overshooting at both the top and bottom boundaries of convection zone. Beyond the boundaries, the convective velocity falls off exponentially, with excitation of internal modes. The resulting mixing greatly affect the evolution of the simulations. Connections with other works of simulation of convection, in which such behavior is found in a different context, are discussed.

The role of hydrodynamic mixing in astrophysics is reviewed, emphasizing connections with laser physics experiments and inertial confinement fusion (ICF). Computer technology now allows two dimensional (2D) simulations, with complex microphysics, of stellar hydrodynamics and evolutionary sequences, and holds the promise for 3D. Careful validation of astrophysical methods, by laboratory experiment, by critical comparison of numerical and analytical methods, and by observation are necessary for the development of simulation methods with reliable predictive capability. Recent and surprising results from isotopic patterns in presolar grains, 2D hydrodynamic simulations of stellar evolution, and laser tests and computer simulations of Richtmeyer-Meshkov and Rayleigh-Taylor instabilities are discussed, and related to stellar evolution and supernovae.

Explosive nucleosynthesis is a combination of the nuclear physics of thermonuclear reactions, and the hydrodynamics of the plasma in which the reactions occur. It depends upon the initial conditions---the stellar evolution up to the explosive instability, and the nature of the explosion mechanism. Some key issues for explosive nucleosynthesis are the interaction of burning with hydrodynamics, the degree of microscopic mixing in convective zones, and the breaking of spherical symmetry by convection and rotation. Recent experiments on high intensity lasers provides new opportunities for laboratory testing of astrophysical hydrodynamic codes. Implications of SN1987A, SN1998bw (GRB980425?), and eta Carina are discussed, as well as the formation of black holes or neutron stars.

The use of supernovae of Type Ia for the determination of accurate distances rests upon the empirical Phillips relation, in which the brightest events are the broadest in time. Implications of new data upon the homogeneity of light curves under the operation of a stretch in time, of the parabolic luminosity increase at the earliest times, and of the time from explosion to maximum light are discussed. The early luminosity is in excellent agreement with predictions (Arnett 1982), and the lack of prominent higher modes of diffusion constrain progenitor and explosion models. Difficulties with reproducing the observed rise time are restricted to radiative transfer models (e.g., Hoeflich and Khokhlov), and probably due to an overestimate of thermal photon escape due to inadequate line lists. Because of the strong dependence of luminosity on 56Ni mass, some simple models can give a Phillips relation of the correct sense.

We present a series of 2-dimensional hydrodynamic simulations of massive disks around protostars. We simulate the same physical problem using both a `Piecewise Parabolic Method' (PPM) code and a `Smoothed Particle Hydrodynamic' (SPH) code, and analyze their differences. The disks studied here range in mass from $0.05 M_*$ to $1.0 M_*$ and in initial minimum Toomre $Q$ value from 1.1 to 3.0. For this problem, the strengths of the codes overlap only in a limited fashion, but similarities exist in their predictions, including spiral arm pattern speeds and morphological features. Our results represent limiting cases (i.e. systems evolved isothermally) rather than true physical systems. Disks become active from the inner regions outward. From the earliest times, their evolution is a strongly dynamic process rather than a smooth progression toward eventual nonlinear behavior. We calculate approximate growth rates for the spiral patterns; the one-armed ($m=1$) spiral arm is not the fastest growing pattern of most disks. In our SPH simulations, disks with initial minimum $Q=1.5$ or lower break up into proto-binary or proto-planetary clumps. However, these simulations cannot follow the physics important for the flow and must be terminated before the system has completely evolved. At their termination, PPM simulations with similar initial conditions show uneven mass distributions within spiral arms, suggesting that clumping behavior might result if they were carried further. Concern that the point-like nature of SPH exaggerates clumping, that our representation of the gravitational potential in PPM is too coarse, and that our physics assumptions are too simple, suggest caution in interpretation of the clumping in both the disk and torus simulations.

By direct hydrodynamic simulation, using the Piecewise Parabolic Method (PPM) code PROMETHEUS, we study the properties of a convective oxygen burning shell in a SN 1987A progenitor star prior to collapse. The convection is too heterogeneous and dynamic to be well approximated by one-dimensional diffusion-like algorithms which have previously been used for this epoch. Qualitatively new phenomena are seen. The simulations are two-dimensional, with good resolution in radius and angle, and use a large (90-degree) slice centered at the equator. The microphysics and the initial model were carefully treated. Many of the qualitative features of previous multi-dimensional simulations of convection are seen, including large kinetic and acoustic energy fluxes, which are not accounted for by mixing length theory. Small but significant amounts of carbon-12 are mixed non-uniformly into the oxygen burning convection zone, resulting in hot spots of nuclear energy production which are more than an order of magnitude more energetic than the oxygen flame itself. Density perturbations (up to 8%) occur at the `edges' of the convective zone and are the result of gravity waves generated by interaction of penetrating flows into the stable region. Perturbations of temperature and electron fraction at the base of the convective zone are of sufficient magnitude to create angular inhomogeneities in explosive nucleosynthesis products, and need to be included in quantitative estimates of yields. Combined with the plume-like velocity structure arising from convection, the perturbations will contribute to the mixing of nickel-56 throughout supernovae envelopes. Runs of different resolution, and angular extent, were performed to test the robustness of these