Title: Minimally-modeled search of higher multipole gravitational-wave radiation in compact binary coalescences
Abstract As the Advanced LIGO and Advanced Virgo interferometers, soon to be joined by the KAGRA interferometer, increase their sensitivity, they detect an ever-larger number of gravitational waves with a significant presence of higher multipoles (HMs) in addition to the dominant (2, 2) multipole. These HMs can be detected with different approaches, such as the minimally-modeled burst search methods, and here we discuss one such approach based on the coherent WaveBurst (cWB) pipeline. During the inspiral phase the HMs produce chirps whose instantaneous frequency is a multiple of the dominant (2, 2) multipole, and here we describe how cWB can be used to detect these spectral features. The search is performed within suitable regions of the time-frequency representation; their shape is determined by optimizing the receiver operating characteristics. This novel method has already been used in the GW190814 discovery paper (Abbott et al 2020 Astrophys. J. Lett. 896 L44) and is very fast and flexible. Here we describe in full detail the procedure used to detect the (3, 3) multipole in GW190814 as well as searches for other HMs during the inspiral phase, and apply it to another event that displays HMs, GW190412, replicating the results obtained with different methods. The procedure described here can be used for the fast analysis of HMs and to support the findings obtained with the model-based Bayesian parameter estimates.  more » « less
Award ID(s):
2110060 1806165
PAR ID:
10344428
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Classical and Quantum Gravity
Volume:
39
Issue:
4
ISSN:
0264-9381
Page Range / eLocation ID:
045001
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract The coherent WaveBurst (cWB) pipeline implements a minimally-modelled search to find a coherent response in the network of gravitational wave detectors of the LIGO-Virgo Col-laboration in the time-frequency domain. In this manuscript, we provide a timely introduction to an extension of the cWB analysis to detect spectral features beyond the main quadrupolar emission of gravitational waves during the inspiral phase of compact binary coalescences; more detailed discussion will be provided in a forthcoming paper [1]. The search is performed by defining specific regions in the time-frequency map to extract the energy of harmonics of main quadrupole mode in the inspiral phase. This method has already been used in the GW190814 discovery paper (Astrophys. J. Lett. 896 L44). Here we show the procedure to detect the (3, 3) multipole in GW190814 within the cWB framework. 
    more » « less
  2. Amplitude and phase of the gravitational waveform from compact binary systems can be decomposed in terms of their mass- and current-type multipole moments. In a modified theory of gravity, one or more of these multipole moments could deviate from general theory of relativity. In this work, we show that a waveform model that parametrizes the amplitude and phase in terms of the multipole moments of the binary can facilitate a novel multiparameter test of general relativity with exquisite precision. Using a network of next-generation gravitational-wave observatories, simultaneous deviation in the leading seven multipoles of a GW190814-like binary can be bounded to within 6%–40% depending on the multipole order, while supermassive black hole mergers observed by the Laser Interferometer Space Antenna achieve a bound of 0.3%–2%. We further argue that bounds from multipoles can be uniquely mapped onto other parametrized tests of general relativity and have the potential to become a downstream analysis from which bounds of other parametric tests of general relativity can be derived. The set of multipole parameters, therefore, provides an excellent basis to carry out precision tests of general relativity. 
    more » « less
  3. We study the possible occurrence of the hadron-quark phase transition (PT) during the merging of neutron star binaries by hydrodynamical simulations employing a set of temperature dependent hybrid equations of state (EoSs). Following previous work we describe an unambiguous and measurable signature of deconfined quark matter in the gravitational-wave (GW) signal of neutron star binary mergers including equal-mass and unequal-mass systems of different total binary mass. The softening of the EoS by the PT at higher densities, i.e. after merging, leads to a characteristic increase of the dominant postmerger GW frequency f_peak relative to the tidal deformability Lambda inferred during the premerger inspiral phase. Hence, measuring such an increase of the postmerger frequency provides evidence for the presence of a strong PT. If the postmerger frequency and the tidal deformability are compatible with results from purely baryonic EoS models yielding very tight relations between f_peak and Lambda, a strong PT can be excluded up to a certain density. We find tight correlations of f_peak and Lambda with the maximum density during the early postmerger remnant evolution. These GW observables thus inform about the density regime which is probed by the remnant and its GW emission. Exploiting such relations we devise a directly applicable, concrete procedure to constrain the onset density of the QCD PT from future GW measurements. We point out two interesting scenarios: if no indications for a PT are inferred from a GW detection, our procedure yields a lower limit on the onset density of the hadron quark PT. On the contrary, if a merger event reveals evidence for the occurrence of deconfined quark matter, the inferred GW parameters set an upper limit on the PT onset density. (abridged) 
    more » « less
  4. The formulation of active pharmaceutical ingredients involves discovering stable crystal packing arrangements or polymorphs, each of which has distinct pharmaceutically relevant properties. Traditional experimental screening techniques utilizing various conditions are commonly supplemented with in silico crystal structure prediction (CSP) to inform the crystallization process and mitigate risk. Predictions are often based on advanced classical force fields or quantum mechanical calculations that model the crystal potential energy landscape but do not fully incorporate temperature, pressure, or solution conditions during the search procedure. This study proposes an innovative alchemical path that utilizes an advanced polarizable atomic multipole force field to predict crystal structures based on direct sampling of the NPT ensemble. The use of alchemical (i.e., nonphysical) intermediates, a novel Monte Carlo barostat, and an orthogonal space tempering bias combine to enhance the sampling efficiency of the deposition/sublimation phase transition. The proposed algorithm was applied to 2-((4-(2-(3,4-dichlorophenyl)ethyl)phenyl)amino)benzoic acid (Cambridge Crystallography Database Centre ID: XAFPAY) as a case study to showcase the algorithm. Each experimentally determined polymorph with one molecule in the asymmetric unit was successfully reproduced via approximately 1000 short 1 ns simulations per space group where each simulation was initiated from random rigid body coordinates and unit cell parameters. Utilizing two threads of a recent Intel CPU (a Xeon Gold 6330 CPU at 2.00 GHz), 1 ns of sampling using the polarizable AMOEBA force field can be acquired in 4 h (equating to more than 300 ns/day using all 112 threads/56 cores of a dual CPU node) within the Force Field X software (https://ffx.biochem.uiowa.edu). These results demonstrate a step forward in the rigorous use of the NPT ensemble during the CSP search process and open the door to future algorithms that incorporate solution conditions using continuum solvation methods. 
    more » « less
  5. Subjected to the tidal field of its companion, each component of a coalescing binary suffers a slow change in its mass (tidal heating) and spin (tidal torquing) during the inspiral and merger. For black holes, these changes are associated with their absorption of energy and angular momentum fluxes. This effect modifies the inspiral rate of the binary, and consequently, the phase and amplitude of its gravitational waveform. Numerical relativity (NR) waveforms contain these effects inherently, whereas analytical approximants for the early inspiral phase have to include them manually in the energy balance equation. In this work, we construct IMRPhenomD_Horizon, a frequency-domain gravitational waveform model that incorporates the effects of tidal heating of black holes. This is achieved by recalibrating the inspiral phase of the waveform model IMRPhenomD to incorporate the phase corrections for tidal heating. We also include corrections to the amplitude, but add them directly to the inspiral amplitude model of IMRPhenomD. First we demonstrate that the inclusion of the corrections, especially in the phase, confers an overall improvement in the phase agreement between the analytical inspiral model (uncalibrated SEOBNRv2) and NR data. The model presented here is faithful, with less than 1% mismatches against a set of hybrid waveforms (except for one outlier that barely breaches this limit). The recalibrated model shows mismatches of up to ∼14% with IMRPhenomD for high mass ratios and spins. Amplitude corrections become less significant for higher mass ratios, whereas the phase corrections leave more impact—suggesting that the former is practically irrelevant for gravitational wave data analysis in Advanced LIGO (aLIGO), Virgo and KAGRA. Comparing with a set of 219 numerical relativity waveforms, we find that the median of mismatches decreases by ∼4% in aLIGO zero-detuned high power noise curve, and by ∼1.5% with a flat noise curve. This implies a modest but notable improvement in waveform accuracy. 
    more » « less