BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Richard O'Shaughnessy

Richard O'Shaughnessy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.COastro-ph.GAastro-ph.IMastro-ph.SRhep-phhep-thnucl-thastro-ph

Catalog footprint

What is connected

45works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

45 published item(s)

preprint2026arXiv

Analysis of GWTC-3 with fully precessing numerical relativity surrogate models

The third Gravitational-Wave Transient Catalog (GWTC-3) contains 90 binary coalescence candidates detected by the LIGO-Virgo-KAGRA Collaboration (LVK). We provide a re-analysis of binary black hole (BBH) events using a recently developed numerical relativity (NR) waveform surrogate model, NRSur7dq4, that includes all $\ell \leq 4$ spin-weighted spherical harmonic modes as well as the complete physical effects of precession. Properties of the remnant black holes' (BH's) mass, spin vector, and kick vector are found using an associated remnant surrogate model NRSur7dq4Remnant. Both NRSur7dq4 and NRSur7dq4Remnant models have errors comparable to numerical relativity simulations and allow for high-accuracy parameter estimates. We restrict our analysis to 47 BBH events that fall within the regime of validity of NRSur7dq4 (mass ratios greater than 1/6 and total masses greater than $60 M_{\odot}$). While for most of these events our results match the LVK analyses that were obtained using the semi-analytical models such as IMRPhenomXPHM and SEOBNRv4PHM, we find that for more than 20\% of events the NRSur7dq4 model recovers noticeably different measurements of black hole properties like the masses and spins, as well as extrinsic properties like the binary inclination and distance. For instance, GW150914_095045 exhibits noticeable differences in spin precession and spin magnitude measurements. Other notable findings include one event (GW191109_010717) that constrains the effective spin $χ_{eff}$ to be negative at a 99.3\% credible level and two events (GW191109_010717 and GW200129_065458) with well-constrained kick velocities. Furthermore, compared to the models used in the LVK analyses, NRSur7dq4 recovers a larger signal-to-noise ratio and/or Bayes factors for several events.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

GW200105: A detailed study of eccentricity in the neutron star-black hole binary

GW200105_162426 is the first neutron star-black hole merger to be confidently confirmed through either gravitational-wave or electromagnetic observations. Although initially analyzed after detection, the event has recently gained renewed attention following a study [Morras et al. arXiv:2503.15393] that employed a post-Newtonian inspiral-only waveform model and reported strong evidence for orbital eccentricity. In this work, we perform a detailed analysis of GW200105 using state-of-the-art effective-one-body waveform models. Importantly, we present the first study of this event utilizing a physically complete model that incorporates both orbital eccentricity and spin precession across the full inspiral, merger, and ringdown stages, along with higher-order gravitational wave modes. Our results support the presence of eccentricity in the signal, with zero eccentricity excluded from the 99% credible interval, but yielding a mass ratio closer to the original LIGO-Virgo-KAGRA analysis, differing from the findings of [Morras et al. arXiv:2503.15393]. Additionally, similar to a previous eccentric-only analysis [de Lluc Planas et al. Astrophys. J. 995, 47 (2025).], we observe a multimodal structure in the eccentricity posterior distribution. We conduct targeted investigations to understand the origin of this multimodality and complement our analysis with numerical relativity simulations to examine how the inclusion of eccentricity impacts the merger dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A Nuclear Equation of State Inferred from Stellar r-process Abundances

Binary neutron star mergers (NSMs) have been confirmed as one source of the heaviest observable elements made by the rapid neutron-capture (r-) process. However, modeling NSM outflows -- from the total ejecta masses to their elemental yields -- depends on the unknown nuclear equation of state (EOS) that governs neutron-star structure. In this work, we derive a phenomenological EOS by assuming that NSMs are the dominant sources of the heavy element material in metal-poor stars with r-process abundance patterns. We start with a population synthesis model to obtain a population of merging neutron star binaries and calculate their EOS-dependent elemental yields. Under the assumption that these mergers were responsible for the majority of r-process elements in the metal-poor stars, we find parameters representing the EOS for which the theoretical NSM yields reproduce the derived abundances from observations of metal-poor stars. For our proof-of-concept assumptions, we find an EOS that is slightly softer than, but still in agreement with, current constraints, e.g., by the Neutron Star Interior Composition Explorer, with $R_{1.4}=12.25\pm 0.03$~km and $M_{\textrm TOV}$ of $2.17\pm 0.03$~M$_\odot$(statistical uncertainties, neglecting modeling systematics).

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cost Minimization in Acquisition for Gravitational Wave Surrogate Modeling

Gravitational wave science is dependent upon expensive numerical simulations, which provide the foundational understanding of binary merger radiation needed to interpret observations of massive binary black holes. The high cost of these simulations limits large-scale campaigns to explore the binary black hole parameter space. Surrogate models have been developed to efficiently interpolate between simulation results, but these models require a sufficiently comprehensive sample to train on. Acquisition functions can be used to identify points in the domain for simulation. We develop a new acquisition function which accounts for the cost of simulating new points. We show that when applied to a 3D domain of binary mass ratio and dimensionless spins, the accumulated cost of simulation is reduced by a factor of about 10.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Implementation of a generalized precession parameter in the RIFT parameter estimation algorithm

Since the initial discovery of gravitational waves in 2015, significant developments have been made towards waveform interpretation and estimation of compact binary source parameters. We present herein an implementation of the generalized precession parameter $\langle χ_p \rangle$, which averages over all angular variations on the precession timescale, within the RIFT parameter estimation framework. Relative to the originally-proposed precession parameter $χ_p$, which characterizes the single largest dynamical spin in a binary, this new parameter $\langle χ_p \rangle$ has a unique domain $1 < \langle χ_p \rangle < 2$, which is exclusive to binaries with two precessing spins. After reviewing the physical differences between these two parameters, we describe how $\langle χ_p \rangle$ was implemented in RIFT and apply it to all 36 events from the second half of the Advanced LIGO and Advanced Virgo third operating run (O3b). In O3b, ten events show significant amounts of precession $\langle χ_p \rangle > 0.5$. Of particular interest is GW191109_010717; we show it has a $\sim28\%$ probability that the originating system necessarily contains two misaligned spins.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Normal Approximate Likelihoods to Gravitational Wave Events

Gravitational wave observations of quasicircular compact binary mergers in principle provide an arbitrarily complex likelihood over eight independent intrinsic parameters: the masses and spins of the two merging objects. In this work, we demonstrate by example that a simple normal approximation over fewer (usually, three) effective dimensions provides a very accurate representation of the likelihood, and allows us to replicate the eight-dimensional posterior over the mass and spin degrees of freedom. Alongside this paper, we provide the parameters for multivariate normal fits for each event published in GWTC-1 and GWTC-2, using the posterior samples from the catalog for each associated release. These normal approximations provide a highly efficient way to characterize gravitational wave observations when combining large numbers of events.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Application of the third RIT binary black hole simulations catalog to parameter estimation of gravitational waves signals from the LIGO-Virgo O1/O2 observational runs

Using exclusively the 777 full numerical waveforms of the third Binary Black Holes RIT catalog, we reanalyze the ten black hole merger signals reported in LIGO/Virgo's O1/O2 observation runs. We obtain binary parameters, extrinsic parameters, and the remnant properties of these gravitational waves events which are consistent with, but not identical to previously presented results. We have also analyzed three additional events (GW170121, GW170304, GW170727) reported in Venumadhav et al. 2019, and found closely matching parameters. We finally assess the accuracy of our waveforms with convergence studies applied to O1/O2 events and found them adequate for current estimation of parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Black Hole Coagulation: Modeling Hierarchical Mergers in Black Hole Populations

Data from the LIGO and Virgo detectors has confirmed that stellar-mass black holes can merge within a Hubble time, leaving behind massive remnant black holes. In some astrophysical environments such as globular clusters and AGN disks, it may be possible for these remnants to take part in further compact-object mergers, producing a population of hierarchically formed black holes. In this work, we present a parameterized framework for describing the population of binary black hole mergers, while self-consistently accounting for hierarchical mergers. The framework casts black holes as particles in a box which can collide based on an effective cross-section, but allows inputs from more detailed astrophysical simulations. Our approach is relevant to any population which is comprised of second or higher generation black holes, such as primordial black holes or dense cluster environments. We describe some possible inputs to this generic model and their effects on the black hole merger populations, and use the model to perform Bayesian inference on the catalog of black holes from LIGO and Virgo's first two observing runs. We find that models with a high rate of hierarchical mergers are disfavored, consistent with previous population analyses. Future gravitational-wave events will further constrain the inputs to this generic hierarchical merger model, enabling a deeper look into the formation environments of binary black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Combining Electromagnetic and Gravitational-Wave Constraints on Neutron-Star Masses and Radii

We perform a joint Bayesian inference of neutron-star mass and radius constraints based on GW170817, observations of quiescent low-mass X-ray binaries (QLMXBs), photospheric radius expansion X-ray bursts (PREs), and X-ray timing observations of J0030+0451. With this data set, the form of the prior distribution still has an impact on the posterior mass-radius (MR) curves and equation of state (EOS), but this impact is smaller than recently obtained when considering QLMXBs alone. We analyze the consistency of the electromagnetic data by including an "intrinsic scattering" contribution to the uncertainties, and find only a slight broadening of the posteriors. This suggests that the gravitational-wave and electromagnetic observations of neutron-star structure are providing a consistent picture of the neutron-star mass-radius curve and the EOS.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Hybrid waveforms for generic precessing binaries for gravitational-wave data analysis

We construct hybrid binary black holes merger waveforms using analytical model waveforms for the early inspiral phase and numerical relativity waveforms for late inspiral to merger and post merger phases. To hybridize analytical and numerical waveforms, we first perform a 3-dimensional rotation to align the instantaneous orbital planes associated with the two waveforms at some fiducial frequency, we then find appropriate phase and time translations that maximize the overlap of the two waveforms in a hybridization interval. We discuss the accuracy and limitations for hybrids constructed by this procedure in the context of LIGO-Virgo-KAGRA observations. Our goal is to hybridize waveforms for more generic precessing binaries and construct longer waveforms that are sufficiently accurate for the parameter estimation techniques for upcoming LIGO observations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Impact of subdominant modes on the interpretation of gravitational-wave signals from heavy binary black hole systems

Over the past year, a handful of new gravitational wave models have been developed to include multiple harmonic modes thereby enabling for the first time fully Bayesian inference studies including higher modes to be performed. Using one recently-developed numerical relativity surrogate model, NRHybSur3dq8, we investigate the importance of higher modes on parameter inference of coalescing massive binary black holes. We focus on examples relevant to the current three-detector network of observatories, with a detector-frame mass set to $120 M_\odot$ and with signal amplitude values that are consistent with plausible candidates for the next few observing runs. We show that for such systems the higher mode content will be important for interpreting coalescing binary black holes, reducing systematic bias, and computing properties of the remnant object. Even for comparable-mass binaries and at low signal amplitude, the omission of higher modes can influence posterior probability distributions. We discuss the impact of our results on source population inference and self-consistency tests of general relativity. Our work can be used to better understand asymmetric binary black hole merger events, such as GW190412. Higher modes are critical for such systems, and their omission usually produces substantial parameter biases.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inferring Parameters of GW170502: The Loudest Intermediate-mass Black Hole Trigger in LIGO's O1/O2 data

Gravitational wave (GW) measurements provide the most robust constraints of the mass of astrophysical black holes. Using state-of-the-art GW signal models and a unique parameter estimation technique, we infer the source parameters of the loudest marginal trigger, GW170502, found by LIGO from 2015 to 2017. If this trigger is assumed to be a binary black hole merger, we find it corresponds to a total mass in the source frame of $157^{+55}_{-41}~\rm{M}_\odot$ at redshift $z=1.37^{+0.93}_{-0.64}$. The primary and secondary black hole masses are constrained to $94^{+44}_{-28}~\rm{M}_{\odot}$ and $62^{+30}_{-25}~\rm{M}_{\odot}$ respectively, with 90\% confidence. Across all signal models, we find $\gtrsim 70\%$ probability for the effective spin parameter $χ_\mathrm{eff}>0.1$. Furthermore, we find that the inclusion of higher-order modes in the analysis narrows the confidence region for the primary black hole mass by 10\%, however, the evidence for these modes in the data remains negligible. The techniques outlined in this study could lead to robust inference of the physical parameters for all intermediate-mass black hole binary candidates $(\gtrsim100~\mathrm{M}_\odot)$ in the current GW network.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inferring the neutron star equation of state simultaneously with the population of merging neutron stars

Observations of the properties of multiple coalescing neutron stars will simultaneously provide insight into neutron star mass and spin distribution, the neutron star merger rate, and the nuclear equation of state. Not all merging binaries containing neutron stars are expected to be identical. Plausible sources of diversity in these coalescing binaries can arise from a broad or multi-peaked NS mass distribution; the effect of different and extreme NS natal spins; the possibility of NS-BH mergers; or even the possibility of phase transitions, allowing for NS with similar mass but strongly divergent radius. In this work, we provide a concrete algorithm to combine all information obtained from GW measurements into a joint constraint on the NS merger rate, the distribution of NS properties, and the nuclear equation of state. Using a concrete example, we show how biased mass distribution inferences can significantly impact the recovered equation of state, even in the small-$N$ limit. With the same concrete example, we show how small-$N$ observations could identify a bimodal mass and spin distribution for merging NS simultaneously with the EOS. Our concordance approach can be immediately generalized to incorporate other observational constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Inspiral-Merger-Ringdown Consistency Tests with Higher Modes on Gravitational Signals from the Second Observing Run of LIGO and Virgo

Current tests of General Relativity are performed using approximations which neglect a key feature of complete solution of Einstein's theory: higher-order modes. Our analysis will reassess these tests, including these higher-order mode effects. We have chosen to perform inspiral-merger-ringdown consistency tests on the gravitational transients detected by LIGO and Virgo during the observing run O2. We use an approximant which includes all higher modes with $\ell \le 4$ (NRSur7dq2) and then, for the most interesting cases, we repeat the tests involving fits on Numerical Relativity simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Compact Binary Merger Rates: Comparison with LIGO/Virgo Upper Limits

We compare evolutionary predictions of double compact object merger rate densities with initial and forthcoming LIGO/Virgo upper limits. We find that: (i) Due to the cosmological reach of advanced detectors, current conversion methods of population synthesis predictions into merger rate densities are insufficient. (ii) Our optimistic models are a factor of 18 below the initial LIGO/Virgo upper limits for BH-BH systems, indicating that a modest increase in observational sensitivity (by a factor of 2.5) may bring the first detections or first gravitational wave constraints on binary evolution. (iii) Stellar-origin massive BH-BH mergers should dominate event rates in advanced LIGO/Virgo and can be detected out to redshift z=2 with templates including inspiral, merger, and ringdown. Normal stars (<150 Msun) can produce such mergers with total redshifted mass up to 400 Msun. (iv) High black hole natal kicks can severely limit the formation of massive BH-BH systems (both in isolated binary and in dynamical dense cluster evolution), and thus would eliminate detection of these systems even at full advanced LIGO/Virgo sensitivity. We find that low and high black hole natal kicks are allowed by current observational electromagnetic constraints. (v) The majority of our models yield detections of all types of mergers (NS-NS, BH-NS, BH-BH) with advanced detectors. Numerous massive BH-BH merger detections will indicate small (if any) natal kicks for massive BHs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Distinguishing black-hole spin-orbit resonances by their gravitational wave signatures. II: Full parameter estimation

Gravitational waves from coalescing binary black holes encode the evolution of their spins prior to merger. In the post-Newtonian regime and on the precession timescale, this evolution has one of three morphologies, with the spins either librating around one of two fixed points ("resonances") or circulating freely. In this work we perform full parameter estimation on resonant binaries with fixed masses and spin magnitudes, changing three parameters: a conserved "projected effective spin" $ξ$ and resonant family $ΔΦ=0,π$ (which uniquely label the source), the inclination $θ_{JN}$ of the binary's total angular momentum with respect to the line of sight (which determines the strength of precessional effects in the waveform), and the signal amplitude. We demonstrate that resonances can be distinguished for a wide range of binaries, except for highly symmetric configurations where precessional effects are suppressed. Motivated by new insight into double-spin evolution, we introduce new variables to characterize precessing black hole binaries which naturally reflects the timescale separation of the system and therefore better encode the dynamical information carried by gravitational waves.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Modeling the source of GW150914 with targeted numerical-relativity simulations

In fall of 2015, the two LIGO detectors measured the gravitational wave signal GW150914, which originated from a pair of merging black holes. In the final 0.2 seconds (about 8 gravitational-wave cycles) before the amplitude reached its maximum, the observed signal swept up in amplitude and frequency, from 35 Hz to 150 Hz. The theoretical gravitational-wave signal for merging black holes, as predicted by general relativity, can be computed only by full numerical relativity, because analytic approximations fail near the time of merger. Moreover, the nearly-equal masses, moderate spins, and small number of orbits of GW150914 are especially straightforward and efficient to simulate with modern numerical-relativity codes. In this paper, we report the modeling of GW150914 with numerical-relativity simulations, using black-hole masses and spins consistent with those inferred from LIGO's measurement. In particular, we employ two independent numerical-relativity codes that use completely different analytical and numerical methods to model the same merging black holes and to compute the emitted gravitational waveform; we find excellent agreement between the waveforms produced by the two independent codes. These results demonstrate the validity, impact, and potential of current and future studies using rapid-response, targeted numerical-relativity simulations for better understanding gravitational-wave observations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effects of host galaxy properties on merging compact binaries detectable by LIGO

Cosmological simulations of galaxy formation can produce present-day galaxies with a large range of assembly and star formation histories. A detailed study of the metallicity evolution and star formation history of such simulations can assist in predicting LIGO-detectable compact object binary mergers. Recent simulations of compact binary evolution suggest the compact object merger rate depends sensitively on the progenitor's metallicity. Rare low-metallicity star formation during galaxy assembly can produce more detected compact binaries than typical star formation. Using detailed simulations of galaxy and chemical evolution, we determine how sensitively the compact binary populations of galaxies with similar present-day appearance depend on the details of their assembly. We also demonstrate by concrete example the extent to which dwarf galaxies overabundantly produce compact binary mergers, particularly binary black holes, relative to more massive galaxies. We discuss the implications for transient multimessenger astronomy with compact binary sources.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The first gravitational-wave source from the isolated evolution of two 40-100 Msun stars

The merger of two massive 30 Msun black holes has been detected in gravitational waves (1,GW150914). This discovery validates recent predictions (2-4) that massive binary black holes would constitute the first detection. However, previous calculations have not sampled the relevant binary black hole progenitors---massive, low-metallicity binary stars---with sufficient accuracy and input physics to enable robust predictions to better than several orders of magnitude (5-10). Here we report a suite of high-precision numerical simulations of binary black hole formation via the evolution of isolated binary stars, providing a framework to interpret GW150914 and predict the properties of subsequent binary black hole gravitational-wave events. Our models imply that these events form in an environment where the metallicity is less than 10 percent of solar; have initial masses of 40-100 Msun; and interact through mass transfer and a common envelope phase. Their progenitors likely form either at 2 Gyr, or somewhat less likely, at 11 Gyr after the Big Bang. Most binary black holes form without supernova explosions, and their spins are nearly unchanged since birth, but do not have to be parallel. The classical field formation of binary black holes proposed in this study, with low natal kicks and restricted common envelope evolution, produces 40 times more binary black holes than dynamical formation channels involving globular clusters (11) and is comparable to the rate from homogeneous evolution channels (12-15). Our calculations predict detections of about 1,000 black hole mergers per year with total mass of 20-80 Msun once second generation ground-based gravitational wave observatories reach full sensitivity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effective potentials and morphological transitions for binary black-hole spin precession

We derive an effective potential for binary black-hole (BBH) spin precession at second post-Newtonian order. This effective potential allows us to solve the orbit-averaged spin-precession equations analytically for arbitrary mass ratios and spins. These solutions are quasiperiodic functions of time: after a fixed period the BBH spins return to their initial relative orientations and jointly precess about the total angular momentum by a fixed angle. Using these solutions, we classify BBH spin precession into three distinct morphologies between which BBHs can transition during their inspiral. We also derive a precession-averaged evolution equation for the total angular momentum that can be integrated on the radiation-reaction time and identify a new class of spin-orbit resonances that can tilt the direction of the total angular momentum during the inspiral. Our new results will help efforts to model and interpret gravitational waves from generic BBH mergers and predict the distributions of final spins and gravitational recoils.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multi-timescale analysis of phase transitions in precessing black-hole binaries

The dynamics of precessing binary black holes (BBHs) in the post-Newtonian regime has a strong timescale hierarchy: the orbital timescale is very short compared to the spin-precession timescale which, in turn, is much shorter than the radiation-reaction timescale on which the orbit is shrinking due to gravitational-wave emission. We exploit this timescale hierarchy to develop a multi-scale analysis of BBH dynamics elaborating on the analysis of Kesden et al. (2015). We solve the spin-precession equations analytically on the precession time and then implement a quasi-adiabatic approach to evolve these solutions on the longer radiation-reaction time. This procedure leads to an innovative "precession-averaged" post-Newtonian approach to studying precessing BBHs. We use our new solutions to classify BBH spin precession into three distinct morphologies, then investigate phase transitions between these morphologies as BBHs inspiral. These precession-averaged post-Newtonian inspirals can be efficiently calculated from arbitrarily large separations, thus making progress towards bridging the gap between astrophysics and numerical relativity.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Precessional instability in binary black holes with aligned spins

Binary black holes on quasicircular orbits with spins aligned with their orbital angular momentum have been testbeds for analytic and numerical relativity for decades, not least because symmetry ensures that such configurations are equilibrium solutions to the spin-precession equations. In this work, we show that these solutions can be unstable when the spin of the higher-mass black hole is aligned with the orbital angular momentum and the spin of the lower-mass black hole is anti-aligned. Spins in these configurations are unstable to precession to large misalignment when the binary separation $r$ is between the values $r_{\rm ud\pm}= (\sqrt{χ_1} \pm \sqrt{q χ_2})^4 (1-q)^{-2} M$, where $M$ is the total mass, $q \equiv m_2/m_1$ is the mass ratio, and $χ_1$ ($χ_2$) is the dimensionless spin of the more (less) massive black hole. This instability exists for a wide range of spin magnitudes and mass ratios and can occur in the strong-field regime near merger. We describe the origin and nature of the instability using recently developed analytical techniques to characterize fully generic spin precession. This instability provides a channel to circumvent astrophysical spin alignment at large binary separations, allowing significant spin precession prior to merger affecting both gravitational-wave and electromagnetic signatures of stellar-mass and supermassive binary black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Rapid gravitational wave parameter estimation with a single spin: Systematic uncertainties in parameter estimation with the SpinTaylorF2 approximation

Reliable low-latency gravitational wave parameter estimation is essential to target limited electromagnetic followup facilities toward astrophysically interesting and electromagnetically relevant sources of gravitational waves. In this study, we examine the tradeoff between speed and accuracy. Specifically, we estimate the astrophysical relevance of systematic errors in the posterior parameter distributions derived using a fast-but-approximate waveform model, SpinTaylorF2 (STF2), in parameter estimation with lalinference_mcmc. Though efficient, the STF2 approximation to compact binary inspiral employs approximate kinematics (e.g., a single spin) and an approximate waveform (e.g., frequency domain versus time domain). More broadly, using a large astrophysically-motivated population of generic compact binary merger signals, we report on the effectualness and limitations of this single-spin approximation as a method to infer parameters of generic compact binary sources. For most low-mass compact binary sources, we find that the STF2 approximation estimates compact binary parameters with biases comparable to systematic uncertainties in the waveform. We illustrate by example the effect these systematic errors have on posterior probabilities most relevant to low-latency electromagnetic followup: whether the secondary is has a mass consistent with a neutron star; whether the masses, spins, and orbit are consistent with that neutron star's tidal disruption; and whether the binary's angular momentum axis is oriented along the line of sight.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Robust parameter estimation for compact binaries with ground-based gravitational-wave observations using the LALInference software library

The Advanced LIGO and Advanced Virgo gravitational wave (GW) detectors will begin operation in the coming years, with compact binary coalescence events a likely source for the first detections. The gravitational waveforms emitted directly encode information about the sources, including the masses and spins of the compact objects. Recovering the physical parameters of the sources from the GW observations is a key analysis task. This work describes the LALInference software library for Bayesian parameter estimation of compact binary signals, which builds on several previous methods to provide a well-tested toolkit which has already been used for several studies. We show that our implementation is able to correctly recover the parameters of compact binary signals from simulated data from the advanced GW detectors. We demonstrate this with a detailed comparison on three compact binary systems: a binary neutron star, a neutron star black hole binary and a binary black hole, where we show a cross-comparison of results obtained using three independent sampling algorithms. These systems were analysed with non-spinning, aligned spin and generic spin configurations respectively, showing that consistent results can be obtained even with the full 15-dimensional parameter space of the generic spin configurations. We also demonstrate statistically that the Bayesian credible intervals we recover correspond to frequentist confidence intervals under correct prior assumptions by analysing a set of 100 signals drawn from the prior. We discuss the computational cost of these algorithms, and describe the general and problem-specific sampling techniques we have used to improve the efficiency of sampling the compact binary coalescence parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Testing General Relativity with Present and Future Astrophysical Observations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A more effective coordinate system for parameter estimation of precessing compact binaries from gravitational waves

Ground-based gravitational wave detectors are sensitive to a narrow range of frequencies, effectively taking a snapshot of merging compact-object binary dynamics just before merger. We demonstrate that by adopting analysis parameters that naturally characterize this 'picture', the physical parameters of the system can be extracted more efficiently from the gravitational wave data, and interpreted more easily. We assess the performance of MCMC parameter estimation in this physically intuitive coordinate system, defined by (a) a frame anchored on the binary's spins and orbital angular momentum and (b) a time at which the detectors are most sensitive to the binary's gravitational wave emission. Using anticipated noise curves for the advanced-generation LIGO and Virgo gravitational wave detectors, we find that this careful choice of reference frame and reference time significantly improves parameter estimation efficiency for BNS, NS-BH, and BBH signals.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Accurate and efficient waveforms for compact binaries on eccentric orbits

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Distinguishing black-hole spin-orbit resonances by their gravitational-wave signatures

If binary black holes form following the successive core collapses of sufficiently massive binary stars, precessional dynamics may align their spins $\mathbf S_1$ and $\mathbf S_2$ and the orbital angular momentum $\mathbf L$ into a plane in which they jointly precess about the total angular momentum $\mathbf J$. These spin orientations are known as spin-orbit resonances since $\mathbf S_1$, $\mathbf S_2$, and $\mathbf L$ all precess at the same frequency to maintain their planar configuration. Two families of such spin-orbit resonances exist, alike in dignity but differentiated by whether the components of the two spins in the orbital plane are either aligned or antialigned. The fraction of binary black holes in each family is determined by the stellar evolution of their progenitors, so if gravitational-wave detectors could measure this fraction they could provide important insights into astrophysical formation scenarios for binary black holes. In this paper, we show that even under the conservative assumption that binary black holes are observed along the direction of $\mathbf J$ (where precession-induced modulations to the gravitational waveforms are minimized), the waveforms of many members of each resonant family can be distinguished from all members of the opposing family in events with signal-to-noise ratios $ρ\simeq 10$, typical of those expected for the first detections with Advanced LIGO/Virgo. We hope that our preliminary findings inspire a greater appreciation of the capability of gravitational-wave detectors to constrain stellar astrophysics and stimulate further studies of the distinguishability of spin-orbit resonant families in more expanded regions of binary black-hole parameter space.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Comparing gravitational waves from nonprecessing and precessing black hole binaries in the corotating frame

Previous analytic and numerical calculations suggest that, at each instant, the emission from a precessing black hole binary closely resembles the emission from a nonprecessing analog. In this paper we quantitatively explore the validity and limitations of that correspondence, extracting the radiation from a large collection of roughly two hundred generic black hole binary merger simulations both in the simulation frame and in a corotating frame that tracks precession. To a first approximation, the corotating-frame waveforms resemble nonprecessing analogs, based on similarity over a band-limited frequency interval defined using a fiducial detector (here, advanced LIGO) and the source's total mass $M$. By restricting attention to masses $M\in 100, 1000 M_\odot$, we insure our comparisons are sensitive only to our simulated late-time inspiral, merger, and ringdown signals. In this mass region, every one of our precessing simulations can be fit by some physically similar member of the \texttt{IMRPhenomB} phenomenological waveform family to better than 95%; most fit significantly better. The best-fit parameters at low and high mass correspond to natural physical limits: the pre-merger orbit and post-merger perturbed black hole. Our results suggest that physically-motivated synthetic signals can be derived by viewing radiation from suitable nonprecessing binaries in a suitable nonintertial reference frame. While a good first approximation, precessing systems have degrees of freedom (i.e., the transverse spins) which a nonprecessing simulation cannot reproduce. We quantify the extent to which these missing degrees of freedom limit the utility of synthetic precessing signals for detection and parameter estimation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Double Compact Objects II: Cosmological Merger Rates

The development of advanced gravitational wave (GW) observatories, such as Advanced LIGO and Advanced Virgo, provides impetus to refine theoretical predictions for what these instruments might detect. In particular, with the range increasing by an order of magnitude, the search for GW sources is extending beyond the "local" Universe and out to cosmological distances. Double compact objects (neutron star-neutron star (NS-NS), black hole-neutron star (BH-NS) and black hole-black hole (BH-BH) systems) are considered to be the most promising gravitational wave sources. In addition, NS-NS and/or BH-NS systems are thought to be the progenitors of gamma ray bursts (GRBs), and may also be associated with kilonovae. In this paper we present the merger event rates of these objects as a function of cosmological redshift. We provide the results for four cases, each one investigating a different important evolution parameter of binary stars. Each case is also presented for two metallicity evolution scenarios. We find that (i) in most cases NS-NS systems dominate the merger rates in the local Universe, while BH-BH mergers dominate at high redshift; (ii) BH-NS mergers are less frequent than other sources per unit volume, for all time; and (iii) natal kicks may alter the observable properties of populations in a significant way, allowing the underlying models of binary evolution and compact object formation to be easily distinguished. This is the second paper in a series of three. The third paper will focus on calculating the detection rates of mergers by gravitational wave telescopes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Reconstructing the sky location of gravitational-wave detected compact binary systems: methodology for testing and comparison

The problem of reconstructing the sky position of compact binary coalescences detected via gravitational waves is a central one for future observations with the ground-based network of gravitational-wave laser interferometers, such as Advanced LIGO and Advanced Virgo. Different techniques for sky localisation have been independently developed. They can be divided in two broad categories: fully coherent Bayesian techniques, which are high-latency and aimed at in-depth studies of all the parameters of a source, including sky position, and "triangulation-based" techniques, which exploit the data products from the search stage of the analysis to provide an almost real-time approximation of the posterior probability density function of the sky location of a detection candidate. These techniques have previously been applied to data collected during the last science runs of gravitational-wave detectors operating in the so-called initial configuration. Here, we develop and analyse methods for assessing the self-consistency of parameter estimation methods and carrying out fair comparisons between different algorithms, addressing issues of efficiency and optimality. These methods are general, and can be applied to parameter estimation problems other than sky localisation. We apply these methods to two existing sky localisation techniques representing the two above-mentioned categories, using a set of simulated inspiral-only signals from compact binary systems with total mass $\le 20\,M_\odot$ and non-spinning components. We compare the relative advantages and costs of the two techniques and show that sky location uncertainties are on average a factor $\approx 20$ smaller for fully coherent techniques than for the specific variant of the "triangulation-based" technique used during the last science runs, at the expense of a factor $\approx 1000$ longer processing time.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Resonant-plane locking and spin alignment in stellar-mass black-hole binaries: a diagnostic of compact-binary formation

We study the influence of astrophysical formation scenarios on the precessional dynamics of spinning black-hole binaries by the time they enter the observational window of second- and third-generation gravitational-wave detectors, such as Advanced LIGO/Virgo, LIGO-India, KAGRA and the Einstein Telescope. Under the plausible assumption that tidal interactions are efficient at aligning the spins of few-solar mass black-hole progenitors with the orbital angular momentum, we find that black-hole spins should be expected to preferentially lie in a plane when they become detectable by gravitational-wave interferometers. This "resonant plane" is identified by the conditions ΔΦ=0° or ΔΦ=+/-180°, where ΔΦis the angle between the components of the black-hole spins in the plane orthogonal to the orbital angular momentum. If the angles ΔΦcan be accurately measured for a large sample of gravitational-wave detections, their distribution will constrain models of compact binary formation. In particular, it will tell us whether tidal interactions are efficient and whether a mechanism such as mass transfer, stellar winds, or supernovae can induce a mass-ratio reversal (so that the heavier black hole is produced by the initially lighter stellar progenitor). Therefore our model offers a concrete observational link between gravitational-wave measurements and astrophysics. We also hope that it will stimulate further studies of precessional dynamics, gravitational-wave template placement and parameter estimation for binaries locked in the resonant plane.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Asymptotic frame selection for binary black hole spacetimes II: Post-Newtonian limit

One way to select a preferred frame from gravitational radiation is via the principal axes of < L L>, an average of the action of rotation group generators on the Weyl tensor at asymptotic infinity. In this paper we evaluate this time-domain average for a quasicircular binary using approximate (post-Newtonian) waveforms. For nonprecessing unequal-mass binaries, we show the dominant eigenvector of this tensor lies along the orbital angular momentum. For precessing binaries, this frame is not generally aligned with either the orbital or total angular momentum, working to leading order in the spins. The difference between these two quantities grows with time, as the binary approaches the end of the inspiral and both precession and higher harmonics become more significant.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Double Compact Objects I: The Significance Of The Common Envelope On Merger Rates

The last decade of observational and theoretical developments in stellar and binary evolution provides an opportunity to incorporate major improvements to the predictions from populations synthesis models. We compute the Galactic merger rates for NS-NS, BH-NS, and BH-BH mergers with the StarTrack code. The most important revisions include: updated wind mass loss rates (allowing for stellar mass black holes up to $80 \msun$), a realistic treatment of the common envelope phase (a process that can affect merger rates by 2--3 orders of magnitude), and a qualitatively new neutron star/black hole mass distribution (consistent with the observed "mass gap"). Our findings include: (i) The binding energy of the envelope plays a pivotal role in determining whether a binary merges within a Hubble time. (ii) Our description of natal kicks from supernovae plays an important role, especially for the formation of BH-BH systems. (iii) The masses of BH-BH systems can be substantially increased in the case of low metallicities or weak winds. (iv) Certain combinations of parameters underpredict the Galactic NS-NS merger rate, and can be ruled out. {\em (v)} Models incorporating delayed supernovae do not agree with the observed NS/BH "mass gap", in accordance with our previous work. This is the first in a series of three papers. The second paper will study the merger rates of double compact objects as a function of redshift, star formation rate, and metallicity. In the third paper we will present the detection rates for gravitational wave observatories, using up-to-date signal waveforms and sensitivity curves.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational waves from BH-NS binaries: Effective Fisher matrices and parameter estimation using higher harmonics

Inspiralling black hole-neutron star (BH-NS) binaries emit a complicated gravitational wave signature, produced by multiple harmonics sourced by their strong local gravitational field and further modulated by the orbital plane's precession. Some features of this complex signal are easily accessible to ground-based interferometers (e.g., the rate of change of frequency); others less so (e.g., the polarization content); and others unavailable (e.g., features of the signal out of band). For this reason, an ambiguity function (a diagnostic of dissimilarity) between two such signals varies on many parameter scales and ranges. In this paper, we present a method for computing an approximate, effective Fisher matrix from variations in the ambiguity function on physically pertinent scales which depend on the relevant signal to noise ratio. As a concrete example, we explore how higher harmonics improve parameter measurement accuracy. As previous studies suggest, for our fiducial BH-NS binaries and for plausible signal amplitudes, we see that higher harmonics at best marginally improve our ability to measure parameters. For non-precessing binaries, these Fisher matrices separate into intrinsic (mass, spin) and extrinsic (geometrical) parameters; higher harmonics principally improve our knowledge about the line of sight. For the precessing binaries, the extra information provided by higher harmonics is distributed across several parameters. We provide concrete estimates for measurement accuracy, using coordinates adapted to the precession cone in the detector's sensitive band.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Nonspinning searches for spinning binaries in ground-based detector data: Amplitude and mismatch predictions in the constant precession cone approximation

Current searches for compact binary mergers by ground-based gravitational-wave detectors assume for simplicity the two bodies are not spinning. If the binary contains compact objects with significant spin, then this can reduce the sensitivity of these searches, particularly for black hole--neutron star binaries. In this paper we investigate the effect of neglecting precession on the sensitivity of searches for spinning binaries using non-spinning waveform models. We demonstrate that in the sensitive band of Advanced LIGO, the angle between the binary's orbital angular momentum and its total angular momentum is approximately constant. Under this \emph{constant precession cone} approximation, we show that the gravitational-wave phasing is modulated in two ways: a secular increase of the gravitational-wave phase due to precession and an oscillation around this secular increase. We show that this secular evolution occurs in precisely three ways, corresponding to physically different apparent evolutions of the binary's precession about the line of sight. We estimate the best possible fitting factor between \emph{any} non-precessing template model and a single precessing signal, in the limit of a constant precession cone. Our closed form estimate of the fitting-factor depends only the geometry of the in-band precession cone; it does not depend explicitly on binary parameters, detector response, or details of either signal model. The precessing black hole--neutron star waveforms least accurately matched by nonspinning waveforms correspond to viewing geometries where the precession cone sweeps the orbital plane repeatedly across the line of sight, in an unfavorable polarization alignment.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Blindly Detecting Merging Supermassive Black Holes with Radio Surveys

Supermassive black holes presumably grow through numerous mergers throughout cosmic time. During each merger, supermassive black hole binaries are surrounded by a circumbinary accretion disk that imposes a significant (~1e4 G for a binary of 1e8 Msun) magnetic field. The motion of the binary through that field will convert the field energy to Poynting flux, with a luminosity ~1e43 erg/s (B/1e4 G)^2 (M/1e8 Msun)^2, some of which may emerge as synchrotron emission at frequencies near 1 GHz where current and planned wide-field radio surveys will operate. We find that the short timescales of many mergers will limit their detectability with most planned blind surveys to <1 per year over the whole sky, independent of the details of the emission process and flux distribution. Including an optimistic estimate for the radio flux makes detection even less likely, with <1 mergers per year over the whole sky. However, wide-field radio instruments may be able to localize systems identified in advance of merger by gravitational waves. Further, radio surveys may be able to detect the weaker emission produced by the binary's motion as it is modulated by spin-orbit precession and inspiral well in advance of merger.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Population Synthesis of Hot Subdwarfs: A Parameter Study

Binaries that contain a hot subdwarf (sdB) star and a main sequence companion may have interacted in the past. This binary population has historically helped determine our understanding of binary stellar evolution. We have computed a grid of binary population synthesis models using different assumptions about the minimum core mass for helium ignition, the envelope binding energy, the common envelope ejection efficiency, the amount of mass and angular momentum lost during stable mass transfer, and the criteria for stable mass transfer on the red giant branch and in the Hertzsprung gap. These parameters separately and together can significantly change the entire predicted population of sdBs. Nonetheless, several different parameter sets can reproduce the observed subpopulation of sdB + white dwarf and sdB + M dwarf binaries, which has been used to constrain these parameters in previous studies. The period distribution of sdB + early F dwarf binaries offers a better test of different mass transfer scenarios for stars that fill their Roche lobes on the red giant branch.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Compact Binary Coalescences in the Band of Ground-based Gravitational-Wave Detectors

As the ground-based gravitational-wave telescopes LIGO, Virgo, and GEO 600 approach the era of first detections, we review the current knowledge of the coalescence rates and the mass and spin distributions of merging neutron-star and black-hole binaries. We emphasize the bi-directional connection between gravitational-wave astronomy and conventional astrophysics. Astrophysical input will make possible informed decisions about optimal detector configurations and search techniques. Meanwhile, rate upper limits, detected merger rates, and the distribution of masses and spins measured by gravitational-wave searches will constrain astrophysical parameters through comparisons with astrophysical models. Future developments necessary to the success of gravitational-wave astronomy are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Intrinsic selection biases of ground-based gravitational wave searches for high-mass BH-BH mergers

The next generation of ground-based gravitational wave detectors may detect a few mergers of comparable-mass M\simeq 100-1000 Msun ("intermediate-mass'', or IMBH) spinning black holes. Black hole spin is known to have a significant impact on the orbit, merger signal, and post-merger ringdown of any binary with non-negligible spin. In particular, the detection volume for spinning binaries depends significantly on the component black hole spins. We provide a fit to the single-detector and isotropic-network detection volume versus (total) mass and arbitrary spin for equal-mass binaries. Our analysis assumes matched filtering to all significant available waveform power (up to l=6 available for fitting, but only l<= 4 significant) estimated by an array of 64 numerical simulations with component spins as large as S_{1,2}/M^2 <= 0.8. We provide a spin-dependent estimate of our uncertainty, up to S_{1,2}/M^2 <= 1. For the initial (advanced) LIGO detector, our fits are reliable for $M\in[100,500]M_\odot$ ($M\in[100,1600]M_\odot$). In the online version of this article, we also provide fits assuming incomplete information, such as the neglect of higher-order harmonics. We briefly discuss how a strong selection bias towards aligned spins influences the interpretation of future gravitational wave detections of IMBH-IMBH mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Pulsar Binary Birthrates with Spin-Opening Angle Correlations

Empirical birthrate estimates for pulsar binaries depend on the fraction of sky subtended by the pulsar beam: the pulsar beaming fraction. This fraction depends on both the pulsar's opening angle and the misalignment angle between its spin and magnetic axes. Previous estimates use the average value for only two pulsars, i.e. PSRs B1913+16 and B1534+12. We explore how birthrate predictions depend on assumptions about opening angle and alignment, using empirically-motivated distributions to define an effective beaming correction factor, f_{b,eff}. For most known pulsars, we expect f_{b,eff} to be less than 6. We also calculate f_{b,eff} for PSRs J0737-3039A and J1141-6545, applying the currently available constraints for their beam geometry. Our median posterior birthrate predictions for tight PSR-NS binaries, wide PSR-NS binaries, and tight PSR-WD binaries are 89/Myr, 0.84/Myr, and 34/Myr, respectively. For pulsars with spin period between 10 ms and 100 ms, we marginalized our posterior birthrate distribution P(R) over a range of plausible beaming correction factors. Our predictions are robust, changing little between several alternative misalignment angle distribution models. [ABRIDGED]

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Unraveling Binary Evolution from Gravitational-Wave Signals and Source Statistics

The next generation of ground-based gravitational-wave detectors are likely to observe gravitational waves from the coalescences of compact-objects binaries. We describe the state of the art for predictions of the rate of compact-binary coalescences and report on initial efforts to develop a framework for converting gravitational-wave observations into improved constraints on astrophysical parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Binary Mergers and Growth of Black Holes in Dense Star Clusters

We model the dynamical evolution of primordial black holes (BHs) in dense star clusters using a simplified treatment of stellar dynamics in which the BHs are assumed to remain concentrated in an inner core, completely decoupled from the background stars. Dynamical interactions involving BH binaries are computed exactly and are generated according to a Monte Carlo prescription. Recoil and ejections lead to complete evaporation of the BH core on a timescale ~10^9 yr for typical globular cluster parameters. Orbital decay driven by gravitational radiation can make binaries merge and, in some cases, successive mergers can lead to significant BH growth. Our highly simplified treatment of the cluster dynamics allows us to study a large number of models and to compute statistical distributions of outcomes, such as the probability of massive BH growth and retention in a cluster. We find that, in most models, there is a significant probability (~20-80%) of BH growth with final masses > 100 M_{\sun}. In one case, a final BH formed with mass ~ 620 M_{\sun}. However, if the typical merger recoil speed (due to asymmetric emission of gravitational radiation) significantly exceeds the cluster escape speed, no growth ever occurs. Independent of the recoil speed, we find that BH-BH mergers enhanced by dynamical interactions in cluster cores present an important source of gravitational waves for ground--based laser interferometers. Under optimistic conditions, the total rate of detections by Advanced LIGO, could be as high as a few tens of events per year from inspiraling BHs from clusters.

0 citations0 reviewsNo code linkedOpen access