skip to main content


Search for: All records

Creators/Authors contains: "Ashton, G."

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. ABSTRACT We introduce Bilby-MCMC, a Markov chain Monte Carlo sampling algorithm tuned for the analysis of gravitational waves from merging compact objects. Bilby-MCMC provides a parallel-tempered ensemble Metropolis-Hastings sampler with access to a block-updating proposal library including problem-specific and machine learning proposals. We demonstrate that learning proposals can produce over a 10-fold improvement in efficiency by reducing the autocorrelation time. Using a variety of standard and problem-specific tests, we validate the ability of the Bilby-MCMC sampler to produce independent posterior samples and estimate the Bayesian evidence. Compared to the widely used Dynesty nested sampling algorithm, Bilby-MCMC is less efficient in producing independent posterior samples and less accurate in its estimation of the evidence. However, we find that posterior samples drawn from the Bilby-MCMC sampler are more robust: never failing to pass our validation tests. Meanwhile, the Dynesty sampler fails the difficult-to-sample Rosenbrock likelihood test, over constraining the posterior. For CBC problems, this highlights the importance of cross-sampler comparisons to ensure results are robust to sampling error. Finally, Bilby-MCMC can be embarrassingly and asynchronously parallelized making it highly suitable for reducing the analysis wall-time using a High Throughput Computing environment. Bilby-MCMC may be a useful tool for the rapid and robust analysis of gravitational-wave signals during the advanced detector era and we expect it to have utility throughout astrophysics. 
    more » « less
  2. ABSTRACT

    Radio pulsar glitches probe far-from-equilibrium processes involving stress accumulation and relaxation in neutron star interiors. Previous studies of glitch rates have focused on individual pulsars with as many recorded glitches as possible. In this work, we analyse glitch rates using all available data including objects that have glitched never or once. We assume the glitch rate follows a homogeneous Poisson process, and therefore exclude pulsars that exhibit quasiperiodic glitching behaviour. Calculating relevant Bayes factors shows that a model in which the glitch rate λ scales as a power of the characteristic age τ is preferred over models that depend arbitrarily on powers of the spin frequency ν and/or its time derivative $\dot{\nu }$. For λ = A(τ/τref)−γ, where τref = 1 yr is a reference time, the posterior distributions are unimodal with $A=0.0066_{-0.002}^{+0.003}\ \rm {yr}^{-1}$ and $\gamma =0.27_{-0.03}^{+0.03}$. Importantly, the data exclude with 99 per cent confidence the possibility γ = 1 canvassed in the literature. When objects with zero-recorded glitches are included, the age-based rate law is still preferred and the posteriors change to give $A=0.0099_{-0.003}^{+0.004}\ \rm {yr}^{-1}$ and $\gamma =0.31_{-0.03}^{+0.03}$. The updated estimates still support increased glitch activity for younger pulsars, while demonstrating that the large number of objects with zero glitches contain important statistical information about the rate, provided that they are part of the same population as opposed to a disjoint population which never glitches for some unknown physical reason.

     
    more » « less
  3. null (Ed.)
  4. Gravitational waves provide a unique tool for observational astronomy. While the first LIGO–Virgo catalogue of gravitational wave transients (GWTC-1) contains 11 signals from black hole and neutron star binaries, the number of observations is increasing rapidly as detector sensitivity improves. To extract information from the observed signals, it is imperative to have fast, flexible, and scalable inference techniques. In a previous paper, we introduced BILBY: a modular and user-friendly Bayesian inference library adapted to address the needs of gravitational-wave inference. In this work, we demonstrate that BILBY produces reliable results for simulated gravitational-wave signals from compact binary mergers, and verify that it accurately reproduces results reported for the 11 GWTC-1 signals. Additionally, we provide configuration and output files for all analyses to allow for easy reproduction, modification, and future use. This work establishes that BILBY is primed and ready to analyse the rapidly growing population of compact binary coalescence gravitational-wave signals. 
    more » « less
  5. ABSTRACT During ϕX174 morphogenesis, 240 copies of the external scaffolding protein D organize 12 pentameric assembly intermediates into procapsids, a reaction reconstituted in vitro . In previous studies, ϕX174 strains resistant to exogenously expressed dominant lethal D genes were experimentally evolved. Resistance was achieved by the stepwise acquisition of coat protein mutations. Once resistance was established, a stimulatory D protein mutation that greatly increased strain fitness arose. In this study, in vitro biophysical and biochemical methods were utilized to elucidate the mechanistic details and evolutionary trade-offs created by the resistance mutations. The kinetics of procapsid formation was analyzed in vitro using wild-type, inhibitory, and experimentally evolved coat and scaffolding proteins. Our data suggest that viral fitness is correlated with in vitro assembly kinetics and demonstrate that in vivo experimental evolution can be analyzed within an in vitro biophysical context. IMPORTANCE Experimental evolution is an extremely valuable tool. Comparisons between ancestral and evolved genotypes suggest hypotheses regarding adaptive mechanisms. However, it is not always possible to rigorously test these hypotheses in vivo . We applied in vitro biophysical and biochemical methods to elucidate the mechanistic details that allowed an experimentally evolved virus to become resistant to an antiviral protein and then evolve a productive use for that protein. Moreover, our results indicate that the respective roles of scaffolding and coat proteins may have been redistributed during the evolution of a two-scaffolding-protein system. In one-scaffolding-protein virus assembly systems, coat proteins promiscuously interact to form heterogeneous aberrant structures in the absence of scaffolding proteins. Thus, the scaffolding protein controls fidelity. During ϕX174 assembly, the external scaffolding protein acts like a coat protein, self-associating into large aberrant spherical structures in the absence of coat protein, whereas the coat protein appears to control fidelity. 
    more » « less
  6. Free, publicly-accessible full text available January 1, 2025
  7. Free, publicly-accessible full text available December 1, 2024
  8. Abstract

    We search for gravitational-wave (GW) transients associated with fast radio bursts (FRBs) detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project, during the first part of the third observing run of Advanced LIGO and Advanced Virgo (2019 April 1 15:00 UTC–2019 October 1 15:00 UTC). Triggers from 22 FRBs were analyzed with a search that targets both binary neutron star (BNS) and neutron star–black hole (NSBH) mergers. A targeted search for generic GW transients was conducted on 40 FRBs. We find no significant evidence for a GW association in either search. Given the large uncertainties in the distances of our FRB sample, we are unable to exclude the possibility of a GW association. Assessing the volumetric event rates of both FRB and binary mergers, an association is limited to 15% of the FRB population for BNS mergers or 1% for NSBH mergers. We report 90% confidence lower bounds on the distance to each FRB for a range of GW progenitor models and set upper limits on the energy emitted through GWs for a range of emission scenarios. We find values of order 1051–1057erg for models with central GW frequencies in the range 70–3560 Hz. At the sensitivity of this search, we find these limits to be above the predicted GW emissions for the models considered. We also find no significant coincident detection of GWs with the repeater, FRB 20200120E, which is the closest known extragalactic FRB.

     
    more » « less
    Free, publicly-accessible full text available September 28, 2024
  9. Abstract The global network of gravitational-wave observatories now includes five detectors, namely LIGO Hanford, LIGO Livingston, Virgo, KAGRA, and GEO 600. These detectors collected data during their third observing run, O3, composed of three phases: O3a starting in 2019 April and lasting six months, O3b starting in 2019 November and lasting five months, and O3GK starting in 2020 April and lasting two weeks. In this paper we describe these data and various other science products that can be freely accessed through the Gravitational Wave Open Science Center at https://gwosc.org . The main data set, consisting of the gravitational-wave strain time series that contains the astrophysical signals, is released together with supporting data useful for their analysis and documentation, tutorials, as well as analysis software packages. 
    more » « less
    Free, publicly-accessible full text available July 28, 2024