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.

Vela X is the prototypical example of a pulsar wind nebula whose morphology and detailed structure have been affected by the interaction with the reverse shock of its host supernova remnant. The resulting complex of filamentary structure and mixed-in ejecta embedded in a nebula that is offset from the pulsar provides the best example we have of this middle-age state that characterizes a significant fraction of composite SNRs, and perhaps all of the large-diameter PWNe seen as TeV sources. Here we report on an XMM-Newton Large Project study of Vela X, supplemented by additional Chandra observations. Through broad spectral modeling as well as detailed spectral investigations of discrete emission regions, we confirm previous studies that report evidence for ejecta material within Vela X, and show that equivalent width variations of O VII and O VIII are consistent with temperature maps within the PWN that show low-temperature regions where the projected SNR emission appears to dominate emission from the ejecta. We identify spectral variations in the nonthermal emission, with hard emission being concentrated near the pulsar. We carry out investigations of the Vela X "cocoon" structure and, with hydrodynamical simulations, show that its overall properties are consistent with structures formed in the late-phase evolution of a composite SNR expanding into a surrounding medium with a density gradient, with ejecta material being swept beyond the pulsar and compressed into an elongated structure in the direction opposite the high external density.

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.

For pulsars born in supernovae, the expansion of the shocked pulsar wind nebula is initially in the freely expanding ejecta of the supernova. While the nebula is in the inner flat part of the ejecta density profile, the swept-up, accelerating shell is subject to the Rayleigh-Taylor instability. We carried out 2 and 3-dimensional simulations showing that the instability gives rise to filamentary structure during this initial phase but does not greatly change the dynamics of the expanding shell. The flow is effectively self-similar. If the shell is powered into the outer steep part of the density profile, the shell is subject to a robust Rayleigh-Taylor instability in which the shell is fragmented and the shocked pulsar wind breaks out through the shell. The flow is not self-similar in this phase. For a wind nebula to reach this phase requires that the deposited pulsar energy be greater than the supernova energy, or that the initial pulsar period be in the ms range for a typical $10^{51}$ erg supernova. These conditions are satisfied by some magnetar models for Type I superluminous supernovae. We also consider the Crab Nebula, which may be associated with a low energy supernova for which this scenario applies.

The interaction between a pulsar wind nebula (PWN) and its host supernova remnant (SNR) can produce a vast array of observable structures. Asymmetry present within these structures derives from the complexity of the composite system, where many factors take turns playing a dominating hand throughout the stages of composite SNR evolution. Of particular interest are systems characterized by blastwave expansion within a nonuniform interstellar medium (ISM), which contain an active pulsar having a substantial "kick" velocity (upward of 300 km s$^{-1}$ ), because these systems tend to produce complex morphologies. We present a numerical model that employs these and several other factors in an effort to generate asymmetry similar to that seen in various X-ray and radio observations. We find that the main parameters driving structure are ISM uniformity and total pulsar spin- down energy, with secondary contributions from factors such as pulsar trajectory and initial spin-down luminosity. We also investigate the dynamics behind PWN "tails," which may form to link active pulsars to a crushed, relic nebula as the reverse shock passes. We find that the directions of such tails are not good indicators of pulsar motion, but direction does reveal the flow of ejecta created by the passage of a reverse shock.

We present the first three-dimensional measurements of the velocity of various ejecta knots in Tycho's supernova remnant, known to result from a Type Ia explosion. Chandra X-ray observations over a 12-year baseline from 2003 to 2015 allow us to measure the proper motion of nearly 60 "tufts" of Si-rich ejecta, giving us the velocity in the plane of the sky. For the line of sight velocity, we use two different methods: a non-equilibrium ionization model fit to the strong Si and S lines in the 1.2-2.8 keV regime, and a fit consisting of a series of Gaussian lines. These methods give consistent results, allowing us to determine the red or blue shift of each of the knots. Assuming a distance of 3.5 kpc, we find total velocities that range from 2400 to 6600 km s$^{-1}$, with a mean of 4430 km s$^{-1}$. We find several regions where the ejecta knots have overtaken the forward shock. These regions have proper motions in excess of 6000 km s$^{-1}$. Some Type Ia supernova explosion models predict a velocity asymmetry in the ejecta. We find no such velocity asymmetries in Tycho, and discuss our findings in light of various explosion models, favoring those delayed detonation models with relatively vigorous and symmetrical deflagrations. Finally, we compare measurements with models of the remnant's evolution that include both smooth and clumpy ejecta profiles, finding that both ejecta profiles can be accommodated by the observations.

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.

We present newly obtained X-ray and radio observations of Tycho's supernova remnant using {\it Chandra} and the Karl G. Jansky Very Large Array in 2015 and 2013/14, respectively. When combined with earlier epoch observations by these instruments, we now have time baselines for expansion measurements of the remnant of 12-15 year in the X-rays and 30 year in the radio. The remnant's large angular size allows for proper motion measurements at many locations around the periphery of the blast wave. We find, consistent with earlier measurements, a clear gradient in the expansion velocity of the remnant, despite its round shape. The proper motions on the western and southwestern sides of the remnant are about a factor of two higher than those in the east and northeast. We showed in an earlier work that this is related to an offset of the explosion site from the geometric center of the remnant due to a density gradient in the ISM, and using our refined measurements reported here, we find that this offset is $\sim 23"$ towards the northeast. An explosion center offset in such a circular remnant has implications for searches for progenitor companions in other remnants.

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.

X-ray astronomy allows study of objects which may be associated with compact objects, i.e. neutron stars or black holes, and also may contain strong magnetic fields. Such objects are categorically non-spherical, and likely non-circular when projected on the sky. Polarization allows study of such geometric effects, and X-ray polarimetry is likely to become feasible for a significant number of sources in the future. Potential targets for future X-ray polarization observations are the high-mass X-ray binaries (HMXBs), which consist of a compact object in orbit with an early type star. In this paper ws show that X-ray polarization from HMXBs has a distinct signature which depends on the source inclination and orbital phase. The presence of the X-ray source displaced from the star creates linear polarization even if the primary wind is spherically symmetric whenever the system is viewed away from conjunction. Direct X-rays dilute this polarization whenever the X-ray source is not eclipsed; at mid-eclipse the net polarization is expected to be small or zero if the wind is circularly symmetric around the line of centers. Resonance line scattering increases the scattering fraction, often by large factors, over the energy band spanned by resonance lines. Real winds are not expected to be spherically symmetric, or circularly symmetric around the line of centers, owing to the combined effects of the compact object gravity and ionization on the wind hydrodynamics. A sample calculation shows that this creates polarization fractions ranging up to tens of percent at mid-eclipse.

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 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.

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.

In an effort to better understand the evolution of composite supernova remnants (SNRs) and the eventual fate of relativistic particles injected by their pulsars, we present a multifaceted investigation of the interaction between a pulsar wind nebula (PWN) and its host SNR G327.1-1.1. Our 350 ks Chandra X-ray observations of SNR G327.1-1.1 reveal a highly complex morphology; a cometary structure resembling a bow shock, prong-like features extending into large arcs in the SNR interior, and thermal emission from the SNR shell. Spectral analysis of the non-thermal emission offers clues about the origin of the PWN structures, while enhanced abundances in the PWN region provide evidence for mixing of supernova ejecta with PWN material. The overall morphology and spectral properties of the SNR suggest that the PWN has undergone an asymmetric interaction with the SNR reverse shock (RS) that can occur as a result of a density gradient in the ambient medium and/or a moving pulsar that displaces the PWN from the center of the remnant. We present hydrodynamical simulations of G327.1-1.1 that show that its morphology and evolution can be described by a ~ 17,000 yr old composite SNR that expanded into a density gradient with an orientation perpendicular to the pulsar's motion. We also show that the RS/PWN interaction scenario can reproduce the broadband spectrum of the PWN from radio to gamma-ray wavelengths. The analysis and modeling presented in this work have important implications for our general understanding of the structure and evolution of composite SNRs.

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 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.

Supergiant High Mass X-ray Binary systems (sgHMXBs) consist of a massive, late type, star and a neutron star. The massive stars exhibits strong, radiatively driven, stellar winds. Wind accretion onto compact object triggers X-ray emission, which alters the stellar wind significantly. Hydrodynamic simulation has been used to study the neutron star - stellar wind interaction it two sgHMXBs: i) A heavily obscured sgHMXB (IGR J17252-3616) discovered by INTEGRAL. To account for observable quantities (i.e., absorbing column density) we have to assume a very slow wind terminal velocity of about 500 km/s and a rather massive neutron star. If confirmed in other obscured systems, this could provide a completely new stellar wind diagnostics. ii) A classical sgHMXB (Vela X-1) has been studied in depth to understand the origin of the off-states observed in this system. Among many models used to account for this observed behavior (clumpy wind, gating mechanism) we propose that self-organized criticality of the accretion stream is the likely reason for the observed behavior. In conclusion, the neutron star, in these two examples, acts very efficiently as a probe to study stellar winds.

{\it Spitzer} images of Tycho's supernova remnant in the mid-infrared reveal limb-brightened emission from the entire periphery of the shell and faint filamentary structures in the interior. As with other young remnants, this emission is produced by dust grains, warmed to $\sim 100$ K in the post-shock environment by collisions with energetic electrons and ions. The ratio of the 70 to 24 $\mu$m fluxes is a diagnostic of the dust temperature, which in turn is a sensitive function of the plasma density. We find significant variations in the 70/24 flux ratio around the periphery of Tycho's forward shock, implying order-of-magnitude variations in density. While some of these are likely localized interactions with dense clumps of the interstellar medium, we find an overall gradient in the ambient density surrounding Tycho, with densities 3-10 times higher in the NE than in the SW. This large density gradient is qualitatively consistent with the variations in the proper motion of the shock observed in radio and X-ray studies. Overall, the mean ISM density around Tycho is quite low ($\sim 0.1-0.2$ cm$^{-3}$), consistent with the lack of thermal X-ray emission observed at the forward shock. We perform two-dimensional hydrodynamic simulations of a Type Ia SN expanding into a density gradient in the ISM, and find that the overall round shape of the remnant is still easily acheivable, even for explosions into significant gradients. However, this leads to an offset of the center of the explosion from the geometric center of the remnant of up to 20%, although lower values of 10% are preferred. The best match with hydrodynamical simulations is achieved if Tycho is located at a large (3-4 kpc) distance in a medium with a mean preshock density of $\sim 0.2$ cm$^{-3}$. Such preshock densities are obtained for highly ($\ga 50$%) porous ISM grains.

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.

Kepler's supernova remnant resulted from a thermonuclear explosion, but is interacting with circumstellar material (CSM) lost from the progenitor system. We describe a statistical technique for isolating X-ray emission due to CSM from that due to shocked ejecta. Shocked CSM coincides well in position with 24 $\mu$m emission seen by {\sl Spitzer}. We find most CSM to be distributed along the bright north rim, but substantial concentrations are also found projected against the center of the remnant, roughly along a diameter with position angle $\sim 100^\circ$. We interpret this as evidence for a disk distribution of CSM before the SN, with the line of sight to the observer roughly in the disk plane. We present 2-D hydrodynamic simulations of this scenario, in qualitative agreement with the observed CSM morphology. Our observations require Kepler to have originated in a close binary system with an AGB star companion.

We explore the morphology of Type Ia supernova remnants (SNRs) using three-dimensional hydrodynamics modeling and an exponential density profile. Our model distinguishes ejecta from the interstellar medium (ISM), and tracks the ionization age of shocked ejecta, both of which allow for additional analysis of the simulated remnants. We also include the adiabatic index as a free parameter, which affects the compressibility of the fluid and emulates the efficiency of cosmic ray acceleration by shock fronts. In addition to generating 3-D images of the simulations, we compute line-of-sight projections through the remnants for comparison against observations of Tycho's SNR and SN 1006. We find that several features observed in these two remnants, such as the separation between the fluid discontinuities and the presence of ejecta knots ahead of the forward shock, can be generated by smooth ejecta without any initial clumpiness. Our results are consistent with SN 1006 being dynamically younger than Tycho's SNR, and more efficiently accelerating cosmic rays at its forward shock. We conclude that clumpiness is not a necessary condition to reproduce many observed features of Type Ia supernova remnants, particularly the radial profiles and the fleecy emission from ejecta at the central region of both remnants.

A population of obscured supergiant High Mass X-ray Binaries (sgHMXBs) has been discovered by INTEGRAL. X-ray wind tomography of IGR J17252-3616 inferred a slow wind velocity to account for the enhanced obscuration. The main goal of this study is to understand under which conditions high obscuration could occur. We have used an hydrodynamical code to simulate the flow of the stellar wind around the neutron star. A grid of simulations was used to study the dependency of the absorbing column density and of the X-ray light-curves on the model parameters. A comparison between the simulation results and the observations of IGR J17252-3616 provides an estimate on these parameters. We have constrained the wind terminal velocity to 500-600 km/s and the neutron star mass to 1.75-2.15 solar masses. We have confirmed that the initial hypothesis of a slow wind velocity with a moderate mass loss rate is valid. The mass of the neutron star can be constrained by studying its impact on the accretion flow.

We investigate the stability of gravitational accretion of an ideal gas onto a compact object moving through a uniform medium at Mach 3. Previous three-dimensional simulations have shown that such accretion is not stable, and that strong rotational 'disk-like' flows are generated and accreted on short time scales. We re-address this problem using overset spherical grids that provide a factor of seven improvement in spatial resolution over previous simulations. With our higher spatial resolution we found these 3D accretion flows remained remarkably axisymmetric. We examined two cases of accretion with different sized accretors. The larger accretor produced very steady flow, with the mass accretion rate varying by less than 0.02% over 30 flow times. The smaller accretor exhibited an axisymmetric breathing mode that modulated the mass accretion rate by a constant 20%. Nonetheless, the flow remained highly axisymmetric with only negligible accretion of angular momentum in both cases.

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 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 report results from a multi-wavelength analysis of the Galactic SNR RCW 86, the proposed remnant of the supernova of 185 A.D. We report new infrared observations from {\it Spitzer} and {\it WISE}, where the entire shell is detected at 24 and 22 $\mu$m. We fit the infrared flux ratios with models of collisionally heated ambient dust, finding post-shock gas densities in the non-radiative shocks of 2.4 and 2.0 cm$^{-3}$ in the SW and NW portions of the remnant, respectively. The Balmer-dominated shocks around the periphery of the shell, large amount of iron in the X-ray emitting ejecta, and lack of a compact remnant support a Type Ia origin for this remnant. From hydrodynamic simulations, the observed characteristics of RCW 86 are successfully reproduced by an off-center explosion in a low-density cavity carved by the progenitor system. This would make RCW 86 the first known case of a Type Ia supernova in a wind-blown bubble. The fast shocks ($> 3000$ km s$^{-1}$) observed in the NE are propagating in the low-density bubble, where the shock is just beginning to encounter the shell, while the slower shocks elsewhere have already encountered the bubble wall. The diffuse nature of the synchrotron emission in the SW and NW is due to electrons that were accelerated early in the lifetime of the remnant, when the shock was still in the bubble. Electrons in a bubble could produce gamma-rays by inverse-Compton scattering. The wind-blown bubble scenario requires a single-degenerate progenitor, which should leave behind a companion star.

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.

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.

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.

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.

Progenitors of Type Ia supernovae (SNe) have been predicted to modify their ambient circumstellar (CSM) and interstellar environments through the action of their powerful winds. While there is X-ray and optical evidence for circumstellar interaction in several remnants of Type Ia SNe, widespread evidence for such interaction in Type Ia SNe themselves has been lacking. We consider prospects for detection of CSM shells that have been predicted to be common around Type Ia SNe. Such shells are most easily detected in Na I absorption lines. Variable (declining) absorption is expected to occur soon after the explosion, primarily during the SN rise time, for shells located within 1 - 10 pc of a SN. The distance of the shell from the SN can be determined by measuring the time scale for line variability.

We revisit the flip-flop instability of two-dimensional planar accretion using high-fidelity numerical simulations. By starting from an initially steady-state axisymmetric solution, we are able to follow the growth of this overstability from small amplitudes. In the small-amplitude limit, before any transient accretion disk is formed, the oscillation period of the accretion shock is comparable to the Keplerian period at the Hoyle-Lyttleton accretion radius (R_a), independent of the size of the accreting object. The growth rate of the overstability increases dramatically with decreasing size of the accretor, but is relatively insensitive to the upstream Mach number of the flow. We confirm that the flip-flop does not require any gradient in the upstream flow. Indeed, a small density gradient as used in the discovery simulations has virtually no influence on the growth rate of the overstability. The ratio of specific heats does influence the overstability, with smaller gamma leading to faster growth of the instability. For a relatively large accretor (a radius of 0.037R_a) planar accretion is unstable for gamma = 4/3, but stable for gamma > 1.6. Planar accretion is unstable even for gamma = 5/3 provided the accretor has a radius of < 0.0025R_a. We also confirm that when the accretor is sufficiently small, the secular evolution is described by sudden jumps between states with counter-rotating quasi-Keplerian accretion disks.

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.

We present initial results of a 750 ks Chandra observation of the remnant of Kepler's supernova of AD 1604. The strength and prominence of iron emission, together with the absence of O-rich ejecta, demonstrate that Kepler resulted from a thermonuclear supernova, even though evidence for circumstellar interaction is also strong. We have analyzed spectra of over 100 small regions, and find that they fall into three classes. (1) The vast majority show Fe L emission between 0.7 and 1 keV and Si and S K alpha emission; we associate these with shocked ejecta. A few of these are found at or beyond the mean blast wave radius. (2) A very few regions show solar O/Fe abundance rations; these we associate with shocked circumstellar medium (CSM). Otherwise O is scarce. (3) A few regions are dominated by continuum, probably synchrotron radiation. Finally, we find no central point source, with a limit about 100 times fainter than the central object in Cas A. The evidence that the blast wave is interacting with CSM may indicate a Ia explosion in a more massive progenitor.

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.

Two-dimensional axisymmetric simulations have shown that the post-bounce accretion shock in core collapse supernovae is subject to the Spherical Accretion Shock Instability, or SASI. Recent three-dimensional simulations have revealed the existence of a non-axisymmetric mode of the SASI as well, where the postshock flow displays a spiral pattern. Here we investigate the growth of these spiral modes using two-dimensional simulations of the post-bounce accretion flow in the equatorial plane of a core-collapse supernova. By perturbing a steady-state model we are able to excite both one, two and three-armed spiral modes that grow exponentially with time, demonstrating that these are linearly unstable modes closely related to the original axisymmetric sloshing modes. By tracking the distribution of angular momentum, we show that these modes are able to efficiently separate the angular momentum of the accretion flow (which maintains a net angular momentum of zero), leading to a significant spin-up of the underlying accreting proto-neutron star.

We calculate the X-ray emission from both constant and time evolving shocked fast winds blown by the central stars of planetary nebulae (PNs) and compare with observations. Using spherically symmetric numerical simulations with radiative cooling, we calculate the flow structure, and the X-ray temperature and luminosity of the hot bubble formed by the shocked fast wind. We find that a constant fast wind gives results that are very close to those obtained from the self-similar solution. We show that in order for a fast shocked wind to explain the observed X-ray properties of PNs, rapid evolution of the wind is essential. More specifically, the mass loss rate of the fast wind should be high early on when the speed is ~300-700 km/s, and then it needs to drop drastically by the time the PN age reaches ~1000 yr. This implies that the central star has a very short pre-PN (post-AGB) phase.

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.

We introduce a million-second observation of the supernova remnant Cassiopeia A with the Chandra X-ray Observatory. The bipolar structure of the Si-rich ejecta (NE jet and SW counterpart) is clearly evident in the new images, and their chemical similarity is confirmed by their spectra. These are most likely due to jets of ejecta as opposed to cavities in the circumstellar medium, since we can reject simple models for the latter. The properties of these jets and the Fe-rich ejecta will provide clues to the explosion of Cas A.

We present a numerical investigation of the development of Rayleigh-Taylor instability at the interface between an expanding Pulsar Wind Nebula and its surrounding Supernova Remnant. These systems have long been thought to be naturally subject to this kind of instability, given their expansion behavior and the density jump at the contact discontinuity. High resolution images of the Crab Nebula at optical frequencies show the presence of a complex network of line-emitting filaments protruding inside the synchrotron nebula. These structures are interpreted as the observational evidence that Rayleigh-Taylor instability is in fact at work. The development of this instability in the regime appropriate to describe Supernova Remnant-Pulsar Wind Nebula systems is non-trivial. The conditions at the interface are likely close to the stability threshold, and the inclusion of the nebular magnetic field, which might play an important role in stabilizing the system, is essential to the modeling. If Rayleigh-Taylor features can grow efficiently a mixing layer in the outer portion of the nebula might form where most of the supernova material is confined. When a magnetic field close to equipartition is included we find that the interface is stable, and that even a weaker magnetic field affects substantially the growth and shape of the fingers.

We present high resolution spherically symmetric relativistic magnetohydrodynamical simulations of the evolution of a pulsar wind nebula inside the free expanding ejecta of the supernova progenitor. The evolution is followed starting from a few years after the supernova explosion and up to an age of the remnant of 1500 years. We consider different values of the pulsar wind magnetization parameter and also different braking indices for the spin-down process. We compare the numerical results with those derived through an approximate semi-analytical approach that allows us to trace the time evolution of the positions of both the pulsar wind termination shock and the contact discontinuity between the nebula and the supernova ejecta. We also discuss, whenever a comparison is possible, to what extent our numerical results agree with former self-similar models, and how these models could be adapted to take into account the temporal evolution of the system. The inferred magnetization of the pulsar wind could be an order of magnitude lower than that derived from time independent analytic models.

Pulsars, formed during supernova explosions, are known to be sources of relativistic magnetized winds whose interaction with the expanding supernova remnants (SNRs) gives rise to a pulsar wind nebula (PWN). We present spherically symmetric relativistic magnetohydrodynamics (RMHD) simulations of the interaction of a pulsar wind with the surrounding SNR, both in particle and magnetically dominated regimes. As shown by previous simulations, the evolution can be divided in three phases: free expansion, a transient phase characterized by the compression and reverberation of the reverse shock, and a final Sedov expansion. The evolution of the contact discontinuity between the PWN and the SNR (and consequently of the SNR itself) is almost independent of the magnetization of the nebula as long as the total (magnetic plus particle) energy is the same. However, a different behaviour of the PWN internal structure is observable during the compression-reverberation phase, depending on the degree of magnetization=2E The simulations were performed using the third order conservative scheme by Del Zanna et al. (2003).

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.

For pulsars similar to the one in the Crab Nebula, most of the energy input to the surrounding wind nebula occurs on a timescale of less than 1000 years; during this time, the nebula expands into freely expanding supernova ejecta. On a timescale 10,000 years, the interaction of the supernova with the surrounding medium drives a reverse shock front toward the center of the remnant, where it crushes the PWN (pulsar wind nebula). One- and two-dimensional, two-fluid simulations of the crushing and re-expansion phases of a PWN show that (1) these phases are subject to Rayleigh-Taylor instabilities that result in the mixing of thermal and nonthermal fluids, and (2) asymmetries in the surrounding interstellar medium give rise to asymmetries in the position of the PWN relative to the pulsar and explosion site. These effects are expected to be observable in the radio emission from evolved PWN because of the long lifetimes of radio emitting electrons. The scenario can explain the chaotic and asymmetric appearance of the Vela X PWN relative to the Vela pulsar without recourse to a directed flow from the vicinity of the pulsar. The displacement of the radio nebulae in G327.1--1.1, MSH15--56 (G326.3--1.8), G0.9+0.1, and W44 relative to the X-ray nebulae may be due to this mechanism. On timescales much greater than the nebular crushing time, the initial PWN may be mixed with thermal gas and become unobservable, so that even the radio emission is dominated by recently injected particles.

We employ hydrodynamic simulations to study the effects of high shock compression ratios, as expected for fast shocks with efficient particle acceleration, on the convective instability of driven waves in supernova remnants. We find that the instability itself does not depend significantly on the compression ratio, but because the width of the interaction region between the forward and reverse shocks can shrink significantly with increasing shock compression, we find that convective instabilities can reach all the way to the forward shock front if compression ratios are high enough.

We compare the predictions of hydrodynamic models of the LMC X-4 X-ray binary system with observations of UV P Cygni lines with the GHRS and STIS spectrographs on the Hubble Space Telescope. The hydrodynamic model determines density and velocity fields of the stellar wind, wind-compressed disk, accretion stream, Keplerian accretion disk, and accretion disk wind. We use a Monte Carlo code to determine the UV P Cygni line profiles by simulating the radiative transfer of UV photons that originate on the star and are scattered in the wind. The qualitative orbital variation predicted is similar to that observed, although the model fails to reproduce the strong orbital asymmetry (the observed absorption is strongest for phi>0.5). The model predicts a mid-eclipse X-ray spectrum, due almost entirely to Compton scattering, with a factor 4 less flux than observed with ASCA. We discuss how the model may need to be altered to explain the spectral variability of the system.

We present here numerical hydrodynamic simulations of line-driven accretion disk winds in cataclysmic variable systems. We calculate wind mass-loss rate, terminal velocities, and line profiles for CIV (1550 A) for various viewing angles. The models are 2.5-dimensional, include an energy balance condition, and calculate the radiation field as a function of position near an optically thick accretion disk. The model results show that centrifugal forces produce collisions of streamlines in the disk wind which in turn generate an enhanced density region, underlining the necessity of two dimensional calculations where these forces may be represented. For disk luminosity Ldisk = Lsun, white dwarf mass Mwd = 0.6 Msun, and white dwarf radii Rwd = 0.01 Rsun, we obtain a wind mass-loss rate of dMwind/dt = 8.0E-12 Msun/yr, and a terminal velocity of ~3000 km/s. The line profiles we obtain are consistent with observations in their general form, in particular in the maximum absorption at roughly half the terminal velocity for the blue-shifted component, in the magnitudes of the wind velocities implied by the absorption components, in the FWHM of the emission components, and in the strong dependence in inclination angle.

Observations of core-collapse supernovae (SNe) have revealed the presence of extensive mixing of radioactive material in SN ejecta. The mixing of radioactive material, mostly freshly synthesized Ni, is not complete, which leads to a two-phase SN ejecta structure. The low-density phase consists of Fe bubbles, created by the energy input from radioactive Co and Ni, surrounded by compressed high-density metal-rich ejecta. We report on the theoretical investigation of supernova remnant (SNR) dynamics with the two-phase SN ejecta. We first present 3-dimensional hydrodynamic simulations of a single Fe bubble immersed in an outer ejecta envelope, and compare the results with previous work on shock-cloud interactions. We then consider randomly distributed Fe bubbles with an average volume filling fraction of 1/2. We find that the presence of Fe bubbles leads to vigorous turbulence and mixing of Fe with other heavy elements and with the ambient normal-abundance gas. The turbulent energy can be an order of magnitude larger than in the case of smooth ejecta. A significant fraction of the shocked ejecta is found in narrow filaments and clumps moving with radial velocities larger than the velocity of the forward shock. Observational consequences of the two-phase ejecta on SNR X-ray spectra and images are briefly mentioned.

The effects of binary tidal forces on transport within an accretion disk are studied with a time-dependent hydrodynamical model of a two-dimensional isothermal accretion disk. Tidal forces quickly truncate the accretion disk to radii of order half the average radius of the Roche lobe, and excite a two-armed spiral wave that remains stationary in the rotating reference frame of the binary system. We measure an effective alpha of order 0.1 near the outer edge of the disk in all of our models, independent of the mass ratio, Mach number, and radial density profile. However, in cold disks with high Mach number, the effective alpha drops rapidly with decreasing radius such that it falls below our threshold of measurement (roughly .001) at a radius of only one third the tidal radius. In warmer disks where the Mach numbers remain below 20, we can measure an effective alpha down to radii 10 times smaller than the maximum size of the disk.

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.

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.

Multiwavelength observations of Type II supernovae have shown evidence for the interaction of supernovae with the dense slow winds from the red supergiant progenitor stars. Observations of planetary nebulae and the nebula around SN 1987A show that the slow winds from extended stars frequently have an axisymme- tric structure with a high density in the equatorial plane. We have carried out numerical calculations of the interaction of a supernova with such an axisymme- tric density distribution. For small values of the angular density gradient at the pole, the asymmetry in the interaction shell is greater than, but close to, that expected from purely radial motion. If the angular density gradient is above a moderate value, the flow qualitatively changes and a protrusion emerges along the axis. For a power-law supernova density profile, the flow approaches a self-similar state in which the protrusion length is $2-4$ times the radius of the main shell. The critical density gradient is larger for steeper density profiles of the ejecta. Most of our calculations are axisymmetric, but we have carried out a 3-dimensional calculation to show that the protrusion is not a numerical artifact along the symmetry axis. For typical supernova parameters, the protrusions take $\gtrsim$ several years to develop. The appearance of the shell with protrusions is similar to that observed in VLBI radio images of the remnant 41.9 +58 in M82 and, possibly, of SN 1986J. We also considered the possibility of asymmetric ejecta and found that it had a relatively small effect on the asymmetry of the interaction region.