The astrophysical foreground from compact binary coalescence signals is expected to be a dominant part of total gravitational-wave (GW) energy density in the frequency band of the third-generation detectors. The detection of any other subdominant stochastic GW background (GWB), especially a primordial GWB, will be disturbed by the astrophysical foreground. Although previous studies have proposed several cleaning methods, the foreground from subthreshold binary neutron stars (BNSs) has been a major obstacle to remove. In this paper, we propose the novel idea to better estimate the unresolved foreground, by using the information about its anisotropies. We simulate the BNS population and compute its angular power spectrum and shot noise. We find that the shot noise from BNSs is too faint to observe after subtracting loud signals due to the limited angular resolution of the third-generation detectors. This justifies the approximation regarding the unresolved foreground as an isotropic component. We also discuss the angular resolution necessary to make our method valid for the foreground subtraction.

Gravitational waves from compact binary coalescences provide unique opportunities to test general relativity (GR) in the strong-field regime. In particular, the merger phase, during which two compact objects finally coalesce, corresponds to the regime of the strongest gravitational fields accessible by direct observation and thus serves as a probe of the nonlinear nature of gravity. In this work, we test GR in the merger phase by analyzing GW150914 using a modified waveform proposed in Watarai et al. [Phys. Rev. D 109, 084058 (2024).], which parametrizes possible deviations from GR during this stage. Within this framework, the inferred deviation parameters can be translated into model-independent constraints on physically meaningful quantities. For GW150914, we find that the additional energy radiated in the merger phase is constrained to be 0.26 & thorn;0.75-0.62 % of the total energy emitted over the entire coalescence predicted by GR, and the peak time shift due to beyond-GR effects is 2.17 & thorn;9.56-9.90 ms, both within the 90% credible interval. These two constraints serve as observational benchmarks for deviations in the nonlinear gravity regime, offering guidance for theoretical investigations of beyond-GR models.

We present a novel search for dark photon dark matter (DM) using terrestrial magnetic field measurements at frequencies below 100 Hz. Coherently oscillating dark photon DM can induce a monochromatic magnetic field via kinetic mixing with ordinary photons. Notably, for dark photon masses mA around 3 & times; 10-14 eV, the signal can be resonantly amplified within a cavity formed by the Earth's surface and the ionosphere. We compute the expected signal incorporating the effect of atmospheric conductivity, and derive new upper limits on the kinetic mixing parameter epsilon from long-term geomagnetic data. These limits improve upon previous ground-based constraints in the mass range of 1 & times; 10-15 eV $ mA $ 2 & times; 10-13 eV.

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.

In tests of gravity using gravitational waves (GWs), GW events analyzed are often selected based on specific criteria, particularly the signal-to-noise ratio. However, such event selection can introduce bias into parameter estimation unless the selection effect is appropriately taken into account in the analysis. In this paper, we investigate how event selection with certain prior information affects parameter inference within the scalar-tensor polarization framework, focusing on the measurement of the scalar mode amplitude parameters. We find that for the Tensor+Scalar(dipole) model, the amplitude of the scalar dipole radiation is overestimated when its true value is nonzero while there is no false deviation in the absence of the scalar mode. The same bias is expected to occur also for the Tensor+Scalar(quadrupole) model. However, the error typically exceeds the bias as the scalar quadrupole mode is difficult to be distinguished from the tensor mode.

We have proposed using laser interferometric gravitational wave detectors to search for ultralight vector and axion dark matter. Vector dark matter can be probed through oscillating forces on suspended mirrors, while axion dark matter can be detected via oscillating polarization rotation of laser beams. This paper reviews these searches with the KAGRA detector in Japan, including the first vector dark matter search with KAGRA's 2020 data and installation of polarization optics for axion dark matter search during the upcoming 2025 observing run.

The displacement-noise-free interferometer (DFI) is designed to eliminate all displacement-induced noise while retaining sensitivity to gravitational wave (GW) signals. Ground-based DFIs suffer from physical arm-length limitations, resulting in poor sensitivity at frequencies below 1 kHz. To address this, previous research introduced a neutron-based DFI, which replaces laser light with neutrons and achieves exceptional sensitivity down to a few hertz. In this study, we conducted a proof-of-principle experiment using a pulsed neutron source at the Japan Proton Accelerator Research Complex (J- PARC). Despite practical constraints that led to deviations from the ideal experimental design, we optimized the setup and developed a novel analysis method that successfully cancels displacement noise while preserving simulated GW signals. This work presents the first successful demonstration of a neutron DFI and a neutron interferometer for GW detection.

The detection of low-frequency gravitational waves on Earth requires the reduction of displacement noise, which dominates the low-frequency band. One method to cancel test mass displacement noise is a neutron displacement-noise-free interferometer (DFI). This paper proposes a new neutron DFI configuration, a Sagnac-type neutron DFI, which uses a Sagnac interferometer in place of the Mach-Zehnder interferometer. We demonstrate that a sensitivity of the Sagnac-type neutron DFI is higher than that of a conventional neutron DFI with the same interferometer scale. This configuration is particularly significant for neutron DFIs with limited space for construction and limited flux from available neutron sources.

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 merger phase of compact binary coalescences is the strongest gravity regime that can be observed. To test the validity of general relativity (GR) in strong gravitational fields, we propose a gravitational waveform parameterized for deviations from GR in the dynamical and nonlinear regime of gravity. Our fundamental idea is that perturbative modifications to a GR waveform can capture possible deviations in the merger phase that are difficult to model in a specific theory of gravity. One of notable points is that our waveform is physically consistent in the sense that the additional radiative losses of energy and angular momentum associated with beyond-GR modifications are included. Our prescription to ensure physical consistency in the whole coalescence process is expected to be applicable to any deviation from the standard model of compact binary coalescence, such as the extended models of gravity or the environmental effects of compact objects, as long as perturbative modifications are considered. Based on the Fisher analysis and the compatibility with Einstein-dilaton Gauss-Bonnet waveforms, we show that our parameterization is a physically-consistent minimal one that captures the deviations in the nonlinear regime.

In nonminimally coupled theories where a scalar field is coupled to the Ricci scalar, neutron stars (NSs) can have scalar charges through an interaction with matter mediated by gravity. On the other hand, the same theories do not give rise to hairy black hole (BH) solutions. The observations of gravitational waves (GWs) emitted from an inspiralling NS-BH binary system allows a possibility of constraining the NS scalar charge. Moreover, the nonminimally coupled scalar-tensor theories generate a breathing scalar mode besides two tensor polarizations. Using the GW200115 data of the coalescence of a BH-NS binary, we place observational constraints on the NS scalar charge as well as the nonminimal coupling strength for a subclass of massless Horndeski theories with a luminal GW propagation. Unlike past related works, we exploit a waveform for a mixture of tensor and scalar polarizations. Taking the breathing mode into account, the scalar charge is more tightly constrained in comparison to the analysis of the tensor GWs alone. In nonminimally coupled theories including Brans-Dicke gravity and spontaneous scalarization scenarios with/without a kinetic screening, we put new bounds on model parameters of each theory.

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.

Magnetars are neutron stars with exceptionally strong dipole magnetic fields which are observed to display a range of x-ray flaring behavior, but the flaring mechanism is not well understood. The third observing run of Advanced LIGO and Virgo extended from April 1, 2019 to March 27, 2020, and contained x-ray flares from known magnetar SGR 1935+2154, as well as the newly-discovered magnetar, Swift J1818-1607. We search for gravitational waves coincident with these magnetar flares with minimally modeled, coherent searches which specifically target both short-duration gravitational waves produced by excited f-modes in the magnetar's core, as well as long-duration gravitational waves motivated by the Quasi-Periodic Oscillations observed in the tails of giant flares. In this paper, we report on the methods and sensitivity estimates of these searches, and the astrophysical implications.

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 time scale of DM.

Gravitational lensing by massive objects along the line of sight to the source causes distortions to gravitational wave (GW) signals; such distortions may reveal information about fundamental physics, cosmology, and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO-Virgo network. We search for repeated signals from strong lensing by (1) performing targeted searches for subthreshold signals, (2) calculating the degree of overlap among the intrinsic parameters and sky location of pairs of signals, (3) comparing the similarities of the spectrograms among pairs of signals, and (4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by (1) frequency-independent phase shifts in strongly lensed images, and (2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the nondetection of GW lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects.

KAGRA is a newly built gravitational wave observatory, a laser interferometer with a 3 km arm length, located in Kamioka, Gifu prefecture, Japan. In this article, we describe the KAGRA data management system, i.e., recording of data, transfer from the KAGRA experiment site to computing resources, as well as data distribution to tier sites, including international sites in Taiwan and Korea. The amount of KAGRA data exceeded 1.0 PiB and increased by about 1.5 TB per day during operation in 2020. Our system has succeeded in data management, and has achieved performance that can withstand observations after 2023, that is, a transfer rate of 20 MB s-1 or more and file storage of sufficient capacity for petabyte class. We also discuss the sharing of data between the global gravitational-wave detector network with other experiments, namely LIGO and Virgo. The latency, which consists of calculation of calibrated strain data and transfer time within the global network, is very important from the view of multi-messenger astronomy using gravitational waves. Real-time calbrated data delivered from the KAGRA detector site and other detectors to our computing system arrive with about 4-15 seconds of latency. These latencies are sufficiently short compared to the time taken for gravitational wave event search computations. We also established a high-latency exchange of offline calibrated data that was aggregated with a better accuracy compared with real-time data.

KAGRA, the kilometer-scale underground gravitational-wave detector, is located at Kamioka, Japan. In April 2020, an astrophysics observation was performed at the KAGRA detector in combination with the GEO 600 detector; this observation operation is called O3GK. The optical configuration in O3GK is based on a power recycled Fabry-Pérot Michelson interferometer; all the mirrors were set at room temperature. The duty factor of the operation was approximately 53%, and the strain sensitivity was $3\times10^{-22}~/\sqrt{\rm{Hz}}$ at 250 Hz. In addition, the binary-neutron-star (BNS) inspiral range was approximately 0.6 Mpc. The contributions of various noise sources to the sensitivity of O3GK were investigated to understand how the observation range could be improved; this study is called a "noise budget". According to our noise budget, the measured sensitivity could be approximated by adding up the effect of each noise. The sensitivity was dominated by noise from the sensors used for local controls of the vibration isolation systems, acoustic noise, shot noise, and laser frequency noise. Further, other noise sources that did not limit the sensitivity were investigated. This paper provides a detailed account of the KAGRA detector in O3GK including interferometer configuration, status, and noise budget. In addition, strategies for future sensitivity improvements such as hardware upgrades, are discussed.

Axions are considered as a candidate of dark matter. The axions form coherent clouds which delay and amplify gravitational waves (GWs) at a resonant frequency. That is, the coherent axion clouds produce secondary GWs following a primary wave. All GWs from compact binary mergers detected so far propagate in the Milky Way halo composed of dark matter so that such GWs are followed by the secondary GWs. Since the properties of the secondary GWs depend on the axion mass and coupling with the parity violating sector of gravity, we can search for a characteristic signal produced by axions. In our previous study, we have developed a search method optimized for the axion signals, and obtained about ten times stronger constraint on the coupling than the previous best one for the axion mass range, $[1.7 \times 10^{-13},\,8.5 \times 10^{-12}]\,\mathrm{eV}$. However, in our previous search, we assume that dark matter is homogeneously distributed in the Milky Way halo, which is unrealistic. In this paper, we extend the dark matter profile to a more realistic one called the core NFW profile, and find that the constraint in our previous study is robust to the dark matter profile.