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 laser interferometer space antenna, 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.

The Impact of Astrophysical Priors on Parameter Inference for GW230529

Abstract We investigate the effects of prior selection on the inferred mass and spin parameters of the neutron star-black hole merger GW230529_181500. Specifically, we explore models motivated by astrophysical considerations, including massive binary and pulsar evolution. We examine mass and spin distributions of neutron stars constrained by radio pulsar observations, alongside black hole spin observations from previous gravitational wave detections. We show that the inferred mass distribution highly depends upon the spin prior. Specifically, under the most restrictive, binary stellar evolution models, we obtain narrower distributions of masses with a black hole mass of $4.3^{+0.1}_{-0.1}\, {\rm M}_{\odot }$and neutron star mass of $1.3^{+0.03}_{-0.03}\, {\rm M}_{\odot }$ where, somewhat surprisingly, it is the prior on component spins which has the greatest impact on the inferred mass distributions. Re-weighting using neutron star mass and spin priors from observations of radio pulsars, with black hole spins from observations of gravitational waves, yields the black hole and the neutron star masses to be $3.8^{+0.5}_{-0.6} \, M_\odot$ and $1.4^{+0.2}_{-0.1} \, M_\odot$ respectively. The sequence of compact object formation — whether the neutron star or the black hole formed first — cannot be determined at the observed signal-to-noise ratio. However, there is no evidence that the black hole was tidally spun up.

Systematic bias due to mismodeling precessing binary black hole ringdown

Accurate waveform modeling is crucial for parameter estimation in gravitational wave astronomy, impacting our understanding of source properties and the testing of general relativity. The precession of orbital and spin angular momenta in binary black hole (BBH) systems with misaligned spins presents a complex challenge for gravitational waveform modeling. Current precessing BBH waveform models employ a coprecessing frame, which precesses along with the binary. In this paper, we investigate a source of bias stemming from the mismodeling of ringdown frequency in the coprecessing frame. We find that this mismodeling of the coprecessing frame ringdown frequency introduces systematic biases in parameter estimation, for high-mass systems in particular, and in the inspiral-merger-ringdown (IMR)-consistency test of general relativity. Employing the waveform model henom, we conduct an IMR-consistency test using a Fisher matrix analysis across parameter space, as well as full injected signal parameter estimation studies. Our results show that this mismodeling particularly affects BBH systems with high-mass ratios, high-spin magnitudes, and highly inclined spins. These findings suggest inconsistencies for all waveform models which do not address this issue. Published by the American Physical Society 2024

Are all models wrong? Falsifying binary formation models in gravitational-wave astronomy using exceptional events

Abstract As the catalogue of gravitational-wave transients grows, several entries appear ‘exceptional’ within the population. Tipping the scales with a total mass of ∼150M⊙, GW190521 likely contained black holes in the pair-instability mass gap. The event GW190814, meanwhile, is unusual for its extreme mass ratio and the mass of its secondary component. A growing model-building industry has emerged to provide explanations for such exceptional events, and Bayesian model selection is frequently used to determine the most informative model. However, Bayesian methods can only take us so far. They provide no answer to the question: does our model provide an adequate explanation for exceptional events in the data? If none of the models we are testing provide an adequate explanation, then it is not enough to simply rank our existing models—we need new ones. In this paper, we introduce a method to answer this question with a frequentist p-value. We apply the method to different models that have been suggested to explain the unusually massive event GW190521: hierarchical mergers in active galactic nuclei and globular clusters. We show that some (but not all) of these models provide adequate explanations for exceptionally massive events like GW190521.

Prograde and retrograde stars in nuclear cluster merger. Evolution of the supermassive black hole binary and the host galactic nucleus

Nuclear star cluster (NSC) mergers, involving the fusion of dense stellar clusters near the centres of galaxies, play a pivotal role in shaping galactic structures. The distribution of stellar orbits has significant effects on the formation and characteristics of extreme mass ratio inspirals (EMRIs). In this study, we address the orbital distribution of stars in merging NSCs and the subsequent effects on supermassive black hole binary (SMBHB) evolution. We ran dedicated direct-summation $N$-body simulations with different initial conditions to do a detailed study of the resulting NSC after their progenitors had merged. Our findings reveal that prograde stars form a flattened structure, while retrograde stars have a more spherical distribution. The axial ratios of the prograde component vary based on the presence and mass ratio of the SMBHs. The fraction of prograde and retrograde stars depends on the merger orbital properties and the SMBH mass ratio. The interactions of retrograde stars with the SMBHB affect the eccentricity and separation evolution of the binary. Our analysis reveals a strong correlation between the angular momentum and eccentricity of the SMBH binary. This relationship could serve as a means to infer information about the stellar dynamics surrounding the binary. We find that prograde orbits are particularly close to the binary of SMBHs, a promising fact regarding EMRI production. Moreover, prograde and retrograde stars have different kinematic structures, with the prograde stars typically rotating faster than the retrograde ones. The line-of-sight velocity and velocity dispersion, as well as the velocity anisotropy of each NSC, depend on the initial merger orbital properties and SMBH mass ratios. The prograde and retrograde stars always show different behaviours. The distribution of stellar orbits and the dynamical properties of each kinematic population can potentially be used as a way to tell the properties of the parent nuclei apart, and has an important impact on expected rates of EMRIs, which will be detected by future gravitational wave observatories such as the Laser Interferometer Space Antenna (LISA).

Characterizing gravitational wave detector networks: from A♯ to cosmic explorer

Abstract Gravitational-wave observations by the Laser Interferometer Gravitational-Wave Observatory (LIGO) and Virgo have provided us a new tool to explore the Universe on all scales from nuclear physics to the cosmos and have the massive potential to further impact fundamental physics, astrophysics, and cosmology for decades to come. In this paper we have studied the science capabilities of a network of LIGO detectors when they reach their best possible sensitivity, called A#, given the infrastructure in which they exist and a new generation of observatories that are factor of 10 to 100 times more sensitive (depending on the frequency), in particular a pair of L-shaped Cosmic Explorer observatories (one 40 km and one 20 km arm length) in the US and the triangular Einstein Telescope with 10 km arms in Europe. We use science metrics derived from the top priorities of several funding agencies to characterize the science capabilities of different networks. The presence of one or two A# observatories in a network containing two or one next generation observatories, respectively, will provide good localization capabilities for facilitating multimessenger astronomy and precision measurement of the Hubble parameter. Two Cosmic Explorer observatories are indispensable for achieving precise localization of binary neutron stars, facilitating detection of electromagnetic counterparts and transforming multimessenger astronomy. Their combined operation is even more important in the detection and localization of high-redshift sources, such as binary neutron stars, beyond the star-formation peak, and primordial black hole mergers, which may occur roughly 100 million years after the Big Bang. The addition of the Einstein Telescope to a network of two Cosmic Explorer observatories is critical for accomplishing all the identified science metrics. For most metrics the triple network of next generation terrestrial observatories are a factor 100 better than what can be accomplished by a network of three A# observatories.

Conversion of 30 W laser light at 1064 nm to 20 W at 2128 nm and comparison of relative power noise

Abstract All current gravitational wave (GW) observatories operate with Nd:YAG lasers with a wavelength of 1064 nm. The sensitivity of future GW observatories could benefit significantly from changing the laser wavelength to approximately 2 μm combined with exchanging the current room temperature test mass mirrors with cryogenically cooled crystalline silicon test masses with mirror coatings from amorphous silicon and amorphous silicon nitride layers. Laser light of the order of ten watts with a low relative power noise would be required. Here we use a laboratory-built degenerate optical parametric oscillator (OPO) to convert the light from a high-power Nd:YAG laser to 2128 nm. With an input power of 30 W, we achieve an output power of 20 W, which corresponds to an external conversion efficiency of approximately 67%. We find that the relative power noise spectrum marginally increases during the wavelength conversion process. Our result is an important step in the development of low-noise light around 2 μm based on existing low-noise Nd:YAG lasers.

Toward 10-watt-level single-frequency fiber laser oscillators

High-power single-frequency laser oscillators are in great demand for some specific applications, such as quantum information and laser interferometer gravitational-wave observatories, where 10-watt or even 100-watt-level narrow-linewidth lasers with extremely low noise levels are required. In this paper, we present numerical investigations on the power scalability of single-frequency distributed Bragg reflector (DBR) ytterbium-doped phosphate fiber lasers and identify the low-quantum-defect operation approach to 10-watt single-frequency laser oscillators with the best figure of merit, defined by the ratio of the laser efficiency to the quantum defect. This research offers valuable insights for the design and development of high-power single-frequency DBR fiber laser oscillators at many other wavelengths.

Crescendo beyond the horizon: more gravitational waves from domain walls bounded by inflated cosmic strings

Abstract Gravitational-wave (GW) signals offer a unique window into the dynamics of the early universe. GWs may be generated by the topological defects produced in the early universe, which contain information on the symmetry of UV physics. We consider the case in which a two-step phase transition produces a network of domain walls bounded by cosmic strings. Specifically, we focus on the case in which there is a hierarchy in the symmetry-breaking scales, and a period of inflation pushes the cosmic string generated in the first phase transition outside the horizon before the second phase transition. We show that the GW signal from the evolution and collapse of this string-wall network has a unique spectrum, and the resulting signal strength can be sizeable. In particular, depending on the model parameters, the resulting signal can show up in a broad range of frequencies and can be discovered by a multitude of future probes, including the pulsar timing arrays and space- and ground-based GW observatories. As an example that naturally gives rise to this scenario, we present a model with the first phase transition followed by a brief period of thermal inflation driven by the field responsible for the second stage of symmetry breaking. The model can be embedded into a supersymmetric setup, which provides a natural realization of this scenario. In this case, the successful detection of the peak of the GW spectrum probes the soft supersymmetry breaking scale and the wall tension.

Frozen stars: Black hole mimickers sourced by a string fluid

The frozen star is a nonsingular, ultracompact object that, to an external observer, looks exactly like a Schwarzschild black hole, but with a different interior geometry and matter composition. The frozen star needs to be sourced by an extremely anisotropic fluid, for which the sum of the radial pressure and energy density is either vanishing or perturbatively small. Here, we show that this matter can be identified with the string fluid resulting from the decay of an unstable D-brane or a brane-antibrane system at the end of open-string tachyon condensation. The string fluid corresponds to flux tubes emanating from the center and ending at the Schwarzschild radius of the star. The effective Lagrangian for this fluid can be recast into a Born-Infeld form. When the fluid Lagrangian is coupled to that of Einstein’s gravity, the static, spherically symmetric solutions of the equations of motion are shown to be the same as those describing the frozen star model. Frozen stars can therefore be viewed as a gravitationally backreacted BIons model or that of gravitationally confined flux tubes. The Born-Infeld Lagrangian provides a complete set of equations that describe the dynamics of the frozen star in a generic state, which is not necessarily static nor spherically symmetric. Our framework therefore should allow a fully dynamical and out-of-equilibrium description of the frozen star model which will be relevant to gravitational waves observations of mergers of astrophysical black holes. Published by the American Physical Society 2024

Mapping eccentricity evolutions between numerical relativity and effective-one-body gravitational waveforms

Orbital eccentricity in compact binaries is considered to be a key tracer of their astrophysical origin, and can be inferred from gravitational-wave observations due to its imprint on the emitted signal. For a robust measurement, accurate waveform models are needed. However, ambiguities in the definition of eccentricity can obfuscate the physical meaning and result in seemingly discrepant measurements. In this work we present a suite of 28 new numerical relativity simulations of eccentric, aligned-spin binary black holes with mass ratios between 1 and 6 and initial post-Newtonian eccentricities between 0.05 and 0.3. We then develop a robust pipeline for measuring the eccentricity evolution as a function of frequency from gravitational-wave observables that is applicable even to signals that span at least ≳7 orbits. We assess the reliability of our procedure and quantify its robustness under different assumptions on the data. Using the eccentricity measured at the first apastron, we initialize effective-one-body waveforms and quantify how the precision in the eccentricity measurement, and therefore the choice of the initial conditions, impacts the agreement with the numerical data. We find that even small deviations in the initial eccentricity can lead to non-negligible differences in the phase and amplitude of the waveforms. However, we demonstrate that we can reliably map the eccentricities between the simulation data and analytic models, which is crucial for robustly building eccentric hybrid waveforms, and to improve the accuracy of eccentric waveform models in the strong-field regime. Published by the American Physical Society 2024

Bounds on the charge of the graviton using gravitational wave observations

If the graviton possesses a non-zero charge qg , gravitational waves (GW) originating from astrophysical sources would experience an additional time delay due to intergalactic magnetic fields. This would result in a modification of the phase evolution of the observed GW signal similar to the effect induced by a massive graviton. As a result, we can reinterpret the most recent upper limits on the graviton's mass as constraints on the joint mass-charge parameter space, finding |qg |/e < 3 × 10-34 where e represents the charge of an electron. Additionally, we illustrate that a charged graviton would introduce a constant phase difference in the gravitational waves detected by two spatially separated GW detectors due to the Aharonov-Bohm effect. Using the non-observation of such a phase difference for the GW event GW190814, we establish a mass-independent constraint |qg |/e < 2 × 10-26. To the best of our knowledge, our results constitute the first-ever bounds on the charge of the graviton. We also discuss various caveats involved in our measurements and prospects for strengthening these bounds with future GW observations.

Physical Models for the Astrophysical Population of Black Holes: Application to the Bump in the Mass Distribution of Gravitational-wave Sources

Abstract Gravitational-wave observations of binary black holes have revealed unexpected structure in the black hole mass distribution. Previous studies employ physically motivated phenomenological models and infer the parameters that control the features of the mass distribution that are allowed in their model, associating the constraints on those parameters with their physical motivations a posteriori. In this work, we take an alternative approach in which we introduce a model parameterizing the underlying stellar and core-collapse physics and obtaining the remnant black hole distribution as a derived by-product. In doing so, we constrain the stellar physics necessary to explain the astrophysical distribution of black hole properties under a given model. We apply this to the mapping between initial mass and remnant black hole mass, accounting for mass-dependent mass loss using a simple parameterized description. Allowing the parameters of the initial mass–remnant mass relationship to evolve with redshift permits correlated and physically reasonable changes to features in the mass function. We find that the current data are consistent with no redshift evolution in the core–remnant mass relationship, but place only weak constraints on the change of these parameters. This procedure can be applied to modeling any physical process underlying the astrophysical distribution. We illustrate this by applying our model to the pulsational pair instability supernova (PPISN) process, previously proposed as an explanation for the observed excess of black holes at ∼35 M ⊙. Placing constraints on the reaction rates necessary to explain the PPISN parameters, we concur with previous results in the literature that the peak observed at ∼35 M ⊙ is unlikely to be a signature from the PPISN process as presently understood.

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.