Gravitational Waveforms Research Articles

One-loop gravitational bremsstrahlung and waveforms from a heavy-mass effective field theory

Using a heavy-mass effective field theory (HEFT), we study gravitational-wave emission in the scattering of two spinless black holes or neutron stars of arbitrary masses at next-to-leading order in the Post-Minkowskian expansion. We compute the contributions to the one-loop scattering amplitude with four scalars and one graviton which are relevant to the calculation of the waveforms, also presenting expressions of classical tree-level amplitudes with four scalars and up to two radiated gravitons. The latter are obtained using a novel on-shell recursion relation for classical amplitudes with four scalars and an arbitrary number of gravitons. Our one-loop five-point amplitude is expressed in terms of a single family of master integrals with the principal value prescription for linearised massive propagators, which we evaluate using differential equations. In our HEFT approach, soft/heavy-mass expansions of complete integrands are avoided, and all hyper-classical iterations and quantum corrections are dropped at the diagrammatic level, thereby computing directly contributions to classical physics. Our result exhibits the expected factorisation of infrared divergences, the correct soft limits, and highly nontrivial cancellations of spurious poles. Finally, using our amplitude result we compute numerically the corresponding next-to-leading corrections to the spectral waveforms and the far-field time-domain waveforms using the Newman-Penrose scalar Ψ4.

Formula omitted]-mode oscillations of compact stars with realistic equations of state in dynamical spacetime

In this study, we perform full three dimensional numerical relativity simulations of non-rotating general relativistic stars. Extending the studies for polytropic equation of state, we investigate the accuracy and robustness of numerical scheme on measuring fundamental (f)-mode frequency for realistic equations of state (EoS). We use various EoS with varying range of stiffness and numerically evolve perturbed stellar models for several mass configurations (in the range of 1.2−2.0M⊙) for each of these EoS. Using the gravitational waveform obtained from the simulations we extract the f-modes of the stars. The obtained results are tested against the pre-existing perturbation methods and find good agreement. Validity and deviation of universal relations have been carried out and are compared with earlier results under Cowling approximation as well as perturbative approaches. We also show that even using perturbed single star simulations can provide good agreement with f-modes extracted from inspiral phase of binary neutron star simulations which are computationally more expensive.

Distinguishing environmental effects on binary black hole gravitational waveforms

Distinguishing environmental effects on binary black hole gravitational waveforms

Post-Keplerian waveform model for binary compact object as sources of space-based gravitational wave detector and its implications

Binary compact objects will be among the important sources for the future space-based gravitational wave detectors. Such binary compact objects include stellar massive binary black hole, binary neutron star, binary white dwarf and mixture of these compact objects. Regarding to the relatively low frequency, the gravitational interaction between the two objects of the binary is weak. Post-Newtonian approximation of general relativity is valid. Previous works about the waveform model for such binaries in the literature consider the dynamics for specific situations which involve detailed complicated matter dynamics between the two objects. We here take a different idea. We adopt the trick used in pulsar timing detection. For any gravity theories and any detailed complicated matter dynamics, the motion of the binary can always be described as a post-Keplerian expansion. And a post-Keplerian gravitational waveform model will be reduced. Instead of object masses, spins, matter’s equation of state parameters and dynamical parameters beyond general relativity, the involved parameters in our post-Keplerian waveform model are the Keplerian orbit elements and their adiabatic variations. Respect to current planning space-based gravitational wave detectors including LISA, Taiji and Tianqin, we find that the involved waveform model parameters can be well determined. And consequently the detail matter dynamics of the binary can be studied then. For binary with purely gravitational interactions, gravity theory can be constrained well.

Third post-Newtonian effective-one-body Hamiltonian in scalar-tensor and Einstein-scalar-Gauss-Bonnet gravity

We build an effective-one-body (EOB) Hamiltonian at third post-Newtonian (3PN) order in scalar-tensor (ST) and Einstein-scalar-Gauss-Bonnet (ESGB) theories of gravity. The latter is an extension of general relativity that predicts scalar hair for black holes. We start from the known two-body Lagrangian at 3PN order, and use order-reduction methods to construct its ordinary Hamiltonian counterpart. We then reduce the conservative two-body dynamics to the (nongeodesic) motion of a test particle in an effective metric by means of canonical transformations. The resulting EOB Hamiltonian is a modification of the general relativistic Hamiltonian, and already at 3PN order, it must account for nonlocal-in-time tail contributions. We include the latter beyond circular orbits and up to sixth order in the binary's orbital eccentricity. We finally calculate the orbital frequency at the innermost stable circular orbit (ISCO) of binary black holes in the shift-symmetric ESGB model. Our work extends F.L. Juli\'e and N. Deruelle [Phys. Rev. D 95, 124054 (2017)], and it is an essential step toward the accurate modeling of gravitational waveforms beyond general relativity.

Combining effective-one-body accuracy and reduced-order-quadrature speed for binary neutron star merger parameter estimation with machine learning

We present mlgw_bns, a gravitational waveform surrogate that allows for a significant improvement in the generation speed of frequency-domain waveforms for binary neutron star mergers, at a negligible cost in accuracy. This improvement is achieved by training a machine-learning model on a dataset of waveforms generated with an accurate but comparatively costlier approximant: the state-of-the-art effective-one-body model TEOBResumSPA. When coupled to a reduced-order scheme, mlgw_bns can accelerate waveform generation up to a factor of $\ensuremath{\sim}35$, outperforming all other approximants of similar accuracy. By analyzing GW170817 in realistic parameter estimation settings with our scheme, we showcase an overall speedup against TEOBResumSPA greater than an order of magnitude. Our methodology will bear a significant impact on the scientific program of next generation detectors by allowing routine usage of accurate effective-one-body models.

Probing phase transitions in neutron stars via the crust-core interfacial mode

Gravitational waves emitted from the binary neutron star (BNS) systems can carry information about the dense matter phase in these compact stars. The crust-core interfacial mode is an oscillation mode in a neutron star, and it depends mostly on the equation of state of the matter in the crust-core transition region. This mode can be resonantly excited by the tidal field of an inspiraling-in BNS system, thereby affecting the emitted gravitational waves and hence could be used to probe the equation of state in the crust-core transition region. In this work, we investigate, in detail, how the first-order phase transition inside the neutron star, if it exists, affects the properties of the crust-core interfacial mode using a Newtonian fluid perturbation theory on a general relativistic background solution of the stellar structure. Two possible types of phase transitions are considered: (1) the phase transitions happen in the fluid core but near the crust-core interface, which results in density discontinuities, and (2) the strong interaction phase transitions in the dense core (as in the conventional hybrid star case). We study how these phase transitions affect the properties of the neutron star oscillation mode excited at the interface, where there exists a shear modulus discontinuity (interfacial mode). In particular, the former phase transition has a minor effect on the mass-radius relation, and the adiabatic tidal deformability, but has the potential to significantly affect the interfacial mode frequency and thereby could be probed using gravitational waves. For the BNS systems, we discuss the possible observational signatures of these phase transitions in the gravitational waveforms and their detectability. Our work enriches the exploration of the physical properties of the crust-core interfacial mode and provides a promising method for probing the phase transition using the seismology of a compact star.

Numerical-relativity-informed effective-one-body model for black-hole–neutron-star mergers with higher modes and spin precession

We present the first effective-one-body (EOB) model for generic-spins quasi-circular black-hole--neutron-star (BHNS) inspiral-merger-ringdown gravitational waveforms (GWs). Our model is based on a new numerical-relativity (NR) informed expression of the BH remnant and its ringdown. It reproduces the NR $(\ell,m)=(2,2)$ waveform with typical phase agreement of ${\lesssim0.5}\,$rad (${\lesssim1}\,$rad) to merger (ringdown). The maximum (minimum) mismatch between the $(2,2)$ and the NR data is 4% (0.6%). Higher modes (HMs) $(2,1)$, $(3,2)$, $(3,3)$, $(4,4)$ and $(5,5)$ are included and their mismatch with the available NR waveforms are up to (down to) a 60% (1%) depending on the inclination. Phase comparison with a 16 orbit precessing simulation shows differences within the NR uncertainties. We demonstrate the applicability of the model in GW parameter estimation by perfoming the first BHNS Bayesian analysis with HMs (and non-precessing spins) of the event GW190814, together with new $(2,2)$-mode analysis of GW200105 and GW200115. For the GW190814 study, the inclusion of HMs gives tighter parameter posteriors. The Bayes factors of our analyses on this event show decisive evidence for the presence of HMs, but no clear preference for a BHNS or a binary black hole (BBH) source. Similarly, we confirm GW200105 and GW200115 show no evidence for tidal effects.

Measuring dark matter spikes around primordial black holes with Einstein Telescope and Cosmic Explorer

Future ground-based gravitational wave observatories will be ideal probes of the environments surrounding black holes with masses $1\ensuremath{-}10{M}_{\ensuremath{\bigodot}}$. Binary black hole mergers with mass ratios of order $q={m}_{2}/{m}_{1}\ensuremath{\lesssim}{10}^{\ensuremath{-}3}$ can remain in the frequency band of such detectors for months or years, enabling precision searches for modifications of their gravitational waveforms with respect to vacuum inspirals. As a concrete example of an environmental effect, we consider here a population of binary primordial black holes which are expected to be embedded in dense cold dark matter spikes. We provide a viable formation scenario for these systems compatible with all observational constraints and predict upper and lower limits on the merger rates of small-mass-ratio pairs. Given a detected signal of one such system by either Einstein Telescope or Cosmic Explorer, we show that the properties of the binary and of the dark matter spike can be measured to excellent precision with one week's worth of data, if the effect of the dark matter spike on the waveform is taken into account. However, we show that there is a risk of biased parameter inference or missing the events entirely if the effect of the predicted dark matter overdensity around these objects is not properly accounted for.

Gravitational waveform and polarizationfrom binary black hole inspiral in dynamical Chern-Simons gravity: from generation to propagation

We calculate the gravitational waveform radiated from spinning black holes (BHs) binary in dynamical Chern-Simons (dCS) gravity. The equation of motion (EOM) of the spinning binary BHs is derived based on the modified Mathisson-Papapetrou-Dixon equation for the spin-aligned circular orbits. The leading-order effects induced by the dCS theory contain spin-spin interaction and monopole-quadrupole interaction, which influences both the EOM of the binary system and corresponding gravitational waveform at the second post-Newtonian (PN) order (i.e., 2PN order). After reporting the waveforms, we investigate the polarization modes of gravitational waves (GWs) in dCS theory. None of the extra modes appears in this theory up to the considered PN order. Moreover, since the time scale of the binary merger is much smaller than that of the cosmological expansion, the parity-violating effect of the dCS theory does not appear in the process of GW generation. However, during the process of GW propagation, amplitude birefringence, a typical parity-violating effect, makes plus and cross modes convert to each other, which modifies the gravitational waveform at 1.5PN order.

Lorentz transformation of three dimensional gravitational wave tensor

Recently there has been more and more interest in the gravitational wave (GW) of moving sources. This paper introduces a Lorentz transformation problem of GWs. Although the Bondi-Metzner-Sachs (BMS) theory has in principle already included the Lorentz transformation of GWs, the transformation of the three-dimensional GW tensor has not been explicitly calculated before. Within four-dimensional spacetime, GWs have the properties of ‘boost weight zero’ and ‘spin weight 2’. This fact makes the Lorentz transformation of GWs difficult to understand. In the current paper, we adopt the traditional three-dimensional tensor description of a GW. Such a transverse-traceless tensor describes the GW freedom directly. We derive the explicit Lorentz transformation of the GW tensor. The transformation is similar to the Lorentz transformation for an electric field vector and a magnetic field vector which are three-dimensional vectors. Based on the deduced Lorentz transformation of the GW three-dimensional tensor, we can construct the gravitational waveform of a moving source with high speed if only the waveform of the corresponding rest waveform is given.

Second release of the CoRe database of binary neutron star merger waveforms

We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity (CoRe) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to . They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

Priorities in gravitational waveforms for future space-borne detectors: vacuum accuracy or environment?

ABSTRACT In preparation for future space-borne gravitational-wave (GW) detectors, should the modelling effort focus on high-precision vacuum templates or on the astrophysical environment of the sources? We perform a systematic comparison of the phase contributions caused by (1) known environmental effects in both gaseous and stellar matter backgrounds, or (2) high-order post-Newtonian (PN) terms in the evolution of mHz GW sources during the inspiral stage of massive binaries. We use the accuracy of currently available analytical waveform models as a benchmark value, finding the following trends: the largest unmodelled phase contributions are likely environmental rather than PN for binaries lighter than ∼107/(1 + z)2 M⊙, where z is the redshift. Binaries heavier than ∼108/(1 + z) M⊙ do not require more accurate inspiral waveforms due to low signal-to-noise ratios (SNRs). For high-SNR sources, environmental phase contributions are relevant at low redshift, while high-order vacuum templates are required at z ≳ 4. Led by these findings, we argue that including environmental effects in waveform models should be prioritized in order to maximize the science yield of future mHz detectors.

Complete gravitational-waveform amplitude modes for quasicircular compact binaries to the 3.5PN order

We complete the computation of the gravitational waveform amplitude for non-spinning compact binaries to the third and a half post-Newtonian (PN) order in the quasi-circular case. This is done by performing a spin-weighted spherical harmonics decomposition of the amplitude. This computation is achieved using the post-Newtonian-Multipolar-post-Minkowskian formalism. The amplitude modes are written in a PN expanded way as well as a factorized form suitable for EOB template building. Combining this work with previous ones, we provide in the supplemental materials the expressions of the modes including all non-spinning and spinning effects to the 3.5PN order for non-precessing quasi-circular binaries.

Nonlinear Effects in Black Hole Ringdown.

We report evidence for nonlinear modes in the ringdown stage of the gravitational waveform produced by the merger of two comparable-mass black holes. We consider both the coalescence of black hole binaries in quasicircular orbits and high-energy, head-on black hole collisions. The presence of nonlinear modes in the numerical simulations confirms that general-relativistic nonlinearities are important and must be considered in gravitational-wave data analysis.

Gravitational Waves and Electromagnetic Radiation from Charged Black Hole Binaries

In this manuscript, we investigate the electromagnetic radiation of a binary system of electrically charged black holes. Using the results of previous works, we compute the analytical expression for the waveform, the phase, and the Fourier transform during the inspiral phase for both the electromagnetic and gravitational radiations. To do so, we consider the quasi-circular approximation and small values for the charge-to-mass ratio in each black hole. In the case of electromagnetic radiation, we focus on the dipole contribution, but we also include the quadrupole term to complete our discussion. We found that the gravitational and electromagnetic waveforms contain two terms, and so does the Fourier transform. However, the behavior is dominated only by one of them. In the frequency-domain waves, for example, the dipole and quadrupole contributions and the gravitational wave are dominated by terms proportional to f−3/6, f−1/6, and f−7/6, respectively. As expected, the gravitational radiation and the quadrupole contribution have the same phase, in contrast to the dipole contribution. Moreover, the electromagnetic wave is more sensitive to changes in the charge-to-mass ratio than the gravitational wave.

Inspiral gravitational waveforms from compact binary systems in Horndeski gravity

In a subclass of Horndeski theories with the speed of gravity equivalent to that of light, we study gravitational radiation emitted during the inspiral phase of compact binary systems. We compute the waveform of scalar perturbations under a post-Newtonian expansion of energy-momentum tensors of pointlike particles that depend on a scalar field. This scalar mode not only gives rise to breathing and longitudinal polarizations of gravitational waves, but it is also responsible for scalar gravitational radiation in addition to energy loss associated with transverse and traceless tensor polarizations. We calculate the Fourier-transformed gravitational waveform of two tensor polarizations under a stationary phase approximation and show that the resulting waveform reduces to the one in a parametrized post-Einsteinian (ppE) formalism. The ppE parameters are directly related to a scalar charge in the Einstein frame, whose existence is crucial to allow the deviation from general relativity (GR). We apply our general framework to several concrete theories and show that a new theory of spontaneous scalarization with a higher-order scalar kinetic term leaves interesting deviations from GR that can be probed by the observations of gravitational waves emitted from neutron star--black hole binaries. If the scalar mass exceeds the order of typical orbital frequencies $\ensuremath{\omega}\ensuremath{\simeq}{10}^{\ensuremath{-}13}\text{ }\text{ }\mathrm{eV}$, which is the case for a recently proposed scalarized neutron star with a self-interacting potential, the gravitational waveform practically reduces to that in GR.

Numerical-relativity surrogate modeling with nearly extremal black-hole spins

Numerical relativity (NR) simulations of binary black hole (BBH) systems provide the most accurate gravitational wave predictions, but at a high computational cost—especially when the black holes have nearly extremal spins (i.e. spins near the theoretical upper limit) or very unequal masses. Recently, the technique of reduced order modeling has enabled the construction of ‘surrogate models’ trained on an existing set of NR waveforms. Surrogate models enable the rapid computation of the gravitational waves emitted by BBHs. Typically these models are used for interpolation to compute gravitational waveforms for BBHs with mass ratios and spins within the bounds of the training set. Because simulations with nearly extremal spins are so technically challenging, surrogate models almost always rely on training sets with only moderate spins. In this paper, we explore how well surrogate models can extrapolate to nearly extremal spins when the training set only includes moderate spins. For simplicity, we focus on one-dimensional surrogate models trained on NR simulations of BBHs with equal masses and equal, aligned spins. We assess the performance of the surrogate models at higher spin magnitudes by calculating the mismatches between extrapolated surrogate model waveforms and NR waveforms, by calculating the differences between extrapolated and NR measurements of the remnant black-hole mass, and by testing how the surrogate model improves as the training set extends to higher spins. We find that while extrapolation in this one-dimensional case is viable for current detector sensitivities, surrogate models for next-generation detectors should use training sets that extend to nearly extremal spins.

Optimizing the placement of numerical relativity simulations using a mismatch predicting neural network

Gravitational wave observations from merging compact objects are becoming commonplace, and as detectors improve and gravitational wave sources become more varied, it is increasingly important to have dense and expansive template banks of predicted gravitational waveforms. Since numerical relativity is the only way to fully solve the nonlinear merger regime of general relativity for comparably massed systems, numerical relativity simulations are critical for gravitational wave detection and analysis. These simulations are computationally expensive, with each simulation placing one point within the high dimensional parameter space of binary black hole coalescences. This makes it important to have a method of placing new simulations in ways that use our computational resources optimally while ensuring sufficient coverage of the parameter space. Accomplishing this requires predicting the impact of a new set of parameters before performing the simulation. To this effect, this paper introduces a neural network to predict the mismatch between the gravitational waves of two binary systems. Using this network, we then show how we can propose new numerical relativity simulations that will provide the most benefit. We also use the network to identify gaps in existing public catalogs and identify degeneracies in the binary black hole parameter space.

Multipole moments on the common horizon in a binary-black-hole simulation

We construct the covariantly defined multipole moments on the common horizon of an equal-mass, nonspinning, quasicircular binary-black-hole system. We see a strong correlation between these multipole moments and the gravitational waveform. We find that the multipole moments are well described by the fundamental quasinormal modes at sufficiently late times. For each nonzero multipole moment with $\ensuremath{\ell}\ensuremath{\le}6$, at least two fundamental quasinormal modes of different $\ensuremath{\ell}$ are detectable in the best model. These models provide faithful estimates of the true mass and spin of the remnant black hole. We also show that by including overtones, the $\ensuremath{\ell}=m=2$ mass multipole moment admits an excellent quasinormal-mode description at all times after the merger. This demonstrates the perhaps surprising power of perturbation theory near the merger.