Source author record

Ben Farr

Ben Farr 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

12works
5topics
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

12 published item(s)

preprint2022arXiv

Ain't No Mountain High Enough: Semi-Parametric Modeling of LIGO-Virgos Binary Black Hole Mass Distribution

We introduce a semi-parametric model for the primary mass distribution of binary black holes (BBHs) observed with gravitational waves (GWs) that applies a cubic-spline perturbation to a power law. We apply this model to the 46 BBHs included in the second gravitational wave transient catalog (GWTC-2). The spline perturbation model recovers a consistent primary mass distribution with previous results, corroborating the existence of a peak at $35\,M_\odot$ ($>97\%$ credibility) found with the \textsc{Powerlaw+Peak} model. The peak could be the result pulsational pair-instability supernovae (PPISNe). The spline perturbation model finds potential signs of additional features in the primary mass distribution at lower masses similar to those previously reported by Tiwari and Fairhurst (2021). However, with fluctuations due to small number statistics, the simpler \textsc{Powerlaw+Peak} and \textsc{BrokenPowerlaw} models are both still perfectly consistent with observations. Our semi-parametric approach serves as a way to bridge the gap between parametric and non-parametric models to more accurately measure the BBH mass distribution. With larger catalogs we will be able to use this model to resolve possible additional features that could be used to perform cosmological measurements, and will build on our understanding of BBH formation, stellar evolution and nuclear astrophysics.

preprint2022arXiv

Constraining Unmodeled Physics with Compact Binary Mergers from GWTC-1

We present a flexible model to describe the effects of generic deviations of observed gravitational wave signals from modeled waveforms in the LIGO and Virgo gravitational wave detectors. With the detection of 11 gravitational wave events from the GWTC-1 catalog, we are able to constrain possible deviations from our modeled waveforms. In this paper we present our coherent spline model that describes the deviations, then choose to validate our model on an example phenomenological and astrophysically motivated departure in waveforms based on extreme spontaneous scalarization. We find that the model is capable of recovering the simulated deviations. By performing model comparisons we observe that the spline model effectively describes the simulated departures better than a normal compact binary coalescence (CBC) model. We analyze the entire GWTC-1 catalog of events with our model and compare it to a normal CBC model, finding that there are no significant departures from the modeled template gravitational waveforms used.

preprint2022arXiv

Poking Holes: Looking for Gaps in LIGO/Virgo's Black Hole Population

Stellar evolution models predict the existence of a gap in the black hole mass spectrum from $\sim55 M_\odot - 120 M_\odot$ due to pair-instability supernovae (PISNe). We investigate the possible existence of such an "upper" mass gap in the second gravitational wave transient catalog (GWTC-2) by hierarchically modeling the astrophysical distribution of black hole masses. We extend the Truncated and Powerlaw+Peak mass distribution families to allow for an explicit gap in the mass distribution, and apply the extended models to GWTC-2. We find that with the Truncated model there is mild evidence favoring an upper mass gap with log Bayes Factor $\ln \mathcal{B} = 2.79$, inferring the lower and upper bounds at $56.12_{-4.38}^{+7.54} M_\odot$, and $103.74_{-6.32}^{+17.01} M_\odot$ respectively. When using the Powerlaw+Peak model, we find no preference for the gap. When imposing tighter priors on the gap bounds centered on the expected PISN gap bounds, the log Bayes factors in favor of a gap mildly increase. These results are however contingent on the parameter inference for the most massive binary, GW190521, for which follow up analyses showed the source may be an intermediate mass ratio merger that has component masses straddling the gap. Using the GW190521 posterior samples from the analysis in Nitz & Capano (2021), we find an increase in Bayes factors in favor of the gap. However, the overall conclusions are unchanged: There is no preference for a gap when using the Powerlaw+Peak model. This work paves the way for constraining the physics of pair-instability and pulsational pair-instability supernovae and high-mass black hole formation.

preprint2022arXiv

SOAR/Goodman Spectroscopic Assessment of Candidate Counterparts of the LIGO-Virgo Event GW190814

On 2019 August 14 at 21:10:39 UTC, the LIGO/Virgo Collaboration (LVC) detected a possible neutron star-black hole merger (NSBH), the first ever identified. An extensive search for an optical counterpart of this event, designated GW190814, was undertaken using the Dark Energy Camera (DECam) on the 4m Victor M. Blanco Telescope at the Cerro Tololo Inter-American Observatory. Target of Opportunity interrupts were issued on 8 separate nights to observe 11 candidates using the 4.1m Southern Astrophysical Research (SOAR) telescope's Goodman High Throughput Spectrograph in order to assess whether any of these transients was likely to be an optical counterpart of the possible NSBH merger. Here, we describe the process of observing with SOAR, the analysis of our spectra, our spectroscopic typing methodology, and our resultant conclusion that none of the candidates corresponded to the gravitational wave merger event but were all instead other transients. Finally, we describe the lessons learned from this effort. Application of these lessons will be critical for a successful community spectroscopic follow-up program for LVC observing run 4 (O4) and beyond.

preprint2022arXiv

Transient glitch mitigation in Advanced LIGO data with $\textit{glitschen}$

"Glitches" -- transient noise artifacts in the data collected by gravitational wave interferometers like LIGO and Virgo -- are an ever-present obstacle for the search and characterization of gravitational wave signals. With some having morphology similar to high mass, high mass-ratio, and extreme-spin binary black hole events, they limit sensitivity to such sources. They can also act as a contaminant for all sources, requiring targeted mitigation before astrophysical inferences can be made. We propose a data driven, parametric model for frequently encountered glitch types using probabilistic principal component analysis. As a noise analog of parameterized gravitational wave signal models, it can be easily incorporated into existing search and detector characterization techniques. We have implemented our approach with the open source glitschen package. Using LIGO's currently most problematic glitch types, the 'blip' and 'tomte', we demonstrate that parametric models of modest dimension can be constructed and used for effective mitigation in both frequentist and Bayesian analyses.

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.

preprint2016arXiv

Parameter estimation on gravitational waves from neutron-star binaries with spinning components

Inspiraling binary neutron stars are expected to be one of the most significant sources of gravitational-wave signals for the new generation of advanced ground-based detectors. We investigate how well we could hope to measure properties of these binaries using the Advanced LIGO detectors, which began operation in September 2015. We study an astrophysically motivated population of sources (binary components with masses $1.2~\mathrm{M}_\odot$--$1.6~\mathrm{M}_\odot$ and spins of less than $0.05$) using the full LIGO analysis pipeline. While this simulated population covers the observed range of potential binary neutron-star sources, we do not exclude the possibility of sources with parameters outside these ranges; given the existing uncertainty in distributions of mass and spin, it is critical that analyses account for the full range of possible mass and spin configurations. We find that conservative prior assumptions on neutron-star mass and spin lead to average fractional uncertainties in component masses of $\sim 16\%$, with little constraint on spins (the median $90\%$ upper limit on the spin of the more massive component is $\sim 0.7$). Stronger prior constraints on neutron-star spins can further constrain mass estimates, but only marginally. However, we find that the sky position and luminosity distance for these sources are not influenced by the inclusion of spin; therefore, if LIGO detects a low-spin population of BNS sources, less computationally expensive results calculated neglecting spin will be sufficient for guiding electromagnetic follow-up.

preprint2016arXiv

Systematic errors in low latency gravitational wave parameter estimation impact electromagnetic follow-up observations

Among the most eagerly anticipated opportunities made possible by Advanced LIGO/Virgo are multimessenger observations of compact mergers. Optical counterparts may be short-lived so rapid characterization of gravitational wave (GW) events is paramount for discovering electromagnetic signatures. One way to meet the demand for rapid GW parameter estimation is to trade off accuracy for speed, using waveform models with simplified treatment of the compact objects' spin. We report on the systematic errors in GW parameter estimation suffered when using different spin approximations to recover generic signals. Component mass measurements can be biased by $>5σ$ using simple-precession waveforms and in excess of $20σ$ when non-spinning templates are employed. This suggests that electromagnetic observing campaigns should not take a strict approach to selecting which LIGO/Virgo candidates warrant follow-up observations based on low-latency mass estimates. For sky localization, we find searched areas are up to a factor of ${\sim}$2 larger for non-spinning analyses, and are systematically larger for any of the simplified waveforms considered in our analysis. Distance biases for the non-precessing waveforms can be in excess of 100\% and are largest when the spin angular momenta are in the orbital plane of the binary. We confirm that spin-aligned waveforms should be used for low-latency parameter estimation at the minimum. Including simple precession, though more computationally costly, mitigates biases except for signals with extreme precession effects. Our results shine a spotlight on the critical need for development of computationally inexpensive precessing waveforms and/or massively parallel algorithms for parameter estimation.

preprint2015arXiv

Parameter estimation for binary neutron-star coalescences with realistic noise during the Advanced LIGO era

Advanced ground-based gravitational-wave (GW) detectors begin operation imminently. Their intended goal is not only to make the first direct detection of GWs, but also to make inferences about the source systems. Binary neutron-star mergers are among the most promising sources. We investigate the performance of the parameter-estimation \edit{(PE)} pipeline that will be used during the first observing run of the Advanced Laser Interferometer Gravitational-wave Observatory (aLIGO) in 2015: we concentrate on the ability to reconstruct the source location on the sky, but also consider the ability to measure masses and the distance. Accurate, rapid sky-localization is necessary to alert electromagnetic (EM) observatories so that they can perform follow-up searches for counterpart transient events. We consider PE accuracy in the presence of \edit{non-stationary}, non-Gaussian noise. We find that the character of the noise makes negligible difference to the PE performance \edit{at a given signal-to-noise ratio}. The source luminosity distance can only be poorly constrained, the median $90\%$ ($50\%$) credible interval scaled with respect to the true distance is $0.85$ ($0.38$). However, the chirp mass is well measured. Our chirp-mass estimates are subject to systematic error because we used gravitational-waveform templates without component spin to carry out inference on signals with moderate spins, but the total error is typically less than $10^{-3} M_\odot$. The median $90\%$ ($50\%$) credible region for sky localization is $\sim600~\mathrm{deg^{2}}$ ($\sim150~\mathrm{deg^{2}}$), with $3\%$ ($30\%$) of detected events localized within $100~\mathrm{deg^{2}}$. Early aLIGO, with only two detectors, will have a sky-localization accuracy for binary neutron stars of hundreds of square degrees; this makes EM follow-up challenging, but not impossible.

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.

preprint2014arXiv

The First Two Years of Electromagnetic Follow-Up with Advanced LIGO and Virgo

We anticipate the first direct detections of gravitational waves (GWs) with Advanced LIGO and Virgo later this decade. Though this groundbreaking technical achievement will be its own reward, a still greater prize could be observations of compact binary mergers in both gravitational and electromagnetic channels simultaneously. During Advanced LIGO and Virgo's first two years of operation, 2015 through 2016, we expect the global GW detector array to improve in sensitivity and livetime and expand from two to three detectors. We model the detection rate and the sky localization accuracy for binary neutron star (BNS) mergers across this transition. We have analyzed a large, astrophysically motivated source population using real-time detection and sky localization codes and higher-latency parameter estimation codes that have been expressly built for operation in the Advanced LIGO/Virgo era. We show that for most BNS events the rapid sky localization, available about a minute after a detection, is as accurate as the full parameter estimation. We demonstrate that Advanced Virgo will play an important role in sky localization, even though it is anticipated to come online with only one-third as much sensitivity as the Advanced LIGO detectors. We find that the median 90% confidence region shrinks from ~500 square degrees in 2015 to ~200 square degrees in 2016. A few distinct scenarios for the first LIGO/Virgo detections emerge from our simulations.

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.