The detection of a stochastic gravitational-wave background (SGWB) is a primary science objective for the Laser Interferometer Space Antenna (LISA). However, extracting these signals is difficult because both the signal and the instrumental noise are stochastic and overlapping in the millihertz band. In this work, we present a Bayesian framework for the joint estimation of LISA noise and SGWB signals. Our approach models the LISA instrumental noise using flexible log-penalized splines, employing a roughness penalty to prevent overfitting while maintaining computational efficiency. For the SGWB, we compare a power-law model with a spline-based model and study how the choice of signal model and noise prior affects signal recovery and detection. Using simulated LISA data, we find that the power-law model gives tighter estimates when the signal follows the assumed shape. However, it fails to recover a localized spectral feature that is not described by a power law, causing the signal to be absorbed by the instrumental-noise spline. The fully spline-based model is less restrictive and successfully recovers such features. We also find that stronger prior information about the test-mass noise helps reduce the degeneracy between the noise and SGWB models at low frequencies, improving signal detection. These results reflect a single trade-off: added model flexibility reduces sensitivity when the assumed signal shape is correct, and prevents bias when it is not.

This paper addresses the important problem of estimating the noise spectral density of next-generation gravitational-wave detectors, such as LISA and the Einstein Telescope (ET), where cross-channel correlations must be accounted for to avoid biased parameter estimation of gravitational-wave signals. Unlike approaches that estimate test-mass and optical-metrology-system noise separately at the single-link level and then map them to the Time-Delay Interferometry (TDI) channels through known transfer functions, we develop a Bayesian nonparametric method that directly estimates the spectral density matrix of the XYZ channels, thereby accommodating additional sources of uncertainty. Our approach combines a flexible matrix-gamma process prior on the matrix-valued coefficients of a Bernstein polynomial basis expansion with a blocked multivariate Whittle likelihood. The prior guarantees Hermitian positive definiteness of the spectral estimate at every frequency. To avoid reversible-jump methods, we use an adaptive Markov chain Monte Carlo (MCMC) algorithm for posterior sampling. The proposed framework can also be used to correct misspecified parametric noise models. Results from a simulation study and simulated correlated-noise data for both LISA and ET demonstrate the effectiveness of the proposed method.

Cross-correlation searches for gravitational-wave backgrounds depend on the geometrical configuration (physical separation and relative orientation) of the detectors comprising the network. Applying standard techniques to a few simple examples, we illustrate how the relative orientation of a pair of Earth-based L-shaped laser interferometers can drastically impact the detectability of both isotropic and anisotropic gravitational-wave backgrounds.

The LIGO-Virgo-KAGRA collaboration has reported the results for searches for sub-solar mass components in compact binary coalescence during observing run O4a. No detection candidates were identified, but the most significant seven triggers were reported. We present the results of Bayesian parameter inference on these triggers. Five of the triggers show agreement between the Bayesian parameter estimation and the search pipeline trigger. Our results show that three of the triggers may contain a possible sub-solar mass component. Parameter estimation indicates that the other two events, if real, would be neutron star - black hole binaries. The remaining two triggers do not provide informative parameter estimation. We also study three O4a compact binary coalescence triggers, and one O3 trigger, that have been noted by three other groups, and our parameter estimation indicates that three of these may contain a sub-solar mass component. We study the data quality associated with these triggers. Finally, we discuss the challenges for parameter estimation on compact binary coalescence events containing a sub-solar mass component: long signal duration, possible small chirp masses, possible small mass ratios, and data quality issues over potentially hundreds of seconds of data.

Flexible and accurate noise characterization is crucial for the precise estimation of gravitational-wave parameters. We introduce a Bayesian method for estimating the power spectral density (PSD) of long, stationary time series, explicitly tailored for Laser Interferometer Space Antenna (LISA) data analysis. Our approach models the PSD as the geometric mean of a parametric and a nonparametric component, combining the knowledge from parametric models with the flexibility to capture deviations from theoretical expectations. The nonparametric component is expressed by a mixture of penalized B splines. Adaptive, data-driven knot placement, performed once at initialization, removes the need for a overfitting. Validation on simulated autoregressive (AR) data of order 4 [AR(4)] demonstrates estimator consistency and shows that well-matched parametric components reduce the integrated absolute error compared to an uninformative baseline, requiring fewer spline knots to achieve comparable accuracy. Applied to one year of simulated LISA X-channel (univariate) noise, our method achieves relative integrated absolute errors of O(10-2), making it suitable for iterative analysis pipelines and multiyear mission data sets.

The LIGO-Virgo-KAGRA collaboration recently reported an exceptional gravitational-wave event, GW231123. This gravitational-wave signal was assumed to be generated from the merger of a binary black hole system, with source frame masses of 137 & thorn;22-17M0 and 103 & thorn;20-52M0 (90% credible intervals). As seen by the two LIGO detectors, the signal has only 5 cycles, between 30 and 80 Hz, over 10 ms. It is of critical importance to confirm the origin of this signal. Here we present the results of a Bayesian model comparison to test whether the gravitational-wave signal was actually generated by a binary black hole merger, or emitted from cusps or kinks on a cosmic string. We find significant evidence for a binary black hole merger versus a cosmic of the

Primordial black holes (PBHs) formed in the early Universe remain a viable dark matter candidate. Since dark matter is expected to be concentrated toward galactic centres, a significant population of PBHs would reside near a supermassive black hole (SMBH), providing a promising environment for gravitational wave emission. We consider PBHs with masses from 10^-15 to 10 solar masses on hyperbolic trajectories past an SMBH, that produce bursts of gravitational radiation during periastron passage. Sagittarius A* is an obvious location for these events, but we show that both M31* and M87* would generate similar signals, albeit at lower frequencies. We assess the detectability of both individual bursts and "popcorn" backgrounds, relative to LISA and the proposed μAres mission. Comparing these results to plausible central halo densities, we find that unbound PBHs are unlikely to yield detectable signals.

This work presents the first demonstration of nonlinear noise regression in the Virgo detector using deep learning techniques. We use DeepClean, a convolutional autoencoder previously shown to be effective in denoising LIGO data, as our tool for modeling and subtracting environmental and technical noise in Virgo. The method uses auxiliary witness channels to learn correlated noise features and remove them from the strain data. For this study, we apply DeepClean to Virgo O3b data, using 225 witness channels selected across 13 targeted frequency bands. Our analysis confirms the presence of nonlinear couplings in the subtracted noise, highlighting the importance of DeepClean-like tools in capturing such effects. We observe up to a 1.3 Mpc improvement in the binary neutron star inspiral range (similar to 2.5% gain), and an average increase of 1.7% in the recovered signal-to-noise ratio for injected binary black hole signals. Parameter estimation studies further confirm that DeepClean does not introduce bias in the recovery of source parameters. These results demonstrate the robustness of DeepClean on Virgo data and support its adoption in real-time noise subtraction frameworks for future observing runs.

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.

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.

Gravitational wave detectors like the Einstein Telescope and LISA generate long multivariate time series, which pose significant challenges in spectral density estimation due to a number of overlapping signals as well as the presence of correlated noise. Addressing both issues is crucial for accurately interpreting the signals detected by these instruments. This paper presents an application of a variational inference spectral density estimation method specifically tailored for dealing with correlated noise in the data. It is flexible in that it does not rely on any specific parametric form for the multivariate spectral density. The method employs a blocked Whittle likelihood approximation for stationary time series and utilizes the Cholesky decomposition of the inverse spectral density matrix to ensure a positive definite estimator. A discounted regularized horseshoe prior is applied to the spline coefficients of each Cholesky factor, and the posterior distribution is computed using a stochastic gradient variational Bayes approach. This method is particularly effective in addressing correlated noise, a significant challenge in the analysis of multivariate data from co-located detectors. The method is demonstrated by analyzing 2000 seconds of simulated Einstein Telescope noise, which shows its ability to produce accurate spectral density estimates and quantify coherence between time series components. This makes it a powerful tool for analyzing correlated noise in gravitational wave data.

Correlated noise sources, particularly magnetic noise, form a risk to future gravitational-wave searches aimed at detecting the gravitational-wave background. To investigate potential noise contamination, we project noise levels based on accurate measurements of how strongly the noise couples to the detector. To make these estimates, we inject, for the first time, broadband, coherent magnetic noise between two gravitational-wave detectors, LIGO Hanford and LIGO Livingston, separated by several thousands of kilometers. We describe the noise injection and its impact on the analysis pipelines, then investigate the accuracy of noise projection techniques used over the past decade. Finally, we present a proof-of-concept demonstration of noise subtraction using Wiener filtering, while also highlighting potential risks associated with this method. This unique dataset with correlated noise caused by magnetic field fluctuations in two gravitational-wave detectors, as well as in an array of witness sensors, provides an excellent testing ground for additional future studies. Ultimately, this study demonstrates that Wiener filtering is effective and can be applied in the eventual detection of the gravitational-wave background by the LIGO-Virgo-KAGRA Collaboration.

Context. Recent studies have revealed that the contribution of extragalactic double white dwarfs (DWDs) to the astrophysical gravitational-wave background could be detectable in the millihertz regime by the LISA space mission. Conversely, the presence of this background could hamper the detection of cosmological backgrounds, which are among the key targets of gravitational-wave astronomy. Aims. We aim to confirm the amplitude and spectrum of the extragalactic DWD background and estimate its detectability with LISA under different assumptions. We also aim to understand the main uncertainties in the amplitude and frequency spectrum and estimate whether the signal could be anisotropic. Methods. We used the population synthesis code COSMIC with several assumptions about binary evolution and initial conditions. We also incorporated a specific treatment to account for the episodes of mass transfer and tidal torques after the formation of the DWDs. Results. Our study is in global agreement with previous studies, although we find a lower contribution at high frequencies, due to a different treatment of mass transfer in stellar binaries. We find that the uncertainties in the amplitude are dominated by the star formation model, and to a lesser degree by the binary evolution model. The inclusion of tidal effects and mass transfer episodes in DWDs can change the amplitude of the estimated background up to a factor of 3 at the highest frequencies. For all the models we consider, we find that this background would be easily detectable, with signal-to-noise-ratio values from 100 to more than 1000 by LISA after 4 years of observations. Under the hypothesis of an homogeneous Universe beyond 200 Mpc, anisotropies associated with the astrophysical population of DWDs will likely not be detectable. We provide phenomenogical fits of the background produced by extragalactic DWDs under different assumptions to be used by the community. Conclusions. We demonstrate significant variability in gravitational-wave background predictions, emphasizing uncertainties due to different astrophysical assumptions. We highlight the importance of determining the position of the knee in the gravitational-wave background spectrum, as it provides insights into mass transfer models. The prediction of this background is of critical importance for LISA in the context of observing other backgrounds.

We report correlations in underground seismic measurements with horizontal separations of several hundreds of meters to a few kilometers in the frequency range 0.01Hz to 40Hz. These seismic correlations could threaten science goals of planned interferometric gravitational-wave detectors such as the Einstein Telescope as well as atom interferometers such as MIGA and ELGAR. We use seismic measurements from four different sites, i.e. the former Homestake mine (USA) as well as two candidate sites for the Einstein Telescope, Sos Enattos (IT) and Euregio Maas-Rhein (NL-BE-DE) and the site housing the MIGA detector, LSBB (FR). At all sites, we observe significant coherence for at least 50% of the time in the majority of the frequency region of interest. Based on the observed correlations in the seismic fields, we predict levels of correlated Newtonian noise from body waves. We project the effect of correlated Newtonian noise from body waves on the capabilities of the triangular design of the Einstein Telescope's to observe an isotropic gravitational-wave background (GWB) and find that, even in case of the most quiet site, its sensitivity will be affected up to $\sim$20Hz. The resolvable amplitude of a GWB signal with a negatively sloped power-law behaviour would be reduced by several orders of magnitude. However, the resolvability of a power-law signal with a slope of e.g. $\alpha=0$ ($\alpha=2/3$) would be more moderately affected by a factor $\sim$ 6-9 ($\sim$3-4) in case of a low noise environment. Furthermore, we bolster confidence in our results by showing that transient noise features have a limited impact on the presented results.

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.

Next-generation gravitational-wave detectors like the Einstein Telescope and Cosmic Explorer, currently in their preparatory phase, have the potential to significantly improve our understanding of astrophysics, cosmology and fundamental physics. We examine how the arm orientations of the proposed detectors influence the sensitivity of a combined Einstein Telescope - Cosmic Explorer network with respect to the sensitivity to the stochastic gravitational-wave background and compact binary coalescences, where measuring both gravitational-wave polarizations is favorable. We present a method to optimize the arm orientations in the network for these two targets, and also demonstrate how to achieve a balanced configuration for both stochastic background and compact binary coalescence searches. For five specific network configurations, we explicitly compare the sensitivity to the stochastic background and binary neutron star mergers. For the latter, we conduct Bayesian parameter estimation on the extrinsic parameters of a reference binary neutron star system to assess sky localization and distance estimation capabilities. These are illustrated through efficiency curves showing the fraction of events meeting sky localization and distance uncertainty criteria as a function of redshift. Our findings suggest that globally coordinating efforts towards the next-generation gravitational-wave detector network is advantageous.

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.

Aims. We present a campaign designed to train the Global Rapid Advanced Network Devoted to the Multi-messenger Addicts (GRANDMA) network and its infrastructure to follow up on transient alerts and detect their early afterglows. In preparation for O4 II campaign, we focused on gamma-ray burst (GRB) alerts since they are expected to serve as the electromagnetic counterpart of gravitational-wave events. Our goal was to improve our response to the alerts and to start prompt observations as soon as possible, so that we may better prepare the GRANDMA network for the fourth observational run of LIGO-Virgo-Kagra (launched at the end of May 2023) and future missions such as SM. Methods. We set up a dedicated infrastructure and organized a rota of follow-up advocates (FAs) to guarantee round-the-clock assistance to our telescope teams, with an aim to receive, manage, and send out observational plans to our partner telescopes. To ensure a large number of observations, we focused on Swift GRBs whose localization errors were generally smaller than the GRANDMA telescopes’ field of view. This allowed us to bypass the transient identification process and focus on the reaction time and efficiency of the network. Results. During the ‘Ready for O4 II’ phase, 11 Swift /INTEGRAL GRB triggers were selected. Of these, nine fields had been observed and three afterglows had been detected (GRB 220403B, GRB 220427A, GRB 220514A) with 17 GRANDMA telescopes and 17 amateur astronomers from the citizen science project Kilonova-Catcher. Here, we highlight the GRB 220427A analysis, where our long-term follow-up of the host galaxy allowed us to obtain a photometric redshift of z = 0.82 ± 0.09 and its lightcurve evolution, as well as to fit the decay slope of the afterglows and study the properties of the host galaxy. Conclusions. During this eight-week-long GRB follow-up campaign, we successfully fulfilled our goal of training telescope teams for O4 and improving the associated technical toolkits. For seven of the GRB alerts, our network was able to start the first observations less than one hour after the GRB trigger time. We also characterized the network efficiency to observe GRB afterglow given the resulting time delay and limiting magnitude, and to its light curve evolution based on the observation of GRB 220427A.

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.

Magnetars are slowly rotating neutron stars that possess the strongest magnetic fields ($10^{14}-10^{15} \mathrm{G}$) known in the cosmos. They display a range of transient high-energy electromagnetic activity. The brightest and most energetic of these events are the gamma-ray bursts (GRBs) known as magnetar giant flares (MGFs), with isotropic energy $E\approx10^{44}-10^{46} \mathrm{erg}$. There are only seven detections identified as MGFs to date: three unambiguous events occurred in our Galaxy and the Magellanic Clouds, and the other four MGF candidates are associated with nearby star-forming galaxies. As all seven identified MGFs are bright at Earth, additional weaker events remain unidentified in archival data. We conducted a search of the Fermi Gamma-ray Burst Monitor (GBM) database for candidate extragalactic MGFs and, when possible, collected localization data from the Interplanetary Network (IPN) satellites. Our search yielded one convincing event, GRB 180128A. IPN localizes this burst with NGC 253, commonly known as the Sculptor Galaxy. This event is the second MGF in modern astronomy to be associated with this galaxy and the first time two bursts are associated with a single galaxy outside our own. Here, we detail the archival search criteria that uncovered this event and its spectral and temporal properties, which are consistent with expectations for a MGF. We also discuss the theoretical implications and finer burst structures resolved from various binning methods. Our analysis provides observational evidence for an eighth identified MGF.