A first-order quantum chromodynamics (QCD) phase transition (PT) may take place in the protocompact star (PCS) produced by a core-collapse supernova (CCSN). In this work, we study the consequences of such a PT in a non-rotating CCSN with axisymmetric hydrodynamic simulations. We find that the PT leads to the collapse of the PCS and results in a loud burst of gravitational waves (GWs). The amplitude of this GW burst is $\sim30$ times larger than the post-bounce GW signal normally found for non-rotating CCSN. It shows a broad peak at high frequencies ($\sim2500-4000$ Hz) in the spectrum, has a duration of $\lesssim5 {\rm ms}$, and carries $\sim3$ orders of magnitude more energy than the other episodes. Also, the peak frequency of the PCS oscillation increases dramatically after the PT-induced collapse. In addition to a second neutrino burst, the GW signal, if detected by the ground-based GW detectors, is decisive evidence of the first-order QCD PT inside CCSNe and provides key information about the structure and dynamics of the PCS.

Radioactive isotopes produced in core-collapse supernovae (CCSNe) provide useful insights into the underlying processes driving the collapse mechanism and the origins of elemental abundances. Their study generates a confluence of major physics research, including experimental measurements of nuclear reaction rates, astrophysical modeling, and $\gamma$-ray observations. Here we identify the key nuclear reaction rates to the nucleosynthesis of observable radioactive isotopes in explosive silicon-burning during CCSNe. Using the nuclear reaction network calculator SkyNet and current REACLIB reaction rates, we evolve temperature-density-time profiles of the innermost $0.45~M_\odot$ ejecta from the core collapse and explosion of a $12~M_\odot$ star. Individually varying 3403 reaction rates by factors of 100, we identify 141 reactions which cause significant differences in the isotopes of interest, namely, $^{43}$K, $^{47}$Ca, $^{44,47}$Sc, $^{44}$Ti, $^{48,51}$Cr, $^{48,49}$V, $^{52,53}$Mn, $^{55,59}$Fe, $^{56,57}$Co, and $^{56,57,59}$Ni. For each of these reactions, we present a novel method to extract the temperature range pertinent to the nucleosynthesis of the relevant isotope; the resulting temperatures lie within the range $T = 0.47$ to $6.15~$GK. Limiting the variations to within $1\sigma$ of STARLIB reaction rate uncertainties further reduces the identified reactions to 48 key rates, which can be used to guide future experimental research. Complete results are presented in tabular form.

With the advent of modern neutrino and gravitational wave detectors, the promise of multi-messenger detections of the next galactic core-collapse supernova has become very real. Such detections will give insight into the core-collapse supernova mechanism, the structure of the progenitor star, and may resolve longstanding questions in fundamental physics. In order to properly interpret these detections, a thorough understanding of the landscape of possible core-collapse supernova events, and their multi-messenger signals, is needed. We present detailed predictions of neutrino and gravitational wave signals from 1D simulations of stellar core collapse, spanning the landscape of core-collapse progenitors from $9-120\,\mathrm{M}_{\odot}$. In order to achieve explosions in 1D, we use the STIR model, which includes the effects of turbulence and convection in 1D supernova simulations to mimic the 3D explosion mechanism. We study the gravitational wave emission from the 1D simulations using an astroseismology analysis of the proto-neutron star. We find that the neutrino and gravitational wave signals are strongly correlated with the structure of the progenitor star and remnant compact object. Using these correlations, future detections of the first few seconds of neutrino and gravitational wave emission from a galactic core-collapse supernova may be able to provide constraints on stellar evolution independent of pre-explosion imaging and the mass of the compact object remnant prior to fallback accretion.

The core collapse of a massive star results in the formation of a proto-neutron star (PNS). If enough material is accreted onto a PNS it will become gravitationally unstable and further collapse into a black-hole (BH). We perform a systematic study of failing core-collapse supernovae in spherical symmetry for a wide range of presupernova progenitor stars and equations of state (EOSs) of nuclear matter. We analyze how variations in progenitor structure and the EOS of dense matter above nuclear saturation density affect the PNS evolution and subsequent BH formation. Comparisons of core-collapse for a given progenitor star and different EOSs show that the path traced by the PNS in mass-entropy phase space $M_{\mathrm{grav}}^{\mathrm{PNS}}-\tilde{s}$ is well correlated with the progenitor compactness and almost EOS independent, apart from the final endpoint. Furthermore, BH formation occurs, to a very good approximation, soon after the PNS overcomes the maximum \textit{gravitational} mass supported by a hot NS with constant entropy equal to $\tilde{s}$. These results show a path to constraining the temperature dependence of the EOS through the detection of neutrinos from a failed galactic supernova.

The core-collapse supernova (CCSN) mechanism is fundamentally three-dimensional with instabilities, convection, and turbulence playing crucial roles in aiding neutrino-driven explosions. Simulations of CCNSe including accurate treatments of neutrino transport and sufficient resolution to capture key instabilities remain amongst the most expensive numerical simulations in astrophysics, prohibiting large parameter studies in 2D and 3D. Studies spanning a large swath of the incredibly varied initial conditions of CCSNe are possible in 1D, though such simulations must be artificially driven to explode. We present a new method for including the most important effects of convection and turbulence in 1D simulations of neutrino-driven CCSNe, called Supernova Turbulence In Reduced-dimensionality, or STIR. Our new approach includes crucial terms resulting from the turbulent and convective motions of the flow. We estimate the strength of convection and turbulence using a modified mixing length theory (MLT) approach introducing a few free parameters to the model which are fit to the results of 3D simulations. For sufficiently large values of the mixing length parameter, turbulence-aided neutrino-driven explosions are obtained. We compare the results of STIR to high-fidelity 3D simulations and perform a parameter study of CCSN explosion using 200 solar-metallicity progenitor models from 9 to 120 $M_\odot$. We find that STIR is a better predictor of which models will explode in multidimensional simulations than other methods of driving explosions in 1D. We also present a preliminary investigation of predicted observable characteristics of the CCSN population from STIR, such as the distributions of explosion energies and remnant masses.

We investigate correlated gravitational wave and neutrino signals from rotating core-collapse supernovae with simulations. Using an improved mode identification procedure based on mode function matching, we show that a linear quadrupolar mode of the core produces a dual imprint on gravitational waves and neutrinos in the early post-bounce phase of the supernova. The angular harmonics of the neutrino emission are consistent with the mode energy around the neutrinospheres, which points to a mechanism for the imprint on neutrinos. Thus, neutrinos carry information about the mode amplitude in the outer region of the core, whereas gravitational waves probe deeper in. We also find that the best-fit mode function has a frequency bounded above by $\sim 420$ Hz, and yet the mode's frequency in our simulations is $\sim 15\%$ higher, due to the use of Newtonian hydrodynamics and a widely used pseudo-Newtonian gravity approximation. This overestimation is particularly important for the analysis of gravitational wave detectability and asteroseismology, pointing to limitations of pseudo-Newtonian approaches for these purposes, possibly even resulting in excitation of incorrect modes. In addition, mode frequency matching (as opposed to mode function matching) could be resulting in mode misidentification in recent work. Lastly, we evaluate the prospects of a multimessenger detection of the mode using current technology. The detection of the imprint on neutrinos is most challenging, with a maximum detection distance of $\sim\!1$ kpc using the IceCube Neutrino Observatory. The maximum distance for detecting the complementary gravitational wave imprint is $\sim\!5$ kpc using Advanced LIGO at design sensitivity.

We explore the influence of progenitor mass and rotation on the gravitational-wave (GW) emission from core-collapse supernovae, during the postbounce, preexplosion, accretion-phase. We present the results from 15 two-dimensional (2D) neutrino radiation-hydrodynamic simulations from initial stellar collapse to $\sim$300 ms after core bounce. We examine the features of the GW signals for four zero-age main sequence (ZAMS) progenitor masses ranging from 12 $M_\odot$ to 60 $M_\odot$ and four core rotation rates from 0 to 3 rad s$^{-1}$. We find that GW strain immediately around core bounce is fairly independent of ZAMS mass and---consistent with previous findings---that it is more heavily dependent on the core angular momentum. At later times, all nonrotating progenitors exhibit loud GW emission, which we attribute to vibrational g-modes of the protoneutron star (PNS) excited by convection in the postshock layer and the standing accretion shock instability (SASI). We find that increasing rotation rates results in muting of the accretion-phase GW signal due to centrifugal effects that inhibit convection in the postshock region, quench the SASI, and slow the rate at which the PNS peak vibrational frequency increases. Additionally, we verify the efficacy of our approximate general relativistic (GR) effective potential treatment of gravity by comparing our core bounce GW strains with the recent 2D GR results of other groups.

The identification of the electromagnetic counterpart of the gravitational wave event, GW170817, and discovery of neutrinos and gamma-rays from TXS 0506+056 heralded the new era of multi-messenger astrophysics. As the number of multi-messenger events rapidly grow over the next decade, the cyberinfrastructure requirements to handle the increase in data rates, data volume, need for event follow up, and analysis across the different messengers will also explosively grow. The cyberinfrastructure requirements to enhance multi-messenger astrophysics will both be a major challenge and opportunity for astronomers, physicists, computer scientists and cyberinfrastructure specialists. Here we outline some of these requirements and argue for a distributed cyberinfrastructure institute for multi-messenger astrophysics to meet these challenges.

As of today, we have directly detected exactly one source in both gravitational waves (GWs) and electromagnetic (EM) radiation, the binary neutron star merger GW170817, its associated gamma-ray burst GRB170817A, and the subsequent kilonova SSS17a/AT 2017gfo. Within ten years, we will detect hundreds of events, including new classes of events such as neutron-star-black-hole mergers, core-collapse supernovae, and almost certainly something completely unexpected. As we build this sample, we will explore exotic astrophysical topics ranging from nucleosynthesis, stellar evolution, general relativity, high-energy astrophysics, nuclear matter, to cosmology. The discovery potential is extraordinary, and investments in this area will yield major scientific breakthroughs. Here we outline some of the most exciting scientific questions that can be answered by combining GW and EM observations.

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.

Shock revival in core-collapse supernovae (CCSNe) may be due to the neutrino mechanism. While it is known that in a neutrino-powered CCSN, explosion begins when the neutrino luminosity of the proto-neutron star exceeds a critical value, the physics of this condition in time-dependent, multidimensional simulations are not fully understood. \citet{Pejcha2012} found that an `antesonic condition' exists for time-steady spherically symmetric models, potentially giving a physical explanation for the critical curve observed in simulations. In this paper, we extend that analysis to time-dependent, spherically symmetric polytropic models. We verify the critical antesonic condition in our simulations, showing that models exceeding it drive transonic winds whereas models below it exhibit steady accretion. In addition, we find that (1) high spatial resolution is needed for accurate determination of the antesonic ratio and shock radius at the critical curve, and that low resolution simulations systematically underpredict these quantities, making explosion more difficult at lower resolution; (2) there is an important physical connection between the critical mass accretion rate at explosion and the mass loss rate of the post-explosion wind: the two are directly proportional at criticality, implying that, at criticality, the wind kinetic power is tied directly to the accretion power; (3) the value of the post-shock adiabatic index $\Gamma$ has a large effect on the length and time scales of the post-bounce evolution of the explosion larger values of $\Gamma$ result in a longer transition from the accretion to wind phases.

Neutrinos play a crucial role in the core-collapse supernova (CCSN) explosion mechanism. The requirement of accurately calculating the transport of neutrinos makes simulations of the CCSN mechanism extremely challenging and computationally expensive. Historically, this stiff challenge has been met by making approximations to the full transport equation. In this work, we compare CCSN simulations in one- and two-dimensions with three approximate neutrino transport schemes, each implemented in the {\tt FLASH} simulation framework. We compare a two moment M1 scheme with an analytic closure (M1), the Isotropic Diffusion Source Approximation (IDSA), and the Advanced Spectral Leakage (ASL) method. We identify and discuss the advantages and disadvantages of each scheme. For each approximate transport scheme, we use identical grid setups, hydrodynamics, and gravity solvers to investigate the transport effects on supernova shock dynamics and neutrino quantities. We find that the transport scheme has a small effect on the evolution of protoneutron star (PNS) radius, PNS mass, and the mass accretion rate. The neutrino luminosities, mean energies, and shock radii have a $\sim$ 10-20\% quantitative difference but the overall qualitative trends are fairly consistent between all three approximations. We find larger differences in the gain region properties, including the gain region mass and the net heating rate in the gain region, as well as the strength of PNS convection in the core. We investigate the progenitor, nuclear equation of state, and stochastic perturbation dependence of our simulations and find similar magnitudes of impact on key quantities. We also compare the computational expense of the various approximations.

In July 2018 an FRIB Theory Alliance program was held on the implications of GW170817 and its associated kilonova for r-process nucleosynthesis. Topics of discussion included the astrophysical and nuclear physics uncertainties in the interpretation of the GW170817 kilonova, what we can learn about the astrophysical site or sites of the r process from this event, and the advances in nuclear experiment and theory most crucial to pursue in light of the new data. Here we compile a selection of scientific contributions to the workshop, broadly representative of progress in r-process studies since the GW170817 event.

The details of the physical mechanism that drives core-collapse supernovae (CCSNe) remain uncertain. While there is an emerging consensus on the qualitative outcome of detailed CCSN mechanism simulations in 2D, only recently have high-fidelity 3D simulations become possible. Here we present the results of an extensive set of 3D CCSN simulations using high-fidelity multidimensional neutrino transport, high-resolution hydrodynamics, and approximate general relativistic gravity. We employ a state-of-the-art 20 solar mass progenitor generated using the Modules for Experiments in Stellar Astrophysics (MESA; Farmer et al. (2016) Paxton et al. (2011, 2013, 2015, 2018) and the SFHo equation of state of Steiner et al. (2013). While none of our 3D CCSN simulations explode within ~500ms after core bounce, we find that the presence of large scale aspherical motion in the Si and O shells surrounding the collapsing iron core aid shock expansion and bring the models closer to the threshold of explosion. We also find some dependence on resolution and geometry (octant vs. full 4$\pi$). As has been noted in other recent works, we find that the post-shock turbulence plays an important role in determining the overall dynamical evolution of our simulations. We find a strong standing accretion shock instability (SASI) that develops at late times during the shock recession epoch. The SASI aids in transient shock expansion phases, but is not enough to result in shock revival. We also report that for a subset of our simulations, we find conclusive evidence for the LESA first reported in Tamborra et al. (2014), but until now, not confirmed by other simulation codes. Both the progenitor asphericities and the SASI-induced transient shock expansion phases generate transient gravitational waves and neutrino signal modulations via perturbations of the protoneutron star by turbulent motions. (abridged)

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

We investigate axisymmetric black hole~(BH) formation and its gravitational wave (GW) and neutrino signals with self-consistent core-collapse supernova simulations of a non-rotating $40~M_\odot$ progenitor star using the isotropic diffusion source approximation for the neutrino transport and a modified gravitational potential for general relativistic effects. We consider four different neutron star (NS) equations of state~(EoS): LS220, SFHo, BHB$\Lambda\phi$ and DD2, and study the impact of the EoS on BH formation dynamics and GW emission. We find that the BH formation time is sensitive to the EoS from 460 to $>1300$~ms and is delayed in multiple dimensions for $\sim~100-250$~ms due to the finite entropy effects. Depending on the EoS, our simulations show the possibility that shock revival can occur along with the collapse of the proto-neutron star~(PNS) to a BH. The gravitational waveforms contain four major features that are similar to previous studies but show extreme values: (1)~a low frequency signal ($\sim~300-500$~Hz) from core-bounce and prompt convection, (2)~a strong signal from the PNS g-mode oscillation among other features, (3)~a high frequency signal from the PNS inner-core convection, and (4)~signals from the standing accretion shock instability and convection. The peak frequency at the onset of BH formation reaches to $\sim~2.3$~kHz. The characteristic amplitude of a 10~kpc object at peak frequency is detectable but close to the noise threshold of the Advanced~LIGO and KAGRA, suggesting that the next generation gravitational wave detector will need to improve the sensitivity at the kHz domain to better observe stellar-mass BH formation from core-collapse supernovae or failed supernovae.

We present results from simulations of core-collapse supernovae in FLASH using a newly-implemented multidimensional neutrino transport scheme and a newly-implemented general relativistic (GR) treatment of gravity. We use a two-moment method with an analytic closure (so-called M1 transport) for the neutrino transport. This transport is multienergy, multispecies, velocity-dependent and truly multidimensional, i.e., we do not assume the commonly used "ray-by-ray" approximation. Our GR gravity is implemented in our Newtonian hydrodynamics simulations via an effective relativistic potential (GREP) that closely reproduces the GR structure of neutron stars and has been shown to match GR simulations of core collapse quite well. In axisymmetry, we simulate core-collapse supernovae with four different progenitor models in both Newtonian and GR gravity. We find that the more compact protoneutron star structure realized in simulations with GR gravity gives higher neutrino luminosities and higher neutrino energies. These differences in turn give higher neutrino heating rates (upwards of $\sim$20-30% over the corresponding Newtonian gravity simulations) that increase the efficacy of the neutrino mechanism. Three of the four models successfully explode in the simulations assuming GREP gravity. In our Newtonian gravity simulations, two of the four models explode, but at times much later than observed in our GR gravity simulations. Our results, both in Newtonian and GR gravity, compare well with several other studies in the literature. These results conclusively show that the approximation of Newtonian gravity for simulating the core-collapse supernova central engine is not acceptable. We also simulate four additional models in GR gravity to highlight the growing disparity between parameterized 1D models of core-collapse supernovae and the current generation of 2D models.

Multidimensional simulations show that non-radial, turbulent, fluid motion is a fundamental component of the core-collapse supernova (CCSN) explosion mechanism. Neutrino-driven convection, the standing accretion shock instability, and relic-perturbations from advanced stages of nuclear burning can all impact the outcome of core collapse in a qualitative and quantitative way. Here, we review the current understanding of these phenomena and their role in the explosion of massive stars. We also discuss the role of protoneutron star convection and of magnetic fields in the context of the delayed neutrino mechanism.

Core-collapse supernova explosions are driven by a central engine that converts a small fraction of the gravitational binding energy released during core collapse to outgoing kinetic energy. The suspected mode for this energy conversion is the neutrino mechanism, where a fraction of the neutrinos emitted from the newly formed protoneutron star are absorbed by and heat the matter behind the supernova shock. Accurate neutrino-matter interaction terms are crucial for simulating these explosions. In this proceedings for IAUS 331, SN 1987A, 30 years later, we explore several corrections to the neutrino-nucleon scattering opacity and demonstrate the effect on the dynamics of the core-collapse supernova central engine via two dimensional neutrino-radiation-hydrodynamics simulations. Our results reveal that the explosion properties are sensitive to corrections to the neutral-current scattering cross section at the 10-20% level, but only for densities at or above $\sim 10^{12}$ g cm$^{-3}$

We explore properties of core-collapse supernova progenitors with respect to the composite uncertainties in the thermonuclear reaction rates by coupling the reaction rate probability density functions provided by the STARLIB reaction rate library with $\texttt{MESA}$ stellar models. We evolve 1000 15 $M_{\odot}$ models from the pre main-sequence to core O-depletion at solar and subsolar metallicities for a total of 2000 Monte Carlo stellar models. For each stellar model, we independently and simultaneously sample 665 thermonuclear reaction rates and use them in a $\texttt{MESA}$ in situ reaction network that follows 127 isotopes from $^{1}$H to $^{64}$Zn. With this framework we survey the core mass, burning lifetime, composition, and structural properties at five different evolutionary epochs. At each epoch we measure the probability distribution function of the variations of each property and calculate Spearman Rank-Order Correlation coefficients for each sampled reaction rate to identify which reaction rate has the largest impact on the variations on each property. We find that uncertainties in $^{14}$N$(p,\gamma)^{15}$O, triple-$\alpha$, $^{12}$C$(\alpha,\gamma)^{16}$O, $^{12}$C($^{12}$C,$p$)$^{23}$Na, $^{12}$C($^{16}$O,$p$)$^{27}$Al, $^{16}$O($^{16}$O,$n$)$^{31}$S, $^{16}$O($^{16}$O,$p$)$^{31}$P, and $^{16}$O($^{16}$O,$\alpha$)$^{28}$Si reaction rates dominate the variations of the properties surveyed. We find that variations induced by uncertainties in nuclear reaction rates grow with each passing phase of evolution, and at core H-, He-depletion are of comparable magnitude to the variations induced by choices of mass resolution and network resolution. However, at core C-, Ne-, and O-depletion, the reaction rate uncertainties can dominate the variation causing uncertainty in various properties of the stellar model in the evolution towards iron core-collapse.

We present results from high-resolution semi-global simulations of neutrino-driven convection in core-collapse supernovae. We employ an idealized setup with parametrized neutrino heating/cooling and nuclear dissociation at the shock front. We study the internal dynamics of neutrino-driven convection and its role in re-distributing energy and momentum through the gain region. We find that even if buoyant plumes are able to locally transfer heat up to the shock, convection is not able to create a net positive energy flux and overcome the downwards transport of energy from the accretion flow. Turbulent convection does, however, provide a significant effective pressure support to the accretion flow as it favors the accumulation of energy, mass and momentum in the gain region. We derive an approximate equation that is able to explain and predict the shock evolution in terms of integrals of quantities such as the turbulent pressure in the gain region or the effects of non-radial motion of the fluid. We use this relation as a way to quantify the role of turbulence in the dynamics of the accretion shock. Finally, we investigate the effects of grid resolution, which we change by a factor 20 between the lowest and highest resolution. Our results show that the shallow slopes of the turbulent kinetic energy spectra reported in previous studies are a numerical artefact. Kolmogorov scaling is progressively recovered as the resolution is increased.

We explore the effects of rotation on convective carbon, oxygen, and silicon shell burning during the late stages of evolution in a 20Msun star. Using the Modules for Experiments in Stellar Astrophysics (MESA) we construct 1D stellar models both with no rotation and with an initial rigid rotation of 50% of critical. At different points during the evolution, we map the 1D models into 2D and follow the multidimensional evolution using the FLASH compressible hydrodynamics code for many convective turnover times until a quasi-steady state is reached. We characterize the strength and scale of convective motions via decomposition of the momentum density into vector spherical harmonics. We find that rotation influences the total power in solenoidal modes, with a slightly larger impact for carbon and oxygen shell burning than for silicon shell burning. Including rotation in one-dimensional (1D) stellar evolution models alters the structure of the star in a manner that has a significant impact on the character of multidimensional convection. Adding modest amounts of rotation to a stellar model that ignores rotation during the evolutionary stage, however, has little impact on the character of resulting convection. Since the spatial scale and strength of convection present at the point of core collapse directly influence the supernova mechanism, our results suggest that rotation could play an important role in setting the stage for massive stellar explosions.

We present the SuperNova Explosion Code (SNEC), an open-source Lagrangian code for the hydrodynamics and equilibrium-diffusion radiation transport in the expanding envelopes of supernovae. Given a model of a progenitor star, an explosion energy, and an amount and distribution of radioactive nickel, SNEC generates the bolometric light curve, as well as the light curves in different broad bands assuming black body emission. As a first application of SNEC, we consider the explosions of a grid of 15 Msun (at zero-age main sequence) stars whose hydrogen envelopes are stripped to different extents and at different points in their evolution. The resulting light curves exhibit plateaus with durations of ~20-100 days if >~1.5-2 Msun of hydrogen-rich material is left and no plateau if less hydrogen-rich material is left. If these shorter plateau lengths are not seen for Type IIP supernovae in nature, it suggests that, at least for zero-age main sequence masses <~ 20 Msun, hydrogen mass loss occurs as an all or nothing process. This perhaps points to the important role binary interactions play in generating the observed mass-stripped supernovae (i.e., Type Ib/c events). These light curves are also unlike what is typically seen for Type IIL supernovae, arguing that simply varying the amount of mass loss cannot explain these events. The most stripped models begin to show double-peaked light curves similar to what is often seen for Type IIb supernovae, confirming previous work that these supernovae can come from progenitors that have a small amount of hydrogen and a radius of ~500 Rsun.

(Abridged) In the implicit large eddy simulation (ILES) paradigm, the dissipative nature of high-resolution shock-capturing schemes is exploited to provide an implicit model of turbulence. Recent 3D simulations suggest that turbulence might play a crucial role in core-collapse supernova explosions, however the fidelity with which turbulence is simulated in these studies is unclear. Especially considering that the accuracy of ILES for the regime of interest in CCSN, weakly compressible and strongly anisotropic, has not been systematically assessed before. In this paper we assess the accuracy of ILES using numerical methods most commonly employed in computational astrophysics by means of a number of local simulations of driven, weakly compressible, anisotropic turbulence. We report a detailed analysis of the way in which the turbulent cascade is influenced by the numerics. Our results suggest that anisotropy and compressibility in CCSN turbulence have little effect on the turbulent kinetic energy spectrum and a Kolmogorov $k^{-5/3}$ scaling is obtained in the inertial range. We find that, on the one hand, the kinetic energy dissipation rate at large scales is correctly captured even at relatively low resolutions, suggesting that very high effective Reynolds number can be achieved at the largest scales of the simulation. On the other hand, the dynamics at intermediate scales appears to be completely dominated by the so-called bottleneck effect, \ie the pile up of kinetic energy close to the dissipation range due to the partial suppression of the energy cascade by numerical viscosity. An inertial range is not recovered until the point where relatively high resolution $\sim 512^3$, which would be difficult to realize in global simulations, is reached. We discuss the consequences for CCSN simulations.

We present the first three dimensional (3D) simulation of the final minutes of iron core growth in a massive star, up to and including the point of core gravitational instability and collapse. We self-consistently capture the development of strong convection driven by violent Si burning in the shell surrounding the iron core. This convective burning builds the iron core to its critical (Chandrasekhar) mass and collapse ensues, driven by electron capture and photodisintegration. The non-spherical structure and motion (turbulent fluctuations) generated by 3D convection is substantial at the point of collapse. We examine the impact of such physically-realistic 3D initial conditions on the core-collapse supernova mechanism using 3D simulations including multispecies neutrino leakage. We conclude that non-spherical progenitor structure should not be ignored, and has a significant and favorable impact on the likelihood for neutrino-driven explosions.

The neutrino-heated "gain layer" immediately behind the stalled shock in a core-collapse supernova is unstable to high-Reynolds-number turbulent convection. We carry out and analyze a new set of 19 high-resolution three-dimensional (3D) simulations with a three-species neutrino leakage/heating scheme and compare with spherically-symmetric (1D) and axisymmetric (2D) simulations carried out with the same methods. We study the postbounce supernova evolution in a $15$-$M_\odot$ progenitor star and vary the local neutrino heating rate, the magnitude and spatial dependence of asphericity from convective burning in the Si/O shell, and spatial resolution. Our simulations suggest that there is a direct correlation between the strength of turbulence in the gain layer and the susceptability to explosion. 2D and 3D simulations explode at much lower neutrino heating rates than 1D simulations. This is commonly explained by the fact that nonradial dynamics allows accreting material to stay longer in the gain layer. We show that this explanation is incomplete. Our results indicate that the effective turbulent ram pressure exerted on the shock plays a crucial role by allowing multi-D models to explode at a lower postshock thermal pressure and thus with less neutrino heating than 1D models. We connect the turbulent ram pressure with turbulent energy at large scales and in this way explain why 2D simulations are erroneously exploding more easily than 3D simulations.

We apply the mathematical formalism of vector spherical harmonics decomposition to convective stellar velocity fields from multi-dimensional hydrodynamics simulations, and show that the resulting power spectra furnish a robust and stable statistical description of stellar convective turbulence. Analysis of the power spectra help identify key physical parameters of the convective process such as the dominant scale of the turbulent motions that influence the structure of massive evolved pre-supernova stars. We introduce the numerical method that can be used to calculate vector spherical harmonics power spectra from 2D and 3D convective shell simulation data. Using this method we study the properties of oxygen shell burning and convection for a 15 Msun star simulated by the hydrodynamics code FLASH in 2D and 3D. We discuss the importance of realistic initial conditions to achieving successful core-collapse supernova explosions in multi-dimensional simulations. We show that the calculated power spectra can be used to generate realizations of the velocity fields of pre-supernova convective shells. We find that the slope of the solenoidal mode power spectrum remains mostly constant throughout the evolution of convection in the oxygen shell in both 2D and 3D simulations. We also find that the characteristic radial scales of the convective elements are smaller in 3D than in 2D while the angular scales are larger in 3D.

Three-dimensional simulations of core-collapse supernovae are granting new insight into the as-yet uncertain mechanism that drives successful explosions. While there is still debate about whether explosions are obtained more easily in 3D than in 2D, it is undeniable that there exist qualitative and quantitative differences between the results of 3D and 2D simulations. We present an extensive set of high-resolution one-, two-, and three-dimensional core-collapse supernova simulations with multispecies neutrino leakage carried out in two different progenitors. Our simulations confirm the results of Couch (2013) indicating that 2D explodes more readily than 3D. We argue that this is due to the inadequacies of 2D to accurately capture important aspects of the three-dimensional dynamics. We find that without artificially enhancing the neutrino heating rate we do not obtain explosions in 3D. We examine the development of neutrino-driven convection and the standing accretion shock instability and find that, in separate regimes, either instability can dominate. We find evidence for growth of the standing accretion shock instability for both 15-$M_\odot$ and 27-$M_\odot$ progenitors, however, it is weaker in 3D exploding models. The growth rate of both instabilities is artificially enhanced along the symmetry axis in 2D as compared with our axis-free 3D Cartesian simulations. Our work highlights the growing consensus that core-collapse supernovae must be studied in 3D if we hope to solve the mystery of how the explosions are powered.

We explore Implicit Monte Carlo (IMC) and Discrete Diffusion Monte Carlo (DDMC) for radiation transport in high-velocity outflows with structured opacity. The IMC method is a stochastic computational technique for nonlinear radiation transport. IMC is partially implicit in time and may suffer in efficiency when tracking Monte Carlo particles through optically thick materials. DDMC accelerates IMC in diffusive domains. Abdikamalov extended IMC and DDMC to multigroup, velocity-dependent transport with the intent of modeling neutrino dynamics in core-collapse supernovae. Densmore has also formulated a multifrequency extension to the originally grey DDMC method. We rigorously formulate IMC and DDMC over a high-velocity Lagrangian grid for possible application to photon transport in the post-explosion phase of Type Ia supernovae. This formulation includes an analysis that yields an additional factor in the standard IMC-to-DDMC spatial interface condition. To our knowledge the new boundary condition is distinct from others presented in prior DDMC literature. The method is suitable for a variety of opacity distributions and may be applied to semi-relativistic radiation transport in simple fluids and geometries. Additionally, we test the code, called SuperNu, using an analytic solution having static material, as well as with a manufactured solution for moving material with structured opacities. Finally, we demonstrate with a simple source and 10 group logarithmic wavelength grid that IMC-DDMC performs better than pure IMC in terms of accuracy and speed when there are large disparities between the magnitudes of opacities in adjacent groups. We also present and test our implementation of the new boundary condition.

Multi-dimensional simulations of advanced nuclear burning stages of massive stars suggest that the Si/O layers of presupernova stars harbor large deviations from the spherical symmetry typically assumed for presupernova stellar structure. We carry out three-dimensional core-collapse supernova simulations with and without aspherical velocity perturbations to assess their potential impact on the supernova hydrodynamics in the stalled shock phase. Our results show that realistic perturbations can qualitatively alter the postbounce evolution, triggering an explosion in a model that fails to explode without them. This finding underlines the need for a multi-dimensional treatment of the presupernova stage of stellar evolution.

Self-gravity computation by multipole expansion is a common approach in problems such as core-collapse and Type Ia supernovae, where single large condensations of mass must be treated. The standard formulation of multipole self-gravity suffers from two significant sources of error, which we correct in the formulation presented in this article. The first source of error is due to the numerical approximation that effectively places grid cell mass at the central point of the cell, then computes the gravitational potential at that point, resulting in a convergence failure of the multipole expansion. We describe a new scheme that avoids this problem by computing gravitational potential at cell faces. The second source of error is due to sub-optimal choice of location for the expansion center, which results in angular power at high multipole $l$ values in the gravitational field, requiring a high --- and expensive --- value of multipole cutoff \lmax. By introducing a global measure of angular power in the gravitational field, we show that the optimal coordinate for the expansion is the square-density-weighted mean location. We subject our new multipole self-gravity algorithm to two rigorous test problems: MacLaurin spheroids for which exact analytic solutions are known, and core-collapse supernovae. We show that key observables of the core-collapse simulations, particularly shock expansion, proto-neutron star motion, and momentum conservation, are extremely sensitive to the accuracy of the multipole gravity, and the accuracy of their computation is greatly improved by our reformulated solver.

We study the effects of rotation on the dynamics, energetics and Ni-56 production of Pair Instability Supernova explosions by performing rotating two-dimensional ("2.5-D") hydrodynamics simulations. We calculate the evolution of eight low metallicity (Z = 10^-3, 10^-4 Zsun) massive (135-245 Msun) PISN progenitors with initial surface rotational velocities 50% that of the critical Keplerian value using the stellar evolution code MESA. We allow for both the inclusion and the omission of the effects of magnetic fields in the angular momentum transport and in chemical mixing, resulting in slowly-rotating and rapidly-rotating final carbon-oxygen cores, respectively. Increased rotation for carbon-oxygen cores of the same mass and chemical stratification leads to less energetic PISN explosions that produce smaller amounts of Ni-56 due to the effect of the angular momentum barrier that develops and slows the dynamical collapse. We find a non-monotonic dependence of Ni-56 production on rotational velocity in situations when smoother composition gradients form at the outer edge of the rotating cores. In these cases, the PISN energetics are determined by the competition of two factors: the extent of chemical mixing in the outer layers of the core due to the effects of rotation in the progenitor evolution and the development of angular momentum support against collapse. Our 2.5-D PISN simulations with rotation are the first presented in the literature. They reveal hydrodynamic instabilities in several regions of the exploding star and increased explosion asymmetries with higher core rotational velocity.

We present 1D, 2D, and 3D hydrodynamical simulations of core-collapse supernovae including a parameterized neutrino heating and cooling scheme in order to investigate the critical core neutrino luminosity (L_crit) required for explosion. In contrast to some previous works, we find that 3D simulations explode later than 2D simulations, and that L_crit at fixed mass accretion rate is somewhat higher in 3D than in 2D. We find, however, that in 2D L_crit increases as the numerical resolution of the simulation increases. In contrast to some previous works, we argue that the average entropy of the gain region is in fact not a good indicator of explosion but is rather a reflection of the greater mass in the gain region in 2D. We compare our simulations to semi-analytic explosion criteria and examine the nature of the convective motions in 2D and 3D. We discuss the balance between neutrino-driven-buoyancy and drag forces. In particular, we show that the drag force will be proportional to a buoyant plume's surface area while the buoyant force is proportional to a plume's volume and, therefore, plumes with greater volume-to-surface area ratios will rise more quickly. We show that buoyant plumes in 2D are inherently larger, with greater volume-to-surface area ratios, than plumes in 3D. In the scenario that the supernova shock expansion is dominated by neutrino-driven buoyancy, this balance between buoyancy and drag forces may explain why 3D simulations explode later than 2D simulations and why L_crit increases with resolution. Finally, we provide a comparison of our results with other calculations in the literature.

We study the dependence of the delayed neutrino-heating mechanism for core-collapse supernovae on the equation of state. Using a simplified treatment of the neutrino physics with a parameterized neutrino luminosity, we explore the relationship between explosion time, mass accretion rate, and neutrino luminosity for a 15 Msun progenitor in 1D and 2D. We test three different equations of state commonly used in core-collapse simulations: the models of Lattimer & Swesty (1991) with incompressibility of 180 MeV and 220 MeV, and the model of Shen et al. (1998), in order of increasing stiffness. We find that for a given neutrino luminosity the time after bounce until explosion increases with the stiffness of the equation of state: the Lattimer & Swesty EOS explode more easily than that of Shen et al. We find this holds in both 1D and 2D, while for all models explosions are obtained more easily in 2D than in 1D. We also discuss the relevance of approximate instability criteria to realistic simulations.

We have searched for optical identifications for 79 Chandra X-ray sources that lie within the half-mass radius of the nearby, core-collapsed globular cluster NGC 6397, using deep Hubble Space Telescope Advanced Camera for Surveys Wide Field Channel imaging in H-alpha, R, and B. Photometry of these images allows us to classify candidate counterparts based on color-magnitude diagram location. In addition to recovering nine previously detected cataclysmic variables (CVs), we have identified six additional faint CV candidates, a total of 42 active binaries (ABs), two millisecond pulsars (MSPs), one candidate active galactic nucleus, and one candidate interacting galaxy pair. Of the 79 sources, 69 have a plausible optical counterpart. The 15 likely and possible CVs in NGC 6397 mostly fall into two groups: a brighter group of six for which the optical emission is dominated by contributions from the secondary and accretion disk, and a fainter group of seven for which the white dwarf dominates the optical emission. There are two possible transitional objects that lie between these groups. The faintest CVs likely lie near the minimum of the CV period distribution, where an accumulation is expected. The spatial distribution of the brighter CVs is much more centrally concentrated than those of the fainter CVs and the active binaries. This may represent the result of an evolutionary process in which CVs are produced by dynamical interactions, such as exchange reactions, near the cluster center and are scattered to larger orbital radii, over their lifetimes, as they age and become fainter.

Shock breakout is the earliest, readily-observable emission from a core-collapse supernova explosion. Observing supernova shock breakout may yield information about the nature of the supernova shock prior to exiting the progenitor and, in turn, about the core-collapse supernova mechanism itself. X-ray Outburst 080109, later associated with SN 2008D, is a very well-observed example of shock breakout from a core-collapse supernova. Despite excellent observational coverage and detailed modeling, fundamental information about the shock breakout, such as the radius of breakout and driver of the light curve time scale, is still uncertain. The models constructed for explaining the shock breakout emission from SN 2008D all assume spherical symmetry. We present a study of the observational characteristics of {\it aspherical} shock breakout from stripped-envelope core-collapse supernovae. We conduct two-dimensional, jet-driven supernova simulations from stripped-envelope progenitors and calculate the resulting shock breakout X-ray spectra and light curves. The X-ray spectra evolve significantly in time as the shocks expand outward and are not well-fit by single-temperature and radius black bodies. The time scale of the X-ray burst light curve of the shock breakout is related to the shock crossing time of the progenitor, not the much shorter light crossing time that sets the light curve time scale in spherical breakouts. This could explain the long shock breakout light curve time scale observed for XRO 080109/SN 2008D.

We present axisymmetric hydrodynamical simulations of the long-term accretion of a rotating GRB progenitor star, a "collapsar," onto the central compact object. The simulations were carried out with the adaptive mesh refinement code FLASH in two spatial dimensions and with an explicit shear viscosity. The evolution of the central accretion rate exhibits phases reminiscent of the long GRB gamma-ray and X-ray light curve, which lends support to the proposal that the luminosity is modulated by the central accretion rate. After a few tens of seconds, an accretion shock sweeps outward through the star. The formation and outward expansion of the accretion shock is accompanied with a sudden and rapid power-law decline in the central accretion rate Mdot ~ t^{-2.8}, which resembles the L_X ~ t^{-3} decline observed in the X-ray light curves. The collapsed, shock-heated stellar envelope settles into a thick, low-mass equatorial disk embedded within a massive, pressure-supported atmosphere. After a few hundred seconds, the inflow of low-angular-momentum material in the axial funnel reverses into an outflow from the surface of the thick disk. Meanwhile, the rapid decline of the accretion rate slows down, or even settles a in steady state with Mdot ~ 5x10^{-5} Msun/s, which resembles the "plateau" phase in the X-ray light curve. While the duration of the "prompt" phase depends on the resolution in our simulations, we provide an analytical model taking into account neutrino losses that estimates the duration to be ~20 s. The model suggests that the steep decline in GRB X-ray light curves is triggered by the circularization of the infalling stellar envelope at radii where the virial temperature is below ~10^{10} K, such that neutrino cooling shuts off and an outward expansion of the accretion shock becomes imminent.

The validity of the hypothesis that the massive black holes in high redshift quasars grew from stellar-sized "seeds" is contingent on a seed's ability to double its mass every few ten million years. This requires that the seed accrete at approximately the Eddington-limited rate. In the specific case of radiatively efficient quasiradial accretion in a metal-poor protogalactic medium, for which the Bondi accretion rate is often prescribed in cosmological simulations of massive black hole formation, we examine the effects of the radiation emitted near the black hole's event horizon on the structure of the surrounding gas flow. We find that the radiation pressure from photoionization significantly reduces the steady-state accretion rate and renders the quasiradial accretion flow unsteady and inefficient. The time-averaged accretion rates are a small fraction of the Eddington-limited accretion rate for Thomson scattering. The pressure of Ly-alpha photons trapped near the HII region surrounding the black hole may further attenuate the inflow. These results suggest that an alternative to quasiradial, radiatively efficient Bondi-like accretion should be sought to explain the rapid growth of quasar-progenitor seed black holes.

We explore the observational characteristics of jet-driven supernovae by simulating bipolar-jet-driven explosions in a red supergiant progenitor. We present results of four models in which we hold the injected kinetic energy at a constant $10^{51}$ ergs across all jet models but vary the specific characteristics of the jets to explore the influence of the nature of jets on the structure of the supernova ejecta. We evolve the explosions past shock-breakout and into quasi-homologous expansion of the supernova envelope into a red supergiant wind. The oppositely-directed, nickel-rich jets give a large-scale asymmetry that may account for the non-spherical excitation and substructure of spectral lines such as H$\alpha$ and He I 10830\AA. Jets with a large fraction of kinetic to thermal energy punch through the progenitor envelope and give rise to explosions that would be observed to be asymmetric from the earliest epochs, inconsistent with spectropolarimetric measurements of Type II supernovae. Jets with higher thermal energy fractions result in explosions that are roughly spherical at large radii but are significantly elongated at smaller radii, deep inside the ejecta, in agreement with the polarimetric observations. We present shock breakout light curves that indicate that strongly aspherical shock breakouts are incompatible with recent {\it GALEX} observations of shock breakout from red supergiant stars. Comparison with observations indicates that jets must deposit their kinetic energy efficiently throughout the ejecta while in the hydrogen envelope. Thermal energy-dominated jets satisfy this criterion and yield many of the observational characteristics of Type II supernovae.

We present the first results from two-dimensional simulations of radiatively-efficient accretion of metal-free gas onto intermediate-mass black holes. We fix the shape of the spectral energy distribution of the radiation produced near the event horizon and study the structure of the irradiated low-angular-momentum accretion flow over three orders of magnitude in radius from the black hole, 10^{14}-10^{17} cm for a 100 M_sun black hole. The luminosity of the central source is made to be proportional to the rate at which gas accretes across the inner boundary, which we set just inside the sonic radius. We find that photoionization heating and radiation pressure modify the structure of the flow. When the ambient gas density is 10^7 cm^{-3}, accretion is intermittent and on average reduced to 32% of the Eddington-limited rate, two orders of magnitude below the "Bondi" rate evaluated ignoring radiation, in agreement with simplified theoretical models. Even if the vicinity of the black hole is supplied with high density gas, accretion is rendered inefficient through heating and radiation pressure.

We investigate the interaction of nonthermal ions (protons and nuclei) accelerated in an ultrarelativistic blastwave with the pre-existing magnetic field of the medium into which the blastwave propagates. While particle acceleration processes such as diffusive shock acceleration can accelerate ions and electrons, the accelerated electrons suffer larger radiative losses. Under certain conditions, the ions can attain higher energies and reach farther ahead of the shock than the electrons, and so the nonthermal particles can be partially charge-separated. To compensate for the charge separation, the upstream plasma develops a return current, which, as it flows across the magnetic field, drives transverse acceleration of the upstream plasma and a growth of density contrast in the shock upstream. If the density contrast is strong by the time the fluid is shocked, vorticity is generated at the shock transition. The resulting turbulence can amplify the post-shock magnetic field to the levels inferred from gamma-ray burst afterglow spectra and light curves. Therefore, since the upstream inhomogeneities are induced by the ions accelerated in the shock, they are generic even if the blastwave propagates into a medium of uniform density. We speculate about the global structure of the shock precursor, and delineate several distinct physical regimes that are classified by an increasing distance from the shock and, correspondingly, a decreasing density of nonthermal particles that reach that distance.

Based on optical, IR and X-ray studies of Cas A, we propose a geometry for the remnant based on a "jet-induced" scenario with significant systematic departures from axial symmetry. In this model, the main jet axis is oriented in the direction of strong blue-shifted motion at an angle of 110 - 120 degrees East of North and about 40 - 50 degrees to the East of the line of sight. Normal to this axis would be an expanding torus as predicted by jet-induced models. In the proposed geometry, iron-peak elements in the main jet-like flow could appear "beyond" the portions of the remnant rich in silicon by projection effects, not the effect of mixing. In the context of the proposed geometry, the displacement of the compact object from the kinematic center of the remnant at a position angle of ~169 degrees can be accommodated if the motion of the compact object is near to, but slightly off from, the direction of the main "jet" axis by of order 30 degrees. In this model, the classical NE "jet," the SW "counter-jet" and other protrusions, particularly the "hole" in the North, are non-asymmetric flows approximately in the equatorial plane, e.g., out through the perimeter of the expanding torus, rather than being associated with the main jet. We explore the spoke-like flow in the equatorial plane in terms of Rayleigh-Taylor, Richtmyer-Meshkov and Kelvin-Helmholz instabilities and illustrate these instabilities with a jet-induced simulation.