We provide a detailed description of the Chimera code, a code developed to model core collapse supernovae in multiple spatial dimensions. The core collapse supernova explosion mechanism remains the subject of intense research. Progress to date demonstrates that it involves a complex interplay of neutrino production, transport, and interaction in the stellar core, three-dimensional stellar core fluid dynamics and its associated instabilities, nuclear burning, and the foundational physics of the neutrino-stellar core weak interactions and the equations of state of all stellar core constituents -particularly, the nuclear equation of state associated with nucleons, both free and bound in nuclei. Chimera, by incorporating detailed neutrino transport, realistic neutrino-matter interactions, three-dimensional hydrodynamics, realistic nuclear, leptonic, and photonic equations of state, and a nuclear reaction network, along with other refinements, can be used to study the role of neutrino radiation, hydrodynamic instabilities, and a variety of input physics in the explosion mechanism itself. It can also be used to compute observables such as neutrino signatures, gravitational radiation, and the products of nucleosynthesis associated with core collapse supernovae. The code contains modules for neutrino transport, multidimensional compressible hydrodynamics, nuclear reactions, a variety of neutrino interactions, equations of state, and modules to provide data for post-processing observables such as the products of nucleosynthesis, and gravitational radiation. Chimera is an evolving code, being updated periodically with improved input physics and numerical refinements. We detail here the current version of the code, from which future improvements will stem, which can in turn be described as needed in future publications.

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.

Building on the framework of Zhang \& Shu \cite{zhangShu_2010a,zhangShu_2010b}, we develop a realizability-preserving method to simulate the transport of particles (fermions) through a background material using a two-moment model that evolves the angular moments of a phase space distribution function $f$. The two-moment model is closed using algebraic moment closures; e.g., as proposed by Cernohorsky \& Bludman \cite{cernohorskyBludman_1994} and Banach \& Larecki \cite{banachLarecki_2017a}. Variations of this model have recently been used to simulate neutrino transport in nuclear astrophysics applications, including core-collapse supernovae and compact binary mergers. We employ the discontinuous Galerkin (DG) method for spatial discretization (in part to capture the asymptotic diffusion limit of the model) combined with implicit-explicit (IMEX) time integration to stably bypass short timescales induced by frequent interactions between particles and the background. Appropriate care is taken to ensure the method preserves strict algebraic bounds on the evolved moments (particle density and flux) as dictated by Pauli's exclusion principle, which demands a bounded distribution function (i.e., $f\in[0,1]$). This realizability-preserving scheme combines a suitable CFL condition, a realizability-enforcing limiter, a closure procedure based on Fermi-Dirac statistics, and an IMEX scheme whose stages can be written as a convex combination of forward Euler steps combined with a backward Euler step. Numerical results demonstrate the realizability-preserving properties of the scheme. We also demonstrate that the use of algebraic moment closures not based on Fermi-Dirac statistics can lead to unphysical moments in the context of fermion transport.

We present a comparison between several simulation codes designed to study the core-collapse supernova mechanism. We pay close attention to controlling the initial conditions and input physics in order to ensure a meaningful and informative comparison. Our goal is three-fold. First, we aim to demonstrate the current level of agreement between various groups studying the core-collapse supernova central engine. Second, we desire to form a strong basis for future simulation codes and methods to compare to. Lastly, we want this work to be a stepping stone for future work exploring more complex simulations of core-collapse supernovae, i.e., simulations in multiple dimensions and simulations with modern neutrino and nuclear physics. We compare the early (first ~500ms after core bounce) spherically-symmetric evolution of a 20 solar mass progenitor star from six different core-collapse supernovae codes: 3DnSNe-IDSA, AGILE-BOLTZTRAN, FLASH, F{\sc{ornax}}, GR1D, and PROMETHEUS-VERTEX. Given the diversity of neutrino transport and hydrodynamic methods employed, we find excellent agreement in many critical quantities, including the shock radius evolution and the amount of neutrino heating. Our results provide an excellent starting point from which to extend this comparison to higher dimensions and compare the development of hydrodynamic instabilities that are crucial to the supernova explosion mechanism, such as turbulence and convection.

In this work we report briefly on the gravitational wave (GW) signal computed in the context of a self-consistent, 3D simulation of a core-collapse supernova (CCSN) explosion of a 15M$_\odot$ progenitor star. We present a short overview of the GW signal, including signal amplitude, frequency distribution, and the energy emitted in the form of GWs for each phase of explosion, along with neutrino luminosities, and discuss correlations between them.

In this Letter, we report on the gravitational wave signal computed in the context of an $ab$ $initio$, three-dimensional simulation of a core collapse supernova explosion, beginning with a 15M$_\odot$ star and using state-of-the-art weak interactions. The simulation was performed with our neutrino hydrodynamics code Chimera. We discuss the potential for detection of our predicted gravitational signal by the current generation of gravitational wave detectors.

The DUNE/LBNF program aims to address key questions in neutrino physics and astroparticle physics. Realizing DUNE's potential to reconstruct low-energy particles in the 10-100 MeV energy range will bring significant benefits for all DUNE's science goals. In neutrino physics, low-energy sensitivity will improve neutrino energy reconstruction in the GeV range relevant for the kinematics of DUNE's long-baseline oscillation program. In astroparticle physics, low-energy capabilities will make DUNE's far detectors the world's best apparatus for studying the electron-neutrino flux from a supernova. This will open a new window to unrivaled studies of the dynamics and neutronization of a star's central core in real time, the potential discovery of the neutrino mass hierarchy, provide new sensitivity to physics beyond the Standard Model, and evidence of neutrino quantum-coherence effects. The same capabilities will also provide new sensitivity to `boosted dark matter' models that are not observable in traditional direct dark matter detectors.

This white paper informs the nuclear astrophysics community and funding agencies about the scientific directions and priorities of the field and provides input from this community for the 2015 Nuclear Science Long Range Plan. It summarizes the outcome of the nuclear astrophysics town meeting that was held on August 21-23, 2014 in College Station at the campus of Texas A&M University in preparation of the NSAC Nuclear Science Long Range Plan. It also reflects the outcome of an earlier town meeting of the nuclear astrophysics community organized by the Joint Institute for Nuclear Astrophysics (JINA) on October 9- 10, 2012 Detroit, Michigan, with the purpose of developing a vision for nuclear astrophysics in light of the recent NRC decadal surveys in nuclear physics (NP2010) and astronomy (ASTRO2010). The white paper is furthermore informed by the town meeting of the Association of Research at University Nuclear Accelerators (ARUNA) that took place at the University of Notre Dame on June 12-13, 2014. In summary we find that nuclear astrophysics is a modern and vibrant field addressing fundamental science questions at the intersection of nuclear physics and astrophysics. These questions relate to the origin of the elements, the nuclear engines that drive life and death of stars, and the properties of dense matter. A broad range of nuclear accelerator facilities, astronomical observatories, theory efforts, and computational capabilities are needed. With the developments outlined in this white paper, answers to long standing key questions are well within reach in the coming decade.

Core-collapse supernovae, the culmination of massive stellar evolution, are spectacular astronomical events and the principle actors in the story of our elemental origins. Our understanding of these events, while still incomplete, centers around a neutrino-driven central engine that is highly hydrodynamically unstable. Increasingly sophisticated simulations reveal a shock that stalls for hundreds of milliseconds before reviving. Though brought back to life by neutrino heating, the development of the supernova explosion is inextricably linked to multi-dimensional fluid flows. In this paper, the outcomes of three-dimensional simulations that include sophisticated nuclear physics and spectral neutrino transport are juxtaposed to learn about the nature of the three dimensional fluid flow that shapes the explosion. Comparison is also made between the results of simulations in spherical symmetry from several groups, to give ourselves confidence in the understanding derived from this juxtaposition.

We present four ab initio axisymmetric core-collapse supernova simulations for 12, 15, 20, and 25 $M_\odot$ progenitors. All of the simulations yield explosions and have been evolved for at least 1.2 seconds after core bounce and 1 second after material first becomes unbound. Simulations were computed with our Chimera code employing spectral neutrino transport, special and general relativistic transport effects, and state-of-the-art neutrino interactions. Continuing the evolution beyond 1 second allows explosions to develop more fully and the processes powering the explosions to become more clearly evident. We compute explosion energy estimates, including the binding energy of the stellar envelope outside the shock, of 0.34, 0.88, 0.38, and 0.70 B ($10^{51}$ ergs) and increasing at 0.03, 0.15, 0.19, and 0.52 B s$^{-1}$, respectively, for the 12, 15, 20, and 25 $M_\odot$ models. Three models developed pronounced prolate shock morphologies, while the 20 $M_\odot$ model, though exhibiting lobes and accretion streams like the other models, develops an approximately spherical, off-center shock as the explosion begins and then becomes moderately prolate $\sim$600 ms after bounce. This reduces the explosion energy relative to the other models by reducing mass accretion during the critical explosion power-up phase. We examine the growth of the explosion energy in our models through detailed analyses of the energy sources and flows. We find that the 12 and 20 $M_\odot$ models have explosion energies comparable to that of the lower range of observed explosion energies while the 15 and 25 $M_\odot$ models are within the range of observed explosion energies, particularly considering the rate at which their explosion energies are increasing. The ejected $^{56}$Ni masses given by our models are all within observational limits as are the proto-neutron star masses and kick velocities. (Truncated)

We present gravitational wave and neutrino signatures obtained in our first principle 3D core-collapse supernova simulation of 15M non-rotating progenitor with Chimera code. Observations of neutrinos emitted by the forming neutron star, and gravitational waves, which are produced by hydrodynamic instabilities is the only way to get direct information about the supernova engine. Both GW and neutrino signals show different phases of supernova evolution.

We present results from an ab initio three-dimensional, multi-physics core collapse supernova simulation for the case of a 15 M progenitor. Our simulation includes multi-frequency neutrino transport with state-of-the-art neutrino interactions in the "ray-by-ray" approximation, and approximate general relativity. Our model exhibits a neutrino-driven explosion. The shock radius begins an outward trajectory at approximately 275 ms after bounce, giving the first indication of a developing explosion in the model. The onset of this shock expansion is delayed relative to our two-dimensional counterpart model, which begins at approximately 200 ms after core bounce. At a time of 441 ms after bounce, the angle-averaged shock radius in our three-dimensional model has reached 751 km. Further quantitative analysis of the outcomes in this model must await further development of the post-bounce dynamics and a simulation that will extend well beyond 1 s after stellar core bounce, based on the results for the same progenitor in the context of our two-dimensional, counterpart model. This more complete analysis will determine whether or not the explosion is robust and whether or not observables such as the explosion energy, 56Ni mass, etc. are in agreement with observations. Nonetheless, the onset of explosion in our ab initio three-dimensional multi-physics model with multi-frequency neutrino transport and general relativity is encouraging.

In this paper, the second in a series, we document the algorithms and solvers for compressible nonrelativistic hydrodynamics implemented in GenASiS (General Astrophysical Simulation System)---a new code being developed initially and primarily, though by no means exclusively, for the simulation of core-collapse supernovae. In the Mathematics division of GenASiS we introduce Solvers, which includes finite-volume updates for generic hyperbolic BalanceEquations and ordinary differential equation integration Steps. We also introduce the Physics division of GenASiS; this extends the Manifolds division of Mathematics into physical Spaces, defines StressEnergies, and combines these into Universes. We benchmark the hydrodynamics capabilities of GenASiS against many standard test problems; the results illustrate the basic competence of our implementation, demonstrate the manifest superiority of the HLLC over the HLL Riemann solver in a number of interesting cases, and provide preliminary indications of the code's ability to scale and to function with cell-by-cell fixed-mesh refinement.

GenASiS (General Astrophysical Simulation System) is a new code being developed initially and primarily, though by no means exclusively, for the simulation of core-collapse supernovae on the world's leading capability supercomputers. This paper---the first in a series---demonstrates a centrally refined coordinate patch suitable for gravitational collapse and documents methods for compressible nonrelativistic hydrodynamics. We benchmark the hydrodynamics capabilities of GenASiS against many standard test problems; the results illustrate the basic competence of our implementation, demonstrate the strengths and limitations of the HLLC relative to the HLL Riemann solver in a number of interesting cases, and provide preliminary indications of the code's ability to scale and to function with cell-by-cell fixed-mesh refinement.

We have performed ab initio neutrino radiation hydrodynamics simulations in three and two spatial dimensions (3D and 2D) of core-collapse supernovae from the same 15 $M_\odot$ progenitor through 440 ms after core bounce. Both 3D and 2D models achieve explosions, however, the onset of explosion (shock revival) is delayed by $\sim$100 ms in 3D relative to the 2D counterpart and the growth of the diagnostic explosion energy is slower. This is consistent with previously reported 3D simulations utilizing iron-core progenitors with dense mantles. In the $\sim$100 ms before the onset of explosion, diagnostics of neutrino heating and turbulent kinetic energy favor earlier explosion in 2D. During the delay, the angular scale of convective plumes reaching the shock surface grows and explosion in 3D is ultimately lead by a single, large-angle plume, giving the expanding shock a directional orientation not dissimilar from those imposed by axial symmetry in 2D simulations. We posit that shock revival and explosion in the 3D simulation may be delayed until sufficiently large plumes form, whereas such plumes form more rapidly in 2D, permitting earlier explosions.

We present the gravitational waveforms computed in ab initio two-dimensional core collapse supernova models evolved with the Chimera code for progenitor masses between 12 and 25 solar masses. All models employ multi-frequency neutrino transport in the ray-by-ray approximation, state-of-the-art weak interaction physics, relativistic transport corrections such as the gravitational redshift of neutrinos, two-dimensional hydrodynamics with the commensurate relativistic corrections, Newtonian self-gravity with a general relativistic monopole correction, and the Lattimer-Swesty equation of state with 220 MeV compressibility, and begin with the most recent Woosley-Heger nonrotating progenitors in this mass range. All of our models exhibit robust explosions. Therefore, our waveforms capture all stages of supernova development: 1) a relatively short and weak prompt signal, 2) a quiescent stage, 3) a strong signal due to convection and SASI activity, 4) termination of active accretion onto the proto-neutron star, and 5) a slowly increasing tail that reaches a saturation value. Fourier decomposition shows that the gravitational wave signals we predict should be observable by AdvLIGO for Galactic events across the range of progenitors considered here. The fundamental limitation of these models is in their imposition of axisymmetry. Further progress will require counterpart three-dimensional models, which are underway.

We have been working within the fundamental paradigm that core collapse supernovae (CCSNe) may be neutrino driven, since the first suggestion of this by Colgate and White nearly five decades ago. Computational models have become increasingly sophisticated, first in one spatial dimension assuming spherical symmetry, then in two spatial dimensions assuming axisymmetry, and now in three spatial dimensions with no imposed symmetries. The increase in the number of spatial dimensions has been accompanied by an increase in the physics included in the models, and an increase in the sophistication with which this physics has been modeled. Computation has played an essential role in the development of CCSN theory, not simply for the obvious reason that such multidimensional, multi-physics, nonlinear events cannot possibly be fully captured analytically, but for its role in discovery. In particular, the discovery of the standing accretion shock instability (SASI) through computation about a decade ago has impacted all simulations performed since then. Today, we appear to be at a threshold, where neutrinos, neutrino-driven convection, and the SASI, working together over time scales significantly longer than had been anticipated in the past, are able to generate explosions, and in some cases, robust explosions, in a number of axisymmetric models. But how will this play out in three dimensions? Early results from the first three-dimensional (3D), multi-physics simulation of the "Oak Ridge" group are promising. I will discuss the essential components of today's models and the requirements of realistic CCSN modeling, present results from our one-, two-, and three-dimensional models, place our models in context with respect to other efforts around the world, and discuss short- and long-term next steps.

We investigate core-collapse supernova (CCSN) nucleosynthesis in polar axisymmetric simulations using the multidimensional radiation hydrodynamics code CHIMERA. Computational costs have traditionally constrained the evolution of the nuclear composition in CCSN models to, at best, a 14-species $\alpha$-network. Such a simplified network limits the ability to accurately evolve detailed composition, neutronization and the nuclear energy generation rate. Lagrangian tracer particles are commonly used to extend the nuclear network evolution by incorporating more realistic networks in post-processing nucleosynthesis calculations. Limitations such as poor spatial resolution of the tracer particles, estimation of the expansion timescales, and determination of the "mass-cut" at the end of the simulation impose uncertainties inherent to this approach. We present a detailed analysis of the impact of these uncertainties on post-processing nucleosynthesis calculations and implications for future models.

We extend the positivity-preserving method of Zhang & Shu (2010, JCP, 229, 3091-3120) to simulate the advection of neutral particles in phase space using curvilinear coordinates. The ability to utilize these coordinates is important for non-equilibrium transport problems in general relativity and also in science and engineering applications with specific geometries. The method achieves high-order accuracy using Discontinuous Galerkin (DG) discretization of phase space and strong stability-preserving, Runge-Kutta (SSP-RK) time integration. Special care in taken to ensure that the method preserves strict bounds for the phase space distribution function $f$; i.e., $f\in[0,1]$. The combination of suitable CFL conditions and the use of the high-order limiter proposed in Zhang & Shu (2010) is sufficient to ensure positivity of the distribution function. However, to ensure that the distribution function satisfies the upper bound, the discretization must, in addition, preserve the divergence-free property of the phase space flow. Proofs that highlight the necessary conditions are presented for general curvilinear coordinates, and the details of these conditions are worked out for some commonly used coordinate systems (i.e., spherical polar spatial coordinates in spherical symmetry and cylindrical spatial coordinates in axial symmetry, both with spherical momentum coordinates). Results from numerical experiments --- including one example in spherical symmetry adopting the Schwarzschild metric --- demonstrate that the method achieves high-order accuracy and that the distribution function satisfies the maximum principle.

We summarize the results of core collapse supernova theory from one-, two-, and three-dimensional models and provide a snapshot of the field at this time. We also present results from the "Oak Ridge" group in this context. Studies in both one and two spatial dimensions define the necessary} physics that must be included in core collapse supernova models: a general relativistic treatment of gravity (at least an approximate one), spectral neutrino transport, including relativistic effects such as gravitational redshift, and a complete set of neutrino weak interactions that includes state-of-the-art electron capture on nuclei and energy-exchanging scattering on electrons and nucleons. Whether or not the necessarily approximate treatment of this physics in current models that include it is sufficient remains to be determined in the context of future models that remove the approximations. We summarize the results of the Oak Ridge group's two-dimensional supernova models. In particular, we demonstrate that robust neutrino-driven explosions can be obtained. We also demonstrate that our predictions of the explosion energies and remnant neutron star masses are in agreement with observations, although a much larger number of models must be developed before more confident conclusions can be made. We provide preliminary results from our ongoing three dimensional model with the same physics. Finally, we speculate on future outcomes and directions.

We present a new derivation of the conservative form of the general relativistic Boltzmann equation and specialize it to the 3+1 metric. The resulting transport equation is intended for use in simulations involving numerical relativity, particularly in the absence of spherical symmetry. The independent variables are lab frame coordinate basis spacetime position components and comoving frame curvilinear momentum space coordinates. With an eye towards astrophysical applications---such as core-collapse supernovae and compact object mergers---in which the fluid includes nuclei and/or nuclear matter at finite temperature, and in which the transported particles are neutrinos, we examine the relationship between lepton number and four-momentum exchange between neutrinos and the fluid.

We present conservative 3+1 general relativistic variable Eddington tensor radiation transport equations, including greater elaboration of the momentum space divergence (that is, the energy derivative term) than in previous work. These equations are intended for use in simulations involving numerical relativity, particularly in the absence of spherical symmetry. The independent variables are the lab frame coordinate basis spacetime position coordinates and the particle energy measured in the comoving frame. With an eye towards astrophysical applications---such as core-collapse supernovae and compact object mergers---in which the fluid includes nuclei and/or nuclear matter at finite temperature, and in which the transported particles are neutrinos, we pay special attention to the consistency of four-momentum and lepton number exchange between neutrinos and the fluid, showing the term-by-term cancellations that must occur for this consistency to be achieved.

We present an overview of four ab initio axisymmetric core-collapse supernova simulations employing detailed spectral neutrino transport computed with our CHIMERA code and initiated from Woosley & Heger (2007) progenitors of mass 12, 15, 20, and 25 M_sol. All four models exhibit shock revival over \sim 200 ms (leading to the possibility of explosion), driven by neutrino energy deposition. Hydrodynamic instabilities that impart substantial asymmetries to the shock aid these revivals, with convection appearing first in the 12 M_sol model and the standing accretion shock instability (SASI) appearing first in the 25 M_sol model. Three of the models have developed pronounced prolate morphologies (the 20 M_sol model has remained approximately spherical). By 500 ms after bounce the mean shock radii in all four models exceed 3,000 km and the diagnostic explosion energies are 0.33, 0.66, 0.65, and 0.70 Bethe (B = $10^{51}$ ergs) for the 12, 15, 20, and 25 M_sol models, respectively, and are increasing. The three least massive of our models are already sufficiently energetic to completely unbind the envelopes of their progenitors (i.e., to explode), as evidenced by our best estimate of their explosion energies, which first become positive at 320, 380, and 440 ms after bounce. By 850 ms the 12 M_sol diagnostic explosion energy has saturated at 0.38 B, and our estimate for the final kinetic energy of the ejecta is \sim 0.3 B, which is comparable to observations for lower-mass progenitors.

Ascertaining the core-collapse supernova mechanism is a complex, and yet unsolved, problem dependent on the interaction of general relativity, hydrodynamics, neutrino transport, neutrino-matter interactions, and nuclear equations of state and reaction kinetics. Ab initio modeling of core-collapse supernovae and their nucleosynthetic outcomes requires care in the coupling and approximations of the physical components. We have built our multi-physics CHIMERA code for supernova modeling in 1-, 2-, and 3-D, using ray-by-ray neutrino transport, approximate general relativity, and detailed neutrino and nuclear physics. We discuss some early results from our current series of exploding 2D simulations and our work to perform computationally tractable simulations in 3D using the "Yin-Yang" grid.

We derive conservative, multidimensional, energy-dependent moment equations for neutrino transport in core-collapse supernovae and related astrophysical systems, with particular attention to the consistency of conservative four-momentum and lepton number transport equations. After taking angular moments of conservative formulations of the general relativistic Boltzmann equation, we specialize to a conformally flat spacetime, which also serves as the basis for four further limits. Two of these---the multidimensional special relativistic case, and a conformally flat formulation of the spherically symmetric general relativistic case---are given in appendices for the sake of comparison with extant literature. The third limit is a weak-field, `pseudo-Newtonian' approach \citep{kim_etal_2009,kim_etal_2012} in which the source of the gravitational potential includes the trace of the stress-energy tensor (rather than just the mass density), and all orders in fluid velocity $v$ are retained. Our primary interest here is in the fourth limit: `$\mathcal{O}(v)$' moment equations for use in conjunction with Newtonian self-gravitating hydrodynamics. We show that the concept of `$\mathcal{O}(v)$' transport requires care when dealing with both conservative four-momentum and conservative lepton number transport, and present two self-consistent options: `$\mathcal{O}(v)$-plus' transport, in which an $\mathcal{O}(v^2)$ energy equation combines with an $\mathcal{O}(v)$ momentum equation to give an $\mathcal{O}(v^2)$ number equation; and `$\mathcal{O}(v)$-minus' transport, in which an $\mathcal{O}(v)$ energy equation combines with an $\mathcal{O}(1)$ momentum equation to give an $\mathcal{O}(v)$ number equation.

We have conducted a series of numerical experiments using spherically symmetric, general relativistic, neutrino radiation hydrodynamics with the code Agile-BOLTZTRAN to examine the effects of modern neutrino opacities on the development of supernova simulations. We test the effects of opacities by removing opacities or by undoing opacity improvements for individual opacities and groups of opacities. We find that improvements to electron capture (EC) on nuclei, namely EC on an ensemble of nuclei using modern nuclear structure models rather than the simpler independent-particle approximation (IPA) for EC on a mean nucleus, plays the most important role during core collapse of all tested neutrino opacities. Low-energy neutrinos emitted by modern nuclear EC preferentially escape during collapse without the energy downscattering on electrons required to enhance neutrino escape and deleptonization for the models with IPA nuclear EC. During shock breakout the primary influence on the emergent neutrinos arises from NIS on electrons. For the accretion phase, non-isoenergetic scattering on free nucleons and pair emission by $e^+e^-$ annihilation have the largest impact on the neutrino emission and shock evolution. Other opacities evaluated, including nucleon--nucleon bremsstrahlung and especially neutrino--positron scattering, have little measurable impact on neutrino emission or shock dynamics. Modern treatments of nuclear electron capture, $e^+e^-$-annihilation pair emission, and non-isoenergetic scattering on electrons and free nucleons are critical elements of core-collapse simulations of all dimensionality.

The stationary accretion shock instability (SASI) plays a central role in modern simulations of the explosion phase of core-collapse supernovae (CCSNe). It may be key to realizing neutrino powered explosions, and possibly links birth properties of pulsars (e.g., kick, spin, and magnetic field) to supernova dynamics. Using high-resolution magnetohydrodynamic simulations, we study the development of turbulence, and subsequent amplification of magnetic fields in a simplified model of the post-bounce core-collapse supernova environment. Turbulence develops from secondary instabilities induced by the SASI. Our simulations suggest that the development of turbulence plays an important role for the subsequent evolution of the SASI. The turbulence also acts to amplify weak magnetic fields via a small-scale dynamo.

We describe the initial implementation of magnetohydrodynamics (MHD) in our astrophysical simulation code \genasis. Then, we present MHD simulations exploring the capacity of the stationary accretion shock instability (SASI) to generate magnetic fields by adding a weak magnetic field to an initially spherically symmetric fluid configuration that models a stalled shock in the post-bounce supernova environment. Upon perturbation and nonlinear SASI development, shear flows associated with the spiral SASI mode contributes to a widespread and turbulent field amplification mechanism. While the SASI may contribute to neutron star magnetization, these simulations do not show qualitatively new features in the global evolution of the shock as a result of SASI-induced magnetic field amplification.

We extend our investigation of magnetic field evolution in three-dimensional flows driven by the stationary accretion shock instability (SASI) with a suite of higher-resolution idealized models of the post-bounce core-collapse supernova environment. Our magnetohydrodynamic simulations vary in initial magnetic field strength, rotation rate, and grid resolution. Vigorous SASI-driven turbulence inside the shock amplifies magnetic fields exponentially; but while the amplified fields reduce the kinetic energy of small-scale flows, they do not seem to affect the global shock dynamics. The growth rate and final magnitude of the magnetic energy are very sensitive to grid resolution, and both are underestimated by the simulations. Nevertheless our simulations suggest that neutron star magnetic fields exceeding $10^{14}$ G can result from dynamics driven by the SASI, \emph{even for non-rotating progenitors}.

We have conducted a series of numerical experiments with the spherically symmetric, general relativistic, neutrino radiation hydrodynamics code Agile-BOLTZTRAN to examine the effects of several approximations used in multidimensional core-collapse supernova simulations. Our code permits us to examine the effects of these approximations quantitatively by removing, or substituting for, the pieces of supernova physics of interest. These approximations include: (1) using Newtonian versus general relativistic gravity, hydrodynamics, and transport; (2) using a reduced set of weak interactions, including the omission of non-isoenergetic neutrino scattering, versus the current state-of-the-art; and (3) omitting the velocity-dependent terms, or observer corrections, from the neutrino Boltzmann kinetic equation. We demonstrate that each of these changes has noticeable effects on the outcomes of our simulations. Of these, we find that the omission of observer corrections is particularly detrimental to the potential for neutrino-driven explosions and exhibits a failure to conserve lepton number. Finally, we discuss the impact of these results on our understanding of current, and the requirements for future, multidimensional models.

Core-collapse supernova models depend on the details of the nuclear and weak interaction physics inputs just as they depend on the details of the macroscopic physics (transport, hydrodynamics, etc.), numerical methods, and progenitors. We present preliminary results from our ongoing comparison studies of nuclear and weak interaction physics inputs to core collapse supernova models using the spherically-symmetric, general relativistic, neutrino radiation hydrodynamics code Agile-Boltztran. We focus on comparisons of the effects of the nuclear EoS and the effects of improving the opacities, particularly neutrino--nucleon interactions.

The CHIMERA code is a multi-dimensional multi-physics engine dedicated primarily to the simulation of core collapse supernova explosions. One of the most important aspects of these explosions is their capacity to produce gravitational radiation that is detectable by Earth-based laser-interferometric gravitational wave observatories such as LIGO and VIRGO. We present here preliminary gravitational signatures of two-dimensional models with non-rotating progenitors. These simulations exhibit explosions, which are followed for more than half a second after stellar core bounce.

Massive stars end their life in an explosion event with kinetic energies of the order 1 Bethe. Immediately after the explosion has been launched, a region of low density and high entropy forms behind the ejecta which is continuously subject to neutrino heating. The neutrinos emitted from the remnant at the center, the protoneutron star (PNS), heat the material above the PNS surface. This heat is partly converted into kinetic energy and the material accelerates to an outflow that is known as the neutrino driven wind. For the first time, we simulate the collapse, bounce, explosion and the neutrino driven wind phases consistently over more than 20 seconds. Our numerical model is based on spherically symmetric general relativistic radiation hydrodynamics using spectral three flavor Boltzmann neutrino transport. In simulations where no explosions are obtained naturally, we model neutrino driven explosions for low and intermediate mass Fe-core progenitor stars by enhancing the charged current reaction rates. In the case of a special progenitor star, the O-Ne-Mg-core, the explosion in spherical symmetry was obtained without enhanced opacities. The post explosion evolution is in qualitative agreement with static steady-state and parametrized dynamic models of the neutrino driven wind. On the other hand, we find generally smaller neutrino luminosities and mean neutrino energies as well as a different evolutionary behavior of the neutrino luminosities and mean neutrino energies. The neutrino driven wind is proton-rich for more than 10 seconds and the contraction of the PNS differs from the assumptions made for the conditions at the inner boundary in previous neutrino driven wind studies. Despite the moderately large entropies per baryon of about 100 and the fast expansion timescale, the conditions found in our model are unlikely to favor...

For finite chemical potential effective models of QCD predict a first order phase transition. In favour for the search of such a phase transition in nature, we construct an equation of state for strange quark matter based on the MIT bag model. We apply this equation of state to highly asymmetric core collapse supernova matter with finite temperatures and large baryon densities. The phase transition is constructed using the general Gibbs conditions, which results in an extended coexistence region between the pure hadronic and pure quark phases in the phase diagram, i.e. the mixed phase. The supernovae are simulated via general relativistic radiation hydrodynamics based on three flavor Boltzmann neutrino transport in spherical symmetry. During the dynamical evolution temperatures above 10 MeV, baryon densities above nuclear saturation density and a proton-to-baryon ratio below 0.2 are obtained. At these conditions the phase transition is triggered which leads to a significant softening of the EoS for matter in the mixed phase. As a direct consequence of the stiffening of the EoS again for matter in the pure quark phase, a shock wave forms at the boundary between the mixed and the pure hadronic phases. This shock is accelerated and propagates outward which releases a burst of neutrinos dominated by electron anti-neutrinos due to the lifted degeneracy of the shock-heated hadronic material. We discuss the radiation-hydrodynamic evolution of the phase transition at the example of several low and intermediate mass Fe-core progenitor stars and illustrate the expected neutrino signal from the phase transition.

We present the gravitational wave signatures for a suite of axisymmetric core collapse supernova models with progenitors masses between 12 and 25 solar masses. These models are distinguished by the fact they explode and contain essential physics (in particular, multi-frequency neutrino transport and general relativity) needed for a more realistic description. Thus, we are able to compute complete waveforms (i.e., through explosion) based on non-parameterized, first-principles models. This is essential if the waveform amplitudes and time scales are to be computed more precisely. Fourier decomposition shows that the gravitational wave signals we predict should be observable by AdvLIGO across the range of progenitors considered here. The fundamental limitation of these models is in their imposition of axisymmetry. Further progress will require counterpart three-dimensional models.

We begin an exploration of the capacity of the stationary accretion shock instability (SASI) to generate magnetic fields by adding a weak, stationary, and radial (but bipolar) magnetic field, and in some cases rotation, to an initially spherically symmetric fluid configuration that models a stalled shock in the post-bounce supernova environment. In axisymmetric simulations we find that cycles of latitudinal flows into and radial flows out of the polar regions amplify the field parallel to the symmetry axis, typically increasing the total magnetic energy by about two orders of magnitude. Nonaxisymmetric calculations result in fundamentally different flows and a larger magnetic energy increase: shearing associated with the SASI spiral mode contributes to a widespread and turbulent field amplification mechanism, boosting the magnetic energy by almost four orders of magnitude (a result which remains very sensitive to the spatial resolution of the numerical simulations). While the SASI may contribute to neutron star magnetization, these simulations do not show qualitatively new features in the global evolution of the shock as a result of SASI-induced magnetic field amplification.

Unraveling the mechanism for core-collapse supernova explosions is an outstanding computational challenge and the problem remains essentially unsolved despite more than four decades of effort. However, much progress in realistic modeling has occurred recently through the availability of multi-teraflop machines and the increasing sophistication of supernova codes. These improvements have led to some key insights which may clarify the picture in the not too distant future. Here we briefly review the current status of the three explosion mechanisms (acoustic, MHD, and neutrino heating) that are currently under active investigation, concentrating on the neutrino heating mechanism as the one most likely responsible for producing explosions from progenitors in the mass range ~10 to ~25 solar masses. We then briefly describe the CHIMERA code, a supernova code we have developed to simulate core-collapse supernovae in 1, 2, and 3 spatial dimensions. We finally describe the results of an ongoing suite of 2D simulations initiated from a 12, 15, 20, and 25 solar mass progenitor. These have all exhibited explosions and are currently in the expanding phase with the shock at between 5,000 and 10,000 km. We finally very briefly describe an ongoing simulation in 3 spatial dimensions initiated from the 15 solar mass progenitor.

Much progress in realistic modeling of core-collapse supernovae has occurred recently through the availability of multi-teraflop machines and the increasing sophistication of supernova codes. These improvements are enabling simulations with enough realism that the explosion mechanism, long a mystery, may soon be delineated. We briefly describe the CHIMERA code, a supernova code we have developed to simulate core-collapse supernovae in 1, 2, and 3 spatial dimensions. We then describe the results of an ongoing suite of 2D simulations initiated from a 12, 15, 20, and 25 solar mass progenitor. These have all exhibited explosions and are currently in the expanding phase with the shock at between 5,000 and 20,000 km. We also briefly describe an ongoing simulation in 3 spatial dimensions initiated from the 15 solar mass progenitor.

The implications of the formation of strange quark matter in neutron stars and in core-collapse supernovae is discussed with special emphasis on the possibility of having a strong first order QCD phase transition at high baryon densities. If strange quark matter is formed in core-collapse supernovae shortly after the bounce, it causes the launch of a second outgoing shock which is energetic enough to lead to a explosion. A signal for the formation of strange quark matter can be read off from the neutrino spectrum, as a second peak in antineutrinos is released when the second shock runs over the neutrinosphere.

We discuss the formation of stellar mass black holes via protoneutron star (PNS) collapse. In the absence of an earlier explosion, the PNS collapses to a black hole due to the continued mass accretion onto the PNS. We present an analysis of the emitted neutrino spectra of all three flavors during the PNS contraction. Special attention is given to the physical conditions which depend on the input physics, e.g. the equation of state (EoS) and the progenitor model. The PNSs are modeled as the central object in core collapse simulations using general relativistic three-flavor Boltzmann neutrino transport in spherical symmetry. The simulations are launched from several massive progenitors of 40 and 50 solar mass. We analyze the electron-neutrino luminosity dependencies and construct a simple approximation for the electron-neutrino luminosity, which depends only on the physical conditions at the electron-neutrinosphere. In addition, we analyze different mu/tau-neutrino pair-reactions separately and compare the differences during the post-bounce phase of failed core collapse supernova explosions of massive progenitors. We also investigate the connection between the increasing mu/tau-neutrino luminosity and the PNS contraction during the accretion phase before black hole formation. Comparing the different post bounce phase of the progenitor models under investigation, we find large differences in the emitted neutrino spectra. These differences and the analysis of the electron-neutrino luminosity indicate a strong progenitor model dependency of the emitted neutrino signal.

The possible role of a first order QCD phase transition at nonvanishing quark chemical potential and temperature for cold neutron stars and for supernovae is delineated. For cold neutron stars, we use the NJL model with nonvanishing color superconducting pairing gaps, which describes the phase transition to the 2SC and the CFL quark matter phases at high baryon densities. We demonstrate that these two phase transitions can both be present in the core of neutron stars and that they lead to the appearance of a third family of solution for compact stars. In particular, a core of CFL quark matter can be present in stable compact star configurations when slightly adjusting the vacuum pressure to the onset of the chiral phase transition from the hadronic model to the NJL model. We show that a strong first order phase transition can have strong impact on the dynamics of core collapse supernovae. If the QCD phase transition sets in shortly after the first bounce, a second outgoing shock wave can be generated which leads to an explosion. The presence of the QCD phase transition can be read off from the neutrino and antineutrino signal of the supernova.

We explore the implications of the QCD phase transition during the postbounce evolution of core-collapse supernovae. Using the MIT bag model for the description of quark matter and assuming small bag constants, we find that the phase transition occurs during the early postbounce accretion phase. This stage of the evolution can be simulated with general relativistic three-flavor Boltzmann neutrino transport. The phase transition produces a second shock wave that triggers a delayed supernova explosion. If such a phase transition happens in a future galactic supernova, its existence and properties should become observable as a second peak in the neutrino signal that is accompanied by significant changes in the energy of the emitted neutrinos. In contrast to the first neutronization burst, this second neutrino burst is dominated by the emission of anti-neutrinos because the electron-degeneracy is lifted when the second shock passes through the previously neutronized matter.

By adding a weak magnetic field to a spherically symmetric fluid configuration that caricatures a stalled shock in the post-bounce supernova environment, we explore the capacity of the stationary accretion shock instability (SASI) to generate magnetic fields. The SASI develops upon perturbation of the initial condition, and the ensuing flow generates--{\em in the absence of rotation}--dynamically significant magnetic fields ($\sim 10^{15}$ G) on a time scale that is relevant for the explosion mechanism of core-collapse supernovae. We describe our model, present some recent results, and discuss their potential relevance for supernova models.

Here we present the results from two sets of simulations, in two and three spatial dimensions. In two dimensions, the simulations include multifrequency flux-limited diffusion neutrino transport in the "ray-by-ray-plus" approximation, two-dimensional self gravity in the Newtonian limit, and nuclear burning through a 14-isotope alpha network. The three-dimensional simulations are model simulations constructed to reflect the post stellar core bounce conditions during neutrino shock reheating at the onset of explosion. They are hydrodynamics-only models that focus on critical aspects of the shock stability and dynamics and their impact on the supernova mechanism and explosion. In two dimensions, we obtain explosions (although in one case weak) for two progenitors (11 and 15 Solar mass models). Moreover, in both cases the explosion is initiated when the inner edge of the oxygen layer accretes through the shock. Thus, the shock is not revived while in the iron core, as previously discussed in the literature. The three-dimensional studies of the development of the stationary accretion shock instability (SASI) demonstrate the fundamentally new dynamics allowed when simulations are performed in three spatial dimensions. The predominant l=1 SASI mode gives way to a stable m=1 mode, which in turn has significant ramifications for the distribution of angular momentum in the region between the shock and proto-neutron star and, ultimately, for the spin of the remnant neutron star. Moreover, the three-dimensional simulations make clear, given the increased number of degrees of freedom, that two-dimensional models are severely limited by artificially imposed symmetries.

The overwhelming evidence that the core collapse supernova mechanism is inherently multidimensional, the complexity of the physical processes involved, and the increasing evidence from simulations that the explosion is marginal presents great computational challenges for the realistic modeling of this event, particularly in 3 spatial dimensions. We have developed a code which is scalable to computations in 3 dimensions which couples PPM Lagrangian with remap hydrodynamics [1], multigroup, flux-limited diffusion neutrino transport [2], with many improvements), and a nuclear network [3]. The neutrino transport is performed in a ray-by-ray plus approximation wherein all the lateral effects of neutrinos are included (e.g., pressure, velocity corrections, advection) except the transport. A moving radial grid option permits the evolution to be carried out from initial core collapse with only modest demands on the number of radial zones. The inner part of the core is evolved after collapse along with the rest of the core and mantle by subcycling the lateral evolution near the center as demanded by the small Courant times. We present results of 2-D simulations of a symmetric and an asymmetric collapse of both a 15 and an 11 M progenitor. In each of these simulations we have discovered that once the oxygen rich material reaches the shock there is a synergistic interplay between the reduced ram pressure, the energy released by the burning of the shock heated oxygen rich material, and the neutrino energy deposition which leads to a revival of the shock and an explosion.

First results from a fully self-consistent, temperature-dependent equation of state that spans the whole density range of neutron stars and supernova cores are presented. The equation of state (EoS) is calculated using a mean-field Hartree-Fock method in three dimensions (3D). The nuclear interaction is represented by the phenomenological Skyrme model in this work, but the EoS can be obtained in our framework for any suitable form of the nucleon-nucleon effective interaction. The scheme we employ naturally allows effects such as (i) neutron drip, which results in an external neutron gas, (ii) the variety of exotic nuclear shapes expected for extremely neutron heavy nuclei, and (iii) the subsequent dissolution of these nuclei into nuclear matter. In this way, the equation of state is calculated across phase transitions without recourse to interpolation techniques between density regimes described by different physical models. EoS tables are calculated in the wide range of densities, temperature and proton/neutron ratios on the ORNL NCCS XT3, using up to 2000 processors simultaneously.

Rotation-powered radio pulsars are born with inferred initial rotation periods of order 300 ms (some as short as 20 ms) in core-collapse supernovae. In the traditional picture, this fast rotation is the result of conservation of angular momentum during the collapse of a rotating stellar core. This leads to the inevitable conclusion that pulsar spin is directly correlated with the rotation of the progenitor star. So far, however, stellar theory has not been able to explain the distribution of pulsar spins, suggesting that the birth rotation is either too slow or too fast. Here we report a robust instability of the stalled accretion shock in core-collapse supernovae that is able to generate a strong rotational flow in the vicinity of the accreting proto-neutron star. Sufficient angular momentum is deposited on the proto-neutron star to generate a final spin period consistent with observations, even beginning with spherically symmetrical initial conditions. This provides a new mechanism for the generation of neutron star spin and weakens, if not breaks, the assumed correlation between the rotational periods of supernova progenitor cores and pulsar spin.

With presently known input physics and computer simulations in 1D, a self-consistent treatment of core collapse supernovae does not yet lead to successful explosions, while 2D models show some promise. Thus, there are strong indications that the delayed neutrino mechanism works combined with a multi-D convection treatment for unstable layers. On the other hand there is a need to provide correct nucleosynthesis abundances for the progressing field of galactic evolution and observations of low metallicity stars. The innermost ejecta is directly affected by the explosion mechanism, i.e. most strongly the yields of Fe-group nuclei for which an induced piston or thermal bomb treatment will not provide the correct yields because the effect of neutrino interactions is not included. We apply parameterized variations to the neutrino scattering cross sections and alternatively, parameterized variations are applied to the neutrino absorption cross sections on nucleons in the ``gain region''. We find that both measures lead to similar results, causing explosions and a Ye>0.5 in the innermost ejected layers, due to the combined effect of a short weak interaction time scale and a negligible electron degeneracy, unveiling the proton-neutron mass difference. We include all weak interactions (electron and positron capture, beta-decay, neutrino and antineutrino capture on nuclei, and neutrino and antineutrino capture on nucleons) and present first nucleosynthesis results for these innermost ejected layers to discuss how they improve predictions for Fe-group nuclei. The proton-rich environment results in enhanced abundances of 45Sc, 49Ti, and 64Zn as requested by chemical evolution studies and observations of low metallicity stars as well as appreciable production of nuclei in the mass range up to A=80.

The computational difficulty of six-dimensional neutrino radiation hydrodynamics has spawned a variety of approximations, provoking a long history of uncertainty in the core-collapse supernova explosion mechanism. Under the auspices of the Terascale Supernova Initiative, we are honoring the physical complexity of supernovae by meeting the computational challenge head-on, undertaking the development of a new adaptive mesh refinement code for self-gravitating, six-dimensional neutrino radiation magnetohydrodynamics. This code--called {\em GenASiS,} for {\em Gen}eral {\em A}strophysical {\em Si}mulation {\em S}ystem--is designed for modularity and extensibility of the physics. Presently in use or under development are capabilities for Newtonian self-gravity, Newtonian and special relativistic magnetohydrodynamics (with `realistic' equation of state), and special relativistic energy- and angle-dependent neutrino transport--including full treatment of the energy and angle dependence of scattering and pair interactions.

The computational difficulty of six-dimensional neutrino radiation hydrodynamics has spawned a variety of approximations, provoking a long history of uncertainty in the core-collapse supernova explosion mechanism. Under the auspices of the Terascale Supernova Initiative, we are honoring the physical complexity of supernovae by meeting the computational challenge head-on, undertaking the development of a new adaptive mesh refinement code for self-gravitating, six-dimensional neutrino radiation magnetohydrodynamics. This code--called {\em GenASiS}, for {\em Gen}eral {\em A}strophysical {\em Si}mulation {\em S}ystem--is designed for modularity and extensibility of the physics. Presently in use or under development are capabilities for Newtonian self-gravity, Newtonian and special relativistic magnetohydrodynamics (with `realistic' equation of state), and special relativistic energy- and angle-dependent neutrino transport--including full treatment of the energy and angle dependence of scattering and pair interactions.

Many astrophysical phenomena exhibit relativistic radiative flows. While velocities in excess of $v \sim 0.1c$ can occur in these systems, it has been common practice to approximate radiative transfer to $\cO(v/c)$. In the case of neutrino transport in core-collapse supernovae, this approximation gives rise to an inconsistency between the lepton number transfer and lab frame energy transfer, which have different $\cO(v/c)$ limits. A solution used in spherically symmetric $\cO(v/c)$ simulations has been to retain, for energy accounting purposes, the $\cO(v^2/c^2)$ terms in the lab frame energy transfer equation that arise from the $\cO(v/c)$ neutrino number transport equation. Avoiding the proliferation of such ``extra'' $\cO(v^2/c^2)$ terms in the absence of spherical symmetry motivates a special relativistic formalism, which we exhibit in coordinates sufficiently general to encompass Cartesian, spherical, and cylindrical coordinate systems.

We use time-dependent, axisymmetric, hydrodynamic simulations to study the linear stability of the stalled, spherical accretion shock that arises in the post-bounce phase of core-collapse supernovae. We show that this accretion shock is stable to radial modes, with decay rates and oscillation frequencies in close agreement with the linear stability analysis of Houck and Chevalier. For non-spherical perturbations we find that the l=1 mode is always unstable for parameters appropriate to core-collapse supernovae. We also find that the l=2 mode is unstable, but typically has a growth rate smaller than that for l=1. Furthermore, the l=1 mode is the only mode found to transition into a nonlinear stage in our simulations. This result provides a possible explanation for the dominance of an l=1 'sloshing' mode seen in many two-dimensional simulations of core-collapse supernovae.

In 2002, Ray Davis and Masatoshi Koshiba were awarded the Nobel Prize in Physics ``for pioneering contributions to astrophysics, in particular for the detection of cosmic neutrinos.'' However, while astronomy has undergone a revolution in understanding by synthesizing data taken at many wavelengths, the universe has only barely been glimpsed in neutrinos, just the Sun and the nearby SN 1987A. An entire universe awaits, and since neutrinos can probe astrophysical objects at densities, energies, and distances that are otherwise inaccessible, the results are expected to be particularly exciting. Similarly, the revolution in quantitative cosmology has heightened the need for very precise tests that depend on the effects of neutrinos, and prominent among them is the search for the effects of neutrino mass, since neutrinos are a small but known component of the dark matter. In this report, we highlight some of the key opportunties for progress in neutrino astrophysics and cosmology, and the implications for other areas of physics.

Accurate neutrino transport has been built into spherically symmetric simulations of stellar core collapse and postbounce evolution. The results of such simulations agree that spherically symmetric models with standard microphysical input fail to explode by the delayed, neutrino-driven mechanism. Independent groups implemented fundamentally different numerical methods to tackle the Boltzmann neutrino transport equation. Here we present a direct and detailed comparison of such neutrino radiation-hydrodynamical simulations for two codes, Agile-Boltztran of the Oak Ridge-Basel group and Vertex of the Garching group. The former solves the Boltzmann equation directly by an implicit, general relativistic discrete angle method on the adaptive grid of a conservative implicit hydrodynamics code with second-order TVD advection. In contrast, the latter couples a variable Eddington factor technique with an explicit, moving-grid, conservative high-order Riemann solver with important relativistic effects treated by an effective gravitational potential. The presented study is meant to test both neutrino radiation-hydrodynamics implementations and to provide a data basis for comparisons and verifications of supernova codes to be developed in the future. Results are discussed for simulations of the core collapse and post-bounce evolution of a 13 solar mass star with Newtonian gravity and a 15 solar mass star with relativistic gravity.

The pursuit of the core collapse supernova explosion mechanism continues. While such efforts have been undertaken over the last four decades, it is only in the last decade that multidimensional models have been developed, and only in the last few years that significant components of supernova models, such as the neutrino transport, have been modeled with sufficient realism. We can now identify what are arguably the important components that must be included in any realistic three-dimensional supernova model. The challenge we now face is to include all of these in the models.

We have developed a formalism for analyzing the stability of a fluid in the presence of neutrinos of all flavors and in the presence of a gravitational field. When applied to an extensive two-dimensional grid of core radii and fluid element sizes for each of several time slices of a number of proto-supernovae, we find no evidence for the neutron finger instability as described by the Livermore group. We find, instead, that the rate of lepton equilibration always exceeds that of thermal equilibration. Furthermore, we find that thermal diffusion driven by a lepton fraction difference and lepton diffusion driven by an entropy difference are both nonzero and that the first of these tends to be large in magnitude. An important consequence of this is the presence of a doubly diffusive instability, which we refer to as ``lepto-entropy fingers,'' in an extensive region below the neutrinosphere where the lepton number, $Y_{\ell}$, is small. This instability is driven by a mechanism very different from that giving rise to neutron fingers, and may play an important role in enhancing the neutrino emission. Deep in the core where the entropy is low and the lepton number higher, our analysis indicates a region unstable to another instability which we refer to as ``lepto-entropy semiconvection.'' These instabilities, particularly lepto-entropy fingers, may have already been seen in some multi-dimensional core collapse simulations described in the literature.

The most important weak nuclear interaction to the dynamics of stellar core collapse is electron capture, primarily on nuclei with masses larger than 60. In prior simulations of core collapse, electron capture on these nuclei has been treated in a highly parameterized fashion, if not ignored. With realistic treatment of electron capture on heavy nuclei come significant changes in the hydrodynamics of core collapse and bounce. We discuss these as well as the ramifications for the post-bounce evolution in core collapse supernovae.

New facilities to measure neutrino-nucleus cross sections, such as those possible in conjunction with spallation neutron sources, could provide an experimental foundation for the many neutrino-nucleus weak interaction rates needed in supernova models. This would enable more realistic supernova models and provide a greatly improved ability to understand the physics fundamental to supernovae by comparison of these models with detailed observations. Charged- and neutral-current neutrino interactions on nuclei in the stellar core play a central role in supernova dynamics and nucleosynthesis as well as being important for supernova neutrino detection. Measurements of these reactions on judiciously chosen targets would provide an invaluable test of the complex theoretical models used to compute the large number of neutrino-nucleus cross sections that are needed.

We present an implicit finite difference representation for general relativistic radiation hydrodynamics in spherical symmetry. Our code, Agile-Boltztran, solves the Boltzmann transport equation for the angular and spectral neutrino distribution functions in self-consistent simulations of stellar core collapse and postbounce evolution. It implements a dynamically adaptive grid in comoving coordinates. Most macroscopically interesting physical quantities are defined by expectation values of the distribution function. We optimize the finite differencing of the microscopic transport equation for a consistent evolution of important expectation values. We test our code in simulations launched from progenitor stars with 13 solar masses and 40 solar masses. ~0.5 s after core collapse and bounce, the protoneutron star in the latter case reaches its maximum mass and collapses further to form a black hole. When the hydrostatic gravitational contraction sets in, we find a transient increase in electron flavor neutrino luminosities due to a change in the accretion rate. The muon- and tauon-neutrino luminosities and rms energies, however, continue to rise because previously shock-heated material with a non-degenerate electron gas starts to replace the cool degenerate material at their production site. We demonstrate this by supplementing the concept of neutrinospheres with a more detailed statistical description of the origin of escaping neutrinos. We compare the evolution of the 13 solar mass progenitor star to simulations with the MGFLD approximation, based on a recently developed flux limiter. We find similar results in the postbounce phase and validate this MGFLD approach for the spherically symmetric case with standard input physics.

Proton capture by 17F plays an important role in the synthesis of nuclei in nova explosions. A revised rate for this reaction, based on a measurement of the 1H(17F,p)17F excitation function using a radioactive 17F beam at ORNL's Holifield Radioactive Ion Beam Facility, is used to calculate the nucleosynthesis in nova outbursts on the surfaces of 1.25 and 1.35 solar mass ONeMg white dwarfs and a 1.00 solar mass CO white dwarf. We find that the new 17F(p,gamma)18Ne reaction rate changes the abundances of some nuclides (e.g., 17O) synthesized in the hottest zones of an explosion on a 1.35 solar mass white dwarf by more than a factor of 10,000 compared to calculations using some previous estimates for this reaction rate, and by more than a factor of 3 when the entire exploding envelope is considered. In a 1.25 solar mass white dwarf nova explosion, this new rate changes the abundances of some nuclides synthesized in the hottest zones by more than a factor of 600, and by more than a factor of 2 when the entire exploding envelope is considered. Calculations for the 1.00 solar mass white dwarf nova show that this new rate changes the abundance of 18Ne by 21%, but has negligible effect on all other nuclides. Comparison of model predictions with observations is also discussed.

Supernova simulations to date have assumed that during core collapse electron captures occur dominantly on free protons, while captures on heavy nuclei are Pauli-blocked and are ignored. We have calculated rates for electron capture on nuclei with mass numbers A=65-112 for the temperatures and densities appropriate for core collapse. We find that these rates are large enough so that, in contrast to previous assumptions, electron capture on nuclei dominates over capture on free protons. This leads to significant changes in core collapse simulations.

Experience with core-collapse supernova simulations shows that accurate accounting of total particle number and 4-momentum can be a challenge for computational radiative transfer. This accurate accounting would be facilitated by the use of particle number and 4-momentum transport equations that allow transparent conversion between volume and surface integrals in both configuration and momentum space. Such conservative formulations of general relativistic kinetic theory in multiple spatial dimensions are presented in this paper, and their relevance to core-collapse supernova simulations is described.

General relativistic multi-group and multi-flavor Boltzmann neutrino transport in spherical symmetry adds a new level of detail to the numerical bridge between microscopic nuclear and weak interaction physics and the macroscopic evolution of the astrophysical object. Although no supernova explosions are obtained, we investigate the neutrino luminosities in various phases of the postbounce evolution for a wide range of progenitor stars between 13 and 40 solar masses. The signal probes the dynamics of material layered in and around the protoneutron star and is, within narrow limits, sensitive to improvements in the weak interaction physics. Only changes that dramatically exceed physical limitations allow experiments with exploding models. We discuss the differences in the neutrino signal and find the electron fraction in the innermost ejecta to exceed 0.5 as a consequence of thermal balance and weak equilibrium at the masscut.

We examine the stability of standing, spherical accretion shocks. Accretion shocks arise in core collapse supernovae (the focus of this paper), star formation, and accreting white dwarfs and neutron stars. We present a simple analytic model and use time-dependent hydrodynamics simulations to show that this solution is stable to radial perturbations. In two dimensions we show that small perturbations to a spherical shock front can lead to rapid growth of turbulence behind the shock, driven by the injection of vorticity from the now non-spherical shock. We discuss the ramifications this instability may have for the supernova mechanism.

Based on general relativistic three-flavor Boltzmann neutrino transport, we compare simulations of the postbounce evolution for five different stellar progenitors with 13 and 20 solar masses (Nomoto & Hashimoto '88), and 15, 25, and 40 solar masses (Woosley & Weaver '95).

With the exception of the Big Bang, responsible for 1,2H, 3,4He, and 7Li, stars act as sources for the composition of the interstellar medium. Cosmic rays are related to the latter and very probably due to acceleration of the mixed interstellar medium by shock waves from supernova remnants. Thus, the understanding of the abundance evolution in the interstellar medium and especially the enrichment of heavy elements, as a function of space and time, is essential. It reflects the history of star formation and the lifetimes of the diverse contributing stellar objects. Therefore, the understanding of the endpoints of stellar evolution is essential as well. These are mainly planetary nebulae and type II/Ib/Ic supernovae as evolutionary endpoints of single stars, but also events in binary systems can contribute, like e.g. supernovae of type Ia, novae and possibly X-ray bursts and neutron star or neutron star - black hole mergers. Despite many efforts, a full and self-consistent understanding of supernovae (the main contributors to nucleosynthesis in galaxies) is not existing, yet. Their fingerprints, however, seen either in spectra, lightcurves, radioactivities/decay gamma-rays or in galactic evolution, can help to constrain the composition of their ejecta and related model uncertainties.

We propose an approximation to full relativity that captures the main gravitational effects of dynamical importance in supernovae. The conceptual link between this formalism and the Newtonian limit is such that it could likely be implemented relatively easily in existing multidimensional Newtonian gravitational hydrodynamics codes employing a Poisson solver. As a test of the formalism's utility, we display results for rapidly rotating (and therefore highly deformed) neutron stars.

We propose an approximation to general relativity that captures the main gravitational effects of dynamical importance in supernovae. The conceptual link between this formalism and the Newtonian limit is such that it could likely be implemented in existing multidimensional Newtonian gravitational hydrodynamics codes employing a Poisson solver. As a test of the formalism's utility, we display results for rapidly rotating (and therefore highly deformed) neutron stars.

We report on the stellar core collapse, bounce, and postbounce evolution of a 13 solar mass star in a self-consistent general relativistic spherically symmetric simulation based on Boltzmann neutrino transport. We conclude that approximations to exact neutrino transport and omission of general relativistic effects were not alone responsible for the failure of numerous preceding attempts to model supernova explosions in spherical symmetry. Compared to simulations in Newtonian gravity, the general relativistic simulation results in a smaller shock radius. We however argue that the higher neutrino luminosities and rms energies in the general relativistic case could lead to a larger supernova explosion energy.

The description of general relativistic radiation hydrodynamics in spherical symmetry is presented in natural coordinate choices. For hydrodynamics, comoving coordinates are chosen, and the momentum phase space for the radiation particles is described in comoving frame four momenta. We also investigate a description of the momentum phase space in terms of particle impact parameter and energy at infinity and derive a simple approximation to the general relativistic Boltzmann equation. Further developed are, however, the exact equations in comoving coordinates, because the description of the interaction between matter and radiation particles is best described in the closely related orthonormal basis comoving with the fluid elements. We achieve a conservative and concise formulation of radiation hydrodynamics that is well suited for numerical implementation by a variety of methods. The contribution of radiation to the general relativistic jump conditions at shock fronts is discussed, and artificial viscosity is consistently included in the derivations in order to support approaches relying on this option.

We present self-consistent general relativistic simulations of stellar core collapse, bounce, and postbounce evolution for 13, 15, and 20 solar mass progenitors in spherical symmetry. Our simulations implement three-flavor Boltzmann neutrino transport and standard nuclear physics. The results are compared to our corresponding simulations with Newtonian hydrodynamics and O(v/c) Boltzmann transport.

Except for 1H, 2H, 3He, 4He, and 7Li, originating from the Big Bang, all heavier elements are made in stellar evolution and stellar explosions. Nuclear physics, and in many cases nuclear structure far from stability, enters in a crucial way. Therefore, we examine in this review the role of nuclear physics in astrophysics in general and in particular how it affects stellar events and the resulting nucleosynthesis. Stellar modeling addresses four major aspects: 1. energy generation and nucleosynthesis, 2. energy transport via conduction, radiation or possibly convection, 3. hydrodynamics/hydrostatics, and finally 4. thermodynamic properties of the matter involved. Nuclear Physics enters via nuclear reaction cross sections and nuclear structure (affecting the composition changes and nuclear energy generation), neutrino-nucleon and neutrino-nucleus cross sections (affecting neutrino opacities and transport), and e.g. the equation of state at and beyond nuclear densities which creates a relation between the nuclear many-body problem and the hydrodynamic response like pressure and entropy. In the following we review these four topics by highlighting the role and impact of nuclear physics in each of these aspects of stellar modeling. The main emphasis is put on the connection to element synthesis.

We apply our recently developed code for spherically symmetric, fully general relativistic (GR) Lagrangian hydrodynamics and multigroup flux-limited diffusion neutrino transport to examine the effects of GR on the hydrodynamics and transport during collapse, bounce, and the critical shock reheating phase of core collapse supernovae. Comparisons of models computed with GR versus Newtonian hydrodynamics show that collapse to bounce takes slightly less time in the GR limit, and that the shock propagates slightly farther out in radius before receding. After a secondary quasistatic rise in the shock radius, the shock radius declines considerably more rapidly in the GR simulations than in the corresponding Newtonian simulations. During the shock reheating phase, core collapse computed with GR hydrodynamics results in a substantially more compact structure from the center out to the stagnated shock. The inflow speed of material behind the shock is also increased. Comparisons also show that the luminosity and rms energy of any neutrino flavor during the shock reheating phase increases when switching from Newtonian to GR hydrodynamics, and decreases when switching from Newtonian to GR transport. This latter decrease in neutrino luminosities and rms energies is less in magnitude than the increase that arise when switching from Newtonian to GR hydrodynamics, with the result that a fully GR simulation gives higher neutrino luminosities and harder neutrino spectra than a fully Newtonian simulation of the same precollapse model.

Early black hole formation in core-collapse supernovae may be triggered by mass accretion or a change in the high-density equation of state. We consider the possibility that black hole formation happens when the flux of neutrinos is still measurably high. If this occurs, then the neutrino signal from the supernova will be terminated abruptly (the transition takes $\lesssim 0.5$ ms). The properties and duration of the signal before the cutoff are important measures of both the physics and astrophysics of the cooling proto-neutron star. For the event rates expected in present and proposed detectors, the cutoff will generally appear sharp, thus allowing model-independent time-of-flight mass tests for the neutrinos after the cutoff. If black hole formation occurs relatively early, within a few ($\sim 1$) seconds after core collapse, then the expected luminosities are of order $L_{BH} = 10^{52}$ erg/s per flavor. In this case, the neutrino mass sensitivity can be extraordinary. For a supernova at a distance $D = 10$ kpc, SuperKamiokande can detect a $\bar{\nu}_e$ mass down to 1.8 eV by comparing the arrival times of the high-energy and low-energy neutrinos in $\bar{\nu}_e + p \to e^+ + n$. This test will also measure the cutoff time, and will thus allow a mass test of $\nu_\mu$ and $\nu_\tau$ relative to $\bar{\nu}_e$. Assuming that $\nu_\mu$ and $\nu_\tau$ are nearly degenerate, as suggested by the atmospheric neutrino results, masses down to about 6 eV can be probed with a proposed lead detector of mass $M_D = 4$ kton (OMNIS). Remarkably, the neutrino mass sensitivity scales as $(D/L_{BH} M_D)^{1/2}$. Therefore, {\it direct} sensitivity to all three neutrino masses in the interesting few-eV range is realistically possible; {\it there are no other known techniques that have this capability}.

The nucleosynthesis and other observable consequences of a nova outburst depend sensitively on the details of the thermonuclear runaway which initiates the outburst. One important source of uncertainty in our current models is the nuclear reaction data used as input for the evolutionary calculations. We present preliminary results of the first analyses of the impact on nova nucleosynthesis of all reaction rate uncertainties considered simultaneously.

With exact three-flavor Boltzmann neutrino transport, we simulate the stellar core collapse, bounce, and postbounce evolution of a 13 solar mass star in spherical symmetry, the Newtonian limit, without invoking convection. In the absence of convection, prior spherically symmetric models, which implemented approximations to Boltzmann transport, failed to produce explosions. We are motivated to consider exact transport to determine if these failures were due to the transport approximations made and to answer remaining fundamental questions in supernova theory. The model presented here is the first in a sequence of models beginning with different progenitors. In this model, a supernova explosion is not obtained. We discuss the ramifications of our results for the supernova mechanism.

In this paper, we discuss the current status of core collapse supernova models and the future developments needed to achieve significant advances in understanding the supernova mechanism and supernova phenomenology, i.e., in developing a supernova standard model.

A new facility to measure neutrino mass differences and mixing angles and neutrino-nucleus cross sections, such as the proposed ORLaND facility at Oak Ridge, would contribute to the experimental determination of vacuum mixing parameters and would provide an experimental foundation for the many neutrino-nucleus weak interaction rates needed in supernova models. This would enable more realistic supernova models and a far greater ability to cull fundamental physics from these models by comparing them with detailed observations. Charged- and neutral-current neutrino interactions on nuclei in the stellar core play a central role in supernova dynamics, nucleosynthesis, and neutrino detection. Measurements of these reactions on select, judiciously chosen targets would provide an invaluable test of the complex theoretical models used to compute the neutrino-nucleus cross sections.

In this paper, we present results from a simulation of stellar core collapse, bounce, and postbounce evolution with Boltzmann neutrino transport. We motivate the development of our Boltzmann solver in light of the sensitivity of the neutrino-heating core collapse supernova paradigm to details in the neutrino transport, particularly near the neutrinospheres, where the neutrinos are neither diffusing nor free streaming and a kinetic description is necessary, and in light of the mixed outcomes and transport approximations used in all prior supernova models in both one and two dimensions. We discuss the implications of our findings for the supernova mechanism and future supernova research. We also present the results of a Boltzmann transport prediction of the early neutrino light curves in the model included here.

Early black hole formation in a core-collapse supernova will abruptly truncate the neutrino fluxes. The sharp cutoff can be used to make model-independent time-of-flight neutrino mass tests. Assuming a neutrino luminosity of $10^{52}$ erg/s per flavor at cutoff and a distance of 10 kpc, SuperKamiokande can detect an electron neutrino mass as small as 1.8 eV, and the proposed OMNIS detector can detect mu and tau neutrino masses as small as 6 eV. This {\it Letter} presents the first technique with direct sensitivity to eV-scale mu and tau neutrino masses.

The observable consequences of a nova outburst depend sensitively on the details of the thermonuclear runaway which initiates the outburst. One of the more important sources of uncertainty is the nuclear reaction data used as input for the evolutionary calculations. A recent paper by Starrfield, Truran, Wiescher, & Sparks (1998) has demonstrated that changes in the reaction rate library used within a nova simulation have significant effects, not just on the production of individual isotopes (which can change by an order of magnitude), but on global observables such as the peak luminosity and the amount of mass ejected. We present preliminary results of systematic analyses of the impact of reaction rate uncertainties on nova nucleosynthesis.

The effects of general relativity (GR) on the hydrodynamics and neutrino transport are examined during the critical shock reheating phase of core collapse supernovae. We find that core collapse computed with GR hydrodynamics results in a substantially more compact core structure out to the shock, the shock radius at stagnation being reduced by a factor of 2. The inflow speed of material behind the shock is also increased by a factor of 2 throughout most of the evolution. We have developed a code for general relativistic multigroup flux-limited diffusion (MGFLD) in static spacetimes and compared the steady-state neutrino distributions for selected time slices of post-bounce models with those computed with Newtonian MGFLD. The GR transport calculations show the expected reductions in neutrino luminosities and rms energies from redshift and curvature effects. Although the effects of GR on the hydrodynamics and neutrino transport seem to work against shock revival, the core configurations are sufficiently different that no firm conclusions can be drawn, except that simulations of core collapse supernovae using Newtonian hydrodynamics and transport are not realistic.

We compare Newtonian three-flavor multigroup Boltzmann (MGBT) and (Bruenn's) multigroup flux-limited diffusion (MGFLD) neutrino transport in postbounce core collapse supernova environments. We focus our study on quantities central to the postbounce neutrino heating mechanism for reviving the stalled shock. Stationary-state three-flavor neutrino distributions are developed in thermally and hydrodynamically frozen time slices obtained from core collapse and bounce simulations that implement Lagrangian hydrodynamics and MGFLD neutrino transport. Most important, we find, for a region above the gain radius, net heating rates for MGBT that are as much as ~2 times the corresponding MGFLD rates, and net cooling rates below the gain radius that are typically ~0.8 times the MGFLD rates. These differences stem from differences in the neutrino luminosities and mean inverse flux factors, which can be as much as 11% and 24%, respectively. They are greatest at earlier postbounce times for a given progenitor mass and, for a given postbounce time, greater for greater progenitor mass. We discuss the ramifications these new results have for the supernova mechanism.

Core-collapse supernovae are, despite their spectacular visual display, neutrino events. Virtually all of the 10^53 ergs of gravitational binding energy released in the formation of the nascent neutron star is carried away in the form of neutrinos and antineutrinos of all three flavors, and these neutrinos are primarily responsible for powering the explosion. This mechanism depends sensitively on the neutrino transport between the neutrinospheres and the shock. In light of this, we have performed a comparison of multigroup Boltzmann neutrino transport (MGBT) and multigroup flux-limited diffusion (MGFLD) in post-core bounce environments. Differences in the mean inverse flux factors, luminosities, and RMS energies translate to heating rates that are up to 2 times larger for Boltzmann transport, with net cooling rates below the gain radius that are typically 0.8 times the MGFLD rates. These differences are greatest at earlier postbounce times for a given progenitor mass, and for a given postbounce time, greater for greater progenitor mass. The increased differences with increased progenitor mass suggest that the net heating enhancement from MGBT is potentially robust and self-regulated.

We investigate neutrino-driven convection in core collapse supernovae and its ramifications for the explosion mechanism, for a 15 solar mass model. Our two-dimensional simulation begins at 12 ms after bounce and proceeds for 500 ms. We couple two-dimensional PPM hydrodynamics to precalculated one-dimensional MGFLD neutrino transport. (The accuracy of this approximation is assessed.) For the moment we sacrifice dimensionality for realism in other aspects of our neutrino transport. MGFLD is an implementation of neutrino transport that simultaneously (a) is multigroup and (b) simulates with greater realism the transport of neutrinos in opaque, semitransparent, and transparent regions. Both are crucial to the accurate determination of postshock neutrino heating, which sensitively depends on the luminosities, spectra, and flux factors of the electron neutrinos and antineutrinos emerging from their respective neutrinospheres. By 212 ms after bounce, we see large-scale, neutrino-driven convection beneath the shock, characterized by higher-entropy, expanding upflows and denser, lower-entropy, finger-like downflows. The radial convection velocities at this time become supersonic just below the shock, reaching magnitudes in excess of 10^{9} cm/sec. Eventually, however, the shock recedes to smaller radii, and at 500 ms after bounce there is no evidence in our simulation of an explosion. While vigorous, neutrino-driven convection in our model does not have a significant impact on the overall shock dynamics. The differences between our results and those of other groups are considered. These most likely result from differences in (1) numerical hydrodynamics methods, (2) initial postbounce models, and most important, (3) neutrino transport approximations.

We couple two-dimensional hydrodynamics to realistic one-dimensional multigroup flux-limited diffusion neutrino transport to investigate protoneutron star convection in core collapse supernovae, and more specifically, the interplay between its development and neutrino transport, for both 15 and 25 solar mass models. In the presence of neutrino transport, protoneutron star convection velocities are too small relative to bulk inflow velocities to result in any significant convective transport of entropy and leptons. A simple analytical model supports our numerical results, indicating that the inclusion of neutrino transport reduces the entropy-driven (lepton-driven) convection growth rates and asymptotic velocities by a factor of 3 (50) at the neutrinosphere and a factor 250 (1000) at a density of 10^{12} g/cm^{3}, for both our 15 and 25 solar mass models. Moreover, when transport is included, the initial postbounce entropy gradient is smoothed out by neutrino diffusion, whereas the initial lepton gradient is maintained by electron capture and neutrino escape near the neutrinosphere. Despite the maintenance of the lepton gradient, protoneutron star convection does not develop over the 100 ms duration typical of all our simulations, except in the instance where ``low-test'' initial conditions are used, which are generated by core collapse and bounce simulations that neglect neutrino-electron scattering and ion-ion screening corrections to neutrino- nucleus elastic scattering.

In one dimension, we present results from comparisons of stationary state multigroup flux-limited diffusion and Boltzmann neutrino transport, focusing on quantities central to the postbounce shock reheating. In two dimensions, we present results from simulations that couple one-dimensional multigroup flux-limited diffusion to two-dimensional PPM hydrodynamics.

We investigate neutrino-driven convection in core collapse supernovae and its ramifications for the explosion mechanism. We begin with an ``optimistic'' 15 solar mass precollapse model, which is representative of the class of stars with compact iron cores. This model is evolved through core collapse and bounce in one dimension using multigroup (neutrino-energy--dependent) flux-limited diffusion (MGFLD) neutrino transport and Lagrangian hydrodynamics, providing realistic initial conditions for the postbounce convection and evolution. Our two-dimensional simulation begins at 106 ms after bounce at a time when there is a well-developed gain region, and proceeds for 400 ms. We couple two-dimensional (PPM) hydrodynamics to one-dimensional MGFLD neutrino transport. At 225 ms after bounce we see large-scale convection behind the shock, characterized by high-entropy, mushroom-like, expanding upflows and dense, low-entropy, finger-like downflows. The upflows reach the shock and distort it from sphericity. The radial convection velocities become supersonic just below the shock, reaching magnitudes in excess of 10^9 cm/sec. Eventually, however, the shock recedes to smaller radii, and at about 500 ms after bounce there is no evidence in our simulation of an explosion or of a developing explosion. Failure in our ``optimistic'' 15 solar mass Newtonian model leads us to conclude that it is unlikely, at least in our approximation, that neutrino-driven convection will lead to explosions for more massive stars with fatter iron cores or in cases in which general relativity is included.

We couple two-dimensional hydrodynamics to detailed one-dimensional multigroup flux-limited diffusion neutrino transport to investigate prompt convection in core collapse supernovae. Our initial conditions, time-dependent boundary conditions, and neutrino distributions for computing neutrino heating, cooling, and deleptonization rates are obtained from one-dimensional simulations that implement multigroup flux-limited diffusion neutrino transport and one-dimensional hydrodynamics. The development and evolution of prompt convection and its ramifications for the shock dynamics are investigated for both 15 and 25 solar mass models, representative of the two classes of stars with compact and extended iron cores, respectively. In the absence of neutrino transport, prompt convection develops and dissipates on a time scale $\sim$15 ms for both models. Prompt convection seeds convection behind the shock, which causes distortions in the shock's sphericity, but on the average, the shock radius is not boosted significantly. In the presence of neutrino transport, prompt convection velocities are too small relative to bulk inflow velocities to result in any significant convective transport of entropy and leptons. A simple analytical model supports our numerical results, indicating that the inclusion of transport reduces the convection growth rates and asymptotic velocities by factors of 4--250.

Modeling the atmospheres of SNe~Ia requires the solution of the NLTE radiative transfer equation. We discuss the formulation of the radiative transfer equation in the co-moving frame. For characteristic velocities larger than about 2000 km/s, the effects of advection on the synthetic spectra are non-negligible, and hence should be included in model calculations. We show that the time-independent or quasi-static approximation is adequate for SNe~Ia near maximum light, as well as for most other astrophysical problems; e.g., hot stars, novae, and other types of supernovae. We examine the use of the Sobolev approximation in modeling moving atmospheres and find that the number of overlapping lines in the co-moving frame make the approximation suspect in models that predict both lines and continua. We briefly discuss the form of the Rosseland mean opacity in the co-moving frame, and present a formula that is easy to implement in radiation hydrodynamics calculations.