Researcher profile

Will M. Farr

Will M. Farr contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

13 published item(s)

preprint2022arXiv

A Standard Siren Cosmological Measurement from the Potential GW190521 Electromagnetic Counterpart ZTF19abanrhr

The identification of the electromagnetic counterpart candidate ZTF19abanrhr to the binary black hole merger GW190521 opens the possibility to infer cosmological parameters from this standard siren with a uniquely identified host galaxy. The distant merger allows for cosmological inference beyond the Hubble constant. Here we show that the three-dimensional spatial location of ZTF19abanrhr calculated from the electromagnetic data remains consistent with the updated sky localization of GW190521 provided by the LIGO-Virgo Collaboration. If ZTF19abanrhr is associated with the GW190521 merger and assuming a flat wCDM model we find that $H_0 =48^{+24}_{-10}$ km/s/Mpc, $Ω_m =0.39^{+0.38}_{-0.29}$, and $w_0 = -1.29^{+0.63}_{-0.50}$ (median and 68% credible interval). If we use the Hubble constant value inferred from another gravitational-wave event, GW170817, as a prior for our analysis, together with assumption of a flat $Λ$CDM and the model-independent constraint on the physical matter density $ω_m$ from Planck, we find $H_0 = 69.1^{8.7}_{-6.0}$ km/s/Mpc.

preprint2022arXiv

Comparing Bayes factors and hierarchical inference for testing general relativity with gravitational waves

In the context of testing general relativity with gravitational waves, constraints obtained with multiple events are typically combined either through a hierarchical formalism or though a combined multiplicative Bayes factor. We show that the well-known dependence of Bayes factors on the analysis priors in regions of the parameter space without likelihood support can lead to strong confidence in favor of incorrect conclusions when one employs the multiplicative Bayes factor. Bayes factors $\mathcal{O}(1)$ are ambivalent as they depend sensitively on the analysis priors, which are rarely set in a principled way; additionally, combined Bayes factors $>\mathcal{O}(10^3)$ can be obtained in favor of the incorrect conclusion depending on the analysis priors when many $\mathcal{O}(1)$ Bayes factors are multiplied, and specifically when the priors are much wider than the underlying population. The hierarchical analysis that instead infers the ensemble distribution of the individual beyond-general-relativity constraints does not suffer from this problem, and generically converges to favor the correct conclusion. Rather than a naive multiplication, a more reliable Bayes factor can be computed from the hierarchical analysis. We present a number of toy models showing that the practice of multiplying Bayes Factors can lead to incorrect conclusions.

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

Evidence of large recoil velocity from a black hole merger signal

The final black hole left behind after a binary black hole merger can attain a recoil velocity, or a "kick", reaching values up to 5000 km/s. This phenomenon has important implications for gravitational wave astronomy, black hole formation scenarios, testing general relativity, and galaxy evolution. We consider the gravitational wave signal from the binary black hole merger GW200129_065458 (henceforth referred to as GW200129), which has been shown to exhibit strong evidence of orbital precession. Using numerical relativity surrogate models, we constrain the kick velocity of GW200129 to $v_f \sim 1542^{+747}_{-1098}$ km/s or $v_f \gtrsim 698$ km/s (one-sided limit), at 90\% credibility. This marks the first identification of a large kick velocity for an individual gravitational wave event. Given the kick velocity of GW200129, we estimate that there is a less than $0.48\%$ ($7.7\%$) probability that the remnant black hole after the merger would be retained by globular (nuclear star) clusters. Finally, we show that kick effects are not expected to cause biases in ringdown tests of general relativity for this event, although this may change in the future with improved detectors.

preprint2022arXiv

Hiding Planets Near and Far: The Parameter Space of Hidden Companions for Known Planetary Systems

Recent ground and space-based observations show that stars with multiple planets are common in the galaxy. Most of these observational methods are biased toward detecting large planets near to their host stars. Because of these observational biases, these systems can hide small, close-in planets or far-orbiting (big or small) companions. These planets can still exert dynamical influence on known planets and have such influence exerted upon them in turn. In certain configurations, this influence can destabilize the system; in others, the star's gravitational influence can instead further stabilize the system. For example, in systems with planets close to the host star, effects arising from general relativity can help to stabilize the configuration. We derive criteria for hidden planets orbiting both beyond and within known planets that quantify how strongly general relativistic effects can stabilize systems that would otherwise be unstable. As a proof-of-concept, we investigate the several planets in a system based on Kepler 56, and show that the outermost planet will not disrupt the system even at high eccentricities, and show that an Earth-radius planet could be stable within this system if it orbits below 0.08 au. Furthermore, we provide specific predictions to known observed systems by constraining the parameter space of possible hidden planets.

preprint2022arXiv

Hints of spin-orbit resonances in the binary black hole population

Binary black hole spin measurements from gravitational wave observations can reveal the binary's evolutionary history. In particular, the spin orientations of the component black holes within the orbital plane, $ϕ_1$ and $ϕ_2$, can be used to identify binaries caught in the so-called spin-orbit resonances. In a companion paper, we demonstrate that $ϕ_1$ and $ϕ_2$ are best measured near the merger of the two black holes. In this work, we use these spin measurements to provide the first constraints on the full six-dimensional spin distribution of merging binary black holes. In particular, we find that there is a preference for $Δϕ= ϕ_1 - ϕ_2 \sim \pm π$ in the population, which can be a signature of spin-orbit resonances. We also find a preference for $ϕ_1 \sim -π/4$ with respect to the line of separation near merger, which has not been predicted for any astrophysical formation channel. However, the strength of these preferences depends on our prior choices, and we are unable to constrain the widths of the $ϕ_1$ and $Δϕ$ distributions. Therefore, more observations are necessary to confirm the features we find. Finally, we derive constraints on the distribution of recoil kicks in the population, and use this to estimate the fraction of merger remnants retained by globular and nuclear star clusters. We make our spin and kick population constraints publicly available.

preprint2022arXiv

Measuring binary black hole orbital-plane spin orientations

Binary black hole spins are among the key observables for gravitational wave astronomy. Among the spin parameters, their orientations within the orbital plane, $ϕ_1$, $ϕ_2$ and $Δϕ=ϕ_1-ϕ_2$, are critical for understanding the prevalence of the spin-orbit resonances and merger recoils in binary black holes. Unfortunately, these angles are particularly hard to measure using current detectors, LIGO and Virgo. Because the spin directions are not constant for precessing binaries, the traditional approach is to measure the spin components at some reference stage in the waveform evolution, typically the point at which the frequency of the detected signal reaches 20 Hz. However, we find that this is a poor choice for the orbital-plane spin angle measurements. Instead, we propose measuring the spins at a fixed dimensionless time or frequency near the merger. This leads to significantly improved measurements for $ϕ_1$ and $ϕ_2$ for several gravitational wave events. Furthermore, using numerical relativity injections, we demonstrate that $Δϕ$ will also be better measured near the merger for louder signals expected in the future. Finally, we show that numerical relativity surrogate models are key for reliably measuring the orbital-plane spin orientations, even at moderate signal-to-noise ratios like $\sim 30-45$.

preprint2021arXiv

Measuring the star formation rate with gravitational waves from binary black holes

A measurement of the history of cosmic star formation is central to understand the origin and evolution of galaxies. The measurement is extremely challenging using electromagnetic radiation: significant modeling is required to convert luminosity to mass, and to properly account for dust attenuation, for example. Here we show how detections of gravitational waves from inspiraling binary black holes made by proposed third-generation detectors can be used to measure the star formation rate (SFR) of massive stars with high precision up to redshifts of ~10. Depending on the time-delay model, the predicted detection rates ranges from ~2310 to ~56,740 per month with the current measurement of local merger rate density. With 30,000 detections, parameters describing the volumetric SFR can be constrained at the few percent level, and the volumetric merger rate can be directly measured to 3% at z ~ 2. Given a parameterized SFR, the characteristic delay time between binary formation and merger can be measured to ~60%.

preprint2020arXiv

A Population-Informed Mass Estimate for Pulsar J0740+6620

Galactic double neutron star systems have a tight mass distribution around $\sim 1.35 M_{\odot}$, but the mass distribution of all known pulsars is broader. Here we reconstruct the Alsing, et al. (2018) bimodal mass distribution of pulsars observed in binary systems, incorporating data from observations of J0740+6620 which were not available at the time of that work. Because J0740+6620 is an outlier in the mass distribution with non-negligible uncertainty in its mass measurement, its mass receives a large correction from the population, becoming $m_{J0740+6620} = 2.03^{+0.10}_{-0.08} \, M_\odot$ (median and 68\% CI). Stochastic samples from our population model, including population-informed pulsar mass estimates, are available at https://github.com/farr/AlsingNSMassReplication .

preprint2020arXiv

A self-consistent method to estimate the rate of compact binary coalescences with a Poisson mixture model

The recently published GWTC-1 - a journal article summarizing the search for gravitational waves (GWs) from coalescing compact binaries in data produced by the LIGO-Virgo network of ground-based detectors during their first and second observing runs - quoted estimates for the rates of binary neutron star, neutron star black hole binary, and binary black hole mergers, as well as assigned probabilities of astrophysical origin for various significant and marginal GW candidate events. In this paper, we delineate the formalism used to compute these rates and probabilities, which assumes that triggers above a low ranking statistic threshold, whether of terrestrial or astrophysical origin, occur as independent Poisson processes. In particular, we include an arbitrary number of astrophysical categories by redistributing, via mass-based template weighting, the foreground probabilities of candidate events, across source classes. We evaluate this formalism on synthetic GW data, and demonstrate that this method works well for the kind of GW signals observed during the first and second observing runs.

preprint2020arXiv

A trend in the effective spin distribution of LIGO binary black holes with mass

Binary black holes (BBHs) detected by gravitational wave (GW) observations could be broadly divided into two formation channels: those formed through field binary evolution and those assembled dynamically in dense stellar systems. Each of these formation channels, and their sub-channels, populate a distinct region in the effective spin-mass ($χ_{\rm eff}-M$) plane. Depending on the branching ratio of different channels, an ensemble of BBHs could show a trend in this plane. Here we fit a mass-dependent distribution for $χ_{\rm eff}$ to the GWTC-1 BBHs from the first and second observing runs of Advanced LIGO and Advanced Virgo. We find a negative correlation between mass and the mean effective spin ($\barχ_{\mathrm{eff}}$), and positive correlation with its dispersion ($σ_{χ_\mathrm{eff}}$) at 75\% and 80\% confidence. This trend is robust against the choice of mass variable, but most pronounced when the mass variable is taken to be the chirp mass of the binary. The result is consistent with significant contributions from both dynamically assembled and field binaries in the GWTC-1 catalog. The upcoming LIGO O3a data release will critically test this interpretation.

preprint2020arXiv

Cosmological Inference using Gravitational Wave Standard Sirens: A Mock Data Challenge

The observation of binary neutron star merger GW170817, along with its optical counterpart, provided the first constraint on the Hubble constant $H_0$ using gravitational wave standard sirens. When no counterpart is identified, a galaxy catalog can be used to provide the necessary redshift information. However, the true host might not be contained in a catalog which is not complete out to the limit of gravitational-wave detectability. These electromagnetic and gravitational-wave selection effects must be accounted for. We describe and implement a method to estimate $H_0$ using both the counterpart and the galaxy catalog standard siren methods. We perform a series of mock data analyses using binary neutron star mergers to confirm our ability to recover an unbiased estimate of $H_0$. Our simulations used a simplified universe with no redshift uncertainties or galaxy clustering, but with different magnitude-limited catalogs and assumed host galaxy properties, to test our treatment of both selection effects. We explore how the incompleteness of catalogs affects the final measurement of $H_0$, as well as the effect of weighting each galaxy's likelihood of being a host by its luminosity. In our most realistic simulation, where the simulated catalog is about three times denser than the density of galaxies in the local universe, we find that a 4.4\% measurement precision can be reached using galaxy catalogs with 50\% completeness and $\sim 250$ binary neutron star detections with sensitivity similar to that of Advanced LIGO's second observing run.

preprint2019arXiv

The Most Massive Binary Black Hole Detections and the Identification of Population Outliers

Advanced LIGO and Virgo detected ten binary black holes (BBHs) in their first two observing runs (O1 and O2). Analysis of these events found strong evidence for a dearth of BBHs with component masses greater than $\sim45 \ M_\odot$, as would be expected from a pair-instability mass gap. Meanwhile, a standalone analysis of the merger GW170729 found its primary mass $m_1 = {51.2^{+16.2}_{-11.0} \ M_\odot}$, with the majority of its posterior support at $m_1 > 45 \ M_\odot$. Although this appears to be in contradiction with the existence of a limit at $\sim45\ M_\odot$, we argue that individual events cannot be evaluated without reference to the entire population. When GW170729 is analyzed jointly with the rest of the detections, as part of a full hierarchical population analysis, its inferred primary mass tightens considerably, to $m_1 = {38.9^{+7.3}_{-4.5} \ M_\odot}$. For a large sample of events in the presence of noise, apparent outliers in the detected distribution are inevitable, even if the underlying population forbids outliers. We discuss methods of distinguishing between statistical fluctuations and population outliers using posterior predictive tests. Applying these tests to the primary mass distribution in O1 and O2, we find that the ten detections are consistent with even the simplest power-law plus maximum-mass model considered by the LVC. This supports the claim that GW170729 is not a population outlier. We also provide non-parametric constraints on the rate of high-mass mergers and conservatively bound the rate of mergers with $m_1 > 45 \ M_\odot$ at $2.8^{+5.4}_{-2.0}\%$ of the total merger rate. After 100 detections like those of O1 and O2 from a population with a maximum primary mass of $45 \, M_\odot$, it would be common for the most massive system to have an observed maximum-likelihood mass $m_1 \gtrsim 70 \, M_\odot$.