Source author record

Duncan A. Brown

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

30 published item(s)

preprint2022arXiv

Advancing the Landscape of Multimessenger Science in the Next Decade

The last decade has brought about a profound transformation in multimessenger science. Ten years ago, facilities had been built or were under construction that would eventually discover the nature of objects in our universe could be detected through multiple messengers. Nonetheless, multimessenger science was hardly more than a dream. The rewards for our foresight were finally realized through IceCube's discovery of the diffuse astrophysical neutrino flux, the first observation of gravitational waves by LIGO, and the first joint detections in gravitational waves and photons and in neutrinos and photons. Today we live in the dawn of the multimessenger era. The successes of the multimessenger campaigns of the last decade have pushed multimessenger science to the forefront of priority science areas in both the particle physics and the astrophysics communities. Multimessenger science provides new methods of testing fundamental theories about the nature of matter and energy, particularly in conditions that are not reproducible on Earth. This white paper will present the science and facilities that will provide opportunities for the particle physics community renew its commitment and maintain its leadership in multimessenger science.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Future Gravitational-Wave Detector Facilities

The next generation of gravitational-wave observatories can explore a wide range of fundamental physics phenomena throughout the history of the universe. These phenomena include access to the universe's binary black hole population throughout cosmic time, to the universe's expansion history independent of the cosmic distance ladders, to stochastic gravitational-waves from early-universe phase transitions, to warped space-time in the strong-field and high-velocity limit, to the equation of state of nuclear matter at neutron star and post-merger densities, and to dark matter candidates through their interaction in extreme astrophysical environments or their interaction with the detector itself. We present the gravitational-wave detector concepts than can drive the future of gravitational-wave astrophysics. We summarize the status of the necessary technology, and the research needed to be able to build these observatories in the 2030s.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Observational Facilities to Study Dark Matter

We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measurement of gravitational waves. Such facilities are a critical component of a multi-pronged experimental program to uncover the nature of dark matter, while often providing complementary measurements of dark energy, neutrino physics, and inflation.

preprint2021arXiv

Reproducing GW150914: the first observation of gravitational waves from a binary black hole merger

In 2016, LIGO and Virgo announced the first observation of gravitational waves from a binary black hole merger, known as GW150914. To establish the confidence of this detection, large-scale scientific workflows were used to measure the event's statistical significance. They used code written by the LIGO/Virgo and were executed on the LIGO Data Grid. The codes are publicly available, but there has not yet been an attempt to directly reproduce the results, although several analyses have replicated the analysis, confirming the detection. We attempt to reproduce the result presented in the GW150914 discovery paper using publicly available code on the Open Science Grid. We show that we can reproduce the main result but we cannot exactly reproduce the LIGO analysis as the original data set used is not public. We discuss the challenges we encountered and make recommendations for scientists who wish to make their work reproducible.

preprint2020arXiv

Constraints on nonlinear tides due to $p$-$g$ mode coupling from the neutron-star merger GW170817

It has been suggested by Weinberg et al. (2013) that an instability due to the nonlinear coupling of a neutron star's tide to its $p$- and $g$-modes could affect the gravitational-wave phase evolution of a neutron-star binary. Weinberg (2016) suggests that this instability can turn on as the gravitational-waves pass through the sensitive band of ground-based detectors, although the size of the effect is not known. The discovery of the binary neutron star merger GW170817 provides an opportunity to look for evidence of nonlinear tides from $p$-$g$ mode coupling. We compute Bayesian evidences that compare waveform models that include the $p$-$g$ mode coupling to models that do not. Using the waveform model and priors of Essick et al. (2016) we find that the observation of GW170817 is consistent with $p$-$g$ mode coupling, in agreement with Abbott, B. P., et al. (2019a). We investigate the properties of the model and find that this consistency is due to degeneracy in a large region of the parameter space between the the model that includes nonlinear tides from $p$-$g$ mode coupling and the standard post-Newtonian model that does not. We investigate the consistency of the GW170817 signal with regions of the parameter space where the effect of nonlinear tides is not degenerate with the standard model. Regions of the nonlinear tide parameter space that have a fitting factor of less than 99$\%$ (98.5$\%$) are disfavored by a Bayes factor of 15 (25). We conclude that the consistency of the GW170817 signal with the model of Essick et al. (2016) is due to parameter degeneracy and that regions where nonlinear tides produce a measurable effect are strongly disfavored.

preprint2020arXiv

DQSEGDB: A time-interval database for storing gravitational wave observatory metadata

The Data Quality Segment Database (DQSEGDB) software is a database service, backend API, frontend graphical web interface, and client package used by the Laser Interferometer Gravitational-Wave Observatory (LIGO), Virgo, GEO600 and the Kamioka Gravitational wave detector for storing and accessing metadata describing the status of their detectors. The DQSEGDB has been used in the analysis of all published detections of gravitational waves in the advanced detector era. The DQSEGDB currently stores roughly 600 million metadata entries and responds to roughly 600,000 queries per day with an average response time of 0.223 ms.

preprint2020arXiv

Measuring the eccentricity of GW170817 and GW190425

Two binary neutron star mergers, GW170817 and GW190425, have been detected by Advanced LIGO and Virgo. These signals were detected by matched-filter searches that assume the star's orbit has circularized by the time their gravitational-wave emission is observable. This suggests that their eccentricity is low, but full parameter estimation of their eccentricity has not yet been performed. We use gravitational-wave observations to measure the eccentricity of GW170817 and GW190425. We find that the eccentricity at a gravitational-wave frequency of 10 Hz is $e \leq 0.024$ and $e \leq 0.048$ for GW170817 and GW190425, respectively (90% confidence). This is consistent with the binaries being formed in the field, as such systems are expected to have circularized to $e \leq 10^{-4}$ by the time they reach the LIGO-Virgo band. Our constraint is a factor of two smaller that an estimate based on GW170817 being detected by searches that neglect eccentricity. However, we caution that we find significant prior dependence in our limits, suggesting that there is limited information in the signals. We note that other techniques used to constrain binary neutron star eccentricity without full parameter estimation may miss degeneracies in the waveform, and that for future signals it will be important to perform full parameter estimation with accurate waveform templates.

preprint2020arXiv

Stringent constraints on neutron-star radii from multimessenger observations and nuclear theory

The properties of neutron stars are determined by the nature of the matter that they contain. These properties can be constrained by measurements of the star's size. We obtain stringent constraints on neutron-star radii by combining multimessenger observations of the binary neutron-star merger GW170817 with nuclear theory that best accounts for density-dependent uncertainties in the equation of state. We construct equations of state constrained by chiral effective field theory and marginalize over these using the gravitational-wave observations. Combining this with the electromagnetic observations of the merger remnant that imply the presence of a short-lived hyper-massive neutron star, we find that the radius of a $1.4\,\rm{M}_\odot$ neutron star is $R_{1.4\,\mathrm{M}_\odot} = 11.0^{+0.9}_{-0.6}~{\rm km}$ (90% credible interval). Using this constraint, we show that neutron stars are unlikely to be disrupted in neutron-star black-hole mergers; subsequently, such events will not produce observable electromagnetic emission.

preprint2019arXiv

Detectability of the subdominant mode in a binary black hole ringdown

The ringdown is the late part of the post-merger signature emitted during the coalescence of two black holes and comprises of a superposition of quasi-normal-modes. Within general relativity, because of the no-hair theorems, the frequencies and damping times of these modes are entirely determined by the mass and angular momentum of the final Kerr black hole. A detection of multiple ringdown modes would potentially allow us to test the no-hair theorem from observational data. The parameters which determine whether sub-dominant ringdown modes can be detected are primarily the overall signal-to-noise ratio present in the ringdown signal, and on the amplitude of the subdominant mode with respect to the dominant mode. In this paper, we use Bayesian inference to determine the detectability of a subdominant mode in a set of simulated analytical ringdown signals. Focusing on the design sensitivity of the Advanced LIGO detectors, we systematically vary the signal-to-noise ratio of the ringdown signal, and the mode amplitude ratio in order to determine what kind of signals are promising for performing black hole spectroscopy.

preprint2019arXiv

Search for Eccentric Binary Neutron Star Mergers in the first and second observing runs of Advanced LIGO

We present a search for gravitational waves from merging binary neutron stars which have non-negligible eccentricity as they enter the LIGO observing band. We use the public Advanced LIGO data which covers the period from 2015 through 2017 and contains $\sim164$ days of LIGO-Hanford and LIGO-Livingston coincident observing time. The search was conducted using matched-filtering using the PyCBC toolkit. We find no significant binary neutron star candidates beyond GW170817, which has previously been reported by searches for binaries in circular orbits. We place a 90% upper limit of $\sim1700$ mergers $\textrm{Gpc}^{-3} \textrm{Yr}^{-1}$ for eccentricities $\lesssim 0.43$ at a dominant-mode gravitational-wave frequency of 10 Hz. The absence of a detection with these data is consistent with theoretical predictions of eccentric binary neutron star merger rates. Using our measured rate we estimate the sensitive volume of future gravitational-wave detectors and compare this to theoretical rate predictions. We find that, in the absence of a prior detection, the rate limits set by six months of Cosmic Explorer observations would constrain all current plausible models of eccentric binary neutron star formation.

preprint2016arXiv

Gravitational waveforms for neutron star binaries from binary black hole simulations

Gravitational waves from binary neutron star (BNS) and black hole/neutron star (BHNS) inspirals are primary sources for detection by the Advanced Laser Interferometer Gravitational-Wave Observatory. The tidal forces acting on the neutron stars induce changes in the phase evolution of the gravitational waveform, and these changes can be used to constrain the nuclear equation of state. Current methods of generating BNS and BHNS waveforms rely on either computationally challenging full 3D hydrodynamical simulations or approximate analytic solutions. We introduce a new method for computing inspiral waveforms for BNS/BHNS systems by adding the post-Newtonian (PN) tidal effects to full numerical simulations of binary black holes (BBHs), effectively replacing the nontidal terms in the PN expansion with BBH results. Comparing a waveform generated with this method against a full hydrodynamical simulation of a BNS inspiral yields a phase difference of $<1$ radian over $\sim 15$ orbits. The numerical phase accuracy required of BNS simulations to measure the accuracy of the method we present here is estimated as a function of the tidal deformability parameter $λ$.

preprint2016arXiv

The PyCBC search for gravitational waves from compact binary coalescence

We describe the PyCBC search for gravitational waves from compact-object binary coalescences in advanced gravitational-wave detector data. The search was used in the first Advanced LIGO observing run and unambiguously identified two black hole binary mergers, GW150914 and GW151226. At its core, the PyCBC search performs a matched-filter search for binary merger signals using a bank of gravitational-wave template waveforms. We provide a complete description of the search pipeline including the steps used to mitigate the effects of noise transients in the data, identify candidate events and measure their statistical significance. The analysis is able to measure false-alarm rates as low as one per million years, required for confident detection of signals. Using data from initial LIGO's sixth science run, we show that the new analysis reduces the background noise in the search, giving a 30% increase in sensitive volume for binary neutron star systems over previous searches.

preprint2015arXiv

Accuracy and precision of gravitational-wave models of inspiraling neutron star -- black hole binaries with spin: comparison with numerical relativity in the low-frequency regime

Coalescing binaries of neutron stars (NS) and black holes (BH) are one of the most important sources of gravitational waves for the upcoming network of ground based detectors. Detection and extraction of astrophysical information from gravitational-wave signals requires accurate waveform models. The Effective-One-Body and other phenomenological models interpolate between analytic results and $10-30$ orbit numerical relativity (NR) merger simulations. In this paper we study the accuracy of these models using new NR simulations that span $36-88$ orbits, with mass-ratios and black hole spins $(q,χ_{BH}) = (7, \pm 0.4), (7, \pm 0.6)$, and $(5, -0.9)$. We find that: (i) the recently published SEOBNRv1 and SEOBNRv2 models of the Effective-One-Body family disagree with each other (mismatches of a few percent) for black hole spins $\geq 0.5$ or $\leq -0.3$, with waveform mismatch accumulating during early inspiral; (ii) comparison with numerical waveforms indicate that this disagreement is due to phasing errors of SEOBNRv1, with SEOBNRv2 in good agreement with all of our simulations; (iii) Phenomenological waveforms disagree with SEOBNRv2 over most of the NSBH binary parameter space; (iv) comparison with NR waveforms shows that most of the model's dephasing accumulates near the frequency interval where it switches to a phenomenological phasing prescription; and finally (v) both SEOBNR and post-Newtonian (PN) models are effectual for NSBH systems, but PN waveforms will give a significant bias in parameter recovery. Our results suggest that future gravitational-wave detection searches and parameter estimation efforts targeted at NSBH systems with $q\lesssim 7$ and $χ_\mathrm{BH} \approx [-0.9, +0.6]$ will benefit from using SEOBNRv2 templates. For larger black hole spins and/or binary mass-ratios, we recommend the models be further investigated as suitable NR simulations become available.

preprint2014arXiv

Accuracy of gravitational waveform models for observing neutron-star--black-hole binaries in Advanced LIGO

Gravitational waves radiated by the coalescence of compact-object binaries containing a neutron star and a black hole are one of the most interesting sources for the ground-based gravitational-wave observatories Advanced LIGO and Advanced Virgo. Advanced LIGO will be sensitive to the inspiral of a $1.4\, M_\odot$ neutron star into a $10\,M_\odot$ black hole to a maximum distance of $\sim 900$ Mpc. Achieving this sensitivity and extracting the physics imprinted in observed signals requires accurate modeling of the binary to construct template waveforms. In a NSBH binary, the black hole may have significant angular momentum (spin), which affects the phase evolution of the emitted gravitational waves. We investigate the ability of post-Newtonian (PN) templates to model the gravitational waves emitted during the inspiral phase of NSBH binaries. We restrict the black hole's spin to be aligned with the orbital angular momentum and compare several approximants. We examine restricted amplitude waveforms that are accurate to 3.5PN order in the orbital dynamics and complete to 2.5PN order in the spin dynamics. We also consider PN waveforms with the recently derived 3.5PN spin-orbit and 3PN spin-orbit tail corrections. We compare these approximants to the effective-one-body model. For all these models, large disagreements start at low to moderate black hole spins, particularly for binaries where the spin is anti-aligned with the orbital angular momentum. We show that this divergence begins in the early inspiral at $v \sim 0.2$ for $χ_{BH} \sim 0.4$. PN spin corrections beyond those currently known will be required for optimal detection searches and to measure the parameters of neutron star--black hole binaries. While this complicates searches, the strong dependence of the gravitational-wave signal on the spin dynamics will make it possible to extract significant astrophysical information.

preprint2014arXiv

Implementing a search for aligned-spin neutron star -- black hole systems with advanced ground based gravitational wave detectors

We study the effect of spins on searches for gravitational waves from compact binary coalescences in realistic simulated early advanced LIGO data. We construct a detection pipeline including matched filtering, signal-based vetoes, a coincidence test between different detectors, and an estimate of the rate of background events. We restrict attention to neutron star--black hole (NS-BH) binary systems, and we compare a search using non-spinning templates to one using templates that include spins aligned with the orbital angular momentum. To run the searches we implement the binary inspiral matched-filter computation in PyCBC, a new software toolkit for gravitational-wave data analysis. We find that the inclusion of aligned-spin effects significantly increases the astrophysical reach of the search. Considering astrophysical NS-BH systems with non-precessing black hole spins, for dimensionless spin components along the orbital angular momentum uniformly distributed in $(-1, 1)$, the sensitive volume of the search with aligned-spin templates is increased by $\sim 50\%$ compared to the non-spinning search; for signals with aligned spins uniformly distributed in the range $(0.7,1)$, the increase in sensitive volume is a factor of $\sim 10$.

preprint2014arXiv

Investigating the effect of precession on searches for neutron-star-black-hole binaries with Advanced LIGO

The first direct detection of neutron-star-black-hole binaries will likely be made with gravitational-wave observatories. Advanced LIGO and Advanced Virgo will be able to observe neutron-star-black-hole mergers at a maximum distance of 900Mpc. To acheive this sensitivity, gravitational-wave searches will rely on using a bank of filter waveforms that accurately model the expected gravitational-wave signal. The angular momentum of the black hole is expected to be comparable to the orbital angular momentum. This angular momentum will affect the dynamics of the inspiralling system and alter the phase evolution of the emitted gravitational-wave signal. In addition, if the black hole's angular momentum is not aligned with the orbital angular momentum it will cause the orbital plane of the system to precess. In this work we demonstrate that if the effect of the black hole's angular momentum is neglected in the waveform models used in gravitational-wave searches, the detection rate of $(10+1.4)M_{\odot}$ neutron-star--black-hole systems would be reduced by $33 - 37%$. The error in this measurement is due to uncertainty in the Post-Newtonian approximations that are used to model the gravitational-wave signal of neutron-star-black-hole inspiralling binaries. We describe a new method for creating a bank of filter waveforms where the black hole has non-zero angular momentum, but is aligned with the orbital angular momentum. With this bank we find that the detection rate of $(10+1.4)M_{\odot}$ neutron-star-black-hole systems would be reduced by $26-33%$. Systems that will not be detected are ones where the precession of the orbital plane causes the gravitational-wave signal to match poorly with non-precessing filter waveforms. We identify the regions of parameter space where such systems occur and suggest methods for searching for highly precessing neutron-star-black-hole binaries.

preprint2013arXiv

Discovery and redshift of an optical afterglow in 71 square degrees: iPTF13bxl and GRB 130702A

We report the discovery of the optical afterglow of the gamma-ray burst (GRB) 130702A, identified upon searching 71 square degrees surrounding the Fermi Gamma-ray Burst Monitor (GBM) localization. Discovered and characterized by the intermediate Palomar Transient Factory (iPTF), iPTF13bxl is the first afterglow discovered solely based on a GBM localization. Real-time image subtraction, machine learning, human vetting, and rapid response multi-wavelength follow-up enabled us to quickly narrow a list of 27,004 optical transient candidates to a single afterglow-like source. Detection of a new, fading X-ray source by Swift and a radio counterpart by CARMA and the VLA confirmed the association between iPTF13bxl and GRB 130702A. Spectroscopy with the Magellan and Palomar 200-inch telescopes showed the afterglow to be at a redshift of z=0.145, placing GRB 130702A among the lowest redshift GRBs detected to date. The prompt gamma-ray energy release and afterglow luminosity are intermediate between typical cosmological GRBs and nearby sub-luminous events such as GRB 980425 and GRB 060218. The bright afterglow and emerging supernova offer an opportunity for extensive panchromatic follow-up. Our discovery of iPTF13bxl demonstrates the first observational proof-of-principle for ~10 Fermi-iPTF localizations annually. Furthermore, it represents an important step towards overcoming the challenges inherent in uncovering faint optical counterparts to comparably localized gravitational wave events in the Advanced LIGO and Virgo era.

preprint2013arXiv

Effect of eccentricity on binary neutron star searches in Advanced LIGO

Binary neutron stars (BNSs) are the primary source of gravitational waves for the Laser Interferometer Gravitational-Wave Observatory (LIGO) and its international partners Virgo and KAGRA. Current BNS searches target field binaries whose orbits will have circularized by radiation reaction before their gravitational waves enter the Advanced LIGO sensitive band at 15 Hz. It has been suggested that a population of BNSs may form by n-body interactions near supermassive black holes or in globular clusters and that these systems may have non-negligible eccentricity in the Advanced LIGO band. We show that for BNS systems with total mass of 2.4 (6.0) solar masses, the effect of eccentricity e < 0.02 (0.05) is negligible and a circular search is effectual for these binaries. For eccentricities up to e = 0.4, we investigate the selection bias caused by neglecting eccentricity in BNS searches. If such high eccentricity systems exist, searches that specifically target eccentric binaries will be needed in Advanced LIGO and Virgo.

preprint2013arXiv

Template Banks for Binary black hole searches with Numerical Relativity waveforms

Gravitational waves (GW) from coalescing stellar-mass black hole binaries (BBH) are expected to be detected by the Advanced Laser Interferometer Gravitational-wave Observatory and Advanced Virgo. Detection searches operate by matched-filtering the detector data using a bank of waveform templates. Traditionally, template banks for BBH are constructed from intermediary analytical waveform models which are calibrated against numerical relativity simulations and which can be aluated for any choice of BBH parameters. This paper explores an alternative to the traditional approach, namely the construction of template banks directly from numerical BBH simulations. Using non-spinning BBH systems as an example, we demonstrate which regions of the mass-parameter plane can be covered with existing numerical BBH waveforms. We estimate the required number and required length of BBH simulations to cover the entire non-spinning BBH parameter plane up to mass-ratio 10, thus illustrating that our approach can be used to guide parameter placement of future numerical simulations. We derive error bounds which are independent of analytical waveform models; therefore, our formalism can be used to independently test the accuracy of such waveform models. The resulting template banks are suitable for advanced LIGO searches.

preprint2013arXiv

When can gravitational-wave observations distinguish between black holes and neutron stars?

Gravitational-wave observations of compact binaries have the potential to uncover the distribution of masses and angular momenta of black holes and neutron stars in the universe. The binary components' physical parameters can be inferred from their effect on the phasing of the gravitational-wave signal, but a partial degeneracy between the components' mass ratio and their angular momenta limits our ability to measure the individual component masses. At the typical signal amplitudes expected by the Advanced Laser Interferometer Gravitational-wave Observatory (signal-to-noise ratios between 10 and 20), we show that it will in many cases be difficult to distinguish whether the components are neutron stars or black holes. We identify when the masses of the binary components could be unambiguously measured outside the range of current observations: a system with a chirp mass $\mathcal{M} \le 0.871 $ M$_\odot$ would unambiguously contain the smallest-mass neutron star observed, and a system with $\mathcal{M} \ge 2.786 \Msun$ must contain a black hole. However, additional information would be needed to distinguish between a binary containing two 1.35 M$_\odot$ neutron stars and an exotic neutron-star--black-hole binary. We also identify those configurations that could be unambiguously identified as black-hole binaries, and show how the observation of an electromagnetic counterpart to a neutron-star--black-hole binary could be used to constrain the black-hole spin.

preprint2012arXiv

Accurate modeling of intermediate-mass-ratio inspirals: Exploring the form of the self-force in the intermediate-mass-ratio regime

In this paper we develop a waveform model that accurately reproduces the dynamical evolution of intermediate-mass-ratio inspirals, as predicted by the effective-one-body (EOB) model introduced in [1], and which enables us to shed some light on the form of the self-force for events with mass-ratio 1:6, 1:10 and 1:100. To complement this study, we make use of self-force results in the extreme-mass-ratio regime, and of predictions of the EOB model introduced in [1], to derive a prescription for the shift of the orbital frequency at the innermost stable circular orbit which consistently captures predictions from the extreme, intermediate and comparable mass-ratio regimes.

preprint2012arXiv

Detecting binary neutron star systems with spin in advanced gravitational-wave detectors

The detection of gravitational waves from binary neutron stars is a major goal of the gravitational-wave observatories Advanced LIGO and Advanced Virgo. Previous searches for binary neutron stars with LIGO and Virgo neglected the component stars' angular momentum (spin). We demonstrate that neglecting spin in matched-filter searches causes advanced detectors to lose more than 3% of the possible signal-to-noise ratio for 59% (6%) of sources, assuming that neutron star dimensionless spins, $c\mathbf{J}/GM^2$, are uniformly distributed with magnitudes between 0 and 0.4 (0.05) and that the neutron stars have isotropically distributed spin orientations. We present a new method for constructing template banks for gravitational wave searches for systems with spin. We present a new metric in a parameter space in which the template placement metric is globally flat. This new method can create template banks of signals with non-zero spins that are (anti-)aligned with the orbital angular momentum. We show that this search loses more than 3% of the maximium signal-to-noise for only 9% (0.2%) of BNS sources with dimensionless spins between 0 and 0.4 (0.05) and isotropic spin orientations. Use of this template bank will prevent selection bias in gravitational-wave searches and allow a more accurate exploration of the distribution of spins in binary neutron stars.

preprint2012arXiv

Importance of including small body spin effects in the modelling of intermediate mass-ratio inspirals. II Accurate parameter extraction of strong sources using higher-order spin effects

We improve the numerical kludge waveform model introduced in [1] in two ways. We extend the equations of motion for spinning black hole binaries derived by Saijo et al. [2] using spin-orbit and spin-spin couplings taken from perturbative and post-Newtonian (PN) calculations at the highest order available. We also include first-order conservative self-force corrections for spin-orbit and spin-spin couplings, which are derived by comparison to PN results. We generate the inspiral evolution using fluxes that include the most recent calculations of small body spin corrections, spin-spin and spin-orbit couplings and higher-order fits to solutions of the Teukolsky equation. Using a simplified version of this model in [1], we found that small body spin effects could be measured through gravitational wave observations from intermediate-mass ratio inspirals (IMRIs) with mass ratio eta ~ 0.001, when both binary components are rapidly rotating. In this paper we study in detail how the spin of the small/big body affects parameter measurement using a variety of mass and spin combinations for typical IMRIs sources. We find that for IMRI events of a moderately rotating intermediate mass black hole (IMBH) of ten thousand solar masses, and a rapidly rotating central supermassive black hole (SMBH) of one million solar masses, gravitational wave observations made with LISA at a fixed signal-to-noise ratio (SNR) of 1000 will be able to determine the inspiralling IMBH mass, the central SMBH mass, the SMBH spin magnitude, and the IMBH spin magnitude to within fractional errors of ~0.001, 0.001, 0.0001, and 9%, respectively. LISA can also determine the location of the source in the sky and the SMBH spin orientation to within ~0.0001 steradians. We show that by including conservative corrections up to 2.5PN order, systematic errors no longer dominate over statistical errors for IMRIs with typical SNR ~1000.

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.

preprint2012arXiv

The NINJA-2 catalog of hybrid post-Newtonian/numerical-relativity waveforms for non-precessing black-hole binaries

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search and parameter-estimation algorithms using numerically generated waveforms, and to foster closer collaboration between the numerical relativity and data analysis communities. The first NINJA project used only a small number of injections of short numerical-relativity waveforms, which limited its ability to draw quantitative conclusions. The goal of the NINJA-2 project is to overcome these limitations with long post-Newtonian - numerical relativity hybrid waveforms, large numbers of injections, and the use of real detector data. We report on the submission requirements for the NINJA-2 project and the construction of the waveform catalog. Eight numerical relativity groups have contributed 63 hybrid waveforms consisting of a numerical portion modelling the late inspiral, merger, and ringdown stitched to a post-Newtonian portion modelling the early inspiral. We summarize the techniques used by each group in constructing their submissions. We also report on the procedures used to validate these submissions, including examination in the time and frequency domains and comparisons of waveforms from different groups against each other. These procedures have so far considered only the $(\ell,m)=(2,2)$ mode. Based on these studies we judge that the hybrid waveforms are suitable for NINJA-2 studies. We note some of the plans for these investigations.

preprint2010arXiv

The Effect of Eccentricity on Searches for Gravitational-Waves from Coalescing Compact Binaries in Ground-based Detectors

Inspiralling compact binaries are expected to circularize before their gravitational-wave signals reach the sensitive frequency band of ground-based detectors. Current searches for gravitational waves from compact binaries using the LIGO and Virgo detectors therefore use circular templates to construct matched filters. Binary formation models have been proposed which suggest that some systems detectable by the LIGO--Virgo network may have non-negligible eccentricity. We investigate the ability of the restricted 3.5 post-Newtonian order TaylorF2 template bank, used by LIGO and Virgo to search for gravitational waves from compact binaries with masses $M \le 35 M_\odot$, to detect binaries with non-zero eccentricity. We model the gravitational waves from eccentric binaries using the $x$-model post-Newtonian formalism proposed by Hinder \emph{et. al.} [I. Hinder, F. Hermann, P. Laguna, and D. Shoemaker, arXiv:0806.1037v1]. We find that small residual eccentricities ($e_0 \lesssim 0.05$ at 40 Hz) do not significantly affect the ability of current LIGO searches to detect gravitational waves from coalescing compact binaries with total mass $2 M_\odot < M < 15 M_\odot$. For eccentricities $e_0 \gtrsim 0.1$, the loss in matched filter signal-to-noise ratio due to eccentricity can be significant and so templates which include eccentric effects will be required to perform optimal searches for such systems.

preprint2009arXiv

Status of NINJA: the Numerical INJection Analysis project

The 2008 NRDA conference introduced the Numerical INJection Analysis project (NINJA), a new collaborative effort between the numerical relativity community and the data analysis community. NINJA focuses on modeling and searching for gravitational wave signatures from the coalescence of binary system of compact objects. We review the scope of this collaboration and the components of the first NINJA project, where numerical relativity groups shared waveforms and data analysis teams applied various techniques to detect them when embedded in colored Gaussian noise.

preprint2009arXiv

Testing gravitational-wave searches with numerical relativity waveforms: Results from the first Numerical INJection Analysis (NINJA) project

The Numerical INJection Analysis (NINJA) project is a collaborative effort between members of the numerical relativity and gravitational-wave data analysis communities. The purpose of NINJA is to study the sensitivity of existing gravitational-wave search algorithms using numerically generated waveforms and to foster closer collaboration between the numerical relativity and data analysis communities. We describe the results of the first NINJA analysis which focused on gravitational waveforms from binary black hole coalescence. Ten numerical relativity groups contributed numerical data which were used to generate a set of gravitational-wave signals. These signals were injected into a simulated data set, designed to mimic the response of the Initial LIGO and Virgo gravitational-wave detectors. Nine groups analysed this data using search and parameter-estimation pipelines. Matched filter algorithms, un-modelled-burst searches and Bayesian parameter-estimation and model-selection algorithms were applied to the data. We report the efficiency of these search methods in detecting the numerical waveforms and measuring their parameters. We describe preliminary comparisons between the different search methods and suggest improvements for future NINJA analyses.

preprint2008arXiv

High-accuracy comparison of numerical relativity simulations with post-Newtonian expansions

Numerical simulations of 15 orbits of an equal-mass binary black hole system are presented. Gravitational waveforms from these simulations, covering more than 30 cycles and ending about 1.5 cycles before merger, are compared with those from quasi-circular zero-spin post-Newtonian (PN) formulae. The cumulative phase uncertainty of these comparisons is about 0.05 radians, dominated by effects arising from the small residual spins of the black holes and the small residual orbital eccentricity in the simulations. Matching numerical results to PN waveforms early in the run yields excellent agreement (within 0.05 radians) over the first $\sim 15$ cycles, thus validating the numerical simulation and establishing a regime where PN theory is accurate. In the last 15 cycles to merger, however, {\em generic} time-domain Taylor approximants build up phase differences of several radians. But, apparently by coincidence, one specific post-Newtonian approximant, TaylorT4 at 3.5PN order, agrees much better with the numerical simulations, with accumulated phase differences of less than 0.05 radians over the 30-cycle waveform. Gravitational-wave amplitude comparisons are also done between numerical simulations and post-Newtonian, and the agreement depends on the post-Newtonian order of the amplitude expansion: the amplitude difference is about 6--7% for zeroth order and becomes smaller for increasing order. A newly derived 3.0PN amplitude correction improves agreement significantly ($<1%$ amplitude difference throughout most of the run, increasing to 4% near merger) over the previously known 2.5PN amplitude terms.

preprint2008arXiv

Rates and Characteristics of Intermediate Mass Ratio Inspirals Detectable by Advanced LIGO

Gravitational waves (GWs) from the inspiral of a neutron star (NS) or stellar-mass black hole (BH) into an intermediate-mass black hole (IMBH) with mass between ~50 and ~350 solar masses may be detectable by the planned advanced generation of ground-based GW interferometers. Such intermediate mass ratio inspirals (IMRIs) are most likely to be found in globular clusters. We analyze four possible IMRI formation mechanisms: (1) hardening of an NS-IMBH or BH-IMBH binary via three-body interactions, (2) hardening via Kozai resonance in a hierarchical triple system, (3) direct capture, and (4) inspiral of a compact object from a tidally captured main-sequence star; we also discuss tidal effects when the inspiraling object is an NS. For each mechanism we predict the typical eccentricities of the resulting IMRIs. We find that IMRIs will have largely circularized by the time they enter the sensitivity band of ground-based detectors. Hardening of a binary via three-body interactions, which is likely to be the dominant mechanism for IMRI formation, yields eccentricities under 10^-4 when the GW frequency reaches 10 Hz. Even among IMRIs formed via direct captures, which can have the highest eccentricities, around 90% will circularize to eccentricities under 0.1 before the GW frequency reaches 10 Hz. We estimate the rate of IMRI coalescences in globular clusters and the sensitivity of a network of three Advanced LIGO detectors to the resulting GWs. We show that this detector network may see up to tens of IMRIs per year, although rates of one to a few per year may be more plausible. We also estimate the loss in signal-to-noise ratio that will result from using circular IMRI templates for data analysis and find that, for the eccentricities we expect, this loss is negligible.