Gravitational Wave Astronomy Research Articles

Rapid Parameter Estimation for Merging Massive Black Hole Binaries Using Continuous Normalizing Flows

Abstract Detecting the coalescences of massive black hole binaries (MBHBs) is one of the primary targets for space-based gravitational wave observatories such as LISA, Taiji, and Tianqin. 
The fast and accurate parameter estimation of merging MBHBs is of great significance for the global fitting of all resolvable sources, as well as the astrophysical interpretation of gravitational wave signals.
However, such analyses usually entail significant computational costs.
To address these challenges, inspired by the latest progress in generative models, we explore the application of continuous normalizing flows (CNFs) on the parameter estimation of MBHBs. 
Specifically, we employ linear interpolation and trig interpolation methods to construct transport paths for training CNFs. 
Additionally, we creatively introduce a parameter transformation method based on the symmetry in the detector's response function. This transformation is integrated within CNFs, allowing us to train the model using a simplified dataset, and then perform parameter estimation on more general data, hence also acting as a crucial factor in improving the training speed.
In conclusion, for the first time, within a comprehensive and reasonable parameter range, we have achieved a complete and unbiased 11-dimensional rapid inference for MBHBs in the presence of astrophysical confusion noise using CNFs. In the experiments based on simulated data, our model produces posterior distributions comparable to those obtained by nested sampling.

Expected insights into Type Ia supernovae from LISA’s gravitational wave observations

The nature of progenitors of Type Ia supernovae has long been debated, primarily due to the elusiveness of the progenitor systems to traditional electromagnetic observation methods. We argue that gravitational wave observations with the upcoming Laser Interferometer Space Antenna (LISA) offer the most promising way to test one of the leading progenitor scenarios – the double-degenerate scenario, which involves a binary system of two white dwarf stars. In this study we review published results, supplementing them with additional calculations for the context of Type Ia supernovae. We discuss the fact that LISA will be able to provide a complete sample of double white dwarf Type Ia supernova progenitors with orbital periods shorter than 16–11 minutes (gravitational wave frequencies above 2–3 millihertz). Such a sample will enable a statistical validation of the double-degenerate scenario by simply counting whether LISA detects enough double white dwarf binaries to account for the measured Type Ia merger rate in Milky Way-like galaxies. Additionally, we illustrate how LISA’s capability to measure the chirp mass will set lower bounds on the primary mass, revealing whether detected double white dwarf binaries will eventually end up as a Type Ia supernova. We estimate that the expected LISA constraints on the Type Ia merger rate for the Milky Way will be 4–9%. We also discuss the potential gravitational wave signal from a Type Ia supernova assuming a double-detonation mechanism and explore how multi-messenger observations could significantly advance our understanding of these transient phenomena.

Analysis of the Principle of Gravitational Wave Detection and State- of-art Results

Over one hundred years after Einstein’s first theoretical prediction of gravitational waves, the recent observations of gravitational waves by interferometric detectors have risen excitement in the science community with the light from a new era of directly observing perturbations in spacetime itself. With the large quantity of new research in the gravitational wave detection area, this study provides a brief summary of the current achievements of the field. This paper first gives a short introduction of the field, then offers a summarize of the theories and principles of gravitational wave detection. Afterwards, it demonstrates modern detectors using the advanced LIGO detectors as an example, and presents its groundbreaking results, signal GW150914 and GW151226. Finally, the limitations and prospects of current observations are proposed. By summarizing fundamental concepts and major achievements in the field of gravitational wave detection, it is hoped to provide future students and researchers with a basic idea of the field and thus promote new ideas and advancements in the rising era of gravitational astronomy.

Modeling of Long-term Afterglow Counterparts to Gravitational Wave Events: The Full View of GRB 170817A

The arrival of gravitational wave astronomy and a growing number of time-domain-focused observatories are set to lead to an increasing number of detections of short gamma-ray bursts (GRBs) launched with a moderate inclination to Earth. Being nearby events, these are also prime candidates for very long-term follow-up campaigns and very long-baseline interferometry, which has implications for multi-messenger modeling, data analysis, and statistical inference methods applied to these sources. Here, we present a comprehensive modeling update that directly incorporates into afterglowpy astrometric observations of the GRB position, Poissonian statistics for faint sources, and modeling of a trans-relativistic population of electrons. We use the revolutionary event GW170817 to demonstrate the impact of these extensions both for the best-fit physics parameters and model selection methods that assess the statistical significance of additional late-time emission components. By including in our analysis the latest Chandra X-ray observations of GRB 170817A, we find only weak evidence (≲2σ) for a new emission component at late times, which makes for a slightly more natural fit to the centroid evolution and prediction for the external medium density.

A geometrically scalable method for manufacturing high quality factor mechanical resonators.

We present what we believe to be a novel, geometrically scalable manufacturing method for creating compact, low-resonance frequency, and high quality factor fused silica resonators. These resonators are intended to be used in inertial sensors for measuring external disturbances of sensitive physics experiments. The novel method uses direct bonding and chemical-mechanical polishing (CMP) in order to overcome the limitations of current subtractive manufacturing methods, which face prohibitive cost and complexity as material removal increases, inherently restricting the design flexibility of the resonator. We demonstrate a prototype with a test mass of only 3 g that reaches a quality factor of Q = 118 000 ± 400 at a resonance frequency of below 20 Hz. This advancement is particularly significant for future gravitational wave observatories, such as the Einstein Telescope.

Searching for triple systems unbound by supernovae

Abstract A large fraction of massive stars are found in higher order systems where the presence of a tertiary may significantly modify the system’s evolution. In particular, it can lead to increased numbers of compact object binaries and accelerate their mergers with important implications for gravitational wave observations. Using Gaia, we constrain the number of Galactic supernovae that produce unbound triples. We do this by searching 8 supernova remnants for stars with consistent Gaia parallaxes and paths intersecting near the center of the supernova remnant at a time consistent with the age of the remnant. We find no candidates for unbound triple systems. Combined with prior work, less than 11.4% of supernovae leave behind unbound triples at a 90% confidence limit. The absence of such systems limits their role in the evolution of massive stars and the formation of merging compact objects.

Test Mass Capture Control for Drag-Free Satellite Based on State-Dependent Riccati Equation Method

The drag-free satellite plays an important role in the space-based gravitational wave observatory. The capture control of test mass after release is a crucial technology that can affect the success of the mission. The test mass must be released to the center of the electrostatic suspension cage accurately. This paper presents a nonlinear dynamic model of drag-free satellites in Lagrange formalism. A capture control scheme for test mass release phase is proposed based on the state-dependent Riccati equation (SDRE) strategy. To deal with the actuator saturation problem, a nonlinear saturation model is introduced to the dynamics of satellite, while the SDRE strategy is applied to the non-affine system. The effectiveness of the proposed methodology is verified by the numerical simulation for the drag-free satellite.

Automated design of digital filters using convolutional neural networks for extracting ringdown gravitational waves

Abstract The observation of gravitational waves is expected to allow new tests of general relativity to be performed. As the gravitational wave signal is hidden by detector noise in observed data, a method to reduce noise is required to analyze the ringdown phase of gravitational wave signals. Recently, some noise reduction methods based on a neural network have been proposed; however, the results of these methods must be considered with caution because the output can contain spurious components. To overcome this limitation, in this study, we developed a neural network– based method to design optimal digital filters for extracting ringdown gravitational wave signals. In this method, no spurious components appear in the output because the digital filters reduce the noise. We conducted simulations with waveforms of gravitational waves from binary black hole coalescence and confirmed that the proposed method designs appropriate filters that reduce detector noise.

Black hole-neutron star mergers in Einstein-scalar-Gauss-Bonnet gravity

Gravitational wave observations of black hole-neutron star binaries, particularly those where the black hole has a lower mass compared to other observed systems, have the potential to place strong constraints on modifications to general relativity that arise at small curvature length scales. Here we study the dynamics of black hole-neutron star mergers in shift-symmetric Einstein-scalar-Gauss-Bonnet gravity, a representative example of such a theory, by numerically evolving the full equations of motion. We consider quasi-circular binaries with different mass-ratios that are consistent with recent gravitational wave observations, including cases with and without tidal disruption of the star, and quantify the impact of varying the coupling controlling deviations from general relativity on the gravitational wave signal and scalar radiation. We find that the main effect on the late inspiral is the accelerated frequency evolution compared to general relativity, and that—even considering Gauss-Bonnet coupling values approaching those where the theory breaks down—the impact on the merger gravitational wave signal is mild, predominately manifesting as a small change in the amplitude of the ringdown. We compare our results to current post-Newtonian calculations and find consistency throughout the inspiral. Published by the American Physical Society 2024

Analyzing One- and Two-bit Data to Reduce Memory Requirements for F -statistic-based Gravitational Wave Searches

Searches for continuous-wave gravitational radiation in data collected by modern long-baseline interferometers, such as the Laser Interferometer Gravitational-wave Observatory (LIGO), the Virgo interferometer, and the Kamioka Gravitational Wave Detector, can be memory intensive. A digitization scheme is described that reduces the 64-bit interferometer output to a one- or two-bit data stream while minimizing distortion and achieving considerable reduction in storage and input/output cost. For the representative example of the coherent, maximum-likelihood matched filter known as the F -statistic, it is found using Monte Carlo simulations that the injected signal only needs to be ≈24% stronger (for one-bit data) and ≈6.4% stronger (for two-bit data with optimal thresholds) than a 64-bit signal in order to be detected with 90% probability in Gaussian noise. The foregoing percentages do not change significantly when the signal frequency decreases secularly, or when the noise statistics are not Gaussian, as verified with LIGO Science Run 6 data.

Gravitational-wave and Gravitational-wave Memory Signatures of Core-collapse Supernovae

In this paper, we calculate the energy, signal-to-noise ratio (SNR), detection range, and angular anisotropy of the matter, matter memory, and neutrino memory gravitational-wave (GW) signatures of 21 three-dimensional initially nonrotating core-collapse supernova (CCSN) models carried to late times. We find that inferred energy, SNR, and detection range are angle-dependent quantities, and that the spread of possible energy, signal to noise, and detection ranges across all viewing angles generally increases with progenitor mass. When examining the low-frequency matter memory and neutrino memory components of the signal, we find that the neutrino memory is the most detectable component of a CCSN GW signal, and that DECIGO is best equipped to detect both matter memory and neutrino memory. Moreover, we find that the polarization angle between the h + and h × strains serves as a unique identifier of matter and neutrino memory. Finally, we develop a Galactic density- and stellar mass-weighted formalism to calculate the rate at which we can expect to detect CCSN GW signals with the Advanced Laser Interferometer Gravitational-Wave Observatory (aLIGO). When considering only the matter component of the signal, the aLIGO detection rate is around 65% of the total Galactic supernova rate, but increases to 90% when incorporating the neutrino memory component. We find that all future detectors (Einstein Telescope, Cosmic Explorer, DECIGO) will be able to detect CCSN GW signals from the entire Galaxy, and for the higher-mass progenitors even into the Local Group of galaxies.

Communicating the gravitational-wave discoveries of the LIGO-Virgo-KAGRA Collaboration

The LIGO-Virgo-KAGRA (LVK) Collaboration has made breakthrough discoveries in gravitational-wave astronomy, a new field that provides a different means of observing our Universe. Gravitational-wave discoveries are possible thanks to the work of thousands of people from across the globe working together. In this article, we discuss the range of engagement activities used to communicate LVK gravitational-wave discoveries and the stories of the people behind the science, using the activities surrounding the release of the third Gravitational-Wave Transient Catalog as a case study.

Nonlinear studies of modifications to general relativity: Comparing different approaches

Studying the dynamical, nonlinear regime of modified theories of gravity remains a theoretical challenge that limits our ability to test general relativity. Here we consider two generally applicable, but approximate, methods for treating modifications to full general relativity that have been used to study binary black hole mergers and other phenomena in this regime, and compare solutions obtained by them to those from solving the full equations of motion. The first method evolves corrections to general relativity order by order in a perturbative expansion, while the second method introduces extra dynamical fields in such a way that strong hyperbolicity is recovered. We use shift-symmetric Einstein-scalar-Gauss-Bonnet gravity as a benchmark theory to illustrate the differences between these methods for several spacetimes of physical interest. We study the formation of scalar hair about initially nonspinning black holes, the collision of black holes with scalar charge, and the inspiral and merger of binary black holes. By directly comparing predictions, we assess the extent to which those from the approximate treatments can be meaningfully confronted with gravitational wave observations. We find that the order-by-order approach cannot faithfully track the solutions when the corrections to general relativity are non-negligible. The second approach, however, can provide consistent solutions, provided the timescale over which the dynamical fields are driven to their target values is made short compared to the physical timescales. Published by the American Physical Society 2024

Probing the speed of scalar induced gravitational waves from observations

The propagation speed of gravitational waves is a fundamental issue in gravitational theory. According to general relativity, gravitational waves propagate at the speed of light. However, alternative theories of gravity propose modifications to general relativity, including variations in the speed of gravitational waves. In this paper, we investigate scalar-induced gravitational waves that propagate at speeds different from the speed of light. First, we analytically calculate the power spectrum of scalar induced gravitational waves based on the speed and spectrum of primordial curvature perturbations. We then explore several scalar power spectra, deriving corresponding fractional energy densities, including monochromatic spectrum, scale-invariant spectrum, and power-law spectrum. Finally, we constrain scalar-induced gravitational waves and evaluate the signatures of their speed from the combination of CMB+BAO and gravitational wave observations. Our numerical results clearly illustrate the influence of the speed of scalar-induced gravitational waves.

Generation of squeezed vacuum state in the millihertz frequency band

The detection of gravitational waves has ushered in a new era of observing the universe. Quantum resource advantages offer significant enhancements to the sensitivity of gravitational wave observatories. While squeezed states for ground-based gravitational wave detection have received marked attention, the generation of squeezed states suitable for mid-to-low-frequency detection has remained unexplored. To address the gap in squeezed state optical fields at ultra-low frequencies, we report on the first direct observation of a squeezed vacuum field until Fourier frequency of 4 millihertz with the quantum noise reduction of up to 8.0 dB, by the employment of a multiple noise suppression scheme. Our work provides quantum resources for future gravitational wave observatories, facilitating the development of quantum precision measurement.

Phase transition and gravitational waves in maximally symmetric composite Higgs model

In this paper we study phase transitions in a maximally symmetric composite Higgs model with next-to-minimal coset, where a pseudoscalar singlet emerges alongside the Higgs doublet. The maximal symmetry guarantees the finiteness of the radiatively generated scalar potential. We explore the scenario involving an explicit source of CP violation in the strong sector, which induces a ℤ2 asymmetric scalar potential, and consequently leads to nonzero vacuum expectation value for the singlet. Current experimental bounds from the LHC are imposed on the masses of the composite resonances, while the CP violating interactions of the pseudo Nambu-Goldstone bosons are tightly constrained from the measurements of the electric dipole moment of the electron. We compute the finite temperature corrections to the potential, incorporating the momentum-dependent form factors in the loop integrals to capture the effect of the strong dynamics. The impact of the resonances from the strong sector on the finite temperature potential are exponentially suppressed. The presence of explicit CP violation leads to strong first-order phase transition from a false vacuum to the electroweak vacuum where the pseudoscalar singlet has a non-zero vacuum expectation value. We illustrate that, as a result of such phase transitions, the production of potentially observable gravitational waves at future detectors will offer a complementary avenue to probe the composite Higgs models, distinct from collider experiments.

Inflationary Butterfly Effect: Nonperturbative Dynamics from Small-Scale Features.

For the first time, we investigate the nonperturbative dynamics of single field inflation with a departure from slow roll. Using simulations, we find that oscillatory features in the potential can drastically alter the course of inflation, with major phenomenological implications. In certain cases, the entire Universe gets trapped in a forever inflating de Sitter state. In others, only some regions get stuck in a false vacuum, offering an alternative channel for primordial black hole formation. Analogous to the flap of a butterfly, these results show that small-scale phenomena can have profound consequences on the evolution of the entire Universe. More generally, our work shows the power of simulations in the exploration of the small-scale physics of inflation, particularly in the regime relevant for gravitational-wave astronomy.

A simulation study on the sub-threshold joint gravitational wave-electromagnetic wave observation on binary neutron star mergers

Abstract The coalescence of binary neutron stars (BNS) is a prolific source of gravitational waves (GWs) and electromagnetic (EM) radiation, offering a dual observational window into the Universe. Lowering the signal-to-noise ratio (S/N) threshold is a simple and cost-effective way to enhance the detection probability of GWs from BNS mergers. In this study, we introduce a metric of the purity of joint GW and EM detections Pjoint, which is in analogue to Pastro in GW only observations. By simulating BNS merger GWs jointly detected by the HLV network and EM counterparts (kilonovae and short Gamma-ray bursts, sGRBs) with an assumed merger rate density of BNS, we generate catalogs of GW events and EM counterparts. Through this simulation, we analyze joint detection pairs, both correct and misidentified. We find the following: 1. For kilonovae, requiring Pjoint > 95% instead of $P_{\rm astro}>95~{{\%}}$ reduces the S/N from 9.2 to 8.5-8.8, allowing 5-13 additional joint detections per year and increasing the GW detection volume by 9-17%; 2. For sGRBs, requiring Pjoint > 95% instead of Pastro reduces the S/N from 9.2 to 8.1-8.5; 3. Increasing kilonova or sGRB detection capability does not improve Pjoint due to a higher rate of misidentifications. We also show that sub-threshold GW and kilonova detections can reduce the uncertainty in measuring the Hubble constant to 89-92% of its original value, and sub-threshold GW and sGRB observations can enhance the precision of constraining the speed of GWs to 88% of previously established values.

Feasibility of loop-gain tuning for general measurement systems inspired by quantum locking for DECIGO

Abstract A series of quantum locking theories have been proposed to enhance the quantum-noise-limited target sensitivity of the DECi-hertz Interferometer Gravitational wave Observatory. The quantum locking that uses a square completion optimizes the sensitivity across all frequencies. However, a substantial amount of data-series must be post-processed since the square completion is a form of signal processing technique. This paper approaches the optimal sensitivity across all frequencies from an alternative perspective: by optimizing the frequency dependence of a servo gain in a feedback loop. The optimal servo gain is formulated by comparing the alternative method with the square completion method for the same optical setup. This will be shown in general noise issues extending the framework of the quantum locking. We find that the optimal servo gain forms a non-feasible filter but has certain characteristics. We also find that the noise of the measurement signal deteriorates proportionally to the noise measured in the feedback loop when the servo gain is slightly imperfect.

Estimating false alarm rates of sub-dominant quasi-normal modes in GW190521

Abstract A major aim of gravitational wave astronomy is to test observationally the Kerr nature of black holes. The strongest such test, with minimal additional assumptions, is provided by observations of multiple ringdown modes, also known as black hole spectroscopy. For the gravitational wave merger event GW190521, we have previously claimed the detection of two ringdown modes emitted by the remnant black hole. In this paper we provide further evidence for the detection of multiple ringdown modes from this event. We analyse the recovery of simulated gravitational wave signals designed to replicate the ringdown properties of GW190521. We quantify how often our detection statistic reports strong evidence for a sub-dominant $(\ell,m,n)=(3,3,0)$ ringdown mode, even when no such mode is present in the simulated signal. We find this only occurs with a probability $\sim 0.02$, which is consistent with a Bayes factor of $56 \pm 1$ (1$\sigma$ uncertainty) found for GW190521. We also quantify our agnostic analysis of GW190521, in which no relationship is assumed between ringdown modes, and find that only 1 in 250 simulated signals without a $(3,3,0)$ mode yields a result as significant as GW190521. Conversely, we verify that when simulated signals do have an observable $(3,3,0)$ mode they consistently yield a strong evidence and significant agnostic results. We also find that constraints on deviations from the $(3,3,0)$ mode on GW190521-like signals with a $(3,3,0)$ mode are consistent with what was obtained from our previous analysis of GW190521. Our results support our previous conclusion that the gravitational wave signal from GW190521 contains an observable sub-dominant $(\ell,m,n)=(3,3,0)$ mode.