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 (
- PAR ID:
- 10343570
- 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
-
Abstract χ 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χ p measurement, 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χ p is a more detectable signature of hierarchical mergers and dynamical formation than antialigned spins. -
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.
-
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
-
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.
-
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