Gravitational-wave astronomy has developed enormously over the last decade with the first detections across different frequency bands, but has yet to access $0.1-10$ $\mathrm{Hz}$ gravitational waves. Gravitational waves in this band are emitted by some of the most enigmatic sources, including intermediate-mass binary black hole mergers, early inspiralling compact binaries, and possibly cosmic inflation. To tap this exciting band, we propose the construction of a detector based on pulsar timing principles, the Artificial Precision Timing Array (APTA). We envision APTA as a solar system array of artificial "pulsars"$-$precision-clock-carrying satellites that emit pulsing electromagnetic signals towards Earth or other centrum. In this fundamental study, we estimate the clock precision needed for APTA to successfully detect gravitational waves. Our results suggest that a clock relative uncertainty of $10^{-17}$, which is currently attainable, would be sufficient for APTA to surpass LISA's sensitivity in the decihertz band and observe $10^3-10^4$ $\mathrm{M}_\odot$ black hole mergers. Future atomic clock technology realistically expected in the next decade would enable the detection of an increasingly diverse set of astrophysical sources, including stellar-mass compact binaries that merge in the LIGO-Virgo-KAGRA band, extreme-mass-ratio inspirals, and Type Ia supernovae. This work opens up a new area of research into designing and constructing artificial gravitational-wave detectors relying on the successful principles of pulsar timing.

We analyzed the 7.92 × 10 ^11 cosmic-ray-induced muon events collected by the IceCube Neutrino Observatory from 2011 May 13, when the fully constructed experiment started to take data, to 2023 May 12. This data set provides an up-to-date cosmic-ray arrival direction distribution in the Southern Hemisphere with unprecedented statistical accuracy covering more than a full period length of a solar cycle. Improvements in Monte Carlo event simulation and better handling of year-to-year differences in data processing significantly reduce systematic uncertainties below the level of statistical fluctuations compared to the previously published results. We confirm the observation of a change in the angular structure of the cosmic-ray anisotropy between 10 TeV and 1 PeV, more specifically in the 100–300 TeV energy range. For the first time, we analyzed the angular power spectrum at different energies. The observed variations of the power spectra with energy suggest relatively reduced large-scale features at high energy compared to those of medium and small scales. The large volume of data enhances the statistical significance at higher energies, up to the PeV scale, and smaller angular scales, down to approximately 6° compared to previous findings.

Mergers of compact objects, binary black holes and mergers including at least one neutron star, are a predicted source of high-energy neutrinos. These astrophysical events are now routinely detected through observation of their gravitational wave signature and, at least in one instance, their electromagnetic counterparts were also detected. Particles accelerated during the coalescence of compact objects may also interact to produce high-energy neutrinos, which have yet to be detected, but observations are ongoing. The LIGO-Virgo-KAGRA Collaboration publicly releases information on candidate gravitational wave events from compact binary coalescences in low latency during the current observing run (O4). To aid the electromagnetic follow-up, using data from the IceCube Neutrino Observatory, we search, in real time, for neutrinos spatially and temporally coincident with these gravitational wave candidate events using a time window of 1000 seconds centered on the merger time. We use two methods, both of which have been previously used to search for neutrino emission from gravitational-wave transients: an unbinned maximum likelihood analysis applied to significant alerts and a Bayesian analysis with astrophysical priors, applied to both significant and low-significance alerts. In addition, we search for long-duration neutrino emission up to 14 days after the merging of binaries containing a neutron star. We report analysis results determined in real time for these searches, and set upper limits on both flux and isotropic-equivalent energy emitted in neutrinos.

While the sources of the diffuse astrophysical neutrino flux detected by the IceCube Neutrino Observatory are still largely unknown, one of the promising methods to improve our understanding of them is investigating the potential temporal and spatial correlations between neutrino alerts and the electromagnetic radiation from blazars. We report on the multiwavelength target-of-opportunity observations of the blazar B3 2247+381, taken in response to an IceCube multiplet alert for a cluster of muon neutrino events compatible with the source location between 2022 May 20 and 2022 November 10. B3 2247+381 was not detected with VERITAS during this time period. The source was found to be in a low-flux state in the optical, ultraviolet, and gamma-ray bands for the time interval corresponding to the neutrino event, but was detected in the hard X-ray band with NuSTAR during this period. We find the multiwavelength spectral energy distribution is described well using a simple one-zone leptonic synchrotron self-Compton radiation model. Moreover, assuming the neutrinos originate from hadronic processes within the jet, the neutrino flux would be accompanied by a photon flux from the cascade emission, and the integrated photon flux required in such a case would significantly exceed the total multiwavelength fluxes and the VERITAS upper limits presented here. The lack of flaring activity observed with VERITAS, combined with the low multiwavelength flux levels, as well as the significance of the neutrino excess being at a 3 σ level (uncorrected for trials), makes B3 2247+381 an unlikely source of the IceCube multiplet. We conclude that the neutrino excess is likely a background fluctuation.

We report a search for high-energy astrophysical neutrino multiplets, detections of multiple neutrino clusters in the same direction within 30 days, based on an analysis of 11.4 yr of IceCube data. A new search method optimized for transient neutrino emission with a monthly timescale is employed, providing a higher sensitivity to neutrino fluxes. This result is sensitive to neutrino transient emission, reaching per-flavor flux of approximately 1 0 − 10 erg cm − 2 s − 1 from the Northern Sky in the energy range E ≳ 50 TeV. The number of doublets and triplets identified in this search is compatible with the atmospheric background hypothesis, which leads us to set limits on the nature of neutrino transient sources with emission timescales of one month.

The recent IceCube detection of TeV neutrino emission from the nearby active galaxy NGC 1068 suggests that active galactic nuclei (AGNs) could make a sizable contribution to the diffuse flux of astrophysical neutrinos. The absence of TeV γ -rays from NGC 1068 indicates neutrino production in the vicinity of the supermassive black hole, where the high radiation density leads to γ -ray attenuation. Therefore, any potential neutrino emission from similar sources is not expected to correlate with high-energy γ -rays. Disk-corona models predict neutrino emission from Seyfert galaxies to correlate with keV X-rays because they are tracers of coronal activity. Using through-going track events from the Northern Sky recorded by IceCube between 2011 and 2021, we report results from a search for individual and aggregated neutrino signals from 27 additional Seyfert galaxies that are contained in the Swift's Burst Alert Telescope AGN Spectroscopic Survey. Besides the generic single power law, we evaluate the spectra predicted by the disk-corona model assuming stochastic acceleration parameters that match the measured flux from NGC 1068. Assuming all sources to be intrinsically similar to NGC 1068, our findings constrain the collective neutrino emission from X-ray bright Seyfert galaxies in the northern sky, but, at the same time, show excesses of neutrinos that could be associated with the objects NGC 4151 and CGCG 420-015. These excesses result in a 2.7 σ significance with respect to background expectations.

Active galactic nuclei (AGN) are promising candidate sources of high-energy astrophysical neutrinos, since they provide environments rich in matter and photon targets where cosmic-ray interactions may lead to the production of gamma rays and neutrinos. We searched for high-energy neutrino emission from AGN using the Swift-BAT Spectroscopic Survey catalog of hard X-ray sources and 12 yr of IceCube muon track data. First, upon performing a stacked search, no significant emission was found. Second, we searched for neutrinos from a list of 43 candidate sources and found an excess from the direction of two sources, the Seyfert galaxies NGC 1068 and NGC 4151. We observed NGC 1068 at flux ϕ ν μ + ν ¯ μ = 4.0 2 − 1.52 + 1.58 × 1 0 − 11 TeV −1 cm −2 s −1 normalized at 1 TeV, with a power-law spectral index γ = 3.10 − 0.22 + 0.26 , consistent with previous IceCube results. The observation of a neutrino excess from the direction of NGC 4151 is at a posttrial significance of 2.9 σ . If interpreted as an astrophysical signal, the excess observed from NGC 4151 corresponds to a flux ϕ ν μ + ν ¯ μ = 1.5 1 − 0.81 + 0.99 × 1 0 − 11 TeV −1 cm −2 s −1 normalized at 1 TeV and γ = 2.83 − 0.28 + 0.35 .

The discovery of joint sources of high-energy neutrinos and gravitational waves has been a primary target for the LIGO, Virgo, KAGRA, and IceCube observatories. The joint detection of high-energy neutrinos and gravitational waves would provide insight into cosmic processes, such as progenitor dynamics and outflows. The joint detection of multiple cosmic messengers can also elevate the significance of the observation when some or all of the constituent messengers are sub-threshold, not significant enough to declare their detection individually. Leveraging data from the LIGO, Virgo, and IceCube observatories, we conducted an archival investigation of sub-threshold multimessenger events. Complementing previous analyses, we used minimal assumptions to search for common sources of sub-threshold gravitational-wave and high-energy neutrino candidates during the third observing run (O3) of the Advanced LIGO and Advanced Virgo detectors. Our search did not identify significant joint sources. We therefore derive constraints on the rate density of joint sources for each compact binary merger population as a function of the energy emitted in neutrinos. Only a fraction of the gravitational-wave sources emit neutrinos, if the neutrino emission has high bolometric energy (>10^52 to 10^54 erg).

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.

Coincident multimessenger observations of cosmic sources can offer numerous benefits, especially when used in the context of synergistic astrophysics. One significant advantage is enhancing the detection significance of separate detectors by correlating their data and assuming joint emission. We have formulated an approach for updating the Bayesian posterior probability of an astrophysical origin, namely pastro, relying on multimessenger coincidences assuming an emission model. The description is applicable to any combination of messengers. We demonstrated the formalism for the gravitational waves and highenergy neutrinos case. Applying our method to the public data of candidate coincident high-energy neutrinos with subthreshold gravitational-wave triggers, we found that in the case of highly energetic neutrino coincidences, pastro can increase from approximately similar to 0.1 to similar to 0.9. The amount of improvement depends on the assumed joint emission model. If models are trusted, the marked improvement makes subthreshold detections much more confident. Moreover, the model dependency can also be used to test the consistency of different models. This work is a crucial step toward the goal of uniting all detectors on equal footing into a statistically integrated, Earth-sized observatory for comprehensive multimessenger astrophysics.

We present the results of a search for 10–1000 GeV neutrinos from 2268 gamma-ray bursts (GRBs) over 8 yr of IceCube-DeepCore data. This work probes burst physics below the photosphere where electromagnetic radiation cannot escape. Neutrinos of tens of giga electronvolts are predicted in sub-photospheric collision of free-streaming neutrons with bulk-jet protons. In a first analysis, we searched for the most significant neutrino-GRB coincidence using six overlapping time windows centered on the prompt phase of each GRB. In a second analysis, we conducted a search for a group of GRBs, each individually too weak to be detectable, but potentially significant when combined. No evidence of neutrino emission is found for either analysis. The most significant neutrino coincidence is for Fermi-GBM GRB bn 140807500, with a p -value of 0.097 corrected for all trials. The binomial test used to search for a group of GRBs had a p -value of 0.65 after all trial corrections. The binomial test found a group consisting only of GRB bn 140807500 and no additional GRBs. The neutrino limits of this work complement those obtained by IceCube at tera electronvolt to peta electronvolt energies. We compare our findings for the large set of GRBs as well as GRB 221009A to the sub-photospheric neutron-proton collision model and find that GRB 221009A provides the most constraining limit on baryon loading. For a jet Lorentz factor of 300 (800), the baryon loading on GRB 221009A is lower than 3.85 (2.13) at a 90% confidence level.

State-of-the-art quantum computers have very limited applicability for accurate calculations. Here, we report the first experimental demonstration of qubit-based matched filtering for a detection of the gravitational-wave signal from a binary black hole merger. With our implementation on noisy superconducting qubits, we obtained a similar signal-to-noise ratio for the binary black hole merger as achievable with classical computation, providing evidence for the utility of qubits for practically relevant tasks. The algorithm we invented for this application is a Monte Carlo algorithm which uses quantum and classical computation together. It provides a quasi-quadratic speed-up for time-domain convolution, similar to achievable with fast Fourier transform.

Neutrino flares in the sky are searched for in data collected by IceCube between 2011 and 2021 May. This data set contains cascade-like events originating from charged-current electron neutrino and tau neutrino interactions and all-flavor neutral-current interactions. IceCube’s previous all-sky searches for neutrino flares used data sets consisting of track-like events originating from charged-current muon neutrino interactions. The cascade data set is statistically independent of the track data sets, and while inferior in angular resolution, the low-background nature makes it competitive and complementary to previous searches. No statistically significant flare of neutrino emission was observed in an all-sky scan. Upper limits are calculated on neutrino flares of varying duration from 1 hr to 100 days. Furthermore, constraints on the contribution of these flares to the diffuse astrophysical neutrino flux are presented, showing that multiple unresolved transient sources may contribute to the diffuse astrophysical neutrino flux.

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 origin of high-energy galactic cosmic rays is yet to be understood, but some galactic cosmic-ray accelerators can accelerate cosmic rays up to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding material or radiation, resulting in the production of gamma-rays and neutrinos. To optimize for the detection of such associated production of gamma-rays and neutrinos for a given source morphology and spectrum, a multimessenger analysis that combines gamma-rays and neutrinos is required. In this study, we use the Multi-Mission Maximum Likelihood framework with IceCube Maximum Likelihood Analysis software and HAWC Accelerated Likelihood to search for a correlation between 22 known gamma-ray sources from the third HAWC gamma-ray catalog and 14 yr of IceCube track-like data. No significant neutrino emission from the direction of the HAWC sources was found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit from the 22 sources. From the neutrino flux limit, we conclude that, for five of the sources, the gamma-ray emission observed by HAWC cannot be produced purely from hadronic interactions. We report the limit for the fraction of gamma-rays produced by hadronic interactions for these five sources.

This Letter presents the result of a 3+1 sterile neutrino search using 10.7 yr of IceCube data. We analyze atmospheric muon neutrinos that traverse the Earth with energies ranging from 0.5 to 100 TeV, incorporating significant improvements in modeling neutrino flux and detector response compared to earlier studies. Notably, for the first time, we categorize data into starting and throughgoing events, distinguishing neutrino interactions with vertices inside or outside the instrumented volume, to improve energy resolution. The best-fit point for a 3+1 model is found to be at sin^{2}(2θ_{24})=0.16 and Δm_{41}^{2}=3.5  eV^{2}, which agrees with previous iterations of this Letter. The result is consistent with the null hypothesis of no sterile neutrinos with a p value of 3.1%.

We present the first three-parameter fit to a 3+1 sterile neutrino model using 7.634 years of data from the IceCube Neutrino Observatory on nu(mu) + (nu) over bar (mu) charged-current interactions in the energy range 500-9976 GeV. Our analysis is sensitive to the mass-squared splitting between the heaviest and lightest mass state (Delta m(41)(2)), the mixing matrix element connecting muon flavor to the fourth mass state (vertical bar U-mu 4 vertical bar(2)), and the element connecting tau flavor to the fourth mass state (vertical bar U-tau 4 vertical bar(2)). Predicted propagation effects in matter enhance the signature through a resonance as atmospheric neutrinos from the Northern Hemisphere traverse the Earth to the IceCube detector at the South Pole. The remaining sterile neutrino matrix elements are left fixed, with vertical bar U-e4 vertical bar(2) = 0 and delta(14) = 0, as they have a negligible effect, and delta(24) = pi is set to give the most conservative limits. The result is consistent with the no-sterile neutrino hypothesis with a probability of 4.3%. Profiling the likelihood of each parameter yields the 90% confidence levels: 2.4 eV(2) < Delta m(41)(2) < 9.6 eV(2), 0.0081 < vertical bar U-mu 4 vertical bar(2) < 0.10, and vertical bar U-tau 4 vertical bar(2) < 0.035, which narrows the allowed parameter-space for vertical bar U-tau 4 vertical bar(2). However, the primary result of this analysis is the first map of the 3+1 parameter space exploring the interdependence of Delta m(41)(2), vertical bar U-mu 4 vertical bar(2), and |U-tau 4 vertical bar(2).