Gravitational Wave Detectors Research Articles

A rapid multi-modal parameter estimation technique for LISA

Abstract The Laser Interferometer Space Antenna (LISA) will observe gravitational-wave signals from a wide range of sources, including massive black hole binaries. Although numerous techniques have been developed to perform Bayesian inference for LISA, they are often computationally expensive; analyses often take at least $\sim 1$ month on a single CPU, even when using accelerated techniques. Not only does this make it difficult to concurrently analyse more than one gravitational-wave signal, it also makes it challenging to rapidly produce parameter estimates for possible electromagnetic follow-up campaigns. simple-pe was recently developed to produce rapid parameter estimates for gravitational-wave signals observed with ground-based gravitational-wave detectors. In this work, we extend simple-pe to produce rapid parameter estimates for LISA sources, including the effects of higher order multipole moments. We show that simple-pe infers the source properties of massive black hole binaries in zero-noise at least $\sim 100\times$ faster than existing techniques; $\sim 12$ hours on a single CPU. We further demonstrate that simple-pe can be applied before existing Bayesian techniques to mitigate biases in multi-modal parameter estimation analyses of MBHBs.

Suppression of Clock-Jitter Noise and Laser Phase Noise in Arm Locking

Abstract Arm-locking technique has been a focus of attention as one of the means to suppress the laser phase noise in space-based gravitational wave (GW) detector. 
 The main idea of the arm-locking technique is to transfer the stability of the detector arm length to laser frequency by introducing a feedback control loop. 
 Generally, laser phase noise will be suppressed by an amount similar to the magnitude of the controller gain. 
 However, on the one hand, clock-jitter noise and optical bench motion noise, as the noise floor of the arm-locking technique, need to be suppressed. 
 On the other hand, limited by the Doppler frequency pulling, the gain of the controller generally cannot be too large.
 It means that even if we do not consider clock-jitter noise and optical bench motion noise, it is difficult to suppress laser phase noise below the noise floor only by arm-locking technique. 
 In this work, we combine self-referenced optical frequency combs (OFC) and arm-locking technique to generate clock signals that are coherently referenced to the closed-loop laser beams, so that the clock-jitter noise is also suppressed by about the level of controller gain.
 We conduct a simulation on the above configuration, and the results show that the performance of the arm-locking is limited by clock-jitter noise in the low-frequency band. 
 To address the issue of insufficient laser phase noise suppression by the arm-locking technique, we further investigate Time-delay interferometry (TDI) combinations under outputs of arbitrary arm-locking configurations. 
 We obtain the equations for eliminating laser phase noise.
 To ensure that the TDI combinations can directly operate in the time domain, we derive a restricted solution space by assuming a specific form for the solutions.

Blazar S 0528+134 is possibly the most powerful emitter in the Universe, including in the range of gravitational waves

A new method of determining the parameters of close binary systems of supermassive black holes (SMBHs) and the level of gravitational waves (GWs) on the Earth's surface are proposed. Data are presented from long-term monitoring of possibly the most powerful emitter in the Universe, S 0528+134, at five frequencies in the radio range from 4.8 GHz to 37 GHz, obtained by the RT-22 radio telescope of the Radio Astronomy Laboratory CrAO in Simeiz, the RT-26 radio telescope at Michigan Observatory, and the 40 m radio telescope of the Owens Valley Radio Observatory. The dynamics of powerful flares that have occurred in the object since its discovery in a sky survey in 1970 were considered. The main physical characteristics of binary SMBHs located in the central regions of the system were obtained. These data were used to find the masses of the SMBH companions, the parameters of their orbits, the energy reserve of the system, and the lifetime of the object before the SMBHs' merger. The level of GWs on the Earth's surface was determined and the possibility of their detection by International Pulsar Timing Array (IPTA) GW detectors was considered.

Strong-lensing cosmography using third-generation gravitational-wave detectors

Abstract We present a detailed exposition of a statistical method for estimating cosmological parameters from the observation of a large number of strongly lensed binary-black-hole (BBH) mergers observable by next (third) generation (XG) gravitational-wave (GW) detectors. This method, first presented in [1], compares the observed number of strongly lensed GW events and their time delay distribution (between lensed images) with observed events to infer cosmological parameters. We show that the precision of the estimation of the cosmological parameters does not have a strong dependance on the assumed BBH redshift distribution model. Using the large number of unlensed mergers, XG detectors are expected to measure the BBH redshift distribution with sufficient precision for the cosmological inference. However, a biased inference of the BBH redshift distribution will bias the estimation of cosmological parameters. An incorrect model for the distribution of lens properties can also lead to a biased cosmological inference. However, Bayesian model selection can assist in selecting the right model from a set of available parametric models for the lens distribution. We also present a way to incorporate the effect of contamination in the data due to the limited efficiency of lensing identification methods, so that it will not bias the cosmological inference.

Swiftly Chasing Gravitational Waves across the Sky in Real Time

We introduce a new capability of the Neil Gehrels Swift Observatory, dubbed “continuous commanding,” that achieves 10 s latency response time on orbit to unscheduled target-of-opportunity requests received on the ground. We show that this will allow Swift to respond to premerger (early-warning) gravitational-wave (GW) detections, rapidly slewing the Burst Alert Telescope (BAT) across the sky to place the GW origin in the BAT field of view at or before merger time. This will dramatically increase the GW/gamma-ray burst (GRB) codetection rate and enable prompt arcminute localization of a neutron star merger. We simulate the full Swift response to a GW early-warning alert, including input sky maps produced at different early-warning times, a complete model of the Swift attitude control system, and a full accounting of the latency between the GW detectors and the spacecraft. 60 s of early warning can double the rate of a prompt GRB detection with arcminute localization, and 140 s guarantees observation anywhere on the unocculted sky, even with localization areas ≫1000 deg2. While 140 s is beyond current GW detector sensitivities, 30–70 s is achievable today. We show that the detection yield is now limited by the latency of LIGO/Virgo cyberinfrastructure and motivate a focus on its reduction. Continuous commanding has been integrated as a general capability of Swift, significantly increasing its versatility in response to the growing demands of time-domain astrophysics. We demonstrate this potential on an externally triggered fast radio burst (FRB), slewing 81° across the sky, and collecting X-ray and UV photons from the source position <150 s after the trigger was received from the Canadian Hydrogen Intensity Mapping Experiment, thereby setting the earliest and deepest such constraints on high-energy activity from nonrepeating FRBs. The Swift Team invites the community to consider and propose novel scientific applications of ultra-low-latency UV, X-ray, and gamma-ray observations.

Evaluating the gravitational wave detectability of globular clusters and the Magellanic Clouds for LISA

We used the stellar evolution code bpass and the gravitational wave (GW) simulation code legwork to simulate populations of compact binaries that may be detected by the future space-based GW detector LISA. Specifically, we simulate the Magellanic Clouds and binary populations mimicking several globular clusters, neglecting dynamical effects. We find that a handful of sources should be detectable in each of the Magellanic Clouds, but for globular clusters the amount of detectable sources will likely be less than one each. We compared our results to earlier research and find that our predicted numbers are several dozen times lower than both the results from calculations that used the stellar evolution code bse and take dynamical effects into account, and results from calculations that used the stellar evolution code SeBa for the Magellanic Clouds. Earlier research that compared bpass models for GW sources in the Galactic disk with bse models found a similarly sized discrepancy. We determine that this discrepancy is caused by differences between the stellar evolution codes, particularly in the treatment of mass transfer and common-envelope events in binaries: in bpass mass transfer is more likely to be stable and tends to lead to less orbital shrinkage in the common-envelope phase than in other codes. This difference results in fewer compact binaries with periods short enough to be detected by LISA existing in the bpass population. For globular clusters, we conclude that the impact of dynamical effects is uncertain based on the literature, but the differences in stellar evolution have an effect of a factor of 20 to 40 on the number of detectable binaries.

Black hole spectroscopy with ground-based atom interferometer and space-based laser interferometer gravitational wave detectors

Gravitational wave (GW) detection allows us to test general relativity in entirely new regimes. A prominent role takes the detection of quasi-normal modes (QNMs), which are emitted after the merger of a binary black hole (BBH) when the highly distorted remnant emits GWs to become a regular Kerr black hole (BH). The BH uniqueness theorems of Kerr black hole solutions in general relativity imply that the frequencies and damping times of QNMs are determined solely by the mass and spin of the remnant BH. Therefore, detecting QNMs offers a unique way to probe the nature of the remnant BH and to test general relativity. We study the detection of a merging BBH in the intermediate-mass range, where the inspiral–merger phase is detected by space-based laser interferometer detectors TianQin and LISA, while the ringdown is detected by the ground-based atom interferometer (AI) observatory AION. The analysis of the ringdown is done using the regular broadband mode of AI detectors as well as the resonant mode optimizing it to the frequencies of the QNMs predicted from the inspiral–merger phase. We find that the regular broadband mode allows constraining the parameters of the BBH with relative errors of the order 10−1 and below from the ringdown. Moreover, for a variety of systems considered, the frequencies and the damping times of the QNMs can be determined with relative errors below 0.1 and 0.2, respectively. We further find that using the resonant mode can improve the parameter estimation for the BBH from the ringdown by a factor of up to three. Utilizing the resonant mode significantly limits the detection of the frequency of the QNMs but improves the detection error of the damping times by around two orders of magnitude.

Listening to the Universe with next generation ground-based gravitational-wave detectors

Listening to the Universe with next generation ground-based gravitational-wave detectors

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.

An upper limit on the spins of merging binary black holes formed through isolated binary evolution

As the sensitivity of ground-based gravitational wave detectors progressively increases, observations of black hole mergers will provide us with the joint distribution of their masses and spins. This will be a critical benchmark to validate different formation scenarios. Merging binary black holes formed through the evolution of isolated binary systems require both components to be stripped of their hydrogen envelopes before core-collapse. The rotation rates of such stripped stars are constrained by the critical rotation limit at their surface, including its deviation from the Keplerian value owing to the outward force provided by radiation. This sets a restriction on their angular momentum content at core-collapse. We aim to determine if this restriction plays a role in the spins of binary black hole mergers. We used detailed calculations of stripped stars with the MESA code at low metallicities ($Z=Z_ $Z_ and $Z_ to determine the dimensionless spins of black holes produced by critically rotating stellar progenitors. To study how such progenitors can arise, we considered their formation through chemically homogeneous evolution (CHE) in binary stars. We used a semi-analytical model to study the physical processes that determine the final angular momentum of CHE binaries, and compared our results against available population synthesis models that rely on detailed binary evolution calculations. We find that above black hole masses of $ 25M_ the dimensionless spin parameter of critically rotating stripped stars ($a=Jc/(GM^2)$) is below unity. This results in an exclusion region at high chirp masses and effective spins that cannot be populated by isolated binary evolution. CHE can produce binaries where both black holes hit this limit, producing a pileup at the boundary of the excluded region. High-spin black holes arise from very low-metallicity CHE systems with short delay times, which merge at higher redshifts. On the other hand, the contribution of CHE to merging binary black holes detected in the third observing run of the LVK collaboration is expected to be dominated by systems with low spins ($ eff &lt;0.5$) that merge near redshift zero. Owing to its higher projected sensitivity and runtime, the fourth observing run of the LVK collaboration can potentially place constraints on the high-spin population and the existence of a limit set by critical rotation.

Uncertainty of the white dwarf astrophysical gravitational wave background

The astrophysical gravitational wave background (AGWB) is a stochastic gravitational wave (GW) signal emitted by different populations of in-spiralling binary systems containing compact objects throughout the Universe . In the frequency range between 10$^ $ and 10$^ $ hertz (Hz), it will be detected by future space-based gravitational wave detectors, such as Laser Interferometer Space Antenna (LISA). In a recent work, we concluded that the white dwarf (WD) contribution to the AGWB dominates that of black holes (BHs) and neutron stars (NSs) We aim to investigate the uncertainties of the WD AGWB that arise from the use of different stellar metallicities, star formation rate density (SFRD) models, and binary evolution models. We used the code we previously developed to determine the WD component of the AGWB. We used a metallicity-dependent SFRD based on an earlier work to construct five different SFRD models. We used four different population models based on a range of common-envelope treatments and six different metallicities for each model. For all possible combinations, the WD component of the AGWB is dominant over other populations of compact objects. The effects of metallicity and population model are less significant than the effect of a (metallicity dependent) SFRD model. We find a range of about a factor of 5 in the level of the WD AGWB around a value of $ WD $ at 1 mHz and a shape that is weakly dependent on the model. We find the uncertainty for the WD component of the AGWB to be about a factor of 5. We note that there are other uncertainties that have an effect on this signal as well. We discuss whether the turnover of the WD AGWB at 10 mHz will be detectable by LISA and find it to be likely. We confirm our previous findings asserting that the WD component of the AGWB dominates over other populations, in particular, BHs.

Comparison of arm cavity optical losses for the two wavelengths of the Einstein telescope gravitational wave detector

Abstract A new generation of gravitational wave detectors is currently being designed with the likely use of a different laser wavelength compared to current instruments. The estimation of the optical losses for this new wavelength is particularly relevant to derive the detector sensitivity and also to anticipate the optical performances of future instruments. In this article, we measured the absorption and angle-resolved scattering of several mirror samples in order to compare optical losses at a wavelength of 1064 and 1550 nm. In addition, we have carried out simulations of the Einstein Telescope arm cavities at 1064 and 1550 nm taking into account losses due to surface low-spatial frequency flatness. Our results suggest that optical losses as measured at 1064 nm are about twice as large as those at 1550 nm as predicted with a simple model.

Measuring neutron star radius with second and third generation gravitational wave detector networks

Abstract The next generation of ground-based interferometric gravitational wave detectors will observe mergers of black holes and neutron stars throughout cosmic time. A large number of the binary neutron star merger events will be observed with extreme high fidelity, and will provide stringent constraints on the equation of state of nuclear matter. In this paper, we investigate the systematic improvement in the measurability of the equation of state with increase in detector sensitivity by combining constraints obtained on the radius of a $1.4 \, \mathrm{M}_{\odot}$ neutron star from a simulated source population. Since the measurability of the equation of state depends on its stiffness, we consider a range of realistic equations of state that span the current observational constraints. We show that a single 40km Cosmic Explorer detector can pin down the neutron star radius for a soft, medium and stiff equation of state to an accuracy of 10m within a decade, whereas the current generation of ground-based detectors like the Advanced LIGO-Virgo network would take $\mathcal{O}(10^5)$ years to do so for a soft equation of state. &amp;#xD;

Approaches of frequency-dependent squeezing for the low frequency detector of the Einstein Telescope

The quantum noise in gravitational-wave detectors can be suppressed in a broadband by frequency-dependent squeezing. It usually requires one large scale filter cavity and even two, for example in the low frequency detector of Einstein Telescope, which is a detuned dual recycling Fabry-Perot Michelson interferometer. In this paper, we study the feasibility of replacing two filter cavities with a coupled-cavity, aiming to reduce the optical losses with less number of optics. It turns out this approach is only theoretically valid, however, the required parameters of the optics do not support practical implementation, which is consistent with the results in Jones [Implications of the quantum noise target for the Einstein Telescope infrastructure design, ]. Furthermore, we investigate the viability of utilizing Einstein-Podolsky-Rosen (EPR) squeezing to eliminate either one or two filter cavities in the Einstein Telescope. It turns out EPR squeezing would allow us to eliminate one filter cavity, and can potentially improve the detector sensitivity with the allowance of higher input squeezing level, benefiting from the longer length of the arm cavity which serves as one of the filter cavities for frequency-dependent squeezing. Published by the American Physical Society 2024

Probing interacting dark sector with the next generation of gravitational-wave detectors

Probing interacting dark sector with the next generation of gravitational-wave detectors

The Virgo Newtonian calibration system for the O4 observing run

Abstract After initial tests performed during previous observing runs, a Newtonian Calibrator (NCal) system was developed and installed on the Virgo gravitational wave detector for the O4 observing run. This system, which is continuously operated, provides the absolute calibration of Virgo for this run. Its 1-σ uncertainty of 0.17% on the amplitudes of the injected signals is better than that obtained with other calibration techniques like the photon calibrator (PCal). This paper presents this NCal system and details the different sources of uncertainties.

Suppression of the Effect of Parametric Instability in the LIGO Voyager Gravitational Wave Detector

The number of unstable combinations of elastic and Stokes optical modes of the LIGO Voyager gravitational wave detector, frequencies and spatial distributions of displacement vectors of elastic modes of mirrors were calculated. The values of overlap factors for unstable modes up to ninth-order optical modes are calculated, taking into account the azimuthal condition of parametric instability. An analysis of the influence of the temperature dependence of the Young’s modulus of the mirror material on the number of unstable modes in the Fabry-Perot resonator was performed.

Projections of the uncertainty on the compact binary population background using popstock

The LIGO-Virgo-KAGRA collaboration has announced the detection to date of almost 100 binary black holes that have been used in several studies to infer the features of the underlying binary black hole population. From these objects it is possible to predict the overall gravitational-wave (GW) fractional energy density contributed by black holes throughout the Universe, and thus estimate the gravitational-wave background (GWB) spectrum emitted in the current GW detector band. These predictions are fundamental in our forecasts for background detection and characterisation, with both present and future instruments. The uncertainties in the inferred population strongly impact the predicted energy spectrum, and in this paper we present a new flexible method to quickly calculate the energy spectrum for varying black hole population features, such as the mass spectrum and redshift distribution. We have implemented this method in an open-access package popstock and extensively tested its capabilities. Using popstock we investigated how uncertainties in these distributions impact our detection capabilities, and present several caveats for background estimation. In particular, we find that the standard assumption that the background signal follows a two-thirds power law at low frequencies is both waveform and mass-model dependent, and that the power-law signal is likely shallower than previously modelled, given the current waveform and population knowledge.

Gravitational-wave detectors perform better with squeezing

Gravitational-wave detectors perform better with squeezing

Online decoupling of the time-varying longitudinal feedback loops for improved performance in Advanced Virgo Plus*

Abstract Gravitational Wave detectors require low-noise sensors combined with high-performance feedback loops to maximize the detector sensitivity in the low-frequency detection range. Some feedback loops in the detector are strongly coupled and their coupling varies over time, which is inherent to the optical configuration and the optical readout scheme used to obtain error signals for the feedback loops. This paper presents a control-based approach to deal with the time-varying interaction among three of the longitudinal degrees of freedom in the Advanced Virgo Plus detector, using the pre-existing infrastructure of the detector. An intuitive Multiple-Input Multiple-Output stability criterion is presented that illustrates how the time-varying coupling affects the stability of the feedback loops, as well as a method that identifies the coupling online and updates a decoupling matrix accordingly. Experiment results of the proposed method on the Advanced Virgo Plus detector illustrate continuous minimization of the coupling over time. Using the derived stability criterion, it is furthermore shown that the coupling is sufficiently minimized such that the interaction terms can be neglected in the control design, resulting in an improved controller performance.