skip to main content


Title: Implementation of a generalized precession parameter in the RIFT parameter estimation algorithm
Abstract Since the initial discovery of gravitational waves in 2015, significant developments have been made towards waveform interpretation and estimation of compact binary source parameters. We present herein an implementation of the generalized precession parameter ⟨ χ p ⟩ [Gerosa et al 2021], which averages over all angular variations on the precession timescale, within the RIFT parameter estimation framework. Relative to the precession parameter χ p , which characterizes the single largest dynamical spin in a binary, ⟨ χ p ⟩ has a unique domain 1 < ⟨ χ p ⟩ < 2, which is exclusive to binaries with two precessing spins. After reviewing the physical differences between these two parameters, we describe how ⟨ χ p ⟩ was implemented in RIFT and apply it to all 36 events from the second half of the Advanced LIGO and Advanced Virgo third observing run (O3b). In O3b, ten events show significant amounts of precession ⟨ χ p ⟩ > 0.5. Of particular interest is GW191109_010717; we show it has a ∼ 28 % probability that the originating system necessarily contains two misaligned spins.  more » « less
Award ID(s):
2110481 1809572 2012057
PAR ID:
10343570
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Classical and Quantum Gravity
Volume:
39
Issue:
12
ISSN:
0264-9381
Page Range / eLocation ID:
125003
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Gravitational-wave observations provide the unique opportunity of studying black hole formation channels and histories—but only if we can identify their origin. One such formation mechanism is the dynamical synthesis of black hole binaries in dense stellar systems. Given the expected isotropic distribution of component spins of binary black holes in gas-free dynamical environments, the presence of antialigned or in-plane spins with respect to the orbital angular momentum is considered a tell-tale sign of a merger’s dynamical origin. Even in the scenario where birth spins of black holes are low, hierarchical mergers attain large component spins due to the orbital angular momentum of the prior merger. However, measuring such spin configurations is difficult. Here, we quantify the efficacy of the spin parameters encoding aligned-spin (χeff) and in-plane spin (χp) at classifying such hierarchical systems. Using Monte Carlo cluster simulations to generate a realistic distribution of hierarchical merger parameters from globular clusters, we can infer mergers’χeffandχp. The cluster populations are simulated using Advanced LIGO-Virgo sensitivity during the detector network’s third observing period and projections for design sensitivity. Using a “likelihood-ratio”-based statistic, we find that ∼2% of the recovered population by the current gravitational-wave detector network has a statistically significantχpmeasurement, whereas noχeffmeasurement was capable of confidently determining a system to be antialigned with the orbital angular momentum at current detector sensitivities. These results indicate that measuring spin-precession throughχpis a more detectable signature of hierarchical mergers and dynamical formation than antialigned spins.

     
    more » « less
  2. ABSTRACT

    We consider the 10 confidently detected gravitational-wave signals in the GWTC-1 catalog, which are consistent with mergers of binary black hole systems, and perform a thorough parameter estimation re-analysis. This is made possible by using computationally efficient waveform models of the current (fourth) generation of the IMRPhenom family of phenomenological waveform models, which consists of the IMRPhenomX frequency-domain models and the IMRPhenomT time-domain models.The analysis is performed with both precessing and non-precessing waveform models with and without subdominant spherical harmonic modes. Results for all events are validated with convergence tests, discussing in particular the events GW170729 and GW151226. For the latter and the other two lowest-mass events, we also compare results between two independent sampling codes, bilbyand lalinference. We find overall consistent results with the original GWTC-1 results, with all Jensen–Shannon divergences between the previous results using IMRPhenomPv2 and our default IMRPhenomXPHM posteriors below 0.045 bits. However, we also discuss cases where subdominant harmonics and/or precession influence the posteriors.

     
    more » « less
  3. Binary black holes with misaligned spins will generically induce both precession and nutation of the orbital angular momentum 𝐋 about the total angular momentum 𝐉. These phenomena modulate the phase and amplitude of the gravitational waves emitted as the binary inspirals to merger. We introduce a “taxonomy” of binary black-hole spin precession that encompasses all the known phenomenology, then present five new phenomenological parameters that describe generic precession and constitute potential building blocks for future gravitational waveform models. These are the precession amplitude ⟨𝜃𝐿⟩, the precession frequency ⟨Ω𝐿⟩, the nutation amplitude Δ⁢𝜃𝐿, the nutation frequency 𝜔, and the precession-frequency variation Δ⁢Ω𝐿. We investigate the evolution of these five parameters during the inspiral and explore their statistical properties for sources with isotropic spins. In particular, we find that nutation of 𝐋 is most prominent for binaries with high spins (𝜒≳0.5) and moderate mass ratios (𝑞∼0.6). 
    more » « less
  4. Emerging lithium-ion battery systems require high-fidelity electrochemical models for advanced control, diagnostics, and design. Accordingly, battery parameter estimation is an active research domain where novel algorithms are being developed to calibrate complex models from input-output data. Amidst these efforts, little focus has been placed on the fundamental mechanisms governing estimation accuracy, spurring the question, why is an estimate accurate or inaccurate? In response, we derive a generalized estimation error equation under the commonly adopted least-squares objective function, which reveals that the error can be represented as a combination of system uncertainties (i.e., in model, measurement, and parameter) and uncertainty-propagating sensitivity structures in the data. We then relate the error equation to conventional error analysis criteria, such as the Fisher information matrix, Cramér-Rao bound, and parameter sensitivity, to assess the benefits and limitations of each. The error equation is validated through several uni- and bivariate estimations of lithium-ion battery electrochemical parameters using experimental data. These results are also analyzed with the error equation to study the error compositions and parameter identifiability under different data. Finally, we show that adding target parameters to the estimation without increasing the amount of data intrinsically reduces the robustness of the results to system uncertainties.

     
    more » « less
  5. Nucleic acid kinetic simulators aim to predict the kinetics of interacting nucleic acid strands. Many simulators model the kinetics of interacting nucleic acid strands as continuous-time Markov chains (CTMCs). States of the CTMCs represent a collection of secondary structures, and transitions between the states correspond to the forming or breaking of base pairs and are determined by a nucleic acid kinetic model. The number of states these CTMCs can form may be exponentially large in the length of the strands, making two important tasks challenging, namely, mean first passage time (MFPT) estimation and parameter estimation for kinetic models based on MFPTs. Gillespie’s stochastic simulation algorithm (SSA) is widely used to analyze nucleic acid folding kinetics, but could be computationally expensive for reactions whose CTMC has a large state space or for slow reactions. It could also be expensive for arbitrary parameter sets that occur in parameter estimation. Our work addresses these two challenging tasks, in the full state space of all non-pseudoknotted secondary structures of each reaction. In the first task, we show how to use a reduced variance stochastic simulation algorithm (RVSSA), which is adapted from SSA, to estimate the MFPT of a reaction’s CTMC. In the second task, we estimate model parameters based on MFPTs. To this end, first, we show how to use a generalized method of moments (GMM) approach, where we minimize a squared norm of moment functions that we formulate based on experimental and estimated MFPTs. Second, to speed up parameter estimation, we introduce a fixed path ensemble inference (FPEI) approach, that we adapt from RVSSA. We implement and evaluate RVSSA and FPEI using the Multistrand kinetic simulator. In our experiments on a dataset of DNA reactions, FPEI speeds up parameter estimation compared to inference using SSA, by more than a factor of three for slow reactions. Also, for reactions with large state spaces, it speeds up parameter estimation by more than a factor of two. 
    more » « less