Researcher profile

Paul D. Lasky

Paul D. Lasky contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

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

19 published item(s)

preprint2023arXiv

Detectability of the gravitational-wave background produced by magnetar giant flares

We study the gravitational-wave background produced by f-mode oscillations of neutron stars triggered by magnetar giant flares. For the gravitational-wave energy, we use analytic formulae obtained via general relativistic magnetohydrodynamic simulations of strongly magnetized neutron stars. Assuming the magnetar giant flare rate is proportional to the star-formation rate, we show the gravitational-wave signal is likely undetectable by third-generation detectors such as the Einstein Telescope and Cosmic Explorer. We calculate the minimum value of the magnetic field and the magnetar giant flare rate necessary for such a signal to be detectable, and discuss these in the context of our current understanding of magnetar flares throughout the Universe.

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

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

Low-efficiency long gamma-ray bursts: A case study with AT2020blt

The Zwicky Transient Facility recently announced the detection of an optical transient AT2020blt at redshift $z=2.9$, consistent with the afterglow of an on-axis gamma-ray burst. However, no prompt emission was observed. We analyse AT2020blt with detailed models, showing the data are best explained as the afterglow of an on-axis long gamma-ray burst, ruling out other hypotheses such as a cocoon and a low-Lorentz factor jet. We search \textit{Fermi} data for prompt emission, setting deeper upper limits on the prompt emission than in the original detection paper. Together with KONUS-\textit{Wind} observations, we show that the gamma-ray efficiency of AT2020blt is $\lesssim 0.3 - 4.5\%$. We speculate that AT2020blt and AT2021any belong to the low-efficiency tail of long gamma-ray burst distributions that are beginning to be readily observed due to the capabilities of new observatories like the Zwicky Transient Facility.

preprint2022arXiv

Multimessenger astronomy with a kHz-band gravitational-wave observatory

Proposed next-generation networks of gravitational-wave observatories include dedicated kilohertz instruments that target neutron star science, such as the proposed Neutron Star Extreme Matter Observatory, NEMO. The original proposal for NEMO highlighted the need for it to exist in a network of gravitational-wave observatories to ensure detection confidence and sky localisation of sources. We show that NEMO-like observatories have significant utility on their own as coincident electromagnetic observations can provide the detection significance and sky localisation. We show that, with a single NEMO-like detector and expected electromagnetic observatories in the late 2020s and early 2030s such as the Vera C. Rubin observatory and SVOM, approximately $40\%$ of all binary neutron star mergers detected with gravitational waves could be confidently identified as coincident multimessenger detections. We show that we expect $2^{+10}_{-1}\rm{yr}^{-1}$ coincident observations of gravitational-wave mergers with gamma-ray burst prompt emission, $13^{+23}_{-10}\rm{yr}^{-1}$ detections with kilonova observations, and $4^{+18}_{-3}\rm{yr}^{-1}$ with broadband afterglows and kilonovae, where the uncertainties are $90\%$ confidence intervals arising from uncertainty in current merger-rate estimates. Combined, this implies a coincident detection rate of $14^{+25}_{-11}\rm{yr}^{-1}$ out to $300\rm{Mpc}$. These numbers indicate significant science potential for a single kilohertz gravitational-wave detector operating without a global network of other gravitational-wave observatories.

preprint2022arXiv

The evolution of binary neutron star post-merger remnants: a review

Two neutron stars merge somewhere in the Universe approximately every 10 seconds, creating violent explosions observable in gravitational waves and across the electromagnetic spectrum. The transformative coincident gravitational-wave and electromagnetic observations of the binary neutron star merger GW170817 gave invaluable insights into these cataclysmic collisions, probing bulk nuclear matter at supranuclear densities, the jet structure of gamma-ray bursts, the speed of gravity, and the cosmological evolution of the local Universe, among other things. Despite the wealth of information, it is still unclear when the remnant of GW170817 collapsed to form a black hole. Evidence from other short gamma-ray bursts indicates a large fraction of mergers may form long-lived neutron stars. We review what is known observationally and theoretically about binary neutron star post-merger remnants. From a theoretical perspective, we review our understanding of the evolution of short- and long-lived merger remnants, including fluid, magnetic-field, and temperature evolution. These considerations impact prospects of detection of gravitational waves from either short- or long-lived neutron star remnants which potentially allows for new probes into the hot nuclear equation of state in conditions inaccessible in terrestrial experiments. We also review prospects for determining post-merger physics from current and future electromagnetic observations, including kilonovae and late-time x-ray and radio afterglow observations.

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

Pitfalls of periodograms: The non-stationarity bias in the analysis of quasi-periodic oscillations

Quasi-periodic oscillations (QPOs) are an important key to understand the dynamic behavior of astrophysical objects during transient events like gamma-ray bursts, solar flares, and magnetar flares. Searches for QPOs often use the periodogram of the time series and perform spectral density estimation using a Whittle likelihood function. However, the Whittle likelihood is only valid if the time series is stationary since the frequency bins are otherwise not statistically independent. We show that if time series are non-stationary, the significance of QPOs can be highly overestimated and estimates of the central frequencies and QPO widths can be overconstrained. The effect occurs if the QPO is only present for a fraction of the time series and the noise level is varying throughout the time series. This can occur for example if background noise from before or after the transient is included in the time series or if the low-frequency noise profile varies strongly over the time series. We confirm the presence of this bias in previously reported results from solar flare data and show that significance can be highly overstated. Finally, we provide some suggestions that help identify if an analysis is affected by this bias.

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

Detection and parameter estimation of binary neutron star merger remnants

Detection and parameter estimation of binary neutron star merger remnants can shed light on the physics of hot matter at supranuclear densities. Here we develop a fast, simple model that can generate gravitational waveforms, and show it can be used for both detection and parameter estimation of post-merger remnants. The model consists of three exponentially-damped sinusoids with a linear frequency-drift term. The median fitting factors between the model waveforms and numerical-relativity simulations exceed 0.90. We detect remnants at a post-merger signal-to-noise ratio of $\ge 7$ using a Bayes-factor detection statistic with a threshold of 3000. We can constrain the primary post-merger frequency to $\pm_{1.2}^{1.4}\%$ at post-merger signal-to-noise ratios of 15 with an increase in precision to $\pm_{0.2}^{0.3}\%$ for post-merger signal-to-noise ratios of 50. The tidal coupling constant can be constrained to $\pm^{9}_{12}\%$ at post-merger signal-to-noise ratios of 15, and $\pm 5\%$ at post-merger signal-to-noise ratios of 50 using a hierarchical inference model.

preprint2020arXiv

Gravitational waves or deconfined quarks: what causes the premature collapse of neutron stars born in short gamma-ray bursts?

We infer the collapse times of long-lived neutron stars into black holes using the X-ray afterglows of 18 short gamma-ray bursts. We then apply hierarchical inference to infer properties of the neutron star equation of state and dominant spin-down mechanism. We measure the maximum non-rotating neutron star mass $M_\mathrm{TOV} = 2.31 ^{+0.36}_{-0.21} M_{\odot}$ and constrain the fraction of remnants spinning down predominantly through gravitational-wave emission to $η= 0.69 ^{+0.21}_{-0.39}$ with $68 \%$ uncertainties. In principle, this method can determine the difference between hadronic and quark equation of states. In practice, however, the data is not yet informative with indications that these neutron stars do not have hadronic equation of states at the $1σ$ level. These inferences all depend on the underlying progenitor mass distribution for short gamma-ray bursts produced by binary neutron star mergers. The recently announced gravitational-wave detection of GW190425 suggests this underlying distribution is different from the locally-measured population of double neutron stars. We show that $M_\mathrm{TOV}$ and $η$ constraints depend on the fraction of binary mergers that form through a distribution consistent with the locally-measured population and a distribution that can explain GW190425. The more binaries that form from the latter distribution, the larger $M_\mathrm{TOV}$ needs to be to satisfy the X-ray observations. Our measurements above are marginalised over this unknown fraction. If instead, we assume GW190425 is not a binary neutron star merger, i.e the underlying mass distribution of double neutron stars is the same as observed locally, we measure $M_\mathrm{TOV} = 2.26 ^{+0.31}_{-0.17} M_{\odot}$.

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

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.

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.