Analytic modelling of tidal effects in the relativistic inspiral of binary neutron stars

To detect the gravitational-wave (GW) signal from binary neutron stars and extract information about the equation of state of matter at nuclear density, it is necessary to match the signal with a bank of accurate templates. We present the two longest (to date) general-relativistic simulations of equal-mass binary neutron stars with different compactnesses, C=0.12 and C=0.14, and compare them with a tidal extension of the effective-one-body (EOB)model. The typical numerical phasing errors over the $\simeq 22$ GW cycles are $\Delta \phi\simeq \pm 0.24$ rad. By calibrating only one parameter (representing a higher-order amplification of tidal effects), the EOB model can reproduce, within the numerical error, the two numerical waveforms essentially up to the merger. By contrast, the third post-Newtonian Taylor-T4 approximant with leading-order tidal corrections dephases with respect to the numerical waveforms by several radians

Relativistic tidal properties of neutron stars

We study the various linear responses of neutron stars to external relativistic tidal fields. We focus on three different tidal responses, associated to three different tidal coefficients: (i) a gravito-electric-type coefficient G\mu_\ell=[length]^{2\ell+1} measuring the \ell^{th}-order mass multipolar moment GM_{a_1... a_\ell} induced in a star by an external \ell^{th}-order gravito-electric tidal field G_{a_1... a_\ell}; (ii) a gravito-magnetic-type coefficient G\sigma_\ell=[length]^{2\ell+1} measuring the \ell^{th} spin multipole moment G S_{a_1... a_\ell} induced in a star by an external \ell^{th}-order gravito-magnetic tidal field H_{a_1... a_\ell}; and (iii) a dimensionless ``shape'' Love number h_\ell measuring the distortion of the shape of the surface of a star by an external \ell^{th}-order gravito-electric tidal field. All the dimensionless tidal coefficients G\mu_\ell/R^{2\ell+1}, G\sigma_\l/R^{2\ell+1} and h_\ell (where R is the radius of the star) are found to have a strong sensitivity to the value of the star's ``compactness'' c\equiv GM/(c_0^2 R) (where we indicate by c_0 the speed of light). In particular, G\mu_\l/R^{2\l+1}\sim k_\ell is found to strongly decrease, as c increases, down to a zero value as c is formally extended to the ``black-hole (BH) limit'' c^{BH}=1/2. The shape Love number h_\ell is also found to significantly decrease as c increases, though it does not vanish in the formal limit c\to c^{BH}. The formal vanishing of \mu_\ell and \sigma_\ell as c\to c^{BH} is a consequence of the no-hair properties of black holes; this suggests, but in no way proves, that the effective action describing the gravitational interactions of black holes may not need to be augmented by nonminimal worldline couplings.Comment: 21 pages, 10 figures. Matches the published versio

Binary black hole merger in the extreme mass ratio limit

We discuss the transition from quasi-circular inspiral to plunge of a system of two nonrotating black holes of masses $m_1$ and $m_2$ in the extreme mass ratio limit $m_1m_2\ll (m_1+m_2)^2$. In the spirit of the Effective One Body (EOB) approach to the general relativistic dynamics of binary systems, the dynamics of the two black hole system is represented in terms of an effective particle of mass $\mu\equiv m_1m_2/(m_1+m_2)$ moving in a (quasi-)Schwarzschild background of mass $M\equiv m_1+m_2$ and submitted to an ${\cal O}(\mu)$ radiation reaction force defined by Pad\'e resumming high-order Post-Newtonian results. We then complete this approach by numerically computing, \`a la Regge-Wheeler-Zerilli, the gravitational radiation emitted by such a particle. Several tests of the numerical procedure are presented. We focus on gravitational waveforms and the related energy and angular momentum losses. We view this work as a contribution to the matching between analytical and numerical methods within an EOB-type framework.Comment: 14 pages, six figures. Revised version. To appear in the CQG special issue based around New Frontiers in Numerical Relativity conference, Golm (Germany), July 17-21 200

Gradient Representations and Affine Structures in AE(n)

We study the indefinite Kac-Moody algebras AE(n), arising in the reduction of Einstein's theory from (n+1) space-time dimensions to one (time) dimension, and their distinguished maximal regular subalgebras sl(n) and affine A_{n-2}^{(1)}. The interplay between these two subalgebras is used, for n=3, to determine the commutation relations of the `gradient generators' within AE(3). The low level truncation of the geodesic sigma-model over the coset space AE(n)/K(AE(n)) is shown to map to a suitably truncated version of the SL(n)/SO(n) non-linear sigma-model resulting from the reduction Einstein's equations in (n+1) dimensions to (1+1) dimensions. A further truncation to diagonal solutions can be exploited to define a one-to-one correspondence between such solutions, and null geodesic trajectories on the infinite-dimensional coset space H/K(H), where H is the (extended) Heisenberg group, and K(H) its maximal compact subgroup. We clarify the relation between H and the corresponding subgroup of the Geroch group.Comment: 43 page

Reduced Hamiltonian for next-to-leading order Spin-Squared Dynamics of General Compact Binaries

Within the post Newtonian framework the fully reduced Hamiltonian (i.e., with eliminated spin supplementary condition) for the next-to-leading order spin-squared dynamics of general compact binaries is presented. The Hamiltonian is applicable to the spin dynamics of all kinds of binaries with self-gravitating components like black holes and/or neutron stars taking into account spin-induced quadrupolar deformation effects in second post-Newtonian order perturbation theory of Einstein's field equations. The corresponding equations of motion for spin, position and momentum variables are given in terms of canonical Poisson brackets. Comparison with a nonreduced potential calculated within the Effective Field Theory approach is made.Comment: 11 pages, minor changes to match published version at CQ

Covariant perturbations of f(R) black holes: the Weyl terms

In this paper we revisit non-spherical perturbations of the Schwarzschild black hole in the context of f(R) gravity. Previous studies were able to demonstrate the stability of the f(R) Schwarzschild black hole against gravitational perturbations in both the even and odd parity sectors. In particular, it was seen that the Regge-Wheeler and Zerilli equations in f(R) gravity obey the same equations as their General Relativity counterparts. More recently, the 1+1+2 semi-tetrad formalism has been used to derive a set of two wave equations: one for transverse, trace-free (tensor) perturbations and one for the additional scalar modes that characterise fourth-order theories of gravitation. The master variable governing tensor perturbations was shown to be a modified Regge-Wheeler tensor obeying the same equation as in General Relativity. However, it is well known that there is a non-uniqueness in the definition of the master variable. In this paper we derive a set of two perturbation variables and their concomitant wave equations that describe gravitational perturbations in a covariant and gauge invariant manner. These variables can be related to the Newman-Penrose (NP) Weyl scalars as well as the master variables from the 2+2 formalism

On the Detection of a Scalar Stochastic Background of Gravitational Waves

In the near future we will witness the coming to a full operational regime of laser interferometers and resonant mass detectors of spherical shape. In this work we study the sensitivity of pairs of such gravitational wave detectors to a scalar stochastic background of gravitational waves. Our computations are carried out both for minimal and non minimal coupling of the scalar fields.Comment: 25 pages, 3 figure

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.Comment: 51 pages, 10 figures; published versio

The Current Status of Binary Black Hole Simulations in Numerical Relativity

Since the breakthroughs in 2005 which have led to long term stable solutions of the binary black hole problem in numerical relativity, much progress has been made. I present here a short summary of the state of the field, including the capabilities of numerical relativity codes, recent physical results obtained from simulations, and improvements to the methods used to evolve and analyse binary black hole spacetimes.Comment: 14 pages; minor changes and corrections in response to referee

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