BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Richard O'Shaughnessy

Richard O'Shaughnessy contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
gr-qcastro-ph.HEastro-ph.COnucl-thastro-ph.IMastro-ph.SR

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
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

16 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&#39;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&#39;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 &#34;intrinsic scattering&#34; 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&#39;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&#39;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
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