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

Sebastiano Bernuzzi

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

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.COastro-ph.SRphysics.comp-ph

Catalog footprint

What is connected

56works
5topics
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

56 published item(s)

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
preprint2022arXiv

Assessment of Effective-One-Body Radiation Reactions for Generic Planar Orbits

In this paper we assess the performances of different analytical prescriptions for the effective-one-body (EOB) radiation reaction along generic planar orbits using exact numerical result in the test-mass limit. We consider three prescriptions put forward in the recent literature: (i) the quasicircular prescription (QC), (ii) the QC with second post-Newtonian (2PN) order noncircular corrections (QC2PN), and (iii) the QC corrected by the noncircular Newtonian prefactor (NCN). The analytical fluxes are then compared against the exact fluxes that are computed by solving the Teukolsky equation with a test-mass source in geodesic motion. We find that the NCN prescription is the most accurate for both eccentric and hyperbolic orbits and it is in robust agreement also for large values of the eccentricity. This result carries over to the comparable masses, as we discuss for a numerical-relativity (NR) case study. We also demonstrate that, while the EOB/NR waveform unfaithfulness is a necessary check for the precision of EOB models, the direct comparison of EOB/NR fluxes is a more stringent and informative test to select the best prescription. Finally, we propose an improved radiation reaction, NCN2PN, that includes noncircular 2PN corrections, in resummed form, as a further multiplicative contribution and that is valid for any mass ratio.

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

Effective one-body model for extreme-mass-ratio spinning binaries on eccentric equatorial orbits: testing radiation reaction and waveform

We provide a systematic analysis of the multipolar gravitational waveform, energy and angular momentum fluxes emitted by a nonspinning test particle orbiting a Kerr black hole along equatorial, eccentric orbits. These quantities are computed by numerically solving the Teukolsky equation in the time domain and are then used to test and improve the radiation reaction (and waveform) of an effective-one-body (EOB) model. Eccentricity is incorporated into EOB by replacing the quasi-circular Newtonian (or leading-order) prefactors in the EOB-factorized multipolar waveform (and fluxes) with their generic counterparts. The comparison between numerical and analytical quantities is carried out over a large portion of the parameter space, notably for orbits close to the separatrix and with high eccentricities. The analytical model agrees to $\sim 1\%$ with the numerical data for orbits with moderate eccentricities ($e\lesssim 0.3$) and moderate spins ($\hat{a}\lesssim 0.5$), although this increases up to $\sim 33\%$ for large, positive, black hole spins ($\sim 0.9$) and large eccentricities ($\sim 0.9$). For moderate eccentricities, the new EOB fluxes can be used to drive the dynamics through the nonadiabatic transition from eccentric inspiral to plunge, merger and ringdown, thus providing accurate an description of the merger. We also show that the same technique can be applied to hyperbolic encounters. In general, our approach to radiation reaction for eccentric inspirals should be seen as a first step toward EOB modelization of extreme-mass-ratio-inspirals waveforms for LISA.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Effective-one-body waveforms for precessing coalescing compact binaries with post-newtonian Twist

Spin precession is a generic feature of compact binary coalescences, which leaves clear imprints in the gravitational waveforms. Building on previous work, we present an efficient time domain inspiral-merger-ringdown effective-one-body model (EOB) for precessing binary black holes, which incorporates subdominant modes beyond $\ell=2$, and the first EOB frequency domain approximant for precessing binary neutron stars. We validate our model against 99 ``short'' numerical relativity precessing waveforms, where we find median mismatches of $5\times 10^{-3}$, $7 \times 10^{-3}$ at inclinations of $0$, $π/3$, and 21 ``long'' waveforms with median mismatches of $4 \times 10^{-3}$ and $5 \times 10^{-3}$ at the same inclinations. Further comparisons against the state-of-the-art $\texttt{NRSur7dq4}$ waveform model yield median mismatches of $4\times 10^{-3}, 1.8 \times 10^{-2}$ at inclinations of $0, π/3$ for 5000 precessing configurations with the precession parameter $χ_p$ up to 0.8 and mass ratios up to 4. To demonstrate the computational efficiency of our model we apply it to parameter estimation and re-analyze the gravitational-wave events GW150914, GW190412, and GW170817.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Entropy-limited higher-order central scheme for neutron star merger simulations

Numerical relativity simulations are the only way to calculate exact gravitational waveforms from binary neutron star mergers and to design templates for gravitational-wave astronomy. The accuracy of these numerical calculations is critical in quantifying tidal effects near merger that are currently one of the main sources of uncertainty in merger waveforms. In this work, we explore the use of an entropy-based flux-limiting scheme for high-order, convergent simulations of neutron star spacetimes. The scheme effectively tracks the stellar surface and physical shocks using the residual of the entropy equation thus allowing the use of unlimited central flux schemes in regions of smooth flow. We perform the first neutron star merger simulations with such a method and demonstrate up to fourth-order convergence in the gravitational waveform phase. The scheme reduces the phase error up to a factor five when compared to state-of-the-art high-order characteristic schemes and can be employed for producing faithful tidal waveforms for gravitational-wave modelling.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Exploiting Newton-factorized, 2PN-accurate, waveform multipoles in effective-one-body models for spin-aligned noncircularized binaries

We present a new approach to factorize and resum the post-Newtonian (PN) waveform for generic equatorial motion to be used within effective-one-body (EOB) based waveform models. The new multipolar waveform factorization improves previous prescriptions in that: (i) the generic Newtonian contribution is factored out from each multipole; (ii) the circular part is factored out and resummed using standard EOB methods and (iii) the residual, 2PN-accurate, noncircular part, and in particular the tail contribution, is additionally resummed using Padé approximants. The resulting waveform is validated in the extreme-mass-ratio limit by comparisons with nine (mostly nonspinning) numerical waveforms either from eccentric inspirals, with eccentricities up to $e=0.9$, or dynamical captures . The resummation of the noncircular tail contribution is found essential to obtain excellent (${\lesssim}0.05$~rad at periastron for $e=0.9$) analytical/numerical agreement and to considerably improve the prescription with just the Newtonian prefactor. In the comparable mass case, the new 2PN waveform shows only a marginal improvement over the previous Newtonian factorization, though yielding maximal unfaithfulness $\simeq 10^{-3}$ with the 28 publicly available numerical relativity simulations with eccentricity up to $\sim 0.3$ (except for a single outlier that grazes $10^{-2}$). We finally use test-particle data to validate the waveform factorization proposed by Khalil et al.~[Phys.~Rev.~104 (2021) 2, 024046] and conclude that its amplitude can be considered reliable (though less accurate, $\sim 6\%$ fractional difference versus $1.5\%$ of our method) only up to eccentricities $\sim 0.3$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Kilohertz Gravitational Waves from Binary Neutron Star Mergers: Inference of Postmerger Signals with the Einstein Telescope

Next-generation detectors are expected to be sensitive to postmerger signals from binary neutron star coalescences and thus to directly probe the remnant dynamics. We investigate the scientific potential of postmerger detections with the Einstein Telescope using full Bayesian analyses with the state-of-the-art waveform model ${\tt NRPMw}$. We find that: (i) Postmerger signals with signal-to-noise ratio (SNR) ${\sim}7$ can be confidently detected with a Bayes' factor of $\log{\cal B}\simeq 5$ ($\rm e$-folded) and the posterior distributions report informative measurements already at SNR ${\sim}6$ for some noise realizations. (ii) The postmerger peak frequency $f_2$ can be confidently identified at SNR $7$ with errors of $O(1~{\rm kHz})$, that decrease below $O(100~{\rm Hz})$ for SNR 10. (iii) The remnant's time of collapse to black hole can be constrained to $O(20~{\rm ms})$ at SNR 10. However, the inference can be biased by noise fluctuationsif the latter exceed the signal's amplitude before collapse. (iv) Violations of the EOS-insentive relations for $f_2$ can be detected at SNR $\gtrsim 8$ if the frequency shifts are $\gtrsim 500~{\rm Hz}$; they can be smoking guns for EOS softening effects at extreme densities. However, the $f_2$ measurement can be significantly biased by subdominant frequency components for short-lived remnants. In these cases, an EOS softening might be better inferred from the remnant's earlier collapse.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Kilohertz Gravitational Waves From Binary Neutron Star Mergers: Numerical-relativity Informed Postmerger Model

We present ${\tt NRPMw}$, an analytical model of gravitational-waves from neutron star merger remnants informed using 618 numerical relativity (NR) simulations. ${\tt NRPMw}$ is designed in the frequency domain using a combination of complex Gaussian wavelets. The wavelet's parameters are calibrated to equations of state (EOS) insensitive relations from NR data. The NR simulations are computed with 21 EOS (7 of which are finite-temperature microphysical models, and 3 of which contain quark phase transitions or hyperonic degrees of freedom) and span total binary masses $M\in[2.4,3.4]~{\rm M}_\odot$, mass ratios up to $q=2$, and (nonprecessing) dimensionless spins magnitudes up to ${0.2}$. The theoretical uncertainties of the EOS-insensitive relations are incorporated in ${\tt NRPMw}$ using recalibration parameters that enhance the flexibility and accuracy of the model. ${\tt NRPMw}$ is NR-faithful with fitting factors ${\gtrsim}0.9$ computed on an independent validation set of 102 simulations.

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

New Avenue for Accurate Analytical Waveforms and Fluxes for Eccentric Compact Binaries

We introduce a new paradigm for constructing accurate analytic waveforms (and fluxes) for eccentric compact binaries. Our recipe builds on the standard Post-Newtonian (PN) approach but (i) retains implicit time-derivatives of the phase space variables in the instantaneous part of the noncircular waveform, and then (ii) suitably factorizes and resums this partly PN-implicit waveform using effective-one-body (EOB) procedures. We test our prescription against the exact results obtained by solving the Teukolsky equation with a test-mass source orbiting a Kerr black hole, and compare the use of the exact vs PN equations of motion for the time derivatives computation. Focusing only on the quadrupole contribution, we find that the use of the exact equations of motion yields an analytical/numerical agreement of the (averaged) angular momentum fluxes that is improved by $40\%$ with respect to previous work, with $4.5\%$ fractional difference for eccentricity $e=0.9$ and black hole dimensionless spin $-0.9\leq \hat{a}\leq +0.9$. We also find a remarkable convergence trend between Newtonian, 1PN and 2PN results. Our approach carries over to the comparable mass case using the resummed EOB equations of motion and paves the way to faithful EOB-based waveform model for long-inspiral eccentric binaries for current and future gravitational wave detectors.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

New Horizons for Fundamental Physics with LISA

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas.

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

Numerical-relativity validation of effective-one-body waveforms in the intermediate-mass-ratio regime

One of the open problems in developing binary black hole (BBH) waveforms for gravitational wave astronomy is to model the intermediate mass ratio regime and connect it to the extreme mass ratio regime. A natural approach is to employ the effective one body (EOB) approach to the two-body dynamics that, by design, can cover the entire mass ratio range and naturally incorporates the extreme mass ratio limit. Here we use recently obtained numerical relativity (NR) data with mass ratios $m_1/m_2=(7,15,\,32,\,64,\,128)$ to test the accuracy of the state-of-the-art EOB model TEOBResumS in the intermediate mass ratio regime. We generally find an excellent EOB/NR consistency around merger and ringdown for all mass ratios and for all available subdominant multipoles, except for the $\ell=m=5$ one. This mode can be crucially improved using the new large-mass ratio NR data of this paper. The EOB/NR inspirals are also consistent with the estimated NR uncertainties. We also use several NR datasets taken by different public catalogs to probe the universal behavior of the multipolar hierarchy of waveform amplitudes at merger, that smoothly connects the equal-mass BBH to the test-mass result. Interestingly, the universal behavior is strengthened if the nonoscillatory memory contribution is included in the NR waveform. Future NR simulations with improved accuracy will be necessary to further probe, and possibly quantitatively refine, the TEOBResumS transition from late inspiral to plunge in the intermediate mass ratio regime.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A hybrid post-Newtonian -- effective-one-body scheme for spin-precessing compact-binary waveforms

We introduce \texttt{TEOBResumSP}: an efficient, accurate hybrid scheme for generating gravitational waveforms from spin-precessing compact binaries. The precessing waveforms are generated via the established technique of Euler rotating the non-precessing \texttt{TEOBResumS} waveforms from a precessing frame to an inertial frame. We obtain the Euler angles by solving the post-Newtonian precession equations expanded to second post-Newtonian order. Current version of \texttt{TEOBResumSP} produces precessing waveforms through the inspiral phase up to the onset of the merger. We compare \texttt{TEOBResumSP} to current state-of-the-art precessing approximants \texttt{NRSur7dq4}, \texttt{SEOBNRv4PHM}, and \texttt{IMRPhenomPv3HM} for 200 cases of precessing compact binary inspirals with orbital inclinations up to 90 degrees, mass ratios up to four, and the effective precession parameter $χ_p$ up to 0.75. We further provide an extended comparison with \texttt{SEOBNRv4PHM} involving 1030 more inspirals with $χ_p\le 1$ and mass ratios up to 10. We find that 91\% of the \texttt{TEOBResumSP}-\texttt{NRSur7dq4} matches, 85\% of the \texttt{TEOBResumSP}-\texttt{SEOBNRv4PHM} matches, and 77\% of the \texttt{TEOBResumSP}-\texttt{IMRPhenomPv3HM} matches are greater than $0.965$. Most disagreements occur for large mass ratios and $χ_p \gtrsim 0.6$. We identify the mismatch of the \emph{non}-precessing $(2,1)$ mode as one of the leading causes of disagreements. We also introduce a new parameter, $χ_{\perp,\text{max}}$, to measure the strength of precession and hint that the mismatch between the above approximants shows an exponential dependence on $χ_{\perp,\text{max}}$ though this requires further study. Our results indicate that \texttt{TEOBResumSP} is on its way to becoming a robust precessing approximant to be employed in the parameter estimation of generic-spin compact binaries.

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

Black-hole remnants from black-hole--neutron-star mergers

Observations of gravitational waves and their electromagnetic counterparts may soon uncover the existence of coalescing compact binary systems formed by a stellar-mass black hole and a neutron star. These mergers result in a remnant black hole, possibly surrounded by an accretion disk. The mass and spin of the remnant black hole depend on the properties of the coalescing binary. We construct a map from the binary components to the remnant black hole using a sample of numerical-relativity simulations of different mass ratios $q$, (anti-)aligned dimensionless spins of the black hole $a_{\rm BH}$, and several neutron star equations of state. Given the binary total mass, the mass and spin of the remnant black hole can therefore be determined from the three parameters $(q,a_{\rm BH},Λ)$, where $Λ$ is the tidal deformability of the neutron star. Our models also incorporate the binary black hole and test-mass limit cases and we discuss a simple extension for generic black hole spins. We combine the remnant characterization with recent population synthesis simulations for various metallicities of the progenitor stars that generated the binary system. We predict that black-hole-neutron-star mergers produce a population of remnant black holes with masses distributed around $7M_\odot$ and $9M_\odot$. For isotropic spin distributions, nonmassive accretion disks are favoured: no bright electromagnetic counterparts are expected in such 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

Machine Learning Gravitational Waves from Binary Black Hole Mergers

We apply machine learning methods to build a time-domain model for gravitational waveforms from binary black hole mergers, called mlgw. The dimensionality of the problem is handled by representing the waveform's amplitude and phase using a principal component analysis. We train mlgw on about $\mathcal{O}(10^3)$ TEOBResumS and SEOBNRv4 effective-one-body waveforms with mass ratios $q\in[1,20]$ and aligned dimensionless spins $s\in[-0.80,0.95]$. The resulting models are faithful to the training sets at the ${\sim}10^{-3}$ level (averaged on the parameter space). The speed up for a single waveform generation is a factor 10 to 50 (depending on the binary mass and initial frequency) for TEOBResumS and approximately an order of magnitude more for SEOBNRv4. Furthermore, mlgw provides a closed form expression for the waveform and its gradient with respect to the orbital parameters; such an information might be useful for future improvements in GW data analysis. As demonstration of the capabilities of mlgw to perform a full parameter estimation, we re-analyze the public data from the first GW transient catalog (GWTC-1). We find broadly consistent results with previous analyses at a fraction of the cost, although the analysis with spin aligned waveforms gives systematic larger values of the effective spins with respect to previous analyses with precessing waveforms. Since the generation time does not depend on the length of the signal, our model is particularly suitable for the analysis of the long signals that are expected to be detected by third-generation detectors. Future applications include the analysis of waveform systematics and model selection in parameter estimation.

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
preprint2016arXiv

Gravitational waveforms for neutron star binaries from binary black hole simulations

Gravitational waves from binary neutron star (BNS) and black hole/neutron star (BHNS) inspirals are primary sources for detection by the Advanced Laser Interferometer Gravitational-Wave Observatory. The tidal forces acting on the neutron stars induce changes in the phase evolution of the gravitational waveform, and these changes can be used to constrain the nuclear equation of state. Current methods of generating BNS and BHNS waveforms rely on either computationally challenging full 3D hydrodynamical simulations or approximate analytic solutions. We introduce a new method for computing inspiral waveforms for BNS/BHNS systems by adding the post-Newtonian (PN) tidal effects to full numerical simulations of binary black holes (BBHs), effectively replacing the nontidal terms in the PN expansion with BBH results. Comparing a waveform generated with this method against a full hydrodynamical simulation of a BNS inspiral yields a phase difference of $<1$ radian over $\sim 15$ orbits. The numerical phase accuracy required of BNS simulations to measure the accuracy of the method we present here is estimated as a function of the tidal deformability parameter $λ$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gravitational waveforms from binary neutron star mergers with high-order WENO schemes in numerical relativity

The theoretical modeling of gravitational waveforms from binary neutron star mergers requires precise numerical relativity simulations. Assessing convergence of the numerical data and building the error budget is currently challenging due to the low accuracy of general-relativistic hydrodynamics schemes and to the grid resolutions that can be employed in $(3+1)$-dimensional simulations. In this work, we explore the use of high-order weighted-essentially-nonoscillatory (WENO) schemes in neutron star merger simulations and investigate the accuracy of the waveforms obtained with such methods. We find that high-order WENO schemes can be robustly employed for simulating the inspiral-merger phase and they significantly improve the assessment of the waveform's error budget with respect to finite-volume methods. High-order WENO schemes can be thus efficiently used for high-quality waveforms production, also in future large-scale investigations of the binary parameter space.

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

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

Binary Neutron Stars with Generic Spin, Eccentricity, Mass ratio, and Compactness - Quasi-equilibrium Sequences and First Evolutions

Information about the last stages of a binary neutron star inspiral and the final merger can be extracted from quasi-equilibrium configurations and dynamical evolutions. In this article, we construct quasi-equilibrium configurations for different spins, eccentricities, mass ratios, compactnesses, and equations of state. For this purpose we employ the SGRID code, which allows us to construct such data in previously inaccessible regions of the parameter space. In particular, we consider spinning neutron stars in isolation and in binary systems; we incorporate new methods to produce highly eccentric and eccentricity reduced data; we present the possibility of computing data for significantly unequal-mass binaries; and we create equal-mass binaries with individual compactness up to 0.23. As a proof of principle, we explore the dynamical evolution of three new configurations. First, we simulate a $q=2.06$ mass ratio which is the highest mass ratio for a binary neutron star evolved in numerical relativity to date. We find that mass transfer from the companion star sets in a few revolutions before merger and a rest mass of $\sim10^{-2}M_\odot$ is transferred between the two stars. This configuration also ejects a large amount of material during merger, imparting a substantial kick to the remnant. Second, we simulate the first merger of a precessing binary neutron star. We present the dominant modes of the gravitational waves for the precessing simulation, where a clear imprint of the precession is visible in the (2,1) mode. Finally, we quantify the effect of an eccentricity reduction procedure on the gravitational waveform. The procedure improves the waveform quality and should be employed in future precision studies, but also other errors, notably truncation errors, need to be reduced in order for the improvement due to eccentricity reduction to be effective. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modeling the complete gravitational wave spectrum of neutron star mergers

In the context of neutron star mergers, we study the gravitational wave spectrum of the merger remnant using numerical relativity simulations. Postmerger spectra are characterized by a main peak frequency $f_2$ related to the particular structure and dynamics of the remnant hot hypermassive neutron star. We show that $f_2$ is correlated with the tidal coupling constant $κ^T_2$ that characterizes the binary tidal interactions during the late-inspiral--merger. The relation $f_2(κ^T_2)$ depends very weakly on the binary total mass, mass-ratio, equation of state, and thermal effects. This observation opens up the possibility of developing a model of the gravitational spectrum of every merger unifying the late-inspiral and postmerger descriptions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modeling the dynamics of tidally-interacting binary neutron stars up to merger

The data analysis of the gravitational wave signals emitted by coalescing neutron star binaries requires the availability of an accurate analytical representation of the dynamics and waveforms of these systems. We propose an effective-one-body (EOB) model that describes the general relativistic dynamics of neutron star binaries from the early inspiral up to merger. Our EOB model incorporates an enhanced attractive tidal potential motivated by recent analytical advances in the post-Newtonian and gravitational self-force description of relativistic tidal interactions. No fitting parameters are introduced for the description of tidal interaction in the late, strong-field dynamics. We compare the model energetics and the gravitational wave phasing with new high-resolution multi-orbit numerical relativity simulations of equal-mass configurations with different equations of state. We find agreement within the uncertainty of the numerical data for all configurations. Our model is the first semi-analytical model which captures the tidal amplification effects close to merger. It thereby provides the most accurate analytical representation of binary neutron star dynamics and waveforms currently available.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Numerical relativity simulations of neutron star merger remnants using conservative mesh refinement

We study equal and unequal-mass neutron star mergers by means of new numerical relativity simulations in which the general relativistic hydrodynamics solver employs an algorithm that guarantees mass conservation across the refinement levels of the computational mesh. We consider eight binary configurations with total mass $M=2.7\,M_\odot$, mass-ratios $q=1$ and $q=1.16$, and four different equation of states (EOSs), and one configuration with a stiff EOS, $M=2.5M_\odot$ and $q=1.5$. We focus on the post-merger dynamics and study the merger remnant, dynamical ejecta and the postmerger gravitational wave spectrum. Although most of the merger remnants form a hypermassive neutron star collapsing to a black hole+disk system on dynamical timescales, stiff EOSs can eventually produce a stable massive neutron star. Ejecta are mostly emitted around the orbital plane; favored by large mass ratios and softer EOS. The postmerger wave spectrum is mainly characterized by non-axisymmetric oscillations of the remnant. The stiff EOS configuration consisting of a $1.5M_\odot$ and a $1.0M_\odot$ neutron star shows a rather peculiar dynamics. During merger the companion star is very deformed; about~$\sim0.03M_\odot$ of rest-mass becomes unbound from the tidal tail due torque; and the merger remnant forms stable neutron star surrounded by a massive accretion disk $\sim0.3M_\odot$. Similar configurations might be particularly interesting for electromagnetic counterparts. Comparing results obtained with and without the conservative mesh refinement algorithm, we find that post-merger simulations can be affected by systematic errors if mass conservation is not enforced in the mesh refinement strategy. However, mass conservation also depends on grid details and on the artificial atmosphere setup. [abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simulations of rotating neutron star collapse with the puncture gauge: end state and gravitational waveforms

We reexamine the gravitational collapse of rotating neutron stars to black holes by new 3+1 numerical relativity simulations employing the Z4c formulation of Einstein equations, the moving puncture gauge conditions, and a conservative mesh refinement scheme for the general relativistic hydrodynamics. The end state of the collapse is compared to the vacuum spacetime resulting from the evolution of spinning puncture initial data. Using a local analysis for the metric fields, we demonstrate that the two spacetimes actually agree. Gravitational waveforms are analyzed in some detail. We connect the emission of radiation to the collapse dynamics using simplified spacetime diagrams, and discuss the similarity of the waveform structure with the one of black hole perturbation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Solving 3D relativistic hydrodynamical problems with WENO discontinuous Galerkin methods

Discontinuous Galerkin (DG) methods coupled to WENO algorithms allow high order convergence for smooth problems and for the simulation of discontinuities and shocks. In this work, we investigate WENO-DG algorithms in the context of numerical general relativity, in particular for general relativistic hydrodynamics. We implement the standard WENO method at different orders, a compact (simple) WENO scheme, as well as an alternative subcell evolution algorithm. To evaluate the performance of the different numerical schemes, we study non-relativistic, special relativistic, and general relativistic testbeds. We present the first three-dimensional simulations of general relativistic hydrodynamics, albeit for a fixed spacetime background, within the framework of WENO-DG methods. The most important testbed is a single TOV-star in three dimensions, showing that long term stable simulations of single isolated neutron stars can be obtained with WENO-DG methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A new gravitational wave generation algorithm for particle perturbations of the Kerr spacetime

We present a new approach to solve the 2+1 Teukolsky equation for gravitational perturbations of a Kerr black hole. Our approach relies on a new horizon penetrating, hyperboloidal foliation of Kerr spacetime and spatial compactification. In particular, we present a framework for waveform generation from point-particle perturbations. Extensive tests of a time domain implementation in the code {\it Teukode} are presented. The code can efficiently deliver waveforms at future null infinity. As a first application of the method, we compute the gravitational waveforms from inspiraling and coalescing black-hole binaries in the large-mass-ratio limit. The smaller mass black hole is modeled as a point particle whose dynamics is driven by an effective-one-body-resummed analytical radiation reaction force. We compare the analytical angular momentum loss to the gravitational wave angular momentum flux. We find that higher-order post-Newtonian corrections are needed to improve the consistency for rapidly spinning binaries. Close to merger, the subdominant multipolar amplitudes (notably the $m=0$ ones) are enhanced for retrograde orbits with respect to prograde ones. We argue that this effect mirrors nonnegligible deviations from circularity of the dynamics during the late-plunge and merger phase. We compute the gravitational wave energy flux flowing into the black hole during the inspiral using a time-domain formalism proposed by Poisson. Finally, a self-consistent, iterative method to compute the gravitational wave fluxes at leading-order in the mass of the particle is presented. For a specific case study with $\hat{a}$=0.9, a simulation that uses the consistent flux differs from one that uses the analytical flux by $\sim35$ gravitational wave cycles over a total of about $250$ cycles. In this case the horizon absorption accounts for about $+5$ gravitational wave cycles.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mergers of binary neutron stars with realistic spin

Simulations of binary neutron stars have seen great advances in terms of physical detail and numerical quality. However, the spin of the neutron stars, one of the simplest global parameters of binaries, remains mostly unstudied. We present the first, fully nonlinear general relativistic dynamical evolutions of the last three orbits for constraint satisfying initial data of spinning neutron star binaries, with astrophysically realistic spins aligned and anti-aligned to the orbital angular momentum. The initial data is computed with the constant rotational velocity approach. The dynamics of the systems is analyzed in terms of gauge-invariant binding energy vs. orbital angular momentum curves. By comparing to a binary black hole configuration we can estimate the different tidal and spin contributions to the binding energy for the first time. First results on the gravitational wave forms are presented. The phase evolution during the orbital motion is significantly affected by spin-orbit interactions, leading to delayed or early mergers. Furthermore, a frequency shift in the main emission mode of the hyper massive neutron star is observed. Our results suggest that a detailed modeling of merger waveforms requires the inclusion of spin, even for the moderate magnitudes observed in binary neutron star systems.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Quasiuniversal properties of neutron star mergers

Binary neutron star mergers are studied using nonlinear 3+1 numerical relativity simulations and the analytical effective-one-body (EOB) model. The EOB model predicts quasiuniversal relations between the mass-rescaled gravitational wave frequency and the binding energy at the moment of merger, and certain dimensionless binary tidal coupling constants depending on the stars Love numbers, compactnesses and the binary mass ratio. These relations are quasiuniversal in the sense that, for a given value of the tidal coupling constant, they depend significantly neither on the equation of state nor on the mass ratio, though they do depend on stars spins. The spin dependence is approximately linear for small spins aligned with the orbital angular momentum. The quasiuniversality is a property of the conservative dynamics; nontrivial relations emerge as the binary interaction becomes tidally dominated. This analytical prediction is qualitatively consistent with new, multi-orbit numerical relativity results for the relevant case of equal-mass irrotational binaries. Universal relations are thus expected to characterize neutron star mergers dynamics. In the context of gravitational wave astronomy, these universal relations may be used to constrain the neutron star equation of state using waveforms that model the merger accurately.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The antikick strikes back: recoil velocities for nearly-extremal binary black hole mergers in the test-mass limit

Gravitational waves emitted from a generic binary black-hole merger carry away linear momentum anisotropically, resulting in a gravitational recoil, or "kick", of the center of mass. For certain merger configurations the time evolution of the magnitude of the kick velocity has a local maximum followed by a sudden drop. Perturbative studies of this "antikick" in a limited range of black hole spins have found that the antikick decreases for retrograde orbits as a function of negative spin. We analyze this problem using a recently developed code to evolve gravitational perturbations from a point-particle in Kerr spacetime driven by an effective-one-body resummed radiation reaction force at linear order in the mass ratio $ν\ll 1$. Extending previous studies to nearly-extremal negative spins, we find that the well-known decrease of the antikick is overturned and, instead of approaching zero, the antikick increases again to reach $Δv/(cν^{2})=3.37\times10^{-3}$ for dimensionless spin $\hat{a}=-0.9999$. The corresponding final kick velocity is $v_{end}/(cν^{2})=0.076$. This result is connected to the nonadiabatic character of the emission of linear momentum during the plunge. We interpret it analytically by means of the quality factor of the flux to capture quantitatively the main properties of the kick velocity. The use of such quality factor of the flux does not require trajectories nor horizon curvature distributions and should therefore be useful both in perturbation theory and numerical relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Eccentric binary neutron star mergers

Neutron star binaries offer a rich phenomenology in terms of gravitational waves and merger remnants. However, most general relativistic studies have been performed for nearly circular binaries, with the exception of head-on collisions. We present the first numerical relativity investigation of mergers of eccentric equal-mass neutron-star binaries that probes the regime between head-on and circular. In addition to gravitational waves generated by the orbital motion, we find that the signal also contains a strong component due to stellar oscillations (f-modes) induced by tidal forces, extending a classical result for Newtonian binaries. The merger can lead to rather massive disks on the order of 10% of the total initial mass.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Error-analysis and comparison to analytical models of numerical waveforms produced by the NRAR Collaboration

The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binary's total mass is ~100-200 solar masses, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios <= 4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Improved effective-one-body description of coalescing nonspinning black-hole binaries and its numerical-relativity completion

We improve the effective-one-body (EOB) description of nonspinning coalescing black hole binaries by incorporating several recent analytical advances, notably: (i) logarithmic contributions to the conservative dynamics; (ii) resummed horizon-absorption contribution to the orbital angular momentum loss; and (iii) a specific radial component of the radiation reaction force implied by consistency with the azimuthal one. We then complete this analytically improved EOB model by comparing it to accurate numerical relativity (NR) simulations performed by the Caltech-Cornell-CITA group for mass ratios $q=(1,2,3,4,6)$. In particular, the comparison to NR data allows us to determine with high-accuracy ($\sim 10^{-4}$) the value of the main EOB radial potential: $A(u;\,ν)$, where $u=GM/(R c^2)$ is the inter-body gravitational potential and $ν=q/(q+1)^2$ is the symmetric mass ratio. We introduce a new technique for extracting from NR data an intrinsic measure of the phase evolution, ($Q_ω(ω)$ diagnostics). Aligning the NR-completed EOB quadrupolar waveform and the NR one at low frequencies, we find that they keep agreeing (in phase and amplitude) within the NR uncertainties throughout the evolution for all mass ratios considered. We also find good agreement for several subdominant multipoles without having to introduce and tune any extra parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical solution of the 2+1 Teukolsky equation on a hyperboloidal and horizon penetrating foliation of Kerr and application to late-time decays

In this work we present a formulation of the Teukolsky equation for generic spin perturbations on the hyperboloidal and horizon penetrating foliation of Kerr recently proposed by Racz and Toth. An additional, spin-dependent rescaling of the field variable can be used to achieve stable, long-term, and accurate time-domain evolutions of generic spin perturbations. As an application (and a severe numerical test), we investigate the late-time decays of electromagnetic and gravitational perturbations at the horizon and future null infinity by means of 2+1 evolutions. As initial data we consider four combinations of (non-)stationary and (non-)compact-support initial data with a pure spin-weighted spherical harmonic profile. We present an extensive study of late time decays of axisymmetric perturbations. We verify the power-law decay rates predicted analytically, together with a certain "splitting" behaviour of the power-law exponent. We also present results for non-axisymmetric perturbations. In particular, our approach allows to study the behaviour of the late time decays of gravitational fields for nearly extremal and extremal black holes. For rapid rotation we observe a very prolonged, weakly damped, quasi-normal-mode phase. For extremal rotation the field at future null infinity shows an oscillatory behaviour decaying as the inverse power of time, while at the horizon it is amplified by several orders of magnitude over long time scales. This behaviour can be understood in terms of the superradiance cavity argument.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Transient Gravitational-Wave Sky

Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Accuracy of numerical relativity waveforms from binary neutron star mergers and their comparison with post-Newtonian waveforms

We present numerical relativity simulations of nine-orbit equal-mass binary neutron star covering the quasicircular late inspiral and merger. The extracted gravitational waveforms are analyzed for convergence and accuracy. Second order convergence is observed up to contact, i.e. about 3-4 cycles to merger; error estimates can be made up to this point. The uncertainties on the phase and the amplitude are dominated by truncation errors and can be minimized to 0.13 rad and less then 1%, respectively, by using several simulations and extrapolating in resolution. In the latter case finite-radius extraction uncertainties become a source of error of the same order and have to be taken into account. The waveforms are tested against accuracy standards for data analysis. The uncertainties on the waveforms are such that accuracy standards are generically not met for signal-to-noise ratios relevant for detection, except for some best cases using extrapolation from several runs. A detailed analysis of the errors is thus imperative for the use of numerical relativity waveforms from binary neutron stars in quantitative studies. The waveforms are compared with the post-Newtonian Taylor T4 approximants both for point-particle and including the analytically known tidal corrections. The T4 approximants accumulate significant phase differences of 2 rad at contact and 4 rad at merger, underestimating the influence of finite size effects. Tidal signatures in the waveforms are thus important at least during the last six orbits of the merger process.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Compact binary evolutions with the Z4c formulation

Numerical relativity simulations of compact binaries with the Z4c and BSSNOK formulations are compared. The Z4c formulation is advantageous in every case considered. In simulations of non-vacuum spacetimes the constraint violations due to truncation errors are between one and three orders of magnitude lower in the Z4c evolutions. Improvements are also found in the accuracy of the computed gravitational radiation. For equal-mass irrotational binary neutron star evolutions we find that the absolute errors in phase and amplitude of the waveforms can be up to a factor of four smaller. The quality of the Z4c numerical data is also demonstrated by a remarkably accurate computation of the ADM mass from surface integrals. For equal-mass non-spinning binary puncture black hole evolutions we find that the absolute errors in phase and amplitude of the waveforms can be up to a factor of two smaller. In the same evolutions we find that away from the punctures the Hamiltonian constraint violation is reduced by between one and two orders of magnitude. Furthermore, the utility of gravitational radiation controlling, constraint preserving boundary conditions for the Z4c formulation is demonstrated. The evolution of spacetimes containing a single compact object confirm earlier results in spherical symmetry. The boundary conditions avoid spurious and non-convergent effects present in high resolution runs with either formulation with a more naive boundary treatment. We conclude that Z4c is preferable to BSSNOK for the numerical solution of the 3+1 Einstein equations with the puncture gauge.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraint damping for the Z4c formulation of general relativity

One possibility for avoiding constraint violation in numerical relativity simulations adopting free-evolution schemes is to modify the continuum evolution equations so that constraint violations are damped away. Gundlach et. al. demonstrated that such a scheme damps low amplitude, high frequency constraint violating modes exponentially for the Z4 formulation of General Relativity. Here we analyze the effect of the damping scheme in numerical applications on a conformal decomposition of Z4. After reproducing the theoretically predicted damping rates of constraint violations in the linear regime, we explore numerical solutions not covered by the theoretical analysis. In particular we examine the effect of the damping scheme on low-frequency and on high-amplitude perturbations of flat spacetime as well and on the long-term dynamics of puncture and compact star initial data in the context of spherical symmetry. We find that the damping scheme is effective provided that the constraint violation is resolved on the numerical grid. On grid noise the combination of artificial dissipation and damping helps to suppress constraint violations. We find that care must be taken in choosing the damping parameter in simulations of puncture black holes. Otherwise the damping scheme can cause undesirable growth of the constraints, and even qualitatively incorrect evolutions. In the numerical evolution of a compact static star we find that the choice of the damping parameter is even more delicate, but may lead to a small decrease of constraint violation. For a large range of values it results in unphysical behavior.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Horizon-absorption effects in coalescing black-hole binaries: An effective-one-body study of the non-spinning case

We study the horizon absorption of gravitational waves in coalescing, circularized, nonspinning black hole binaries. The horizon absorbed fluxes of a binary with a large mass ratio (q=1000) obtained by numerical perturbative simulations are compared with an analytical, effective-one-body (EOB) resummed expression recently proposed. The perturbative method employs an analytical, linear in the mass ratio, effective-one-body (EOB) resummed radiation reaction, and the Regge-Wheeler-Zerilli (RWZ) formalism for wave extraction. Hyperboloidal (transmitting) layers are employed for the numerical solution of the RWZ equations to accurately compute horizon fluxes up to the late plunge phase. The horizon fluxes from perturbative simulations and the EOB-resummed expression agree at the level of a few percent down to the late plunge. An upgrade of the EOB model for nonspinning binaries that includes horizon absorption of angular momentum as an additional term in the resummed radiation reaction is then discussed. The effect of this term on the waveform phasing for binaries with mass ratios spanning 1 to 1000 is investigated. We confirm that for comparable and intermediate-mass-ratio binaries horizon absorbtion is practically negligible for detection with advanced LIGO and the Einstein Telescope (faithfulness greater than or equal to 0.997).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Tidal effects in binary neutron star coalescence

We compare dynamics and waveforms from binary neutron star coalescence as computed by new long-term ($\sim 10 $ orbits) numerical relativity simulations and by the tidal effective-one-body (EOB) model including analytical tidal corrections up to second post-Newtonian order (2PN). The current analytical knowledge encoded in the tidal EOB model is found to be sufficient to reproduce the numerical data up to contact and within their uncertainties. Remarkably, no calibration of any tidal EOB free parameters is required, beside those already fitted to binary black holes data. The inclusion of 2PN tidal corrections minimizes the differences with the numerical data, but it is not possible to significantly distinguish them from the leading-order tidal contribution. The presence of a relevant amplification of tidal effects is likely to be excluded, although it can appear as a consequence of numerical inaccuracies. We conclude that the tidally-completed effective-one-body model provides nowadays the most advanced and accurate tool for modelling gravitational waveforms from binary neutron star inspiral up to contact. This work also points out the importance of extensive tests to assess the uncertainties of the numerical data, and the potential need of new numerical strategies to perform accurate simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Binary black hole coalescence in the extreme-mass-ratio limit: testing and improving the effective-one-body multipolar waveform

We discuss the properties of the effective-one-body (EOB) multipolar gravitational waveform emitted by nonspinning black-hole binaries of masses $μ$ and $M$ in the extreme-mass-ratio limit, $μ/M=ν\ll 1$. We focus on the transition from quasicircular inspiral to plunge, merger and ringdown.We compare the EOB waveform to a Regge-Wheeler-Zerilli (RWZ) waveform computed using the hyperboloidal layer method and extracted at null infinity. Because the EOB waveform keeps track analytically of most phase differences in the early inspiral, we do not allow for any arbitrary time or phase shift between the waveforms. The dynamics of the particle, common to both wave-generation formalisms, is driven by leading-order ${\cal O}(ν)$ analytically--resummed radiation reaction. The EOB and the RWZ waveforms have an initial dephasing of about $5\times 10^{-4}$ rad and maintain then a remarkably accurate phase coherence during the long inspiral ($\sim 33$ orbits), accumulating only about $-2\times 10^{-3}$ rad until the last stable orbit, i.e. $Δϕ/ϕ\sim -5.95\times 10^{-6}$. We obtain such accuracy without calibrating the analytically-resummed EOB waveform to numerical data, which indicates the aptitude of the EOB waveform for LISA-oriented studies. We then improve the behavior of the EOB waveform around merger by introducing and tuning next-to-quasi-circular corrections both in the gravitational wave amplitude and phase. For each multipole we tune only four next-to-quasi-circular parameters by requiring compatibility between EOB and RWZ waveforms at the light-ring. The resulting phase difference around merger time is as small as $\pm 0.015$ rad, with a fractional amplitude agreement of 2.5%. This suggest that next-to-quasi-circular corrections to the phase can be a useful ingredient in comparisons between EOB and numerical relativity waveforms.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Binary black hole coalescence in the large-mass-ratio limit: the hyperboloidal layer method and waveforms at null infinity

We compute and analyze the gravitational waveform emitted to future null infinity by a system of two black holes in the large mass ratio limit. We consider the transition from the quasi-adiabatic inspiral to plunge, merger, and ringdown. The relative dynamics is driven by a leading order in the mass ratio, 5PN-resummed, effective-one-body (EOB), analytic radiation reaction. To compute the waveforms we solve the Regge-Wheeler-Zerilli equations in the time-domain on a spacelike foliation which coincides with the standard Schwarzschild foliation in the region including the motion of the small black hole, and is globally hyperboloidal, allowing us to include future null infinity in the computational domain by compactification. This method is called the hyperboloidal layer method, and is discussed here for the first time in a study of the gravitational radiation emitted by black hole binaries. We consider binaries characterized by five mass ratios, $ν=10^{-2,-3,-4,-5,-6}$, that are primary targets of space-based or third-generation gravitational wave detectors. We show significative phase differences between finite-radius and null-infinity waveforms. We test, in our context, the reliability of the extrapolation procedure routinely applied to numerical relativity waveforms. We present an updated calculation of the gravitational recoil imparted to the merger remnant by the gravitational wave emission. As a self consistency test of the method, we show an excellent fractional agreement (even during the plunge) between the 5PN EOB-resummed mechanical angular momentum loss and the gravitational wave angular momentum flux computed at null infinity. New results concerning the radiation emitted from unstable circular orbits are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Constraint preserving boundary conditions for the Z4c formulation of general relativity

We discuss high order absorbing constraint preserving boundary conditions for the Z4c formulation of general relativity coupled to the moving puncture family of gauges. We are primarily concerned with the constraint preservation and absorption properties of these conditions. In the frozen coefficient approximation, with an appropriate first order pseudo-differential reduction, we show that the constraint subsystem is boundary stable on a four dimensional compact manifold. We analyze the remainder of the initial boundary value problem for a spherical reduction of the Z4c formulation with a particular choice of the puncture gauge. Numerical evidence for the efficacy of the conditions is presented in spherical symmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational-wave extraction from neutron-star oscillations

We compare different gravitational-wave extraction methods used in three-dimensional nonlinear simulations against linear simulations of perturbations of spherical spacetimes with matter. We present results from fully general-relativistic simulations of a system composed by an oscillating and non-rotating star emitting gravitational radiation. Results about the onset of non-linear effects are also shown.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical relativity simulations of binary neutron stars

We present a new numerical relativity code designed for simulations of compact binaries involving matter. The code is an upgrade of the BAM code to include general relativistic hydrodynamics and implements state-of-the-art high-resolution-shock-capturing schemes on a hierarchy of mesh refined Cartesian grids with moving boxes. We test and validate the code in a series of standard experiments involving single neutron star spacetimes. We present test evolutions of quasi-equilibrium equal-mass irrotational binary neutron star configurations in quasi-circular orbits which describe the late inspiral to merger phases. Neutron star matter is modeled as a zero-temperature fluid; thermal effects can be included by means of a simple ideal-gas prescription. We analyze the impact that the use of different values of damping parameter in the Gamma-driver shift condition has on the dynamics of the system. The use of different reconstruction schemes and their impact in the post-merger dynamics is investigated. We compute and characterize the gravitational radiation emitted by the system. Self-convergence of the waves is tested, and we consistently estimate error-bars on the numerically generated waveforms in the inspiral phase.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The trumpet solution from spherical gravitational collapse with puncture gauges

We investigate the stationary end-state obtained by evolving a collapsing spherical star with the gauges routinely adopted to study puncture black holes. We compare the end-state of the collapse with the trumpet solution found in the evolution of a single wormhole slice and show that the two solutions closely agree. We demonstrate that the agreement is caused by the use of the Gamma-driver shift condition, which allows the matter to fall inwards into a region of spacetime that is not resolved by the numerical grid, and which simultaneously finds the stationary coordinates of the trumpet outside the matter.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Binary black hole merger in the extreme-mass-ratio limit: a multipolar analysis

Building up on previous work, we present a new calculation of the gravitational wave (GW) emission generated during the transition from quasi-circular inspiral to plunge, merger and ringdown by a binary system of nonspinning black holes, of masses $m_1$ and $m_2$, in the extreme mass ratio limit, $m_1 m_2\ll(m_1+m_2)^2$. The relative dynamics of the system is computed {\it without making any adiabatic approximation} by using an effective one body (EOB) description, namely by representing the binary by an effective particle of mass $μ=m_1 m_2/(m_1+m_2)$ moving in a (quasi-)Schwarzschild background of mass $M=m_1+m_2$ and submitted to an $Ø(ν)$ 5PN-resummed analytical radiation reaction force, with $ν=μ/M$. The gravitational wave emission is calculated via a multipolar Regge-Wheeler-Zerilli type perturbative approach (valid in the limit $ν\ll 1$). We consider three mass ratios, $ν={10^{-2},10^{-3},10^{-4}}$,and we compute the multipolar waveform up to $\ell=8$. We estimate energy and angular momentum losses during the quasi-universal and quasi-geodesic part of the plunge phase and we analyze the structure of the ringdown. We calculate the gravitational recoil, or "kick", imparted to the merger remnant by the gravitational wave emission and we emphasize the importance of higher multipoles to get a final value of the recoil $v/(cν^2)=0.0446$. We finally show that there is an {\it excellent fractional agreement} ($\sim 10^{-3}$) (even during the plunge) between the 5PN EOB analytically-resummed radiation reaction flux and the numerically computed gravitational wave angular momentum flux. This is a further confirmation of the aptitude of the EOB formalism to accurately model extreme-mass-ratio inspirals, as needed for the future space-based LISA gravitational wave detector.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Constraint violation in free evolution schemes: comparing BSSNOK with a conformal decomposition of Z4

We compare numerical evolutions performed with the BSSNOK formulation and a conformal decomposition of a Z4-like formulation of General Relativity. The important difference between the two formulations is that the Z4 formulation has a propagating Hamiltonian constraint, whereas BSSNOK has a zero-speed characteristic variable in the constraint subsystem. In spherical symmetry we evolve both puncture and neutron star initial data. We demonstrate that the propagating nature of the Z4 constraints leads to results that compare favorably with BSSNOK evolutions, especially when matter is present in the spacetime. From the point of view of implementation the new system is a simple modification of BSSNOK.

0 citations0 reviewsNo code linkedOpen access