Search for: All records

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 3D simulations of core-collapse supernovae, we find that in the early stages of the PNS evolution, typically 0.4 s 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 s 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.

We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity (CoRe) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to$\ell =m=4$. They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425.CoResimulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

We investigate r-process nucleosynthesis and kilonova emission resulting from binary neutron star (BNS) mergers based on a three-dimensional (3D) general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. The simulation includes a microphysical finite-temperature equation of state (EOS) and neutrino emission and absorption effects via a leakage scheme. We track the thermodynamic properties of the ejecta using Lagrangian tracer particles and determine its composition using the nuclear reaction network SkyNet. We investigate the impact of neutrinos on the nucleosynthetic yields by varying the neutrino luminosities during post-processing. The ejecta show a broad distribution with respect to their electron fraction Ye, peaking between ∼0.25–0.4 depending on the neutrino luminosity employed. We find that the resulting r-process abundance patterns differ from solar, with no significant production of material beyond the second r-process peak when using luminosities recorded by the tracer particles. We also map the HMNS outflows to the radiation hydrodynamics code SNEC and predict the evolution of the bolometric luminosity as well as broadband light curves of the kilonova. The bolometric light curve peaks on the timescale of a day and the brightest emission is seen in the infrared bands. This is the first direct calculation of the r-process yields and kilonova signal expected from HMNS winds based on 3D GRMHD simulations. For longer-lived remnants, these winds may be the dominant ejecta component producing the kilonova emission.

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 localization area too broad is still an open question. We simulate 28 BNS mergers with the chirp mass of GW190425 and mass ratio 1 ≤ q ≤ 1.67, using numerical-relativity simulations with finite-temperature, composition dependent equations of state (EOS) and neutrino radiation. The energy emitted in GWs is $\lesssim 0.083\mathrm{\, M_\odot }c^2$ with peak luminosity of 1.1–$2.4\times ~10^{58}/(1+q)^2\, {\rm {erg \, s^{-1}}}$. Dynamical ejecta and disc mass range between 5 × 10−6–10−3 and 10−5–$0.1 \mathrm{\, M_\odot }$, respectively. Asymmetric mergers, especially with stiff EOSs, unbind more matter and form heavier discs compared to equal mass binaries. The angular momentum of the disc is 8–$10\mathrm{\, M_\odot }~GM_{\rm {disc}}/c$ over three orders of magnitude in Mdisc. While the nucleosynthesis shows no peculiarity, the simulated kilonovae are relatively dim compared with GW170817. For distances compatible with GW190425, AB magnitudes are always dimmer than ∼20 mag for the B, r, and K bands, with brighter kilonovae associated to more asymmetric binaries and stiffer EOSs. 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.

We develop a method to compute synthetic kilonova light curves that combine numerical relativity simulations of neutron star mergers and the SNEC radiation–hydrodynamics code. We describe our implementation of initial and boundary conditions, r-process heating, and opacities for kilonova simulations. We validate our approach by carefully checking that energy conservation is satisfied and by comparing the SNEC results with those of two semi-analytic light-curve models. We apply our code to the calculation of colour light curves for three binaries having different mass ratios (equal and unequal mass) and different merger outcome (short-lived and long-lived remnants). We study the sensitivity of our results to hydrodynamic effects, nuclear physics uncertainties in the heating rates, and duration of the merger simulations. We find that hydrodynamics effects are typically negligible and that homologous expansion is a good approximation in most cases. However, pressure forces can amplify the impact of uncertainties in the radioactive heating rates. We also study the impact of shocks possibly launched into the outflows by a relativistic jet. None of our models match AT2017gfo, the kilonova in GW170817. This points to possible deficiencies in our merger simulations and kilonova models that neglect non-LTE effects and possible additional energy injection from the merger remnant and to the need to go beyond the assumption of spherical symmetry adopted in this work.

Neutrinos are copiously emitted by neutron star mergers, due to the high temperatures reached by dense matter during the merger and its aftermath. Neutrinos influence the merger dynamics and shape the properties of the ejecta, including the resultingr-process nucleosynthesis and kilonova emission. In this work, we analyse neutrino emission from a large sample of binary neutron star merger simulations in Numerical Relativity, covering a broad range of initial masses, nuclear equation of state and viscosity treatments. We extract neutrino luminosities and mean energies, and compute quantities of interest such as the peak values, peak broadnesses, time averages and decrease time scales. We provide a systematic description of such quantities, including their dependence on the initial parameters of the system. We find that for equal-mass systems the total neutrino luminosity (several$$10^{53}{\hbox {erg}~{\hbox {s}}^{-1}}$$${10}^{53}\text{erg}{\text{s}}^{-1}$) decreases as the reduced tidal deformability increases, as a consequence of the less violent merger dynamics. Similarly, tidal disruption in asymmetric mergers leads to systematically smaller luminosities. Peak luminosities can be twice as large as the average ones. Electron antineutrino luminosities dominate (initially by a factor of 2-3) over electron neutrino ones, while electron neutrinos and heavy flavour neutrinos have similar luminosities. Mean energies are nearly constant in time and independent on the binary parameters. Their values reflect the different decoupling temperature inside the merger remnant. Despite present uncertainties in neutrino modelling, our results provide a broad and physically grounded characterisation of neutrino emission, and they can serve as a reference point to develop more sophisticated neutrino transport schemes.

Strong magnetic fields play an important role in powering the emission of neutron stars. Nevertheless, a full understanding of the interior configuration of the field remains elusive. In this work, we present general relativistic magnetohydrodynamics (MHD) simulations of the magnetic field evolution in neutron stars lasting ${\sim } {880}\,$ms (∼6.5 Alfvén crossing periods) and up to resolutions of $0.1155\,$km using Athena++. We explore two different initial conditions, one with purely poloidal magnetic field and the other with a dominant toroidal component, and study the poloidal and toroidal field energies, the growth times of the various instability-driven oscillation modes, and turbulence. We find that the purely poloidal setup generates a toroidal field, which later decays exponentially reaching $1{{\ \rm per\ cent}}$ of the total magnetic energy, showing no evidence of reaching equilibrium. The initially stronger toroidal field setup, on the other hand, loses up to 20 per cent of toroidal energy and maintains this state till the end of our simulation. We also explore the hypothesis, drawn from previous MHD simulations, that turbulence plays an important role in the quasi-equilibrium state. An analysis of the spectra in our higher resolution setups reveals, however, that in most cases we are not observing turbulence at small scales, but rather a noisy velocity field inside the star. We also observe that the majority of the magnetic energy gets dissipated as heat increasing the internal energy of the star, while a small fraction gets radiated away as electromagnetic radiation.

We present a new moment-based energy-integrated neutrino transport code for neutron star merger simulations in general relativity. In the merger context, ours is the first code to include Doppler effects at all orders in υ/c, retaining all non-linear neutrino–matter coupling terms. The code is validated with a stringent series of tests. We show that the inclusion of full neutrino–matter coupling terms is necessary to correctly capture the trapping of neutrinos in relativistically moving media, such as in differentially rotating merger remnants. We perform preliminary simulations proving the robustness of the scheme in simulating ab-initio mergers to black hole collapse and long-term neutron star remnants up to ${\sim }70\,$ ms. The latter is the longest dynamical space-time, 3D, general relativistic simulations with full neutrino transport to date. We compare results obtained at different resolutions and using two different closures for the moment scheme. We do not find evidences of significant out-of-thermodynamic equilibrium effects, such as bulk viscosity, on the post-merger dynamics or gravitational wave emission. Neutrino luminosities and average energies are in good agreement with theory expectations and previous simulations by other groups using similar schemes. We compare dynamical and early wind ejecta properties obtained with M1 and with our older neutrino treatment. We find that the M1 results have systematically larger proton fractions. However, the differences in the nucleosynthesis yields are modest. This work sets the basis for future detailed studies spanning a wider set of neutrino reactions, binaries, and equations of state.