Search for: All records

We perform the first 3D ab-initio general-relativistic neutrino-radiation hydrodynamics of a long-lived neutron star merger remnant spanning a fraction of its cooling timescale. We find that neutrino cooling becomes the dominant energy loss mechanism after the gravitational-wave dominated phase (∼20 ms postmerger). Electron flavor antineutrino luminosity dominates over electron flavor neutrino luminosity at early times, resulting in a secular increase of the electron fraction in the outer layers of the remnant. However, the two luminosities become comparable ∼20–40 ms postmerger. A dense gas of electron antineutrinos is formed in the outer core of the remnant at densities ∼10^14.5g cm⁻³, corresponding to temperature hot spots. The neutrinos account for ∼10% of the lepton number in this region. Despite the negative radial temperature gradient, the radial entropy gradient remains positive, and the remnant is stably stratified according to the Ledoux criterion for convection. A massive accretion disk is formed from the material squeezed out of the collisional interface between the stars. The disk carries a large fraction of the angular momentum of the system, allowing the remnant massive neutron star to settle to a quasi-steady equilibrium within the region of possible, stable, rigidly rotating configurations. The remnant is differentially rotating, but it is stable against the magnetorotational instability. Other MHD mechanisms operating on longer timescales are likely responsible for the removal of the differential rotation. Our results indicate the remnant massive neutron star is thus qualitatively different from a protoneutron stars formed in core-collapse supernovae.

 

It is widely believed that the binary neutron star merger GW190425 produced a black hole promptly upon merger. Motivated by the potential association with the fast radio burst FRB 20190425A, which took place 2.5 h after the merger, we revisit the question of the outcome of GW190425 by means of numerical relativity simulations. We show that current laboratory and astrophysical constraints on the equation of state of dense matter do not rule out the formation of a long-lived remnant. However, the formation of a stable remnant would have produced a bright kilonova, in tension with upper limits by ZTF at the location and time of FRB 20190425A. Moreover, the ejecta would have been optically thick to radio emission for days to months, preventing a putative FRB from propagating out. The predicted dispersion measure is also several orders of magnitude larger than that observed for FRB 20190425A. Our results indicate that FRB 20190425A and GW190425 are not associated. However, we cannot completely rule out the formation of a long-lived remnant, due to the incomplete coverage of the relevant sky regions. More observations of GW190425-like events, including potential upper limit, have the potential to constrain nuclear physics. To this aim, it is important that follow-up observational campaigns of gravitational wave events are informed by the properties of the source, such as their chirp mass, and we urge the LIGO-Virgo-KAGRA collaboration to promptly release them publicly.

 

We present a 3D general-relativistic magnetohydrodynamic simulation of a short-lived neutron star remnant formed in the aftermath of a binary neutron star merger. The simulation uses an M1 neutrino transport scheme to track neutrino–matter interactions and is well suited to studying the resulting nucleosynthesis and kilonova emission. A magnetized wind is driven from the remnant and ejects neutron-rich material at a quasi-steady-state rate of 0.8 × 10⁻¹M_⊙s⁻¹. We find that the ejecta in our simulations underproducer-process abundances beyond the secondr-process peak. For sufficiently long-lived remnants, these outflowsalonecan produce blue kilonovae, including the blue kilonova component observed for AT2017gfo.

 

We present fully general relativistic 3D numerical simulations of accretion-induced collapse (AIC) of white dwarfs (WDs). We evolve three different WD models (non-rotating, rotating at 80 per cent and 99 per cent of the Keplerian mass shedding limit) that collapse due to electron capture. For each of these models, we provide a detailed analysis of their gravitational waves (GWs), neutrinos, and electromagnetic counterpart and discuss their detectability. Our results suggest that fast rotating AICs could be detectable up to a distance of 8 Mpc with third-generation GW observatories, and up to 1 Mpc with LIGO. AIC progenitors are expected to have large angular momentum due to their accretion history, which is a determining factor for their stronger GW emission compared to core-collapse supernovae (CCSNe). Regarding neutrino emission, we found no significant difference between AICs and CCSNe. In the electromagnetic spectrum, we find that AICs are two orders of magnitude fainter than type Ia supernovae. Our work places AICs as realistic targets for future multimessenger searches with third generation ground-based GW detectors.

 

We study the impact of finite-temperature effects in numerical-relativity simulations of binary neutron star mergers with microphysical equations of state and neutrino transport in which we vary the effective nucleon masses in a controlled way. We find that, as the specific heat is increased, the merger remnants become colder and more compact due to the reduced thermal pressure support. Using a full Bayesian analysis, we demonstrate that this effect will be measurable in the postmerger gravitational wave signal with next-generation observatories at signal-to-noise ratios of 15.

 

ABSTRACT We present a systematic numerical relativity study of the impact of different physics input and grid resolution in binary neutron star mergers. We compare simulations employing a neutrino leakage scheme, leakage plus M0 scheme, the M1 transport scheme, and pure hydrodynamics. Additionally, we examine the effect of a sub-grid scheme for turbulent viscosity. We find that the overall dynamics and thermodynamics of the remnant core are robust, implying that the maximum remnant density could be inferred from gravitational wave observations. Black hole collapse instead depends significantly on viscosity and grid resolution. Differently from recent work, we identify possible signatures of neutrino effects in the gravitational waves only at the highest resolutions considered; new high-resolution simulations will be thus required to build accurate gravitational wave templates to observe these effects. Different neutrino transport schemes impact significantly mass, geometry, and composition of the remnant’s disc and ejecta; M1 simulations show systematically larger proton fractions, reaching maximum values larger than 0.4. r-process nucleosynthesis yields reflect the different ejecta compositions; they are in agreement and reproduce residual solar abundances only if M0 or M1 neutrino transport schemes are adopted. We compute kilonova light curves using spherically-symmetric radiation-hydrodynamics evolutions up to 15 d post-merger, finding that they are mostly sensitive to the ejecta mass and electron fraction; accounting for multiple ejecta components appears necessary for reliable light curve predictions. We conclude that advanced neutrino schemes and resolutions higher than current standards are essential for robust long-term evolutions and detailed astrophysical predictions.