Search for: All records

Abstract Observations of neutron star mergers have the potential to unveil detailed physics of matter and gravity in regimes inaccessible by other experiments. Quantitative comparisons to theory and parameter estimation require nonlinear numerical simulations. However, the detailed physics of energy and momentum transfer between different scales, and the formation and interaction of small scale structures, which can be probed by detectors, are not captured by current simulations. This is where turbulence enters neutron star modelling. This review will outline the theory and current status of turbulence modelling for relativistic neutron star merger simulations. 

Abstract Quasiperiodic oscillations (QPOs) have been recently discovered in the short gamma-ray bursts (GRBs) 910711 and 931101B. Their frequencies are consistent with those of the quasiradial and quadrupolar oscillations of binary neutron star (BNS) merger remnants, as obtained in numerical relativity simulations. These simulations reveal quasi-universal relations between the remnant oscillation frequencies and the tidal coupling constant of the binaries. Under the assumption that the observed QPOs are due to these postmerger oscillations, we use the frequency–tide relations in a Bayesian framework to infer the source redshift, as well as the chirp mass and the binary tidal deformability of the BNS progenitors for GRBs 910711 and 931101B. We further use this inference to estimate bounds on the mass–radius relation for neutron stars. By combining the estimates from the two GRBs, we find a 68% credible range $R_{1.4}=12.48_{-0.40}^{+0.41}$km for the radius of a neutron star with massM= 1.4M_⊙, which is one of the tightest bounds to date. 

Abstract We present the first seconds-long 2D general relativistic neutrino magnetohydrodynamic simulations of accretion-induced collapse (AIC) in rapidly rotating, strongly magnetized white dwarfs (WDs), which might originate as remnants of double-WD mergers. This study examines extreme combinations of magnetic fields and rotation rates, motivated both by the need to address the limitations of 2D axisymmetric simulations and to explore the physics of AIC under rare conditions that, while yet to be observationally confirmed, may be consistent with current theoretical models and account for unusual events. Under these assumptions, our results demonstrate that, if realizable, such systems can generate relativistic jets and neutron-rich outflows with properties consistent with long gamma-ray bursts (LGRBs) accompanied by kilonovae, such as GRB 211211A and GRB 230307A. These findings highlight the potential role of AIC in heavyr-process element production and offer a framework for understanding rare LGRBs associated with kilonova emission. Longer-duration 3D simulations are needed to fully capture magnetic field amplification, resolve instabilities, and determine the fate of the energy retained by the magnetar at the end of the simulations. 

Abstract We study mass ejection from a binary neutron star merger producing a long-lived massive neutron star remnant with general-relativistic neutrino-radiation hydrodynamics simulations. In addition to outflows generated by shocks and tidal torques during and shortly after the merger, we observe the appearance of a wind driven by spiral density waves in the disk. This spiral-wave-driven outflow is predominantly located close to the disk orbital plane and have a broad distribution of electron fractions. At higher latitudes, a high electron-fraction wind is driven by neutrino radiation. The combined nucleosynthesis yields from all the ejecta components is in good agreement with Solar abundance measurements. 

Abstract We study the ringdown signal of black holes formed in prompt-collapse binary neutron star mergers. We analyze data from 47 numerical relativity simulations. We show that the $(\ell =2,m=2)$and $(\ell =2,m=1)$multipoles of the gravitational wave signal are well fitted by decaying damped exponentials, as predicted by black-hole perturbation theory. We show that the ratio of the amplitude in the two modes depends on the progenitor binary mass ratioqand reduced tidal parameter $\tilde{\mathrm{\Lambda }}$. Unfortunately, the numerical uncertainty in our data is too large to fully quantify this dependency. If confirmed, these results will enable novel tests of general relativity in the presence of matter with next-generation gravitational-wave observatories. 

ABSTRACT 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. 

Understanding the details ofr-process nucleosynthesis in binary neutron star merger (BNSM) ejecta is key to interpreting kilonova observations and identifying the role of BNSMs in the origin of heavy elements. We present a self-consistent, two-dimensional, ray-by-ray radiation-hydrodynamic evolution of BNSM ejecta with an online nuclear network (NN) up to a timescale of days. For the first time, an initial numerical relativity ejecta profile composed of the dynamical component and spiral-wave and disk winds is evolved including detailedr-process reactions and nuclear heating effects. A simple model for the jet energy deposition is also included. Our simulation highlights that the common approach of relating in postprocessing the final nucleosynthesis yields to the initial thermodynamic profile of the ejecta can lead to inaccurate predictions. Moreover, we find that neglecting the details of the radiation-hydrodynamic evolution of the ejecta in nuclear calculations can introduce deviations of up to 1 order of magnitude in the final abundances of several elements, including very light and secondr-process peak elements. The presence of a jet affects element production only in the innermost part of the polar ejecta, and it does not alter the global nucleosynthesis results. Overall, our analysis shows that employing an online NN improves the reliability of nucleosynthesis and kilonova light-curve predictions. 

We present a Bayesian framework for joint and coherent analyses of multimessenger binary neutron star signals. The method, implemented in our bajes infrastructure, incorporates a joint likelihood for multiple datasets, support for various semi-analytical kilonova models and numerical-relativity (NR) informed relations for the mass ejecta, as well as a technique to include and marginalize over modeling uncertainties. As a first application, we analyze the gravitational-wave GW170817 and the kilonova AT2017gfo data. These results are then combined with the most recent X-ray pulsars analyses of PSR J0030+0451 and PSR J0740+6620 to obtain EOS constraints.Various constraints on the mass-radius diagram and neutron star properties are then obtained by resampling over a set of ten million parametrized EOS built under minimal assumptions. We find that a joint and coherent approach improves the inference of the extrinsic parameters (distance) and, among the instrinc parameters, the mass ratio. The inclusion of NR informed relations strongly improves over the case of using an agnostic prior on the intrinsic parameters. Comparing Bayes factors, we find that the two observations are better explained by the common source hypothesis only by assuming NR-informed relations. These relations break some of the degeneracies in the employed kN models. The EOS inference folding-in PSR J0952-0607 minimum-maximum mass, PSR J0030+0451 and PSR J0740+6620 data constrains, among other quantities, the neutron star radius to R1.4=12.30−0.56+0.81R1.4​=12.30−0.56+0.81​ km (R1.4=13.20−0.90+0.91R1.4​=13.20−0.90+0.91​ km) and the maximum mass to Mmax=2.28−0.17+0.25 M⊙Mmax​=2.28−0.17+0.25​ M⊙​ (Mmax=2.32−0.19+0.30 M⊙Mmax​=2.32−0.19+0.30​ M⊙​) where the ST+PDT (PDT-U) analysis of Vinciguerra et a (2023) for PSR J0030+0451 is employed. Hence, the systematics on PSR J0030+0451 data reduction currently dominate the mass-radius diagram constraints.