Based on our recent three-dimensional core-collapse supernova (CCSN) simulations including both exploding and non-exploding models, we study the detailed neutrino signals in representative terrestrial neutrino observatories, Super-Kamiokande (Hyper-Kamiokande), DUNE, JUNO, and IceCube. We find that the physical origin of difference in the neutrino signals between 1D and 3D is mainly proto-neutron-star (PNS) convection. We study the temporal and angular variations of the neutrino signals and discuss the detectability of the time variations driven by the spiral Standing Accretion Shock Instability (spiral SASI) when it emerges for non-exploding models. In addition, we determine that there can be a large angular asymmetry in the event rate ($\gtrsim 50 \%$), but that the time-integrated signal has a relatively modest asymmetry ($\lesssim 20 \%$). Both features are associated with the lepton-number emission self-sustained asymmetry (LESA) and the spiral SASI. Moreover, our analysis suggests that there is an interesting correlation between the total neutrino energy (TONE) and the cumulative number of neutrino events in each detector, a correlation that can facilitate data analyses of real observations. We demonstrate the retrieval of neutrino energy spectra for all flavors of neutrino by applying a novel spectrum reconstruction technique to the data from multiple detectors. We find that this new method is capable of estimating the TONE within the error of $\sim$20\% if the distance to the CCSN is $\lesssim 6$ kpc.

The gravitational-wave astronomical revolution began in 2015 with LIGO's observation of the coalescence of two stellar-mass black holes. Over the coming decades, ground-based detectors like LIGO will extend their reach, discovering thousands of stellar-mass binaries. In the 2030s, the space-based LISA will enable gravitational-wave observations of the massive black holes in galactic centres. Between LISA and ground-based observatories lies the unexplored decihertz gravitational-wave frequency band. Here, we propose a Decihertz Observatory to cover this band, and complement observations made by other gravitational-wave observatories. The decihertz band is uniquely suited to observation of intermediate-mass ($\sim 10^2$-$10^4 M_\odot$) black holes, which may form the missing link between stellar-mass and massive black holes, offering a unique opportunity to measure their properties. Decihertz observations will be able to detect stellar-mass binaries days to years before they merge and are observed by ground-based detectors, providing early warning of nearby binary neutron star mergers, and enabling measurements of the eccentricity of binary black holes, providing revealing insights into their formation. Observing decihertz gravitational-waves also opens the possibility of testing fundamental physics in a new laboratory, permitting unique tests of general relativity and the Standard Model of particle physics. Overall, a Decihertz Observatory will answer key questions about how black holes form and evolve across cosmic time, open new avenues for multimessenger astronomy, and advance our understanding of gravitation, particle physics and cosmology.

This paper presents the first systematic study of proto-neutron star (PNS) convection in three dimensions (3D) based on our latest numerical Fornax models of core-collapse supernova (CCSN). We confirm that PNS convection commonly occurs, and then quantify the basic physical characteristics of the convection. By virtue of the large number of long-term models, the diversity of PNS convective behavior emerges. We find that the vigor of PNS convection is not correlated with CCSN dynamics at large radii, but rather with the mass of PNS $-$ heavier masses are associated with stronger PNS convection. We find that PNS convection boosts the luminosities of $\nu_{\mu}$, $\nu_{\tau}$, $\bar{\nu}_{\mu}$, and $\bar{\nu}_{\tau}$ neutrinos, while the impact on other species is complex due to a competition of factors. Finally, we assess the consequent impact on CCSN dynamics and the potential for PNS convection to generate pulsar magnetic fields.

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.

Using our new state-of-the-art core-collapse supernova (CCSN) code Fornax, we explore the dependence upon spatial resolution of the outcome and character of three-dimensional (3D) supernova simulations. For the same 19-M$_{\odot}$ progenitor star, energy and radial binning, neutrino microphysics, and nuclear equation of state, changing only the number of angular bins in the $\theta$ and $\phi$ directions, we witness that our lowest resolution 3D simulation does not explode. However, when jumping progressively up in resolution by factors of two in each angular direction on our spherical-polar grid, models then explode, and explode slightly more vigorously with increasing resolution. This suggests that there can be a qualitative dependence of the outcome of 3D CCSN simulations upon spatial resolution. The critical aspect of higher spatial resolution is the adequate capturing of the physics of neutrino-driven turbulence, in particular its Reynolds stress. The greater numerical viscosity of lower-resolution simulations results in greater drag on the turbulent eddies that embody turbulent stress, and, hence, in a diminution of their vigor. Turbulent stress not only pushes the temporarily stalled shock further out, but bootstraps a concomitant increase in the deposited neutrino power. Both effects together lie at the core of the resolution dependence we observe.

We provide the time series and angular distributions of the neutrino and gravitational-wave emissions of eleven state-of-the-art three-dimensional non-rotating core-collapse supernova models and explore correlations between these signatures and the real-time dynamics of the shock and the proto-neutron-star core. The neutrino emissions are roughly isotropic on average, with instantaneous excursions about the mean inferred luminosity of as much as $\pm$20%. The deviation from isotropy is least for the "$\nu_{\mu}$"-type neutrinos and the lowest-mass progenitors. Instantaneous temporal luminosity variations along a given direction for exploding models average $\sim$2$-$4%, but can be as high as $\sim$10%. For non-exploding models, they can achieve $\sim$25%. The temporal variations in the neutrino emissions correlate with the temporal and angular variations in the mass accretion rate. We witness the LESA phenomenon in all our models and find that the vector direction of the LESA dipole and that of the inner Y$_\mathrm{e}$ distribution are highly correlated. For our entire set of 3D models, we find strong connections between the cumulative neutrino energy losses, the radius of the proto-neutron star, and the $f$-mode frequency of the gravitational wave emissions. When physically normalized, the progenitor-to-progenitor variation in any of these quantities is no more than $\sim$10%. Moreover, the reduced $f$-mode frequency is independent of time after bounce to better than $\sim$10%. Therefore, simultaneous measurement of gravitational waves and neutrinos from a given supernova event can be used synergistically to extract real physical quantities of the supernova core.

We optimize the third-generation gravitational-wave detector to maximize the range to detect core-collapse supernovae. Based on three-dimensional simulations for core-collapse and the corresponding gravitational-wave waveform emitted, the corresponding detection range for these waveforms is limited to within our galaxy even in the era of third-generation detectors. The corresponding event rate is two per century. We find from the waveforms that to detect core-collapse supernovae with an event rate of one per year, the gravitational-wave detectors need a strain sensitivity of 3$\times10^{-27}~$Hz$^{-1/2}$ in a frequency range from 100~Hz to 1500~Hz. We also explore detector configurations technologically beyond the scope of third-generation detectors. We find with these improvements, the event rate for gravitational-wave observations from CCSN is still low, but is improved to one in twenty years.

We study the gravitational wave signal from eight new 3D core-collapse supernova simulations. We show that the signal is dominated by $f$- and $g$-mode oscillations of the protoneutron star and its frequency evolution encodes the contraction rate of the latter, which, in turn, is known to depend on the star's mass, on the equation of state, and on transport properties in warm nuclear matter. A lower-frequency component of the signal, associated with the standing accretion shock instability, is found in only one of our models. Finally, we show that the energy radiated in gravitational waves is proportional to the amount of turbulent energy accreted by the protoneutron star.

Using the new state-of-the-art core-collapse supernova (CCSN) code F{\sc{ornax}}, we have simulated the three-dimensional dynamical evolution of the cores of 9-, 10-, 11-, 12-, and 13-M$_{\odot}$ stars from the onset of collapse. Stars from 8-M$_{\odot}$ to 13-M$_{\odot}$ constitute roughly 50% of all massive stars, so the explosive potential for this mass range is important to the overall theory of CCSNe. We find that the 9-, 10-, 11-, and 12-M$_{\odot}$ models explode in 3D easily, but that the 13-M$_{\odot}$ model does not. From these findings, and the fact that slightly more massive progenitors seem to explode \citep{vartanyan2019}, we suggest that there is a gap in explodability near 12-M$_{\odot}$ to 14-M$_{\odot}$ for non-rotating progenitor stars. Factors conducive to explosion are turbulence behind the stalled shock, energy transfer due to neutrino-matter absorption and neutrino-matter scattering, many-body corrections to the neutrino-nucleon scattering rate, and the presence of a sharp silicon-oxygen interface in the progenitor. Our 3D exploding models frequently have a dipolar structure, with the two asymmetrical exploding lobes separated by a pinched waist where matter temporarily continues to accrete. This process maintains the driving neutrino luminosty, while partially shunting matter out of the way of the expanding lobes, thereby modestly facilitating explosion. The morphology of all 3D explosions is characterized by multiple bubble structures with a range of low-order harmonic modes. Though much remains to be done in CCSN theory, these and other results in the literature suggest that, at least for these lower-mass progenitors, supernova theory is converging on a credible solution.

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.

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.

In this paper, we present the results of our three-dimensional, multi-group, multi-neutrino-species radiation/hydrodynamic simulation using the state-of-the-art code F{\sc{ornax}} of the terminal dynamics of the core of a non-rotating 16-M$_{\odot}$ stellar progenitor. The calculation incorporates redistribution by inelastic scattering, a correction for the effect of many-body interactions on the neutrino-nucleon scattering rates, approximate general relativity (including the effects of gravitational redshifts), velocity-dependent frequency advection, and an implementation of initial perturbations in the progenitor core. The model explodes within $\sim$100 milliseconds of bounce (near when the silicon-oxygen interface is accreted through the temporarily-stalled shock) and by the end of the simulation (here, $\sim$677 milliseconds after bounce) is accumulating explosion energy at a rate of $\sim$2.5$\times$10$^{50}$ ergs s$^{-1}$. The supernova explosion resembles an asymmetrical multi-plume structure, with one hemisphere predominating. The gravitational mass of the residual proto-neutron star at $\sim$677 milliseconds is $\sim$1.42 M$_{\odot}$. Even at the end of the simulation, explosion in most of the solid angle is accompanied by some accretion in an annular fraction at the wasp-like waist of the debris field. The ejecta electron fraction (Y$_e$) is distributed from $\sim$0.48 to $\sim$0.56, with most of the ejecta mass proton-rich. This may have implications for supernova nucleosynthesis, and could have a bearing on the p- and $\nu$p-processes and on the site of the first peak of the r-process. The ejecta spatial distributions of both Y$_e$ and mass density are predominantly in wide-angle plumes and large-scale structures, but are nevertheless quite patchy.

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 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 study gravitational waves (GWs) from a set of two-dimensional multi-group neutrino radiation hydrodynamic simulations of core-collapse supernovae (CCSNe). Our goal is to systematize the current knowledge about the post-bounce CCSN GW signal and recognize the templatable features that could be used by the ground-based laser interferometers. We demonstrate that starting from ~400ms after core bounce the dominant GW signal represents the fundamental quadrupole (l=2) oscillation mode (f-mode) of the proto-neutron star (PNS), which can be accurately reproduced by a linear perturbation analysis of the angle-averaged PNS profile. Before that, in the time interval between ~200 and ~400ms after bounce, the dominant mode has two radial nodes and represents a g-mode. We associate the high-frequency noise in the GW spectrograms above the main signal with p-modes, while below the dominant frequency there is a region with very little power. The collection of models presented here summarizes the dependence of the CCSN GW signal on the progenitor mass, equation of state, many-body corrections to the neutrino opacity, and rotation. Weak dependence of the dominant GW frequency on the progenitor mass motivates us to provide a simple fit for it as a function of time, which can be used as a prior when looking for CCSN candidates in the LIGO data.

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.

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.

Recent discoveries of circumbinary planets by Kepler mission provide motivation for understanding their birthplaces - protoplanetary disks around stellar binaries with separations <1 AU. We explore properties and evolution of such circumbinary disks focusing on modification of their structure caused by tidal coupling to the binary. We develop a set of analytical scaling relations describing viscous evolution of the disk properties, which are verified and calibrated using 1D numerical calculations with realistic inputs. Injection of angular momentum by the central binary suppresses mass accretion onto the binary and causes radial distribution of the viscous angular momentum flux F_J to be different from that in a standard accretion disk around a single star with no torque at the center. Disks with no mass accretion at the center develop F_J profile which is flat in radius. Radial profiles of temperature and surface density are also quite different from those in disks around single stars. Damping of the density waves driven by the binary and viscous dissipation dominate heating of the inner disk (within 1-2 AU), pushing the iceline beyond 3-5 AU, depending on disk mass and age. Irradiation by the binary governs disk thermodynamics beyond ~10 AU. However, self-shadowing by the hot inner disk may render central illumination irrelevant out to ~20 AU. Spectral energy distribution of a circumbinary disk exhibits a distinctive bump around 10 micron, which may facilitate identification of such disks around unresolved binaries. Efficient tidal coupling to the disk drives orbital inspiral of the binary and may cause low-mass and compact binaries to merge into a single star within the disk lifetime. We generally find that circumbinary disks present favorable sites for planet formation (despite wider zone of volatile depletion), in agreement with the statistics of Kepler circumbinary planets.