Gravitational Wave Astronomy Research Articles

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 centre 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 per cent of supernovae leave behind unbound triples at a 90 per cent 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

Stable configuration design for libration point gravitational wave observatory

The Sun-Earth L2 libration point configuration is one of the options for space-based gravitational wave detection. Long-term configuration stability is crucial for high-precision measurements, challenged by the strong nonlinear dynamics of the Sun-Earth three-body system. This paper proposes an efficient design method and determines the feasible parameter domain for the libration point gravitational wave observatory. First, the dynamic model for the libration point configuration is established, and the stability indexes are defined. The sensitive parameters that affect the relative geometric configuration are discussed and the phase angle is found to be the key factor. Then, an efficient design method is proposed, and the procedure is divided into two steps. The phase angle of the Earth phase offset orbit and the libration point configuration are optimized successively. Finally, the proposed method is applied to the LAGRANGE mission concept. The results show that the three stability indexes decrease by 59%, 42% and 23%, respectively. Moreover, a mapping between configuration parameters and stability indexes is established. The feasible parameter domain for the stable libration point configuration is discussed. The feasible amplitudes domain in the x and z directions of the libration point orbit should be less than 6200 km and 42000 km, respectively, to guarantee configuration stability. This research could provide a reference for the stable design and implementation of gravitational wave detection missions utilizing libration point configuration in the future.

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.

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.

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.

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.

Concept study of a storage ring-based gravitational wave observatory: Gravitational wave strain and synchrotron radiation noise

This work for the first time addresses the feasibility of measuring millihertz gravitational waves (mHz GWs) with a storage ring-based detector. While this overall challenge consists of several partial problems, here we focus solely on quantifying design limitations imposed by the kinetic energy and radiated power of circulating ions at relativistic velocities. We propose an experiment based on the measurement of the time-of-flight signal of an ion chain. One of the dominant noise sources inherent to the measurement principle for such a GW detector is the shot noise of the emitted synchrotron radiation. We compute the noise amplitude of arrival time signals obtained by analytical estimates and simulations of ions with different masses and velocities circulating in a storage ring with the circumference of the Large Hadron Collider. Thereby, we show that our experiment design could reduce the noise amplitude due to the synchrotron radiation in the frequency range 10−4–10−2 Hz to one or two orders of magnitude below the expected GW signals from of astrophysical sources, such as supermassive binary black holes or extreme mass-ratio inspirals. Other key requirements for building a working storage ring-based GW detector include the generation and acceleration of heavy ion chains with the required energy resolution, their injection and continued storage, as well as the detection method to be used for the determination of the particle arrival time. However, these are not the focus of the work presented here, in which we instead concentrate on the definition of a working principle in terms of ion type, kinetic energy, and ring design, which will later serve as a starting point when addressing a more complete experimental setup. Published by the American Physical Society 2024

Black hole—neutron star binary mergers: the impact of stellar compactness

Abstract Recent gravitational wave (GW) observations include possible detections of black hole—neutron star binary mergers. As with binary black hole mergers, numerical simulations help characterize the sources. For binary systems with neutron star components, the simulations help to predict the imprint of tidal deformations and disruptions on the GW signals. In a previous study, we investigated how the mass of the black hole has an impact on the disruption of the neutron star and, as a consequence, on the shape of the GWs emitted. We extend these results to study the effects of varying the compactness of the neutron star. We consider neutron star compactness in the 0.113–0.2 range for binaries with mass ratios of 3 and 5. As the compactness and the mass ratio increase, the binary system behaves during the late inspiral and merger more like a black hole binary. For the cases with the least compact neutron star, the GWs emitted, in terms of mismatches, are the most distinguishable from those by a binary black hole. The disruption of the star significantly suppresses the kicks on the final black hole. The disruption also affects, although not dramatically, the spin of the final black hole. Lastly, for neutron stars with low compactness, the quasi-normal ringing of the black hole after the merger does not show a clean quasi-normal ringing because of the late accretion of debris from the neutron star.

Formation of GW230529 from Isolated Binary Evolution

In this paper, we explore the formation of the mass-gap black hole-neutron star (mgBHNS) merger detected in gravitational wave (GW) event, i.e., GW230529, from the isolated binary evolution channel, and study potential signatures of its electromagnetic counterparts. By adopting the “delayed” supernova prescription and reasonable model realizations, our population synthesis simulation results can simultaneously match the rate densities of mgBHNS and total BHNS mergers inferred from the population analyses, along with the population distribution of the BH mass in BHNS mergers reported by the LIGO–Virgo–KAGRA Collaboration. Because GW230529 contributes significantly to the inferred mgBHNS rate densities, we suggest that GW230529 can be explained through the isolated binary evolution channel. Considering the AP4 (DD2) equation of state, the probability that GW230529 can make tidal disruption is 12.8% (63.2%). If GW230529 is a disrupted event, its kilonova peak apparent magnitude is predicted ∼23–24 mag, and hence, can be detected by the present survey projects and Large Synoptic Survey Telescope. Since GW230529 could be an off-axis event inferred from the GW observation, its associated gamma-ray burst (GRB) might be too dim to be observed by γ-ray detectors, interpreting the lack of GRB observations. Our study suggests the existence of mgBHNS mergers formed through the isolated binary evolution channel due to the discovery of GW230529, indicating that BHNS mergers are still likely to be multimessenger sources that emit GWs, GRBs, and kilonovae. Although mgBHNS mergers account for ∼50% of the cosmological BHNS population, we find that ≳90% of disrupted BHNS mergers are expected to originate from mgBHNS mergers.