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.

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

Massive stars with a core-halo structure are interesting objects for stellar physics and hydrodynamics. Using simulations for stellar evolution, radiation hydrodynamics, and radiative transfer, we study the explosion of stars with an extended and tenuous envelope (i.e., stars in which 95% of the mass is contained within 10% of the surface radius or less). We consider both H-rich supergiant and He-giant progenitors resulting from close-binary evolution and dying with a final mass of 2.8-5Msun. An extended envelope causes the supernova (SN) shock to brake and a reverse shock to form, sweeping core material into a dense shell. The shock deposited energy, which suffers little degradation from expansion, is trapped in ejecta layers of moderate optical depth, thereby enhancing the SN luminosity at early times. With the delayed 56Ni heating, we find that the resulting optical and near-IR light curves all exhibit a double-peak morphology. We show how an extended progenitor can explain the blue and featureless optical spectra of some Type IIb and Ib SNe. The dense shell formed by the reverse shock leads to line profiles with a smaller and near-constant width. This ejecta property can explain the statistically narrower profiles of Type IIb compared to Type Ib SNe, as well as the peculiar Halpha profile seen in SN1993J. At early times, our He-giant star explosion model shows a high luminosity, a blue color, and featureless spectra reminiscent of the Type Ib SN2008D, suggesting a low-mass progenitor.

The progenitors of low-luminosity Type II-Plateau supernovae (SNe II-P) are believed to be red supergiant (RSG) stars, but there is much disparity in the literature concerning their mass at core collapse and therefore on the main sequence. Here, we model the SN radiation arising from the low-energy explosion of RSG stars of 12, 25, and 27 M$_{\odot}$ on the main sequence and formed through single star evolution. Despite the narrow range in ejecta kinetic energy (2.5$-$4.2$\times$10$^{50}$ erg) in our model set, the SN observables from our three models are significantly distinct, reflecting the differences in progenitor structure (e.g., surface radius, H-rich envelope mass, He-core mass). Our higher mass RSG stars give rise to Type II SNe that tend to have bluer colors at early times, a shorter photospheric phase, and a faster declining $V$-band light curve (LC) more typical of Type II-linear SNe, in conflict with the LC plateau observed for low-luminosity SNe II. The complete fallback of the CO core in the low-energy explosions of our high mass RSG stars prevents the ejection of any ${}^{56}$Ni (nor any core O or Si), in contrast to low-luminosity SNe II-P, which eject at least 0.001 M$_{\odot}$ of ${}^{56}$Ni. In contrast to observations, type II SN models from higher mass RSGs tend to show an H$\alpha$ absorption that remains broad at late times (due to a larger velocity at the base of the H-rich envelope). In agreement with the analyses of pre-explosion photometry, we conclude that low-luminosity SNe II-P likely arise from low-mass rather than high-mass RSG stars.

Using 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations, we study the ejecta properties required to match the early and late-time photometric and spectroscopic properties of supernovae (SNe) associated with long-duration gamma-ray bursts (LGRBs). To match the short rise time, narrow light curve peak, and extremely broad spectral lines of SN1998bw requires a model with <3Msun ejecta but a high explosion energy of a few 1e52erg and 0.5Msun of Ni56. However the relatively high luminosity, the presence of narrow spectral lines of intermediate mass elements, and the low ionization at the nebular stage are matched with a more standard C-rich Wolf-Rayet (WR) star explosion, with an ejecta of >10Msun, an explosion energy >1e51erg, and only 0.1Msun of Ni56. As the two models are mutually exclusive, the breaking of spherical symmetry is essential to match the early/late photometric/spectroscopic properties of SN1998bw. This conclusion confirms the notion that the ejecta of SN1998bw is aspherical on large scales. More generally, with asphericity, the energetics and Ni56 mass of LGRB/SNe are reduced and their ejecta mass is increased, favoring a massive fast-rotating Wolf-Rayet star progenitor. Contrary to persisting claims in favor of the proto-magnetar model for LGRB/SNe, such progenitor/ejecta properties are compatible with collapsar formation. Ejecta properties of LGRB/SNe inferred from 1D radiative-transfer modeling are fundamentally flawed.

Supernova (SN) 2008bk is a well observed low-luminosity Type II event visually associated with a low-mass red-supergiant progenitor. To model SN 2008bk, we evolve a 12Msun star from the main sequence until core collapse, when it has a total mass of 9.88Msun, a He-core mass of 3.22Msun, and a radius of 502Rsun. We then artificially trigger an explosion that produces 8.29Msun of ejecta with a total energy of 2.5x10^50erg and ~0.009Msun of 56Ni. We model the subsequent evolution of the ejecta with non-Local-Thermodynamic-Equilibrium time-dependent radiative transfer. Although somewhat too luminous and energetic, this model reproduces satisfactorily the multi-band light curves and multi-epoch spectra of SN 2008bk, confirming the suitability of a low-mass massive star progenitor. As in other low-luminosity SNe II, the structured Halpha profile at the end of the plateau phase is probably caused by BaII 6496.9A rather than asphericity. We discuss the sensitivity of our results to changes in progenitor radius and mass, as well as chemical mixing. A 15% increase in progenitor radius causes a 15% increase in luminosity and a 0.2mag V-band brightening of the plateau but leaves its length unaffected. An increase in ejecta mass by 10% lengthens the plateau by ~10d. Chemical mixing introduces slight changes to the bolometric light curve, limited to the end of the plateau, but has a large impact on colours and spectra at nebular times.

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for a large grid of supernovae (SNe) IIb/Ib/Ic that result from the terminal explosion of the mass donor in a close-binary system. Our sample covers ejecta masses $M_{\rm e}$ of 1.7$-$5.2M$_\odot$, kinetic energies $E_{\rm kin}$ of 0.6$-$5.0$\times$10$^{51}$erg, and $^{56}$Ni masses of 0.05$-$0.30M$_\odot$. We find a strong correlation between the $^{56}$Ni mass and the photometric properties at maximum, and between the rise time to bolometric maximum and the post-maximum decline rate. We confirm the small scatter in ($V-R$) at 10d past $R$-band maximum. The quantity $V_{\rm m} \equiv \sqrt{2E_{\rm kin}/M_{\rm e}}$ is comparable to the Doppler velocity measured from HeI 5875\AA\ at maximum in SNe IIb/Ib, although some scatter arises from the uncertain level of chemical mixing. The OI7772\AA\ line may be used for SNe Ic, but the correspondence deteriorates with higher ejecta mass/energy. We identify a temporal reversal of the Doppler velocity at maximum absorption in the $\sim$1.05$\mu$m feature in all models. The reversal is due to HeI alone and could serve as a test for the presence of helium in SNe Ic. Because of variations in composition and ionisation, the ejecta opacity shows substantial variations with both velocity and time. This is in part the origin of the offset between our model light curves and the predictions from the Arnett model.

Using radiation-hydrodynamics and radiative-transfer simulations, we explore the origin of the spectral diversity of interacting supernovae (SNe) of type IIn. We revisit SN1994W and investigate the dynamical configurations that can give rise to spectra with narrow lines at all times. We find that a standard ~10Msun 10^51erg SN ejecta ramming into a 0.4Msun dense CSM is inadequate for SN1994W, as it leads to the appearance of broad lines at late times. This structure, however, generates spectra that exhibit the key morphological changes seen in SN1998S. For SN1994W, we consider a completely different configuration, which involves the interaction at a large radius of a low mass inner shell with a high mass outer shell. Such a structure may arise in an 8-12Msun star from a nuclear flash (e.g., of Ne) followed within a few years by core collapse. Our simulations show that the large mass of the outer shell leads to the complete braking of the inner shell material, the formation of a slow dense shell, and the powering of a luminous SN IIn, even for a low inner shell energy. Early on, our model line profiles are typical of SNe IIn, exhibiting narrow cores and broad electron-scattering wings. As observed in SN1994W, they also remain narrow at late times. Our SN1994W model invokes two low energy ejections, both atypical of observed massive stars, and illustrates the diversity of configurations leading to SNe IIn. These results also highlight the importance of spectra to constrain the dynamical properties and understand the origin of SNe IIn.

We present 1-D non-Local-Thermodynamic-Equilibrium time-dependent radiative-transfer simulations for supernovae (SNe) of type IIb, Ib, and Ic that result from the terminal explosion of the mass donor in a close-binary system. Here, we select three ejecta with a total kinetic energy of ~1.2e51erg, but characterised by different ejecta masses (2-5Msun), composition, and chemical mixing. The type IIb/Ib models correspond to the progenitors that have retained their He-rich shell at the time of explosion. The type Ic model arises from a progenitor that has lost its helium shell, but retains 0.32Msun of helium in a CO-rich core of 5.11Msun. We discuss their photometric and spectroscopic properties during the first 2-3 months after explosion, and connect these to their progenitor and ejecta properties including chemical stratification. For these three models, Arnett's rule overestimates the 56Ni mass by ~50% while the procedure of Katz et al., based on an energy argument, yields a more reliable estimate. The presence of strong CI lines around 9000A prior to maximum is an indicator that the pre-SN star was under-abundant in helium. As noted by others, the 1.08micron feature is a complex blend of CI, MgII, and HeI lines, which makes the identification of He uncertain in SNe Ibc unless other HeI lines can be identified. Our models show little scatter in (V-R) colour 10d after R-band maximum. We also address a number of radiative transfer properties of SNe Ibc, including the notion of a photosphere, the inference of a representative ejecta expansion rate, spectrum formation, blackbody fits and "correction factors".

Type Ia supernovae (SNe Ia), thermonuclear explosions of carbon-oxygen white dwarfs (CO-WDs), are currently the best cosmological "standard candles", but the triggering mechanism of the explosion is unknown. It was recently shown that the rate of head-on collisions of typical field CO-WDs in triple systems may be comparable to the SNe Ia rate. Here we provide evidence supporting a scenario in which the majority of SNe Ia are the result of such head-on collisions of CO-WDs. In this case, the nuclear detonation is due to a well understood shock ignition, devoid of commonly introduced free parameters such as the deflagration velocity or transition to detonation criteria. By using two-dimensional hydrodynamical simulations with a fully resolved ignition process, we show that zero-impact-parameter collisions of typical CO-WDs with masses $0.5-1\,M_{\odot}$ result in explosions that synthesize $^{56}$Ni masses in the range of $\sim0.1-1\,M_{\odot}$, spanning the wide distribution of yields observed for the majority of SNe Ia. All collision models yield the same late-time ($>60$ days since explosion) bolometric light curve when normalized by $^{56}$Ni masses (to better than $30\%$), in agreement with observations. The calculated widths of the $^{56}$Ni-mass-weighted-line-of-sight velocity distributions are correlated with the calculated $^{56}$Ni yield, agreeing with the observed correlation. The strong correlation, shown here for the first time, between $^{56}$Ni yield and total mass of the colliding CO-WDs (insensitive to their mass ratio), is suggestive as the source for the continuous distribution of observed SN Ia features, possibly including the Philips relation.

We explore the properties of Type II-Plateau (II-P) supernovae (SNe) together with their red-supergiant (RSG) star progenitors. Using MESA STAR, we modulate the parameters (e.g., mixing length, overshoot, rotation, metallicity) that control the evolution of a 15Msun main-sequence star to produce a variety of physical pre-SN models and SN II-P ejecta. We extend previous modeling of SN II-P radiation to include photospheric and nebular phases, as well as multi-band light curves and spectra. Our treatment does not assume local thermodynamic equilibrium, is time dependent, treats explicitly the effects of line blanketing, and incorporates non-thermal processes. We find that the color properties of SNe II-P require large model atoms for FeI and FeII, much larger than adopted in Dessart & Hillier (2011). The color properties also imply RSG progenitors of limited extent (~500Rsun) --- larger progenitor stars produce a SN II-P radiation that remains too blue for too long. This finding calls for a reduction of RSG radii, perhaps through a strengthening of convective energy transport in RSG envelopes. Increased overshoot and rotation reduce the ratio of ejecta to helium-core mass, similarly to an increase in main-sequence mass, and thus complicate the inference ofprogenitor masses. In contrast to the great sensitivity on progenitor radius, SN II-P color evolution appears insensitive to variations in explosion energy. Finally, we document the numerous SN II-P signatures that vary with progenitor metallicity, revealing their potential for metallicity determinations in the nearby and distant Universe.

We present non-LTE time-dependent radiative-transfer simulations of pair-instability supernovae (PISNe) stemming from red-supergiant (RSG), blue-supergiant (BSG) and Wolf-Rayet (WR) star rotation-free progenitors born in the mass range 160-230Msun, at 10^-4 Zsun. Although subject to uncertainties in convection and stellar mass-loss rates, our initial conditions come from physically-consistent models that treat evolution from the main-sequence, the onset of the pair-production instability, and the explosion phase. With our set of input models characterized by large 56Ni and ejecta masses, and large kinetic energies, we recover qualitatively the Type II-Plateau, II-peculiar, and Ib/c light-curve morphologies, although they have larger peak bolometric luminosities (~10^9 to 10^10 Lsun) and a longer duration (~200d). We discuss the spectral properties for each model during the photospheric and nebular phases, including Balmer lines in II-P and II-pec at early times, the dominance of lines from intermediate-mass-elements (IMEs) near the bolometric maximum, and the strengthening of metal line blanketing thereafter. Having similar He-core properties, all models exhibit similar post-peak spectra that are strongly blanketed by FeII and FeI lines, characterized by red colors, and that arise from photospheres/ejecta with a temperature of <4000K. Combined with the modest line widths after bolometric peak, these properties contrast with those of known super-luminous SNe suggesting that PISNe are yet to be discovered. Being reddish, PISNe will be difficult to observe at high redshift except when they stem from RSG explosions, in which case they could be used as metallicity probes and distance indicators.

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

Much uncertainty surrounds the origin of super-luminous supernovae (SNe). Motivated by the discovery of the Type Ic SN2007bi, we study its proposed association with a pair-instability SN (PISN). We compute stellar-evolution models for primordial ~200Msun stars, simulating the implosion/explosion due to the pair-production instability, and use them as inputs for detailed non-LTE time-dependent radiative-transfer simulations that include non-local energy deposition and non-thermal processes. We retrieve the basic morphology of PISN light curves from red-supergiant, blue-supergiant, and Wolf-Rayet (WR) star progenitors. Although we confirm that a progenitor 100Msun helium core (PISN model He100) fits well the SN2007bi light curve, the low ratios of its kinetic energy and 56Ni mass to the ejecta mass, similar to standard core-collapse SNe, conspire to produce cool photospheres, red spectra subject to strong line blanketing, and narrow line profiles, all conflicting with SN2007bi observations. He-core models of increasing 56Ni-to-ejecta mass ratio have bluer spectra, but still too red to match SN2007bi, even for model He125 -- the effect of 56Ni heating is offset by the associated increase in blanketing. In contrast, the delayed injection of energy by a magnetar represents a more attractive alternative to reproduce the blue, weakly-blanketed, and broad-lined spectra of super-luminous SNe. The extra heat source is free of blanketing and is not explicitly tied to the ejecta. Experimenting with a ~9Msun WR-star progenitor, initially exploded to yield a ~1.6B SN Ib/c ejecta but later influenced by tunable magnetar-like radiation, we produce a diversity of blue spectral morphologies reminiscent of SN2007bi, the peculiar Type Ib SN2005bf, and super-luminous SN2005ap-like events.

We show that converging spherical and cylindrical shock waves may ignite a detonation wave in a combustible medium, provided the radius at which the shocks become strong exceeds a critical radius, R_c. An approximate analytic expression for R_c is derived for an ideal gas equation of state and a simple (power-law-Arrhenius) reaction law, and shown to reproduce the results of numerical solutions. For typical acetylene-air experiments we find R_c~0.1 mm (spherical) and R_c~1 mm (cylindrical). We suggest that the deflagration to detonation transition (DDT) observed in these systems may be due to converging shocks produced by the turbulent deflagration flow, which reaches sub (but near) sonic velocities on scales >>R_c. Our suggested mechanism differs from that proposed by Zel'dovich et al., in which a fine-tuned spatial gradient in the chemical induction time is required to be maintained within the turbulent deflagration flow. Our analysis may be readily extended to more complicated equations of state and reaction laws. An order of magnitude estimate of R_c within a white dwarf at the pre-detonation conditions believed to lead to Type Ia supernova explosions is 0.1 km, suggesting that our proposed mechanism may be relevant for DDT initiation in these systems. The relevance of our proposed ignition mechanism to DDT initiation may be tested by both experiments and numerical simulations.

A unique feature of deflagration-to-detonation (DDT) white dwarf explosion models of SNe of type Ia is the presence of a strong shock wave propagating through the outer envelope. We consider the early emission expected in such models, which is produced by the expanding shock-heated outer part of the ejecta and precedes the emission driven by radioactive decay. We expand on earlier analyses by considering the modification of the pre-detonation density profile by the weak-shocks generated during the deflagration phase, the time evolution of the opacity, and the deviation of the post-shock equation of state from that obtained for radiation pressure domination. A simple analytic model is presented and shown to provide an acceptable approximation to the results of 1D numerical DDT simulations. Our analysis predicts a thousand second long UV/optical flash with a luminosity of ~1 to 3*1e39 erg/s. Lower luminosity corresponds to faster (turbulent) deflagration velocity. The predicted luminosity of the UV flash is an order of magnitude lower than that of earlier estimates, and is expected to be strongly suppressed at times longer than an hour due to the deviation from pure radiation domination.

We present non-LTE time-dependent radiative-transfer simulations of supernova (SN) IIb/Ib/Ic spectra and light curves, based on ~1B-energy piston-driven ejecta, with and without 56Ni, produced from single and binary Wolf-Rayet (W-R) stars evolved at solar and sub-solar metallicities. Our bolometric light curves show a 10-day long post-breakout plateau with a luminosity of 1-5x10^7Lsun. In our 56Ni-rich models, with ~3Msun ejecta masses, this plateau precedes a 20-30-day long re-brightening phase initiated by the outward-diffusing heat wave powered by radioactive decay at depth. In low ejecta-mass models with moderate mixing, Gamma-ray leakage starts as early as ~50d after explosion and causes the nebular luminosity to steeply decline by ~0.02mag/d. Such signatures, which are observed in standard SNe IIb/Ib/Ic, are consistent with low-mass progenitors derived from a binary-star population. We propose that the majority of stars with an initial mass ~<20Msun yield SNe II-P if 'effectively" single, SNe IIb/Ib/Ic if part of a close binary system, and SN-less black holes if more massive. Our ejecta, with outer hydrogen mass fractions as low as ~>0.01 and a total hydrogen mass of ~>0.001Msun, yield the characteristic SN IIb spectral morphology at early times. However, by ~15d after the explosion, only Halpha may remain as a weak absorption feature. Our binary models, characterised by helium surface mass fractions of ~>0.85, systematically show HeI lines during the post-breakout plateau, irrespective of the 56Ni abundance. Synthetic spectra show a strong sensitivity to metallicity, which offers the possibility to constrain it directly from SN spectroscopic modelling.

We present new results from the only 2D multi-group, multi-angle calculations of core-collapse supernova evolution. The first set of results from these calculations was published in Ott et al. (2008). We have followed a nonrotating and a rapidly rotating 20 solar mass model for ~400 ms after bounce. We show that the radiation fields vary much less with angle than the matter quantities in the region of net neutrino heating. This obtains because most neutrinos are emitted from inner radiative regions and because the specific intensity is an integral over sources from many angles at depth. The latter effect can only be captured by multi-angle transport. We then compute the phase relationship between dipolar oscillations in the shock radius and in matter and radiation quantities throughout the postshock region. We demonstrate a connection between variations in neutrino flux and the hydrodynamical shock oscillations, and use a variant of the Rayleigh test to estimate the detectability of these neutrino fluctuations in IceCube and Super-K. Neglecting flavor oscillations, fluctuations in our nonrotating model would be detectable to ~10 kpc in IceCube, and a detailed power spectrum could be measured out to ~5 kpc. These distances are considerably lower in our rapidly rotating model or with significant flavor oscillations. Finally, we measure the impact of rapid rotation on detectable neutrino signals. Our rapidly rotating model has strong, species-dependent asymmetries in both its peak neutrino flux and its light curves. The peak flux and decline rate show pole-equator ratios of up to ~3 and ~2, respectively.

The collapse of a massive star's core, followed by a neutrino-driven, asymmetric supernova explosion, can naturally lead to pulsar recoils and neutron star kicks. Here, we present a two-dimensional, radiation-hydrodynamic simulation in which core collapse leads to significant acceleration of a fully-formed, nascent neutron star (NS) via an induced, neutrino-driven explosion. During the explosion, a ~10% anisotropy in the low-mass, high-velocity ejecta lead to recoil of the high-mass neutron star. At the end of our simulation, the NS has achieved a velocity of ~150 km s$^{-1}$ and is accelerating at ~350 km s$^{-2}$, but has yet to reach the ballistic regime. The recoil is due almost entirely to hydrodynamical processes, with anisotropic neutrino emission contributing less than 2% to the overall kick magnitude. Since the observed distribution of neutron star kick velocities peaks at ~300-400 km s$^{-1}$, recoil due to anisotropic core-collapse supernovae provides a natural, non-exotic mechanism with which to obtain neutron star kicks.

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

We present radiation-hydrodynamics simulations of core-collapse supernova (SN) explosions, artificially generated by driving a piston at the base of the envelope of a rotating or non-rotating red-supergiant progenitor star. We search for trends in ejecta kinematics in the resulting Type II-Plateau (II-P) SN, exploring dependencies with explosion energy and pre-SN stellar-evolution model. We recover the trivial result that larger explosion energies yield larger ejecta velocities in a given progenitor. However, we emphasise that for a given explosion energy, the increasing helium-core mass with main-sequence mass of such Type II-P SN progenitors leads to ejection of core-embedded oxygen-rich material at larger velocities. We find that the photospheric velocity at 15d after shock breakout is a good and simple indicator of the explosion energy in our selected set of pre-SN models. This measurement, combined with the width of the nebular-phase OI6303-6363A line, can be used to place an upper-limit on the progenitor main-sequence mass. Using the results from our simulations, we find that the current, but remarkably scant, late-time spectra of Type II-P SNe support progenitor main-sequence masses inferior to ~20Msun and thus, corroborate the inferences based on the direct, but difficult, progenitor identification in pre-explosion images. The narrow width of OI6303-6363A in Type II-P SNe with nebular spectra does not support high-mass progenitors in the range 25-30Msun. Combined with quantitative spectroscopic modelling, such diagnostics offer a means to constrain the main-sequence mass of the progenitor, the mass fraction of the core ejected, and thus, the mass of the compact remnant formed.

Observations of transient phenomena in the Universe reveal a spectrum of mass-ejection properties associated with massive stars, covering from Type II/Ib/Ic core-collapse supernovae (SNe) to giant eruptions of Luminous Blue Variables (LBV) and optical transients. Here, we hypothesize that a fraction of these phenomena may have an explosive origin, the distinguishing ingredient being the ratio of the prompt energy release E_dep to the envelope binding energy E_binding. Using one-dimensional one-group radiation hydrodynamics and a set of 10-25Msun, massive-star models, we explore the dynamical response of a stellar envelope subject to a strong, sudden, and deeply-rooted energy release. Following energy deposition, a shock systematically forms, crosses the progenitor envelope on a day timescale, and breaks-out with a signal of hour-to-days duration and a 10^5-10^11 Lsun luminosity. For E_dep > E_binding, full envelope ejection results with a SN-like bolometric luminosity and kinetic energy, modulations being commensurate to the energy deposited and echoing the diversity of Type II-Plateau SNe. For E_dep ~ E_binding, partial envelope ejection results with a small expansion speed, and a more modest but year-long luminosity plateau, reminiscent of LBV eruptions or so-called SN impostors. For E_dep < E_binding, we obtain a "puffed-up" star, secularly relaxing back to thermal equilibrium. In parallel with gravitational collapse and Type II SNe, we argue that the thermonuclear combustion of merely a few 0.01Msun of C/O could power a wide range of explosions/eruptions in loosely-bound massive stars, as those in the 8-12Msun range, or in more massive ones owing to their proximity to the Eddington limit and/or critical rotation.

Core-collapse supernovae are among Nature's most energetic events. They mark the end of massive star evolution and pollute the interstellar medium with the life-enabling ashes of thermonuclear burning. Despite their importance for the evolution of galaxies and life in the universe, the details of the core-collapse supernova explosion mechanism remain in the dark and pose a daunting computational challenge. We outline the multi-dimensional, multi-scale, and multi-physics nature of the core-collapse supernova problem and discuss computational strategies and requirements for its solution. Specifically, we highlight the axisymmetric (2D) radiation-MHD code VULCAN/2D and present results obtained from the first full-2D angle-dependent neutrino radiation-hydrodynamics simulations of the post-core-bounce supernova evolution. We then go on to discuss the new code Zelmani which is based on the open-source HPC Cactus framework and provides a scalable AMR approach for 3D fully general-relativistic modeling of stellar collapse, core-collapse supernovae and black hole formation on current and future massively-parallel HPC systems. We show Zelmani's scaling properties to more than 16,000 compute cores and discuss first 3D general-relativistic core-collapse results.

We present VULCAN/2D multi-group flux-limited-diffusion radiation hydrodynamics simulations of binary neutron star (BNS) mergers, using the Shen equation of state, covering ~100 ms, and starting from azimuthal-averaged 2D slices obtained from 3D SPH simulations of Rosswog & Price for 1.4 Msun (baryonic) neutron stars with no initial spins, co-rotating spins, and counter-rotating spins. Snapshots are post-processed at 10 ms intervals with a multi-angle neutrino-transport solver. We find polar-enhanced neutrino luminosities, dominated by $\bar{\nu}_e$ and ``$\nu_\mu$'' neutrinos at peak, although $\nu_e$ emission may be stronger at late times. We obtain typical peak neutrino energies for $\nu_e$, $\bar{\nu}_e$, and ``$\nu_\mu$'' of ~12, ~16, and ~22 MeV. The super-massive neutron star (SMNS) formed from the merger has a cooling timescale of ~1 s. Charge-current neutrino reactions lead to the formation of a thermally-driven bipolar wind with <$\dot{M}$> ~10$^{-3}$ Msun/s, baryon-loading the polar regions, and preventing any production of a GRB prior to black-hole formation. The large budget of rotational free energy suggests magneto-rotational effects could produce a much greater polar mass loss. We estimate that ~10$^{-4}$ Msun of material with electron fraction in the range 0.1-0.2 become unbound during this SMNS phase as a result of neutrino heating. We present a new formalism to compute the $\nu_i\bar{\nu}_i$ annihilation rate based on moments of the neutrino specific intensity computed with our multi-angle solver. Cumulative annihilation rates, which decay as $t^{-1.8}$, decrease over our 100 ms window from a few 10$^{50}$ to ~10$^{49}$ erg/s, equivalent to a few 10$^{54}$ to ~10$^{53}$ $e^-e^+$ pairs per second.

We perform axisymmetric (2D) multi-angle, multi-group neutrino radiation-hydrodynamic calculations of the postbounce phase of core-collapse supernovae using a genuinely 2D discrete-ordinate (S_n) method. We follow the long-term postbounce evolution of the cores of one nonrotating and one rapidly-rotating 20-solar-mass stellar model for ~400 milliseconds from 160 ms to ~550 ms after bounce. We present a multi-D analysis of the multi-angle neutrino radiation fields and compare in detail with counterpart simulations carried out in the 2D multi-group flux-limited diffusion (MGFLD) approximation to neutrino transport. We find that 2D multi-angle transport is superior in capturing the global and local radiation-field variations associated with rotation-induced and SASI-induced aspherical hydrodynamic configurations. In the rotating model, multi-angle transport predicts much larger asymptotic neutrino flux asymmetries with pole to equator ratios of up to ~2.5, while MGFLD tends to sphericize the radiation fields already in the optically semi-transparent postshock regions. Along the poles, the multi-angle calculation predicts a dramatic enhancement of the neutrino heating by up to a factor of 3, which alters the postbounce evolution and results in greater polar shock radii and an earlier onset of the initially rotationally weakened SASI. In the nonrotating model, differences between multi-angle and MGFLD calculations remain small at early times when the postshock region does not depart significantly from spherical symmetry. At later times, however, the growing SASI leads to large-scale asymmetries and the multi-angle calculation predicts up to 30% higher average integral neutrino energy deposition rates than MGFLD.

Recent stellar evolutionary calculations of low-metallicity massive fast-rotating main-sequence stars yield iron cores at collapse endowed with high angular momentum. It is thought that high angular momentum and black hole formation are critical ingredients of the collapsar model of long-soft gamma-ray bursts (GRBs). Here, we present 2D multi-group, flux-limited-diffusion MHD simulations of the collapse, bounce, and immediate post-bounce phases of a 35-Msun collapsar-candidate model of Woosley & Heger. We find that, provided the magneto-rotational instability (MRI) operates in the differentially-rotating surface layers of the millisecond-period neutron star, a magnetically-driven explosion ensues during the proto-neutron star phase, in the form of a baryon-loaded non-relativistic jet, and that a black hole, central to the collapsar model, does not form. Paradoxically, and although much uncertainty surrounds stellar mass loss, angular momentum transport, magnetic fields, and the MRI, current models of chemically homogeneous evolution at low metallicity yield massive stars with iron cores that may have too much angular momentum to avoid a magnetically-driven, hypernova-like, explosion in the immediate post-bounce phase. We surmise that fast rotation in the iron core may inhibit, rather than enable, collapsar formation, which requires a large angular momentum not in the core but above it. Variations in the angular momentum distribution of massive stars at core collapse might explain both the diversity of Type Ic supernovae/hypernovae and their possible association with a GRB. A corollary might be that, rather than the progenitor mass, the angular momentum distribution, through its effect on magnetic field amplification, distinguishes these outcomes.

We present 2D multi-group flux-limited diffusion magnetohydrodynamics (MHD) simulations of the Accretion-Induced Collapse (AIC) of a rapidly-rotating white dwarf. We focus on the dynamical role of MHD processes after the formation of a millisecond-period protoneutron star. We find that including magnetic fields and stresses can lead to a powerful explosion with an energy of a few Bethe, rather than a weak one of at most 0.1 Bethe, with an associated ejecta mass of ~0.1Msun, instead of a few 0.001Msun. The core is spun down by ~30% within 500ms after bounce, and the rotational energy extracted from the core is channeled into magnetic energy that generates a strong magnetically-driven wind, rather than a weak neutrino-driven wind. Baryon loading of the ejecta, while this wind prevails, precludes it from becoming relativistic. This suggests that a GRB is not expected to emerge from such AICs during the early protoneutron star phase, except in the unlikely event that the massive white dwarf has sufficient mass to lead to black hole formation. In addition, we predict both negligible 56Ni-production (that should result in an optically-dark, adiabatically-cooled explosion) and the ejection of 0.1Msun of material with an electron fraction of 0.1-0.2. Such pollution by neutron-rich nuclei puts strong constraints on the possible rate of such AICs. Moreover, being free from ``fallback,'' such highly-magnetized millisecond-period protoneutron stars may later become magnetars, and the magnetically-driven winds may later transition to Poynting-flux-dominated, relativistic winds, eventually detectable as GRBs at cosmological distances. However, the low expected event rate of AICs will constrain them to be, at best, a small subset of GRB and/or magnetar progenitors.

We present here the first 2D rotating, multi-group, radiation magnetohydrodynamics (RMHD) simulations of supernova core collapse, bounce, and explosion. In the context of rapid rotation, we focus on the dynamical effects of magnetic stresses and the creation and propagation of MHD jets. We find that a quasi-steady state can be quickly established after bounce, during which a well-collimated MHD jet is maintained by continuous pumping of power from the differentially rotating core. If the initial spin period of the progenitor core is $\sles$ 2 seconds, the free energy reservoir in the secularly evolving protoneutron star is adequate to power a supernova explosion, and may be enough for a hypernova. The jets are well collimated by the infalling material and magnetic hoop stresses, and maintain a small opening angle. We see evidence of sausage instabilities in the emerging jet stream. Neutrino heating is sub-dominant in the rapidly rotating models we explore, but can contribute 10$-$25% to the final explosion energy. Our simulations suggest that even in the case of modest or slow rotation, a supernova explosion might be followed by a secondary, weak MHD jet explosion, which, because of its weakness may to date have gone unnoticed in supernova debris. Furthermore, we suggest that the generation of a non-relativistic MHD precursor jet during the early protoneutron star/supernova phase is implicit in both the collapsar and "millisecond magnetar" models of GRBs. The multi-D, multi-group, rapidly rotating RMHD simulations we describe here are a start along the path towards more realistic simulations of the possible role of magnetic fields in some of Nature's most dramatic events.

We report a new finite-difference scheme for two-dimensional magnetohydrodynamics (MHD) simulations, with and without rotation, in unstructured grids with quadrilateral cells. The new scheme is implemented within the code VULCAN/2D, which already includes radiation-hydrodynamics in various approximations and can be used with arbitrarily moving meshes (ALE). The MHD scheme, which consists of cell-centered magnetic field variables, preserves the nodal finite difference representation of $div(\bB)$ by construction, and therefore any initially divergence-free field remains divergence-free through the simulation. In this paper, we describe the new scheme in detail and present comparisons of VULCAN/2D results with those of the code ZEUS/2D for several one-dimensional and two-dimensional test problems. The code now enables two-dimensional simulations of the collapse and explosion of the rotating, magnetic cores of massive stars. Moreover, it can be used to simulate the very wide variety of astrophysical problems for which multi-D radiation-magnetohydrodynamics (RMHD) is relevant.

In this paper, we bring together various of our published and unpublished findings from our recent 2D multi-group, flux-limited radiation hydrodynamic simulations of the collapse and explosion of the cores of massive stars. Aided by 2D and 3D graphical renditions, we motivate the acoustic mechanism of core-collapse supernova explosions and explain, as best we currently can, the phases and phenomena that attend this mechanism. Two major foci of our presentation are the outer shock instability and the inner core g-mode oscillations. The former sets the stage for the latter, which damp by the generation of sound. This sound propagates outward to energize the explosion and is relevant only if the core has not exploded earlier by some other means. Hence, it is a more delayed mechanism than the traditional neutrino mechanism that has been studied for the last twenty years since it was championed by Bethe and Wilson. We discuss protoneutron star convection, accretion-induced-collapse, gravitational wave emissions, pulsar kicks, the angular anisotropy of the neutrino emissions, a subset of numerical issues, and a new code we are designing that should supercede our current supernova code VULCAN/2D. Whatever ideas last from this current generation of numerical results, and whatever the eventual mechanism(s), we conclude that the breaking of spherical symmetry will survive as one of the crucial keys to the supernova puzzle.

In the context of 2D, axisymmetric, multi-group, radiation/hydrodynamic simulations of core-collapse supernovae over the full 180$^{\circ}$ domain, we present an exploration of the progenitor dependence of the acoustic mechanism of explosion. All progenitor models we have tested with our Newtonian code explode. We investigate the roles of the Standing-Accretion-Shock-Instability (SASI), the excitation of core g-modes, the generation of core acoustic power, the ejection of matter with r-process potential, the wind-like character of the explosion, and the fundamental anisotropy of the blasts. We find that the breaking of spherical symmetry is central to the supernova phenomenon and the blasts, when top-bottom asymmetric, are self-collimating. We see indications that the initial explosion energies are larger for the more massive progenitors, and smaller for the less massive progenitors, and that the neutrino contribution to the explosion energy may be an increasing function of progenitor mass. The degree of explosion asymmetry we obtain is completely consistent with that inferred from the polarization measurements of Type Ic supernovae. Furthermore, we calculate for the first time the magnitude and sign of the net impulse on the core due to anisotropic neutrino emission and suggest that hydrodynamic and neutrino recoils in the context of our asymmetric explosions afford a natural mechanism for observed pulsar proper motions. [abridged]

We present a new theory for the gravitational wave signatures of core-collapse supernovae. Previous studies identified axisymmetric rotating core collapse, core bounce, postbounce convection, and anisotropic neutrino emission as the primary processes and phases for the radiation of gravitational waves. Our results, which are based on axisymmetric, Newtonian radiation-hydrodynamics supernova simulations (Burrows et al. 2006), indicate that the dominant emission process of gravitational waves in core-collapse supernovae may be the oscillations of the protoneutron star core. The oscillations are predominantly of g-mode character, are excited hundreds of milliseconds after bounce, and typically last for several hundred milliseconds. Our results suggest that even nonrotating core-collapse supernovae should be visible to current LIGO-class detectors throughout the Galaxy, and depending on progenitor structure, possibly out to Megaparsec distances.

We present 2.5D radiation-hydrodynamics simulations of the accretion-induced collapse (AIC) of white dwarfs, starting from 2D rotational equilibrium configurations of a 1.46-Msun and a 1.92-Msun model. Electron capture leads to the collapse to nuclear densities of these cores within a few tens of milliseconds. The shock generated at bounce moves slowly, but steadily, outwards. Within 50-100ms, the stalled shock breaks out of the white dwarf along the poles. The blast is followed by a neutrino-driven wind that develops within the white dwarf, in a cone of ~40deg opening angle about the poles, with a mass loss rate of 5-8 x 10^{-3} Msun/yr. The ejecta have an entropy on the order of 20-50 k_B/baryon, and an electron fraction distribution that is bimodal. By the end of the simulations, at >600ms after bounce, the explosion energy has reached 3-4 x 10^{49}erg and the total ejecta mass has reached a few times 0.001Msun. We estimate the asymptotic explosion energies to be lower than 10^{50}erg, significantly lower than those inferred for standard core collapse. The AIC of white dwarfs thus represents one instance where a neutrino mechanism leads undoubtedly to a successful, albeit weak, explosion. We document in detail the numerous effects of the fast rotation of the progenitors: The neutron stars are aspherical; the ``nu_mu'' and anti-nu_e neutrino luminosities are reduced compared to the nu_e neutrino luminosity; the deleptonized region has a butterfly shape; the neutrino flux and electron fraction depend strongly upon latitude (a la von Zeipel); and a quasi-Keplerian 0.1-0.5-Msun accretion disk is formed.

Based on multi-dimensional, multi-group, flux-limited-diffusion hydrodynamic simulations of core-collapse supernovae with the VULCAN/2D code, we study the physical conditions within and in the vicinity of the nascent protoneutron star (PNS). Our numerical study follows the evolution of the collapsing envelope of the 11-Msun model of Woosley & Weaver, from ~200 ms before bounce to ~300ms after bounce, on a spatial grid that switches from Cartesian at the PNS center to spherical above a 10km radius. As has been shown previously, we do not see any large-scale overturn of the inner PNS material. Convection, directly connected to the PNS, is found to occur in two distinct regions: between 10 and 20km, coincident with the region of negative lepton gradient, and exterior to the PNS above 50km. Separating these two regions, an interface, with no sizable outward or inward motion, is the site of gravity waves, emerging at ~200-300ms after core bounce, excited by the convection in the outer convective zone. In the PNS, convection is always confined within the neutrinospheric radii for all neutrino energies above just a few MeV. We find that such convective motions do not appreciably enhance the nu_e neutrino luminosity, and that they enhance the anti-nu_e and "nu_mu" luminosities modestly, by ~15% and ~30%, respectively, during the first post-bounce 100--200ms. Moreover, we see no evidence of doubly-diffusive instabilities in the PNS, expected to operate on diffusion timescales of at least a second, much longer than the millisecond timescale associated with PNS convection. PNS convection is thus found to be a secondary feature of the core-collapse phenomenon, rather than a decisive ingredient for a successful explosion.

We present results from an extensive set of one- and two-dimensional radiation-hydrodynamic simulations of the supernova core collapse, bounce, and postbounce phases, and focus on the protoneutron star (PNS) spin periods and rotational profiles as a function of initial iron core angular velocity, degree of differential rotation, and progenitor mass. For the models considered, we find a roughly linear mapping between initial iron core rotation rate and PNS spin. The results indicate that the magnitude of the precollapse iron core angular velocities is the single most important factor in determining the PNS spin. Differences in progenitor mass and degree of differential rotation lead only to small variations in the PNS rotational period and profile. Based on our calculated PNS spins, at ~ 200-300 milliseconds after bounce, and assuming angular momentum conservation, we estimate final neutron star rotation periods. We find periods of one millisecond and shorter for initial central iron core periods of below ~ 10 s. This is appreciably shorter than what previous studies have predicted and is in disagreement with current observational data from pulsar astronomy. After considering possible spindown mechanisms that could lead to longer periods we conclude that there is no mechanism that can robustly spin down a neutron star from ~ 1 ms periods to the "injection" periods of tens to hundreds of milliseconds observed for young pulsars. Our results indicate that, given current knowledge of the limitations of neutron star spindown mechanisms, precollapse iron cores must rotate with periods around 50-100 seconds to form neutron stars with periods generically near those inferred for the radio pulsar population.

We present a new mechanism for core-collapse supernova explosions that relies upon acoustic power generated in the inner core as the driver. In our simulation using an 11-solar-mass progenitor, a strong advective-acoustic oscillation a la Foglizzo with a period of ~25-30 milliseconds (ms) arises ~200 ms after bounce. Its growth saturates due to the generation of secondary shocks, and kinks in the resulting shock structure funnel and regulate subsequent accretion onto the inner core. However, this instability is not the primary agent of explosion. Rather, it is the acoustic power generated in the inner turbulent region and most importantly by the excitation and sonic damping of core g-mode oscillations. An l=1 mode with a period of ~3 ms grows to be prominent around ~500 ms after bounce. The accreting protoneutron star is a self-excited oscillator. The associated acoustic power seen in our 11-solar-mass simulation is sufficient to drive the explosion. The angular distribution of the emitted sound is fundamentally aspherical. The sound pulses radiated from the core steepen into shock waves that merge as they propagate into the outer mantle and deposit their energy and momentum with high efficiency. The core oscillation acts like a transducer to convert accretion energy into sound. An advantage of the acoustic mechanism is that acoustic power does not abate until accretion subsides, so that it is available as long as it may be needed to explode the star. [abridged]

We analyze the sensitivity of the flame propagation in a Chandrasekhar mass white dwarf to initial conditions during the subsonic burning phase, using 2D simulations of the full WD. Results are presented for a wide variety of initial flame distributions including central and off-center single point and multi-point, simultaneous and non-simultaneous, ignitions. We also examine the effects of convective velocity field which should exist at the core before the thermo-nuclear runaway. Our main conclusion suggests that the amounts of burning products and their distributions through the deflagration phase are extremely sensitive to initial conditions, much more sensitive than presented in previous studies. In particular, we find that more complex configurations such as even slight off-center ignitions, non-simultaneous multi-point ignitions and velocity fields tend to favor solutions in which individual plumes rise faster than the bulk of a typical Rayleigh-Taylor driven, unstable burning front. The difference to previous calculations for an octant of a WD may be understood as a consequence of the suppression of l=1,2 modes. Our results are consistent with full star calculations by the Chicago group. Moreover, the total amount of nuclear burning during the phase of subsonic burning depends sensitively on the initial conditions and may cause the WD to pulsate or to become unbound. We discuss the implications of the results on current models for Type Ia SNe, limitations imposed by the 2-D nature of our study, and suggest directions for further study.

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

Using the 2D multi-group, flux-limited diffusion version of the code VULCAN/2D, that also incorporates rotation, we have calculated the collapse, bounce, shock formation, and early post-bounce evolutionary phases of a core-collapse supernova for a variety of initial rotation rates. This is the first series of such multi-group calculations undertaken in supernova theory with fully multi-D tools. We find that though rotation generates pole-to-equator angular anisotropies in the neutrino radiation fields, the magnitude of the asymmetries is not as large as previously estimated. Moreover, we find that the radiation field is always more spherically symmetric than the matter distribution, with its plumes and convective eddies. We present the dependence of the angular anisotropy of the neutrino fields on neutrino species, neutrino energy, and initial rotation rate. Only for our most rapidly rotating model do we start to see qualitatively different hydrodynamics, but for the lower rates consistent with the pre-collapse rotational profiles derived in the literature the anisotropies, though interesting, are modest. This does not mean that rotation does not play a key role in supernova dynamics. The decrease in the effective gravity due to the centripetal effect can be quite important. Rather, it means that when a realistic mapping between initial and final rotational profiles and 2D multi-group radiation-hydrodynamics are incorporated into collapse simulations the anisotropy of the radiation fields may be only a secondary, not a pivotal factor, in the supernova mechanism.

We review some of the reasons for believing that the generic core-collapse supernova is neutrino-driven, not MHD-jet driven. We include a discussion of the possible role of rotation in supernova blast energetics and morphology, and speculate on the origin of Cas A's and SN1987A's ejecta fields. Two "explosive" phenomena may be associated with most core collapses, the neutrino-driven supernova itself and an underenergetic jet-like ejection that follows. The latter may be a magnetic wind that easily penetrates the debris created by the much more energetic supernova. We predict that many core-collapse supernova remnants should have sub-dominant jet-like features. In Cas A, we associate this sub-dominant collimated wind with the "jet/counter-jet" structure observed. We suggest that the actual Cas A explosion itself is at nearly right angles to this jet, along a rotation axis that coincides with the bulk of the ejecta, the iron lobes, and the putative direction of motion of the point source. It may be that when rotation becomes sufficiently rapid that the strong-neutrino-driven-supernova/weak-jet duality switches to a strong-MHD-jet scenario that might be associated with hypernovae, and in some cases GRBs. Finally, we present a calculation using a new 2D multi-group, flux-limited radiation/hydrodynamics code we have recently developed for the simulation of core-collapse supernovae. We discuss the rotation-induced anisotropy in the neutrino radiation field, neutrino heating, and the neutrino flux vectors and speculate on rotation's possible role in the supernova mechanism and in the overall supernova phenomenon.

We have developed a time-dependent, multi-energy-group, and multi-angle (S$_n$) Boltzmann transport scheme for radiation hydrodynamics simulations, in one and two spatial dimensions. The implicit transport is coupled to both 1D (spherically-symmetric) and 2D (axially-symmetric) versions of the explicit Newtonian hydrodynamics code VULCAN. The 2D variant, VULCAN/2D, can be operated in general structured or unstructured grids and though the code can address many problems in astrophysics it was constructed specifically to study the core-collapse supernova problem. Furthermore, VULCAN/2D can simulate the radiation/hydrodynamic evolution of differentially rotating bodies. We summarize the equations solved and methods incorporated into the algorithm and present results of a time-dependent 2D test calculation. A more complete description of the algorithm is postponed to another paper. We highlight a 2D test run that follows for 22 milliseconds the immediate post-bounce evolution of a collapsed core. We present the relationship between the anisotropies of the overturning matter field and the distribution of the corresponding flux vectors, as a function of energy group. This is the first 2D multi-group, multi-angle, time-dependent radiation/hydro calculation ever performed in core collapse studies. Though the transport module of the code is not gray and does not use flux limiters (however, there is a flux-limited variant of VULCAN/2D), it still does not include energy redistribution and most velocity-dependent terms.

We have carried out an extensive set of two-dimensional, axisymmetric, purely-hydrodynamic calculations of rotational stellar core collapse with a realistic, finite-temperature nuclear equation of state and realistic massive star progenitor models. For each of the total number of 72 different simulations we performed, the gravitational wave signature was extracted via the quadrupole formula in the slow-motion, weak-field approximation. We investigate the consequences of variation in the initial ratio of rotational kinetic energy to gravitational potential energy and in the initial degree of differential rotation. Furthermore, we include in our model suite progenitors from recent evolutionary calculations that take into account the effects of rotation and magnetic torques. For each model, we calculate gravitational radiation wave forms, characteristic wave strain spectra, energy spectra, final rotational profiles, and total radiated energy. In addition, we compare our model signals with the anticipated sensitivities of the 1st- and 2nd-generation LIGO detectors coming on line. We find that most of our models are detectable by LIGO from anywhere in the Milky Way.

Delayed detonation in an exploding white dwarf, which propagates from an off-center transition point, rather than from a spherical transition shell, is described and simulated. The differences between the results of 2D simulations and the 1D case are presented and discussed. The two dimensional effects become significant in transition density below 3.e7 g/cm^3, where the energetics, the production of Fe group elements and the symmetry of the explosion are all affected. In the 2D case the explosion is less energetic and less Ni is produced in the detonation phase of the explosion. For low transition density the reduction in Ni mass can reach 20-30 percent. The asymmetry in abundances between regions close to the transition point and regions far from that point is large, and could be a source to polarization patterns in the emitted light. We conclude that the spatial and temporal distribution of transition locations, is an important parameter which must be included in delayed detonation models for Type Ia supernovae. \