Search

In this Letter, we present the design and performance of the frequency-dependent squeezed vacuum source that will be used for the broadband quantum noise reduction of the Advanced Virgo Plus gravitational-wave detector in the upcoming observation run. The frequency-dependent squeezed field is generated by a phase rotation of a frequency-independent squeezed state through a 285 m long, high-finesse, near-detuned optical resonator. With about 8.5 dB of generated squeezing, up to 5.6 dB of quantum noise suppression has been measured at high frequency while close to the filter cavity resonance frequency, the intracavity losses limit this value to about 2 dB. Frequency-dependent squeezing is produced with a rotation frequency stability of about 6 Hz rms, which is maintained over the long term. The achieved results fulfill the frequency dependent squeezed vacuum source requirements for Advanced Virgo Plus. With the current squeezing source, considering also the estimated squeezing degradation induced by the interferometer, we expect a reduction of the quantum shot noise and radiation pressure noise of up to 4.5 dB and 2 dB, respectively.

We present the first directed searches for long-transient and continuous gravitational waves from ultralight vector boson clouds around known black holes (BHs). We use LIGO data from the first part of the fourth LIGO-Virgo-KAGRA observing run. The searches target two distinct types of BHs and use two new semicoherent methods: hidden Markov model (HMM) tracking for the remnant BHs of the mergers GW230814_230901 and GW231123_135430 (referred to as GW230814 and GW231123 in this study) and a dedicated method using the Band Sampled Data (BSD) framework for the galactic BH in the Cygnus X-1 binary system. Without finding evidence of a signal from vector bosons in the data, we estimate the mass range that can be constrained. For the HMM searches targeting the remnants from GW231123 and GW230814, we disfavor vector boson masses in the ranges [0.94, 1.08] and $[2.75,3.28]\times {10}^{-13}\mathrm{eV}$, respectively, at 30% confidence, assuming a 1% false alarm probability. Although these searches are only marginally sensitive to signals from merger remnants at relatively large distances, future observations are expected to yield more stringent constraints with high confidence. For the BSD search targeting the BH in Cygnus X-1, we exclude vector boson masses in the range $[0.85,1.59]\times {10}^{-13}\mathrm{eV}$ at 95% confidence, assuming an initial BH spin larger than 0.5.

The angular distribution of gravitational-wave power from persistent sources may exhibit anisotropies arising from the large-scale structure of the Universe. This motivates directional searches for astrophysical and cosmological gravitational-wave backgrounds, as well as continuous-wave emitters. We present results of such a search using data from the first observing run through the first portion of the fourth observing run of the LIGO-Virgo-KAGRA Collaborations. We apply gravitational-wave radiometer techniques to generate skymaps and search for both narrowband and broadband persistent gravitational-wave sources. Additionally, we use spherical harmonic decomposition to probe spatially extended sources. No evidence of persistent gravitational-wave signals is found, and we set the most stringent constraints to date on such emissions. For narrowband point sources, our sensitivity estimate to effective strain amplitude lies in the range $(0.03–8.4)\times {10}^{-24}$ across all-sky and frequency range ($20-160$) Hz. For targeted sources—Scorpius X-1, SN 1987A, the Galactic Center, Terzan 5, and NGC 6397—we constrain the strain amplitude with best limits ranging from $\sim 1.1\times {10}^{-25}$ to $6.5\times {10}^{-24}$. For persistent broadband sources, we constrain the gravitational-wave flux $F_{\alpha ,\widehat{n}}^{95\%,\mathrm{UL}}(25\mathrm{Hz})<(0.008–5.5)\times {10}^{-8}\mathrm{erg}{\mathrm{cm}}^{-2}{\mathrm{s}}^{-1}{\mathrm{Hz}}^{-1}$, depending on the sky direction $\widehat{n}$ and spectral index $\alpha =0,2/3$, 3. Finally, for extended sources, we place upper limits on the angular power spectrum $C_{\ell }^{1/2}<(0.63–17)\times {10}^{-10}{\mathrm{sr}}^{-1}$.

We present a search for gravitational waves from inspiraling, planetary-mass ultracompact binaries using data from the first part of the fourth observing run of LIGO, Virgo, and KAGRA. Finding no evidence of such systems, we determine the maximum distance reach for such objects and their merger rate densities. Then, we identify classes of primordial black hole mass distributions for which these rate limits can be translated into relevant constraints on the mass distribution of primordial black holes, assuming that they are all composed of dark matter, in the mass range $[{10}^{-6},{10}^{-3}]M_{\odot }$. Our constraints for the class of primordial black hole mass distribution functions $f(m)$ we consider reach down to $f(m)\sim 0.1$. They are consistent with existing microlensing results in the planetary-mass range, provide a complementary probe to subsolar mass objects, and are publicly available [1].

The recent discovery by Advanced LIGO and Advanced Virgo of a gravitational wave signal from a binary neutron star inspiral has enabled tests of general relativity (GR) with this new type of source. This source, for the first time, permits tests of strong-field dynamics of compact binaries in the presence of matter. In this Letter, we place constraints on the dipole radiation and possible deviations from GR in the post-Newtonian coefficients that govern the inspiral regime. Bounds on modified dispersion of gravitational waves are obtained; in combination with information from the observed electromagnetic counterpart we can also constrain effects due to large extra dimensions. Finally, the polarization content of the gravitational wave signal is studied. The results of all tests performed here show good agreement with GR.

Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1{)}_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1{)}_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation timescale of DM.

The gravitational-wave signal GW250114 was observed by the two LIGO detectors with a network matched-filter signal-to-noise ratio of 80. The signal was emitted by the coalescence of two black holes with near-equal masses $m_1={33.6}_{-0.8}^{+1.2}{M}_{\odot }$ and $m_2={32.2}_{-1.3}^{+0.8}{M}_{\odot }$, and small spins ${\chi }_{1,2}\le 0.26$ (90% credibility) and negligible eccentricity $e\le 0.03$. Postmerger data excluding the peak region are consistent with the dominant quadrupolar $(\ell =|m|=2)$ mode of a Kerr black hole and its first overtone. We constrain the modes’ frequencies to $\pm 30\%$ of the Kerr spectrum, providing a test of the remnant’s Kerr nature. We also examine Hawking’s area law, also known as the second law of black hole mechanics, which states that the total area of the black hole event horizons cannot decrease with time. A range of analyses that exclude up to five of the strongest merger cycles confirm that the remnant area is larger than the sum of the initial areas to high credibility.

We present a search for short-duration gravitational-wave transients in data from the first eight months of Advanced LIGO-Virgo-KAGRA’s fourth observing run, denoted O4a. We use four analyses which are sensitive to a wide range of potential signals lasting up to a few seconds in the 16–4096 Hz band. Excluding binary black hole merger candidates that were already identified by low-latency analyses, we find no statistically significant evidence for other gravitational-wave transients. We measure the sensitivity of the search for representative signals, including sine-Gaussians, Gaussian pulses, and white-noise bursts with different frequencies and durations, adopting a false alarm rate of 1 per 100 years as detection threshold. Depending on signal type, we find improvements over previous searches by factors of 2 to 10 in terms of sensitivity to strain amplitude and of 90% confidence upper limit on the rate density of sources. We also evaluate a variety of core-collapse supernova models and find that, for some models, the search could have detected gravitational waves from stellar core-collapse throughout the Milky Way. Finally, we consider neutron star $f$-modes associated with pulsar glitches and find that, assuming a source similar to the Vela Pulsar, the search could have detected a gravitational-wave signal from a glitch with fractional frequency change as small as $\sim 2$ to $6\times {10}^{-5}$ depending on the neutron star mass.

The binary black hole signal GW250114, the loudest gravitational wave detected to date, offers a unique opportunity to test Einstein’s general relativity (GR) in the high-velocity, strong-gravity regime and probe whether the remnant conforms to the Kerr metric. Upon perturbation, black holes emit a spectrum of damped sinusoids with specific, complex frequencies. Our analysis of the postmerger signal shows that at least two quasinormal modes are required to explain the data, with the most damped remaining statistically significant for about one cycle. We probe the remnant’s Kerr nature by constraining the spectroscopic pattern of the dominant quadrupolar ($\ell =m=2$) mode and its first overtone to match the Kerr prediction to tens of percent at multiple postpeak times. The measured mode amplitudes and phases agree with a numerical-relativity simulation having parameters close to GW250114. By fitting a parametrized waveform that incorporates the full inspiral-merger-ringdown sequence, we constrain the fundamental $(\ell =m=4)$ mode to tens of percent and bound the quadrupolar frequency to within a few percent of the GR prediction. We perform a suite of tests—spanning inspiral, merger, and ringdown—finding constraints that are comparable to, and in some cases 2–3 times more stringent than those obtained by combining dozens of events in the fourth Gravitational-Wave Transient Catalog. These results constitute the most stringent single-event verification of GR and the Kerr nature of black holes to date, and outline the power of black-hole spectroscopy for future gravitational-wave observations.

We introduce a new approach to modified gravity which generalizes the recently proposed hybrid metric-Palatini gravity. The gravitational action is taken to depend on a general function of both the metric and Palatini curvature scalars. The dynamical equivalence with a nonminimally coupled biscalar field gravitational theory is proved. The evolution of cosmological solutions is studied using dynamical systems techniques.

We present results from the search for an isotropic gravitational-wave background using Advanced LIGO and Advanced Virgo data from O1 through O4a, the first part of the fourth observing run. This background is the accumulated signal from unresolved sources throughout cosmic history and encodes information about the merger history of compact binaries throughout the Universe, as well as exotic physics and potentially primordial processes from the early cosmos. Our cross-correlation analysis reveals no statistically significant background signal, enabling us to constrain several theoretical scenarios. For compact binary coalescences which approximately follow a $2/3$ power-law spectrum, we constrain the fractional energy density to ${\mathrm{\Omega }}_{\mathrm{GW}}(25\mathrm{Hz})\le 2.0\times {10}^{-9}$ (95% credibility), a factor of 1.7 improvement over previous results. Scale-invariant backgrounds are constrained to ${\mathrm{\Omega }}_{\mathrm{GW}}(25\mathrm{Hz})\le 2.8\times {10}^{-9}$, representing a $2.1\times$ sensitivity gain. We also place new limits on alternative gravity theories predicting nonstandard polarization modes and confirm that terrestrial magnetic noise sources remain below our detection threshold. Combining these spectral limits with population models for Gravitational-Wave Transient Catalog (GWTC)-4, the latest gravitational-wave event catalog, we find our constraints remain above predicted merger backgrounds but are approaching detectability. The joint analysis combining the background limits shown here with GWTC-4 enables improved inference of the binary black hole merger rate evolution across cosmic time. Furthermore, employing GWTC-4 inference results and standard modeling choices, we estimate that the total background arising from compact binary coalescences is ${\mathrm{\Omega }}_{\mathrm{CBC}}(25\mathrm{Hz})={0.9}_{-0.5}^{+1.1}\times {10}^{-9}$ at 90% confidence, where the largest contribution is due to binary black holes only, ${\mathrm{\Omega }}_{\mathrm{BBH}}(25\mathrm{Hz})={0.8}_{-0.5}^{+1.1}\times {10}^{-9}$.

We present a search of data from the first part of the fourth observing run of LIGO-Virgo-KAGRA for three kinds of dark matter—dilatons (spin-0), dark photons (spin-1), and tensor bosons (spin-2)—using three independent methods, each extended to search for all three candidates. Each dark matter candidate could interact with different standard-model particles in the instruments, causing unique differential strains on the interferometers. While we do not find any evidence for a signal, we place the most stringent upper limits to date on the couplings of each dark matter model to the interferometers. For scalars with masses between $[4\times {10}^{-14},1.5\times {10}^{-13}]\mathrm{eV}$ that couple to photons or electrons, our constraints on the scalar-to-electron or photon coupling improve upon those from the third observing run by 1 order of magnitude, with the tightest limit of $\sim {10}^{-20}{\mathrm{GeV}}^{-1}$ at a mass of $\sim 2\times {10}^{-13}\mathrm{eV}$. For vectors with masses between $[7\times {10}^{-13},8.47\times {10}^{-12}]\mathrm{eV}$ that couple to baryons, our constraints on the dark-photon or baryon coupling supersede those from MICROSCOPE and Eöt-Wash by 1 to 2 orders of magnitude, reaching a minimum of $\sim 5\times {10}^{-24}$ at a mass of $\sim {10}^{-12}\mathrm{eV}$. For tensors with masses of $[4\times {10}^{-14},8.47\times {10}^{-12}]\mathrm{eV}$ (the full mass range analyzed) that couple via a Yukawa interaction, our constraints, which are the first obtained from gravitational-wave interferometers, surpass those from fifth-force experiments by 4 to 5 orders of magnitude, achieving a limit as low as $\sim 8\times {10}^{-9}$ at $\sim 2\times {10}^{-13}\mathrm{eV}$. Our results show that each method has substantially different sensitivity to each dark matter candidate, highlighting the importance of a multimethod approach to search for new physics. Moreover, they demonstrate that gravitational-wave interferometers can simultaneously probe multiple ultralight dark matter models and place unprecedented constraints on tensor fields arising in a broad range of modified-gravity theories.

We present an all-sky search for long-duration gravitational waves (GWs) from the first part of the LIGO-Virgo-KAGRA fourth observing run (O4), called O4a and comprising data taken between May 24, 2023, and January 16, 2024. The GW signals targeted by this search are the so-called “long-duration” ($\gtrsim 1\mathrm{s}$) transients expected from a variety of astrophysical processes, including nonaxisymmetric deformations in magnetars or eccentric binary coalescences. We make minimal assumptions on the emitted GW waveforms in terms of morphologies and durations. Overall, our search targets signals with durations of $\sim 1–1000\mathrm{s}$ and frequency content in the range 16–2048 Hz. In the absence of significant detections, we report the sensitivity limits of our search in terms of root-sum-square signal amplitude ($h_{\mathrm{rss}}$) of reference waveforms. These limits improve upon the results from the third LIGO-Virgo-KAGRA observing run (O3) by about 30% on average. Moreover, this analysis demonstrates substantial progress in our ability to search for long-duration GW signals owing to enhancements in pipeline detection efficiencies. As detector sensitivities continue to advance and observational runs grow longer, unmodeled long-duration searches will increasingly be able to explore a range of compelling astrophysical scenarios involving neutron stars and black holes.

GW240925 and GW250207 are two loud gravitational-wave signals from binary black hole coalescences observed with network signal-to-noise ratios $\sim 32$ and $\sim 69$, respectively, by the LIGO Hanford-LIGO Livingston-Virgo network. Gravitational-wave signals from coalescing binaries have characteristic phase and amplitude evolution predicted by general relativity. These signal waveforms, together with measured instrumental calibration uncertainties, are used to infer source parameters. However, for sufficiently loud detections, it is possible to constrain the calibration of the detectors directly using the signals themselves. We present the first informative astrophysical measurements of gravitational-wave detector calibration. For GW240925, we verify the inference of Hanford calibration from the astrophysical signal through cross-checks with known calibration errors obtained from in situ measurements. At the time of GW250207, the Hanford detector was not fully stabilized, leading to elevated calibration uncertainties; thus, astrophysical calibration is essential to obtain accurate data and to enable source localization. These well-localized, high signal-to-noise observations have the potential to offer precise measurements of source properties, stringent tests of general relativity, and informative dark siren measurements, provided that calibration uncertainties are properly incorporated. As detector sensitivity improves, astrophysical calibration will become an increasingly valuable complement to in situ calibration measurements. Obtaining accurate calibration will be essential for precision gravitational-wave science.

A new variational approach for general relativity and modified theories of gravity is presented. In addition to the metric tensor, two independent affine connections enter the action as dynamical variables. In the matter action the dependence upon one of the connections is left completely unspecified. When the variation is applied to the Einstein-Hilbert action the Einstein field equations are recovered. However when applied to $f(R)$ and scalar-tensor theories, it yields gravitational field equations which differ from their equivalents obtained with a metric or Palatini variation and reduce to the former ones only when no connections appear in the matter action.

A significant fraction of stars are members of gravitationally bound hierarchies containing three or more components. Almost all low-mass stars in binaries with periods shorter than three days are part of a hierarchical system. We therefore anticipate that a large fraction of compact galactic binaries detected by the Laser Interferometer Space Antenna will be members of hierarchical triple or quadruple system. The acceleration imparted by the hierarchical companions can be detected in the gravitational wave signal for outer periods as large as 100 years. For systems with periods that are shorter than, or comparable to, the mission lifetime, it will be possible to measure the period and eccentricity of the outer orbit. Laser Interferometer Space Antenna observations of hierarchical stellar systems will provide insight into stellar evolution, including the role that Lidov-Kozai oscillations play in driving systems toward merger.

Multifrequency gravitational wave (GW) observations are useful probes of the formation processes of coalescing stellar-mass binary black holes (BBHs). We discuss the phase drift in the GW inspiral waveform of the merging BBH caused by its center-of-mass acceleration. The acceleration strongly depends on the location where a BBH forms within a galaxy, allowing observations of the early inspiral phase of Laser Interferometer Gravitational Wave Observatory (LIGO)-like BBH mergers by the Laser Interferometer Space Antenna (LISA) to test the formation mechanism. In particular, BBHs formed in dense nuclear star clusters or via compact accretion disks around a nuclear supermassive black hole in active galactic nuclei would suffer strong acceleration, and produce large phase drifts measurable by LISA. The host galaxies of the coalescing BBHs in these scenarios can also be uniquely identified in the LISA error volume, without electromagnetic counterparts. A nondetection of phase drifts would rule out or constrain the contribution of the nuclear formation channels to the stellar-mass BBH population.

In the first-order formalism of gravitational theories, the spacetime connection is considered as an independent variable to vary together with the metric. However, the metric still generates its Levi-Civita connection that turns out to determine the geodesics of matter. Recently, “hybrid” gravity theories have been introduced by constructing actions involving both the independent Palatini connection and the metric Levi-Civita connection. In this study a method is developed to analyze the field content of such theories, in particular to determine whether the propagating degrees of freedom are ghosts or tachyons. New types of second-, fourth- and sixth-order derivative gravity theories are investigated and the so-called $f(X)$ theories are singled out as a viable class of “hybrid” extensions of general relativity. These are the theories in which the corrections to Einstein’s theory are written in terms of the deviations from the usual trace equation $X=R-{\kappa }^{2}T$.

Space-borne gravitational-wave (GW) detectors observing at millihertz and decihertz frequencies are expected to detect large numbers of quasimonochromatic signals. The first and second time derivative of the GW frequency (${\dot{f}}_0$ and ${\ddot{f}}_0$) can be measured for the most favorable sources and used to look for negative post-Newtonian corrections, which can be induced by the source’s environment or modifications of general relativity. We present an analytical, Fisher-matrix-based approach to estimate how precisely such corrections can be constrained. We use this method to estimate the bounds attainable on the time evolution of the gravitational constant $G(t)$ with different classes of quasimonochromatic sources observable with LISA and DECIGO, two representative space-borne detectors for millihertz and decihertz GW frequencies. We find that the most constraining source among a simulated population of LISA galactic binaries could yield $\dot{G}/G_0\lesssim {10}^{-6}{\mathrm{yr}}^{-1}$, while the best currently known verification binary will reach $\dot{G}/G_0\lesssim {10}^{-4}{\mathrm{yr}}^{-1}$. We also perform Monte Carlo simulations using quasimonochromatic waveforms to check the validity of our Fisher-matrix approach, as well as inspiralling waveforms to analyse binaries that do not satisfy the quasimonochromatic assumption. We find that our analytical Fisher matrix produces good order-of-magnitude constraints even for sources well beyond its regime of validity. Monte Carlo investigations also show that chirping stellar-mass compact binaries detected by DECIGO-like detectors at cosmological distances of tens of Mpc can yield constraints as tight as $\dot{G}/G_0\lesssim {10}^{-11}{\mathrm{yr}}^{-1}$.

Laser Interferometer Space Antenna (LISA) observations of massive black hole binaries (MBHBs) will provide long duration inspiral signals with high signal-to-noise ratio (SNR) data, ideal for testing general relativity (GR) in the strong-field and relativistic regime regime. We present an extension of the flexible theory-independent (FTI) framework, adapted to gravitational waves (GWs) from MBHBs observed with LISA, to perform parametrized inspiral tests of GR. This approach introduces generic deviations to the post-Newtonian (PN) coefficients of the frequency-domain GW phase while accounting for the time- and frequency-dependent instrument response, thus effectively identifying potential deviations from GR by constraining modifications to the PN phasing formula. Complementary analyses using Fisher matrix and full Bayesian approaches confirm that LISA observations could improve constraints on agnostic, scale-independent deviations from GR by at least two orders of magnitude compared to the most recent LIGO-Virgo-KAGRA measurements. Since LISA’s sensitivity to different GW phases—inspiral, merger, and ringdown—varies across the MBHB parameter space with masses between ${10}^4$ and ${10}^7{M}_{\odot }$, the optimal regime for testing agnostic deviations is not known a priori. Our results illustrate how the strength of these constraints depends significantly on both the total mass and the SNR, reflecting the trade-off between inspiral and merger-ringdown contributions to the observed signal. We also investigate the interplay between inspiral-only versus inspiral-merger-ringdown analyses in constraining these inspiral deviation parameters. This work contributes to the development of robust tests of GR with LISA, enhancing our ability to probe the nature of gravity and BHs with GW observations.