results for au:Johnson_McDaniel_N in:gr-qc

- Nov 21 2017 gr-qc astro-ph.HE arXiv:1711.06843v1We present the results of a search for long-duration gravitational wave transients in the data of the LIGO Hanford and LIGO Livingston second generation detectors between September 2015 and January 2016, with a total observational time of 49 days. The search targets gravitational wave transients of \unit[10 -- 500]s duration in a frequency band of \unit[24 -- 2048]Hz, with minimal assumptions about the signal waveform, polarization, source direction, or time of occurrence. No significant events were observed. %All candidate triggers were consistent with the expected background, As a result we set 90\% confidence upper limits on the rate of long-duration gravitational wave transients for different types of gravitational wave signals. We also show that the search is sensitive to sources in the Galaxy emitting at least $\sim$ \unit[$10^{-8}$]$\mathrm{M_{\odot} c^2}$ in gravitational waves.
- Nov 16 2017 astro-ph.HE gr-qc arXiv:1711.05578v1On June 8, 2017 at 02:01:16.49 UTC, a gravitational-wave signal from the merger of two stellar-mass black holes was observed by the two Advanced LIGO detectors with a network signal-to-noise ratio of 13. This system is the lightest black hole binary so far observed, with component masses $12^{+7}_{-2}\,M_\odot$ and $7^{+2}_{-2}\,M_\odot$ (90% credible intervals). These lie in the range of measured black hole masses in low-mass X-ray binaries, thus allowing us to compare black holes detected through gravitational waves with electromagnetic observations. The source's luminosity distance is $340^{+140}_{-140}$ Mpc, corresponding to redshift $0.07^{+0.03}_{-0.03}$. We verify that the signal waveform is consistent with the predictions of general relativity.
- Sep 28 2017 gr-qc astro-ph.HE arXiv:1709.09660v3On August 14, 2017 at 10:30:43 UTC, the Advanced Virgo detector and the two Advanced LIGO detectors coherently observed a transient gravitational-wave signal produced by the coalescence of two stellar mass black holes, with a false-alarm-rate of $\lesssim$ 1 in 27000 years. The signal was observed with a three-detector network matched-filter signal-to-noise ratio of 18. The inferred masses of the initial black holes are $30.5_{-3.0}^{+5.7}$ Msun and $25.3_{-4.2}^{+2.8}$ Msun (at the 90% credible level). The luminosity distance of the source is $540_{-210}^{+130}~\mathrm{Mpc}$, corresponding to a redshift of $z=0.11_{-0.04}^{+0.03}$. A network of three detectors improves the sky localization of the source, reducing the area of the 90% credible region from 1160 deg$^2$ using only the two LIGO detectors to 60 deg$^2$ using all three detectors. For the first time, we can test the nature of gravitational wave polarizations from the antenna response of the LIGO-Virgo network, thus enabling a new class of phenomenological tests of gravity.
- Jun 07 2017 gr-qc astro-ph.HE arXiv:1706.01812v1We describe the observation of GW170104, a gravitational-wave signal produced by the coalescence of a pair of stellar-mass black holes. The signal was measured on January 4, 2017 at 10:11:58.6 UTC by the twin advanced detectors of the Laser Interferometer Gravitational-Wave Observatory during their second observing run, with a network signal-to-noise ratio of 13 and a false alarm rate less than 1 in 70,000 years. The inferred component black hole masses are $31.2^{+8.4}_{-6.0}\,M_\odot$ and $19.4^{+5.3}_{-5.9}\,M_\odot$ (at the 90% credible level). The black hole spins are best constrained through measurement of the effective inspiral spin parameter, a mass-weighted combination of the spin components perpendicular to the orbital plane, $\chi_\mathrm{eff} = -0.12^{+0.21}_{-0.30}.$ This result implies that spin configurations with both component spins positively aligned with the orbital angular momentum are disfavored. The source luminosity distance is $880^{+450}_{-390}~\mathrm{Mpc}$ corresponding to a redshift of $z = 0.18^{+0.08}_{-0.07}$. We constrain the magnitude of modifications to the gravitational-wave dispersion relation and perform null tests of general relativity. Assuming that gravitons are dispersed in vacuum like massive particles, we bound the graviton mass to $m_g \le 7.7 \times 10^{-23}~\mathrm{eV}/c^2$. In all cases, we find that GW170104 is consistent with general relativity.
- Apr 25 2017 gr-qc arXiv:1704.06784v1Advanced LIGO's recent observations of gravitational waves (GWs) from merging binary black holes have opened up a unique laboratory to test general relativity (GR) in the highly relativistic regime. One of the tests used to establish the consistency of the first LIGO event with a binary black hole merger predicted by GR was the inspiral-merger-ringdown consistency test. This involves inferring the mass and spin of the remnant black hole from the inspiral (low-frequency) part of the observed signal and checking for the consistency of the inferred parameters with the same estimated from the post-inspiral (high-frequency) part of the signal. Based on the observed rate of binary black hole mergers, we expect the advanced GW observatories to observe hundreds of binary black hole mergers every year when operating at their design sensitivities, most of them with modest signal to noise ratios (SNRs). Anticipating such observations, this paper shows how constraints from a large number of events with modest SNRs can be combined to produce strong constraints on deviations from GR. Using kludge modified GR waveforms, we demonstrate how this test could identify certain types of deviations from GR if such deviations are present in the signal waveforms. We also study the robustness of this test against reasonable variations of a variety of different analysis parameters.
- Jun 16 2016 gr-qc astro-ph.CO arXiv:1606.04856v3The first observational run of the Advanced LIGO detectors, from September 12, 2015 to January 19, 2016, saw the first detections of gravitational waves from binary black hole mergers. In this paper we present full results from a search for binary black hole merger signals with total masses up to $100 M_\odot$ and detailed implications from our observations of these systems. Our search, based on general-relativistic models of gravitational wave signals from binary black hole systems, unambiguously identified two signals, GW150914 and GW151226, with a significance of greater than $5\sigma$ over the observing period. It also identified a third possible signal, LVT151012, with substantially lower significance, and with an 87% probability of being of astrophysical origin. We provide detailed estimates of the parameters of the observed systems. Both GW150914 and GW151226 provide an unprecedented opportunity to study the two-body motion of a compact-object binary in the large velocity, highly nonlinear regime. We do not observe any deviations from general relativity, and place improved empirical bounds on several high-order post-Newtonian coefficients. From our observations we infer stellar-mass binary black hole merger rates lying in the range $9-240 \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$. These observations are beginning to inform astrophysical predictions of binary black hole formation rates, and indicate that future observing runs of the Advanced detector network will yield many more gravitational wave detections.
- Jun 06 2016 gr-qc astro-ph.HE arXiv:1606.01210v1This paper presents updated estimates of source parameters for GW150914, a binary black-hole coalescence event detected by the Laser Interferometer Gravitational-wave Observatory (LIGO) on September 14, 2015 [1]. Reference presented parameter estimation [2] of the source using a 13-dimensional, phenomenological precessing-spin model (precessing IMRPhenom) and a 11-dimensional nonprecessing effective-one-body (EOB) model calibrated to numerical-relativity simulations, which forces spin alignment (nonprecessing EOBNR). Here we present new results that include a 15-dimensional precessing-spin waveform model (precessing EOBNR) developed within the EOB formalism. We find good agreement with the parameters estimated previously [2], and we quote updated component masses of $35^{+5}_{-3}\mathrm{M}_\odot$ and $30^{+3}_{-4}\mathrm{M}_\odot$ (where errors correspond to 90% symmetric credible intervals). We also present slightly tighter constraints on the dimensionless spin magnitudes of the two black holes, with a primary spin estimate $0.65$ and a secondary spin estimate $0.75$ at 90% probability. Reference [2] estimated the systematic parameter-extraction errors due to waveform-model uncertainty by combining the posterior probability densities of precessing IMRPhenom and nonprecessing EOBNR. Here we find that the two precessing-spin models are in closer agreement, suggesting that these systematic errors are smaller than previously quoted.
- Feb 12 2016 gr-qc astro-ph.HE arXiv:1602.03840v2On September 14, 2015, the Laser Interferometer Gravitational-wave Observatory (LIGO) detected a gravitational-wave transient (GW150914); we characterize the properties of the source and its parameters. The data around the time of the event were analyzed coherently across the LIGO network using a suite of accurate waveform models that describe gravitational waves from a compact binary system in general relativity. GW150914 was produced by a nearly equal mass binary black hole of $36^{+5}_{-4} M_\odot$ and $29^{+4}_{-4} M_\odot$; for each parameter we report the median value and the range of the 90% credible interval. The dimensionless spin magnitude of the more massive black hole is bound to be $<0.7$ (at 90% probability). The luminosity distance to the source is $410^{+160}_{-180}$ Mpc, corresponding to a redshift $0.09^{+0.03}_{-0.04}$ assuming standard cosmology. The source location is constrained to an annulus section of $610$ deg$^2$, primarily in the southern hemisphere. The binary merges into a black hole of $62^{+4}_{-4} M_\odot$ and spin $0.67^{+0.05}_{-0.07}$. This black hole is significantly more massive than any other inferred from electromagnetic observations in the stellar-mass regime.
- Feb 09 2016 gr-qc astro-ph.HE arXiv:1602.02453v2The coalescences of stellar-mass black-hole binaries through their inspiral, merger, and ringdown are among the most promising sources for ground-based gravitational-wave (GW) detectors. If a GW signal is observed with sufficient signal-to-noise ratio, the masses and spins of the black holes can be estimated from just the inspiral part of the signal. Using these estimates of the initial parameters of the binary, the mass and spin of the final black hole can be uniquely predicted making use of general-relativistic numerical simulations. In addition, the mass and spin of the final black hole can be independently estimated from the merger--ringdown part of the signal. If the binary black hole dynamics is correctly described by general relativity (GR), these independent estimates have to be consistent with each other. We present a Bayesian implementation of such a test of general relativity, which allows us to combine the constraints from multiple observations. Using kludge modified GR waveforms, we demonstrate that this test can detect sufficiently large deviations from GR, and outline the expected constraints from upcoming GW observations using the second-generation of ground-based GW detectors.
- Sep 23 2015 gr-qc arXiv:1509.06366v1In the adiabatic post-Newtonian (PN) approximation, the phase evolution of gravitational waves (GWs) from inspiralling compact binaries in quasicircular orbits is computed by equating the change in binding energy with the GW flux. This energy balance equation can be solved in different ways, which result in multiple approximants of the PN waveforms. Due to the poor convergence of the PN expansion, these approximants tend to differ from each other during the late inspiral. Which of these approximants should be chosen as templates for detection and parameter estimation of GWs from inspiraling compact binaries is not obvious. In this paper, we present estimates of the effective higher order (beyond the currently available 4PN and 3.5PN) non-spinning terms in the PN expansion of the binding energy and the GW flux that minimize the difference of multiple PN approximants (TaylorT1, TaylorT2, TaylorT4, TaylorF2) with effective one body waveforms calibrated to numerical relativity (EOBNR). We show that PN approximants constructed using the effective higher order terms show significantly better agreement (as compared to 3.5PN) with the inspiral part of the EOBNR. For non-spinning binaries with component masses $m_{1,2} \in [1.4 M_\odot, 15 M_\odot]$, most of the approximants have a match (faithfulness) of better than 99% with both EOBNR and each other.
- Information about the last stages of a binary neutron star inspiral and the final merger can be extracted from quasi-equilibrium configurations and dynamical evolutions. In this article, we construct quasi-equilibrium configurations for different spins, eccentricities, mass ratios, compactnesses, and equations of state. For this purpose we employ the SGRID code, which allows us to construct such data in previously inaccessible regions of the parameter space. In particular, we consider spinning neutron stars in isolation and in binary systems; we incorporate new methods to produce highly eccentric and eccentricity reduced data; we present the possibility of computing data for significantly unequal-mass binaries; and we create equal-mass binaries with individual compactness up to 0.23. As a proof of principle, we explore the dynamical evolution of three new configurations. First, we simulate a $q=2.06$ mass ratio which is the highest mass ratio for a binary neutron star evolved in numerical relativity to date. We find that mass transfer from the companion star sets in a few revolutions before merger and a rest mass of $\sim10^{-2}M_\odot$ is transferred between the two stars. This configuration also ejects a large amount of material during merger, imparting a substantial kick to the remnant. Second, we simulate the first merger of a precessing binary neutron star. We present the dominant modes of the gravitational waves for the precessing simulation, where a clear imprint of the precession is visible in the (2,1) mode. Finally, we quantify the effect of an eccentricity reduction procedure on the gravitational waveform. The procedure improves the waveform quality and should be employed in future precision studies, but also other errors, notably truncation errors, need to be reduced in order for the improvement due to eccentricity reduction to be effective. [abridged]
- Mar 10 2015 gr-qc arXiv:1503.02638v2It is now possible to compute linear in mass-ratio terms in the post-Newtonian (PN) expansion for compact binaries to very high orders using black hole perturbation theory applied to various invariants. For instance, a computation of the redshift invariant of a point particle in a circular orbit about a black hole in linear perturbation theory gives the linear-in-mass-ratio portion of the binding energy of a circular binary with arbitrary mass ratio. This binding energy, in turn, encodes the system's conservative dynamics. We give a method for extracting the analytic forms of these PN coefficients from high-accuracy numerical data using experimental mathematics techniques, notably an integer relation algorithm. Such methods should be particularly important when the calculations progress to the considerably more difficult case of perturbations of the Kerr metric. As an example, we apply this method to the redshift invariant in Schwarzschild. Here we obtain analytic coefficients to 12.5PN, and higher-order terms in mixed analytic-numerical form to 21.5PN, including analytic forms for the complete 13.5PN coefficient, and all the logarithmic terms at 13PN. At these high orders, an individual coefficient can have over 30 terms, including a wide variety of transcendental numbers, when written out in full. We are still able to obtain analytic forms for such coefficients from the numerical data through a careful study of the structure of the expansion. The structure we find also allows us to predict certain "leading logarithm"-type contributions to all orders. The additional terms in the expansion we obtain improve the accuracy of the PN series for the redshift observable, even in the very strong-field regime inside the innermost stable circular orbit, particularly when combined with exponential resummation.
- Aug 20 2014 gr-qc astro-ph.SR arXiv:1408.4136v2Binary neutron stars in circular orbits can be modeled as helically symmetric, i.e., stationary in a rotating frame. This symmetry gives rise to a first integral of the Euler equation, often employed for constructing equilibrium solutions via iteration. For eccentric orbits, however, the lack of helical symmetry has prevented the use of this method, and the numerical relativity community has often resorted to constructing initial data by superimposing boosted spherical stars without solving the Euler equation. The spuriously excited neutron star oscillations seen in evolutions of such data arise because such configurations lack the appropriate tidal deformations and are stationary in a linearly comoving---rather than rotating---frame. We consider eccentric configurations at apoapsis that are instantaneously stationary in a rotating frame. We extend the notion of helical symmetry to eccentric orbits, by approximating the elliptical orbit of each companion as instantaneously circular, using the ellipse's inscribed circle. The two inscribed helical symmetry vectors give rise to approximate instantaneous first integrals of the Euler equation throughout each companion. We use these integrals as the basis of a self-consistent iteration of the Einstein constraints to construct conformal thin-sandwich initial data for eccentric binaries. We find that the spurious stellar oscillations are reduced by at least an order of magnitude, compared with those found in evolutions of superposed initial data. The tidally induced oscillations, however, are physical and qualitatively similar to earlier evolutions. Finally, we show how to incorporate radial velocity due to radiation reaction in our inscribed helical symmetry vectors, which would allow one to obtain truly non-eccentric initial data when our eccentricity parameter $e$ is set to zero.
- May 08 2014 gr-qc arXiv:1405.1572v2(Abridged) High-order terms in the post-Newtonian (PN) expansions of various quantities for compact binaries exhibit a combinatorial increase in complexity, including ever-increasing numbers of transcendentals. Here we consider the gravitational wave energy flux at infinity from a point particle in a circular orbit around a Schwarzschild black hole, which is known to 22PN beyond the lowest-order Newtonian prediction, at which point each order has over 1000 terms. We introduce a factorization that considerably simplifies the spherical harmonic modes of the energy flux (and thus also the amplitudes of the spherical harmonic modes of the gravitational waves); it is likely that much of the complexity this factorization removes is due to curved-space wave propagation (e.g., tail effects). For the modes with azimuthal number l of 7 or greater, this factorization reduces the expressions for the modes that enter the 22PN total energy flux to pure integer PN series with rational coefficients, which amounts to a reduction of up to a factor of ~150 in the total number of terms in a given mode. The reduction in complexity becomes less dramatic for smaller l, due to the structure of the expansion, though the factorization is still able to remove all the half-integer PN terms. For the 22PN l = 2 modes, this factorization still reduces the total number of terms (and size) by a factor of ~10 and gives purely rational coefficients through 8PN. This factorization also improves the convergence of the series, though we find the exponential resummation introduced for the full energy flux by Isoyama et al. to be even more effective at improving the convergence of the individual modes, producing improvements of over four orders of magnitude over the original series for some modes. However, the exponential resummation is not as effective at simplifying the series, particularly for the higher-order modes.
- Jul 22 2013 gr-qc arXiv:1307.5307v3The Numerical-Relativity-Analytical-Relativity (NRAR) collaboration is a joint effort between members of the numerical relativity, analytical relativity and gravitational-wave data analysis communities. The goal of the NRAR collaboration is to produce numerical-relativity simulations of compact binaries and use them to develop accurate analytical templates for the LIGO/Virgo Collaboration to use in detecting gravitational-wave signals and extracting astrophysical information from them. We describe the results of the first stage of the NRAR project, which focused on producing an initial set of numerical waveforms from binary black holes with moderate mass ratios and spins, as well as one non-spinning binary configuration which has a mass ratio of 10. All of the numerical waveforms are analysed in a uniform and consistent manner, with numerical errors evaluated using an analysis code created by members of the NRAR collaboration. We compare previously-calibrated, non-precessing analytical waveforms, notably the effective-one-body (EOB) and phenomenological template families, to the newly-produced numerical waveforms. We find that when the binary's total mass is ~100-200 solar masses, current EOB and phenomenological models of spinning, non-precessing binary waveforms have overlaps above 99% (for advanced LIGO) with all of the non-precessing-binary numerical waveforms with mass ratios <= 4, when maximizing over binary parameters. This implies that the loss of event rate due to modelling error is below 3%. Moreover, the non-spinning EOB waveforms previously calibrated to five non-spinning waveforms with mass ratio smaller than 6 have overlaps above 99.7% with the numerical waveform with a mass ratio of 10, without even maximizing on the binary parameters.
- Interferometric detectors will very soon give us an unprecedented view of the gravitational-wave sky, and in particular of the explosive and transient Universe. Now is the time to challenge our theoretical understanding of short-duration gravitational-wave signatures from cataclysmic events, their connection to more traditional electromagnetic and particle astrophysics, and the data analysis techniques that will make the observations a reality. This paper summarizes the state of the art, future science opportunities, and current challenges in understanding gravitational-wave transients.
- Mar 14 2013 astro-ph.SR gr-qc arXiv:1303.3259v2We show that there is a direct relation between upper limits on (or potential future measurements of) the m = 2 quadrupole moments of slowly rotating neutron stars and the l = m = 2 deformation of the star's surface, in full general relativity, to first order in the perturbation. This relation only depends on the star's structure through its mass and radius. All one has to assume about the star's constituents is that the stress-energy tensor at its surface is that of a perfect fluid, which will be true with good accuracy in almost all the situations of interest, and that the magnetic field configuration there is force-free, which is likely to be a good approximation. We then apply this relation to the stars which have direct LIGO/Virgo bounds on their m = 2 quadrupole moment, below the spin-down limit, and compare with the expected surface deformations due to rotation. In particular, we find that LIGO observations have constrained the Crab pulsar's l = m = 2 surface deformation to be smaller than its l = 2, m = 0 deformation due to rotation, for all the causal equations of state we consider, a statement that would not have been able to be made just using the upper bounds on the l = m = 2 deformation from electromagnetic observations.
- Aug 28 2012 astro-ph.SR gr-qc arXiv:1208.5227v2We present a method for calculating the maximum elastic quadrupolar deformations of relativistic stars, generalizing the previous Newtonian, Cowling approximation integral given by [G. Ushomirsky et al., Mon. Not. R. Astron. Soc. 319, 902 (2000)]. (We also present a method for Newtonian gravity with no Cowling approximation.) We apply these methods to the m = 2 quadrupoles most relevant for gravitational radiation in three cases: crustal deformations, deformations of crystalline cores of hadron-quark hybrid stars, and deformations of entirely crystalline color superconducting quark stars. In all cases, we find suppressions of the quadrupole due to relativity compared to the Newtonian Cowling approximation, particularly for compact stars. For the crust these suppressions are up to a factor ~6, for hybrid stars they are up to ~4, and for solid quark stars they are at most ~2, with slight enhancements instead for low mass stars. We also explore ranges of masses and equations of state more than in previous work, and find that for some parameters the maximum quadrupoles can still be very large. Even with the relativistic suppressions, we find that 1.4 solar mass stars can sustain crustal quadrupoles of a few times 10^39 g cm^2 for the SLy equation of state or close to 10^40 g cm^2 for equations of state that produce less compact stars. Solid quark stars of 1.4 solar masses can sustain quadrupoles of around 10^44 g cm^2. Hybrid stars typically do not have solid cores at 1.4 solar masses, but the most massive ones (~2 solar masses) can sustain quadrupoles of a few times 10^41 g cm^2 for typical microphysical parameters and a few times 10^42 g cm^2 for extreme ones. All of these quadrupoles assume a breaking strain of 0.1 and can be divided by 10^45 g cm^2 to yield the fiducial "ellipticities" quoted elsewhere.
- Jul 07 2009 gr-qc arXiv:0907.0891v1(Abridged) By asymptotically matching a post-Newtonian (PN) metric to two tidally perturbed Schwarzschild metrics, we generate approximate initial data (in the form of a 4-metric) for a nonspinning black hole binary in a circular orbit. We carry out this matching through O(v^4) in the binary's orbital velocity v, so the resulting data are conformally curved. Far from the holes, we use the appropriate PN metric that accounts for retardation, which we construct using the highest-order PN expressions available to compute the binary's past history. The data set's uncontrolled remainders are thus O(v^5) throughout the timeslice; we also generate an extension to the data set that has uncontrolled remainders of O(v^6) in the purely PN portion of the timeslice (i.e., not too close to the holes). The resulting data are smooth, since we join all the metrics together by smoothly interpolating between them. We perform this interpolation using transition functions constructed to avoid introducing excessive additional constraint violations. Due to their inclusion of tidal deformations and outgoing radiation, these data should substantially reduce the initial spurious ("junk") radiation observed in current simulations that use conformally flat initial data. Such reductions in the nonphysical components of the initial data will be necessary for simulations to achieve the accuracy required to supply Advanced LIGO and LISA with the templates necessary for parameter estimation.