Search

The future detection of gravitational waves (GWs) from a galactic core-collapse supernova will provide information on the physics inside protoneutron stars (PNS). In this work, we apply three different classification methods for the PNS non-radial oscillation modes: Cowling classification, Generalized Cowling Nomenclature (GCN), and a Classification Based on Modal Properties (CBMP). Using PNS models from $3$D simulations of core-collapse supernovae, we find that in the early stages of the PNS evolution, typically before $0.4$ seconds after the bounce, the Cowling classification is inconsistent, but the GCN and the CBMP provide complementary information that helps to understand the evolution of the modes. In the GCN, we note several avoided crossings as the mode frequencies evolve at early times, while the CBMP tracks the modes across the avoided crossings. We verify that the strongest emission of GWs by the PNS corresponds to the $f$-mode in the GCN, indicating that the mode trapping region alternates between the core and the envelope at each avoided crossing. At later times, approximately $0.4$ seconds after the bounce, the three classification methods present a similar description of the mode spectrum. We use our results to test universal relations for the PNS modes according to their classification and find that the behaviour of the universal relations for $f$- and $p$-modes is remarkably simple in the CBMP., Comment: 11 pages, 8 figures. Matches the version accepted on MNRAS

Nuclear Astrophysics is a field at the intersection of nuclear physics and astrophysics, which seeks to understand the nuclear engines of astronomical objects and the origin of the chemical elements. This white paper summarizes progress and status of the field, the new open questions that have emerged, and the tremendous scientific opportunities that have opened up with major advances in capabilities across an ever growing number of disciplines and subfields that need to be integrated. We take a holistic view of the field discussing the unique challenges and opportunities in nuclear astrophysics in regards to science, diversity, education, and the interdisciplinarity and breadth of the field. Clearly nuclear astrophysics is a dynamic field with a bright future that is entering a new era of discovery opportunities., Comment: 96 pages. Submitted to Journal of Physics G

GW190425 was the second gravitational wave (GW) signal compatible with a binary neutron star (BNS) merger detected by the Advanced LIGO and Advanced Virgo detectors. Since no electromagnetic counterpart was identified, whether the associated kilonova was too dim or the localisation area too broad is still an open question. We simulate 28 BNS mergers with the chirp mass of GW190425 and mass ratio $1 \leq q \leq 1.67$, using numerical-relativity simulations with finite temperature, composition dependent nuclear equation of state (EOS) and neutrino radiation. The energy emitted in GWs is $\lesssim 0.083 M_{\odot} c^2$ with peak luminosity of $1.1-2.4 \times 10^{58} {\rm erg~s^{-1}}/(1+q)^2$. Dynamical ejecta and disc mass range between $5 \times 10^{-6} - \sim 10^{-3}~M_{\odot}$ and $10^{-5} - 0.1~M_{\odot}$, respectively. Asymmetric mergers, especially with stiff EOS, unbind more matter and to form heavier discs compared to equal mass binaries. The angular momentum of the disc is $8-10 M_{\odot}~GM_{\rm disc}/c$ over three orders of magnitude in $M_{\rm disc}$. While the nucleosynthesis shows no peculiarity, the simulated kilonovae are relatively dim compared with the GW170817 event. For distances compatible with GW190425, AB magnitudes are always dimmer than $\sim20~{\rm mag}$ for the $B$, $r$ and $K$ bands, with brighter kilonovae associated to more asymmetric binaries and stiffer EOS. We suggest that, even assuming a good coverage of GW190425's sky location, the kilonova could hardly have been detected by present wide-field surveys and no firm constraints on the binary parameters or EOS can be argued from the lack of the detection., Comment: 22 pages, 17 figures

Using 250 neutron star merger simulations with microphysics, we explore for the first time the role of nuclear incompressibility in the prompt collapse threshold for binaries with different mass ratios. We demonstrate that observations of prompt collapse thresholds, either from binaries with two different mass ratios or with one mass ratio but combined with the knowledge of the maximum neutron star mass or compactness, will constrain the incompressibility at the maximum neutron star density, $K_{\rm max}$, to within tens of percent. This, otherwise inaccessible, measure of $K_{\rm max}$ can potentially reveal the presence of hyperons or quarks inside neutron stars., Comment: 7 pages, 4 figures; 2 tables and 4 figures in the appendix

For the first $\sim3$ years after the binary neutron star merger event GW 170817 the radio and X-ray radiation has been dominated by emission from a structured relativistic off-axis jet propagating into a low-density medium with n $< 0.01\,\rm{cm^{-3}}$. We report on observational evidence for an excess of X-ray emission at $\delta t>900$ days after the merger. With $L_x\approx5\times 10^{38}\,\rm{erg\,s^{-1}}$ at 1234 days, the recently detected X-ray emission represents a $\ge 3.2\,\sigma$ (Gaussian equivalent) deviation from the universal post jet-break model that best fits the multi-wavelength afterglow at earlier times. In the context of JetFit afterglow models, current data represent a departure with statistical significance $\ge 3.1\,\sigma$, depending on the fireball collimation, with the most realistic models showing excesses at the level of $\ge 3.7\,\sigma$. A lack of detectable 3 GHz radio emission suggests a harder broad-band spectrum than the jet afterglow. These properties are consistent with the emergence of a new emission component such as synchrotron radiation from a mildly relativistic shock generated by the expanding merger ejecta, i.e. a kilonova afterglow. In this context, we present a set of ab-initio numerical-relativity BNS merger simulations that show that an X-ray excess supports the presence of a high-velocity tail in the merger ejecta, and argues against the prompt collapse of the merger remnant into a black hole. Radiation from accretion processes on the compact-object remnant represents a viable alternative. Neither a kilonova afterglow nor accretion-powered emission have been observed before, as detections of BNS mergers at this phase of evolution are unprecedented., Comment: 30 pages, 13 figures, Accepted for publication in ApJL

We study the production of very light elements ($Z < 20$) in the dynamical and spiral-wave wind ejecta of binary neutron star mergers by combining detailed nucleosynthesis calculations with the outcome of numerical relativity merger simulations. All our models are targeted to GW170817 and include neutrino radiation. We explore different finite-temperature, composition dependent nuclear equations of state and binary mass ratios, and find that hydrogen and helium are the most abundant light elements. For both elements, the decay of free neutrons is the driving nuclear reaction. In particular, $\sim 0.5-2 \times 10^{-6} M_{\odot}$ of hydrogen are produced in the fast expanding tail of the dynamical ejecta, while $\sim 1.5-11 \times 10^{-6} M_{\odot}$ of Helium are synthesized in the bulk of the dynamical ejecta, usually in association with heavy r-process elements. By computing synthetic spectra, we find that the possibility of detecting hydrogen and helium features in kilonova spectra is very unlikely for fiducial masses and luminosities, even when including non local thermodynamics equilibrium effects. The latter could be crucial to observe He lines a few days after merger for faint kilonovae or for luminous kilonovae ejecting large masses of helium. Finally, we compute the amount of strontium synthesized in the dynamical and spiral-wave wind ejecta, and find that it is consistent with (or even larger than, in the case of a long lived remnant) the one required to explain early spectral features in the kilonova of GW170817., Comment: 21 pages, 9 figures, Accepted for publication in The Astrophysical Journal

We present Chandra and VLA observations of GW170817 at ~521-743 days post merger, and a homogeneous analysis of the entire Chandra data set. We find that the late-time non-thermal emission follows the expected evolution from an off-axis relativistic jet, with a steep temporal decay $F_{\nu}\propto t^{-1.95\pm0.15}$ and a simple power-law spectrum $F_{\nu}\propto \nu^{-0.575\pm0.007}$. We present a new method to constrain the merger environment density based on diffuse X-ray emission from hot plasma in the host galaxy and we find $n\le 9.6 \times 10^{-3}\,\rm{cm^{-3}}$. This measurement is independent from inferences based on the jet afterglow modeling and allows us to partially solve for model degeneracies. The updated best-fitting model parameters with this density constraint are a fireball kinetic energy $E_0 = 1.5_{-1.1}^{+3.6}\times 10^{49}\,\rm{erg}$ ($E_{iso}= 2.1_{-1.5}^{+6.4}\times10^{52}\, \rm{erg}$), jet opening angle $\theta_{0}= 5.9^{+1.0}_{-0.7}\,\rm{deg}$ with characteristic Lorentz factor $\Gamma_j = 163_{-43}^{+23}$, expanding in a low-density medium with $n_0 = 2.5_{-1.9}^{+4.1} \times 10^{-3}\, \rm{cm^{-3}}$ and viewed $\theta_{obs} = 30.4^{+4.0}_{-3.4}\, \rm{deg}$ off-axis. The synchrotron emission originates from a power-law distribution of electrons with $p=2.15^{+0.01}_{-0.02}$. The shock microphysics parameters are constrained to $\epsilon_{\rm{e}} = 0.18_{-0.13}^{+0.30}$ and $\epsilon_{\rm{B}}=2.3_{-2.2}^{+16.0} \times 10^{-3}$. We investigate the presence of X-ray flares and find no statistically significant evidence of $\ge2.5\sigma$ of temporal variability at any time. Finally, we use our observations to constrain the properties of synchrotron emission from the deceleration of the fastest kilonova ejecta with energy $E_k^{KN}\propto (\Gamma\beta)^{-\alpha}$ into the environment, finding that shallow stratification indexes $\alpha\le6$ are disfavored., Comment: version accepted for publication in ApJL, 13 pages, 6 figures

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., Comment: submitted to Astro2020

We conduct a series of numerical experiments into the nature of three-dimensional (3D) hydrodynamics in the postbounce stalled-shock phase of core-collapse supernovae using 3D general-relativistic hydrodynamic simulations of a $27$-$M_\odot$ progenitor star with a neutrino leakage/heating scheme. We vary the strength of neutrino heating and find three cases of 3D dynamics: (1) neutrino-driven convection, (2) initially neutrino-driven convection and subsequent development of the standing accretion shock instability (SASI), (3) SASI dominated evolution. This confirms previous 3D results of Hanke et al. 2013, ApJ 770, 66 and Couch & Connor 2014, ApJ 785, 123. We carry out simulations with resolutions differing by up to a factor of $\sim$4 and demonstrate that low resolution is artificially favorable for explosion in the 3D convection-dominated case, since it decreases the efficiency of energy transport to small scales. Low resolution results in higher radial convective fluxes of energy and enthalpy, more fully buoyant mass, and stronger neutrino heating. In the SASI-dominated case, lower resolution damps SASI oscillations. In the convection-dominated case, a quasi-stationary angular kinetic energy spectrum $E(\ell)$ develops in the heating layer. Like other 3D studies, we find $E(\ell) \propto \ell^{-1}$ in the "inertial range," while theory and local simulations argue for $E(\ell) \propto \ell^{-5/3}$. We argue that current 3D simulations do not resolve the inertial range of turbulence and are affected by numerical viscosity up to the energy containing scale, creating a "bottleneck" that prevents an efficient turbulent cascade., Comment: 24 pages, 15 figures. Accepted for publication in The Astrophysical Journal. Added one figure and made minor modifications to text according to suggestions from the referee