Researcher profile

Eric Thrane

Eric Thrane contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
27works
0followers
4topics
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

27 published item(s)

preprint2022arXiv

A Rosetta Stone for eccentric gravitational waveform models

Orbital eccentricity is a key signature of dynamical binary black hole formation. The gravitational waves from a coalescing binary contain information about its orbital eccentricity, which may be measured if the binary retains sufficient eccentricity near merger. Dedicated waveforms are required to measure eccentricity. Several models have been put forward, and show good agreement with numerical relativity at the level of a few percent or better. However, there are multiple ways to define eccentricity for inspiralling systems, and different models internally use different definitions of eccentricity, making it difficult to directly compare eccentricity measurements. In this work, we systematically compare two eccentric waveform models, $\texttt{SEOBNRE}$ and $\texttt{TEOBResumS}$, by developing a framework to translate between different definitions of eccentricity. This mapping is constructed by minimizing the relative mismatch between the two models over eccentricity and reference frequency, before evolving the eccentricity of one model to the same reference frequency as the other model. We show that for a given value of eccentricity passed to $\texttt{SEOBNRE}$, one must input a $20$-$50\%$ smaller value of eccentricity to $\texttt{TEOBResumS}$ in order to obtain a waveform with the same empirical eccentricity. We verify this mapping by repeating our analysis for eccentric numerical relativity simulations, demonstrating that $\texttt{TEOBResumS}$ reports a correspondingly smaller value of eccentricity than $\texttt{SEOBNRE}$.

preprint2022arXiv

Building better spin models for merging binary black holes: Evidence for non-spinning and rapidly spinning nearly aligned sub-populations

Recent work paints a conflicting portrait of the distribution of black hole spins in merging binaries measured with gravitational waves. Some analyses find that a significant fraction of merging binaries contain at least one black hole with a spin tilt $>90^\circ$ with respect to the orbital angular momentum vector, which has been interpreted as a signature for dynamical assembly. Other analyses find the data are consistent with a bimodal population in which some binaries contain black holes with negligible spin while the rest contain black holes with spin vectors preferentially aligned with the orbital angular momentum vector. In this work, we scrutinize models for the distribution of black hole spins to pinpoint possible failure modes in which the model yields a faulty conclusion. We reanalyze data from the second LIGO--Virgo gravitational-wave transient catalog (GWTC-2) using a revised spin model, which allows for a sub-population of black holes with negligible spins. In agreement with recent results by Roulet et al., we show that the GWTC-2 detections are consistent with two distinct sub-populations. We estimate that $29-75\%$ (90\% credible interval) of merging binaries contain black holes with negligible spin $χ\approx 0$. The remaining binaries are part of a second sub-population in which the spin vectors are preferentially (but not exactly) aligned to the orbital angular momentum. The black holes in this second sub-population are characterized by spins of $χ\sim0.45$. We suggest that the inferred spin distribution is consistent with the hypothesis that all merging binaries form via the field formation scenario.

preprint2022arXiv

Consistency of the Parkes Pulsar Timing Array Signal with a Nanohertz Gravitational-wave Background

Pulsar timing array experiments have recently reported strong evidence for a common-spectrum stochastic process with a strain spectral index consistent with that expected of a nanohertz-frequency gravitational-wave background, but with negligible yet non-zero evidence for spatial correlations required for a definitive detection. However, it was pointed out by the Parkes Pulsar Timing Array (PPTA) collaboration that the same models used in recent analyses resulted in strong evidence for a common-spectrum process in simulations where none is present. In this work, we introduce a methodology to distinguish pulsar power spectra with the same amplitude from noise power spectra of similar but distinct amplitudes. The former is the signature of a spatially uncorrelated pulsar term of a nanohertz gravitational-wave background, whereas the latter could represent ensemble pulsar noise properties. We test the methodology on simulated data sets. We find that the reported common process in PPTA pulsars is indeed consistent with the spectral feature of a pulsar term. We recommend this methodology as one of the validity tests that the real astrophysical and cosmological backgrounds should pass, as well as for inferences about the spatially uncorrelated component of the background.

preprint2022arXiv

Fast, flexible, and accurate evaluation of gravitational-wave Malmquist bias with machine learning

Many astronomical surveys are limited by the brightness of the sources, and gravitational-wave searches are no exception. The detectability of gravitational waves from merging binaries is affected by the mass and spin of the constituent compact objects. To perform unbiased inference on the distribution of compact binaries, it is necessary to account for this selection effect, which is known as Malmquist bias. Since systematic error from selection effects grows with the number of events, it will be increasingly important over the coming years to accurately estimate the observational selection function for gravitational-wave astronomy. We employ density estimation methods to accurately and efficiently compute the compact binary coalescence selection function. We introduce a simple pre-processing method, which significantly reduces the complexity of the required machine learning models. We demonstrate that our method has smaller statistical errors at comparable computational cost than the method currently most widely used allowing us to probe narrower distributions of spin magnitudes. The currently used method leaves $10-50\%$ of the interesting black hole spin models inaccessible; our new method can probe $>99\%$ of the models and has a lower uncertainty for $>80\%$ of the models.

preprint2022arXiv

Linking the rates of neutron star binaries and short gamma-ray bursts

Short gamma-ray bursts are believed to be produced by both binary neutron star (BNS) and neutron star-black hole (NSBH) mergers. We use current estimates for the BNS and NSBH merger rates to calculate the fraction of observable short gamma-ray bursts produced through each channel. This allows us to constrain merger rates of BNS to $\mathcal{R}_{\rm{BNS}}=384^{+431}_{-213}{\rm{Gpc}^{-3} \rm{yr}^{-1}}$ ($90\%$ credible interval), a $16\%$ decrease in the rate uncertainties from the second LIGO--Virgo Gravitational-Wave Transient Catalog, GWTC-2. Assuming a top-hat emission profile with a large Lorentz factor, we constrain the average opening angle of gamma-ray burst jets produced in BNS mergers to $\approx 15^\circ$. We also measure the fraction of BNS and NSBH mergers that produce an observable short gamma-ray burst to be $0.02^{+0.02}_{-0.01}$ and $0.01 \pm 0.01$, respectively and find that $\gtrsim 40\%$ of BNS mergers launch jets (90\% confidence). We forecast constraints for future gravitational-wave detections given different modelling assumptions, including the possibility that BNS and NSBH jets are different. With $24$ BNS and $55$ NSBH observations, expected within six months of the LIGO-Virgo-KAGRA network operating at design sensitivity, it will be possible to constrain the fraction of BNS and NSBH mergers that launch jets with $10\%$ precision. Within a year of observations, we can determine whether the jets launched in NSBH mergers have a different structure than those launched in BNS mergers and rule out whether $\gtrsim 80\%$ of binary neutron star mergers launch jets. We discuss the implications of future constraints on understanding the physics of short gamma-ray bursts and binary evolution.

preprint2022arXiv

Orbital dynamics and extreme scattering event properties from long-term scintillation observations of PSR J1603-7202

We model long-term variations in the scintillation of binary pulsar PSR J1603$-$7202, observed by the 64 m Parkes radio telescope (Murriyang) between 2004 and 2016. We find that the time variation in the scintillation arc curvature is well-modelled by scattering from an anisotropic thin screen of plasma between the Earth and the pulsar. Using our scintillation model, we measure the inclination angle and longitude of ascending node of the orbit, yielding a significant improvement over the constraints from pulsar timing. From our measurement of the inclination angle, we place a lower bound on the mass of J1603$-$7202's companion of $\gtrsim 0.5\,\text{M}_\odot$ assuming a pulsar mass of $\gtrsim1.2\,\text{M}_\odot$. We find that the scintillation arcs are most pronounced when the electron column density along the line of sight is increased, and that arcs are present during a known extreme scattering event. We measure the distance to the interstellar plasma and its velocity, and we discuss some structures seen in individual scintillation arcs within the context of our model.

preprint2022arXiv

When models fail: an introduction to posterior predictive checks and model misspecification in gravitational-wave astronomy

Bayesian inference is a powerful tool in gravitational-wave astronomy. It enables us to deduce the properties of merging compact-object binaries and to determine how these mergers are distributed as a population according to mass, spin, and redshift. As key results are increasingly derived using Bayesian inference, there is increasing scrutiny on Bayesian methods. In this review, we discuss the phenomenon of \textit{model misspecification}, in which results obtained with Bayesian inference are misleading because of deficiencies in the assumed model(s). Such deficiencies can impede our inferences of the true parameters describing physical systems. They can also reduce our ability to distinguish the "best fitting" model: it can be misleading to say that Model~A is preferred over Model~B if both models are manifestly poor descriptions of reality. Broadly speaking, there are two ways in which models fail: models that fail to adequately describe the data (either the signal or the noise) have misspecified likelihoods. Population models -- designed, for example, to describe the distribution of black hole masses -- may fail to adequately describe the true population due to a misspecified prior. We recommend tests and checks that are useful for spotting misspecified models using examples inspired by gravitational-wave astronomy. We include companion python notebooks to illustrate essential concepts.

preprint2021arXiv

Black-hole spectroscopy, the no-hair theorem and GW150914: Kerr vs. Occam

The &#34;no-hair&#34; theorem states that astrophysical black holes are fully characterised by just two numbers: their mass and spin. The gravitational-wave emission from a perturbed black-hole consists of a superposition of damped sinusoids, known as \textit{quasi-normal modes}. Quasi-normal modes are specified by three integers $(\ell,m,n)$: the $(\ell, m)$ integers describe the angular properties and $(n)$ specifies the (over)tone. If the no-hair theorem holds, the frequencies and damping times of quasi-normal modes are determined uniquely by the mass and spin of the black hole, while phases and amplitudes depend on the particular perturbation. Current tests of the no-hair theorem, attempt to identify these modes in a semi-agnostic way, without imposing priors on the source of the perturbation. This is usually known as \textit{black-hole spectroscopy}. Applying this framework to GW150914, the measurement of the first overtone led to the confirmation of the theorem to $20\%$ level. We show, however, that such semi-agnostic tests cannot provide strong evidence in favour of the no-hair theorem, even for extremely loud signals, given the increasing number of overtones (and free parameters) needed to fit the data. This can be solved by imposing prior assumptions on the origin of the perturbed black hole that can further constrain the explored parameters: in particular, our knowledge that the ringdown is sourced by a binary black hole merger. Applying this strategy to GW150914 we find a natural log Bayes factor of $\sim 6.5$ in favour of the Kerr nature of its remnant, indicating that the hairy object hypothesis is disfavoured with $<1:600$ with respect to the Kerr black-hole one.

preprint2021arXiv

Probing Extremal Gravitational-Wave Events with Coarse-Grained Likelihoods

As catalogs of gravitational-wave transients grow, new records are set for the most extreme systems observed to date. The most massive observed black holes probe the physics of pair instability supernovae while providing clues about the environments in which binary black hole systems are assembled. The least massive black holes, meanwhile, allow us to investigate the purported neutron star-black hole mass gap, and binaries with unusually asymmetric mass ratios or large spins inform our understanding of binary and stellar evolution. Existing outlier tests generally implement leave-one-out analyses, but these do not account for the fact that the event being left out was by definition an extreme member of the population. This results in a bias in the evaluation of outliers. We correct for this bias by introducing a coarse-graining framework to investigate whether these extremal events are true outliers or whether they are consistent with the rest of the observed population. Our method enables us to study extremal events while testing for population model misspecification. We show that this ameliorates biases present in the leave-one-out analyses commonly used within the gravitational-wave community. Applying our method to results from the second LIGO--Virgo transient catalog, we find qualitative agreement with the conclusions of Abbott et al, ApJL 913 L7 (2021). GW190814 is an outlier because of its small secondary mass. We find that neither GW190412 nor GW190521 are outliers.

preprint2021arXiv

Temperature dependent appearance of exotic matter makes nascent neutron stars spin faster

Neutron stars offer the opportunity to study the behaviour of matter at densities and temperatures inaccessible to terrestrial experiments. Gravitational-wave observations of binary neutron star coalescences can constrain the neutron-star equation of state before and after merger. After the neutron star binary merges, hyperons can form in the remnant, changing the behaviour of the neutron-star equation of state. In this study, we use finite-entropy equations of state to show that a post-merger remnant can spin up due to cooling. The magnitude of the spin-up depends on the neutron-star equation of state. If hyperons are present, the post-merger spin-up changes the peak gravitational-wave frequency by $\sim 540$ Hz, when the entropy per baryon drops from $s=2$ $k_B$ to $s=0$ $k_B$. If hyperons are not present, the post-merger spin-up changes by $\sim 360$ Hz, providing a gravitational-wave signature for exotic matter. We expect the same qualitative behaviour whenever temperature dependent phase transitions are triggered.

preprint2021arXiv

The imprint of superradiance on hierarchical black hole mergers

Ultralight bosons are a proposed solution to outstanding problems in cosmology and particle physics: they provide a dark-matter candidate while potentially explaining the strong charge-parity problem. If they exist, ultralight bosons can interact with black holes through the superradiant instability. In this work we explore the consequences of this instability on the evolution of hierarchical black holes within dense stellar clusters. By reducing the spin of individual black holes, superradiance reduce the recoil velocity of merging binary black holes, which, in turn, increases the retention fraction of hierarchical merger remnants. We show that the existence of ultralight bosons with mass $ 2\times10^{-14}\lesssim μ/\textrm{eV} \lesssim2\times10^{-13}$ would lead to an increased rate of hierarchical black hole mergers in nuclear star clusters. An ultralight boson in this energy range would result in up to $\approx60\%$ more present-day nuclear star clusters supporting hierarchical growth. The presence of an ultralight boson can also double the rate of intermediate mass black hole mergers to $\approx0.08$\,Gpc$^{-3}$\,yr$^{-1}$ in the local Universe. These results imply that a select range of ultralight boson mass can have far-reaching consequences for the population of black holes in dense stellar environments. Future studies into black hole cluster populations and the spin distribution of hierarchically formed black holes will test this scenario.

preprint2020arXiv

Constraining short gamma-ray burst jet properties with gravitational waves and gamma rays

Gamma-ray burst (GRB) prompt emission is highly beamed, and understanding the jet geometry and beaming configuration can provide information on the poorly understood central engine and circum-burst environment. Prior to the advent of gravitational-wave astronomy, astronomers relied on observations of jet breaks in the multi-wavelength afterglow to determine the GRB opening angle, since the observer&#39;s viewing angle relative to the system cannot be determined from the electromagnetic data alone. Gravitational-wave observations, however, provide an independent measurement of the viewing angle. We describe a Bayesian method for determining the geometry of short GRBs using coincident electromagnetic and gravitational-wave observations. We demonstrate how an ensemble of multi-messenger detections can be used to measure the distributions of the jet energy, opening angle, Lorentz factor, and angular profile of short GRBs; we find that for a population of 100 such observations, we can constrain the mean of the opening angle distribution to within 10 degrees regardless of the angular emission profile. Conversely, the constraint on the energy distribution depends on the shape of the profile, which can be distinguished.

preprint2020arXiv

Constraining temperature distribution inside LIGO test masses from frequencies of their vibrational modes

Thermal distortion of test masses, as well as thermal drift of their vibrational mode frequencies, present a major challenge for operation of the Advanced LIGO and Advanced VIRGO interferometers, reducing optical efficiency, which limits sensitivity and potentially causing instabilities which reduce duty-cycle. In this paper, we demonstrate that test-mass vibrational mode frequency data can be used to overcome some of these difficulties. First, we derive a general expression for the change in a mode frequency as a function of temperature distribution inside the test mass. Then we show how the mode frequency dependence on temperature distribution can be used to identify the wavefunction of observed vibrational modes. We then show how monitoring the frequencies of multiple vibrational modes allows the temperature distribution inside the test mass to be strongly constrained. Finally, we demonstrate using simulations, the potential to improve the thermal model of the test mass, providing independent and improved estimates of important parameters such as the coating absorption coefficient and the location of point absorbers.

preprint2020arXiv

Constraining the Gravitational-Wave Afterglow From a Binary Neutron Star Coalescence

Binary neutron star mergers are rich laboratories for physics, accessible with ground-based interferometric gravitational-wave detectors such as the Advanced LIGO and Advanced Virgo. If a neutron star remnant survives the merger, it can emit gravitational waves that might be detectable with the current or next generation detectors. The physics of the long-lived post-merger phase is not well understood and makes modelling difficult. In particular the phase of the gravitational-wave signal is not well modelled. In this paper, we explore methods for using long duration post-merger gravitational-wave signals to constrain the parameters and the properties of the remnant. We develop a phase-agnostic likelihood model that uses only the spectral content for parameter estimation and demonstrate the calculation of a Bayesian upper limit in the absence of a signal. With the millisecond magnetar model, we show that for an event like GW170817, the ellipticity of a long-lived remnant can be constrained to less than about 0.5 in the parameter space used.

preprint2020arXiv

Gravitational-wave astronomy with a physical calibration model

We carry out astrophysical inference for compact binary merger events in LIGO-Virgo&#39;s first gravitational-wave transient catalog (GWTC-1) using a physically motivated calibration model. We demonstrate that importance sampling can be used to reduce the cost of what would otherwise be a computationally challenging analysis. We show that including the physical estimate for the calibration error distribution has negligible impact on the inference of parameters for the events in GWTC-1. Studying a simulated signal with matched filter signal-to-noise ratio $\text{SNR}=200$, we project that a calibration error estimate typical of GWTC-1 is likely to be negligible for the current generation of gravitational-wave detectors. We argue that other sources of systematic error---from waveforms, prior distributions, and noise modelling---are likely to be more important. Finally, using the events in GWTC-1 as standard sirens, we infer an astrophysically-informed improvement on the estimate of the calibration error in the LIGO interferometers.

preprint2020arXiv

Gravitational-wave astronomy with an uncertain noise power spectral density

In order to extract information about the properties of compact binaries, we must estimate the noise power spectral density of gravitational-wave data, which depends on the properties of the gravitational-wave detector. In practice, it is not possible to know this perfectly, only to estimate it from the data. Multiple estimation methods are commonly used and each has a corresponding statistical uncertainty. However, this uncertainty is widely ignored when measuring the physical parameters describing compact binary coalescences, and the appropriate likelihoods which account for the uncertainty are not well known. In order to perform increasingly precise astrophysical inference and model selection, it will be essential to account for this uncertainty. In this work, we derive the correct likelihood for one of the most widely used estimation methods in gravitational-wave transient analysis, the median average. We demonstrate that simulated Gaussian noise follows the predicted distributions. We then examine real gravitational-wave data at and around the time of GW151012, a relatively low-significance binary black hole merger event. We show that the data are well described by stationary-Gaussian noise and explore the impact of different noise power spectral density estimation methods on the astrophysical inferences we draw about GW151012.

preprint2020arXiv

Inferring the population properties of binary black holes from unresolved gravitational waves

The vast majority of compact binary mergers in the Universe produce gravitational waves that are too weak to yield unambiguous detections; they are unresolved. We present a method to infer the population properties of compact binaries -- such as their merger rates, mass spectrum, and spin distribution -- using both resolved and unresolved gravitational waves. By eliminating entirely the distinction between resolved and unresolved signals, we eliminate bias from selection effects. To demonstrate this method, we carry out a Monte Carlo study using an astrophysically motivated population of binary black holes. We show that some population properties of compact binaries are well constrained by unresolved signals after about one week of observation with Advanced LIGO at design sensitivity.

preprint2020arXiv

Is there a spectral turnover in the spin noise of millisecond pulsars?

Pulsar timing arrays provide a unique means to detect nanohertz gravitational waves through long-term measurements of pulse arrival times from an ensemble of millisecond pulsars. After years of observations, some timing array pulsars have been shown to be dominated by low-frequency red noise, including spin noise that might be associated with pulsar rotational irregularities. The power spectral density of pulsar timing red noise is usually modeled with a power law or a power law with a turnover frequency below which the noise power spectrum plateaus. If there is a turnover in the spin noise of millisecond pulsars, residing within the observation band of current and/or future pulsar timing measurements, it may be easier than projected to resolve the gravitational-wave background from supermassive binary black holes. Additionally, the spectral turnover can provide valuable insights on neutron star physics. In the recent study by Melatos and Link, the authors provided a derivation of the model for power spectral density of spin noise from superfluid turbulence in the core of a neutron star, from first principles. The model features a spectral turnover, which depends on the dynamical response time of the superfluid and the steady-state angular velocity lag between the crust and the core of the star. In this work, we search for a spectral turnover in spin noise using the first data release of the International Pulsar Timing Array. Through Bayesian model selection, we find no evidence of a spectral turnover. Our analysis also shows that data from pulsars J1939+2134, J1024$-$0719 and J1713+0747 prefers the power-law model to the superfluid turbulence model.

preprint2020arXiv

Measuring the neutron star equation of state with gravitational waves: the first forty binary neutron star mergers

Gravitational waves from binary neutron star coalescences contain rich information about matter at supranuclear densities encoded by the neutron star equation of state. We can measure the equation of state by analyzing the tidal interactions between neutron stars, which is quantified by the tidal deformability. Multiple merger events are required to probe the equation of state over a range of neutron star masses. The more events included in the analysis, the stronger the constraints on the equation of state. In this paper, we build on previous work to explore the constraints that LIGO and Virgo are likely to place on the neutron star equation of state by combining the first forty binary neutron star detections, a milestone we project to be reached during the first year of accumulated design-sensitivity data. We carry out Bayesian inference on a realistic mock dataset of binaries to obtain posterior distributions for neutron star tidal parameters. In order to combine posterior samples from multiple observations, we employ a random forest regressor, which allows us to efficiently interpolate the likelihood distribution. Assuming a merger rate of 1540 Gpc$^{-3}$ yr$^{-1}$ and a LIGO-Virgo detector network operating for one year at the sensitivity of the third-observation run, plus an additional eight months of design sensitivity, we find that the radius of a 1.4 $M_\odot$ neutron star can be constrained to $\sim 10$% at 90% confidence. At the same time, the pressure at twice the nuclear saturation density can be constrained to $\sim 45$ % at 90% confidence. Finally, we add an appendix following publication of the paper in the journal, showing the posterior distribution of the maximum neutron star mass allowed by the equation of state. We find that the maximum mass can be constrained to $\sim 0.3$ $M_\odot$ at 90% confidence.

preprint2020arXiv

Memory Effect or Cosmic String? Classifying Gravitational-Wave Bursts with Bayesian Inference

A variety of gravitational-wave transient sources can be modeled in the Fourier domain using a power law. This simple power-law model provides a reasonable approximation for gravitational-wave bursts from cosmic string cusps, cosmic string kinks, and the memory effect. Each of these sources is described using a different spectral index. In this work, we simulate interferometer strain data with injections of power-law and memory bursts to demonstrate parameter estimation, signal detection, and model selection. We show how Bayesian inference can be used to measure the power-law spectral index, thereby distinguishing between different astrophysical scenarios.

preprint2020arXiv

On the origin of GW190425

The LIGO/Virgo collaborations recently announced the detection of a likely binary neutron star merger, GW190425. The total mass of GW190425 is significantly larger than the masses of Galactic double neutron stars known through radio astronomy. This suggests that GW190425 formed differently from Galactic double neutron stars. We hypothesize that GW190425 formed via unstable &#34;case BB&#34; mass transfer. According to this hypothesis, the progenitor of GW190425 was a binary consisting of a neutron star and a ${\sim} 4-5 {M_\odot}$ helium star, which underwent a common-envelope process. Following the supernovae of the helium star core, a tight, eccentric, double neutron star was formed, which merged in ${\lesssim}$ 10 Myr. The helium star progenitor may explain the unusually large mass of GW190425, while the short time to merger may explain why we do not see similar systems in radio. In order to test this hypothesis, we measure the eccentricity of GW190425 using publicly available LIGO/Virgo data. We constrain the eccentricity at 10 Hz to be $e \leq 0.007$ with $90\%$ confidence. This result provides no evidence for or against the unstable mass transfer scenario because the binary is likely to have circularized to $e\lesssim10^{-4}$ by the time it entered the LIGO/Virgo band. Future detectors operating in lower frequency bands will enable us to discern the formation channel of mergers similar to GW190425 using eccentricity measurements.

preprint2020arXiv

Standard-siren cosmology using gravitational waves from binary black holes

Gravitational-wave astronomy provides a unique new way to study the expansion history of the Universe. In this work, we investigate the impact future gravitational-wave observatories will have on cosmology. Third-generation observatories like the Einstein Telescope and Cosmic Explorer will be sensitive to essentially all of the binary black hole coalescence events in the Universe. Recent work by \cite{farr2019future} points out that features in the stellar-mass black hole population break the mass-redshift degeneracy, facilitating precise determination of the Hubble parameter without electromagnetic counterparts or host galaxy catalogues. Using a hierarchical Bayesian inference model, we show that with one year of observations by the Einstein Telescope, the Hubble constant will be measured to $\lesssim 1\%$. We also show that this method can be used to perform Bayesian model selection between cosmological models. As an illustrative example, we find that a decisive statement can be made comparing the $Λ$CDM and RHCT cosmological models using two weeks of data from the Einstein Telescope.

preprint2020arXiv

The astrophysical odds of GW151216

The gravitational-wave candidate GW151216 is a proposed binary black hole event from the first observing run of the Advanced LIGO detectors. Not identified as a bona fide signal by the LIGO--Virgo collaboration, there is disagreement as to its authenticity, which is quantified by $p_\text{astro}$, the probability that the event is astrophysical in origin. Previous estimates of $p_\text{astro}$ from different groups range from 0.18 to 0.71, making it unclear whether this event should be included in population analyses, which typically require $p_\text{astro}>0.5$. Whether GW151216 is an astrophysical signal or not has implications for the population properties of stellar-mass black holes and hence the evolution of massive stars. Using the astrophysical odds, a Bayesian method which uses the signal coherence between detectors and a parameterised model of non-astrophysical detector noise, we find that $p_\text{astro}=0.03$, suggesting that GW151216 is unlikely to be a genuine signal. We also analyse GW150914 (the first gravitational-wave detection) and GW151012 (initially considered to be an ambiguous detection) and find $p_\text{astro}$ values of 1 and 0.997 respectively. We argue that the astrophysical odds presented here improve upon traditional methods for distinguishing signals from noise.

preprint2020arXiv

Toward the unambiguous identification of supermassive binary black holes through Bayesian inference

Supermassive binary black holes at sub-parsec orbital separations have yet to be discovered, with the possible exception of blazar OJ~287. In parallel to the global hunt for nanohertz gravitational waves from supermassive binaries using pulsar timing arrays, there has been a growing sample of candidates reported from electromagnetic surveys, particularly searches for periodic variations in optical light curves of quasars. However, the periodicity search is prone to false positives from quasar red noise and quasi-periodic oscillations from the accretion disc of a single supermassive black hole---especially when the data span fewer than a few signal cycles. We present a Bayesian method for the detection of quasar (quasi-)periodicity in the presence of red noise. We apply this method to the binary candidate PG1302$-$102, and show that a) there is very strong support (Bayes factor $>10^6$) for quasi-periodicity, and b) the data slightly favour a quasi-periodic oscillation over a sinusoidal signal, which we interpret as modest evidence against the binary black hole hypothesis. We also find that the prevalent damped random walk red-noise model is disfavored with more than 99.9\% credibility. Finally, we outline future work that may enable the unambiguous identification of supermassive binary black holes.

preprint2019arXiv

Searching for Eccentricity: Signatures of Dynamical Formation in the First Gravitational-Wave Transient Catalogue of LIGO and Virgo

Binary black holes are thought to form primarily via two channels: isolated evolution and dynamical formation. The component masses, spins, and eccentricity of a binary black hole system provide clues to its formation history. We focus on eccentricity, which can be a signature of dynamical formation. Employing the spin-aligned eccentric waveform model SEOBNRE, we perform Bayesian inference to measure the eccentricity of binary black hole merger events in the first Gravitational-Wave Transient Catalogue of LIGO and Virgo. We find that all of these events are consistent with zero eccentricity. We set upper limits on eccentricity ranging from 0.02 to 0.05 with 90% confidence at a reference frequency of 10 Hz. These upper limits do not significantly constrain the fraction of LIGO-Virgo events formed dynamically in globular clusters, because only ~ 5% are expected to merge with measurable eccentricity. However, with the Gravitational-Wave Transient Catalogue set to expand dramatically over the coming months, it may soon be possible to significantly constrain the fraction of mergers taking place in globular clusters using eccentricity measurements.

preprint2019arXiv

Thanks for the memory: measuring gravitational-wave memory in the first LIGO/Virgo gravitational-wave transient catalog

Gravitational-wave memory, a strong-field effect of general relativity, manifests itself as a permanent displacement in spacetime. We develop a Bayesian framework to detect gravitational-wave memory with the Advanced LIGO/Virgo detector network. We apply this algorithm on the ten binary black hole mergers in LIGO/Virgo&#39;s first transient gravitational-wave catalog. We find no evidence of memory, which is consistent with expectations. In order to estimate when memory will be detected, we use the best current population estimates to construct a realistic sample of binary black hole observations for LIGO/Virgo at design sensitivity. We show that an ensemble of $\mathcal{O}(2000)$ binary black hole observations can be used to find definitive evidence for gravitational-wave memory. We conclude that memory is likely to be detected in the early A+/Virgo+ era.

preprint2019arXiv

Ultra-relativistic astrophysics using multi-messenger observations of double neutron stars with LISA and the SKA

Recent work highlights that tens of Galactic double neutron stars are likely to be detectable in the millihertz band of the space-based gravitational-wave observatory, LISA. Kyutoku and Nishino point out that some of these binaries might be detectable as radio pulsars using the Square Kilometer Array (SKA). We point out that the joint LISA+SKA detection of a $f_\text{gw}\gtrsim$1 mHz binary, corresponding to a binary period of $\lesssim$400 s, would enable precision measurements of ultra-relativistic phenomena. We show that, given plausible assumptions, multi-messenger observations of ultra-relativistic binaries can be used to constrain the neutron star equation of state with remarkable fidelity. It may be possible to measure the mass-radius relation with a precision of $\approx$0.2% after 10 yr of observations with the SKA. Such a measurement would be roughly an order of magnitude more precise than possible with other proposed observations. We summarize some of the other remarkable science made possible with multi-messenger observations of millihertz binaries, and discuss the prospects for the detection of such objects.