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.Comment: 9 pages, 9 figure

Closed-form tidal approximants for binary neutron star gravitational waveforms constructed from high-resolution numerical relativity simulations

We construct closed-form gravitational waveforms (GWs) with tidal effects for the coalescence and merger of binary neutron stars. The method relies on a new set of eccentricity-reduced and high-resolution numerical relativity (NR) simulations and is composed of three steps. First, tidal contributions to the GW phase are extracted from the time-domain NR data. Second, those contributions are employed to fix high-order coefficients in an effective and resummed post-Newtonian expression. Third, frequency-domain tidal approximants are built using the stationary phase approximation. Our tidal approximants are valid from the low frequencies to the strong-field regime and up to merger. They can be analytically added to any binary black hole GW model to obtain a binary neutron star waveform, either in the time or in the frequency domain. This work provides simple, flexible, and accurate models ready to be used in both searches and parameter estimation of binary neutron star events

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 = v « 1. We focus on the transition from quasicircular inspiral to plunge, merger, and ringdown. We compare the EOB waveform to a Regge-Wheeler-Zerilli 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 a leading-order O(v) analytically resummed radiation reaction. The EOB and the Regge-Wheeler-Zerilli waveforms have an initial dephasing of about 5 X 10^(-4) rad and maintain then a remarkably accurate phase coherence during the long inspiral (~33 orbits), accumulating only about -2 X 10^(-3) rad until the last stable orbit, i.e. ΔØ/Ø~-5.95 X 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 studies concerning the Laser Interferometer Space Antenna. We then improve the behavior of the EOB waveform around merger by introducing and tuning next-to-quasicircular corrections in both the gravitational wave amplitude and phase. For each multipole we tune only four next-to-quasicircular parameters by requiring compatibility between EOB and Regge-Wheeler-Zerilli waveforms at the light ring. The resulting phase difference around the merger time is as small as ±0.015 rad, with a fractional amplitude agreement of 2.5%. This suggest that next-to-quasicircular corrections to the phase can be a useful ingredient in comparisons between EOB and numerical-relativity waveforms

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;\,\nu)$, where $u=GM/(R c^2)$ is the inter-body gravitational potential and $\nu=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_\omega(\omega)$ 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.Comment: 42 pages, 22 figures. Improved version, to appear in Phys. Rev. D. The EOB code will be freely available at eob.ihes.f

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

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.Comment: 13 pages, 24 figure

Quasi-5.5PN TaylorF2 approximant for compact binaries: point-mass phasing and impact on the tidal polarizability inference

We derive a point-mass (nonspinning) frequency-domain TaylorF2 phasing approximant at quasi-5.5 post-Newtonian (PN) accuracy for the gravitational wave from coalescing compact binaries. The new approximant is obtained by Taylor-expanding the effective-one-body (EOB) resummed energy and and angular momentum flux along circular orbits with all the known test-particle information up to 5.5PN. The -- yet uncalculated -- terms at 4PN order and beyond entering both the energy flux and the energy are taken into account as free parameters and then set to zero. We compare the quasi-5.5PN and 3.5PN approximants against full EOB waveforms using gauge-invariant phasing diagnostics $Q_\omega=\hat\omega^2/\dot{\hat\omega}$, where $\hat\omega$ is the dimensionless gravitational-wave frequency. The quasi-5.5PN phasing is found to be systematically closer to the EOB one than the 3.5PN one. Notably, the quasi-5.5PN (3.5PN) approximant accumulates a EOB$-$PN dephasing of $\Delta\Psi^{\rm EOBPN}\sim10^{-3}$rad ($0.13$rad) up to frequency $\hat\omega \simeq 0.06$, 6 orbits to merger, ($\hat\omega \simeq 0.086$, 2 orbits to merger) for a fiducial binary neutron star system. We explore the performance of the quasi-5.5PN approximant on the measurement of the tidal polarizability parameter $\tilde\Lambda$ using injections of EOB waveforms hybridized with numerical relativity merger waveforms. We prove that the quasi-5.5PN point-mass approximant augmented with 6PN-accurate tidal terms allows one to reduce (and in many cases even eliminate) the biases in the measurement of $\tilde\Lambda$ that are instead found when the standard 3.5PN point-mass baseline is used. Methodologically, we demonstrate that the combined use of $Q_\omega$ analysis and of the Bayesian parameter estimation offers a new tool to investigate the impact of systematics on gravitational-wave inference.Comment: 13 pages, 7 figures, submitted to Phys. Rev.

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)

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