We have conducted nineteen state-of-the-art 3D core-collapse supernova simulations spanning a broad range of progenitor masses. This is the largest collection of sophisticated 3D supernova simulations ever performed. We have found that while the majority of these models explode, not all do, and that even models in the middle of the available progenitor mass range may be less explodable. This does not mean that those models for which we did not witness explosion would not explode in Nature, but that they are less prone to explosion than others. One consequence is that the "compactness" measure is not a metric for explodability. We find that lower-mass massive star progenitors likely experience lower-energy explosions, while the higher-mass massive stars likely experience higher-energy explosions. Moreover, most 3D explosions have a dominant dipole morphology, have a pinched, wasp-waist structure, and experience simultaneous accretion and explosion. We reproduce the general range of residual neutron-star masses inferred for the galactic neutron-star population. The most massive progenitor models, however, in particular vis \`a vis explosion energy, need to be continued for longer physical times to asymptote to their final states. We find that while the majority of the inner ejecta have Y$_e = 0.5$, there is a substantial proton-rich tail. This result has important implications for the nucleosynthetic yields as a function of progenitor. Finally, we find that the non-exploding models eventually evolve into compact inner configurations that experience a quasi-periodic spiral SASI mode. We otherwise see little evidence of the SASI in the exploding models.

This paper describes the design and implementation of our new multi-group, multi-dimensional radiation hydrodynamics (RHD) code Fornax and provides a suite of code tests to validate its application in a wide range of physical regimes. Instead of focusing exclusively on tests of neutrino radiation hydrodynamics relevant to the core-collapse supernova problem for which Fornax is primarily intended, we present here classical and rigorous demonstrations of code performance relevant to a broad range of multi-dimensional hydrodynamic and multi-group radiation hydrodynamic problems. Our code solves the comoving-frame radiation moment equations using the M1 closure, utilizes conservative high-order reconstruction, employs semi-explicit matter and radiation transport via a high-order time stepping scheme, and is suitable for application to a wide range of astrophysical problems. To this end, we first describe the philosophy, algorithms, and methodologies of Fornax and then perform numerous stringent code tests, that collectively and vigorously exercise the code, demonstrate the excellent numerical fidelity with which it captures the many physical effects of radiation hydrodynamics, and show excellent strong scaling well above 100k MPI tasks.

For a suite of fourteen core-collapse models during the dynamical first second after bounce, we calculate the detailed neutrino "light" curves expected in the underground neutrino observatories Super-Kamiokande, DUNE, JUNO, and IceCube. These results are given as a function of neutrino-oscillation modality (normal or inverted hierarchy) and progenitor mass (specifically, post-bounce accretion history), and illuminate the differences between the light curves for 1D (spherical) models that don't explode with the corresponding 2D (axisymmetric) models that do. We are able to identify clear signatures of explosion (or non-explosion), the post-bounce accretion phase, and the accretion of the silicon/oxygen interface. In addition, we are able to estimate the supernova detection ranges for various physical diagnostics and the distances out to which various temporal features embedded in the light curves might be discerned. We find that the progenitor mass density profile and supernova dynamics during the dynamical explosion stage should be identifiable for a supernova throughout most of the galaxy in all the facilities studied and that detection by any one of them, but in particular more than one in concert, will speak volumes about the internal dynamics of supernovae.

We present results of 2D axisymmetric core-collapse supernova simulations, employing the FORNAX code, of nine progenitor models spanning 12 to 25 M$_{\odot}$ and evolved over a 20,000-km grid. We find that four of the nine models explode with inelastic scattering off electrons and neutrons as well as the many-body correction to neutrino-nucleon scattering opacities. We show that these four models feature sharp Si-O interfaces in their density profiles, and that the corresponding dip in density reduces the accretion rate around the stalled shock and prompts explosion. The non-exploding models lack such a steep feature, suggesting that Si-O interface is one key to explosion. Furthermore, we show that all of the non-exploding models can be nudged to explosion with modest changes to macrophysical inputs, including moderate rotation and perturbations to infall velocities, as well as to microphysical inputs, including changes to neutrino-nucleon interaction rates, suggesting that all the models are perhaps close to criticality. Exploding models have energies of few $\times$10$^{50}$ ergs at the end of our simulation, and are rising, suggesting the need to continue these simulations over larger grids and for longer times to reproduce the energies seen in Nature. We find that the morphology of the explosion contributes to the explosion energy, with more isotropic ejecta producing larger explosion energies. However, we do not find evidence for the Lepton-number Emission Self-Sustained Asymmetry. Finally, we look at PNS properties and explore the role of dimension in our simulations. We find that convection in the proto-neutron star (PNS) produces larger PNS radii as well as greater "$\nu_\mu$" luminosities in 2D compared to 1D.

We explore with self-consistent 2D F{\sc{ornax}} simulations the dependence of the outcome of collapse on many-body corrections to neutrino-nucleon cross sections, the nucleon-nucleon bremsstrahlung rate, electron capture on heavy nuclei, pre-collapse seed perturbations, and inelastic neutrino-electron and neutrino-nucleon scattering. Importantly, proximity to criticality amplifies the role of even small changes in the neutrino-matter couplings, and such changes can together add to produce outsized effects. When close to the critical condition the cumulative result of a few small effects (including seeds) that individually have only modest consequence can convert an anemic into a robust explosion, or even a dud into a blast. Such sensitivity is not seen in one dimension and may explain the apparent heterogeneity in the outcomes of detailed simulations performed internationally. A natural conclusion is that the different groups collectively are closer to a realistic understanding of the mechanism of core-collapse supernovae than might have seemed apparent.

Momentum deposition by radiation pressure from young, massive stars may help to destroy molecular clouds and unbind stellar clusters by driving large-scale outflows. We extend our previous numerical radiation hydrodynamic study of turbulent, star-forming clouds to analyze the detailed interaction between non-ionizing UV radiation and the cloud material. Our simulations trace the evolution of gas and star particles through self-gravitating collapse, star formation, and cloud destruction via radiation-driven outflows. These models are idealized in that we include only radiation feedback and adopt an isothermal equation of state. Turbulence creates a structure of dense filaments and large holes through which radiation escapes, such that only ~50% of the radiation is (cumulatively) absorbed by the end of star formation. The surface density distribution of gas by mass as seen by the central cluster is roughly lognormal with sigma_ln(Sigma) = 1.3-1.7, similar to the externally-projected surface density distribution. This allows low surface density regions to be driven outwards to nearly 10 times their initial escape speed v_esc. Although the velocity distribution of outflows is broadened by the lognormal surface density distribution, the overall efficiency of momentum injection to the gas cloud is reduced because much of the radiation escapes. The mean outflow velocity is approximately twice the escape speed from the initial cloud radius. Our results are also informative for understanding galactic-scale wind driving by radiation, in particular the relationship between velocity and surface density for individual outflow structures, and the resulting velocity and mass distributions arising from turbulent sources.

We present an implementation of an adaptive ray tracing (ART) module in the Athena hydrodynamics code that accurately and efficiently handles the radiative transfer involving multiple point sources on a three-dimensional Cartesian grid. We adopt a recently proposed parallel algorithm that uses non-blocking, asynchronous MPI communications to accelerate transport of rays across the computational domain. We validate our implementation through several standard test problems including the propagation of radiation in vacuum and the expansions of various types of HII regions. Additionally, scaling tests show that the cost of a full ray trace per source remains comparable to that of the hydrodynamics update on up to $\sim 10^3$ processors. To demonstrate application of our ART implementation, we perform a simulation of star cluster formation in a marginally bound, turbulent cloud, finding that its star formation efficiency is $12\%$ when both radiation pressure forces and photoionization by UV radiation are treated. We directly compare the radiation forces computed from the ART scheme with that from the M1 closure relation. Although the ART and M1 schemes yield similar results on large scales, the latter is unable to resolve the radiation field accurately near individual point sources.

We present new 1D (spherical) and 2D (axisymmetric) simulations of electron-capture (EC) and low-mass iron-core-collapse supernovae (SN). We consider six progenitor models: the ECSN progenitor from Nomoto (1984, 1987); two ECSN-like low-mass low-metallicity iron core progenitors from Heger (private communication); and the 9-, 10-, and 11-$M_\odot$ (zero-age main sequence) progenitors from Sukhbold et al. (2016). We confirm that the ECSN and ESCN-like progenitors explode easily even in 1D with explosion energies of up to a 0.15 Bethes ($1 {\rm B} \equiv 10^{51}\ {\rm erg}$), and are a viable mechanism for the production of very low-mass neutron stars. However, the 9-, 10-, and 11-$M_\odot$ progenitors do not explode in 1D and are not even necessarily easier to explode than higher-mass progenitor stars in 2D. We study the effect of perturbations and of changes to the microphysics and we find that relatively small changes can result in qualitatively different outcomes, even in 1D, for models sufficiently close to the explosion threshold. Finally, we revisit the impact of convection below the protoneutron star (PNS) surface. We analyze, 1D and 2D evolutions of PNSs subject to the same boundary conditions. We find that the impact of PNS convection has been underestimated in previous studies and could result in an increase of the neutrino luminosity by up to factors of two.

We perform the first self-consistent, time-dependent, multi-group calculations in two dimensions (2D) to address the consequences of using the ray-by-ray+ transport simplification in core-collapse supernova simulations. Such a dimensional reduction is employed by many researchers to facilitate their resource-intensive calculations. Our new code (F{\sc{ornax}}) implements multi-D transport, and can, by zeroing out transverse flux terms, emulate the ray-by-ray+ scheme. Using the same microphysics, initial models, resolution, and code, we compare the results of simulating 12-, 15-, 20-, and 25-M$_{\odot}$ progenitor models using these two transport methods. Our findings call into question the wisdom of the pervasive use of the ray-by-ray+ approach. Employing it leads to maximum post-bounce/pre-explosion shock radii that are almost universally larger by tens of kilometers than those derived using the more accurate scheme, typically leaving the post-bounce matter less bound and artificially more "explodable." In fact, for our 25-M$_{\odot}$ progenitor, the ray-by-ray+ model explodes, while the corresponding multi-D transport model does not. Therefore, in two dimensions the combination of ray-by-ray+ with the axial sloshing hydrodynamics that is a feature of 2D supernova dynamics can result in quantitatively, and perhaps qualitatively, incorrect results.

Radiation feedback from stellar clusters is expected to play a key role in setting the rate and efficiency of star formation in giant molecular clouds (GMCs). To investigate how radiation forces influence realistic turbulent systems, we have conducted a series of numerical simulations employing the {\it Hyperion} radiation hydrodynamics solver, considering the regime that is optically thick to ultraviolet (UV) and optically thin to infrared (IR) radiation. Our model clouds cover initial surface densities between $\Sigma_{\rm cl,0} \sim 10-300~M_{\odot}~{\rm pc^{-2}}$, with varying initial turbulence. We follow them through turbulent, self-gravitating collapse, formation of star clusters, and cloud dispersal by stellar radiation. All our models display a lognormal distribution of gas surface density $\Sigma$; for an initial virial parameter $\alpha_{\rm vir,0} = 2$, the lognormal standard deviation is $\sigma_{\rm ln \Sigma} = 1-1.5$ and the star formation rate coefficient $\varepsilon_{\rm ff,\bar\rho} = 0.3-0.5$, both of which are sensitive to turbulence but not radiation feedback. The net star formation efficiency $\varepsilon_\mathrm{final}$ increases with $\Sigma_{\rm cl,0}$ and decreases with $\alpha_{\rm vir,0}$. We interpret these results via a simple conceptual framework, whereby steady star formation increases the radiation force, such that local gas patches at successively higher $\Sigma$ become unbound. Based on this formalism (with fixed $\sigma_{\rm ln \Sigma}$), we provide an analytic upper bound on $\varepsilon_\mathrm{final}$, which is in good agreement with our numerical results. The final star formation efficiency depends on the distribution of Eddington ratios in the cloud and is strongly increased by turbulent compression of gas.

Radiation feedback from young star clusters embedded in giant molecular clouds (GMCs) is believed to be important to the control of star formation. For the most massive and dense clouds, including those in which super star clusters (SSCs) are born, pressure from reprocessed radiation exerted on dust grains may disperse a significant portion of the cloud mass back into the interstellar medium (ISM). Using our radiaton hydrodynamics (RHD) code, Hyperion, we conduct a series of numerical simulations to test this idea. Our models follow the evolution of self-gravitating, strongly turbulent clouds in which collapsing regions are replaced by radiating sink particles representing stellar clusters. We evaluate the dependence of the star formation efficiency (SFE) on the size and mass of the cloud and $\kappa$, the opacity of the gas to infrared (IR) radiation. We find that the single most important parameter determining the evolutionary outcome is $\kappa$, with $\kappa \gtrsim 15 \text{ cm}^2 \text{ g}^{-1}$ needed to disrupt clouds. For $\kappa = 20-40 \text{ cm}^2 \text{ g}^{-1}$, the resulting SFE=50-70% is similar to empirical estimates for some SSC-forming clouds. The opacities required for GMC disruption likely apply only in dust-enriched environments. We find that the subgrid model approach of boosting the direct radiation force $L/c$ by a "trapping factor" equal to a cloud's mean IR optical depth can overestimate the true radiation force by factors of $\sim 4-5$. We conclude that feedback from reprocessed IR radiation alone is unlikely to significantly reduce star formation within GMCs unless their dust abundances or cluster light-to-mass ratios are enhanced.

HD 166191 has been identified by several studies as hosting a rare and extremely bright warm debris disc with an additional outer cool disc component. However, an alternative interpretation is that the star hosts a disc that is currently in transition between a full gas disc and a largely gas-free debris disc. With the help of new optical to mid-IR spectra and Herschel imaging, we argue that the latter interpretation is supported in several ways: i) we show that HD 166191 is co-moving with the ~4 Myr-old Herbig Ae star HD 163296, suggesting that the two have the same age, ii) the disc spectrum of HD 166191 is well matched by a standard radiative transfer model of a gaseous protoplanetary disc with an inner hole, and iii) the HD 166191 mid-IR silicate feature is more consistent with similarly primordial objects. We note some potential issues with the debris disc interpretation that should be considered for such extreme objects, whose lifetime at the current brightness is mush shorter than the stellar age, or in the case of the outer component requires a mass comparable to the solid component of the Solar nebula. These aspects individually and collectively argue that HD 166191 is a 4-5 Myr old star that hosts a gaseous transition disc. Though it does not argue in favour of either scenario, we find strong evidence for 3-5 um disc variability. We place HD 166191 in context with discs at different evolutionary stages, showing that it is a potentially important object for understanding the protoplanetary to debris disc transition.

We describe a module for the Athena code that solves the gray equations of radiation hydrodynamics (RHD), based on the first two moments of the radiative transfer equation. We use a combination of explicit Godunov methods to advance the gas and radiation variables including the non-stiff source terms, and a local implicit method to integrate the stiff source terms. We adopt the M1 closure relation and include all leading source terms. We employ the reduced speed of light approximation (RSLA) with subcycling of the radiation variables in order to reduce computational costs. Our code is dimensionally unsplit in one, two, and three space dimensions and is parallelized using MPI. The streaming and diffusion limits are well-described by the M1 closure model, and our implementation shows excellent behavior for a problem with a concentrated radiation source containing both regimes simultaneously. Our operator-split method is ideally suited for problems with a slowly varying radiation field and dynamical gas flows, in which the effect of the RSLA is minimal. We present an analysis of the dispersion relation of RHD linear waves highlighting the conditions of applicability for the RSLA. To demonstrate the accuracy of our method, we utilize a suite of radiation and RHD tests covering a broad range of regimes, including RHD waves, shocks, and equilibria, which show second-order convergence in most cases. As an application, we investigate radiation-driven ejection of a dusty, optically thick shell in the interstellar medium (ISM). Finally, we compare the timing of our method with other well-known iterative schemes for the RHD equations. Our code implementation, Hyperion, is suitable for a wide variety of astrophysical applications and will be made freely available on the Web.

The optical afterglow of long-duration GRB 071003 is among the brightest yet to be detected from any GRB, with R ~ 12 mag in KAIT observations starting 42 s after the GRB trigger, including filtered detections during prompt emission. However, our high S/N ratio afterglow spectrum displays only extremely weak absorption lines at what we argue is the host redshift of z = 1.60435 - in contrast to the three other, much stronger Mg II absorption systems observed at lower redshifts. Together with Keck adaptive optics observations which fail to reveal a host galaxy coincident with the burst position, our observations suggest a halo progenitor and offer a cautionary tale about the use of Mg II for GRB redshift determination. We present early through late-time observations spanning the electromagnetic spectrum, constrain the connection between the prompt emission and early variations in the light curve (we observe no correlation), and discuss possible origins for an unusual, marked rebrightening that occurs a few hours after the burst: likely either a late-time refreshed shock or a wide-angle secondary jet. Analysis of the late-time afterglow is most consistent with a wind environment, suggesting a massive star progenitor. Together with GRB 070125, this may indicate that a small but significant portion of star formation in the early universe occurred far outside what we consider a normal galactic disk.

The robotic ROTSE-III telescope network detected prompt optical emission contemporaneous with the gamma-ray emission of Swift events GRB051109A and GRB051111. Both datasets have continuous coverage at high signal-to-noise levels from the prompt phase onwards, thus the early observations are readily compared to the Swift XRT and BAT high energy detections. In both cases, the optical afterglow is established, declining steadily during the prompt emission. For GRB051111, there is evidence of an excess optical component during the prompt emission. The component is consistent with the flux spectrally extrapolated from the gamma-rays, using the gamma-ray spectral index. A compilation of spectral information from previous prompt detections shows that such a component is unusual. The existence of two prompt optical components - one connected to the high-energy emission, the other to separate afterglow flux, as indicated in GRB051111 - is not compatible with a simple ``external-external'' shock model for the GRB and its afterglow.