post-Newtonian Equations Research Articles

Stellar dynamics of extreme-mass-ratio inspirals

Inspiral of compact stellar remnants into massive black holes (MBHs) is accompanied by the emission of gravitational waves at frequencies that are potentially detectable by space-based interferometers. Event rates computed from statistical (Fokker-Planck, Monte-Carlo) approaches span a wide range due to uncertaintities about the rate coefficients. Here we present results from direct integration of the post-Newtonian N-body equations of motion descrbing dense clusters of compact stars around Schwarzschild MBHs. These simulations embody an essentially exact (at the post-Newtonian level) treatment of the interplay between stellar dynamical relaxation, relativistic precession, and gravitational-wave energy loss. The rate of capture of stars by the MBH is found to be greatly reduced by relativistic precession, which limits the ability of torques from the stellar potential to change orbital angular momenta. Penetration of this "Schwarzschild barrier" does occasionally occur, resulting in capture of stars onto orbits that gradually inspiral due to gravitational wave emission; we discuss two mechanisms for barrier penetration and find evidence for both in the simulations. We derive an approximate formula for the capture rate, which predicts that captures would be strongly disfavored from orbits with semi-major axes below a certain value; this prediction, as well as the predicted rate, are verified in the N-body integrations. We discuss the implications of our results for the detection of extreme-mass-ratio inspirals from galactic nuclei with a range of physical properties.

Post-Newtonian cosmological dynamics of plane-parallel perturbations and back-reaction

We study the general relativistic non-linear dynamics of self-gravitating irrotational dust in a cosmological setting, adopting the comoving and synchronous gauge, where all the equations can be written in terms of the metric tensor of spatial hyper-surfaces orthogonal to the fluid flow. Performing an expansion in inverse powers of the speed of light, we obtain the post-Newtonian equations, which yield the lowest-order relativistic effects arising during the non-linear evolution. We then specialize our analysis to globally plane-parallel configurations, i.e. to the case where the initial perturbation field depends on a single coordinate. The leading order of our expansion, corresponding to the ``Newtonian background'', is the Zel'dovich approximation, which, for plane-parallel perturbations in the Newtonian limit, represents an exact solution. This allows us to find the exact analytical form for the post-Newtonian metric, thereby providing the post-Newtonian extension of the Zel'dovich solution: this accounts for some relativistic effects, such as the non-Gaussianity of primordial perturbations. An application of our solution in the context of the back-reaction proposal is eventually given, providing a post-Newtonian estimation of kinematical back-reaction, mean spatial curvature, average scale-factor and expansion rate.

Resonant trapping of stars by merging massive black hole binaries

A massive black hole binary might resonantly trap a star (e.g. a white dwarf) and the gas released by its tidal disruption might emit electromagnetic wave signals around the coalescence of the binary. With post-Newtonian equations of motion including gravitational radiation reaction, we numerically studied resonant trappings by black hole binaries with mass ratio 1/100. It is found that 2:1 (and simultaneously 4:2) mean motion resonances of the binaries would be strong and could, in principle, draw small third objects deep into relativistic regimes (e.g. ~10 Schwarzschild radii). The inclinations of the trapped objects could increase significantly and, in some cases, retrograde orbits could be realized eventually.

Post-Newtonian effects on Lagrange’s equilateral triangular solution for the three-body problem

Continuing work initiated in earlier publications [Yamada, Asada, Phys. Rev. D 82, 104019 (2010), 83, 024040 (2011)], we investigate the post-Newtonian effects on Lagrange's equilateral triangular solution for the three-body problem. For three finite masses, it is found that the equilateral triangular configuration satisfies the post-Newtonian equation of motion in general relativity, if and only if all three masses are equal. When a test mass is included, the equilateral configuration is possible for two cases: (1) one mass is finite and the other two are zero, or (2) two of the masses are finite and equal, and the third one is zero, namely a symmetric binary with a test mass. The angular velocity of the post-Newtonian equilateral triangular configuration is always smaller than the Newtonian one, provided that the masses and the side length are the same.

Characterization of the gravitational wave emission of three black holes

We study the gravitational wave emission of three compact objects using post-Newtonian (PN) equations of motion derived from the Arnowitt-Deser-Misner Hamiltonian formulation, where we include (for the first time in this context) terms up to 2.5 PN order. We perform numerical simulations of a hierarchical configuration of three compact bodies in which a binary system is perturbed by a third, lighter body initially far away from the binary. The relative importance of the different PN orders is examined. We compute the waveform in the linear regime considering mass quadrupole, current quadrupole and mass octupole contributions. Performing a spherical harmonic decomposition of the waveforms we find that from the l = 3 modes it is possible to extract information about the third body, in particular, the period, eccentricity of its orbit, and the inclination angle between the inner and outer binary orbits.

Conservative corrections to the innermost stable circular orbit (ISCO) of a Kerr black hole: A new gauge-invariant post-Newtonian ISCO condition, and the ISCO shift due to test-particle spin and the gravitational self-force

The innermost stable circular orbit (ISCO) delimits the transition from circular orbits to those that plunge into a black hole. In the test-mass limit, well-defined ISCO conditions exist for the Kerr and Schwarzschild spacetimes. In the finite-mass case, there are a large variety of ways to define an ISCO in a post-Newtonian (PN) context. Here I generalize the gauge-invariant ISCO condition of Blanchet & Iyer (2003) to the case of spinning (nonprecessing) binaries. The Blanchet-Iyer ISCO condition has two desirable and unexpected properties: (1) it exactly reproduces the Schwarzschild ISCO in the test-mass limit, and (2) it accurately approximates the recently-calculated shift in the Schwarzschild ISCO frequency due to the conservative-piece of the gravitational self-force [Barack & Sago (2009)]. The generalization of this ISCO condition to spinning binaries has the property that it also exactly reproduces the Kerr ISCO in the test-mass limit (up to the order at which PN spin corrections are currently known). The shift in the ISCO due to the spin of the test-particle is also calculated. Remarkably, the gauge-invariant PN ISCO condition exactly reproduces the ISCO shift predicted by the Papapetrou equations for a fully-relativistic spinning particle. It is surprising that an analysis of the stability of the standard PN equations of motion is able (without any form of "resummation") to accurately describe strong-field effects of the Kerr spacetime. The ISCO frequency shift due to the conservative self-force in Kerr is also calculated from this new ISCO condition, as well as from the effective-one-body Hamiltonian of Barausse & Buonanno (2010). These results serve as a useful point-of-comparison for future gravitational self-force calculations in the Kerr spacetime.

Formulation of the Post-Newtonian Equations of Motion of the Restricted Three Body Problem

In the present work the geodesic equation represents the equations of motion of the particles along the geodesics is derived. The deviation of the curved space-time metric tensor from that of the Minkowski tensor is considered as a perturbation. The quantities is expanded in powers of c-2. The equations of motion of the relativistic three body problem in the PN formalism are obtained.

Motion and gravitational wave forms of eccentric compact binaries with orbital-angular-momentum-aligned spins under next-to-leading order in spin–orbit and leading order in spin(1)–spin(2) and spin-squared couplings

A quasi-Keplerian parameterization for the solutions of second post-Newtonian (PN) accurate equations of motion for spinning compact binaries is obtained including leading order spin–spin and next-to-leading order spin–orbit interactions. Rotational deformation of the compact objects is incorporated. For arbitrary mass ratios the spin orientations are taken to be parallel or anti-parallel to the orbital angular momentum vector. The emitted gravitational wave forms are given in analytic form up to 2PN point particle, 1.5PN spin–orbit and 1PN spin–spin contributions, whereby the spins are assumed to be of 0PN order.

Symplectic integration of post-Newtonian equations of motion with spin

We present a noncanonically symplectic integration scheme tailored to numerically computing the post-Newtonian motion of a spinning black-hole binary. Using a splitting approach we combine the flows of orbital and spin contributions. In the context of the splitting, it is possible to integrate the individual terms of the spin-orbit and spin-spin Hamiltonians analytically, exploiting the special structure of the underlying equations of motion. The outcome is a symplectic, time-reversible integrator, which can be raised to arbitrary order by composition. A fourth-order version is shown to give excellent behavior concerning error growth and conservation of energy and angular momentum in long-term simulations. Favorable properties of the integrator are retained in the presence of weak dissipative forces due to radiation damping in the full post-Newtonian equations.

Final spins from the merger of precessing binary black holes

The inspiral of binary black holes is governed by gravitational radiation reaction at binary separations r < 1000 M, yet it is too computationally expensive to begin numerical-relativity simulations with initial separations r > 10 M. Fortunately, binary evolution between these separations is well described by post-Newtonian equations of motion. We examine how this post-Newtonian evolution affects the distribution of spin orientations at separations r ~ 10 M where numerical-relativity simulations typically begin. Although isotropic spin distributions at r ~ 1000 M remain isotropic at r ~ 10 M, distributions that are initially partially aligned with the orbital angular momentum can be significantly distorted during the post-Newtonian inspiral. Spin precession tends to align (anti-align) the binary black hole spins with each other if the spin of the more massive black hole is initially partially aligned (anti-aligned) with the orbital angular momentum, thus increasing (decreasing) the average final spin. Spin precession is stronger for comparable-mass binaries, and could produce significant spin alignment before merger for both supermassive and stellar-mass black hole binaries. We also point out that precession induces an intrinsic accuracy limitation (< 0.03 in the dimensionless spin magnitude, < 20 degrees in the direction) in predicting the final spin resulting from the merger of widely separated binaries.

Measuring black hole spin in OJ287

We model the binary black hole system OJ287 as a spinning primary and a non-spinning secondary. It is assumed that the primary has an accretion disk which is impacted by the secondary at specific times. These times are identified as major outbursts in the light curve of OJ287. This identification allows an exact solution of the orbit, with very tight error limits. Nine outbursts from both the historical photographic records as well as from recent photometric measurements have been used as fixed points of the solution: 1913, 1947, 1957, 1973, 1983, 1984, 1995, 2005 and 2007 outbursts. This allows the determination of eight parameters of the orbit. Most interesting of these are the primary mass of $${1.84\cdot 10^{10}M_\odot}$$ , the secondary mass $${1.40\cdot 10^{8} M_\odot}$$ , major axis precession rate 39°.1 per period, and the eccentricity of the orbit 0.70. The dimensionless spin parameter is 0.28 ± 0.01 (1 sigma). The last parameter will be more tightly constrained in 2015 when the next outburst is due. The outburst should begin on 15 December 2015 if the spin value is in the middle of this range, on 3 January 2016 if the spin is 0.25, and on 26 November 2015 if the spin is 0.31. We have also tested the possibility that the quadrupole term in the post-Newtonian equations of motion does not exactly follow Einstein’s theory: a parameter q is introduced as one of the 8 parameters. Its value is within 30% (1 sigma) of the Einstein’s value q = 1. This supports the no-hair theorem of black holes within the achievable precision. We have also measured the loss of orbital energy due to gravitational waves. The loss rate is found to agree with Einstein’s value with the accuracy of 2% (1 sigma). There is a possibility of improving the accuracy of both quantities using the exact timing of the outburst on 21 July 2019. Because of closeness of OJ287 to the Sun (8–12°), the observations would be best carried out by a telescope in space.

Third-and-a-half order post-Newtonian equations of motion for relativistic compact binaries using the strong field point particle limit

We report our rederivation of the equations of motion for relativistic compact binaries through the third-and-a-half post-Newtonian (3.5 PN) order approximation to general relativity using the strong field point particle limit to describe self-gravitating stars instead of the Dirac delta functional. The computation is done in harmonic coordinates. Our equations of motion describe the orbital motion of the binary consisting of spherically symmetric non-rotating stars. The resulting equations of motion fully agree with the 3.5 PN equations of motion derived in the previous works. We also show that the locally defined energy of the star has a simple relation with its mass up to the 3.5 PN order.

Effective field theory calculation of second post-Newtonian binary dynamics

We use the effective field theory for gravitational bound states, proposed by Goldberger and Rothstein, to compute the interaction Lagrangian of a binary system at the second post-Newtonian order. Throughout the calculation, we use a metric parametrization based on a temporal Kaluza-Klein decomposition and test the claim by Kol and Smolkin that this parametrization provides important calculational advantages. We demonstrate how to use the effective field theory method efficiently in precision calculations, and we reproduce known results for the second post-Newtonian order equations of motion in harmonic gauge in a straightforward manner.

The 2nd-order Post-Newtonian Orbit Equation of Light

Based on the 2nd-order post-Newtonian approximation under the DSX frame of the general relativity theory, the 2nd-order post-Newtonian orbital equation of light in the axis-symmetrical stationary spacetime is derived, and from this, the angle of deflection of light propagating in the equatorial plane is derived. The obtained results are consistent with those of the Schwarzchild and Kerr metrics within the limits of measuring precision.

Summary of sessions B1/B2 and B2: relativistic astrophysics and numerical relativity

The numerical relativity session at GR18 was dominated by physics results on binary black hole mergers. Several groups can now simulate these from a time when the post-Newtonian equations of motion are still applicable, through several orbits and the merger to the ringdown phase, obtaining plausible gravitational waves at infinity, and showing some evidence of convergence with resolution. The results of different groups roughly agree. This new-won confidence has been used by these groups to begin mapping out the (finite dimensional) initial data space of the problem, with a particular focus on the effect of black hole spins, and the acceleration by gravitational wave recoil to hundreds of km s−1 of the final merged black hole. Other work was presented on a variety of topics, such as evolutions with matter, extreme mass ratio inspirals and technical issues such as gauge choices.

Reducing eccentricity in black-hole binary evolutions with initial parameters from post-Newtonian inspiral

Standard choices of quasicircular orbit parameters for black-hole binary evolutions result in eccentric inspiral. We introduce a conceptually simple method, which is to integrate the post-Newtonian equations of motion through hundreds of orbits, and read off the values of the momenta at the separation at which we wish to start a fully general relativistic numerical evolution. For the particular case of nonspinning equal-mass inspiral with an initial coordinate separation of $D=11M$ we show that this approach reduces the eccentricity by at least a factor of 5 to $e&lt;0.002$ as compared to using standard quasicircular initial parameters.

Post-Newtonian diagnosis of quasiequilibrium configurations of neutron star–neutron star and neutron star–black hole binaries

We use a post-Newtonian diagnostic tool to examine numerically generated quasiequilibrium initial data sets for nonspinning double neutron star and neutron star-black hole binary systems. The post-Newtonian equations include the effects of tidal interactions, parametrized by the compactness of the neutron stars and by suitable values of 'apsidal' constants, which measure the degree of distortion of stars subjected to tidal forces. We find that the post-Newtonian diagnostic agrees well with the double neutron star initial data, typically to better than half a percent except where tidal distortions are becoming extreme. We show that the differences could be interpreted as representing small residual eccentricity in the initial orbits. In comparing the diagnostic with preliminary numerical data on neutron star-black hole binaries, we find less agreement.

On derivation of EIH (Einstein–Infeld–Hoffman) equations of motion from the linearized metric of general relativity theory

The half-century old idea of Infeld to use the variational principle of the general relativity field equations is reminded to show that the commonly employed EIH (Einstein–Infeld–Hoffman) equations of motion may be derived from the linearized (weak-field) metric alone. Based on it, the linearized metric might be sufficient for post-Newtonian celestial mechanics and astrometry enabling one to derive the post-Newtonian equations of motion and rotation of celestial bodies as well as the post-Newtonian equations of light propagation within the general relativity framework.

Strong field effects on binary systems in Einstein-aether theory

``Einstein--aether'' theory is a generally covariant theory of gravity containing a dynamical preferred frame. This article continues an examination of effects on the motion of binary pulsar systems in this theory, by incorporating effects due to strong fields in the vicinity of neutron star pulsars. These effects are included through an effective approach, by treating the compact bodies as point particles with nonstandard, velocity dependent interactions parametrized by dimensionless sensitivities. Effective post-Newtonian equations of motion for the bodies and the radiation damping rate are determined. More work is needed to calculate values of the sensitivities for a given fluid source; therefore, precise constraints on the theory's coupling constants cannot yet be stated. It is shown, however, that strong field effects will be negligible given current observational uncertainties if the dimensionless couplings are less than roughly 0.1 and two conditions that match the PPN parameters to those of pure general relativity are imposed. In this case, weak field results suffice. There then exists a one-parameter family of Einstein--aether theories with ``small-enough'' couplings that passes all current observational tests. No conclusion can be reached for larger couplings until the sensitivities for a given source can be calculated.

Inspiral, merger, and ringdown of unequal mass black hole binaries: A multipolar analysis

We study the inspiral, merger, and ringdown of unequal mass black hole binaries by analyzing a catalogue of numerical simulations for seven different values of the mass ratio (from q=M2/M1=1 to q=4). We compare numerical and post-Newtonian results by projecting the waveforms onto spin-weighted spherical harmonics, characterized by angular indices (l,m). We find that the post-Newtonian equations predict remarkably well the relation between the wave amplitude and the orbital frequency for each (l,m), and that the convergence of the post-Newtonian series to the numerical results is nonmonotonic. To leading order, the total energy emitted in the merger phase scales like η2 and the spin of the final black hole scales like η, where η=q/(1+q)2 is the symmetric mass ratio. We study the multipolar distribution of the radiation, finding that odd-l multipoles are suppressed in the equal mass limit. Higher multipoles carry a larger fraction of the total energy as q increases. We introduce and compare three different definitions for the ringdown starting time. Applying linear-estimation methods (the so-called Prony methods) to the ringdown phase, we find resolution-dependent time variations in the fitted parameters of the final black hole. By cross correlating information from different multipoles, we show that ringdown fits can be used to obtain precise estimates of the mass and spin of the final black hole, which are in remarkable agreement with energy and angular momentum balance calculations.