Showing 1–46 of 46 results for author: Stein, L C

Waveform Modelling for the Laser Interferometer Space Antenna

Authors: LISA Consortium Waveform Working Group, Niayesh Afshordi, Sarp Akçay, Pau Amaro Seoane, Andrea Antonelli, Josu C. Aurrekoetxea, Leor Barack, Enrico Barausse, Robert Benkel, Laura Bernard, Sebastiano Bernuzzi, Emanuele Berti, Matteo Bonetti, Béatrice Bonga, Gabriele Bozzola, Richard Brito, Alessandra Buonanno, Alejandro Cárdenas-Avendaño, Marc Casals, David F. Chernoff, Alvin J. K. Chua, Katy Clough, Marta Colleoni, Mekhi Dhesi, Adrien Druart , et al. (121 additional authors not shown)

Abstract: LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmologic… ▽ More LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome. △ Less

Submitted 20 December, 2023; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: 239 pages, 11 figures, white paper from the LISA Consortium Waveform Working Group, invited for submission to Living Reviews in Relativity, updated with comments from community

Closed-form solutions of spinning, eccentric binary black holes at 1.5 post-Newtonian order

Authors: Rickmoy Samanta, Sashwat Tanay, Leo C. Stein

Abstract: The closed-form solution of the 1.5 post-Newtonian (PN) accurate binary black hole (BBH) Hamiltonian system has proven to be difficult to obtain for a long time since its introduction in 1966. Closed-form solutions of the PN BBH systems with arbitrary parameters (masses, spins, eccentricity) are required for modeling the gravitational waves (GWs) emitted by them. Accurate models of GWs are crucial… ▽ More The closed-form solution of the 1.5 post-Newtonian (PN) accurate binary black hole (BBH) Hamiltonian system has proven to be difficult to obtain for a long time since its introduction in 1966. Closed-form solutions of the PN BBH systems with arbitrary parameters (masses, spins, eccentricity) are required for modeling the gravitational waves (GWs) emitted by them. Accurate models of GWs are crucial for their detection by LIGO/Virgo and LISA. Only recently, two solution methods for solving the BBH dynamics were proposed in arXiv:1908.02927 (without using action-angle variables), and arXiv:2012.06586, arXiv:2110.15351 (action-angle based). This paper combines the ideas laid out in the above articles, fills the missing gaps and provides the two solutions which are fully 1.5PN accurate. We also present a public Mathematica package BBHpnToolkit which implements these two solutions and compares them with a fully numerical treatment. The level of agreement between these solutions provides a numerical verification for all the five actions constructed in arXiv:2012.06586, and arXiv:2110.15351. This paper hence serves as a stepping stone for pushing the action-angle-based solution to 2PN order via canonical perturbation theory. △ Less

Submitted 14 December, 2023; v1 submitted 4 October, 2022; originally announced October 2022.

Comments: 14 pages, 3 figures

Journal ref: Phys. Rev. D 108, 124039; Published 14 December 2023

Nonlinearities in Black Hole Ringdowns

Authors: Keefe Mitman, Macarena Lagos, Leo C. Stein, Sizheng Ma, Lam Hui, Yanbei Chen, Nils Deppe, François Hébert, Lawrence E. Kidder, Jordan Moxon, Mark A. Scheel, Saul A. Teukolsky, William Throwe, Nils L. Vu

Abstract: The gravitational wave strain emitted by a perturbed black hole (BH) ringing down is typically modeled analytically using first-order BH perturbation theory. In this Letter we show that second-order effects are necessary for modeling ringdowns from BH merger simulations. Focusing on the strain's $(\ell,m)=(4,4)$ angular harmonic, we show the presence of a quadratic effect across a range of binary… ▽ More The gravitational wave strain emitted by a perturbed black hole (BH) ringing down is typically modeled analytically using first-order BH perturbation theory. In this Letter we show that second-order effects are necessary for modeling ringdowns from BH merger simulations. Focusing on the strain's $(\ell,m)=(4,4)$ angular harmonic, we show the presence of a quadratic effect across a range of binary BH mass ratios that agrees with theoretical expectations. We find that the quadratic $(4,4)$ mode's amplitude exhibits quadratic scaling with the fundamental $(2,2)$ mode -- its parent mode. The nonlinear mode's amplitude is comparable to or even larger than that of the linear $(4,4)$ mode. Therefore, correctly modeling the ringdown of higher harmonics -- improving mode mismatches by up to 2 orders of magnitude -- requires the inclusion of nonlinear effects. △ Less

Submitted 22 February, 2023; v1 submitted 15 August, 2022; originally announced August 2022.

Comments: 6+2 pages, 4 figures, 1 table. Matches PRL version

Journal ref: Phys. Rev. Lett. 130, 081402 (2023)

Tidally-induced nonlinear resonances in EMRIs with an analogue model

Authors: David Bronicki, Alejandro Cárdenas-Avendaño, Leo C. Stein

Abstract: One of the important targets for the future space-based gravitational wave observatory LISA is extreme mass ratio inspirals (EMRIs), where long and accurate waveform modeling is necessary for detection and characterization. Modeling EMRI dynamics requires accounting for effects such as the ones induced by an external tidal field, which can break integrability at resonances and cause significant de… ▽ More One of the important targets for the future space-based gravitational wave observatory LISA is extreme mass ratio inspirals (EMRIs), where long and accurate waveform modeling is necessary for detection and characterization. Modeling EMRI dynamics requires accounting for effects such as the ones induced by an external tidal field, which can break integrability at resonances and cause significant dephasing. In this paper, we use a Newtonian analogue of a Kerr black hole to study the effect of an external tidal field on the dynamics and the gravitational waveform. We have developed a numerical framework that takes advantage of the integrability of the background system to evolve it with a symplectic splitting integrator, and compute approximate gravitational waveforms to estimate the timescale over which the perturbation affects the dynamics. Comparing this timescale with the characteristic time under radiation reaction at resonance, we introduce a tool for quantifying the regime in which tidal effects might be included when modeling EMRI gravitational waves. As an application of this framework, we perform a detailed analysis of the dynamics at one resonance to show how different entry points into the resonance in phase-space can produce substantially different dynamics, and how one can estimate bounds for the parameter space where tidal effects may become dominant. Such bounds will scale as $\varepsilon \gtrsim C \, q$, where $\varepsilon$ measures the strength of the external tidal field, $q$ is the mass ratio, and $C$ is a number which depends on the resonance and the shape of the tide. We demonstrate how to estimate $C$ using our framework for the 2:3 radial to polar frequency resonance in our model system. This framework can serve as a proxy for proper modeling of the tidal perturbation in the fully relativistic case. △ Less

Submitted 10 October, 2023; v1 submitted 16 March, 2022; originally announced March 2022.

Comments: V2. Updated to match published version

Journal ref: Class. Quantum Grav. 40 215015, 2023

High Precision Ringdown Modeling: Multimode Fits and BMS Frames

Authors: Lorena Magaña Zertuche, Keefe Mitman, Neev Khera, Leo C. Stein, Michael Boyle, Nils Deppe, François Hébert, Dante A. B. Iozzo, Lawrence E. Kidder, Jordan Moxon, Harald P. Pfeiffer, Mark A. Scheel, Saul A. Teukolsky, William Throwe, Nils Vu

Abstract: Quasi-normal mode (QNM) modeling is an invaluable tool for characterizing remnant black holes, studying strong gravity, and testing GR. Only recently have QNM studies begun to focus on multimode fitting to numerical relativity (NR) strain waveforms. As GW observatories become even more sensitive they will be able to resolve higher-order modes. Consequently, multimode QNM fits will be critically im… ▽ More Quasi-normal mode (QNM) modeling is an invaluable tool for characterizing remnant black holes, studying strong gravity, and testing GR. Only recently have QNM studies begun to focus on multimode fitting to numerical relativity (NR) strain waveforms. As GW observatories become even more sensitive they will be able to resolve higher-order modes. Consequently, multimode QNM fits will be critically important, and in turn require a more thorough treatment of the asymptotic frame at $\mathscr{I}^+$. The first main result of this work is a method for systematically fitting a QNM model containing many modes to a numerical waveform produced using Cauchy-characteristic extraction (CCE), an extraction technique which is known to resolve memory effects. We choose the modes to model based on their power contribution to the residual between numerical and model waveforms. We show that the all-mode strain mismatch improves by a factor of $\sim10^5$ when using multimode fitting as opposed to only fitting the $(2,\pm2,n)$ modes. Our most significant result addresses a critical point that has been overlooked in the QNM literature: the importance of matching the Bondi-van der Burg-Metzner-Sachs (BMS) frame of the numerical waveform to that of the QNM model. We show that by mapping the numerical waveforms$-$which exhibit the memory effect$-$to a BMS frame known as the super rest frame, there is an improvement of $\sim10^5$ in the all-mode strain mismatch compared to using a strain waveform whose BMS frame is not fixed. Furthermore, we find that by mapping CCE waveforms to the super rest frame, we can obtain all-mode mismatches that are, on average, a factor of $\sim4$ better than using the publicly-available extrapolated waveforms. We illustrate the effectiveness of these modeling enhancements by applying them to families of waveforms produced by NR and comparing our results to previous QNM studies. △ Less

Submitted 11 May, 2022; v1 submitted 29 October, 2021; originally announced October 2021.

Comments: 17 + 2 pages, 11 figures, 2 tables

Journal ref: Phys. Rev. D 105, 104015 (2022)

Action-angle variables of a binary black-hole with arbitrary eccentricity, spins, and masses at 1.5 post-Newtonian order

Authors: Sashwat Tanay, Gihyuk Cho, Leo C. Stein

Abstract: Accurate and efficient modeling of the dynamics of binary black holes (BBHs) is crucial to their detection through gravitational waves (GWs), with LIGO/Virgo/KAGRA, and LISA in the future. Solving the dynamics of a BBH system with arbitrary parameters without simplifications (like orbit- or precession-averaging) in closed form is one of the most challenging problems for the GW community. One poten… ▽ More Accurate and efficient modeling of the dynamics of binary black holes (BBHs) is crucial to their detection through gravitational waves (GWs), with LIGO/Virgo/KAGRA, and LISA in the future. Solving the dynamics of a BBH system with arbitrary parameters without simplifications (like orbit- or precession-averaging) in closed form is one of the most challenging problems for the GW community. One potential approach is using canonical perturbation theory which constructs perturbed action-angle variables from the unperturbed ones of an integrable Hamiltonian system. Having action-angle variables of the integrable 1.5 post-Newtonian (PN) BBH system is therefore imperative. In this paper, we continue the work initiated by two of us in arXiv:2012.06586, where we presented four out of five actions of a BBH system with arbitrary eccentricity, masses, and spins, at 1.5PN order. Here we compute the remaining fifth action using a novel method of extending the phase space by introducing unmeasurable phase space coordinates. We detail how to compute all the frequencies, and sketch how to explicitly transform from the action-angle variables to the usual positions and momenta. This analytically solves the dynamics at 1.5PN. This lays the groundwork to analytically solve the conservative dynamics of the BBH system with arbitrary masses, spins, and eccentricity, at higher PN order, by using canonical perturbation theory. △ Less

Submitted 23 May, 2023; v1 submitted 28 October, 2021; originally announced October 2021.

Comments: 26 pages, 6 figures

Journal ref: Phys. Rev. D 107, 103040 (2023)

Constraining gravitational wave amplitude birefringence and Chern-Simons gravity with GWTC-2

Authors: Maria Okounkova, Will M. Farr, Maximiliano Isi, Leo C. Stein

Abstract: We perform a new test of general relativity (GR) with signals from GWTC-2, the LIGO and Virgo catalog of gravitational wave detections. We search for the presence of amplitude birefringence, in which left versus right circularly polarized modes of gravitational waves are exponentially enhanced and suppressed during propagation. Such an effect is present in various beyond-GR theories but is absent… ▽ More We perform a new test of general relativity (GR) with signals from GWTC-2, the LIGO and Virgo catalog of gravitational wave detections. We search for the presence of amplitude birefringence, in which left versus right circularly polarized modes of gravitational waves are exponentially enhanced and suppressed during propagation. Such an effect is present in various beyond-GR theories but is absent in GR. We constrain the amount of amplitude birefringence consistent with the data through an opacity parameter $κ$, which we bound to be $κ\lesssim 0.74 \textrm{ Gpc}^{-1}$. We then use these theory-agnostic results to constrain Chern-Simons gravity, a beyond-GR theory with motivations in quantum gravity. We bound the canonical Chern-Simons lengthscale to be $\ell_0 \lesssim 1.0 \times 10^3$ km, improving on previous long-distance measurement results by a factor of two. △ Less

Submitted 5 October, 2022; v1 submitted 26 January, 2021; originally announced January 2021.

Comments: Updated to match published version, 8 pages (+6 pages of appendices), 7 figures

Journal ref: Phys. Rev. D 106, 044067 (2022)

Integrability of eccentric, spinning black hole binaries up to second post-Newtonian order

Authors: Sashwat Tanay, Leo C. Stein, José T. Gálvez Ghersi

Abstract: Accurate and efficient modeling of the dynamics of binary black holes (BBHs) is crucial to their detection and parameter estimation through gravitational waves, both with LIGO/Virgo and LISA. General BBH configurations will have misaligned spins and eccentric orbits, eccentricity being particularly relevant at early times. Modeling these systems is both analytically and numerically challenging. Ev… ▽ More Accurate and efficient modeling of the dynamics of binary black holes (BBHs) is crucial to their detection and parameter estimation through gravitational waves, both with LIGO/Virgo and LISA. General BBH configurations will have misaligned spins and eccentric orbits, eccentricity being particularly relevant at early times. Modeling these systems is both analytically and numerically challenging. Even though the 1.5 post-Newtonian (PN) order is Liouville integrable, numerical work has demonstrated chaos at 2PN order, which impedes the existence of an analytic solution. In this article we revisit integrability at both 1.5PN and 2PN orders. At 1.5PN, we construct four (out of five) action integrals. At 2PN, we show that the system is indeed integrable - but in a perturbative sense - by explicitly constructing five mutually-commuting constants of motion. Because of the KAM theorem, this is consistent with the past numerical demonstration of chaos. Our method extends to higher PN orders, opening the door for a fully analytical solution to the generic eccentric, spinning BBH problem. △ Less

Submitted 19 February, 2021; v1 submitted 11 December, 2020; originally announced December 2020.

Comments: 11+2 pages, 2 figures, 1 ancillary Mathematica file; v2: Updated to match version accepted by PRD

Journal ref: Phys. Rev. D 103, 064066 (2021)

Brans-Dicke theory in Bondi-Sachs form: Asymptotically flat solutions, asymptotic symmetries and gravitational-wave memory effects

Authors: Shammi Tahura, David A. Nichols, Alexander Saffer, Leo C. Stein, Kent Yagi

Abstract: Gravitational-wave memory effects are identified by their distinctive effects on families of freely falling observers: after a burst of waves pass by their locations, memory effects can cause lasting relative displacements of the observers. These effects are closely related to the infrared properties of gravity and other massless field theories, including their asymptotic symmetries and conserved… ▽ More Gravitational-wave memory effects are identified by their distinctive effects on families of freely falling observers: after a burst of waves pass by their locations, memory effects can cause lasting relative displacements of the observers. These effects are closely related to the infrared properties of gravity and other massless field theories, including their asymptotic symmetries and conserved quantities. In this paper, we investigate the connection between memory effects, symmetries, and conserved quantities in Brans-Dicke theory. We compute the field equations in Bondi coordinates, and we define a set of boundary conditions that represent asymptotically flat solutions in this context. Next, we derive the asymptotic symmetry group of these spacetimes, and we find that it is the same as the Bondi-Metzner-Sachs group in general relativity. Because there is an additional polarization of gravitational waves in Brans-Dicke theory, we compute the memory effects associated with this extra polarization (the so-called "breathing" mode). This breathing mode produces a uniform expansion (or contraction) of a ring of freely falling observers. After these breathing gravitational waves pass by the observers' locations, there are two additional memory effects that depend on their initial displacements and relative velocities. Neither of these additional memory effects seems to be related to asymptotic symmetries or conserved quantities; rather, they are determined by the properties of the nonradiative region before and after the bursts of the scalar field and the gravitational waves. We discuss the properties of these regions necessary to support nontrivial breathing-mode-type memory effects. △ Less

Submitted 5 July, 2021; v1 submitted 27 July, 2020; originally announced July 2020.

Comments: Appendix added, section 2B modified, minor modifications to the introduction, matches published version

Journal ref: Phys. Rev. D 103, 104026 (2021)

A fixed point for black hole distributions

Authors: José T. Gálvez Ghersi, Leo C. Stein

Abstract: Understanding distributions of black holes is crucial to both astrophysics and quantum gravity. Studying astrophysical population statistics has even been suggested as a channel to constrain black hole formation from the quantum vacuum. Here we propose a Gedankenexperiment to show that the non-linear properties of binary mergers (simulated with accurate surrogate models) generate an attractor in t… ▽ More Understanding distributions of black holes is crucial to both astrophysics and quantum gravity. Studying astrophysical population statistics has even been suggested as a channel to constrain black hole formation from the quantum vacuum. Here we propose a Gedankenexperiment to show that the non-linear properties of binary mergers (simulated with accurate surrogate models) generate an attractor in the space of distributions. Our results show that the joint distribution of spin magnitude and fractional mass loss evolves to a fixed point, converging in a few generations. The features of this fixed point distribution do not depend on the choice of initial distribution. Since a black hole merger is irreversible it produces entropy - possibly the largest source of entropy in the universe. The fixed-point distributions are neither isothermal nor isentropic, and initially thermodynamic states evolve away from thermality. We finally evaluate the specific entropy production rate per merger from initially thermal and non-thermal distributions, which converges to a constant. △ Less

Submitted 26 October, 2020; v1 submitted 22 July, 2020; originally announced July 2020.

Comments: Major revision and results updated. New appendix included. Tested convergence of the joint 2D-3D PDFs and independence of the kick velocity. 9+2 pages. 12 Figures. Comments are welcome

The location of the last stable orbit in Kerr spacetime

Authors: Leo C. Stein, Niels Warburton

Abstract: Black hole spacetimes, like the Kerr spacetime, admit both stable and plunging orbits, separated in parameter space by the separatrix. Determining the location of the separatrix is of fundamental interest in understanding black holes, and is of crucial importance for modeling extreme mass-ratio inspirals. Previous numerical approaches to locating the Kerr separatrix were not always efficient or st… ▽ More Black hole spacetimes, like the Kerr spacetime, admit both stable and plunging orbits, separated in parameter space by the separatrix. Determining the location of the separatrix is of fundamental interest in understanding black holes, and is of crucial importance for modeling extreme mass-ratio inspirals. Previous numerical approaches to locating the Kerr separatrix were not always efficient or stable across all of parameter space. In this paper we show that the Kerr separatrix is the zero set of a single polynomial in parameter space. This gives two main results. First, we thoroughly analyze special cases (extreme Kerr, polar orbits, etc.), finding strict bounds on the limits of roots, and unifying a number of results in the literature. Second, we pose a stable numerical method which is guaranteed to quickly and robustly converge to the separatrix. This new approach is implemented in the Black Hole Perturbation Toolkit, and results in a ~45x speedup over the prior robust approach. △ Less

Submitted 5 March, 2020; v1 submitted 16 December, 2019; originally announced December 2019.

Comments: 11+4 pages, 6 figures, 1 ancillary Mathematica file. v2: Minor edits, matched version published in PRD. Code available at https://bhptoolkit.org/

Journal ref: Phys. Rev. D 101, 064007 (2020)

Numerical relativity simulation of GW150914 beyond general relativity

Authors: Maria Okounkova, Leo C. Stein, Jordan Moxon, Mark A. Scheel, Saul A. Teukolsky

Abstract: We produce the first astrophysically-relevant numerical binary black hole gravitational waveform in a higher-curvature theory of gravity beyond general relativity. We simulate a system with parameters consistent with GW150914, the first LIGO detection, in order-reduced dynamical Chern-Simons gravity, a theory with motivations in string theory and loop quantum gravity. We present results for the le… ▽ More We produce the first astrophysically-relevant numerical binary black hole gravitational waveform in a higher-curvature theory of gravity beyond general relativity. We simulate a system with parameters consistent with GW150914, the first LIGO detection, in order-reduced dynamical Chern-Simons gravity, a theory with motivations in string theory and loop quantum gravity. We present results for the leading-order corrections to the merger and ringdown waveforms, as well as the ringdown quasi-normal mode spectrum. We estimate that such corrections may be discriminated in detections with signal to noise ratio $\gtrsim 180-240$, with the precise value depending on the dimension of the GR waveform family used in data analysis. △ Less

Submitted 7 May, 2020; v1 submitted 6 November, 2019; originally announced November 2019.

Comments: 7 pages + appendices, 8 figures, Updated to match Phys. D. Rev article

Journal ref: Phys. Rev. D 101, 104016 (2020)

qnm: A Python package for calculating Kerr quasinormal modes, separation constants, and spherical-spheroidal mixing coefficients

Authors: Leo C. Stein

Abstract: $\mathtt{qnm}$ is an open-source Python package for computing the Kerr quasinormal mode frequencies, angular separation constants, and spherical-spheroidal mixing coefficients. The $\mathtt{qnm}$ package includes a Leaver solver with the Cook-Zalutskiy spectral approach to the angular sector, and a caching mechanism to avoid repeating calculations. We provide a large cache of low $\ell, m, n$ mode… ▽ More $\mathtt{qnm}$ is an open-source Python package for computing the Kerr quasinormal mode frequencies, angular separation constants, and spherical-spheroidal mixing coefficients. The $\mathtt{qnm}$ package includes a Leaver solver with the Cook-Zalutskiy spectral approach to the angular sector, and a caching mechanism to avoid repeating calculations. We provide a large cache of low $\ell, m, n$ modes, which can be downloaded and installed with a single function call, and interpolated to provide good initial guess for root-polishing at new values of spin. △ Less

Submitted 1 October, 2019; v1 submitted 27 August, 2019; originally announced August 2019.

Comments: 3 pages, 0 figures. Matches version accepted by the Journal of Open Source Software. The source repository: https://github.com/duetosymmetry/qnm , PyPI entry: https://pypi.org/project/qnm/ , documentation: https://qnm.readthedocs.io . Suggestions and pull requests welcome

Journal ref: J. Open Source Softw., 4(42), 1683 (2019)

Numerical binary black hole collisions in dynamical Chern-Simons gravity

Authors: Maria Okounkova, Leo C. Stein, Mark A. Scheel, Saul A. Teukolsky

Abstract: We produce the first numerical relativity binary black hole gravitational waveforms in a higher-curvature theory beyond general relativity. In particular, we study head-on collisions of binary black holes in order-reduced dynamical Chern-Simons gravity. This is a precursor to producing beyond-general-relativity waveforms for inspiraling binary black hole systems that are useful for gravitational w… ▽ More We produce the first numerical relativity binary black hole gravitational waveforms in a higher-curvature theory beyond general relativity. In particular, we study head-on collisions of binary black holes in order-reduced dynamical Chern-Simons gravity. This is a precursor to producing beyond-general-relativity waveforms for inspiraling binary black hole systems that are useful for gravitational wave detection. Head-on collisions are interesting in their own right, however, as they cleanly probe the quasi-normal mode spectrum of the final black hole. We thus compute the leading-order dynamical Chern-Simons modifications to the complex frequencies of the post-merger gravitational radiation. We consider equal-mass systems, with equal spins oriented along the axis of collision, resulting in remnant black holes with spin. We find modifications to the complex frequencies of the quasi-normal mode spectrum that behave as a power law with the spin of the remnant, and that are not degenerate with the frequencies associated with a Kerr black hole of any mass and spin. We discuss these results in the context of testing general relativity with gravitational wave observations. △ Less

Submitted 15 November, 2019; v1 submitted 20 June, 2019; originally announced June 2019.

Comments: 13 pages (+ appendices), 15 figures

Journal ref: Phys. Rev. D 100, 104026 (2019)

Surrogate models for precessing binary black hole simulations with unequal masses

Authors: Vijay Varma, Scott E. Field, Mark A. Scheel, Jonathan Blackman, Davide Gerosa, Leo C. Stein, Lawrence E. Kidder, Harald P. Pfeiffer

Abstract: Only numerical relativity simulations can capture the full complexities of binary black hole mergers. These simulations, however, are prohibitively expensive for direct data analysis applications such as parameter estimation. We present two new fast and accurate surrogate models for the outputs of these simulations: the first model, NRSur7dq4, predicts the gravitational waveform and the second mod… ▽ More Only numerical relativity simulations can capture the full complexities of binary black hole mergers. These simulations, however, are prohibitively expensive for direct data analysis applications such as parameter estimation. We present two new fast and accurate surrogate models for the outputs of these simulations: the first model, NRSur7dq4, predicts the gravitational waveform and the second model, \RemnantModel, predicts the properties of the remnant black hole. These models extend previous 7-dimensional, non-eccentric precessing models to higher mass ratios, and have been trained against 1528 simulations with mass ratios $q\leq4$ and spin magnitudes $χ_1,χ_2 \leq 0.8$, with generic spin directions. The waveform model, NRSur7dq4, which begins about 20 orbits before merger, includes all $\ell \leq 4$ spin-weighted spherical harmonic modes, as well as the precession frame dynamics and spin evolution of the black holes. The final black hole model, \RemnantModel, models the mass, spin, and recoil kick velocity of the remnant black hole. In their training parameter range, both models are shown to be more accurate than existing models by at least an order of magnitude, with errors comparable to the estimated errors in the numerical relativity simulations. We also show that the surrogate models work well even when extrapolated outside their training parameter space range, up to mass ratios $q=6$. △ Less

Submitted 10 October, 2019; v1 submitted 22 May, 2019; originally announced May 2019.

Comments: Matches published version. Models publicly available at https://zenodo.org/record/3455886#.XZ9s1-dKjBI and https://pypi.org/project/surfinBH

Journal ref: Phys. Rev. Research 1, 033015 (2019)

The SXS Collaboration catalog of binary black hole simulations

Authors: Michael Boyle, Daniel Hemberger, Dante A. B. Iozzo, Geoffrey Lovelace, Serguei Ossokine, Harald P. Pfeiffer, Mark A. Scheel, Leo C. Stein, Charles J. Woodford, Aaron B. Zimmerman, Nousha Afshari, Kevin Barkett, Jonathan Blackman, Katerina Chatziioannou, Tony Chu, Nicholas Demos, Nils Deppe, Scott E. Field, Nils L. Fischer, Evan Foley, Heather Fong, Alyssa Garcia, Matthew Giesler, Francois Hebert, Ian Hinder , et al. (15 additional authors not shown)

Abstract: Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of… ▽ More Accurate models of gravitational waves from merging black holes are necessary for detectors to observe as many events as possible while extracting the maximum science. Near the time of merger, the gravitational waves from merging black holes can be computed only using numerical relativity. In this paper, we present a major update of the Simulating eXtreme Spacetimes (SXS) Collaboration catalog of numerical simulations for merging black holes. The catalog contains 2018 distinct configurations (a factor of 11 increase compared to the 2013 SXS catalog), including 1426 spin-precessing configurations, with mass ratios between 1 and 10, and spin magnitudes up to 0.998. The median length of a waveform in the catalog is 39 cycles of the dominant $\ell=m=2$ gravitational-wave mode, with the shortest waveform containing 7.0 cycles and the longest 351.3 cycles. We discuss improvements such as correcting for moving centers of mass and extended coverage of the parameter space. We also present a thorough analysis of numerical errors, finding typical truncation errors corresponding to a waveform mismatch of $\sim 10^{-4}$. The simulations provide remnant masses and spins with uncertainties of 0.03% and 0.1% ($90^{\text{th}}$ percentile), about an order of magnitude better than analytical models for remnant properties. The full catalog is publicly available at https://www.black-holes.org/waveforms . △ Less

Submitted 11 September, 2019; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: 33+18 pages, 13 figures, 4 tables, 2,018 binaries. Catalog metadata in ancillary JSON file. v2: Matches version accepted by CQG. Catalog available at https://www.black-holes.org/waveforms

Journal ref: Class. Quantum Grav. 36, 195006 (2019)

The binary black hole explorer: on-the-fly visualizations of precessing binary black holes

Authors: Vijay Varma, Leo C. Stein, Davide Gerosa

Abstract: Binary black hole mergers are of great interest to the astrophysics community, not least because of their promise to test general relativity in the highly dynamic, strong field regime. Detections of gravitational waves from these sources by LIGO and Virgo have garnered widespread media and public attention. Among these sources, precessing systems (with misaligned black-hole spin/orbital angular mo… ▽ More Binary black hole mergers are of great interest to the astrophysics community, not least because of their promise to test general relativity in the highly dynamic, strong field regime. Detections of gravitational waves from these sources by LIGO and Virgo have garnered widespread media and public attention. Among these sources, precessing systems (with misaligned black-hole spin/orbital angular momentum) are of particular interest because of the rich dynamics they offer. However, these systems are, in turn, more complex compared to nonprecessing systems, making them harder to model or develop intuition about. Visualizations of numerical simulations of precessing systems provide a means to understand and gain insights about these systems. However, since these simulations are very expensive, they can only be performed at a small number of points in parameter space. We present binaryBHexp, a tool that makes use of surrogate models of numerical simulations to generate on-the-fly interactive visualizations of precessing binary black holes. These visualizations can be generated in a few seconds, and at any point in the 7-dimensional parameter space of the underlying surrogate models. With illustrative examples, we demonstrate how this tool can be used to learn about precessing binary black hole systems. △ Less

Submitted 8 April, 2019; v1 submitted 15 November, 2018; originally announced November 2018.

Comments: Matches CQG version. Visualizations available at https://vijayvarma392.github.io/binaryBHexp

High-accuracy mass, spin, and recoil predictions of generic black-hole merger remnants

Authors: Vijay Varma, Davide Gerosa, Leo C. Stein, François Hébert, Hao Zhang

Abstract: We present accurate fits for the remnant properties of generically precessing binary black holes, trained on large banks of numerical-relativity simulations. We use Gaussian process regression to interpolate the remnant mass, spin, and recoil velocity in the 7-dimensional parameter space of precessing black-hole binaries with mass ratios $q\leq2$, and spin magnitudes $χ_1,χ_2\leq0.8$. For precessi… ▽ More We present accurate fits for the remnant properties of generically precessing binary black holes, trained on large banks of numerical-relativity simulations. We use Gaussian process regression to interpolate the remnant mass, spin, and recoil velocity in the 7-dimensional parameter space of precessing black-hole binaries with mass ratios $q\leq2$, and spin magnitudes $χ_1,χ_2\leq0.8$. For precessing systems, our errors in estimating the remnant mass, spin magnitude, and kick magnitude are lower than those of existing fitting formulae by at least an order of magnitude (improvement is also reported in the extrapolated region at high mass ratios and spins). In addition, we also model the remnant spin and kick directions. Being trained directly on precessing simulations, our fits are free from ambiguities regarding the initial frequency at which precessing quantities are defined. We also construct a model for remnant properties of aligned-spin systems with mass ratios $q\leq8$, and spin magnitudes $χ_1,χ_2\leq0.8$. As a byproduct, we also provide error estimates for all fitted quantities, which can be consistently incorporated into current and future gravitational-wave parameter-estimation analyses. Our model(s) are made publicly available through a fast and easy-to-use Python module called surfinBH. △ Less

Submitted 10 January, 2019; v1 submitted 24 September, 2018; originally announced September 2018.

Comments: 6+5 pages. Matches PRL version. Python implementation available at https://pypi.org/project/surfinBH/

Journal ref: Phys. Rev. Lett. 122, 011101 (2019)

Measuring stochastic gravitational-wave energy beyond general relativity

Authors: Maximiliano Isi, Leo C. Stein

Abstract: Gravity theories beyond general relativity (GR) can change the properties of gravitational waves: their polarizations, dispersion, speed, and, importantly, energy content are all heavily theory- dependent. All these corrections can potentially be probed by measuring the stochastic gravitational- wave background. However, most existing treatments of this background beyond GR overlook modifications… ▽ More Gravity theories beyond general relativity (GR) can change the properties of gravitational waves: their polarizations, dispersion, speed, and, importantly, energy content are all heavily theory- dependent. All these corrections can potentially be probed by measuring the stochastic gravitational- wave background. However, most existing treatments of this background beyond GR overlook modifications to the energy carried by gravitational waves, or rely on GR assumptions that are invalid in other theories. This may lead to mistranslation between the observable cross-correlation of detector outputs and gravitational-wave energy density, and thus to errors when deriving observational constraints on theories. In this article, we lay out a generic formalism for stochastic gravitational- wave searches, applicable to a large family of theories beyond GR. We explicitly state the (often tacit) assumptions that go into these searches, evaluating their generic applicability, or lack thereof. Examples of problematic assumptions are: statistical independence of linear polarization amplitudes; which polarizations satisfy equipartition; and which polarizations have well-defined phase velocities. We also show how to correctly infer the value of the stochastic energy density in the context of any given theory. We demonstrate with specific theories in which some of the traditional assumptions break down: Chern-Simons gravity, scalar-tensor theory, and Fierz-Pauli massive gravity. In each theory, we show how to properly include the beyond-GR corrections, and how to interpret observational results. △ Less

Submitted 5 July, 2018; originally announced July 2018.

Comments: 18 pages (plus appendices), 1 figure

Report number: LIGO-P1700234

Journal ref: Phys. Rev. D 98, 104025 (2018)

Black holes, gravitational waves and fundamental physics: a roadmap

Authors: Leor Barack, Vitor Cardoso, Samaya Nissanke, Thomas P. Sotiriou, Abbas Askar, Krzysztof Belczynski, Gianfranco Bertone, Edi Bon, Diego Blas, Richard Brito, Tomasz Bulik, Clare Burrage, Christian T. Byrnes, Chiara Caprini, Masha Chernyakova, Piotr Chrusciel, Monica Colpi, Valeria Ferrari, Daniele Gaggero, Jonathan Gair, Juan Garcia-Bellido, S. F. Hassan, Lavinia Heisenberg, Martin Hendry, Ik Siong Heng , et al. (181 additional authors not shown)

Abstract: The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horiz… ▽ More The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress. △ Less

Submitted 1 February, 2019; v1 submitted 13 June, 2018; originally announced June 2018.

Comments: White Paper for the COST action "Gravitational Waves, Black Holes, and Fundamental Physics", 272 pages, 12 figures; v4: updated references and author list. Overall improvements and corrections. To appear in Classical and Quantum Gravity

Black-hole kicks from numerical-relativity surrogate models

Authors: Davide Gerosa, François Hébert, Leo C. Stein

Abstract: Binary black holes radiate linear momentum in gravitational waves as they merge. Recoils imparted to the black-hole remnant can reach thousands of km/s, thus ejecting black holes from their host galaxies. We exploit recent advances in gravitational waveform modeling to quickly and reliably extract recoils imparted to generic, precessing, black hole binaries. Our procedure uses a numerical-relativi… ▽ More Binary black holes radiate linear momentum in gravitational waves as they merge. Recoils imparted to the black-hole remnant can reach thousands of km/s, thus ejecting black holes from their host galaxies. We exploit recent advances in gravitational waveform modeling to quickly and reliably extract recoils imparted to generic, precessing, black hole binaries. Our procedure uses a numerical-relativity surrogate model to obtain the gravitational waveform given a set of binary parameters, then from this waveform we directly integrate the gravitational-wave linear momentum flux. This entirely bypasses the need of fitting formulae which are typically used to model black-hole recoils in astrophysical contexts. We provide a thorough exploration of the black-hole kick phenomenology in the parameter space, summarizing and extending previous numerical results on the topic. Our extraction procedure is made publicly available as a module for the Python programming language named SURRKICK. Kick evaluations take ~0.1s on a standard off-the-shelf machine, thus making our code ideal to be ported to large-scale astrophysical studies. △ Less

Submitted 28 May, 2018; v1 submitted 12 February, 2018; originally announced February 2018.

Comments: More: https://davidegerosa.com/surrkick - Source: https://github.com/dgerosa/surrkick - pypi: https://pypi.python.org/pypi/surrkick - Published in PRD

Journal ref: Phys. Rev. D 97, 104049 (2018)

Separating metric perturbations in near-horizon extremal Kerr

Authors: Baoyi Chen, Leo C. Stein

Abstract: Linear perturbation theory is a powerful toolkit for studying black hole spacetimes. However, the perturbation equations are hard to solve unless we can use separation of variables. In the Kerr spacetime, metric perturbations do not separate, but curvature perturbations do. The cost of curvature perturbations is a very complicated metric-reconstruction procedure. This procedure can be avoided usin… ▽ More Linear perturbation theory is a powerful toolkit for studying black hole spacetimes. However, the perturbation equations are hard to solve unless we can use separation of variables. In the Kerr spacetime, metric perturbations do not separate, but curvature perturbations do. The cost of curvature perturbations is a very complicated metric-reconstruction procedure. This procedure can be avoided using a symmetry-adapted choice of basis functions in highly symmetric spacetimes, such as near-horizon extremal Kerr. In this paper, we focus on this spacetime, and (i) construct the symmetry-adapted basis functions; (ii) show their orthogonality; and (iii) show that they lead to separation of variables of the scalar, Maxwell, and metric perturbation equations. This separation turns the system of partial differential equations into one of ordinary differential equations over a compact domain, the polar angle. △ Less

Submitted 13 September, 2017; v1 submitted 17 July, 2017; originally announced July 2017.

Comments: 9+9 pages, 0 figures, 6 tables; 3 ancillary Mathematica files. Matches published version

Journal ref: Phys. Rev. D 96, 064017 (2017)

Numerical binary black hole mergers in dynamical Chern-Simons: I. Scalar field

Authors: Maria Okounkova, Leo C. Stein, Mark A. Scheel, Daniel A. Hemberger

Abstract: Testing general relativity in the non-linear, dynamical, strong-field regime of gravity is one of the major goals of gravitational wave astrophysics. Performing precision tests of general relativity (GR) requires numerical inspiral, merger, and ringdown waveforms for binary black hole (BBH) systems in theories beyond GR. Currently, GR and scalar-tensor gravity are the only theories amenable to num… ▽ More Testing general relativity in the non-linear, dynamical, strong-field regime of gravity is one of the major goals of gravitational wave astrophysics. Performing precision tests of general relativity (GR) requires numerical inspiral, merger, and ringdown waveforms for binary black hole (BBH) systems in theories beyond GR. Currently, GR and scalar-tensor gravity are the only theories amenable to numerical simulations. In this article, we present a well-posed perturbation scheme for numerically integrating beyond-GR theories that have a continuous limit to GR. We demonstrate this scheme by simulating BBH mergers in dynamical Chern-Simons gravity (dCS), to linear order in the perturbation parameter. We present mode waveforms and energy fluxes of the dCS pseudoscalar field from our numerical simulations. We find good agreement with analytic predictions at early times, including the absence of pseudoscalar dipole radiation. We discover new phenomenology only accessible through numerics: a burst of dipole radiation during merger. We also quantify the self-consistency of the perturbation scheme. Finally, we estimate bounds that GR-consistent LIGO detections could place on the new dCS length scale, approximately $\ell \lesssim \mathcal{O}(10)~\mathrm{km}$. △ Less

Submitted 17 August, 2017; v1 submitted 22 May, 2017; originally announced May 2017.

Comments: 14+4 pages, 8 figures, 1 table; visualization available at http://www.youtube.com/watch?v=WQH-1b_XUM4 . Matches published version

Journal ref: Phys. Rev. D 96, 044020 (2017)

Black Hole Based Tests of General Relativity

Authors: Kent Yagi, Leo C. Stein

Abstract: General relativity has passed all solar system experiments and neutron star based tests, such as binary pulsar observations, with flying colors. A more exotic arena for testing general relativity is in systems that contain one or more black holes. Black holes are the most compact objects in the universe, providing probes of the strongest-possible gravitational fields. We are motivated to study str… ▽ More General relativity has passed all solar system experiments and neutron star based tests, such as binary pulsar observations, with flying colors. A more exotic arena for testing general relativity is in systems that contain one or more black holes. Black holes are the most compact objects in the universe, providing probes of the strongest-possible gravitational fields. We are motivated to study strong-field gravity since many theories give large deviations from general relativity only at large field strengths, while recovering the weak-field behavior. In this article, we review how one can probe general relativity and various alternative theories of gravity by using electromagnetic waves from a black hole with an accretion disk, and gravitational waves from black hole binaries. We first review model-independent ways of testing gravity with electromagnetic/gravitational waves from a black hole system. We then focus on selected examples of theories that extend general relativity in rather simple ways. Some important characteristics of general relativity include (but are not limited to) (i) only tensor gravitational degrees of freedom, (ii) the graviton is massless, (iii) no quadratic or higher curvatures in the action, and (iv) the theory is 4 dimensional. Altering a characteristic leads to a different extension of general relativity: (i) scalar-tensor theories, (ii) massive gravity theories, (iii) quadratic gravity, and (iv) theories with large extra dimensions. Within each theory, we describe black hole solutions, their properties, and current and projected constraints on each theory using black hole-based tests of gravity. We close this review by listing some of the open problems in model-independent tests and within each specific theory. △ Less

Submitted 18 February, 2016; v1 submitted 7 February, 2016; originally announced February 2016.

Comments: 95 pages, 11 figures; a nomenclature for the braneworld BH conjecture corrected, references added and updated; Invited review article to the Focus Issue on "Black holes and fundamental fields" published in Classical and Quantum Gravity

Journal ref: Class. Quantum Grav. 33 054001 (2016)

Extremal Black Holes in Dynamical Chern-Simons Gravity

Authors: Robert McNees, Leo C. Stein, Nicolás Yunes

Abstract: Rapidly rotating black hole solutions in theories beyond general relativity play a key role in experimental gravity, as they allow us to compute observables in extreme spacetimes that deviate from the predictions of general relativity. Such solutions are often difficult to find in beyond-general-relativity theories due to the inclusion of additional fields that couple to the metric non-linearly an… ▽ More Rapidly rotating black hole solutions in theories beyond general relativity play a key role in experimental gravity, as they allow us to compute observables in extreme spacetimes that deviate from the predictions of general relativity. Such solutions are often difficult to find in beyond-general-relativity theories due to the inclusion of additional fields that couple to the metric non-linearly and non-minimally. In this paper, we consider rotating black hole solutions in one such theory, dynamical Chern-Simons gravity, where the Einstein-Hilbert action is modified by the introduction of a dynamical scalar field that couples to the metric through the Pontryagin density. We treat dynamical Chern-Simons gravity as an effective field theory and work in the decoupling limit, where corrections are treated as small perturbations from general relativity. We perturb about the maximally-rotating Kerr solution, the so-called extremal limit, and develop mathematical insight into the analysis techniques needed to construct solutions for generic spin. First we find closed-form, analytic expressions for the extremal scalar field, and then determine the trace of the metric perturbation, giving both in terms of Legendre decompositions. Retaining only the first three and four modes in the Legendre representation of the scalar field and the trace, respectively, suffices to ensure a fidelity of over 99% relative to full numerical solutions. The leading-order mode in the Legendre expansion of the trace of the metric perturbation contains a logarithmic divergence at the extremal Kerr horizon, which is likely to be unimportant as it occurs inside the perturbed dynamical Chern-Simons horizon. The techniques employed here should enable the construction of analytic, closed-form expressions for the scalar field and metric perturbations on a background with arbitrary rotation. △ Less

Submitted 26 November, 2016; v1 submitted 16 December, 2015; originally announced December 2015.

Comments: 25+9 pages (single column), 10 figures, 1 table; matches published version

Journal ref: Class. Quantum Grav. 33, 235013 (2016)

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

Authors: Kent Yagi, Leo C. Stein, Nicolas Yunes

Abstract: 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 t… ▽ More 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. △ Less

Submitted 29 December, 2015; v1 submitted 7 October, 2015; originally announced October 2015.

Comments: 29 pages, 6 figures; references updated, matches PRD version

Journal ref: Phys. Rev. D 93 (2016) 2, 024010

Slimplectic Integrators: Variational Integrators for General Nonconservative Systems

Authors: David Tsang, Chad R. Galley, Leo C. Stein, Alec Turner

Abstract: Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative interactions. In this Letter, we develop the "slimplectic" integrator, a new type of numerical integrator that shares many of the benefits of traditional symple… ▽ More Symplectic integrators are widely used for long-term integration of conservative astrophysical problems due to their ability to preserve the constants of motion; however, they cannot in general be applied in the presence of nonconservative interactions. In this Letter, we develop the "slimplectic" integrator, a new type of numerical integrator that shares many of the benefits of traditional symplectic integrators yet is applicable to general nonconservative systems. We utilize a fixed time-step variational integrator formalism applied to the principle of stationary nonconservative action developed in Galley, 2013; Galley, Tsang & Stein, 2014. As a result, the generalized momenta and energy (Noether current) evolutions are well-tracked. We discuss several example systems, including damped harmonic oscillators, Poynting-Robertson drag, and gravitational radiation reaction, by utilizing our new publicly available code to demonstrate the slimplectic integrator algorithm. Slimplectic integrators are well-suited for integrations of systems where nonconservative effects play an important role in the long-term dynamical evolution. As such they are particularly appropriate for cosmological or celestial N-body dynamics problems where nonconservative interactions, e.g. gas interactions or dissipative tides, can play an important role. △ Less

Submitted 22 July, 2015; v1 submitted 28 June, 2015; originally announced June 2015.

Comments: 6 pages, 5 Figures; Accepted to ApJL; code repository at http://github.com/davtsang/slimplectic (typo fixed in repository URL)

Journal ref: 2015 ApJ 809 L9

Testing General Relativity with Present and Future Astrophysical Observations

Authors: Emanuele Berti, Enrico Barausse, Vitor Cardoso, Leonardo Gualtieri, Paolo Pani, Ulrich Sperhake, Leo C. Stein, Norbert Wex, Kent Yagi, Tessa Baker, C. P. Burgess, Flávio S. Coelho, Daniela Doneva, Antonio De Felice, Pedro G. Ferreira, Paulo C. C. Freire, James Healy, Carlos Herdeiro, Michael Horbatsch, Burkhard Kleihaus, Antoine Klein, Kostas Kokkotas, Jutta Kunz, Pablo Laguna, Ryan N. Lang , et al. (28 additional authors not shown)

Abstract: One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The… ▽ More One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime. △ Less

Submitted 1 December, 2015; v1 submitted 28 January, 2015; originally announced January 2015.

Comments: 188 pages, 46 figures, 6 tables, 903 references. Matches version published in Classical and Quantum Gravity. Supplementary data files available at http://www.phy.olemiss.edu/~berti/research/ and http://centra.tecnico.ulisboa.pt/network/grit/files/

Journal ref: Class. Quantum Grav. 32, 243001 (2015)

The principle of stationary nonconservative action for classical mechanics and field theories

Authors: Chad R. Galley, David Tsang, Leo C. Stein

Abstract: We further develop a recently introduced variational principle of stationary action for problems in nonconservative classical mechanics and extend it to classical field theories. The variational calculus used is consistent with an initial value formulation of physical problems and allows for time-irreversible processes, such as dissipation, to be included at the level of the action. In this formal… ▽ More We further develop a recently introduced variational principle of stationary action for problems in nonconservative classical mechanics and extend it to classical field theories. The variational calculus used is consistent with an initial value formulation of physical problems and allows for time-irreversible processes, such as dissipation, to be included at the level of the action. In this formalism, the equations of motion are generated by extremizing a nonconservative action $\mathcal{S}$, which is a functional of a doubled set of degrees of freedom. The corresponding nonconservative Lagrangian contains a potential $K$ which generates nonconservative forces and interactions. Such a nonconservative potential can arise in several ways, including from an open system interacting with inaccessible degrees of freedom or from integrating out or coarse-graining a subset of variables in closed systems. We generalize Noether's theorem to show how Noether currents are modified and no longer conserved when $K$ is non-vanishing. Consequently, the nonconservative aspects of a physical system are derived solely from $K$. We show how to use the formalism with examples of nonconservative actions for discrete systems including forced damped harmonic oscillators, radiation reaction on an accelerated charge, and RLC circuits. We present examples for nonconservative classical field theories. Our approach naturally allows for irreversible thermodynamic processes to be included in an unconstrained variational principle. We present the nonconservative action for a Navier-Stokes fluid including the effects of viscous dissipation and heat diffusion, as well as an action that generates the Maxwell model for viscoelastic materials, which can be easily generalized to more realistic rheological models. We show that the nonconservative action can be derived as the classical limit of a more complete quantum theory. △ Less

Submitted 9 December, 2014; originally announced December 2014.

Comments: 42 pages, 4 figures, 1 table

Rapidly rotating black holes in dynamical Chern-Simons gravity: Decoupling limit solutions and breakdown

Authors: Leo C. Stein

Abstract: Rapidly rotating black holes are a prime arena for understanding corrections to Einstein's theory of general relativity (GR). We construct solutions for rapidly rotating black holes in dynamical Chern-Simons (dCS) gravity, a useful and motivated example of a post-GR correction. We treat dCS as an effective theory and thus work in the decoupling limit, where we apply a perturbation scheme using the… ▽ More Rapidly rotating black holes are a prime arena for understanding corrections to Einstein's theory of general relativity (GR). We construct solutions for rapidly rotating black holes in dynamical Chern-Simons (dCS) gravity, a useful and motivated example of a post-GR correction. We treat dCS as an effective theory and thus work in the decoupling limit, where we apply a perturbation scheme using the Kerr metric as the background solution. Using the solutions to the scalar field and the trace of the metric perturbation, we determine the regime of validity of our perturbative approach. We find that the maximal spin limit may be divergent, and the decoupling limit is strongly restricted for rapid rotation. Rapidly-rotating stellar-mass BHs can potentially be used to place strong bounds on the coupling parameter $\ell$ of dCS. In order for the black hole observed in GRO J1655-40 to be within the decoupling limit we need $\ell \lesssim 22$ km, a value 7 orders of magnitude smaller than present Solar System bounds on dynamical Chern-Simons gravity. △ Less

Submitted 22 August, 2014; v1 submitted 9 July, 2014; originally announced July 2014.

Comments: 8 pages, 4 figures; matches version appearing in PRD

Journal ref: Phys. Rev. D 90, 044061 (2014)

Why I-Love-Q

Authors: Kent Yagi, Leo C. Stein, George Pappas, Nicolas Yunes, Theocharis A. Apostolatos

Abstract: 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,… ▽ More 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. △ Less

Submitted 8 September, 2014; v1 submitted 29 June, 2014; originally announced June 2014.

Comments: 23 pages, 22 figures; matches version appearing in PRD

Journal ref: Phys. Rev. D 90, 063010 (2014)

Three-Hair Relations for Rotating Stars: Nonrelativistic Limit

Authors: Leo C. Stein, Kent Yagi, Nicolas Yunes

Abstract: 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 completel… ▽ More 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. △ Less

Submitted 17 May, 2014; v1 submitted 16 December, 2013; originally announced December 2013.

Comments: 5 pages, 3 figures. Updated to match version published in ApJ

Journal ref: 2014 ApJ 788 15

Parametrizing and constraining scalar corrections to general relativity

Authors: Leo C. Stein, Kent Yagi

Abstract: We parametrize a large class of corrections to general relativity which include a long-ranged gravitational scalar field as a dynamical degree of freedom in two ways: parametrizing the structure of the correction to the action, and parametrizing the scalar hair (multipole structure) that compact objects and black holes attain. The presence of this scalar hair violates the no-hair theorems present… ▽ More We parametrize a large class of corrections to general relativity which include a long-ranged gravitational scalar field as a dynamical degree of freedom in two ways: parametrizing the structure of the correction to the action, and parametrizing the scalar hair (multipole structure) that compact objects and black holes attain. The presence of this scalar hair violates the no-hair theorems present in general relativity, which leads to several important effects. The effects we consider are (i) the interaction between an isolated body and an external scalar field, (ii) the scalar multipole-multipole interaction between two bodies in a compact binary, (iii) the additional pericenter precession of a binary, (iv) the scalar radiation from a binary, and (v) the modification to the gravitational wave phase from a binary. We apply this framework to example theories including Einstein-dilaton-Gauss-Bonnet gravity and dynamical Chern-Simons gravity, and estimate the size of the effects. Finally, we estimate the bounds that can be placed on parameters of the theories from the precession of pulsar binaries and from gravitational waves. △ Less

Submitted 24 February, 2014; v1 submitted 24 October, 2013; originally announced October 2013.

Comments: 18 pages, 3 figures, 1 table; matches version appearing in PRD

Journal ref: Phys. Rev. D 89, 044026 (2014)

Isolated and Binary Neutron Stars in Dynamical Chern-Simons Gravity

Authors: Kent Yagi, Leo C. Stein, Nicolas Yunes, Takahiro Tanaka

Abstract: 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… ▽ More 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. △ Less

Submitted 12 April, 2016; v1 submitted 7 February, 2013; originally announced February 2013.

Comments: 28 pages, 7 figures, 3 tables; erratum added, main results unchanged

Journal ref: Phys. Rev. D 87, 084058 (2013)

Search for gravitational waves associated with gamma-ray bursts during LIGO science run 6 and Virgo science runs 2 and 3

Authors: The LIGO Scientific Collaboration, Virgo Collaboration, J. Abadie, B. P. Abbott, R. Abbott, T. D. Abbott, M. Abernathy, T. Accadia, F. Acernese, C. Adams, R. Adhikari, C. Affeldt, M. Agathos, K. Agatsuma, P. Ajith, B. Allen, E. Amador Ceron, D. Amariutei, S. B. Anderson, W. G. Anderson, K. Arai, M. A. Arain, M. C. Araya, S. M. Aston, P. Astone , et al. (785 additional authors not shown)

Abstract: We present the results of a search for gravitational waves associated with 154 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments in 2009-2010, during the sixth LIGO science run and the second and third Virgo science runs. We perform two distinct searches: a modeled search for coalescences of either two neutron stars or a neutron star and black hole; and a search f… ▽ More We present the results of a search for gravitational waves associated with 154 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments in 2009-2010, during the sixth LIGO science run and the second and third Virgo science runs. We perform two distinct searches: a modeled search for coalescences of either two neutron stars or a neutron star and black hole; and a search for generic, unmodeled gravitational-wave bursts. We find no evidence for gravitational-wave counterparts, either with any individual GRB in this sample or with the population as a whole. For all GRBs we place lower bounds on the distance to the progenitor, under the optimistic assumption of a gravitational-wave emission energy of 10^-2 M c^2 at 150 Hz, with a median limit of 17 Mpc. For short hard GRBs we place exclusion distances on binary neutron star and neutron star-black hole progenitors, using astrophysically motivated priors on the source parameters, with median values of 16 Mpc and 28 Mpc respectively. These distance limits, while significantly larger than for a search that is not aided by GRB satellite observations, are not large enough to expect a coincidence with a GRB. However, projecting these exclusions to the sensitivities of Advanced LIGO and Virgo, which should begin operation in 2015, we find that the detection of gravitational waves associated with GRBs will become quite possible. △ Less

Submitted 24 September, 2012; v1 submitted 10 May, 2012; originally announced May 2012.

Comments: Accepted by ApJ. 19 pages, 8 figures, 2 tables, science summary page at http://www.ligo.org/science/Publication-S6GRB/index.php . Public access area to figures, tables at https://dcc.ligo.org/cgi-bin/DocDB/ShowDocument?docid=p1000121

Report number: LIGO-P1000121

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

Authors: Kent Yagi, Leo C. Stein, Nicolas Yunes, Takahiro Tanaka

Abstract: 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… ▽ More 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. △ Less

Submitted 2 December, 2015; v1 submitted 26 October, 2011; originally announced October 2011.

Comments: 26 pages, 3 figures, 2 tables (matches version published in PRD); v3 prepended an erratum

Journal ref: Phys. Rev. D 85, 064022 (2012); Phys. Rev. D 93, 029902 (2016)

Directional limits on persistent gravitational waves using LIGO S5 science data

Authors: B. P. Abbott, R. Abbott, M. Abernathy, T. Accadia, F. Acernese, C. Adams, R. Adhikari, P. Ajith, B. Allen, G. S. Allen, E. Amador Ceron, R. S. Amin, S. B. Anderson, W. G. Anderson, F. Antonucci, M. A. Arain, M. C. Araya, M. Aronsson, K. G. Arun, Y. Aso, S. M. Aston, P. Astone, D. Atkinson, P. Aufmuth, C. Aulbert , et al. (689 additional authors not shown)

Abstract: The gravitational-wave (GW) sky may include nearby pointlike sources as well as astrophysical and cosmological stochastic backgrounds. Since the relative strength and angular distribution of the many possible sources of GWs are not well constrained, searches for GW signals must be performed in a model-independent way. To that end we perform two directional searches for persistent GWs using data fr… ▽ More The gravitational-wave (GW) sky may include nearby pointlike sources as well as astrophysical and cosmological stochastic backgrounds. Since the relative strength and angular distribution of the many possible sources of GWs are not well constrained, searches for GW signals must be performed in a model-independent way. To that end we perform two directional searches for persistent GWs using data from the LIGO S5 science run: one optimized for pointlike sources and one for arbitrary extended sources. The latter result is the first of its kind. Finding no evidence to support the detection of GWs, we present 90% confidence level (CL) upper-limit maps of GW strain power with typical values between 2-20x10^-50 strain^2 Hz^-1 and 5-35x10^-49 strain^2 Hz^-1 sr^-1 for pointlike and extended sources respectively. The limits on pointlike sources constitute a factor of 30 improvement over the previous best limits. We also set 90% CL limits on the narrow-band root-mean-square GW strain from interesting targets including Sco X-1, SN1987A and the Galactic Center as low as ~7x10^-25 in the most sensitive frequency range near 160 Hz. These limits are the most constraining to date and constitute a factor of 5 improvement over the previous best limits. △ Less

Submitted 9 September, 2011; v1 submitted 8 September, 2011; originally announced September 2011.

Comments: 10 pages, 4 figures

Journal ref: Phys. Rev. Lett. 107, 271102 (2011)

Non-Spinning Black Holes in Alternative Theories of Gravity

Authors: Nicolas Yunes, Leo C. Stein

Abstract: 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.… ▽ More 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. △ Less

Submitted 8 March, 2011; v1 submitted 14 January, 2011; originally announced January 2011.

Comments: 7 pages, submitted to PRD

Journal ref: Phys.Rev.D83:104002,2011

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

Authors: Leo C. Stein, Nicolas Yunes

Abstract: 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… ▽ More 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. △ Less

Submitted 21 January, 2011; v1 submitted 14 December, 2010; originally announced December 2010.

Comments: 19 pages. v3: Clarified averaging properties and gauge choice. Version accepted to PRD

Journal ref: Phys.Rev.D83:064038,2011

A search for gravitational waves associated with the August 2006 timing glitch of the Vela pulsar

Authors: The LIGO Scientific Collaboration, J. Abadie, B. P. Abbott, R. Abbott, R. Adhikari, P. Ajith, B. Allen, G. Allen, E. Amador Ceron, R. S. Amin, S. B. Anderson, W. G. Anderson, M. A. Arain, M. Araya, Y. Aso, S. Aston, P. Aufmuth, C. Aulbert, S. Babak, P. Baker, S. Ballmer, D. Barker, B. Barr, P. Barriga, L. Barsotti , et al. (477 additional authors not shown)

Abstract: The physical mechanisms responsible for pulsar timing glitches are thought to excite quasi-normal mode oscillations in their parent neutron star that couple to gravitational wave emission. In August 2006, a timing glitch was observed in the radio emission of PSR B0833-45, the Vela pulsar. At the time of the glitch, the two co-located Hanford gravitational wave detectors of the Laser Interferometer… ▽ More The physical mechanisms responsible for pulsar timing glitches are thought to excite quasi-normal mode oscillations in their parent neutron star that couple to gravitational wave emission. In August 2006, a timing glitch was observed in the radio emission of PSR B0833-45, the Vela pulsar. At the time of the glitch, the two co-located Hanford gravitational wave detectors of the Laser Interferometer Gravitational-wave observatory (LIGO) were operational and taking data as part of the fifth LIGO science run (S5). We present the first direct search for the gravitational wave emission associated with oscillations of the fundamental quadrupole mode excited by a pulsar timing glitch. No gravitational wave detection candidate was found. We place Bayesian 90% confidence upper limits of 6.3e-21 to 1.4e-20 on the peak intrinsic strain amplitude of gravitational wave ring-down signals, depending on which spherical harmonic mode is excited. The corresponding range of energy upper limits is 5.0e44 to 1.3e45 erg. △ Less

Submitted 23 November, 2010; v1 submitted 5 November, 2010; originally announced November 2010.

Comments: 13 pages, 5 figures

Journal ref: Phys.Rev.D83:042001,2011; Publisher-note D83:069902,2011; Phys.Rev.D83:069902,2011

First search for gravitational waves from the youngest known neutron star

Authors: LIGO Scientific Collaboration, J. Abadie, B. P. Abbott, R. Abbott, M. Abernathy, C. Adams, R. Adhikari, P. Ajith, B. Allen, G. Allen, E. Amador Ceron, R. S. Amin, S. B. Anderson, W. G. Anderson, M. A. Arain, M. Araya, M. Aronsson, Y. Aso, S. Aston, D. E. Atkinson, P. Aufmuth, C. Aulbert, S. Babak, P. Baker, S. Ballmer , et al. (515 additional authors not shown)

Abstract: We present a search for periodic gravitational waves from the neutron star in the supernova remnant Cassiopeia A. The search coherently analyzes data in a 12-day interval taken from the fifth science run of the Laser Interferometer Gravitational-Wave Observatory. It searches gravitational wave frequencies from 100 to 300 Hz, and covers a wide range of first and second frequency derivatives appropr… ▽ More We present a search for periodic gravitational waves from the neutron star in the supernova remnant Cassiopeia A. The search coherently analyzes data in a 12-day interval taken from the fifth science run of the Laser Interferometer Gravitational-Wave Observatory. It searches gravitational wave frequencies from 100 to 300 Hz, and covers a wide range of first and second frequency derivatives appropriate for the age of the remnant and for different spin-down mechanisms. No gravitational wave signal was detected. Within the range of search frequencies, we set 95% confidence upper limits of 0.7--1.2e-24 on the intrinsic gravitational wave strain, 0.4--4e-4 on the equatorial ellipticity of the neutron star, and 0.005--0.14 on the amplitude of r-mode oscillations of the neutron star. These direct upper limits beat indirect limits derived from energy conservation and enter the range of theoretical predictions involving crystalline exotic matter or runaway r-modes. This is the first gravitational wave search to present upper limits on r-modes. △ Less

Submitted 9 September, 2010; v1 submitted 13 June, 2010; originally announced June 2010.

Comments: 26 pages, 5 figures

Report number: LIGO-P1000028-v7

Journal ref: Astrophys.J.722:1504-1513,2010

Predictions for the Rates of Compact Binary Coalescences Observable by Ground-based Gravitational-wave Detectors

Authors: LIGO Scientific Collaboration, Virgo Collaboration, J. Abadie, B. P. Abbott, R. Abbott, M Abernathy, T. Accadia, F. Acernese, C. Adams, R. Adhikari, P. Ajith, B. Allen, G. Allen, E. Amador Ceron, R. S. Amin, S. B. Anderson, W. G. Anderson, F. Antonucci, S. Aoudia, M. A. Arain, M. Araya, M. Aronsson, K. G. Arun, Y. Aso, S. Aston , et al. (687 additional authors not shown)

Abstract: We present an up-to-date, comprehensive summary of the rates for all types of compact binary coalescence sources detectable by the Initial and Advanced versions of the ground-based gravitational-wave detectors LIGO and Virgo. Astrophysical estimates for compact-binary coalescence rates depend on a number of assumptions and unknown model parameters, and are still uncertain. The most confident amo… ▽ More We present an up-to-date, comprehensive summary of the rates for all types of compact binary coalescence sources detectable by the Initial and Advanced versions of the ground-based gravitational-wave detectors LIGO and Virgo. Astrophysical estimates for compact-binary coalescence rates depend on a number of assumptions and unknown model parameters, and are still uncertain. The most confident among these estimates are the rate predictions for coalescing binary neutron stars which are based on extrapolations from observed binary pulsars in our Galaxy. These yield a likely coalescence rate of 100 per Myr per Milky Way Equivalent Galaxy (MWEG), although the rate could plausibly range from 1 per Myr per MWEG to 1000 per Myr per MWEG. We convert coalescence rates into detection rates based on data from the LIGO S5 and Virgo VSR2 science runs and projected sensitivities for our Advanced detectors. Using the detector sensitivities derived from these data, we find a likely detection rate of 0.02 per year for Initial LIGO-Virgo interferometers, with a plausible range between 0.0002 and 0.2 per year. The likely binary neutron-star detection rate for the Advanced LIGO-Virgo network increases to 40 events per year, with a range between 0.4 and 400 per year. △ Less

Submitted 26 March, 2010; v1 submitted 12 March, 2010; originally announced March 2010.

Journal ref: Class.Quant.Grav.27:173001,2010

Search for gravitational-wave inspiral signals associated with short Gamma-Ray Bursts during LIGO's fifth and Virgo's first science run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, J. Abadie, B. P. Abbott, R. Abbott, T. Accadia, F. Acernese, R. Adhikari, P. Ajith, B. Allen, G. Allen, E. Amador Ceron, R. S. Amin, S. B. Anderson, W. G. Anderson, F. Antonucci, S. Aoudia, M. A. Arain, M. Araya, K. G. Arun, Y. Aso, S. Aston, P. Astone, P. Aufmuth, C. Aulbert , et al. (643 additional authors not shown)

Abstract: Progenitor scenarios for short gamma-ray bursts (short GRBs) include coalescenses of two neutron stars or a neutron star and black hole, which would necessarily be accompanied by the emission of strong gravitational waves. We present a search for these known gravitational-wave signatures in temporal and directional coincidence with 22 GRBs that had sufficient gravitational-wave data available in… ▽ More Progenitor scenarios for short gamma-ray bursts (short GRBs) include coalescenses of two neutron stars or a neutron star and black hole, which would necessarily be accompanied by the emission of strong gravitational waves. We present a search for these known gravitational-wave signatures in temporal and directional coincidence with 22 GRBs that had sufficient gravitational-wave data available in multiple instruments during LIGO's fifth science run, S5, and Virgo's first science run, VSR1. We find no statistically significant gravitational-wave candidates within a [-5, +1) s window around the trigger time of any GRB. Using the Wilcoxon-Mann-Whitney U test, we find no evidence for an excess of weak gravitational-wave signals in our sample of GRBs. We exclude neutron star-black hole progenitors to a median 90% CL exclusion distance of 6.7 Mpc. △ Less

Submitted 3 March, 2010; v1 submitted 4 January, 2010; originally announced January 2010.

Comments: 8 pages, 2 figures, Accepted by ApJ

Journal ref: Astrophys.J.715:1453-1461,2010

Searches for gravitational waves from known pulsars with S5 LIGO data

Authors: The LIGO Scientific Collaboration, The Virgo Collaboration, B. P. Abbott, R. Abbott, F. Acernese, R. Adhikari, P. Ajith, B. Allen, G. Allen, M. Alshourbagy, R. S. Amin, S. B. Anderson, W. G. Anderson, F. Antonucci, S. Aoudia, M. A. Arain, M. Araya, H. Armandula, P. Armor, K. G. Arun, Y. Aso, S. Aston, P. Astone, P. Aufmuth, C. Aulbert , et al. (656 additional authors not shown)

Abstract: We present a search for gravitational waves from 116 known millisecond and young pulsars using data from the fifth science run of the LIGO detectors. For this search ephemerides overlapping the run period were obtained for all pulsars using radio and X-ray observations. We demonstrate an updated search method that allows for small uncertainties in the pulsar phase parameters to be included in th… ▽ More We present a search for gravitational waves from 116 known millisecond and young pulsars using data from the fifth science run of the LIGO detectors. For this search ephemerides overlapping the run period were obtained for all pulsars using radio and X-ray observations. We demonstrate an updated search method that allows for small uncertainties in the pulsar phase parameters to be included in the search. We report no signal detection from any of the targets and therefore interpret our results as upper limits on the gravitational wave signal strength. The most interesting limits are those for young pulsars. We present updated limits on gravitational radiation from the Crab pulsar, where the measured limit is now a factor of seven below the spin-down limit. This limits the power radiated via gravitational waves to be less than ~2% of the available spin-down power. For the X-ray pulsar J0537-6910 we reach the spin-down limit under the assumption that any gravitational wave signal from it stays phase locked to the X-ray pulses over timing glitches, and for pulsars J1913+1011 and J1952+3252 we are only a factor of a few above the spin-down limit. Of the recycled millisecond pulsars several of the measured upper limits are only about an order of magnitude above their spin-down limits. For these our best (lowest) upper limit on gravitational wave amplitude is 2.3x10^-26 for J1603-7202 and our best (lowest) limit on the inferred pulsar ellipticity is 7.0x10^-8 for J2124-3358. △ Less

Submitted 26 February, 2010; v1 submitted 19 September, 2009; originally announced September 2009.

Comments: 39 pages, 5 figures, Accepted in The Astrophysical Journal

Report number: LIGO-P080112-v6

Journal ref: Astrophys.J.713:671-685,2010

Search for gravitational-wave bursts associated with gamma-ray bursts using data from LIGO Science Run 5 and Virgo Science Run 1

Authors: LIGO Scientific Collaboration, Virgo Collaboration, B. P. Abbott, R. Abbott, F. Acernese, R. Adhikari, P. Ajith, B. Allen, G. Allen, M. Alshourbagy, R. S. Amin, S. B. Anderson, W. G. Anderson, F. Antonucci, S. Aoudia, M. A. Arain, M. Araya, H. Armandula, P. Armor, K. G. Arun, Y. Aso, S. Aston, P. Astone, P. Aufmuth, C. Aulbert , et al. (643 additional authors not shown)

Abstract: We present the results of a search for gravitational-wave bursts associated with 137 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments during the fifth LIGO science run and first Virgo science run. The data used in this analysis were collected from 2005 November 4 to 2007 October 1, and most of the GRB triggers were from the Swift satellite. The search uses a co… ▽ More We present the results of a search for gravitational-wave bursts associated with 137 gamma-ray bursts (GRBs) that were detected by satellite-based gamma-ray experiments during the fifth LIGO science run and first Virgo science run. The data used in this analysis were collected from 2005 November 4 to 2007 October 1, and most of the GRB triggers were from the Swift satellite. The search uses a coherent network analysis method that takes into account the different locations and orientations of the interferometers at the three LIGO-Virgo sites. We find no evidence for gravitational-wave burst signals associated with this sample of GRBs. Using simulated short-duration (<1 s) waveforms, we set upper limits on the amplitude of gravitational waves associated with each GRB. We also place lower bounds on the distance to each GRB under the assumption of a fixed energy emission in gravitational waves, with typical limits of D ~ 15 Mpc (E_GW^iso / 0.01 M_o c^2)^1/2 for emission at frequencies around 150 Hz, where the LIGO-Virgo detector network has best sensitivity. We present astrophysical interpretations and implications of these results, and prospects for corresponding searches during future LIGO-Virgo runs. △ Less

Submitted 7 April, 2010; v1 submitted 26 August, 2009; originally announced August 2009.

Comments: 16 pages, 3 figures. Updated references. To appear in ApJ.

Report number: LIGO-P0900023-v16

Journal ref: Astrophys.J.715:1438-1452, 2010

Beating the spin-down limit on gravitational wave emission from the Crab pulsar

Authors: The LIGO Scientific Collaboration, B. Abbott, R. Abbott, R. Adhikari, P. Ajith, B. Allen, G. Allen, R. Amin, S. B. Anderson, W. G. Anderson, M. A. Arain, M. Araya, H. Armandula, P. Armor, Y. Aso, S. Aston, P. Aufmuth, C. Aulbert, S. Babak, S. Ballmer, H. Bantilan, B. C. Barish, C. Barker, D. Barker, B. Barr , et al. (419 additional authors not shown)

Abstract: We present direct upper limits on gravitational wave emission from the Crab pulsar using data from the first nine months of the fifth science run of the Laser Interferometer Gravitational-wave Observatory (LIGO). These limits are based on two searches. In the first we assume that the gravitational wave emission follows the observed radio timing, giving an upper limit on gravitational wave emissi… ▽ More We present direct upper limits on gravitational wave emission from the Crab pulsar using data from the first nine months of the fifth science run of the Laser Interferometer Gravitational-wave Observatory (LIGO). These limits are based on two searches. In the first we assume that the gravitational wave emission follows the observed radio timing, giving an upper limit on gravitational wave emission that beats indirect limits inferred from the spin-down and braking index of the pulsar and the energetics of the nebula. In the second we allow for a small mismatch between the gravitational and radio signal frequencies and interpret our results in the context of two possible gravitational wave emission mechanisms. △ Less

Submitted 22 July, 2008; v1 submitted 30 May, 2008; originally announced May 2008.

Comments: Accepted for Ap. J. Lett. Minor changes in results due to calibration correction

Report number: LIGO-P070118-00-Z

Journal ref: Astrophys.J.683:L45-L50,2008; Erratum-ibid 706:L203-L204,2009