BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Paul D. Lasky

Paul D. Lasky appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
gr-qcastro-ph.HEastro-ph.COastro-ph.SRastro-ph.IMastro-phastro-ph.GAnucl-th

Catalog footprint

What is connected

40works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

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

0 citations0 reviewsNo code linkedOpen access
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}$.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

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

The "no-hair" 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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}$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Gravitational-wave astronomy with a physical calibration model

We carry out astrophysical inference for compact binary merger events in LIGO-Virgo'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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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'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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A neutron star progenitor for FRBs? Insights from polarisation measurements

Fast Radio Bursts (FRBs) are intense, millisecond-duration broadband radio transients, the emission mechanisms of which are not understood. Masui et al. recently presented Green Bank Telescope observations of FRB 110523, which displayed temporal variation of the linear polarisation position angle (PA). This effect is commonly seen in radio pulsars and is attributed to a changing projected magnetic field orientation in the emission region as the star rotates. If a neutron star is the progenitor of this FRB, and the emission mechanism is pulsar-like, we show that the progenitor is either rapidly rotating, or the emission originates from a region of complex magnetic field geometry. The observed PA variation could also be caused by propagation effects within a neutron-star magnetosphere, or by spatially varying magnetic fields if the progenitor lies in a dense, highly magnetised environment. Although we urge caution in generalising results from FRB 110523 to the broader FRB population, our analysis serves as a guide to interpreting future polarisation measurements of FRBs, and presents another means of elucidating the origins of these enigmatic ephemera.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detectability of Gravitational Waves from High-Redshift Binaries

Recent non-detection of gravitational-wave backgrounds from pulsar timing arrays casts further uncertainty on the evolution of supermassive black hole binaries. We study the capabilities of current gravitational-wave observatories to detect individual binaries and demonstrate that, contrary to conventional wisdom, some are in principle detectable throughout the Universe. In particular, a binary with rest-frame mass $\gtrsim10^{10}\,M_\odot$ can be detected by current timing arrays at arbitrarily high redshifts. The same claim will apply for less massive binaries with more sensitive future arrays. As a consequence, future searches for nanohertz gravitational waves could be expanded to target evolving high-redshift binaries. We calculate the maximum distance at which binaries can be observed with pulsar timing arrays and other detectors, properly accounting for redshift and using realistic binary waveforms.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Detecting gravitational-wave memory with LIGO: implications of GW150914

It may soon be possible for Advanced LIGO to detect hundreds of binary black hole mergers per year. We show how the accumulation of many such measurements will allow for the detection of gravitational-wave memory: a permanent displacement of spacetime that comes from strong-field, general relativistic effects. We estimate that Advanced LIGO operating at design sensitivity may be able to make a signal-to-noise ratio 3(5) detection of memory with ~35 (90) events with masses and distance similar to GW150914. We highlight the importance of incorporating higher-order gravitational-wave modes for parameter estimation of binary black hole mergers, and describe how our methods can also be used to detect higher-order modes themselves before Advanced LIGO reaches design sensitivity.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Gravitational-wave cosmology across 29 decades in frequency

Quantum fluctuations of the gravitational field in the early Universe, amplified by inflation, produce a primordial gravitational-wave background across a broad frequency band. We derive constraints on the spectrum of this gravitational radiation, and hence on theories of the early Universe, by combining experiments that cover 29 orders of magnitude in frequency. These include Planck observations of cosmic microwave background temperature and polarization power spectra and lensing, together with baryon acoustic oscillations and big bang nucleosynthesis measurements, as well as new pulsar timing array and ground-based interferometer limits. While individual experiments constrain the gravitational-wave energy density in specific frequency bands, the combination of experiments allows us to constrain cosmological parameters, including the inflationary spectral index, $n_t$, and the tensor-to-scalar ratio, $r$. Results from individual experiments include the most stringent nanohertz limit of the primordial background to date from the Parkes Pulsar Timing Array, $Ω_{\rm gw}(f)<2.3\times10^{-10}$. Observations of the cosmic microwave background alone limit the gravitational-wave spectral index at 95\% confidence to $n_t\lesssim5$ for a tensor-to-scalar ratio of $r = 0.11$. However, the combination of all the above experiments limits $n_t<0.36$. Future Advanced LIGO observations are expected to further constrain $n_t<0.34$ by 2020. When cosmic microwave background experiments detect a non-zero $r$, our results will imply even more stringent constraints on $n_t$ and hence theories of the early Universe.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The freedom to choose neutron star magnetic field equilibria

Our ability to interpret and glean useful information from the large body of observations of strongly magnetised neutron stars rests largely on our theoretical understanding of magnetic field equilibria. We answer the following question: is one free to arbitrarily prescribe magnetic equilibria such that fluid degrees of freedom can balance the equilibrium equations? We examine this question for various models for neutron star matter; from the simplest single-fluid barotrope to more realistic non-barotropic multifluid models with superfluid/superconducting components, muons and entropy. We do this for both axi- and non-axisymmetric equilibria, and in Newtonian gravity and general relativity. We show that, in axisymmetry, the most realistic model allows complete freedom in choosing a magnetic field equilibrium whereas non-axisymmetric equilibria are never completely arbitrary.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gravitational Waves from Neutron Stars: A Review

Neutron stars are excellent emitters of gravitational waves. Squeezing matter beyond nuclear densities invites exotic physical processes, many of which violently transfer large amounts of mass at relativistic velocities, disrupting spacetime and generating copious quantities of gravitational radiation. I review mechanisms for generating gravitational waves with neutron stars. This includes gravitational waves from radio and millisecond pulsars, magnetars, accreting systems and newly born neutron stars, with mechanisms including magnetic and thermoelastic deformations, various stellar oscillation modes and core superfluid turbulence. I also focus on what physics can be learnt from a gravitational wave detection, and where additional research is required to fully understand the dominant physical processes at play.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pulsar timing noise and the minimum observation time to detect gravitational waves with pulsar timing arrays

The sensitivity of pulsar timing arrays to gravitational waves is, at some level, limited by timing noise. Red timing noise - the stochastic wandering of pulse arrival times with a red spectrum - is prevalent in slow-spinning pulsars and has been identified in many millisecond pulsars. Phenomenological models of timing noise, such as from superfluid turbulence, suggest that the timing noise spectrum plateaus below some critical frequency, $f_c$, potentially aiding the hunt for gravitational waves. We examine this effect for individual pulsars by calculating minimum observation times, $T_{\rm min}(f_c)$, over which the gravitational wave signal becomes larger than the timing noise plateau. We do this in two ways: 1) in a model-independent manner, and 2) by using the superfluid turbulence model for timing noise as an example to illustrate how neutron star parameters can be constrained. We show that the superfluid turbulence model can reproduce the data qualitatively from a number of pulsars observed as part of the Parkes Pulsar Timing Array. We further show how a value of $f_c$, derived either through observations or theory, can be related to $T_{\rm min}$. This provides a diagnostic whereby the usefulness of timing array pulsars for gravitational-wave detection can be quantified.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gravitational waves from rapidly rotating neutron stars

Rapidly rotating neutron stars in Low Mass X-ray Binaries have been proposed as an interesting source of gravitational waves. In this chapter we present estimates of the gravitational wave emission for various scenarios, given the (electromagnetically) observed characteristics of these systems. First of all we focus on the r-mode instability and show that a 'minimal' neutron star model (which does not incorporate exotica in the core, dynamically important magnetic fields or superfluid degrees of freedom), is not consistent with observations. We then present estimates of both thermally induced and magnetically sustained mountains in the crust. In general magnetic mountains are likely to be detectable only if the buried magnetic field of the star is of the order of $B\approx 10^{12}$ G. In the thermal mountain case we find that gravitational wave emission from persistent systems may be detected by ground based interferometers. Finally we re-asses the idea that gravitational wave emission may be balancing the accretion torque in these systems, and show that in most cases the disc/magnetosphere interaction can account for the observed spin periods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nuclear Equation of State from Observations of Short Gamma-Ray Burst Remnants

The favoured progenitor model for short $γ$-ray bursts (SGRBs) is the merger of two neutron stars that triggers an explosion with a burst of collimated $γ$-rays. Following the initial prompt emission, some SGRBs exhibit a plateau phase in their $X$-ray light curves that indicates additional energy injection from a central engine, believed to be a rapidly rotating, highly magnetised neutron star. The collapse of this `protomagnetar' to a black hole is likely to be responsible for a steep decay in $X$-ray flux observed at the end of the plateau. In this letter, we show that these observations can be used to effectively constrain the equation of state of dense matter. In particular, we show that the known distribution of masses in binary neutron star systems, together with fits to the $X$-ray light curves, provide constraints that exclude the softest and stiffest plausible equations of state. We further illustrate how a future gravitational wave observation with Advanced LIGO/Virgo can place tight constraints on the equation of state, by adding into the picture a measurement of the chirp mass of the SGRB progenitor.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The birth of black holes: neutron star collapse times, gamma-ray bursts and fast radio bursts

Recent observations of short gamma-ray bursts (SGRBs) suggest that binary neutron star (NS) mergers can create highly magnetised, millisecond NSs. Sharp cut-offs in X-ray afterglow plateaus of some SGRBs hint at the gravitational collapse of these remnant NSs to black holes. The collapse of such `supramassive' NSs also describes the blitzar model, a leading candidate for the progenitors of fast radio bursts (FRBs). The observation of an FRB associated with an SGRB would provide compelling evidence for the blitzar model and the binary NS merger scenario of SGRBs, and lead to interesting constraints on the NS equation of state. We predict the collapse times of supramassive NSs created in binary NS mergers, finding that such stars collapse $\sim10\,{\rm s}$ -- $4.4\times10^{4}\,{\rm s}$ (95% confidence) after the merger. This directly impacts observations targeting NS remnants of binary NS mergers, providing the optimal window for high time resolution radio and X-ray follow-up of SGRBs and gravitational wave bursts.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Neutron star deformation due to multipolar magnetic fields

Certain multi-wavelength observations of neutron stars, such as intermittent radio emissions from rotation-powered pulsars beyond the pair-cascade death line, the pulse profile of the magnetar SGR 1900+14 after its 1998 August 27 giant flare, and X-ray spectral features of PSR J0821-4300 and SGR 0418+5729, suggest that the magnetic fields of non-accreting neutron stars are not purely dipolar and may contain higher-order multipoles. Here, we calculate the ellipticity of a non-barotropic neutron star with (i) a quadrupole poloidal-toroidal field, and (ii) a purely poloidal field containing arbitrary multipoles, deriving the relation between the ellipticity and the multipole amplitudes. We present, as a worked example, a purely poloidal field comprising dipole, quadrupole, and octupole components. We show the correlation between field energy and ellipticity for each multipole, that the l=4 multipole has the lowest energy, and that l=5 has the lowest ellipticity. We show how a mixed multipolar field creates an observationally testable mismatch between the principal axes of inertia (to be inferred from gravitational wave data) and the magnetic inclination angle. Strong quadrupole and octupole components (with amplitudes ~ 10^2 times higher than the dipole) in SGR 0418+5729 still yield ellipticity ~ 10^{-8}, consistent with current gravitational wave upper limits. The existence of higher multipoles in fast-rotating objects (e.g., newborn magnetars) has interesting implications for the braking law and hence phase tracking during coherent gravitational wave searches.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stochastic gravitational wave background from hydrodynamic turbulence in differentially rotating neutron stars

Hydrodynamic turbulence driven by crust-core differential rotation imposes a fundamental noise floor on gravitational wave observations of neutron stars. The gravitational wave emission peaks at the Kolmogorov decoherence frequency which, for reasonable values of the crust-core shear, ΔΩ, occurs near the most sensitive part of the frequency band for ground-based, long-baseline interferometers. We calculate the energy density spectrum of the stochastic gravitational wave background from a cosmological population of turbulent neutron stars generalising previous calculations for individual sources. The spectrum resembles a piecewise power law, Ω_{gw}(ν)=Ω_αν^α, with α=-1 and 7 above and below the decoherence frequency respectively, and its normalisation scales as Ω_α\propto(ΔΩ)^{7}. Non-detection of a stochastic signal by Initial LIGO implies an upper limit on ΔΩand hence by implication on the internal relaxation time-scale for the crust and core to come into co-rotation, τ_{d}=ΔΩ/\dotΩ, where \dotΩ is the observed electromagnetic spin-down rate, with τ_{d}\lesssim 10^{7} yr for accreting millisecond pulsars and τ_{d}\lesssim 10^{5} yr for radio-loud pulsars. Target limits on τ_{d} are also estimated for future detectors, namely Advanced LIGO and the Einstein Telescope, and are found to be astrophysically interesting.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Tilted torus magnetic fields in neutron stars and their gravitational wave signatures

Gravitational-wave diagnostics are developed for discriminating between varieties of mixed poloidal-toroidal magnetic fields in neutron stars, with particular emphasis on differentially rotating protoneutron stars. It is shown that tilted torus magnetic fields, defined as the sum of an internal/external poloidal component, whose axis of symmetry is tilted with respect to the rotation axis, and an internal toroidal component, whose axis of symmetry is aligned with the rotation axis, deform the star triaxially, unlike twisted torus fields, which deform the star biaxially. Utilizing an analytic tilted torus example, we show that these two topologies can be distinguished by their gravitational wave spectrum and polarization phase portraits. For example, the relative amplitudes and frequencies of the spectral peaks allows one to infer the relative strengths of the toroidal and poloidal components of the field, and the magnetic inclination angle. Finally, we show how a tilted torus field arises naturally from magnetohydrodynamic simulations of differentially rotating neutron stars, and how the gravitational wave spectrum evolves as the internal toroidal field winds up. These results point to the sorts of experiments that may become possible once gravitational wave interferometers detect core-collapse supernovae routinely.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravitational Waves and Hydromagnetic Instabilities in Rotating Magnetized Neutron Stars

We perform nonlinear general relativistic ideal magnetohydrodynamic simulations of poloidal magnetic fields in rotating polytropic neutron stars. We have three primary goals: i) to understand the nature of magnetohydrodynamic instabilities inherent to poloidal magnetic fields in non-rotating and rotating neutron stars, ii) to explore the possible space of stable equilibrium configurations and iii) to understand gravitational wave emissions caused by the catastrophic reconfiguration of magnetic fields associated with giant magnetar flares. Our key physical contributions can be summarized as follows: i) gravitational waves from f-modes caused by magnetar flares are unlikely to be detected in the current or near-future generation of gravitational waves observatories, ii) gravitational waves from Alfven waves propagating inside the neutron star are more likely candidates, although this interpretation relies on the unknown damping time of these modes, iii) any magnetic field equilibria derived from our simulations are characterized as non-axisymmetric, with approximately 65% of their magnetic energy in the poloidal field, iv) rotation acts to separate the timescales of different instabilities in our system, with the varicose mode playing a more major role due to a delayed kink instability and v) despite the slowing growth rate of the kink mode, it is always present in our simulations, even for models where the rotational period is of the same order as the Alfven timescale.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Are gravitational waves from giant magnetar flares observable?

Are giant flares in magnetars viable sources of gravitational radiation? Few theoretical studies have been concerned with this problem, with the small number using either highly idealized models or assuming a magnetic field orders of magnitude beyond what is supported by observations. We perform nonlinear general-relativistic magnetohydrodynamics simulations of large-scale hydromagnetic instabilities in magnetar models. We utilise these models to find gravitational wave emissions over a wide range of energies, from 10^40 to 10^47 erg. This allows us to derive a systematic relationship between the surface field strength and the gravitational wave strain, which we find to be highly nonlinear. In particular, for typical magnetar fields of a few times 10^15 G, we conclude that a direct observation of f-modes excited by global magnetic field reconfigurations is unlikely with present or near-future gravitational wave observatories, though we also discuss the possibility that modes in a low-frequency band up to 100 Hz could be sufficiently excited to be relevant for observation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravitational Lensing with Three-Dimensional Ray Tracing

High redshift sources suffer from magnification or demagnification due to weak gravitational lensing by large scale structure. One consequence of this is that the distance-redshift relation, in wide use for cosmological tests, suffers lensing-induced scatter which can be quantified by the magnification probability distribution. Predicting this distribution generally requires a method for ray-tracing through cosmological N-body simulations. However, standard methods tend to apply the multiple thin-lens approximation. In an effort to quantify the accuracy of these methods, we develop an innovative code that performs ray-tracing without the use of this approximation. The efficiency and accuracy of this computationally challenging approach can be improved by careful choices of numerical parameters; therefore, the results are analysed for the behaviour of the ray-tracing code in the vicinity of Schwarzschild and Navarro-Frenk-White lenses. Preliminary comparisons are drawn with the multiple lens-plane ray-bundle method in the context of cosmological mass distributions for a source redshift of $z_{s}=0.5$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Hydromagnetic Instabilities in Neutron Stars

We model the non-linear ideal magnetohydrodynamics of poloidal magnetic fields in neutron stars in general relativity assuming a polytropic equation of state. We identify familiar hydromagnetic modes, in particular the 'sausage/varicose' mode and 'kink' instability inherent to poloidal magnetic fields. The evolution is dominated by the kink instability, which causes a cataclysmic reconfiguration of the magnetic field. The system subsequently evolves to new, non-axisymmetric, quasi-equilibrium end-states. The existence of this branch of stable quasi-equilibria may have consequences for magnetar physics, including flare generation mechanisms and interpretations of quasi-periodic oscillations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Shape, shear and flexion II - Quantifying the flexion formalism for extended sources with the ray-bundle method

Flexion-based weak gravitational lensing analysis is proving to be a useful adjunct to traditional shear-based techniques. As flexion arises from gradients across an image, analytic and numerical techniques are required to investigate flexion predictions for extended image/source pairs. Using the Schwarzschild lens model, we demonstrate that the ray-bundle method for gravitational lensing can be used to accurately recover second flexion, and is consistent with recovery of zero first flexion. Using lens plane to source plane bundle propagation, we find that second flexion can be recovered with an error no worse than 1% for bundle radii smaller than Δθ = 0.01 θ_E and lens plane impact pararameters greater than θ_E + Δθ, where θ_E is the angular Einstein radius. Using source plane to lens plane bundle propagation, we demonstrate the existence of a preferred flexion zone. For images at radii closer to the lens than the inner boundary of this zone, indicative of the true strong lensing regime, the flexion formalism should be used with caution (errors greater than 5% for extended image/source pairs). We also define a shear zone boundary, beyond which image shapes are essentially indistinguishable from ellipses (1% error in ellipticity). While suggestive that a traditional weak lensing analysis is satisfactory beyond this boundary, a potentially detectable non-zero flexion signal remains.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stability and Quasinormal Modes of Black holes in Tensor-Vector-Scalar theory: Scalar Field Perturbations

The imminent detection of gravitational waves will trigger precision tests of gravity through observations of quasinormal ringing of black holes. While General Relativity predicts just two polarizations of gravitational waves, the so-called plus and cross polarizations, numerous alternative theories of gravity predict up to six different polarizations which will potentially be observed in current and future generations of gravitational wave detectors. Bekenstein's Tensor-Vector-Scalar (TeVeS) theory and its generalization fall into one such class of theory that predict the full gamut of six polarizations of gravitational waves. In this paper we begin the study of quasinormal modes (QNMs) in TeVeS by studying perturbations of the scalar field in a spherically symmetric background. We show that, at least in the case where superluminal propagation of perturbations is not present, black holes are generically stable to this kind of perturbation. We also make a unique prediction that, as the limit of the various coupling parameters of the theory tend to zero, the QNM spectrum tends to $1/\sqrt{2}$ times the QNM spectrum induced by scalar perturbations of a Schwarzschild black hole in General Relativity due to the intrinsic presence of the background vector field. We further show that the QNM spectrum does not vary significantly from this value for small values of the theory's coupling parameters, however can vary by as much as a few percent for larger, but still physically relevant parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The effect of pressure gradients on luminosity distance - redshift relations

Inhomogeneous cosmological models have had significant success in explaining cosmological observations without the need for dark energy. Generally, these models imply inhomogeneous matter distributions alter the observable relations that are taken for granted when assuming the Universe evolves according to the standard Friedmann equations. Moreover, it has recently been shown that both inhomogeneous matter and pressure distributions are required in both early and late stages of cosmological evolution. These associated pressure gradients are required in the early Universe to sufficiently describe void formation, whilst late-stage pressure gradients stop the appearance of anomalous singularities. In this paper we investigate the effect of pressure gradients on cosmological observations by deriving the luminosity distance - redshift relations in spherically symmetric, inhomogeneous spacetimes endowed with a perfect fluid. By applying this to a specific example for the energy density distribution and using various equations of state, we are able to explicitly show that pressure gradients may have a non-negligble effect on cosmological observations. In particular, we show that a non-zero pressure gradient can imply significantly different residual Hubble diagrams for $z\lesssim1$ compared to when the pressure is ignored. This paper therefore highlights the need to properly consider pressure gradients when interpreting cosmological observations.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Structure of Neutron Stars in Tensor-Vector-Scalar Theory

Bekenstein's Tensor-Vector-Scalar (TeVeS) theory has had considerable success in explaining various phenomena without the need for dark matter. However, it is difficult to observationally discern the differences between TeVeS and predictions made within the Lambda-cold dark matter concordance model. This implies that alternative tests are required that independently verify which theory is correct. For this we turn to the strong-field regime of TeVeS. In particular, we solve the spherically symmetric equations of hydrostatic equilibrium for a perfect fluid with a realistic equation of state to build models of neutron stars in TeVeS. We show that causality within the neutron star is only maintained for certain cosmological values of the scalar field, which allows us to put constraints on this value independently of cosmological observations. We also discuss in detail the internal structure of neutron stars and how each of the free parameters in the theory effects the overall size and mass of the neutron stars. In particular, the radii of neutron stars in TeVeS can significantly differ from those in General Relativity for certain values of the vector field coupling, which allows us to also place extra constraints on this parameter. Finally, we discuss future observations of neutron stars using both the electromagnetic and gravitational wave spectrums that will allow for tests of the appropriate theory of gravity.

0 citations0 reviewsNo code linkedOpen access