Researcher profile

David Radice

David Radice contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
20works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2026arXiv

Impact of Neutrino Flavor Conversions on Neutron Star Merger Dynamics, Ejecta, Nucleosynthesis, and Multi-Messenger Signals

We present numerical relativity simulations of binary neutron star mergers incorporating neutrino flavor transformations triggered by fast flavor instability, quantum many-body effects, or potential beyond standard model physics. In both long-lived and short-lived remnant scenarios, neutrino flavor conversions modify species-dependent neutrino luminosities and mean energies, and drive the matter towards more neutron rich conditions. They produce up to $300\%$ more neutron rich ejecta and significantly boost the r-process yields, especially in low-density, near-equatorial outflows. We identify regions unstable to fast flavor instabilities and find that these instabilities persist despite flavor conversions. We further test the sensitivity to the equilibration timescale of the flavor conversions, finding that slower flavor conversions can interact with thermodynamic equilibration, and increase the neutron richness of the ejecta. Flavor conversions may also contribute to stronger gravitational wave and neutrino emissions, pointing to a correlation between neutrino transport and merger dynamics. These results highlight the potential impact of flavor conversions while motivating future work to improve on theoretical understanding of flavor instabilities in global simulations.

preprint2023arXiv

Detectability of QCD phase transitions in binary neutron star mergers: Bayesian inference with the next generation gravitational wave detectors

We study the detectability of postmerger QCD phase transitions in neutron star binaries with next-generation gravitational-wave detectors Cosmic Explorer and Einstein Telescope. We perform numerical relativity simulations of neutron star mergers with equations of state that include a quark deconfinement phase transition through either a Gibbs or Maxwell construction. These are followed by Bayesian parameter estimation of the associated gravitational-wave signals using the $\tt{NRPMw}$ waveform model, with priors inferred from the analysis of the inspiral signal. We assess the ability of the model to measure the postmerger peak frequency $f_2^{\rm peak}$ and identify aspects that should be improved in the model. We show that, even at postmerger signal to noise ratios as low as 10, the model can distinguish (at the 90% level) $f_2^{\rm peak}$ between binaries with and without a phase transition in most cases. Phase-transition induced deviations in the $f_2^{\rm peak}$ from the predictions of equation-of-state insensitive relations can also be detected if they exceed $1.6\,σ$. Our results suggest that next-generation gravitational wave detectors can measure phase transition effects in binary neutron star mergers. However, unless the phase transition is ``strong'', disentangling it from other hadronic physics uncertainties will require significant theory improvements.

preprint2023arXiv

Revealing Phase Transition in Dense Matter with Gravitational Wave Spectroscopy of Binary Neutron Star Mergers

We use numerical relativity simulations of binary neutron star mergers to show that high density deconfinement phase transitions (PTs) to quark matter can be probed using multimodal postmerger gravitational wave (GW) spectroscopy. Hadron-quark PTs suppress the one-armed spiral instability in the remnant. This is manifested in an anti-correlation between the energy carried in the $l=2, m=1$ GW mode and energy density gap which separates the two phases. Consequently, a single measurement of the signal-to-noise ratios of the $l=2, m=1$ and $l=2, m=2$ GW modes could constrain the energy density gap of the PT.

preprint2022arXiv

Binary neutron star merger simulations with neutrino transport and turbulent viscosity: impact of different schemes and grid resolution

We present a systematic numerical relativity study of the impact of different treatment of microphysics and grid resolution in binary neutron star mergers. We consider series of simulations at multiple resolutions comparing hydrodynamics, neutrino leakage scheme, leakage augmented with the M0 scheme and the more consistent M1 transport scheme. Additionally, we consider the impact of a sub-grid scheme for turbulent viscosity. We find that viscosity helps to stabilise the remnant against gravitational collapse but grid resolution has a larger impact than microphysics on the remnant's stability. The gravitational wave (GW) energy correlates with the maximum remnant density, that can be thus inferred from GW observations. M1 simulations shows the emergence of a neutrino trapped gas that locally decreases the temperature a few percent when compared to the other simulation series. This out-of-thermodynamics equilibrium effect does not alter the GW emission at the typical resolutions considered for mergers. Different microphysics treatments impact significantly mass, geometry and composition of the remnant's disc and ejecta. M1 simulations show systematically larger proton fractions. The different ejecta compositions reflect into the nucleosynthesis yields, that are robust only if both neutrino emission and absorption are simulated. Synthetic kilonova light curves calculated by means of spherically-symmetric radiation-hydrodynamics evolutions up to 15 days post-merger are mostly sensitive to ejecta's mass and composition; they can be reliably predicted only including the various ejecta components. We conclude that advanced microphysics in combination with resolutions higher than current standards appear essential for robust long-term evolutions and astrophysical predictions.

preprint2022arXiv

Constraints on the maximum densities of neutron stars from postmerger gravitational waves with third-generation observations

Using data from 289 numerical relativity simulations of merging binary neutron stars, we identify, for the first time, a robust quasi-universal relation connecting the postmerger peak gravitational-wave frequency and the value of the density at the center of the maximum mass nonrotating neutron star. This relation offers a new possibility for precision equation-of-state constraints with next-generation ground-based gravitational-wave interferometers. Mock Einstein Telescope observations of fiducial events indicate that Bayesian inferences can constrain the maximum density to ${\sim}15\%$ ($90\%$ confidence level) for a single signal at the minimum sensitivity threshold for a detection. If the postmerger signal is included in a full-spectrum (inspiral-merger-postmerger) analysis of such signal, the pressure-density function can be tightly constrained up to the maximum density, and the maximum neutron star mass can be measured with an accuracy better than $12\%$ ($90\%$ confidence level).

preprint2022arXiv

Exotic Compact Objects: The Dark White Dwarf

Several dark matter models allow for the intriguing possibility of exotic compact object formation. These objects might have unique characteristics that set them apart from their baryonic counterparts. Furthermore, gravitational wave observations of their mergers may provide the only direct window on a potentially entirely hidden sector. Here we discuss dark white dwarfs, starting with an overview of the microphysical model and analytic scaling relations of macroscopic properties derived from the non-relativistic limit. We use the full relativistic formalism to confirm these scaling relations and demonstrate that dark white dwarfs, if they exist, would have masses and tidal deformabilities that are very different from those of baryonic compact objects. Further, and most importantly, we demonstrate that dark white dwarf mergers would be detectable by current or planned gravitational observatories across several orders of magnitude in the particle-mass parameter space. Lastly, we find universal relations analogous to the compactness-Love and binary Love relations in neutron star literature. Using these results, we show that gravitational wave observations would constrain the properties of the dark matter particles constituting these objects.

preprint2022arXiv

Gravitational Waves from Core-Collapse Supernovae

We summarize our current understanding of gravitational wave emission from core-collapse supernovae. We review the established results from multi-dimensional simulations and, wherever possible, provide back-of-the-envelope calculations to highlight the underlying physical principles. The gravitational waves are predominantly emitted by protoneutron star oscillations. In slowly rotating cases, which represent the most common type of the supernovae, the oscillations are excited by multi-dimensional hydrodynamic instabilities, while in rare rapidly rotating cases, the protoneutron star is born with an oblate deformation due to the centrifugal force. The gravitational wave signal may be marginally visible with current detectors for a source within our galaxy, while future third-generation instruments will enable more robust and detailed observations. The rapidly rotating models that develop non-axisymmetric instabilities may be visible up to a megaparsec distance with the third-generation detectors. Finally, we discuss strategies for multi-messenger observations of supernovae.

preprint2022arXiv

Long-term GRMHD simulation of magnetic field in isolated neutron stars

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 simulations of the magnetic field evolution in neutron stars lasting 500 ms (5 Alfven crossing times) and up to resolutions of 0.231 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% 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% 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 reveal, 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.

preprint2022arXiv

Neutrino emission from binary neutron star mergers: characterizing light curves and mean energies

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 resulting $r$-process nucleosynthesis and kilonova emission. In this work, we analyze neutrino emission from a large sample of merger radiation hydrodynamics 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}{\rm erg~s^{-1}}$) decreases for increasing reduced tidal deformability, 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 characterization of neutrino emission, and they can serve as a reference point to develop more sophisticated neutrino transport schemes.

preprint2022arXiv

Numerical relativity simulations of prompt collapse mergers: threshold mass and phenomenological constraints on neutron star properties after GW170817

We determine the threshold mass for prompt (no bounce) black hole formation in equal-mass neutron star (NS) mergers using a new set of 227 numerical relativity simulations. We consider 23 phenomenological and microphysical finite temperature equations of state (EOS), including models with hyperons and first-order phase transitions to deconfined quarks. We confirm the existence of EOS-insensitive relations between the threshold mass, the binary tidal parameter at the threshold ($Λ_{th}$), the maximum mass of nonrotating NSs, and the radii of reference mass NSs. We correct the systematic errors in previously reported fitting coefficients that were obtained with approximate general-relativity simulations. We combine the EOS-insensitive relations, phenomenological constraints on NS properties and observational data from GW170817 to derive an improved lower limit on radii of maximum mass and 1.6 M$_\odot$ NS of 9.81 km and 10.90 km, respectively. We also constrain the radius and quadrupolar tidal deformability ($Λ$) of a 1.4 $M_\odot$ NS to be larger than 10.74 km and 172, respectively. We consider uncertainties in all independent parameters -- fitting coefficients as well as GW170817 masses while reporting the range of radii constraints. We introduce new methods to constrain the upper as well as lower limit of NS maximum mass using future BNS detections and their identification as prompt or delayed collapse. With future observations it will be possible to derive even tighter constraints on the properties of matter at and above nuclear density using the method proposed in this work.

preprint2022arXiv

Radiation hydrodynamics modeling of kilonovae with SNEC

We develop a method to compute synthetic kilonova light curves that combines numerical relativity simulations of neutron star mergers and the $\texttt{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 $\texttt{SNEC}$ results with those of two semi-analytic light curve models. We apply our code to the calculation of color 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 which 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.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Numerical relativity for next-generation gravitational-wave probes of fundamental physics

The next generation of gravitational-wave detectors, conceived to begin operations in the 2030s, will probe fundamental physics with exquisite sensitivity. These observations will measure the equation of state of dense nuclear matter in the most extreme environments in the universe, reveal with exquisite fidelity the nonlinear dynamics of warped spacetime, put general relativity to the strictest test, and perhaps use black holes as cosmic particle detectors. Achieving each of these goals will require a new generation of numerical relativity simulations that will run at scale on the supercomputers of the 2030s to achieve the necessary accuracy, which far exceeds the capabilities of numerical relativity and high-performance computing infrastructures available today.

preprint2020arXiv

A systematic study of proto-neutron star convection in three-dimensional core-collapse supernova simulations

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 $ν_μ$, $ν_τ$, $\barν_μ$, and $\barν_τ$ 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.

preprint2020arXiv

Accretion-induced prompt black hole formation in asymmetric neutron star mergers, dynamical ejecta and kilonova signals

We present new numerical relativity results of neutron star mergers with chirp mass $1.188M_\odot$ and mass ratios $q=1.67$ and $q=1.8$ using finite-temperature equations of state (EOS), approximate neutrino transport and a subgrid model for magnetohydrodynamics-induced turbulent viscosity. The EOS are compatible with nuclear and astrophysical constraints and include a new microphysical model derived from ab-initio calculations based on the Brueckner-Hartree-Fock approach. We report for the first time evidence for accretion-induced prompt collapse in high-mass-ratio mergers, in which the tidal disruption of the companion and its accretion onto the primary star determine prompt black hole formation. As a result of the tidal disruption, an accretion disc of neutron-rich and cold matter forms with baryon masses ${\sim}0.15M_\odot$, and it is significantly heavier than the remnant discs in equal-masses prompt collapse mergers. Massive dynamical ejecta of order ${\sim}0.01M_\odot$ also originate from the tidal disruption. They are neutron rich and expand from the orbital plane with a crescent-like geometry. Consequently, bright, red and temporally extended kilonova emission is predicted from these mergers. Our results show that prompt black hole mergers can power bright electromagnetic counterparts for high-mass-ratio binaries, and that the binary mass ratio can be in principle constrained from multimessenger observations.

preprint2020arXiv

Binary Neutron Star Merger Simulations with a Calibrated Turbulence Model

Magnetohydrodynamic (MHD) turbulence in neutron star (NS) merger remnants can impact their evolution and multimessenger signatures, complicating the interpretation of present and future observations. Due to the high Reynolds numbers and the large computational costs of numerical relativity simulations, resolving all the relevant scales of the turbulence will be impossible for the foreseeable future. Here, we adopt a method to include subgrid-scale turbulence in moderate resolution simulations by extending the large-eddy simulation (LES) method to general relativity (GR). We calibrate our subgrid turbulence model with results from very-high-resolution GRMHD simulations, and we use it to perform NS merger simulations and study the impact of turbulence. We find that turbulence has a quantitative, but not qualitative impact on the evolution of NS merger remnants, on their gravitational wave signatures, and on the outflows generated in binary NS mergers. Our approach provides a viable path to quantify uncertainties due to turbulence in NS mergers.

preprint2020arXiv

kiloHertz gravitational waves from binary neutron star remnants: time-domain model and constraints on extreme matter

The remnant star of a neutron star merger is an anticipated loud source of kiloHertz gravitational waves that conveys unique information on the equation of state of hot matter at extreme densities. Observations of such signals are hampered by the photon shot noise of ground-based interferometers and pose a challenge for gravitational-wave astronomy. We develop an analytical time-domain waveform model for postmerger signals informed by numerical relativity simulations. The model completes effective-one-body waveforms for quasi-circular nonspinning binaries in the kiloHertz regime. We show that a template-based analysis can detect postmerger signals with a minimal signal-to-noise ratios (SNR) of 8, corresponding to GW170817-like events for third-generation interferometers. Using Bayesian model selection and the complete inspiral-merger-postmerger waveform model it is possible to infer whether the merger outcome is a prompt collapse to a black hole or a remnant star. In the latter case, the radius of the maximum mass (most compact) nonrotating neutron star can be determined to kilometer precision. We demonstrate the feasibility of inferring the stiffness of the equation of state at extreme densities using the quasiuniversal relations deduced from numerical-relativity simulations.

preprint2020arXiv

Numerical Relativity Simulations of the Neutron Star Merger GW170817: Long-Term Remnant Evolutions, Winds, Remnant Disks, and Nucleosynthesis

We present a systematic numerical-relativity study of the dynamical ejecta, winds and nucleosynthesis in neutron star merger remnants. Binaries with the chirp mass compatible with GW170817, different mass ratios, and five microphysical equations of state (EOS) are simulated with an approximate neutrino transport and a subgrid model for magnetohydrodynamics turbulence up to 100 milliseconds postmerger. Spiral density waves propagating from the neutron star remnant to the disk trigger a wind with mass flux ${\sim}0.1{-}0.5\,{\rm M_\odot/s}$ persisting for the entire simulation as long as the remnant does not collapse to black hole. This wind has average electron fraction $\gtrsim 0.3$ and average velocity ${\sim}0.1-0.17\,$c and thus is a site for the production of weak $r$-process elements (mass number $A<195$). Disks around long-lived remnants have masses ${\sim}0.1{-}0.2\,{\rm M_\odot}$, temperatures peaking at $\lesssim10\,$MeV near the inner edge, and a characteristic double-peak distribution in entropy resulting from shocks propagating through the disk. The dynamical and spiral-wave ejecta computed in our targeted simulations are not compatible with those inferred from AT2017gfo using two-components kilonova models. Rather, they indicate that multi-component kilonova models including disk winds are necessary to interpret AT2017gfo. The nucleosynthesis in the combined dynamical ejecta and spiral-wave wind in the comparable-mass long-lived mergers robustly accounts for all the $r$-process peaks, from mass number ${\sim}75$ to actinides in terms of solar abundances. Total abundandes are weakly dependent on the EOS, while the mass ratio affect the production of first peak elements.

preprint2020arXiv

The Dynamics of Binary Neutron Star Mergers and of GW170817

With the first observation of a binary neutron star merger through gravitational waves and light GW170817, compact binary mergers have now taken the center stage in nuclear astrophysics. They are thought to be one of the main astrophysical sites of production of r-process elements, and merger observations have become a fundamental tool to constrain the properties of matter. Here, we review our current understanding of the dynamics of neutron star mergers, in general, and of GW170817 in particular. We discuss the physical processes governing the inspiral, merger, and postmerger evolution, and we highlight the connections between these processes, the dynamics, and the multimessenger observables. Finally, we discuss open questions and issues in the field and the need to address them through a combination of better theoretical models and new observations.

preprint2019arXiv

Inferring prompt black-hole formation in neutron star mergers from gravitational-wave data

The gravitational-wave GW170817 is associated to the inspiral phase of a binary neutron star coalescence event. The LIGO-Virgo detectors sensitivity at high frequencies was not sufficient to detect the signal corresponding to the merger and post-merger phases. Hence, the question whether the merger outcome was a prompt black hole formation or not must be answered using either the pre-merger gravitational wave signal or electromagnetic counterparts. In this work we present two methods to infer the probability of prompt black hole formation, using the analysis of the inspiral gravitational-wave signal. Both methods combine the posterior distribution from the gravitational-wave data analysis with numerical relativity results. One method relies on the use of phenomenological models for the equation of state and on the estimate of the collapse threshold mass. The other is based on the estimate of the tidal polarizability parameter $\tildeΛ$ that is correlated in an equation-of-state agnostic way with the prompt BH formation. We analyze GW170817 data and find that the two methods consistently predict a probability of ~ 50-70% for prompt black-hole formation, which however may significantly decrease below 10% if the maximum mass constraint from PSR J0348+0432 or PSR J0740+6620 is imposed.

preprint2019arXiv

The Overarching Framework of Core-Collapse Supernova Explosions as Revealed by 3D Fornax Simulations

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 &#34;compactness&#34; 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 à 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.