BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David Radice

David Radice appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.HEgr-qcastro-ph.SRphysics.flu-dynastro-ph.COmath-phmath.MPphysics.comp-phastro-ph.IMhep-ph

Catalog footprint

What is connected

36works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

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

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 "compactness" 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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dynamical Mass Ejection from Binary Neutron Star Mergers

We present fully general-relativistic simulations of binary neutron star mergers with a temperature and composition dependent nuclear equation of state. We study the dynamical mass ejection from both quasi-circular and dynamical-capture eccentric mergers. We systematically vary the level of our treatment of the microphysics to isolate the effects of neutrino cooling and heating and we compute the nucleosynthetic yields of the ejecta. We find that eccentric binaries can eject significantly more material than quasi-circular binaries and generate bright infrared and radio emission. In all our simulations the outflow is composed of a combination of tidally- and shock-driven ejecta, mostly distributed over a broad $\sim 60^\circ$ angle from the orbital plane, and, to a lesser extent, by thermally driven winds at high latitudes. Ejecta from eccentric mergers are typically more neutron rich than those of quasi-circular mergers. We find neutrino cooling and heating to affect, quantitatively and qualitatively, composition, morphology, and total mass of the outflows. This is also reflected in the infrared and radio signatures of the binary. The final nucleosynthetic yields of the ejecta are robust and insensitive to input physics or merger type in the regions of the second and third r-process peaks. The yields for elements on the first peak vary between our simulations, but none of our models is able to explain the Solar abundances of first-peak elements without invoking additional first-peak contributions from either neutrino and viscously-driven winds operating on longer timescales after the mergers, or from core-collapse supernovae.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

How loud are neutron star mergers?

We present results from the first large parameter study of neutron star mergers using fully general relativistic simulations with finite-temperature microphysical equations of state and neutrino cooling. We consider equal and unequal-mass binaries drawn from the galactic population and simulate each binary with three different equations of state. Our focus is on the emission of energy and angular momentum in gravitational waves in the postmerger phase. We find that the emitted gravitational-wave energy in the first $\sim$$10\,\mathrm{ms}$ of the life of the resulting hypermassive neutron star (HMNS) is about twice the energy emitted over the entire inspiral history of the binary. The total radiated energy per binary mass is comparable to or larger than that of nonspinning black hole inspiral-mergers. About $0.8-2.5\%$ of the binary mass-energy is emitted at kHz frequencies in the early HMNS evolution. We find a clear dependence of the postmerger GW emission on binary configuration and equation of state and show that it can be encoded as a broad function of the binary tidal coupling constant $κ^T_2$. Our results also demonstrate that the dimensionless spin of black holes resulting from subsequent HMNS collapse are limited to $\lesssim0.7-0.8$. This may significantly impact the neutrino pair annihilation mechanism for powering short gamma-ray bursts (sGRB).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Neutrino-Driven Convection in Core-Collapse Supernovae: High-Resolution Simulations

We present results from high-resolution semi-global simulations of neutrino-driven convection in core-collapse supernovae. We employ an idealized setup with parametrized neutrino heating/cooling and nuclear dissociation at the shock front. We study the internal dynamics of neutrino-driven convection and its role in re-distributing energy and momentum through the gain region. We find that even if buoyant plumes are able to locally transfer heat up to the shock, convection is not able to create a net positive energy flux and overcome the downwards transport of energy from the accretion flow. Turbulent convection does, however, provide a significant effective pressure support to the accretion flow as it favors the accumulation of energy, mass and momentum in the gain region. We derive an approximate equation that is able to explain and predict the shock evolution in terms of integrals of quantities such as the turbulent pressure in the gain region or the effects of non-radial motion of the fluid. We use this relation as a way to quantify the role of turbulence in the dynamics of the accretion shock. Finally, we investigate the effects of grid resolution, which we change by a factor 20 between the lowest and highest resolution. Our results show that the shallow slopes of the turbulent kinetic energy spectra reported in previous studies are a numerical artefact. Kolmogorov scaling is progressively recovered as the resolution is increased.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Shock-Turbulence Interaction in Core-Collapse Supernovae

Nuclear shell burning in the final stages of the lives of massive stars is accompanied by strong turbulent convection. The resulting fluctuations aid supernova explosion by amplifying the non-radial flow in the post-shock region. In this work, we investigate the physical mechanism behind this amplification using a linear perturbation theory. We model the shock wave as a one-dimensional planar discontinuity and consider its interaction with vorticity and entropy perturbations in the upstream flow. We find that, as the perturbations cross the shock, their total turbulent kinetic energy is amplified by a factor of $\sim\!2$, while the average linear size of turbulent eddies decreases by about the same factor. These values are not sensitive to the parameters of the upstream turbulence and the nuclear dissociation efficiency at the shock. Finally, we discuss the implication of our results for the supernova explosion mechanism. We show that the upstream perturbations can decrease the critical neutrino luminosity for producing explosion by several percent.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The One-Armed Spiral Instability in Neutron Star Mergers and its Detectability in Gravitational Waves

We study the development and saturation of the $m=1$ one-armed spiral instability in remnants of binary neutron star mergers by means of high-resolution long-term numerical relativity simulations. Our results suggest that this instability is a generic outcome of neutron stars mergers in astrophysically relevant configurations; including both "stiff" and "soft" nuclear equations of state. We find that, once seeded at merger, the $m=1$ mode saturates within $\sim 10\ \mathrm{ms}$ and persists over secular timescales. Gravitational waves emitted by the $m=1$ instability have a peak frequency around $1-2\ \mathrm{kHz}$ and, if detected, could be used to constrain the equation of state of neutron stars. We construct hybrid waveforms spanning the entire Advanced LIGO band by combining our high-resolution numerical data with state-of-the-art effective-one-body waveforms including tidal effects. We use the complete hybrid waveforms to study the detectability of the one-armed spiral instability for both Advanced LIGO and the Einstein Telescope. We conclude that the one-armed spiral instability is not an efficient gravitational wave emitter. Its observation by current generation detectors is unlikely and will require third-generation interferometers.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-Order Numerical-Relativity Simulations of Binary Neutron Stars

We report simulations of the inspiral and merger of binary neutron stars performed with \texttt{WhiskyTHC}, the first of a new generation of numerical relativity codes employing higher than second-order methods for both the spacetime and the hydrodynamic evolution. We find that the use of higher-order schemes improves substantially the quality of the gravitational waveforms extracted from the simulations when compared to those computed using traditional second-order schemes. The reduced de-phasing and the faster convergence rate allow us to estimate the phase evolution of the gravitational waves emitted, as well as the magnitude of finite-resolution effects, without the need of phase- or time-alignments or rescalings of the waves, as sometimes done in other works. Furthermore, by using an additional unpublished simulation at very high resolution, we confirm the robustness of our high convergence order of $3.2$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Implicit large eddy simulations of anisotropic weakly compressible turbulence with application to core-collapse supernovae

(Abridged) In the implicit large eddy simulation (ILES) paradigm, the dissipative nature of high-resolution shock-capturing schemes is exploited to provide an implicit model of turbulence. Recent 3D simulations suggest that turbulence might play a crucial role in core-collapse supernova explosions, however the fidelity with which turbulence is simulated in these studies is unclear. Especially considering that the accuracy of ILES for the regime of interest in CCSN, weakly compressible and strongly anisotropic, has not been systematically assessed before. In this paper we assess the accuracy of ILES using numerical methods most commonly employed in computational astrophysics by means of a number of local simulations of driven, weakly compressible, anisotropic turbulence. We report a detailed analysis of the way in which the turbulent cascade is influenced by the numerics. Our results suggest that anisotropy and compressibility in CCSN turbulence have little effect on the turbulent kinetic energy spectrum and a Kolmogorov $k^{-5/3}$ scaling is obtained in the inertial range. We find that, on the one hand, the kinetic energy dissipation rate at large scales is correctly captured even at relatively low resolutions, suggesting that very high effective Reynolds number can be achieved at the largest scales of the simulation. On the other hand, the dynamics at intermediate scales appears to be completely dominated by the so-called bottleneck effect, \ie the pile up of kinetic energy close to the dissipation range due to the partial suppression of the energy cascade by numerical viscosity. An inertial range is not recovered until the point where relatively high resolution $\sim 512^3$, which would be difficult to realize in global simulations, is reached. We discuss the consequences for CCSN simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

High-Order Fully General-Relativistic Hydrodynamics: new Approaches and Tests

We present a new approach for achieving high-order convergence in fully general-relativistic hydrodynamic simulations. The approach is implemented in WhiskyTHC, a new code that makes use of state-of-the-art numerical schemes and was key in achieving, for the first time, higher than second-order convergence in the calculation of the gravitational radiation from inspiraling binary neutron stars Radice et al. (2013). Here, we give a detailed description of the algorithms employed and present results obtained for a series of classical tests involving isolated neutron stars. In addition, using the gravitational-wave emission from the late inspiral and merger of binary neutron stars, we make a detailed comparison between the results obtained with the new code and those obtained when using standard second-order schemes commonly employed for matter simulations in numerical relativity. We find that even at moderate resolutions and for binaries with large compactness, the phase accuracy is improved by a factor 50 or more.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A New Spherical Harmonics Scheme for Multi-Dimensional Radiation Transport I: Static Matter Configurations

Recent work by McClarren & Hauck [29] suggests that the filtered spherical harmonics method represents an efficient, robust, and accurate method for radiation transport, at least in the two-dimensional (2D) case. We extend their work to the three-dimensional (3D) case and find that all of the advantages of the filtering approach identified in 2D are present also in the 3D case. We reformulate the filter operation in a way that is independent of the timestep and of the spatial discretization. We also explore different second- and fourth-order filters and find that the second-order ones yield significantly better results. Overall, our findings suggest that the filtered spherical harmonics approach represents a very promising method for 3D radiation transport calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Beyond second-order convergence in simulations of binary neutron stars in full general-relativity

Despite the recent rapid progress in numerical relativity, a convergence order less than the second has so far plagued codes solving the Einstein-Euler system of equations. We report simulations of the inspiral of binary neutron stars in quasi-circular orbits computed with a new code employing high-order, high-resolution shock-capturing, finite-differencing schemes that, for the first time, go beyond the second-order barrier. In particular, without any tuning or alignment, we measure a convergence order above three both in the phase and in the amplitude of the gravitational waves. Because the new code is able to calculate waveforms with very small phase errors already at modest resolutions, we are able to obtain accurate estimates of tidal effects in the inspiral that are essentially free from the large numerical viscosity typical of lower-order methods, and even for the challenging large compactness and small-deformability binary considered here. We find a remarkable agreement between our Richardson-extrapolated waveform and the one from the tidally corrected post-Newtonian (PN) Taylor-T4 model, with a de-phasing smaller than 0.2 radians during the seven orbits of the inspiral and up to the contact point. Because our results can be used reliably to assess the validity of the PN or other approximations at frequencies significantly larger than those considered so far in the literature, they seem to exclude at these compactnesses significant tidal amplifications from next-to-next-to-leading--order terms in the PN expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Universality and intermittency in relativistic turbulent flows of a hot plasma

With the aim of determining the statistical properties of relativistic turbulence and unveiling novel and non-classical features, we resent the results of direct numerical simulations of driven turbulence in an ultrarelativistic hot plasma using high-order numerical schemes. We study the statistical properties of flows with average Mach number ranging from $\sim 0.4$ to $\sim 1.7$ and with average Lorentz factors up to $\sim 1.7$. We find that flow quantities, such as the energy density or the local Lorentz factor, show large spatial variance even in the subsonic case as compressibility is enhanced by relativistic effects. The velocity field is highly intermittent, but its power-spectrum is found to be in good agreement with the predictions of the classical theory of Kolmogorov. Overall, our results indicate that relativistic effects are able to significantly enhance the intermittency of the flow and affect the high-order statistics of the velocity field, while leaving unchanged the low-order statistics, which instead appear to be universal and in good agreement with the classical Kolmogorov theory. To the best of our knowledge, these are the most accurate simulations of driven relativistic turbulence to date.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

THC: a new high-order finite-difference high-resolution shock-capturing code for special-relativistic hydrodynamics

We present THC: a new high-order flux-vector-splitting code for Newtonian and special-relativistic hydrodynamics designed for direct numerical simulations of turbulent flows. Our code implements a variety of different reconstruction algorithms, such as the popular weighted essentially non oscillatory and monotonicity-preserving schemes, or the more specialised bandwidth-optimised WENO scheme that has been specifically designed for the study of compressible turbulence. We show the first systematic comparison of these schemes in Newtonian physics as well as for special-relativistic flows. In particular we will present the results obtained in simulations of grid-aligned and oblique shock waves and nonlinear, large-amplitude, smooth adiabatic waves. We will also discuss the results obtained in classical benchmarks such as the double-Mach shock reflection test in Newtonian physics or the linear and nonlinear development of the relativistic Kelvin-Helmholtz instability in two and three dimensions. Finally, we study the turbulent flow induced by the Kelvin-Helmholtz instability and we show that our code is able to obtain well-converged velocity spectra, from which we benchmark the effective resolution of the different schemes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Universality and intermittency in relativistic turbulent flows of a hot gas

With the aim of determining the statistical properties of relativistic turbulence and unveiling novel and non-classical features, we present the results of direct numerical simulations of driven turbulence in an ultrarelativistic hot plasma using high-order numerical schemes. We study the statistical properties of flows with average Mach number ranging from $\sim 0.4$ to $\sim 1.7$ and with average Lorentz factors up to $\sim 1.7$. We find that flow quantities, such as the energy density or the local Lorentz factor, show large spatial variance even in the subsonic case as compressibility is enhanced by relativistic effects. The velocity field is highly intermittent, but its power-spectrum is found to be in good agreement with the predictions of the classical theory of Kolmogorov.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Discontinuous Galerkin methods for general-relativistic hydrodynamics: formulation and application to spherically symmetric spacetimes

We have developed the formalism necessary to employ the discontinuous-Galerkin approach in general-relativistic hydrodynamics. The formalism is firstly presented in a general 4-dimensional setting and then specialized to the case of spherical symmetry within a 3+1 splitting of spacetime. As a direct application, we have constructed a one-dimensional code, EDGES, which has been used to asses the viability of these methods via a series of tests involving highly relativistic flows in strong gravity. Our results show that discontinuous Galerkin methods are able not only to handle strong relativistic shock waves but, at the same time, to attain very high orders of accuracy and exponential convergence rates in smooth regions of the flow. Given these promising prospects and their affinity with a pseudospectral solution of the Einstein equations, discontinuous Galerkin methods could represent a new paradigm for the accurate numerical modelling in relativistic astrophysics.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Critical Phenomena in Neutron Stars I: Linearly Unstable Nonrotating Models

We consider the evolution in full general relativity of a family of linearly unstable isolated spherical neutron stars under the effects of very small, perturbations as induced by the truncation error. Using a simple ideal-fluid equation of state we find that this system exhibits a type-I critical behaviour, thus confirming the conclusions reached by Liebling et al. [1] for rotating magnetized stars. Exploiting the relative simplicity of our system, we are able carry out a more in-depth study providing solid evidences of the criticality of this phenomenon and also to give a simple interpretation of the putative critical solution as a spherical solution with the unstable mode being the fundamental F-mode. Hence for any choice of the polytropic constant, the critical solution will distinguish the set of subcritical models migrating to the stable branch of the models of equilibrium from the set of subcritical models collapsing to a black hole. Finally, we study how the dynamics changes when the numerically perturbation is replaced by a finite-size, resolution independent velocity perturbation and show that in such cases a nearly-critical solution can be changed into either a sub or supercritical. The work reported here also lays the basis for the analysis carried in a companion paper, where the critical behaviour in the the head-on collision of two neutron stars is instead considered [2].

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Critical Phenomena in Neutron Stars II: Head-on Collisions

We consider the head-on collision of equal-mass neutron stars boosted towards each other and we study the behavior of such systems near the threshold of black-hole formation. In particular, we confirm the previous findings by [1] that a type-I critical phenomenon can be observed by fine-tuning the initial mass of the two neutron stars. At the same time, we argue against the interpretation that the critical solution is not a perturbed spherical star and show instead that the metastable star corresponds to a (perturbed) equilibrium solution on the unstable branch of the equilibrium configurations. As a result, the head-on collision of two neutron stars near the critical threshold can be seen as a transition in the space of configurations from an initial stable solution over to a critical metastable one which can either migrate to a stable solution or collapse to a black hole. The critical exponent for this process shows a fine structure which was already observed in the case of the critical collapse of scalar fields but never before for perfect fluids.

0 citations0 reviewsNo code linkedOpen access