We present a complete pipeline for detecting and characterising gravitational waves (GWs) produced by the inspiral of stellar-mass binary black holes in data from the Laser Interferometer Space Antenna (LISA). The analysis framework relies on an efficient time-frequency implementation of an adaptive semi-coherent detection statistic, which we show to be robust against non-stationary noise and the presence of gaps of varying duration and cadence. The search is able to detect signals down to a coherent signal-to-noise ratio $\approx$ 11 over the full parameter space of black holes with spins aligned to the orbital angular momentum and orbital eccentricity $\leq$ 0.01 when deployed on the 2-year-long LISA Data Challenge Yorsh. The search can be run within a day using $\approx$ 40 GPUs. The techniques presented here have wider applications in GW astronomy, in particular the search for extreme-mass-ratio inspirals in LISA data.

We propose a gravitational wave (GW) detector based on ultrastable optical cavities enabling the detection of GW signals in the mostly unexplored 10-5-1 Hz frequency band. We illustrate the working principle of the detector and discuss that several classes of GW sources, both of astrophysical and cosmological origin, may be within the detection range of this instrument. Our work suggests that terrestrial GW detection in the milli-Hz frequency range is potentially within reach with current technology.

Evidence has emerged for a stochastic signal correlated among 67 pulsars within the 15 yr pulsar-timing data set compiled by the NANOGrav collaboration. Similar signals have been found in data from the European, Indian, Parkes, and Chinese pulsar timing arrays. This signal has been interpreted as indicative of the presence of a nanohertz stochastic gravitational-wave background (GWB). To explore the internal consistency of this result, we investigate how the recovered signal strength changes as we remove the pulsars one by one from the data set. We calculate the signal strength using the (noise-marginalized) optimal statistic, a frequentist metric designed to measure the correlated excess power in the residuals of the arrival times of the radio pulses. We identify several features emerging from this analysis that were initially unexpected. The significance of these features, however, can only be assessed by comparing the real data to synthetic data sets. After conducting identical analyses on simulated data sets, we do not find anything inconsistent with the presence of a stochastic GWB in the NANOGrav 15 yr data. The methodologies developed here can offer additional tools for application to future, more sensitive data sets. While this analysis provides an internal consistency check of the NANOGrav results, it does not eliminate the necessity for additional investigations that could identify potential systematics or uncover unmodeled physical phenomena in the data.

Progress in gravitational-wave (GW) astronomy depends upon having sensitive detectors with good data quality. Since the end of the Laser Interferometer Gravitational-Wave Observatory-Virgo-KAGRA third Observing run in March 2020, detector-characterization efforts have lead to increased sensitivity of the detectors, swifter validation of GW candidates and improved tools used for data-quality products. In this article, we discuss these efforts in detail and their impact on our ability to detect and study GWs. These include the multiple instrumental investigations that led to reduction in transient noise, along with the work to improve software tools used to examine the detectors data-quality. We end with a brief discussion on the role and requirements of detector characterization as the sensitivity of our detectors further improves in the future Observing runs.

The Laser Interferometer Space Antenna (LISA) will detect 100 galactic binary systems comprised of black holes (BHs) and neutron stars (NSs). Identifying the nature of the constituents of these binaries as BHs or NSs, and distinguishing them from ∼ 10^4 detected double white dwarfs will be challenging. In the absence of any other information, the inferred values of the component masses can be used to classify the nature of these objects. However, short-period galactic binaries ∼ 10^7 - 10^3 yr from coalescence produce a quasi-monochromatic signal which carries little information about their masses. We generate synthetic LISA data sets containing gravitational waves (GWs) from galactic binary BHs, binary NSs and BHNSs drawn from an astrophysically realistic population produced through the isolated binary evolution channel. We process the data with an end-to-end Bayesian inference pipeline to explore the accuracy with which the individual component masses can be measured. We find that for ≈ 10% - 50% of the detected systems LISA will be able to measure the individual component masses by measuring the orbital eccentricity, periapse precession frequency, and GW induced frequency derivative. Typical fractional mass errors are ≈ 1% - 100% (depending on the specific value of the source parameters), which will enable in many circumstances the classification of the objects as BHs or NSs. For these binaries, LISA will also be able to determine their 3-dimensional position in the Milky Way. The LISA-detected sample of these double compact objects will provide new information about the galactic population of BHs and NSs, the star formation history of the Milky Way and the astrophysical processes leading to the formation of these systems. However, for a significant fraction of the LISA-detected binaries the nature of their constituent objects may remain unclear.

We present the results of a search for gravitational-wave transients associated with core-collapse supernova SN 2023ixf, which was observed in the galaxy Messier 101 via optical emission on 2023 May 19, during the LIGO–Virgo–KAGRA 15th Engineering Run. We define a five-day on-source window during which an accompanying gravitational-wave signal may have occurred. No gravitational waves have been identified in data when at least two gravitational-wave observatories were operating, which covered ∼14% of this five-day window. We report the search detection efficiency for various possible gravitational-wave emission models. Considering the distance to M101 (6.7 Mpc), we derive constraints on the gravitational-wave emission mechanism of core-collapse supernovae across a broad frequency spectrum, ranging from 50 Hz to 2 kHz, where we assume the gravitational-wave emission occurred when coincident data are available in the on-source window. Considering an ellipsoid model for a rotating proto-neutron star, our search is sensitive to gravitational-wave energy 1 × 10 −4 M ⊙ c 2 and luminosity 2.6 × 10 −4 M ⊙ c 2 s −1 for a source emitting at 82 Hz. These constraints are around an order of magnitude more stringent than those obtained so far with gravitational-wave data. The constraint on the ellipticity of the proto-neutron star that is formed is as low as 1.08, at frequencies above 1200 Hz, surpassing past results.

We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received with low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wave Transient Catalogs (GWTC-3). Targeted searches were carried out on the entire GW sample using the maximum-likelihood Non-imaging Transient Reconstruction and Temporal Search pipeline on the BAT data made available via the GUANO infrastructure. We do not detect any significant electromagnetic emission that is temporally and spatially coincident with any of the GW candidates. We report flux upper limits in the 15-350 keV band as a function of sky position for all the catalog candidates. For GW candidates where the Swift-BAT false alarm rate is less than 10(-3) Hz, we compute the GW-BAT joint false alarm rate. Finally, the derived Swift-BAT upper limits are used to infer constraints on the putative electromagnetic emission associated with binary black hole mergers.

Compact objects observed in gravitational-wave astronomy so far always come in pairs and never individually. Identifying the two components of a binary system is a delicate operation that is often taken for granted. The labeling procedure (i.e., which is object "1" and which is object "2") effectively acts as systematics, or, equivalently, an unspecified prior, in gravitational-wave data inference. The common approach is to label the objects solely by their masses, on a sample-by-sample basis. We show that object identification can instead be tackled using the posterior distribution as a whole. We frame the problem in terms of constrained clustering-a flavor of semisupervised machine learning-and find that unfolding the labeling systematics can significantly impact, and arguably improve, our interpretation of the data. In particular, the precision of black-hole spin measurements improves by up to 50%, multimodalities and tails tend to disappear, posteriors become closer to Gaussian distributions, and the identification of the nature of the object (i.e., black hole vs neutron star) is facilitated. We estimate that about 10% of the LIGO/Virgo posterior samples are affected by this relabeling, i.e., they might have been attributed to the other compact object in the observed binaries.

We demonstrate an end-to-end technique for observing and characterizing massive black hole binary signals before they merge with the LISA space-based gravitational-wave observatory. Our method uses a zero-latency whitening filter, originally designed for rapidly observing compact binary mergers in ground-based observatories, to be able to observe signals with no additional latency due to filter length. We show that with minimal computational cost, we are able to reliably observe signals as early as 14 days premerger as long as the signal has accrued a signal-to-noise ratio of at least 8 in the LISA data. We also demonstrate that this method can be used to characterize the source properties, providing early estimates of the source's merger time, chirp mass, and sky localization. Early observation and characterization of massive black holes is crucial to enable the possibility of rapid multimessenger observations, and to ensure that LISA can enter a protected operating period when the merger signal arrives.

Ultralight axionlike particles (ALPs) can be a viable solution to the dark matter problem. The scalar field associated with ALPs, coupled to the electromagnetic field acts as an active birefringent medium, altering the polarization properties of light through which it propagates. In particular, oscillations of the axionic field induce monochromatic variations of the plane of linearly polarized radiation of astrophysical signals. The radio emission of millisecond pulsars provides an excellent tool to search for such manifestations, given their high fractional linear polarization and negligible fluctuations of their polarization properties. We have searched for evidence of ALPs in polarimetry measurements of pulsars collected and preprocessed for the European Pulsar Timing Array (EPTA) campaign. Focusing on the twelve brightest sources in linear polarization, we searched for an astrophysical signal from axions using both frequentist and Bayesian statistical frameworks. For the frequentist analysis, which uses Lomb-Scargle periodograms at its core, no statistically significant signal has been found. The model used for the Bayesian analysis has been adjusted to accommodate multiple deterministic systematics that may be present in the data. A statistically significant signal has been found in the dataset of multiple pulsars with common frequencies between 10-8 and 2 x 10-8 Hz, which can most likely be explained by the residual Faraday rotation in the terrestrial ionosphere. Strong bounds on the coupling constant ga gamma, in the same ballpark as other searches, have been obtained in the mass range between 6 x 10-24 and 5 x 10-21 eV. We conclude by discussing the problems that can limit the sensitivity of our search for ultralight axions in the polarimetry data of pulsars, and possible ways to resolve them.

High-quality test mass substrates play essential roles in laser interferometric gravitational wave detectors. Inhomogeneous birefringence distribution in test mass substrates, however, can degrade the sensitivity of the detector by introducing the optical loss and disturbing the interferometer controls. In this paper, we present a heterodyne polarimetry method that enables in situ birefringence characterizations, hence diagnosing the gravitational wave interferometer. We experimentally demonstrate the proposed method with a tabletop setup. We also discuss its applicability to current and future gravitational wave detectors and the detectable limit.

On May 24, 2023, the Advanced Laser Interferometer Gravitational-Wave Observatory (LIGO), joined by the Advanced Virgo and KAGRA detectors, began the fourth observing run for a two-year-long dedicated search for gravitational waves. The LIGO Hanford and Livingston detectors have achieved an unprecedented sensitivity to gravitational waves, with an angle-averaged median range to binary neutron star mergers of 152 and 160 Mpc, and duty cycles of 65.0% and 71.2%, respectively, with a coincident duty cycle of 52.6%. The maximum range achieved by the LIGO Hanford detector is 165 Mpc and the LIGO Livingston detector 177 Mpc, both achieved during the second part of the fourth observing run. For the fourth run, the quantum-limited sensitivity of the detectors was increased significantly due to the higher intracavity power from laser system upgrades and replacement of core optics, and from the addition of a 300 m filter cavity to provide the squeezed light with a frequency-dependent squeezing angle, part of the Athorn upgrade program. Altogether, the Athorn upgrades led to reduced detector-wide losses for the squeezed vacuum states of light which, alongside the filter cavity, enabled broadband quantum noise reduction of up to 5.2 dB at the Hanford observatory and 6.1 dB at the Livingston observatory. Improvements to sensors and actuators as well as significant controls commissioning increased low frequency sensitivity. This paper details these instrumental upgrades, analyzes the noise sources that limit detector sensitivity, and describes the commissioning challenges of the fourth observing run.

On 2022 February 15–16, multiple spacecraft measured one of the most intense solar energetic particle (SEP) events observed so far in Solar Cycle 25. This study provides an overview of interesting observations made by multiple spacecraft during this event. Parker Solar Probe (PSP) and BepiColombo were close to each other at 0.34–0.37 au (a radial separation of ∼0.03 au) as they were impacted by the flank of the associated coronal mass ejection (CME). At about 100° in the retrograde direction and 1.5 au away from the Sun, the radiation detector on board the Curiosity surface rover observed the largest ground-level enhancement on Mars since surface measurements began. At intermediate distances (0.7–1.0 au), the presence of stream interaction regions (SIRs) during the SEP arrival time provides additional complexities regarding the analysis of the distinct contributions of CME-driven versus SIR-driven events in observations by spacecraft such as Solar Orbiter and STEREO-A, and by near-Earth spacecraft like ACE, SOHO, and WIND. The proximity of PSP and BepiColombo also enables us to directly compare their measurements and perform cross-calibration for the energetic particle instruments on board the two spacecraft. Our analysis indicates that energetic proton measurements from BepiColombo and PSP are in reasonable agreement with each other to within a factor of ∼1.35. Finally, this study introduces the various ongoing efforts that will collectively improve our understanding of this impactful, widespread SEP event.

We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers.

The Heisenberg uncertainty principle dictates that the position and momentum of an object cannot be simultaneously measured with arbitrary precision, giving rise to an apparent limitation known as the standard quantum limit (SQL). Gravitational-wave detectors use photons to continuously measure the positions of freely falling mirrors and so are affected by the SQL. We investigated the performance of the Laser Interferometer Gravitational-Wave Observatory (LIGO) after the experimental realization of frequency-dependent squeezing designed to surpass the SQL. For the LIGO Livingston detector, we found that the upgrade reduces quantum noise below the SQL by a maximum of three decibels between 35 and 75 hertz while achieving a broadband sensitivity improvement, increasing the overall detector sensitivity during astrophysical observations.

The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling supermassive black hole binaries (SMBHBs); inflation, phase transitions, cosmic strings and tensor mode generation by non-linear evolution of scalar perturbations in the early Universe; oscillations of the Galactic potential in the presence of ultra-light dark matter (ULDM). At the current stage of emerging evidence, it is impossible to discriminate among the different origins. Therefore, in this paper, we consider each process separately, and investigate the implications of the signal under the hypothesis that it is generated by that specific process. We find that the signal is consistent with a cosmic population of inspiralling SMBHBs, and its relatively high amplitude can be used to place constraints on binary merger timescales and the SMBH-host galaxy scaling relations. If this origin is confirmed, this is the first direct evidence that SMBHBs merge in nature, adding an important observational piece to the puzzle of structure formation and galaxy evolution. As for early Universe processes, the measurement would place tight constraints on the cosmic string tension and on the level of turbulence developed by first-order phase transitions. Other processes would require non-standard scenarios, such as a blue-tilted inflationary spectrum or an excess in the primordial spectrum of scalar perturbations at large wavenumbers. Finally, a ULDM origin of the detected signal is disfavoured, which leads to direct constraints on the abundance of ULDM in our Galaxy.

Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level.

We present the results of a search for continuous gravitational wave signals (CGWs) in the second data release (DR2) of the European Pulsar Timing Array (EPTA) Collaboration. The most significant candidate event from this search has a gravitational wave frequency of 4-5 nHz. Such a signal could be generated by a supermassive black hole binary (SMBHB) in the local Universe. We present the results of a follow-up analysis of this candidate using both Bayesian and frequentist methods. The Bayesian analysis gives a Bayes factor of 4 in favour of the presence of the CGW over a common uncorrelated noise process. In contrast, the frequentist analysis estimates the p-value of the candidate to be < 1%, also assuming the presence of common uncorrelated red noise. However, comparing a model that includes both a CGW and a gravitational wave background (GWB) to a GWB only, the Bayes factor in favour of the CGW model is only 0.7. Therefore, we cannot conclusively determine the origin of the observed feature, nor can we rule it out as a CGW source. We present results of simulations that demonstrate that data containing a weak gravitational wave background can be misinterpreted as data including a CGW and vice versa, providing two plausible explanations for the EPTA DR2 data. Further investigations combining data from all PTA collaborations will be needed to reveal the true origin of this feature.

The magnetar SGR 1935+2154 is the only known Galactic source of fast radio bursts (FRBs). FRBs from SGR 1935+2154 were first detected by the Canadian Hydrogen Intensity Mapping Experiment (CHIME)/FRB and the Survey for Transient Astronomical Radio Emission 2 in 2020 April, after the conclusion of the LIGO, Virgo, and KAGRA Collaborations' O3 observing run. Here, we analyze four periods of gravitational wave (GW) data from the GEO600 detector coincident with four periods of FRB activity detected by CHIME/FRB, as well as X-ray glitches and X-ray bursts detected by NICER and NuSTAR close to the time of one of the FRBs. We do not detect any significant GW emission from any of the events. Instead, using a short-duration GW search (for bursts <= 1 s) we derive 50% (90%) upper limits of 10(48) (10(49)) erg for GWs at 300 Hz and 10(49) (10(50)) erg at 2 kHz, and constrain the GW-to-radio energy ratio to <= 10(14)-10(16). We also derive upper limits from a long-duration search for bursts with durations between 1 and 10 s. These represent the strictest upper limits on concurrent GW emission from FRBs.

We report the observation of a coalescing compact binary with component masses 2.5-4.5 M-circle dot and 1.2-2.0 M-circle dot (all measurements quoted at the 90% credible level). The gravitational-wave signal GW230529_181500 was observed during the fourth observing run of the LIGO-Virgo-KAGRA detector network on 2023 May 29 by the LIGO Livingston observatory. The primary component of the source has a mass less than 5 M-circle dot at 99% credibility. We cannot definitively determine from gravitational-wave data alone whether either component of the source is a neutron star or a black hole. However, given existing estimates of the maximum neutron star mass, we find the most probable interpretation of the source to be the coalescence of a neutron star with a black hole that has a mass between the most massive neutron stars and the least massive black holes observed in the Galaxy. We provisionally estimate a merger rate density of 55(-47)(+127) Gpc-3yr(-1) for compact binary coalescences with properties similar to the source of GW230529_181500; assuming that the source is a neutron star-black hole merger, GW230529_181500-like sources may make up the majority of neutron star-black hole coalescences. The discovery of this system implies an increase in the expected rate of neutron star-black hole mergers with electromagnetic counterparts and provides further evidence for compact objects existing within the purported lower mass gap.