Source author record

Nicolas Yunes

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

89works
12topics
4close collaborators

Actions

Connect this record

Log in to claim

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

89 published item(s)

preprint2022arXiv

Are Parametrized Tests of General Relativity with Gravitational Waves Robust to Unknown Higher Post-Newtonian Order Effects?

Gravitational wave observations have great potential to reveal new information about the fundamental nature of gravity, but extracting that information can be difficult. One popular technique is the parametrized inspiral test of general relativity (a realization of the parametrized post-Einsteinian framework), where the gravitational waveform, as calculated in Einstein's theory as a series expansion in the orbital velocity, is parametrically deformed at a given set of orders in velocity. However, most current approaches usually only analyze the data while considering a single, specific modification at a time. Are then constraints placed with a single modification robust to our ignorance of higher post-Newtonian order corrections? We show here that for a wide class of theories, specifically those that admit a post-Newtonian expansion, single-parameter tests are indeed robust. In particular, through a series of full Bayesian parameter estimation studies on several different sets of synthetic data, we show that single-parameter constraints are not degraded but rather are improved by the inclusion of multiple parameters, provided one includes information about the mathematical structure of the series. We then exemplify this with a specific theory of gravity, shift-symmetric scalar Gauss-Bonnet theory, where the waveform has been calculated to higher post-Newtonian orders than leading. We show that the inclusion of these higher order terms strengthens single-parameter constraints, instead of weakening them, and that the strengthening is very mild. This analysis therefore provides strong evidence that single-parameter post-Einsteinian tests of general relativity are robust to ignorance of high post-Newtonian order terms in the general relativistic deformations.

preprint2022arXiv

Chaos in Quadratic Gravity

While recent gravitational wave observations by LIGO and Virgo allow for tests of general relativity in the extreme gravity regime, these observations are still blind to a large swath of phenomena outside these instruments' sensitivity curves. Future gravitational-wave detectors, such as LISA, will enable probes of longer-duration and lower-frequency events. In particular, LISA will enable the characterization of the non-linear dynamics of extreme mass-ratio inspirals, when a small compact object falls into a supermassive black hole. In this paper, we study the motion of test particles around spinning black holes in two quadratic gravity theories: scalar Gauss-Bonnet and dynamical Chern-Simons gravity. We show that geodesic trajectories around rotating black holes in these theories are likely to not have a fourth constant of the motion. In particular, we show that Poincaré sections of the orbital phase space present chaotic features that will affect the inspiral of small compact objects into supermassive black holes in these theories. Nevertheless, the characteristic size of these chaotic features is tiny and their location in parameter space is very close to the event horizon of the supermassive black hole. Therefore, the detection of such chaotic features with LISA is likely very challenging, at best.

preprint2022arXiv

Gravitational-Wave and X-ray Probes of the Neutron Star Equation of State

Neutron stars are a remarkable marriage of Einstein's theory of general relativity with nuclear physics. Their interiors harbor extreme matter that cannot be probed in the laboratory. At such high densities and pressures, their cores may consist predominantly of exotic matter such as free quarks or hyperons. Gravitational wave observations from the Laser Interferometer Gravitational-wave Observatory (LIGO) and from other interferometers, and X-ray observations from the Neutron Star Interior Composition Explorer (NICER), are beginning to pierce through the veil. These observations provide information about neutron star cores, and therefore, about the physics that makes such objects possible. In this review, we discuss what we have learned about the physics of neutron stars from gravitational wave and X-ray observations. We focus on what has been observed with certainty and what should be observable in the near future, with an eye out for the physics that these new observations will teach us.

preprint2022arXiv

How Do Spherical Black Holes Grow Monopole Hair?

Black holes in certain modified gravity theories that contain a scalar field coupled to curvature invariants are known to possess (monopole) scalar hair while non-black-hole spacetimes (like neutron stars) do not. Therefore, as a neutron star collapses to a black hole, scalar hair must grow until it settles to the stationary black hole solution with (monopole) hair. In this paper, we study this process in detail and show that the growth of scalar hair is tied to the appearance and growth of the event horizon (before an apparent horizon forms), which forces scalar modes that would otherwise (in the future) become divergent to be radiated away. We prove this result rigorously in general first for a large class of modified theories, and then we exemplify the results by studying the temporal evolution of the scalar field in scalar Gauss-Bonnet gravity in two backgrounds: (i) a collapsing Oppenheimer-Snyder background, and (ii) a collapsing neutron star background. In case (i), we find an exact scalar field solution analytically, while in case (ii) we solve for the temporal evolution of the scalar field numerically, with both cases supporting the conclusion presented above. Our results suggest that the emission of a burst of scalar field radiation is a necessary condition for black hole formation in a large class of modified theories of gravity.

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.

preprint2022arXiv

Parity Violation in Spin-Precessing Binaries: Gravitational Waves from the Inspiral of Black Holes in Dynamical Chern-Simons Gravity

Spin precession in compact binaries is intricately tuned to the multipole structure of the underlying bodies. For black holes, violations of the no-hair theorems induced by modifications to general relativity correct the precession dynamics, which in turn imprints onto the amplitude and phase modulations of the gravitational waves emitted by the binary. Recently, the spin precession equations were derived up to second order in spin for dynamical Chern-Simons gravity, a parity violating modified theory of gravity. We here solve these equations and construct, for the first time, analytic expressions for the time- and frequency-domain gravitational waves emitted in the quasi-circular inspiral of spin-precessing black hole binaries in a modified theory of gravity using the post-Newtonian approximation. Working within the small coupling approximation and using multiple scale analysis, we show that the corrections to the nutation phase enter at relative 1PN order, and the corrections to the precession angle and Thomas phase enter at relative 0PN order. Making use of the stationary phase approximation and shifted uniform asymptotics, we find that the Fourier phase of the waveform is characterized by three modifications, two due to the back-reaction of the precession dynamics onto the spin-orbit and spin-spin couplings that enter at 1.5PN and 2PN orders, and a 2PN modification due to the emission of dipole radiation. We also find that back-reaction of the precession dynamics forces the dCS corrections to the Fourier amplitude to enter at 0PN order, as opposed to 2PN order, as expected for spin-aligned binaries. Our work lays the first foundational stones to build an inspiral-merger-ringdown phenomenological model for spin-precessing binaries in a modified theory of gravity.

preprint2022arXiv

Quasinormal modes of slowly-rotating black holes in dynamical Chern-Simons gravity

The detection of gravitational waves from compact binary mergers by the LIGO/Virgo collaboration has, for the first time, allowed us to test relativistic gravity in its strong, dynamical and nonlinear regime, thus opening a new arena to confront general relativity (and modifications thereof) against observations. We consider a theory which modifies general relativity by introducing a scalar field coupled to a parity-violating curvature term known as dynamical Chern-Simons gravity. In this theory, spinning black holes are different from their general relativistic counterparts and can thus serve as probes to this theory. We study linear gravito-scalar perturbations of black holes in dynamical Chern-Simons gravity at leading-order in spin and (i) obtain the perturbed field equations describing the evolution of the perturbed gravitational and scalar fields, (ii) numerically solve these equations by direct integration to calculate the quasinormal mode frequencies for the dominant and higher multipoles and tabulate them, (iii) find strong evidence that these rotating black holes are linearly stable, and (iv) present general fitting functions for different multipoles for gravitational and scalar quasinormal mode frequencies in terms of spin and Chern-Simons coupling parameter. Our results can be used to validate the ringdown of small-spin remnants of numerical relativity simulations of black hole binaries in dynamical Chern-Simons gravity and pave the way towards future tests of this theory with gravitational wave ringdown observations.

preprint2021arXiv

The slope, the hill, the drop, and the swoosh: Learning about the nuclear matter equation of state from the binary Love relations

Analyses that connect astrophysical observations of neutron stars with nuclear matter properties sometimes rely on equation-of-state insensitive relations. We show that the slope of the binary Love relations (i.e.~between the tidal deformabilities of binary neutron stars) encodes the rate of change of the nuclear matter speed of sound below three times nuclear saturation density. Twin stars lead to relations that present a signature ''hill'', ''drop'', and ''swoosh'' due to the second (mass-radius) stable branch, requiring a new description of the binary love relations.

preprint2020arXiv

Gravitational Waves from Quasi-Circular Black Hole Binaries in Dynamical Chern-Simons Gravity

Dynamical Chern-Simons gravity cannot be strongly constrained with current experiments because it reduces to General Relativity in the weak-field limit. This theory, however, introduces modifications in the non-linear, dynamical regime, and thus, it could be greatly constrained with gravitational waves from the late inspiral of black hole binaries. We complete the first self-consistent calculation of such gravitational waves in this theory. For favorable spin-orientations, advanced ground-based detectors may improve existing solar-system constraints by 6 orders of magnitude.

preprint2020arXiv

Gravitational-wave versus X-ray tests of strong-field gravity

Electromagnetic observations of the radiation emitted by an accretion disk around a black hole, as well as gravitational wave observations of coalescing binaries, can be used to probe strong-field gravity. We here compare the constraints that these two types of observations can impose on theory-agnostic, parametric deviations from the Schwarzschild metric. On the gravitational wave side, we begin by computing the leading-order deviation to the Hamiltonian of a binary system in a quasi-circular orbit within the post-Newtonian approximation, given a parametric deformation of the Schwarzschild metrics of each binary component. We then compute the leading-order deviation to the gravitational waves emitted by such a binary in the frequency domain, assuming purely Einsteinian radiation-reaction. We compare this model to the LIGO-Virgo collaboration gravitational wave detections and place constraints on the metric deformation parameters, concluding with an estimate of the constraining power of aLIGO at design sensitivity. On the electromagnetic side, we first simulate observations with current and future X-ray instruments of an X-ray binary with a parametrically-deformed Schwarzschild black hole, and we then estimate constraints on the deformation parameters using these observations.We find that current gravitational wave observations have already placed constraints on the metric deformation parameters than are slightly more stringent than what can be achieved with current X-ray instruments. Moreover, future gravitational wave observations with aLIGO at design sensitivity by 2026 will be even more stringent, becoming stronger than constraints achievable with future ATHENA X-ray observations before it flies in 2034.

preprint2020arXiv

Prospects for Fundamental Physics with LISA

In this paper, which is of programmatic rather than quantitative nature, we aim to further delineate and sharpen the future potential of the LISA mission in the area of fundamental physics. Given the very broad range of topics that might be relevant to LISA, we present here a sample of what we view as particularly promising directions, based in part on the current research interests of the LISA scientific community in the area of fundamental physics. We organize these directions through a "science-first" approach that allows us to classify how LISA data can inform theoretical physics in a variety of areas. For each of these theoretical physics classes, we identify the sources that are currently expected to provide the principal contribution to our knowledge, and the areas that need further development. The classification presented here should not be thought of as cast in stone, but rather as a fluid framework that is amenable to change with the flow of new insights in theoretical physics.

preprint2019arXiv

Gravitational Instability of Exotic Compact Objects

Exotic compact objects with physical surfaces a Planckian distance away from where the horizon would have been are inspired in quantum gravity. Most of these objects are defined by a classical spacetime metric, such as boson stars, gravastars and wormholes. We show that these classical objects are gravitationally unstable because accretion of ordinary and dark matter, and gravitational waves forces them to collapse to a black hole by the Hoop conjecture. To avoid collapse, their surface must be a macroscopic distance away from the horizon or they must violate the null energy condition.

preprint2016arXiv

Binary Love Relations

When in a tight binary, the mutual tidal deformations of neutron stars imprint onto observables, encoding information about their internal structure at supranuclear densities and gravity in the extreme-gravity regime. Gravitational wave observations of their late binary inspiral may serve as a tool to extract the individual tidal deformabilities, but this is made difficult by degeneracies between them in the gravitational wave model. We here resolve this problem by discovering approximately universal relations between dimensionless combinations of the individual tidal deformabilities. We show that these relations break degeneracies in the gravitational wave model, allowing for the accurate extraction of both deformabilities. Such measurements can be used to better differentiate between equation-of-state models, and improve tests of General Relativity and cosmology.

preprint2016arXiv

Can the Slow-Rotation Approximation be used in Electromagnetic Observations of Black Holes?

Future electromagnetic observations of black holes may allow us to test General Relativity in the strong-field regime. Such tests, however, require knowledge of rotating black hole solutions in modified gravity theories, a class of which does not admit the Kerr metric as a solution. Several rotating black hole solutions in modified theories have only been found in the slow-rotation approximation (i.e. assuming the spin angular momentum is much smaller than the mass squared). We here investigate whether the systematic error due to the approximate nature of these black hole metrics is small enough relative to the observational error to allow their use in electromagnetic observations to constrain deviations from General Relativity. We address this by considering whether electromagnetic observables constructed from a slow-rotation approximation to the Kerr metric can fit observables constructed from the full Kerr metric with systematic errors smaller than current observational errors. We focus on black hole shadow and continuum spectrum observations, as these are the least influenced by accretion disk physics, with current observational errors of about 10%. We find that the fractional systematic error introduced by using a second-order, slowly-rotating Kerr metric is at most 2% for shadows created by black holes with dimensionless spins $χ\leq0.6$. We also find that the systematic error introduced by using the slowly-rotating Kerr metric as an exact metric when constructing continuum spectrum observables is negligible for black holes with dimensionless spins of $χ\lesssim 0.9$. Our results suggest that the modified gravity solutions found in the slow-rotation approximation may be used to constrain realistic deviations from General Relativity with continuum spectrum and black hole shadow observations.

preprint2016arXiv

I-Love-Q Relations: From Compact Stars to Black Holes

The relations between most observables associated with a compact star, such as the mass and radius of a neutron star or a quark star, typically depend strongly on their unknown internal structure. The I-Love-Q relations (between the moment of inertia, the tidal deformability and the quadrupole moment) are however approximately insensitive to this structure. These relations become exact for stationary black holes in General Relativity as shown by the no-hair theorems. In this paper, we take the first steps toward studying how the approximate I-Love-Q relations become exact in the limit as compact stars become black holes. To do so, we consider a toy model, i.e. incompressible stars with anisotropic pressure, which allows us to model an equilibrium sequence of stars with their compactness approaching the black hole limit arbitrarily closely. We extract the I-Love-Q trio by numerically constructing such a sequence in the slow-rotation and small-tide approximations. We find that the I-Love-Q relations approach the black hole limit in a nontrivial way, with the quadrupole moment and the tidal deformability changing sign as the compactness and the amount of anisotropy are increased. Generalizing Maclaurin spheroids to anisotropic stars, we show that the multipole moments also change sign in the Newtonian limit as the amount of anisotropy is increased. We also prove analytically that the stellar moment of inertia reaches the black hole limit as the compactness reaches the black hole value in the strongly anisotropic limit. Modeling the black hole limit through a sequence of anisotropic stars, however, fails when considering other theories of gravity. We calculate the scalar dipole charge and the moment of inertia in a parity-violating modified theory and find that these quantities do not tend to their black hole counterparts as the anisotropic stellar sequence approaches the black hole limit.

preprint2016arXiv

Isolated and Binary Neutron Stars in Dynamical Chern-Simons Gravity

We study isolated and binary neutron stars in dynamical Chern-Simons gravity. This theory modifies the Einstein-Hilbert action through the introduction of a dynamical scalar field coupled to the Pontryagin density. We here treat this theory as an effective model, working to leading order in the Chern-Simons coupling. We first construct isolated neutron star solutions in the slow-rotation expansion to quadratic order in spin. We find that isolated neutron stars acquire a scalar dipole charge that corrects its spin angular momentum to linear order in spin and corrects its mass and quadrupole moment to quadratic order in spin, as measured by an observer at spatial infinity. We then consider neutron stars binaries that are widely separated and solve for their orbital evolution in this modified theory. We find that the evolution of post-Keplerian parameters is modified, with the rate of periastron advance being the dominant correction at first post-Newtonian order. We conclude by applying these results to observed pulsars with the aim to place constraints on dynamical Chern-Simons gravity. We find that the modifications to the observed mass are degenerate with the neutron star equation of state, which prevents us from testing the theory with the inferred mass of the millisecond pulsar J1614-2230. We also find that the corrections to the post-Keplerian parameters are too small to be observable today even with data from the double binary pulsar J0737-3039. Our results suggest that pulsar observations are not currently capable of constraining dynamical Chern-Simons gravity, and thus, gravitational-wave observations may be the only path to a stringent constraint of this theory in the imminent future.

preprint2016arXiv

Musings on Lorentz Violation Given the Recent Gravitational-Wave Observations of Coalescing Binary Black Holes

The recent observation of gravitational waves by the LIGO/Virgo collaboration provides a unique opportunity to probe the extreme gravity of coalescing binary black holes. In this regime, the gravitational interaction is not only strong, but the spacetime curvature is large, characteristic velocities are a non-negligible fraction of the speed of light, and the time scale on which the curvature and gravity change is small. This contribution discusses some consequences of these observations on modifications to General Relativity, with a special emphasis on Lorentz-violating theories.

preprint2016arXiv

The Effect of Cosmological Evolution on Solar System Constraints and on the Scalarization of Neutron Stars in Massless Scalar-Tensor Theories

Certain scalar-tensor theories of gravity that generalize Jordan-Fierz-Brans-Dicke theory are known to predict non-trivial phenomenology for neutron stars. In these theories, first proposed by Damour and Esposito-Farèse, the scalar field has a standard kinetic term, and couples conformally to the matter fields. The weak equivalence principle is therefore satisfied, but scalar effects may arise in strong-field regimes, e.g. allowing for violations of the strong equivalence principle in neutron stars ("spontaneous scalarization") or in sufficiently tight binary neutron-star systems ("dynamical/induced scalarization"). The original scalar-tensor theory proposed by Damour and Esposito-Farèse is in tension with solar-system constraints (for couplings that lead to scalarization), if one accounts for cosmological evolution of the scalar field and no mass term is included in the action. We here extend the conformal coupling of that theory, in order to ascertain if, in this way, solar-system tests can be passed, while retaining a non-trivial phenomenology for neutron stars. We find that even with this generalized conformal coupling, it is impossible to construct a theory that passes both Big-Bang nucleosynthesis and solar-system constraints, while simultaneously allowing for scalarization in isolated/binary neutron stars.

preprint2016arXiv

Theoretical Physics Implications of the Binary Black-Hole Mergers GW150914 and GW151226

The gravitational wave observations GW150914 and GW151226 by Advanced LIGO provide the first opportunity to learn about physics in the extreme gravity environment of coalescing binary black holes. The LIGO Scientific Collaboration and the Virgo Collaboration have verified that this observation is consistent with General Relativity. This paper expands their analysis to a larger class of anomalies, highlighting the inferences that can be drawn on non-standard theoretical physics mechanisms. We find that these events constrain a plethora of mechanisms associated with the generation and propagation of gravitational waves, including the activation of scalar fields, gravitational leakage into large extra dimensions, the variability of Newton's constant, a modified dispersion relation, gravitational Lorentz violation and the strong equivalence principle. Though other observations limit many of these mechanisms already, GW150914 and GW151226 are unique in that they are direct probes of dynamical strong-field gravity and of gravitational wave propagation. We also show that GW150914 constrains inferred properties of exotic compact object alternatives to Kerr black holes. We argue, however, that the true potential for GW150914 to both rule out exotic objects and constrain physics beyond General Relativity is severely limited by the lack of understanding of the merger regime in almost all relevant modified gravity theories. This event thus significantly raises the bar that these theories have to pass, both in terms of having a sound theoretical underpinning, and being able to solve the equations of motion for binary merger events. We conclude with a discussion of the additional inferences that can be drawn if the lower-confidence observation of an electromagnetic counterpart to GW150914 holds true; this would provide dramatic constraints on the speed of gravity and gravitational Lorentz violation.

preprint2016arXiv

Theory-Agnostic Constraints on Black-Hole Dipole Radiation with Multi-Band Gravitational-Wave Astrophysics

The aLIGO detection of the black-hole binary GW150914 opened a new era for probing extreme gravity. Many gravity theories predict the emission of dipole gravitational radiation by binaries. This is excluded to high accuracy in binary pulsars, but entire classes of theories predict this effect predominantly (or only) in binaries involving black holes. Joint observations of GW150914-like systems by aLIGO and eLISA will improve bounds on dipole emission from black-hole binaries by six orders of magnitude relative to current constraints, provided that eLISA is not dramatically descoped.

preprint2015arXiv

Challenging the Presence of Scalar Charge and Dipolar Radiation in Binary Pulsars

Corrections to general relativity that introduce long-ranged scalar fields which are non-minimally coupled to curvature typically predict that neutron stars possess a non-trivial scalar field profile. An observer far from a star is most sensitive to the spherically-symmetric piece of this profile that decays linearly with the inverse of the distance, the so-called scalar charge, which is related to the emission of dipolar radiation from compact binaries. The presence of dipolar radiation has the potential to very strongly constrain extended theories of gravity. These facts may lead people to believe that gravitational theories with long-ranged scalar fields have already been constrained strongly from binary pulsar observations. Here we challenge this "lore" by investigating the decoupling limit of Gauss-Bonnet gravity as an example, in which the scalar field couples linearly to the Gauss-Bonnet density in the action. We prove a theorem that neutron stars in this theory cannot possess a scalar charge. Thus Gauss-Bonnet gravity evades the strong binary pulsar constraints on dipole radiation. We discuss the astrophysical systems yielding the best constraints on Gauss-Bonnet gravity and related quadratic gravity theories. To achieve this we explicitly compute the scalar charge for slowly-rotating neutron stars in quadratic gravity theories. In generic case, either neutron star-binary or neutron star-black hole systems can be used to constrain the theory, but Gauss-Bonnet gravity evades the neutron star-binary constraints. However, black holes in Gauss-Bonnet gravity do anchor scalar charge. The best constraints on Gauss-Bonnet gravity will thus come from black hole observations, for example via pulsar-black hole binaries. We estimate these constraints to be ten times better than the current estimated bound, and also include estimated constraints on generic quadratic gravity theories.

preprint2015arXiv

Four-Hair Relations for Differentially Rotating Neutron Stars in the Weak-Field Limit

The opportunity to study physics at supra-nuclear densities through X-ray observations of neutron stars has led to in-depth investigations of certain approximately universal relations that can remove degeneracies in pulse profile models. One such set of relations determines all of the multipole moments of a neutron star just from the first three (the mass monopole, the current dipole and the mass quadrupole moment) approximately independently of the equation of state. These three-hair relations were found to hold in neutron stars that rotate rigidly, as is the case in old pulsars, but neutron stars can also rotate differentially, as is the case for proto-neutron stars and hypermassive transient remnants of binary mergers. We here extend the three-hair relations to differentially rotating stars for the first time with a generic rotation law using two approximations: a weak-field scheme (an expansion in powers of the neutron star compactness) and a perturbative differential rotation scheme (an expansion about rigid rotation). These approximations allow us to analytically derive approximately universal relations that allow us to determine all of the multipole moments of a (perturbative) differentially rotating star in terms of only the first four moments. These new four-hair relations for differentially rotating neutron stars are found to be approximately independent of the equation of state to a higher degree than the three-hair relations for uniformly rotating stars. Our results can be instrumental in the development of four-hair relations for rapidly differentially rotating stars in full General Relativity using numerical simulations.

preprint2015arXiv

Improved Universality in the Neutron Star Three-Hair Relations

No-hair like relations between the multipole moments of the exterior gravitational field of neutron stars have recently been found to be approximately independent of the star's internal structure. This approximate, equation-of-state universality arises after one adimensionalizes the multipole moments appropriately, which then begs the question of whether there are better ways to adimensionalize the moments to obtain stronger universality. We here investigate this question in detail by considering slowly-rotating neutron stars both in the non-relativistic limit and in full General Relativity. We find that there exist normalizations that lead to stronger equation-of-state universality in the relations among the moment of inertia and the quadrupole, octopole and hexadecapole moments of neutron stars. We determine the optimal normalization that minimizes the equation-of-state dependence in these relations. The results found here may have applications in the modeling of X-ray pulses and atomic line profiles from millisecond pulsars with NICER and LOFT.

preprint2015arXiv

Post-Newtonian, Quasi-Circular Binary Inspirals in Quadratic Modified Gravity

We consider a general class of quantum gravity-inspired, modified gravity theories, where the Einstein-Hilbert action is extended through the addition of all terms quadratic in the curvature tensor coupled to scalar fields with standard kinetic energy. This class of theories includes Einstein-Dilaton-Gauss-Bonnet and Chern-Simons modified gravity as special cases. We analytically derive and solve the coupled field equations in the post-Newtonian approximation, assuming a comparable-mass, spinning black hole binary source in a quasi-circular, weak-field/slow-motion orbit. We find that a naive subtraction of divergent piece associated with the point-particle approximation is ill-suited to represent compact objects in these theories. Instead, we model them by appropriate effective sources built so that known strong-field solutions are reproduced in the far-field limit. In doing so, we prove that black holes in Einstein-Dilaton-Gauss-Bonnet and Chern-Simons theory can have hair, while neutron stars have no scalar monopole charge, in diametrical opposition to results in scalar-tensor theories. We then employ techniques similar to the direct integration of the relaxed Einstein equations to obtain analytic expressions for the scalar field, metric perturbation, and the associated gravitational wave luminosity measured at infinity. We find that scalar field emission mainly dominates the energy flux budget, sourcing electric-type (even-parity) dipole scalar radiation and magnetic-type (odd-parity) quadrupole scalar radiation, correcting the General Relativistic prediction at relative -1PN and 2PN orders. Such modifications lead to corrections in the emitted gravitational waves that can be mapped to the parameterized post-Einsteinian framework. Such modifications could be strongly constrained with gravitational wave observations.

preprint2015arXiv

Probing the Internal Composition of Neutron Stars with Gravitational Waves

Gravitational waves from neutron star binary inspirals contain information about the equation of state of supranuclear matter. In the absence of definitive experimental evidence that determines the correct equation of state, a number of diverse models that give the pressure in a neutron star as function of its density have been proposed. These models differ not only in the approximations and techniques they use to solve the many-body Schrödinger equation, but also in the neutron star composition they assume. We study whether gravitational wave observations of neutron star binaries in quasicircular inspirals will allow us to distinguish between equations of state of differing internal composition, thereby providing important information about the properties of extremely high density matter. We carry out a Bayesian model selection analysis, and find that second generation gravitational wave detectors can heavily constrain equations of state that contain only quark matter, but hybrid stars containing both normal and quark matter are harder to distinguish from normal matter stars. A gravitational wave detection with a signal-to-noise ratio of 30 and masses around $1.4M_{\odot}$ could either detect or rule out strange quark stars with a 20 to 1 confidence. The presence of kaon condensates or hyperons in neutron star inner cores cannot be easily confirmed. For example, for the equations of state studied in this paper, even a gravitational wave signal with a signal-to-noise ratio as high as 60 would not allow us to claim a detection of kaon condensates or hyperons with confidence greater than 5 to 1. On the other hand, if kaon condensates and hyperons do not form in neutron stars, a gravitational wave signal with similar signal-to-noise ratio would be able to constrain their existence with an 11 to 1 confidence for high-mass systems.

preprint2015arXiv

Relating Follicly-Challenged Compact Stars to Bald Black Holes

Compact stars satisfy certain no-hair relations through which their multipole moments are given by their mass, spin and quadrupole moment. These relations are approximately independent of their equation of state, relating pressure to density. Such relations are similar to the black hole no-hair theorems, but these possess event horizons inside which information that led to their formation can hide. Compact stars do not possess horizons, so whether their no-hair relations are related to the black hole ones is unclear. We investigate how the two relations are related by studying relations among multipole moments for compact stars with anisotropic pressure as a toy model, which allows such stars to be more compact than those with isotropic pressure. We here show numerically that the compact star no-hair relations approach the black hole ones as the compactness approaches that of a black hole. We also prove analytically that the current dipole moment exactly reaches the black hole limit quadratically in compactness as strongly-anisotropic stars approach the black hole limit. We moreover show that compact stars become progressively oblate in this limit, even if prolate at low compactness due to strong anisotropies.

preprint2015arXiv

Slowly-Rotating Black Holes in Einstein-Dilaton-Gauss-Bonnet Gravity: Quadratic Order in Spin Solutions

We derive a stationary and axisymmetric black hole solution in Einstein-Dilaton-Gauss-Bonnet gravity to quadratic order in the ratio of the spin angular momentum to the black hole mass squared. This solution introduces new corrections to previously found nonspinning and linear-in-spin solutions. The location of the event horizon and the ergosphere are modified, as well as the quadrupole moment. The new solution is of Petrov type I, although lower order in spin solutions are of Petrov type D. There are no closed timelike curves or spacetime regions that violate causality outside of the event horizon in the new solution. We calculate the modifications to the binding energy, Kepler's third law, and properties of the innermost stable circular orbit. These modifications are important for determining how the electromagnetic properties of accretion disks around supermassive black holes are changed from those expected in general relativity.

preprint2014arXiv

A Workshop that Works

The main goal of a scientific workshop is to bring together experts in a specific field or related fields to collaborate, to discuss, and to creatively make progress in a particular area. The organizational aspects of such a meeting play a critical role in achieving these goals. We here present suggestions from scientists to scientists that hopefully help in organizing a successful scientific workshop that maximizes collaboration and creativity.

preprint2014arXiv

Accurate and efficient waveforms for compact binaries on eccentric orbits

Compact binaries that emit gravitational waves in the sensitivity band of ground-based detectors can have non-negligible eccentricities just prior to merger, depending on the formation scenario. We develop a purely analytic, frequency-domain model for gravitational waves emitted by compact binaries on orbits with small eccentricity, which reduces to the quasi-circular post-Newtonian approximant TaylorF2 at zero eccentricity and to the post-circular approximation of Yunes et al. (2009) at small eccentricity. Our model uses a spectral approximation to the (post-Newtonian) Kepler problem to model the orbital phase as a function of frequency, accounting for eccentricity effects up to ${\cal{O}}(e^8)$ at each post-Newtonian order. Our approach accurately reproduces an alternative time-domain eccentric waveform model for eccentricities $e\in [0, 0.4]$ and binaries with total mass less than 12 solar masses. As an application, we evaluate the signal amplitude that eccentric binaries produce in different networks of existing and forthcoming gravitational waves detectors. Assuming a population of eccentric systems containing black holes and neutron stars that is uniformly distributed in co-moving volume, we estimate that second generation detectors like Advanced LIGO could detect approximately 0.1-10 events per year out to redshift $z\sim 0.2$, while an array of Einstein Telescope detectors could detect hundreds of events per year to redshift $z \sim 2.3$.

preprint2014arXiv

Applicability of the Newman-Janis Algorithm to Black Hole Solutions of Modified Gravity Theories

The Newman-Janis algorithm has been widely used to construct rotating black hole solutions from non-rotating counterparts. While this algorithm was developed within General Relativity, it has more recently been applied to non-rotating solutions in modified gravity theories. We find that the application of the Newman-Janis algorithm to an arbitrary non-GR spherically-symmetric solution introduces pathologies in the resulting axially-symmetric metric. This then establishes that, in general, the Newman-Janis algorithm should not used to construct rotating black hole solutions outside of General Relativity.

preprint2014arXiv

Constraining Lorentz-violating, Modified Dispersion Relations with Gravitational Waves

Modified gravity theories generically predict a violation of Lorentz invariance, which may lead to a modified dispersion relation for propagating modes of gravitational waves. We construct a parametrized dispersion relation that can reproduce a range of known Lorentz-violating predictions and investigate their impact on the propagation of gravitational waves. A modified dispersion relation forces different wavelengths of the gravitational wave train to travel at slightly different velocities, leading to a modified phase evolution observed at a gravitational-wave detector. We show how such corrections map to the waveform observable and to the parametrized post-Einsteinian framework, proposed to model a range of deviations from General Relativity. Given a gravitational-wave detection, the lack of evidence for such corrections could then be used to place a constraint on Lorentz violation. The constraints we obtain are tightest for dispersion relations that scale with small power of the graviton's momentum and deteriorate for a steeper scaling.

preprint2014arXiv

Detection and Parameter Estimation of Gravitational Waves from Compact Binary Inspirals with Analytical Double-Precessing Templates

We study the performance of various analytical frequency-domain templates for detection and parameter estimation of gravitational waves from spin-precessing, quasi-circular, compact binary inspirals. We begin by assessing the extent to which non-spinning, spin-aligned, and the new (analytical, frequency-domain, small-spin) double-precessing frequency-domain templates can be used to detect signals from such systems. For effective, dimensionless spin values above $0.2$, the use of non-spinning or spin-aligned templates for detection purposes will result in a loss of up to $30%$ of all events, while in the case of the double-precessing model, this never exceeds $6%$. Moreover, even for signals from systems with small spins, non-spinning and spin-aligned templates introduce large biases in the extracted masses and spins. The use of a model that encodes spin-induced precession effects, such as the double-precessing model, improves the mass and spin extraction by up to an order of magnitude. The additional information encoded in the spin-orbit interaction is invaluable if one wishes to extract the maximum amount of information from gravitational wave signals.

preprint2014arXiv

Effective No-Hair Relations for Neutron Stars and Quark Stars: Relativistic Results

Astrophysical charge-free black holes are known to satisfy no-hair relations through which all multipole moments can be specified in terms of just their mass and spin angular momentum. We here investigate the possible existence of no-hair-like relations among multipole moments for neutron stars and quark stars that are independent of their equation of state. We calculate the multipole moments of these stars up to hexadecapole order by constructing uniformly-rotating and unmagnetized stellar solutions to the Einstein equations. For slowly-rotating stars, we construct stellar solutions to quartic order in spin in a slow-rotation expansion, while for rapidly-rotating stars, we solve the Einstein equations numerically with the LORENE and RNS codes. We find that the multipole moments extracted from these numerical solutions are consistent with each other. We confirm that the current-dipole is related to the mass-quadrupole in an approximately equation of state independent fashion, which does not break for rapidly rotating neutron stars or quark stars. We further find that the current-octupole and the mass-hexadecapole moments are related to the mass-quadrupole in an approximately equation of state independent way to $\sim 10%$, worsening in the hexadecapole case. All of our findings are in good agreement with previous work that considered stellar solutions to leading-order in a weak-field expansion. The quartic in spin, slowly-rotating solutions found here allow us to estimate the systematic errors in the measurement of the neutron star's mass and radius with future X-ray observations, such as NICER and LOFT. We find that the effect of these quartic-in-spin terms on the quadrupole and hexadecapole moments and stellar eccentricity may dominate the error budget for very rapidly-rotating neutron stars. The new universal relations found here should help to reduce such systematic errors.

preprint2014arXiv

Linear Stability Analysis of Dynamical Quadratic Gravity

We perform a linear stability analysis of dynamical, quadratic gravity in the high-frequency, geometric optics approximation. This analysis is based on a study of gravitational and scalar modes propagating on spherically-symmetric and axially-symmetric, vacuum solutions of the theory. We find dispersion relations that do no lead to exponential growth of the propagating modes, suggesting the theory is linearly stable on these backgrounds. The modes are found to propagate at subluminal and superluminal speeds, depending on the propagating modes' direction relative to the background geometry, just as in dynamical Chern-Simons gravity.

preprint2014arXiv

Love can be Tough to Measure

The waveform phase for a neutron star binary can be split into point-particle terms and finite-size terms (characterized by the Love number) that account for equation of state effects. The latter first enter at 5 post-Newtonian (PN) order (i.e. proportional to the tenth power of the orbital velocity), but the former are only known completely to 3.5 PN order, with higher order terms only known to leading-order in the mass-ratio. We here find that not including point-particle terms at 4PN order to leading- and first-order in the mass ratio in the template model can severely deteriorate our ability to measure the equation of state. This problem can be solved if one uses numerical waveforms once their own systematic errors are under control.

preprint2014arXiv

Metric of a tidally perturbed spinning black hole

We explicitly construct the metric of a Kerr black hole that is tidally perturbed by the external universe in the slow-motion approximation. This approximation assumes that the external universe changes slowly relative to the rotation rate of the hole, thus allowing the parameterization of the Newman-Penrose scalar $ψ_0$ by time-dependent electric and magnetic tidal tensors. This approximation, however, does not constrain how big the spin of the background hole can be and, in principle, the perturbed metric can model rapidly spinning holes. We first generate a potential by acting with a differential operator on $ψ_0$. From this potential we arrive at the metric perturbation by use of the Chrzanowski procedure in the ingoing radiation gauge. We provide explicit analytic formulae for this metric perturbation in spherical Kerr-Schild coordinates, where the perturbation is finite at the horizon. This perturbation is parametrized by the mass and Kerr spin parameter of the background hole together with the electric and magnetic tidal tensors that describe the time evolution of the perturbation produced by the external universe. In order to take the metric accurate far away from the hole, these tidal tensors should be determined by asymptotically matching this metric to another one valid far from the hole. The tidally perturbed metric constructed here could be useful in initial data constructions to describe the metric near the horizons of a binary system of spinning holes. This perturbed metric could also be used to construct waveforms and study the absorption of mass and angular momentum by a Kerr black hole when external processes generate gravitational radiation.

preprint2014arXiv

Mis-Modelling in Gravitational Wave Astronomy: The Trouble With Templates

Waveform templates are a powerful tool for extracting and characterizing gravitational wave signals, acting as highly restrictive priors on the signal morphologies that allow us to extract weak events buried deep in the instrumental noise. The templates map the waveform shapes to physical parameters, thus allowing us to produce posterior probability distributions for these parameters. However, there are attendant dangers in using highly restrictive signal priors. If strong field gravity is not accurately described by General Relativity (GR), then using GR templates may result in fundamental bias in the recovered parameters, or even worse, a complete failure to detect signals. Here we study such dangers, concentrating on three distinct possibilities. First, we show that there exist modified theories compatible with all existing tests that would fail to be detected by the LIGO/Virgo network using searches based on GR templates, but which would be detected using a one parameter post-Einsteinian extension. Second, we study modified theories that produce departures from GR that turn on suddenly at a critical frequency, producing waveforms that do not naively fit into the simplest parameterized post-Einsteinian (ppE) scheme. We show that even the simplest ppE templates are still capable of picking up these strange signals and diagnosing a departure from GR. Third, we study whether using inspiral-only ppE waveforms for signals that include merger and ringdown can lead to problems in misidentifying a GR departure. We present an easy technique that allows us to self-consistently identify the inspiral portion of the signal, and thus remove these potential biases, allowing GR tests to be performed on higher mass signals that merge within the detector band. We close by studying a parameterized waveform model that may allow us to test GR using the full inspiral-merger-ringdown signal.

preprint2014arXiv

Post-Circular Expansion of Eccentric Binary Inspirals: Fourier-Domain Waveforms in the Stationary Phase Approximation

We lay the foundations for the construction of analytic expressions for Fourier-domain gravitational waveforms produced by eccentric, inspiraling compact binaries in a post-circular or small-eccentricity approximation. The time-dependent, "plus" and "cross" polarizations are expanded in Bessel functions, which are then self-consistently re-expanded in a power series about zero initial eccentricity to eighth order. The stationary phase approximation is then employed to obtain explicit analytic expressions for the Fourier transform of the post-circular expanded, time-domain signal. We exemplify this framework by considering Newtonian-accurate waveforms, which in the post-circular scheme give rise to higher harmonics of the orbital phase and amplitude corrections both to the amplitude and the phase of the Fourier domain waveform. Such higher harmonics lead to an effective increase in the inspiral mass reach of a detector as a function of the binary's eccentricity e_0 at the time when the binary enters the detector sensitivity band. Using the largest initial eccentricity allowed by our approximations (e_0 < 0.4), the mass reach is found to be enhanced up to factors of approximately 5 relative to that of circular binaries for Advanced LIGO, LISA, and the proposed Einstein Telescope at a signal-to-noise ratio of ten. A post-Newtonian generalization of the post circular scheme is also discussed, which holds the promise to provide "ready-to-use" Fourier-domain waveforms for data analysis of eccentric inspirals.

preprint2014arXiv

Projected Constraints on Lorentz-Violating Gravity with Gravitational Waves

Gravitational waves are excellent tools to probe the foundations of General Relativity in the strongly dynamical and non-linear regime. One such foundation is Lorentz symmetry, which can be broken in the gravitational sector by the existence of a preferred time direction, and thus, a preferred frame at each spacetime point. This leads to a modification in the orbital decay rate of binary systems, and also in the generation and chirping of their associated gravitational waves. We here study whether waves emitted in the late, quasi-circular inspiral of non-spinning, neutron star binaries can place competitive constraints on two proxies of gravitational Lorentz-violation: Einstein-Æther theory and khronometric gravity. We model the waves in the small-coupling (or decoupling) limit and in the post-Newtonian approximation, by perturbatively solving the field equations in small deformations from General Relativity and in the small-velocity/weak-gravity approximation. We assume a gravitational wave consistent with General Relativity has been detected with second- and third-generation, ground-based detectors, and with the proposed space-based mission, DECIGO, with and without coincident electromagnetic counterparts. Without a counterpart, a detection consistent with General Relativity of neutron star binaries can only place competitive constraints on gravitational Lorentz violation when using future, third-generation or space-based instruments. On the other hand, a single counterpart is enough to place constraints that are 10 orders of magnitude more stringent than current binary pulsar bounds, even when using second-generation detectors. This is because Lorentz violation forces the group velocity of gravitational waves to be different from that of light, and this difference can be very accurately constrained with coincident observations.

preprint2014arXiv

Projected Constraints on Scalarization with Gravitational Waves from Neutron Star Binaries

Certain scalar-tensor theories have the property of endowing stars with scalar hair, sourced either by the star's own compactness (spontaneous scalarization) or, for binary systems, by the companion's scalar hair (induced scalarization) or by the orbital binding energy (dynamical scalarization). Scalarized stars in binaries present different conservative dynamics than in General Relativity, and can also excite a scalar mode in the metric perturbation that carries away dipolar radiation. As a result, the binary orbit shrinks faster than predicted in General Relativity, modifying the rate of decay of the orbital period. In spite of this, scalar-tensor theories can pass existing binary pulsar tests, because observed pulsars may not be compact enough or sufficiently orbitally bound to activate scalarization. Gravitational waves emitted during the last stages of compact binary inspirals are thus ideal probes of scalarization effects. For the standard projected sensitivity of advanced LIGO, we here show that, if neutron stars are sufficiently compact to enter the detector's sensitivity band already scalarized, then gravitational waves could place constraints at least comparable to binary pulsars. If the stars dynamically scalarize while inspiraling in band, then constraints are still possible provided the scalarization occurs sufficiently early in the inspiral, roughly below an orbital frequency of 50Hz. In performing these studies, we derive an easy-to-calculate data analysis measure, an integrated phase difference between a General Relativistic and a modified signal, that maps directly to the Bayes factor so as to determine whether a modified gravity effect is detectable. Finally, we find that custom-made templates are equally effective as model-independent, parameterized post-Einsteinian waveforms at detecting such modified gravity effects at realistic signal-to-noise ratios.

preprint2014arXiv

Spin-Precession: Breaking the Black Hole--Neutron Star Degeneracy

Mergers of compact stellar remnant are prime targets for the LIGO/Virgo gravitational wave detectors. One hopes that the gravitational wave signals from these merger events can be used to study the mass and spin distribution of stellar remnants, and provide information about black hole horizons and the material properties of neutron stars. However, it has been suggested that degeneracies in the way that the star's mass and spin are imprinted in the waveforms may make it impossible to distinguish between black holes and neutron stars. Here we show that the precession of the orbital plane due to spin-orbit coupling breaks the mass-spin degeneracy, and allows us to distinguish between standard neutron stars and alternative possibilities, such as black holes or exotic neutron stars with large masses and spins.

preprint2014arXiv

Strong Binary Pulsar Constraints on Lorentz Violation in Gravity

Binary pulsars are excellent laboratories to test the building blocks of Einstein's theory of General Relativity. One of these is Lorentz symmetry which states that physical phenomena appear the same for all inertially moving observers. We study the effect of violations of Lorentz symmetry in the orbital evolution of binary pulsars and find that it induces a much more rapid decay of the binary's orbital period due to the emission of dipolar radiation. The absence of such behavior in recent observations allows us to place the most stringent constraints on Lorentz violation in gravity, thus verifying one of the cornerstones of Einstein's theory much more accurately than any previous gravitational observation.

preprint2014arXiv

Three-Hair Relations for Rotating Stars: Nonrelativistic Limit

The gravitational field outside of astrophysical black holes is completely described by their mass and spin frequency, as expressed by the no-hair theorems. These theorems assume vacuum spacetimes, and thus they apply only to black holes and not to stars. Despite this, we analytically find that the gravitational potential of arbitrarily rapid rigidly rotating stars can still be described completely by only their mass, spin angular momentum, and quadrupole moment. Although these results are obtained in the nonrelativistic limit (to leading order in a weak-field expansion of general relativity, GR), they are also consistent with fully relativistic numerical calculations of rotating neutron stars. This description of the gravitational potential outside the source in terms of just three quantities is approximately universal (independent of equation of state). Such universality may be used to break degeneracies in pulsar and future gravitational wave observations to extract more physics and test GR in the strong-field regime.

preprint2014arXiv

Toward Realistic and Practical No-Hair Relations for Neutron Stars in the Non-Relativistic Limit

The gravitational properties of astrophysical objects depend sensitively on their internal structure. In Newtonian theory, the gravitational potential of a rotating star can be fully described by an infinite number of multipole moments of its mass distribution. Recently, this infinite number of moments for uniformly-rotating stars were shown semi-analytically to be expressible in terms of just the first three: the mass, the spin, and the quadrupole moment of the star. The relations between the various lower multipole moments were additionally shown to depend weakly on the equation of state, when considering neutron stars and assuming single polytropic equations of state. Here we extend this result in two ways. First, we show that the universality also holds for realistic equations of state, thus relaxing the need to use single polytropes. Second, we derive purely analytical universal relations by perturbing the equations of structure about an $n=0$ polytrope that reproduce semi-analytic results to $\mathcal{O}(1\%)$. We also find that the linear-order perturbation vanishes in some cases, which provides further evidence and a deeper understanding of the universality.

preprint2014arXiv

Why I-Love-Q

Black holes are said to have no hair because all of their multipole moments can be expressed in terms of just their mass, charge and spin angular momentum. The recent discovery of approximately equation-of-state-independent relations among certain multipole moments in neutron stars suggests that they are also approximately bald. We here explore the yet unknown origin for this universality. First, we investigate which region of the neutron star's interior and of the equation of state is most responsible for the universality. We find that the universal relation between the moment of inertia and the quadrupole moment is dominated by the star's outer-core, a shell of width (50-95)% of the total radius, which corresponds to the density range $(10^{14}$-$10^{15})$g/cm$^3$. Second, we study the impact on the universality of approximating stellar isodensity contours as self-similar ellipsoids. An analytical calculation in the Newtonian limit reveals that the shape of the ellipsoids does not affect the universal relations, but relaxing the self-similarity assumption can completely destroy it. Third, we investigate the eccentricity profiles of rotating relativistic stars and find that the ellipticity is roughly constant, with variations of roughly (20-30)% in the region that matters to the universal relations. Fourth, we repeat the above analysis for non-compact, regular stars and find that the ellipticity is not constant, with variations that easily exceed 100% and universality is lost. These findings suggest that universality arises as an emergent approximate symmetry: as one flows in the stellar-structure phase space from non-compact star region to the relativistic star region, the eccentricity variation inside stars decreases, leading to approximate self-similarity in their isodensity contours, which then leads to the universal behavior observed in their exterior multipole moments.

preprint2013arXiv

A Rosetta Stone for Parameterized Tests of Gravity

Several model-independent parameterizations of deviations from General Relativity have been developed to test Einstein's theory. Although these different parameterizations were developed for different gravitational observables, they ultimately all test the same underlying physics. In this paper, we develop connections between the parameterized post-Newtonian, parameterized post-Keplerian, and the parameterized post-Einsteinian frameworks, developed to carry out tests of General Relativity with Solar System, binary pulsar, and gravitational wave observations respectively. These connections allow us to use knowledge gained from one framework to inform and guide tests using the others. Relating these parameterizations and combining the results from each approach strengthens our tests of General Relativity.

preprint2013arXiv

Gravitational Wave Tests of General Relativity with Ground-Based Detectors and Pulsar Timing Arrays

This review is focused on tests of Einstein's theory of General Relativity with gravitational waves that are detectable by ground-based interferometers and pulsar timing experiments. Einstein's theory has been greatly constrained in the quasi-linear, quasi-stationary regime, where gravity is weak and velocities are small. Gravitational waves will allow us to probe a complimentary, yet previously unexplored regime: the non-linear and dynamical strong-field regime. Such a regime is, for example, applicable to compact binaries coalescing, where characteristic velocities can reach fifty percent the speed of light and compactnesses can reach a half. This review begins with the theoretical basis and the predicted gravitational wave observables of modified gravity theories. The review continues with a brief description of the detectors, including both gravitational wave interferometers and pulsar timing arrays, leading to a discussion of the data analysis formalism that is applicable for such tests. The review ends with a discussion of gravitational wave tests for compact binary systems.

preprint2013arXiv

Gravitational Wave Tests of Strong Field General Relativity with Binary Inspirals: Realistic Injections and Optimal Model Selection

We study generic tests of strong-field General Relativity using gravitational waves emitted during the inspiral of compact binaries. Previous studies have considered simple extensions to the standard post-Newtonian waveforms that differ by a single term in the phase. Here we improve on these studies by (i) increasing the realism of injections and (ii) determining the optimal waveform families for detecting and characterizing such signals. We construct waveforms that deviate from those in General Relativity through a series of post-Newtonian terms, and find that these higher-order terms can affect our ability to test General Relativity, in some cases by making it easier to detect a deviation, and in some cases by making it more difficult. We find that simple single-phase post-Einsteinian waveforms are sufficient for detecting deviations from General Relativity, and there is little to be gained from using more complicated models with multiple phase terms. The results found here will help guide future attempts to test General Relativity with advanced ground-based detectors.

preprint2013arXiv

Gravitational Waveforms for Precessing, Quasicircular Compact Binaries with Multiple Scale Analysis: Small Spin Expansion

We obtain analytical gravitational waveforms in the frequency-domain for precessing, quasi-circular compact binaries with small spins, applicable, for example, to binary neutron star inspirals. We begin by calculating an analytic solution to the precession equations, obtained by expanding in the dimensionless spin parameters and using multiple-scale analysis to separate timescales. We proceed by analytically computing the Fourier transform of time-domain waveform through the stationary phase approximation. We show that the latter is valid for systems with small spins. Finally, we show that these waveforms have a high overlap with numerical waveforms obtained through direct integration of the precession equations and discrete Fourier transformations. The resulting, analytic waveform family is ideal for detection and parameter estimation of gravitational waves emitted by inspiraling binary neutron stars with ground-based detectors.

preprint2013arXiv

I-Love-Q

Neutron stars and quark stars are not only characterized by their mass and radius, but also by how fast they spin, through their moment of inertia, and how much they can be deformed, through their Love number and quadrupole moment. These depend sensitively on the star's internal structure, and thus on unknown nuclear physics. We find universal relations between the moment of inertia, the Love number and the quadrupole moment that are independent of the neutron star's and quark star's internal structure. These can be used to learn about the deformability of these compact objects through observations of the moment of inertia, break degeneracies in gravitational wave detection to measure spin in binary inspirals and test General Relativity in a nuclear-structure independent fashion.

preprint2013arXiv

Linear Stability Analysis and the Speed of Gravitational Waves in Dynamical Chern-Simons Modified Gravity

We perform a linear stability analysis of dynamical Chern-Simons modified gravity in the geometric optics approximation and find that it is linearly stable on the backgrounds considered. Our analysis also reveals that gravitational waves in the modified theory travel at the speed of light in Minkowski spacetime. However, on a Schwarzschild background the characteristic speed of propagation along a given direction splits into two modes, one subluminal and one superluminal. The width of the splitting depends on the azimuthal components of the propagation vector, is linearly proportional to the mass of the black hole, and decreases with the third inverse power of the distance from the black hole. Radial propagation is unaffected, implying that as probed by gravitational waves the location of the event horizon of the spacetime is unaltered. The analysis further reveals that when a high frequency, pure gravitational wave is scattered from a black hole, a scalar wave of comparable amplitude is excited, and vice-versa.

preprint2013arXiv

Resonant Post-Newtonian Eccentricity Excitation in Hierarchical Three-body Systems

We study the secular, hierarchical three-body problem to first-order in a post-Newtonian expansion of General Relativity. We expand the first-order post-Newtonian Hamiltonian to leading-order in the ratio of the semi-major axis of the two orbits. In addition to the well-known terms that correspond to the GR precession of the inner and outer orbits, we find a new secular post-Newtonian interaction term that can affect the long-term evolution of the triple. We explore the parameter space for highly inclined and eccentric systems, where the Kozai-Lidov mechanism can produce large-amplitude oscillations in the eccentricities. The standard lore, i.e., that General Relativity effects suppress eccentricity, is only consistent with the parts of phase space where the General Relativity timescales are several orders of magnitude shorter than the secular Newtonian one. In other parts of phase space, however, post-Newtonian corrections combined with the three body ones, can excite eccentricities. In particular, for systems where the General Relativity timescale is comparable to the secular Newtonian timescales, the three-body interactions give rise to a resonant-like eccentricity excitation. Furthermore, for triples with a comparable-mass inner binary, where the eccentric Kozai-Lidov mechanism is suppressed, post-Newtonian corrections can further increase the eccentricity and lead to orbital flips even when the timescale of the former is much longer than the timescale of the secular Kozai-Lidov quadrupole perturbations.

preprint2013arXiv

Spontaneous Generation of Angular Momentum in Holographic Theories

The Schwarzschild black two-brane in four-dimensional anti-de Sitter space is dual to a finite temperature state in three-dimensional conformal field theory. We show that the solution acquires a nonzero angular momentum density when a gravitational Chern-Simons coupling is turned on in the bulk, even though the solution is not modified. A similar phenomenon is found for the Reissner-Nordstrom black two-brane with axionic coupling to the gauge field. We discuss interpretation of this phenomenon from the point of view of the boundary three-dimensional conformal field theory.

preprint2013arXiv

Stealth Bias in Gravitational-Wave Parameter Estimation

Inspiraling binaries of compact objects are primary targets for current and future gravitational-wave observatories. Waveforms computed in General Relativity are used to search for these sources, and will probably be used to extract source parameters from detected signals. However, if a different theory of gravity happens to be correct in the strong-field regime, source-parameter estimation may be affected by a fundamental bias: that is, by systematic errors induced due to the use of waveforms derived in the incorrect theory. If the deviations from General Relativity are not large enough to be detectable on their own and yet these systematic errors remain significant (i.e., larger than the statistical uncertainties in parameter estimation), fundamental bias cannot be corrected in a single observation, and becomes stealth bias. In this article we develop a scheme to determine in which cases stealth bias could be present in gravitational-wave astronomy. For a given observation, the answer depends on the detection signal-to-noise ratio and on the strength of the modified-gravity correction. As an example, we study three representative stellar-mass binary systems that will be detectable with second-generation ground-based observatories. We find that significant systematic bias can occur whether or not modified gravity can be positively detected, for correction strengths that are not currently excluded by any other experiment. Thus, stealth bias may be a generic feature of gravitational-wave detections, and it should be considered and characterized, using expanded models such as the parametrized post-Einstein framework, when interpreting the results of parameter-estimation analyses.

preprint2013arXiv

Tidal heating and torquing of a Kerr black hole to next-to-leading order in the tidal coupling

We calculate the linear vacuum perturbations of a Kerr black hole surrounded by a slowly-varying external spacetime to third order in the ratio of the black-hole mass to the radius of curvature of the external spacetime. This expansion applies to two relevant physical scenarios: (i) a small Kerr black hole immersed in the gravitational field of a much larger external black hole, and (ii) a Kerr black hole moving slowly around another external black hole of comparable mass. This small-hole/slow-motion approximation allows us to parametrize the perturbation through slowly-varying, time-dependent electric and magnetic tidal tensors, which then enables us to separate the Teukolsky equation and compute the Newman-Penrose scalar analytically to third order in our expansion parameter. We obtain generic expressions for the mass and angular momentum flux through the perturbed black hole horizon, as well as the rate of change of the horizon surface area, in terms of certain invariants constructed from the electric and magnetic tidal tensors. We conclude by applying these results to the two scenarios described above.

preprint2012arXiv

Approximate Waveforms for Extreme-Mass-Ratio Inspirals: The Chimera Scheme

We describe a new kludge scheme to model the dynamics of generic extreme-mass-ratio inspirals (EMRIs; stellar compact objects spiraling into a spinning supermassive black hole) and their gravitational-wave emission. The Chimera scheme is a hybrid method that combines tools from different approximation techniques in General Relativity: (i) A multipolar, post-Minkowskian expansion for the far-zone metric perturbation (the gravitational waveforms) and for the local prescription of the self-force; (ii) a post-Newtonian expansion for the computation of the multipole moments in terms of the trajectories; and (iii) a BH perturbation theory expansion when treating the trajectories as a sequence of self-adjusting Kerr geodesics. The EMRI trajectory is made out of Kerr geodesic fragments joined via the method of osculating elements as dictated by the multipolar post-Minkowskian radiation-reaction prescription. We implemented the proper coordinate mapping between Boyer-Lindquist coordinates, associated with the Kerr geodesics, and harmonic coordinates, associated with the multipolar post-Minkowskian decomposition. The Chimera scheme is thus a combination of approximations that can be used to model generic inspirals of systems with extreme to intermediate mass ratios, and hence, it can provide valuable information for future space-based gravitational-wave observatories, like LISA, and even for advanced ground detectors. The local character in time of our multipolar post-Minkowskian self-force makes this scheme amenable to study the possible appearance of transient resonances in generic inspirals.

preprint2012arXiv

Asymptotically Matched Spacetime Metric for Non-Precessing, Spinning Black Hole Binaries

We construct a closed-form, fully analytical 4-metric that approximately represents the spacetime evolution of non-precessing, spinning black hole binaries from infinite separations up to a few orbits prior to merger. We employ the technique of asymptotic matching to join a perturbed Kerr metric in the neighborhood of each spinning black hole to a near-zone, post-Newtonian metric farther out. The latter is already naturally matched to a far-zone, post-Minkowskian metric that accounts for full temporal retardation. The result is a 4-metric that is approximately valid everywhere in space and in a small bundle of spatial hypersurfaces. We here restrict our attention to quasi- circular orbits, but the method is valid for any orbital motion or physical scenario, provided an overlapping region of validity or buffer zone exists. A simple extension of such a metric will allow for future studies of the accretion disk and jet dynamics around spinning back hole binaries.

preprint2012arXiv

Gravitational Waves from Quasicircular Extreme Mass-Ratio Inspirals as Probes of Scalar-Tensor Theories

A stellar-mass compact object spiraling into a supermassive black hole, an extreme-mass-ratio inspiral (EMRI), is one of the targets of future gravitational-wave detectors and it offers a unique opportunity to test General Relativity (GR) in the strong-field. We study whether generic scalar-tensor (ST) theories can be further constrained with EMRIs. We show that in the EMRI limit, all such theories universally reduce to massive or massless Brans-Dicke theory and that black holes do not emit dipolar radiation to all orders in post-Newtonian (PN) theory. For massless theories, we calculate the scalar energy flux in the Teukolsky formalism to all orders in PN theory and fit it to a high-order PN expansion. We derive the PN ST corrections to the Fourier transform of the gravitational wave response and map it to the parameterized post-Einsteinian framework. We use the effective-one-body framework adapted to EMRIs to calculate the ST modifications to the gravitational waveform. We find that such corrections are smaller than those induced in the early inspiral of comparable-mass binaries, leading to projected bounds on the coupling that are worse than current Solar System ones. Brans-Dicke theory modifies the weak-field, with deviations in the energy flux that are largest at small velocities. For massive theories, superradiance can lead to resonances in the scalar energy flux that can lead to floating orbits outside the innermost stable circular orbit and that last until the supermassive black hole loses enough mass and spin-angular momentum. If such floating orbits occur in the frequency band of LISA, they would lead to a large dephasing (~1e6 rads), preventing detection with GR templates. A detection that is consistent with GR would then rule out floating resonances at frequencies lower than the lowest observed frequency, allowing for the strongest constraints yet on massive ST theories.

preprint2012arXiv

NR/HEP: roadmap for the future

Physics in curved spacetime describes a multitude of phenomena, ranging from astrophysics to high energy physics. The last few years have witnessed further progress on several fronts, including the accurate numerical evolution of the gravitational field equations, which now allows highly nonlinear phenomena to be tamed. Numerical relativity simulations, originally developed to understand strong field astrophysical processes, could prove extremely useful to understand high-energy physics processes like trans-Planckian scattering and gauge-gravity dualities. We present a concise and comprehensive overview of the state-of-the-art and important open problems in the field(s), along with guidelines for the next years. This writeup is a summary of the "NR/HEP Workshop" held in Madeira, Portugal from August 31st to September 3rd 2011.

preprint2011arXiv

Accuracy of the Post-Newtonian Approximation: Optimal Asymptotic Expansion for Quasi-Circular, Extreme-Mass Ratio Inspirals

We study the accuracy of the post-Newtonian (PN) approximation and its formal region of validity, by investigating its optimal asymptotic expansion for the quasi-circular, adiabatic inspiral of a point particle into a Schwarzschild black hole. By comparing the PN expansion of the energy flux to numerical calculations in the perturbative Teukolsky formalism, we show that (i) the inclusion of higher multipoles is necessary to establish the accuracy of high-order PN terms, and (ii) the region of validity of PN theory is largest at relative O(1/c^6) (3PN order). The latter result suggests that the series diverges beyond 3PN order, at least in the extreme mass-ratio limit, probably due to the appearance of logarithmic terms in the energy flux. The study presented here is a first formal attempt to determine the region of validity of the PN approximation using asymptotic analysis. Therefore, it should serve as a template to perform similar studies on other systems, such as comparable-mass quasi-circular inspirals computed by high-accuracy numerical relativistic simulations.

preprint2011arXiv

Accuracy of the post-Newtonian approximation. II. Optimal asymptotic expansion of the energy flux for quasicircular, extreme mass-ratio inspirals into a Kerr black hole

We study the effect of black hole spin on the accuracy of the post-Newtonian approximation. We focus on the gravitational energy flux for the quasicircular, equatorial, extreme mass-ratio inspiral of a compact object into a Kerr black hole of mass M and spin J. For a given dimensionless spin a=J/M^2 (in geometrical units), the energy flux depends only on the orbital velocity v or (equivalently) on the Boyer-Lindquist orbital radius r. We investigate the formal region of validity of the Taylor post-Newtonian expansion of the energy flux (which is known up to order v^8 beyond the quadrupole formula), generalizing previous work by two of us. The "error function" used to determine the region of validity of the post-Newtonian expansion can have two qualitatively different kinds of behavior, and we deal with these two cases separately. We find that, at any fixed post-Newtonian order, the edge of the region of validity (as measured by v/v_{ISCO}, where v_{ISCO} is the orbital velocity at the innermost stable circular orbit) is only weakly dependent on a. Unlike in the nonspinning case, the lack of sufficiently high order terms does not allow us to determine if there is a convergent to divergent transition at order v^6. Independently of a, the inclusion of angular multipoles up to and including l=5 in the numerical flux is necessary to achieve the level of accuracy of the best-known (N=8) PN expansion of the energy flux.

preprint2011arXiv

Approximate Waveforms for Extreme-Mass-Ratio Inspirals in Modified Gravity Spacetimes

Extreme-mass-ratio inspirals, in which a stellar-mass compact object spirals into a supermassive black hole, are prime candidates for detection with space-borne milliHertz gravitational wave detectors, similar to the Laser Interferometer Space Antenna. The gravitational waves generated during such inspirals encode information about the background in which the small object is moving, providing a tracer of the spacetime geometry and a probe of strong-field physics. In this paper, we construct approximate, "analytic-kludge" waveforms for such inspirals with parameterized post-Einsteinian corrections that allow for generic, model-independent deformations of the supermassive black hole background away from the Kerr metric. These approximate waveforms include all of the qualitative features of true waveforms for generic inspirals, including orbital eccentricity and relativistic precession. The deformations of the Kerr metric are modeled using a recently proposed, modified gravity bumpy metric, which parametrically deforms the Kerr spacetime while ensuring that three approximate constants of the motion remain for geodesic orbits: a conserved energy, azimuthal angular momentum and Carter constant. The deformations represent modified gravity effects and have been analytically mapped to several modified gravity black hole solutions in four dimensions. In the analytic kludge waveforms, the conservative motion is modeled by a post-Newtonian expansion of the geodesic equations in the deformed spacetimes, which in turn induce modifications to the radiation-reaction force. These analytic-kludge waveforms serve as a first step toward complete and model-independent tests of General Relativity with extreme mass-ratio inspirals.

preprint2011arXiv

Effective Gravitational Wave Stress-energy Tensor in Alternative Theories of Gravity

The inspiral of binary systems in vacuum is controlled by the stress-energy of gravitational radiation and any other propagating degrees of freedom. For gravitational waves, the dominant contribution is characterized by an effective stress-energy tensor at future null infinity. We employ perturbation theory and the short-wavelength approximation to compute this stress-energy tensor in a wide class of alternative theories. We find that this tensor is generally a modification of that first computed by Isaacson, where the corrections can dominate over the general relativistic term. In a wide class of theories, however, these corrections identically vanish at asymptotically flat, future, null infinity, reducing the stress-energy tensor to Isaacson's. We exemplify this phenomenon by first considering dynamical Chern-Simons modified gravity, which corrects the action via a scalar field and the contraction of the Riemann tensor and its dual. We then consider a wide class of theories with dynamical scalar fields coupled to higher-order curvature invariants, and show that the gravitational wave stress-energy tensor still reduces to Isaacson's. The calculations presented in this paper are crucial to perform systematic tests of such modified gravity theories through the orbital decay of binary pulsars or through gravitational wave observations.

preprint2011arXiv

Gravitational Waves from Compact Binaries as Probes of the Universe

The future detection of gravitational wave forces us to consider the many ways in which astrophysics, gravitational wave theory and fundamental theory will interact. In this paper, I summarize some recent work done to develop such an interface. In particular, I concentrate on how non-vacuum astrophysical environments can modify the gravitational wave signal emitted by compact binary inspirals, and whether signatures from the former are detectable by current and future gravitational wave detectors. I also describe the interface between gravitational wave modeling and fundamental theory, focusing on the status of the parameterized post-Einsteinian framework (a general framework to detect deviations away from General Relativity in future gravitational wave data) and its current data analysis implementation.

preprint2011arXiv

Non-Spinning Black Holes in Alternative Theories of Gravity

We study two large classes of alternative theories, modifying the action through algebraic, quadratic curvature invariants coupled to scalar fields. We find one class that admits solutions that solve the vacuum Einstein equations and another that does not. In the latter, we find a deformation to the Schwarzschild metric that solves the modified field equations in the small coupling approximation. We calculate the event horizon shift, the innermost stable circular orbit shift, and corrections to gravitational waves, mapping them to the parametrized post-Einsteinian framework.

preprint2011arXiv

Observable Signatures of EMRI Black Hole Binaries Embedded in Thin Accretion Disks

We examine the electromagnetic (EM) and gravitational wave (GW) signatures of stellar-mass compact objects (COs) spiraling into a supermassive black hole (extreme mass-ratio inspirals or EMRIs), embedded in a thin, radiation-pressure dominated, accretion disk. At large separations, the tidal effect of the secondary CO clears a gap. We show that the gap refills during the late GW-driven phase of the inspiral, leading to a sudden EM brightening of the source. The accretion disk leaves an imprint on the GW through its angular momentum exchange with the binary, the mass increase of the binary members due to accretion, and its gravity. We compute the disk-modified GWs both in an analytical Newtonian approximation and in a numerical effective-one-body approach. We find that disk-induced migration provides the dominant perturbation to the inspiral, with weaker effects from the mass accretion onto the CO and hydrodynamic drag. Depending on whether a gap is present, the perturbation of the GW phase is between 10 and 1000 radians per year, detectable with the future Laser Interferometer Space Antenna (LISA) at high significance. The Fourier transform of the disk-modified GW in the stationary phase approximation is sensitive to disk parameters with a frequency trend different from post-Newtonian vacuum corrections. Our results suggest that observations of EMRIs may place new sensitive constraints on the physics of accretion disks.

preprint2011arXiv

Resonances in Extreme Mass-Ratio Inspirals: Asymptotic and Hyperasymptotic Analysis

An expected source of gravitational waves for future detectors in space are the inspirals of small compact objects into much more massive black holes. These sources have the potential to provide a wealth of information about astronomy and fundamental physics. On short timescales the orbit of the small object is approximately geodesic. Generic geodesics for a Kerr black hole spacetime have a complete set of integrals and can be characterized by three frequencies of the motion. Over the course of an inspiral, a typical system will pass through resonances where two of these frequencies become commensurate. The effect of the resonance will be to alter significantly the rate of inspiral for the duration of the resonance. Understanding the impact of these resonances on gravitational wave phasing is important to detect and exploit these signals for astrophysics and fundamental physics. Two differential equations that might describe the passage of an inspiral through such a resonance are investigated. These differ depending on whether it is the phase or the frequency components of a Fourier expansion of the motion that are taken to be continuous through the resonance. Asymptotic and hyperasymptotic analysis are used to find the late-time analytic behavior of the solution for a system that has passed through a resonance. Linearly growing (weak resonances) or linearly decaying (strong resonances) solutions are found depending on the strength of the resonance. In the weak-resonance case, frequency resonances leave an imprint (a resonant memory) on the gravitational frequency evolution. The transition between weak and strong resonances is characterized by a square-root singularity, and as one approaches this transition from above, the solutions to the frequency resonance equation bunch up into families exponentially fast.

preprint2011arXiv

The Effect of Massive Perturbers on Extreme Mass-Ratio Inspiral Waveforms

Extreme mass ratio inspirals, in which a stellar-mass object merges with a supermassive black hole, are prime sources for space-based gravitational wave detectors because they will facilitate tests of strong gravity and probe the spacetime around rotating compact objects. In the last few years of such inspirals, the total phase is in the millions of radians and details of the waveforms are sensitive to small perturbations. We show that one potentially detectable perturbation is the presence of a second supermassive black hole within a few tenths of a parsec. The acceleration produced by the perturber on the extreme mass-ratio system produces a steady drift that causes the waveform to deviate systematically from that of an isolated system. If the perturber is a few tenths of a parsec from the extreme-mass ratio system (plausible in as many as a few percent of cases) higher derivatives of motion might also be detectable. In that case, the mass and distance of the perturber can be derived independently, which would allow a new probe of merger dynamics.

preprint2010arXiv

A Tale of Two Jets

One of the most interesting high-energy, astrophysical phenomena are relativistic jets emitted from highly localized sky location. Such jets are common in Nature, observed to high redshift and in a range of wavelengths. Their precise generation mechanism remains a bit of a mystery, but they are generically believed to be powered by black holes. We here summarize the recent simulations of Palenzuela, Lehner and Liebling that shed light on the jet generation mechanism. These authors studied the merger of two non-spinning black holes in the presence of a magnetic field, perpendicular to the orbital plane and anchored by a circumbinary accretion disk, in the "force-free" approximation. They found that each black hole essentially acts as a "straw" that stirs the magnetic field lines around the center of mass as the black holes inspiral. The twisting of the magnetic field lines then generates jets around each black hole, even though these are not spinning. Their simulations show the formation of such a dual jet geometry and how it transitions to a single jet one, as the black holes merge due to gravitational wave emission.

preprint2010arXiv

Binary Pulsar Constraints on the Parameterized post-Einsteinian Framework

We constrain the parameterized post-Einsteinian framework with binary pulsar observations of orbital period decay due to gravitational wave emission. This framework proposes to enhance the amplitude and phase of gravitational waveform templates through post-Einsteinian parameters to search for generic deviations from General Relativity in gravitational wave data. Such enhancements interpolate between General Relativity and alternative theory predictions, but their magnitude must be such as to satisfy all current experiments and observations. The data that currently constrains the parameterized post-Einsteinian framework the most is the orbital period decay of binary pulsars. We use such observations to place upper limits on the magnitude of post-Einsteinian parameters, which will be critical when gravitational waves are detected and this framework is implemented.

preprint2010arXiv

Constraining Parity Violation in Gravity with Measurements of Neutron-Star Moments of Inertia

Neutron stars are sensitive laboratories for testing general relativity, especially when considering deviations where velocities are relativistic and gravitational fields are strong. One such deviation is described by dynamical, Chern-Simons modified gravity, where the Einstein-Hilbert action is modified through the addition of the gravitational parity-violating Pontryagin density coupled to a field. This four-dimensional effective theory arises naturally both in perturbative and non-perturbative string theory, loop quantum gravity, and generic effective field theory expansions. We calculate here Chern-Simons modifications to the properties and gravitational fields of slowly spinning neutron stars. We find that the Chern-Simons correction affects only the gravitomagnetic sector of the metric to leading order, thus introducing modifications to the moment of inertia but not to the mass-radius relation. We show that an observational determination of the moment of inertia to an accuracy of 10%, as is expected from near-future observations of the double pulsar, will place a constraint on the Chern-Simons coupling constant of ξ^{1/4} < 5 km, which is at least three-orders of magnitude stronger than the previous strongest bound.

preprint2010arXiv

Constraining the evolutionary history of Newton's constant with gravitational wave observations

Space-borne gravitational wave detectors, such as the proposed Laser Interferometer Space Antenna, are expected to observe black hole coalescences to high redshift and with large signal-to-noise ratios, rendering their gravitational waves ideal probes of fundamental physics. The promotion of Newton's constant to a time-function introduces modifications to the binary's binding energy and the gravitational wave luminosity, leading to corrections in the chirping frequency. Such corrections propagate into the response function and, given a gravitational wave observation, they allow for constraints on the first time-derivative of Newton's constant at the time of merger. We find that space-borne detectors could indeed place interesting constraints on this quantity as a function of sky position and redshift, providing a {\emph{constraint map}} over the entire range of redshifts where binary black hole mergers are expected to occur. A LISA observation of an equal-mass inspiral event with total redshifted mass of 10^5 solar masses for three years should be able to measure $\dot{G}/G$ at the time of merger to better than 10^(-11)/yr.

preprint2010arXiv

Extreme Mass-Ratio Inspirals in the Effective-One-Body Approach: Quasi-Circular, Equatorial Orbits around a Spinning Black Hole

We construct effective-one-body waveform models suitable for data analysis with LISA for extreme-mass ratio inspirals in quasi-circular, equatorial orbits about a spinning supermassive black hole. The accuracy of our model is established through comparisons against frequency-domain, Teukolsky-based waveforms in the radiative approximation. The calibration of eight high-order post-Newtonian parameters in the energy flux suffices to obtain a phase and fractional amplitude agreement of better than 1 radian and 1 % respectively over a period between 2 and 6 months depending on the system considered. This agreement translates into matches higher than 97 % over a period between 4 and 9 months, depending on the system. Better agreements can be obtained if a larger number of calibration parameters are included. Higher-order mass ratio terms in the effective-one-body Hamiltonian and radiation-reaction introduce phase corrections of at most 30 radians in a one year evolution. These corrections are usually one order of magnitude larger than those introduced by the spin of the small object in a one year evolution. These results suggest that the effective-one-body approach for extreme mass ratio inspirals is a good compromise between accuracy and computational price for LISA data analysis purposes.

preprint2010arXiv

Gravitational Waves Notes, Issue #2 : "A probe of spacetime and astrophysics: EMRIs"

GW Notes was born from the need for a journal where the distinct communities involved in gravitation wave research might gather. While these three communities - Astrophysics, General Relativity and Data Analysis - have made significant collaborative progress over recent years, we believe that it is indispensable to future advancement that they draw closer, and that they speak a common idiom. For this GW Notes issue we have approached Nicolás Yunes (Princeton University) to extend in high detail his recent work on EMRI waveforms for our highlight article.

preprint2010arXiv

Integrated Sachs-Wolfe Effect for Gravitational Radiation

Gravitational waves are messengers carrying valuable information about their sources. For sources at cosmological distances, the waves will contain also the imprint left by the intervening matter. The situation is in close analogy with cosmic microwave photons, for which the large-scale structures the photons traverse contribute to the observed temperature anisotropies, in a process known as the integrated Sachs-Wolfe effect. We derive the gravitational wave counterpart of this effect for waves propagating on a Friedman-Robertson-Walker background with scalar perturbations. We find that the phase, frequency and amplitude of the gravitational waves experience Sachs-Wolfe type integrated effects, this in addition to the magnification effects on the amplitude from gravitational lensing. We show that for supermassive black hole binaries, the integrated effects could account for measurable changes on the frequency, chirp mass and luminosity distance of the binary, thus unveiling the presence of inhomogeneities, and potentially dark energy, in the Universe.

preprint2010arXiv

Modeling Extreme Mass Ratio Inspirals within the Effective-One-Body Approach

We present the first models of extreme-mass-ratio inspirals within the effective-one-body (EOB) formalism, focusing on quasi-circular orbits into non-rotating black holes. We show that the phase difference and (Newtonian normalized) amplitude difference between analytical EOB and numerical Teukolsky-based gravitational waveforms can be reduced to less than 10^(-1) rad and less than 2 x 10^(-3), respectively, after a 2-year evolution. The inclusion of post-Newtonian self-force terms in the EOB approach leads to a phase disagreement of roughly 6-27 rad after a 2-year evolution. Such inclusion could also allow for the EOB modeling of waveforms from intermediate-mass ratio, quasi-circular inspirals.

preprint2010arXiv

Semianalytical estimates of scattering thresholds and gravitational radiation in ultrarelativistic black hole encounters

Ultrarelativistic collisions of black holes are ideal gedanken experiments to study the nonlinearities of general relativity. In this paper we use semianalytical tools to better understand the nature of these collisions and the emitted gravitational radiation. We explain many features of the energy spectra extracted from numerical relativity simulations using two complementary semianalytical calculations. In the first calculation we estimate the radiation by a "zero-frequency limit" analysis of the collision of two point particles with finite impact parameter. In the second calculation we replace one of the black holes by a point particle plunging with arbitrary energy and impact parameter into a Schwarzschild black hole, and we explore the multipolar structure of the radiation paying particular attention to the near-critical regime. We also use a geodesic analogy to provide qualitative estimates of the dependence of the scattering threshold on the black hole spin and on the dimensionality of the spacetime.

preprint2010arXiv

Superkicks in ultrarelativistic encounters of spinning black holes

We study ultrarelativistic encounters of two spinning, equal-mass black holes through simulations in full numerical relativity. Two initial data sequences are studied in detail: one that leads to scattering and one that leads to a grazing collision and merger. In all cases, the initial black hole spins lie in the orbital plane, a configuration that leads to the so-called "superkicks". In astrophysical, quasicircular inspirals, such kicks can be as large as ~3,000 km/s; here, we find configurations that exceed ~15,000 km/s. We find that the maximum recoil is to a good approximation proportional to the total amount of energy radiated in gravitational waves, but largely independent of whether a merger occurs or not. This shows that the mechanism predominantly responsible for the superkick is not related to merger dynamics. Rather, a consistent explanation is that the "bobbing" motion of the orbit causes an asymmetric beaming of the radiation produced by the in-plane orbital motion of the binary, and the net asymmetry is balanced by a recoil. We use our results to formulate some conjectures on the ultimate kick achievable in any black hole encounter.

preprint2010arXiv

Testing gravitational parity violation with coincident gravitational waves and short gamma-ray bursts

Gravitational parity violation is a possibility motivated by particle physics, string theory and loop quantum gravity. One effect of it is amplitude birefringence of gravitational waves, whereby left and right circularly-polarized waves propagate at the same speed but with different amplitude evolution. Here we propose a test of this effect through coincident observations of gravitational waves and short gamma-ray bursts from binary mergers involving neutron stars. Such gravitational waves are highly left or right circularly-polarized due to the geometry of the merger. Using localization information from the gamma-ray burst, ground-based gravitational wave detectors can measure the distance to the source with reasonable accuracy. An electromagnetic determination of the redshift from an afterglow or host galaxy yields an independent measure of this distance. Gravitational parity violation would manifest itself as a discrepancy between these two distance measurements. We exemplify such a test by considering one specific effective theory that leads to such gravitational parity-violation, Chern-Simons gravity. We show that the advanced LIGO-Virgo network and all-sky gamma-ray telescopes can be sensitive to the propagating sector of Chern-Simons gravitational parity violation to a level roughly two orders of magnitude better than current stationary constraints from the LAGEOS satellites.

preprint2010arXiv

Testing Modified Gravity with Gravitational Wave Astronomy

The emergent area of gravitational wave astronomy promises to provide revolutionary discoveries in the areas of astrophysics, cosmology, and fundamental physics. One of the most exciting possibilities is to use gravitational-wave observations to test alternative theories of gravity. In this contribution we describe how to use observations of extreme-mass-ratio inspirals by the future Laser Interferometer Space Antenna to test a particular class of theories: Chern-Simons modified gravity.

preprint2010arXiv

Towards Tests of Alternative Theories of Gravity with LISA

The inspiral of stellar compact objects into massive black holes, usually known as extreme-mass-ratio inspirals (EMRIs), is one of the most important sources of gravitational-waves for the future Laser Interferometer Space Antenna (LISA). Intermediate-mass-ratio inspirals (IMRIs are also of interest to advance ground-based gravitational-wave observatories. We discuss here how modifications to the gravitational interaction can affect the signals emitted by these systems and their detectability by LISA. We concentrate in particular on Chern-Simons modified gravity, a theory that emerges in different quantum gravitational approaches.

preprint2009arXiv

Chern-Simons Modified General Relativity

Chern-Simons modified gravity is an effective extension of general relativity that captures leading-order, gravitational parity violation. Such an effective theory is motivated by anomaly cancelation in particle physics and string theory. In this review, we begin by providing a pedagogical derivation of the three distinct ways such an extension arises: (1) in particle physics, (2) from string theory and (3) geometrically. We then review many exact and approximate, vacuum solutions of the modified theory, and discuss possible matter couplings. Following this, we review the myriad astrophysical, solar system, gravitational wave and cosmological probes that bound Chern-Simons modified gravity, including discussions of cosmic baryon asymmetry and inflation. The review closes with a discussion of possible future directions in which to test and study gravitational parity violation.

preprint2009arXiv

Conformally curved binary black hole initial data including tidal deformations and outgoing radiation

(Abridged) By asymptotically matching a post-Newtonian (PN) metric to two tidally perturbed Schwarzschild metrics, we generate approximate initial data (in the form of a 4-metric) for a nonspinning black hole binary in a circular orbit. We carry out this matching through O(v^4) in the binary's orbital velocity v, so the resulting data are conformally curved. Far from the holes, we use the appropriate PN metric that accounts for retardation, which we construct using the highest-order PN expressions available to compute the binary's past history. The data set's uncontrolled remainders are thus O(v^5) throughout the timeslice; we also generate an extension to the data set that has uncontrolled remainders of O(v^6) in the purely PN portion of the timeslice (i.e., not too close to the holes). The resulting data are smooth, since we join all the metrics together by smoothly interpolating between them. We perform this interpolation using transition functions constructed to avoid introducing excessive additional constraint violations. Due to their inclusion of tidal deformations and outgoing radiation, these data should substantially reduce the initial spurious ("junk") radiation observed in current simulations that use conformally flat initial data. Such reductions in the nonphysical components of the initial data will be necessary for simulations to achieve the accuracy required to supply Advanced LIGO and LISA with the templates necessary for parameter estimation.

preprint2009arXiv

Double Binary Pulsar Test of Dynamical Chern-Simons Modified Gravity

Chern-Simons modified gravity is a string-theory and loop-quantum-gravity inspired effective theory that modifies General Relativity by adding a parity-violating Pontryagin density to the Einstein-Hilbert action multiplied by a coupling scalar. We strongly constrain non-dynamical Chern-Simons modified gravity with a timelike Chern-Simons scalar through observations of the double binary pulsar PSR J0737-3039A/B. We first calculate Chern-Simons corrections to the orbital evolution of binary systems. We find that the ratio of the correction to periastron precession to the general relativistic prediction scales quadratically with the semi-major axis and inversely with the square of the object's radius. Binary pulsar systems are thus ideal to test this theory, since periastron precession can be measured with sub-degree accuracies and the semi-major axis is millions of times larger than the stellar radius. Using data from PSR J0737-3039A/B we dramatically constrain the non-dynamical Chern-Simons coupling to $M_{\CS} := 1/|\dot{\vartheta}| > 33 {\textrm{meV}}$, approximately a hundred billion times better than current Solar System tests.

preprint2009arXiv

Testing Effective Quantum Gravity with Gravitational Waves from Extreme-Mass-Ratio Inspirals

Testing deviation of GR is one of the main goals of the proposed {\emph{Laser Interferometer Space Antenna}}, a space-based gravitational-wave observatory. For the first time, we consistently compute the generation of gravitational waves from extreme-mass ratio inspirals (stellar compact objects into supermassive black holes) in a well-motivated alternative theory of gravity, that to date remains weakly constrained by double binary pulsar observations. The theory we concentrate on is Chern-Simons (CS) modified gravity, a 4-D, effective theory that is motivated both from string theory and loop-quantum gravity, and which enhances the Einstein-Hilbert action through the addition of a dynamical scalar field and the parity-violating Pontryagin density. We show that although point particles continue to follow geodesics in the modified theory, the background about which they inspiral is a modification to the Kerr metric, which imprints a CS correction on the gravitational waves emitted. CS modified gravitational waves are sufficiently different from the General Relativistic expectation that they lead to significant dephasing after 3 weeks of evolution, but such dephasing will probably not prevent detection of these signals, but instead lead to a systematic error in the determination of parameters. We end with a study of radiation-reaction in the modified theory and show that, to leading-order, energy-momentum emission is not CS modified, except possibly for the subdominant effect of scalar-field emission. The inclusion of radiation-reaction will allow for tests of CS modified gravity with space-borne detectors that might be two orders of magnitude larger than current binary pulsar bounds.

preprint2008arXiv

Constraining effective quantum gravity with LISA

All modern routes leading to a quantum theory of gravity -- i.e., perturbative quantum gravitational one-loop exact correction to the global chiral current in the standard model, string theory, and perhaps even loop quantum gravity -- require supplementing the Einstein-Hilbert action with a parity-violating Chern-Simons term. Such a term leads to amplitude-birefringent gravitational wave propagation: i.e., one (circular) polarization state amplified with propagation while the other is attenuated. The proposed Laser Interferometer Space Antenna (LISA) is capable of observing gravitational wave sources at cosmological distances, suggesting the possibility that LISA observations may place a strong bound on this manifestation of quantum gravity. Here we report on a calculation of the effect that spacetime amplitude birefringence has on the signal LISA is capable of observing from inspiraling supermassive black hole binaries at large redshift. We find that the birefringence manifests itself in the observations as an anomalous precession of the binary's orbital angular momentum as it evolves toward coalescence, whose magnitude depends on the integrated history of the Chern-Simons coupling over the worldline of radiation wavefront. We estimate that LISA could place bounds on Chern-Simons modified gravity that are several orders of magnitude stronger than the present Solar System constraints, thus providing a probe of the quantum structure of spacetime.

preprint2007arXiv

Frankenstein's Glue: Transition functions for approximate solutions

Approximations are commonly employed to find approximate solutions to the Einstein equations. These solutions, however, are usually only valid in some specific spacetime region. A global solution can be constructed by gluing approximate solutions together, but this procedure is difficult because discontinuities can arise, leading to large violations of the Einstein equations. In this paper, we provide an attempt to formalize this gluing scheme by studying transition functions that join approximate solutions together. In particular, we propose certain sufficient conditions on these functions and proof that these conditions guarantee that the joined solution still satisfies the Einstein equations to the same order as the approximate ones. An example is also provided for a binary system of non-spinning black holes, where the approximate solutions are taken to be given by a post-Newtonian expansion and a perturbed Schwarzschild solution. For this specific case, we show that if the transition functions satisfy the proposed conditions, the joined solution does not contain any violations to the Einstein equations larger than those already inherent in the approximations. We further show that if these functions violate the proposed conditions, then the matter content of the spacetime is modified by the introduction of a matter-shell, whose stress-energy tensor depends on derivatives of these functions.

preprint2007arXiv

Gravitational Recoil from Binary Black Hole Mergers: the Close-Limit Approximation

[abridged] The coalescence of a binary black hole system is one of the main sources of gravitational waves that present and future detectors will study. Apart from the energy and angular momentum that these waves carry, for unequal-mass binaries there is also a net flux of linear momentum that implies a recoil velocity of the resulting final black hole in the opposite direction. We present a computation of the recoil velocity based on the close-limit approximation scheme, which gives excellent results for head-on and grazing collisions of black holes when compared to full numerical relativistic calculations. We obtain a maximum recoil velocity of ~ 64 km/s for a symmetric mass ratio eta = M_1 M_2/(M_1+M_2)^2 ~ 0.19 and an initial proper separation of 4 M, where M is the total ADM mass of the system. This separation is the maximum at which the close-limit approximation is expected to provide accurate results. If we supplement this estimate with PN calculations up to the innermost stable circular orbit, we obtain a lower bound for the recoil velocity, with a maximum around 84 km/s. This is a lower bound because it neglects the initial merger phase. We can however obtain a rough estimate by using PN methods or the close-limit approximation. Since both methods are known to overestimate the amount of radiation, we obtain in this way an upper bound for the recoil with maxima in the range of 220-265 km/s. We also provide non-linear fits to these estimated upper and lower bounds. These estimates are subject to uncertainties related to issues such as the choice of initial data and higher effects in perturbation theory. Nonetheless, our estimates are consistent with previous results in the literature and suggest a narrower range of possible recoil velocities.