Showing 1–50 of 526 results for author: Coughlin, M

Learned proposals in trans-dimensional inference are optimal at equilibrium, not during assembly

Authors: Argyro Sasli, Nikolaos Karnesis, Minas Karamanis, Michael L. Katz, Dimitrios Kourtesis, Michael W. Coughlin, Vuk Mandic, Nikolaos Stergioulas

Abstract: Inferring the dimension of a model - the number of components needed to explain data - jointly with the parameters is a pervasive problem, from counting sources in an image to mixture modeling, and reversible-jump Markov chain Monte Carlo solves it exactly but mixes slowly. Learned proposals are well established at fixed dimension, but whether they can accelerate the dimension-changing moves thems… ▽ More Inferring the dimension of a model - the number of components needed to explain data - jointly with the parameters is a pervasive problem, from counting sources in an image to mixture modeling, and reversible-jump Markov chain Monte Carlo solves it exactly but mixes slowly. Learned proposals are well established at fixed dimension, but whether they can accelerate the dimension-changing moves themselves has remained largely untested. We show that the answer has a structural origin: the optimal proposal for the dimension-changing birth move is a different object in different phases of the run. While the fit is being assembled it must match the current residual - a state-dependent quantity no state-independent network can represent - but at equilibrium it degenerates to the posterior's single-component marginal, which is exactly the distribution an adaptive normalizing flow learns from the sampler's own history. A learned state-independent birth proposal is therefore useless in one phase and optimal in the other. Controlled experiments confirm the attribution: applied with an exact Metropolis--Hastings correction that leaves the target invariant for any network, the learned births leave acceptance rates unchanged yet accelerate model-order mixing - in a ten-seed benchmark they meet a pre-specified stopping rule in six of ten runs, typically several times sooner, where a strong hand-tuned baseline meets it in one (one-sided p=0.03) - and an isolation experiment shows the same flow deployed within-model buys nothing. Making no domain-specific assumptions, the same sampler counts sources in a noisy image and reconstructs signals across scientific domains, including gravitational waves from ground- and space-based detectors and a scalp EEG recording. We release the method as HyperWave, an open-source package. △ Less

Submitted 4 August, 2026; originally announced August 2026.

Comments: 23 pages, 5 figures, 2 tables

LIGO A$^\sharp$: Detector Design and Science Prospects Beyond A+

Authors: L. Sun, K. Kuns, B. J. J. Slagmolen, P. Fritschel, P. Schmidt, B. T. Lantz, S. S. Y. Chua, Divyajyoti, S. W. Ballmer, M. A. Barton, A. V. Cumming, K. L. Dooley, J. C. Driggers, A. Effler, M. Evans, B. Farr, G. González, N. Lu, D. J. Ottaway, C. Palomba, O. J. Piccinni, G. Pratten, S. Raja, A. P. Subhash, P. J. Sutton , et al. (1131 additional authors not shown)

Abstract: We present the LIGO A$^\sharp$ detector concept, an upgrade for the LIGO observatories based on room-temperature interferometers beyond the fifth observing run (O5). Building on the A+ sensitivity, A$^\sharp$ targets broadband sensitivity improvements through heavier test masses, improved suspensions and seismic isolation, increased arm-cavity power, enhanced frequency-dependent squeezing, reduced… ▽ More We present the LIGO A$^\sharp$ detector concept, an upgrade for the LIGO observatories based on room-temperature interferometers beyond the fifth observing run (O5). Building on the A+ sensitivity, A$^\sharp$ targets broadband sensitivity improvements through heavier test masses, improved suspensions and seismic isolation, increased arm-cavity power, enhanced frequency-dependent squeezing, reduced coating thermal noise considering two scenarios, and improved control of mechanical motion and optical modes. We describe the principal design choices, projected noise performance, and corresponding astrophysical prospects. LIGO A$^\sharp$ substantially increases compact-binary detection rates, strengthens population inference, and improves both early-warning times and localization for binary neutron star mergers. The improved sensitivity enables more detailed studies of compact-binary coalescences, including higher-order multipoles, intermediate-mass black holes, remnant black hole ringdown, and the neutron star equation of state. It also broadens the discovery potential for new gravitational-wave sources such as continuous waves and bursts, should enable detection of the stochastic background from compact binary mergers if it remains undetected after O5, and strengthens the role of gravitational-wave detectors as probes of fundamental physics. We discuss key technical challenges and the role of A$^\sharp$ as both a major scientific upgrade for the 2030s and a technology pathfinder for next-generation gravitational-wave observatories, such as Cosmic Explorer. △ Less

Submitted 12 August, 2026; originally announced August 2026.

Comments: 78 pages, 22 figures

Report number: LIGO-P2600307

Constraints on ultralight bosons from merging binary and remnant black holes observed during the second and third parts of the fourth LIGO-Virgo-KAGRA observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1786 additional authors not shown)

Abstract: We present constraints on ultralight bosons using binary black hole mergers observed in the second and third parts of the fourth LIGO-Virgo-KAGRA observing run. Directed searches are conducted for long-transient gravitational waves from ultralight vector boson clouds around merger remnants, using a hidden-Markov-model (HMM) tracking scheme. We target the remnant black holes formed in the binary co… ▽ More We present constraints on ultralight bosons using binary black hole mergers observed in the second and third parts of the fourth LIGO-Virgo-KAGRA observing run. Directed searches are conducted for long-transient gravitational waves from ultralight vector boson clouds around merger remnants, using a hidden-Markov-model (HMM) tracking scheme. We target the remnant black holes formed in the binary coalescences that produced GW250114 and GW250207. We find no evidence for such signals from either target. Estimating our search sensitivity at a threshold corresponding to a 1% false alarm probability, we thus disfavor vector boson masses in the range of $[2.80, 3.95]\times 10^{-13}$ eV with greater than 90% confidence. In addition, we derive constraints on ultralight scalar and vector bosons from the inferred high spins of the constituent black holes in three binaries, using events GW240515, GW241113, and GW241225_08. The excluded mass ranges in this approach depend on the assumed black-hole ages. At $10^5$ years, corresponding to typical dynamically formed binaries, we exclude scalar and vector bosons in the ranges $[1.39, 6.94]\times 10^{-13}$ eV and $[0.32, 14.4]\times 10^{-13}$ eV at 90% confidence, respectively. △ Less

Submitted 11 August, 2026; originally announced August 2026.

Comments: 27 pages (12 pages author list, 15 pages main paper), 5 figures

Report number: LIGO-P2600218

Binary neutron stars in the next-generation era: Multi-messenger detection prospects and constraints on the equation of state, mass distribution, and cosmology

Authors: Hauke Koehn, Thibeau Wouters, Gilad Sadeh, Peter T. H. Pang, Mattia Bulla, Chris Van Den Broeck, Michael W. Coughlin, Tim Dietrich

Abstract: Next-generation gravitational-wave (GW) observatories will provide crucial insights into the nature of neutron star (NS) matter and the cosmological expansion history. We estimate the number of multi-messenger detections from binary neutron stars (BNS) with the Einstein Telescope (ET) and Cosmic Explorer (CE), and project the resulting constraints on the equation of state (EOS), BNS mass distribut… ▽ More Next-generation gravitational-wave (GW) observatories will provide crucial insights into the nature of neutron star (NS) matter and the cosmological expansion history. We estimate the number of multi-messenger detections from binary neutron stars (BNS) with the Einstein Telescope (ET) and Cosmic Explorer (CE), and project the resulting constraints on the equation of state (EOS), BNS mass distribution, and cosmology via joint hierarchical Bayesian inference. Assuming a local merger rate of 106.6 Gpc$^{-3}$ yr$^{-1}$ and considering two different mass functions, a narrow one centred around 1.4 $M_\odot$ and a wide one ranging between 1.1--2 $M_\odot$, we find that for ET, our mock follow-up algorithm results in at least $\sim40$ and up to $\sim100$ successfully identified electromagnetic counterparts per year, depending on the detector layout and mass distribution. In a joint network with CE, the number of multi-messenger detections can range from $\sim 200$ to $\sim500$. Additionally, several more afterglows from gamma-ray bursts or KNe could be found with dedicated late-time observations. Based on the identified multi-messenger events, we perform an injection campaign to hierarchically constrain the EOS, mass distribution, and cosmology in a fully Bayesian framework. Focussing on ET alone, we show how in an ideal scenario, GW signals, KNe, and host galaxy redshifts can constrain the canonical NS radius $R_{1.4}$ within $\sim 0.2$ km and the Hubble constant $H_0$ within $\sim 1$ km s$^{-1}$ Mpc$^{-1}$, while recovering the essential features of the mass distribution. By comparing inference results that rely solely on GW data and those that incorporate light curve information, we find that while KN light-curve posteriors have a negligible impact on the EOS constraints, they can benefit the inference of cosmological parameters. △ Less

Submitted 30 July, 2026; originally announced July 2026.

Comments: 30 pages, 14 figures, comments welcome

4th TDAMM Workshop White Paper

Authors: Daniel Kocevski, Raffaella Margutti, Michael Coughlin, Christopher Fryer, Kate Alexander, Igor Lorenzo Andreoni, Jean-Luc Atteia, Elias Aydi, Matthew Baring, Joe Bright, Eric Burns, Brad Cenko, Sanjana Curtis, Fabio DeColle, Courey Elliott, Corinne Fletcher, Anna Franckowiak, Carla Fr"ohlich, Adam Goldstein, Alec Habig, Erica Hammerstein, Anna Y. Q. Ho, D. Andrew Howell, Chin-Ping Hu, Xiaoshan Huang , et al. (24 additional authors not shown)

Abstract: Time-Domain and Multi-Messenger Astrophysics (TDAMM) is entering a new era in which the rate and diversity of transient discoveries will grow rapidly across electromagnetic, gravitational-wave, neutrino, and cosmic-ray facilities. The scientific return from these investments will increasingly depend not on discovery alone, but on the ability to identify, prioritize, and coordinate follow-up observ… ▽ More Time-Domain and Multi-Messenger Astrophysics (TDAMM) is entering a new era in which the rate and diversity of transient discoveries will grow rapidly across electromagnetic, gravitational-wave, neutrino, and cosmic-ray facilities. The scientific return from these investments will increasingly depend not on discovery alone, but on the ability to identify, prioritize, and coordinate follow-up observations across a heterogeneous and globally distributed network of observatories. This white paper summarizes the outcomes of the Fourth TDAMM Workshop and assesses the near-term discovery landscape, the infrastructure and tools that support coordinated observations, and the technical, policy, and capability gaps that may limit future progress. The workshop identified three principal challenges: insufficiently scalable and interoperable alert and coordination infrastructure, policies that impede rapid multi-facility observations and rare-event science, and the potential loss of critical high-energy, rapid-response, and spectroscopic capabilities. The white paper identifies the need for sustained support for alert distribution, brokers, standardized observatory metadata, cross-facility coordination platforms, and unified follow-up repositories; expanded joint observing opportunities and funding mechanisms for coordinated analysis; and strategic investment in future TDAMM facilities. The white paper also present a framework for community observing plans that would establish pre-coordinated responses to rare, high-impact events, supported by transparent governance, immediate public data release, and regular community revision. Science overviews and detailed observing strategies are provided for gamma-ray bursts, tidal disruption events, X-ray binaries, novae, supernovae, magnetars, compact binary mergers, and high-energy neutrino sources. △ Less

Submitted 24 July, 2026; originally announced July 2026.

Comments: 140 pages, 1 figure

GWTC-5.0: Tests of General Relativity

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1800 additional authors not shown)

Abstract: The signals from the LIGO-Virgo-KAGRA network of gravitational-wave (GW) detectors allow us to perform sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. We present the results of seven tests of GR using the observed binary signals in the fifth GW Transient Catalog (GWTC-5.0), i.e., up to and including the second part of the fourth observing run (O4b).… ▽ More The signals from the LIGO-Virgo-KAGRA network of gravitational-wave (GW) detectors allow us to perform sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. We present the results of seven tests of GR using the observed binary signals in the fifth GW Transient Catalog (GWTC-5.0), i.e., up to and including the second part of the fourth observing run (O4b). We restrict our analysis to the confident signals, henceforth called events, observed by at least two detectors that have estimated false alarm rates $\le 10^{-3} \ \rm{yr}^{-1}$. These include 72 events from O4b and five events from the first part of the fourth observing run that are now analyzed due to their increased significance from updated search results, bringing the total number of events for tests of GR in the cumulative GWTC to 168. After subtracting the best-fit waveforms, we find the residuals are consistent with detector noise for all events considered. We also find no strong evidence for additional polarizations beyond those predicted by GR. We perform tests of GW generation, improving the constraints on deviations from the GR post-Newtonian coefficients by factors of 1.2-2.6. Finally, we find overall consistency of the remnants with GR using both time- and frequency-domain methods. For GW240621_195059, postmerger data are consistent with the dominant quadrupolar ($\ell=|m|=2$) mode of a Kerr black hole and its first overtone, with spurious high-frequency content preventing a spectroscopic constraint of GR. In the frequency-domain ringdown analysis, the GR prediction lies in the tails of the combined results, possibly due to the limited catalog size. However, the combined results indicate improved consistency with GR over GWTC-4.0, owing to the contribution of GW250114 with a network matched-filter signal-to-noise ratio of 76.9. Overall, we find no evidence for physics beyond GR. △ Less

Submitted 21 July, 2026; originally announced July 2026.

Comments: 44 pages and 10 figures

Report number: LIGO-P2500781

Can We Find the Emission Mechanism Behind the Extremely Bright GRB 230812B?

Authors: Utkarsh Pathak, Sameer K. Patil, Hitesh Tanenia, Yashowardhan Rai, Viswajeet Swain, Anuraag Arya, Gaurav Waratkar, Anirudh Salgundi, Mehul Goyal, Varun Bhalerao, Igor Andreoni, Sarah Antier, G. C. Anupama, Sudhanshu Barway, Dipankar Bhattacharya, Tanmoy Chattopadhyay, Michael Coughlin, Anirban Dutta, Christoffer Fremling, Nidhal Guessoum, Mansi Kasliwal, Oliver J. Roberts, Gokul P. Srinivasaragavan, Rishabh Singh Teja, Santosh Vadawale , et al. (1 additional authors not shown)

Abstract: GRB 230812B is a bright long-duration GRB with a luminous, long-lived afterglow and an AstroSat/CZTI polarization measurement during the prompt phase, enabling a joint study of its prompt spectral evolution, polarization, and broadband afterglow. Time-resolved spectroscopy of the prompt emission shows that during the rising phase, the low-energy Band-function index exceeds the synchrotron line of… ▽ More GRB 230812B is a bright long-duration GRB with a luminous, long-lived afterglow and an AstroSat/CZTI polarization measurement during the prompt phase, enabling a joint study of its prompt spectral evolution, polarization, and broadband afterglow. Time-resolved spectroscopy of the prompt emission shows that during the rising phase, the low-energy Band-function index exceeds the synchrotron line of death, favoring the presence of an additional thermal component. At later times, from $T_0+2$ s to $T_0+32$ s, the prompt spectra are consistent with predominantly non-thermal emission. Polarization analysis of the prompt emission in the $300$-$600$ keV band yields a marginal lower limit on the polarization fraction of $Π\gtrsim 50\%$ at the $1σ$ level. The long X-ray monitoring of the afterglow shows no jet break over the observed baseline. Multiwavelength afterglow modeling favors a wide jet with an inferred half-opening angle of $θ_j = 15^{+6}_{-4}$ degrees observed close to the jet axis with a viewing angle of $θ_v = 0.9^{+1.8}_{-0.6}$ degrees. The inferred circumburst density is low, $n_0 = 1.2^{+0.3}_{-0.1}\times10^{-4}\,\textrm{cm}^{-3}$, and the isotropic-equivalent kinetic energy of the jet is $E_{{\rm k}, iso} = 4.0^{+1.5}_{-0.8} \times 10^{53}$ erg. Taken together, the prompt spectral evolution favors an early phase with a thermal contribution followed by a later phase dominated by non-thermal emission. The polarization constraint in the late prompt phase is consistent with synchrotron emission, although a higher-significance polarization measurement will be required to robustly constrain the magnetic-field geometry and the relative contribution of photospheric emission. △ Less

Submitted 20 July, 2026; originally announced July 2026.

Comments: 16 pages, 7 figures, 3 Tables. Submitted to ApJ

Catching Disguised Transients with ASTRANet: Anomaly-Aware Spectroscopic Classification and Conformal Calibration

Authors: Argyro Sasli, Maojie Xu, Alexandra Junell, Hailey Markoff, Avyukt Raghuvanshi, Felipe F. Nunes, Theophile Jegou Du Laz, Jesper Sollerman, Christoffer Fremling, Drew Oldag, Antoine Le Calloch, Sushant Sharma Chaudhary, Sneha Maharjan, Maxine West, Benny Border, Nabeel Rehemtulla, Richard Dekany, Joahan Castaneda Jaimes, Russ R. Laher, Reed Riddle, Mansi M. Kasliwal, Matthew J. Graham, Ashish A. Mahabal, Michael W. Coughlin

Abstract: Time-domain surveys discover thousands of transients per year, but the spectroscopic identification of rare and physically peculiar objects remains rate-limited by closed-set classifiers that confidently assign every input to a known class -- including spectra that genuinely belong to no known class. We present the \texttt{ASTRANet} framework, a confidence-aware infrastructure for spectroscopic tr… ▽ More Time-domain surveys discover thousands of transients per year, but the spectroscopic identification of rare and physically peculiar objects remains rate-limited by closed-set classifiers that confidently assign every input to a known class -- including spectra that genuinely belong to no known class. We present the \texttt{ASTRANet} framework, a confidence-aware infrastructure for spectroscopic transient classification built around three coupled modules: a hierarchical spectral classifier that operates directly on observer-frame spectra without requiring host-galaxy redshift or spectral phase as inputs; an anomaly detection layer (\texttt{ASTRANet-Sentinel}) that non-linearly combines $16$ embedding-space anomaly scores spanning four physically motivated families; and a conformal uncertainty quantification layer (\texttt{ASTRANet-CP}). We validate the framework on a held-out evaluation set of $289$ rare and out-of-taxonomy transients spanning $11$ classes deliberately excluded from training, chosen to span the full physical diversity of the rare-anomaly population: AGN-related outliers, GRB-related events, gap transients, novae, and peculiar supernovae. Through five astrophysically distinct failure modes of closed-set classifiers, we show that classifier-internal uncertainty and embedding-based anomaly detection are structurally complementary axes of confidence rather than alternative implementations of the same estimator. We further introduce AD-stratified Mondrian conformal prediction (AD-MCP) within \texttt{ASTRANet-CP}, achieving uniform conditional coverage across anomaly-score strata where vanilla Mondrian under-covers in the operational regime. This establishes the methodological infrastructure for confidence-aware spectroscopic discovery in the Vera C.\ Rubin Observatory era. △ Less

Submitted 8 July, 2026; originally announced July 2026.

Comments: 31 pages, 6 figures, 19 tables

Sub-Torque-Balance Upper Limits on Continuous Gravitational Waves from Scorpius X-1

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, the Precision Ephemerides for Gravitational-Wave Searches, Project, :, A. G. Abac, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, C. Adamcewicz, S. Adhicary, D. Adhikari, N. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend , et al. (1814 additional authors not shown)

Abstract: We present the results of a search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO data from the first part of the fourth LIGO-Virgo-KAGRA observing run. By applying the resampling version of the cross-correlation pipeline to search for signal frequencies $f_0$ between $25$ and $200\un{Hz}$ (corresponding to neutron star spin frequencies of $12.5$ to… ▽ More We present the results of a search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO data from the first part of the fourth LIGO-Virgo-KAGRA observing run. By applying the resampling version of the cross-correlation pipeline to search for signal frequencies $f_0$ between $25$ and $200\un{Hz}$ (corresponding to neutron star spin frequencies of $12.5$ to $100\un{Hz}$ for GW due to triaxiality, or $\sim15-20$ to $\sim120-150\un{Hz}$ for GW due to $r$-modes), we set upper limits below the standard torque balance level, independent of neutron star spin inclination, for $50\un{Hz}\lesssim f_0\lesssim200\un{Hz}$. While uncertainties in the modelling of torque and equation of state limit the strength of our inference, our results nonetheless argue against torque balance in this spin range for a neutron star described by a hadronic equation of state. The most sensitive upper limits on the gravitational wave amplitude $h_0$, at the upper end of the frequency band searched, approach $5\times10^{-26}$ marginalized over inclination angle and $2\times10^{-26}$ assuming the most favorable inclination. The marginalized upper limits correspond to a sensitivity depth of $70-75\un{Hz}^{-1/2}$, improving sensitivity considerably over previous searches. Expressed as constraints on the triaxial deformation of the neutron star, the limits correspond to an ellipticity of $3\times10^{-5}$ if the GW frequency $f_0$ is $75\un{Hz}$ and $3\times10^{-6}$ if $f_0=200\un{Hz}$, approaching deformations which could be supported by ordinary nuclear matter. Outliers from the search were ruled out as potential signals by a combination of hierarchical followup and analysis of additional data from later in the observing run. △ Less

Submitted 8 July, 2026; originally announced July 2026.

Comments: 30 pages, 7 figures

Report number: LIGO-P2500260-v10

nmma: An extended Bayesian framework for Nuclear Multimessenger Astronomy in the Era of Next-Generation Detectors

Authors: Henrik Rose, Hauke Koehn, Thibeau Wouters, Peter T. H. Pang, Mattia Bulla, Michael W. Coughlin, Tim Dietrich

Abstract: Context. The joint analysis of different (astro-)physical messengers, in particular gravitational-wave data and electromagnetic follow-up observations, allows us to establish, explore and deepen links between different physical fields. As new survey capacities and improved detection methods will lead to a significant increase in the number of multimessenger detections in the upcoming decades, effi… ▽ More Context. The joint analysis of different (astro-)physical messengers, in particular gravitational-wave data and electromagnetic follow-up observations, allows us to establish, explore and deepen links between different physical fields. As new survey capacities and improved detection methods will lead to a significant increase in the number of multimessenger detections in the upcoming decades, efficient and versatile software frameworks are essential to maximise the scientific outcome of such multimessenger studies. Aims. We present a major upgrade to the Nuclear Multimessenger Astronomy (nmma) framework, incorporating various recent developments in theoretical modelling and machine learning in a modularised and easily extendable Bayesian framework. For the first time, this allows direct sampling on nuclear parameters alongside gravitational-wave, kilonova and afterglow parameters. Methods. We combine fast surrogate models for electromagnetic transients with speed-ups from emulators that map nuclear parameters to macroscopic neutron-star properties. Additional acceleration methods for the evaluation of state-of-the-art waveform approximants enable full Bayesian analyses of multimessenger events at the speed required in the era of next-generation detectors. Results. We demonstrate the capabilities of the upgraded nmma framework through a series of representative applications. Reanalysing the 2017 multi-messenger detection of a neutron-star merger, we achieve 20- to 60-fold speed-ups while using more detailed physical models compared to previous studies. Moreover, we demonstrate for a hypothetical future detection how we can simultaneously constrain nuclear parameters and the Hubble parameter with robustly quantified uncertainties. △ Less

Submitted 3 July, 2026; originally announced July 2026.

Comments: 12(+8) pages, 6(+4) figures. To be submitted to A&A by next week, comments welcome

AI-enabled gravitational-waves searches for binary neutron stars at optimal sensitivity

Authors: Bhavya Gupta, Deep Chatterjee, William Benoit, Ethan Marx, Christina Reissel, Seiya Tsukamoto, Kyungseop Yoon, Michael W. Coughlin, Philip Harris, Erik Katsavounidis

Abstract: Gravitational Waves (GWs) represent the newest window of astronomy, furthering our understanding of compact objects like black holes and neutron stars in the Universe. The signal from two merging neutron stars is especially interesting since it brings the prospect of concordant electromagnetic and neutrino emissions. Such multi-messenger observations have a transformational impact on fundamental p… ▽ More Gravitational Waves (GWs) represent the newest window of astronomy, furthering our understanding of compact objects like black holes and neutron stars in the Universe. The signal from two merging neutron stars is especially interesting since it brings the prospect of concordant electromagnetic and neutrino emissions. Such multi-messenger observations have a transformational impact on fundamental physics, nuclear matter, astrophysics, and gravity. It was first witnessed in 2017 with the detection of the binary neutron star (BNS) merger GW170817. However, searching for BNS signals in real-time in the LIGO-Virgo-KAGRA (LVK) GW detectors presents a computational challenge, as the data streaming out must be matched against $\sim$ million reference waveforms, which requires up to a thousand CPU cores. We present a different approach using neural networks to learn the presence of a signal in the data. Our algorithm, called Aframe, was deployed in the LVK's fourth observing run and was the first artificial intelligence (AI)-enabled search to detect multiple binary black holes (BBHs) live. In this work, we demonstrate that the approach extends to the lower-mass BNS regime, and is the first AI-enabled search that achieves sensitivity comparable to matched-filter pipelines at lower computational and latency costs. The challenge of the longer-duration BNS signals is addressed by heterodyning the data, following which the network architecture used for BBHs is sufficient to distinguish signal versus background. We also show that this analysis requires a single non-flagship GPU for online deployment. Furthermore, the design and adoption of inference-as-a-service tools allow rapid offline analysis using a distributed pool of GPU resources. Hence, aside from the use case of rapid online data analysis, we also establish the use of Aframe for efficient archival data analysis. △ Less

Submitted 1 July, 2026; originally announced July 2026.

Leveraging Multimodality for Real-Time Classification of Transients and Variables found by the Zwicky Transient Facility

Authors: Ved G. Shah, Nabeel Rehemtulla, Adam A. Miller, Sushant Sharma Chaudhary, Michael W. Coughlin, Antoine Le Calloch, Matthew J. Graham, Joahan Castaneda Jaimes, Theophile Jegou du Laz, Ashish A. Mahabal, Frank J. Masci, Josiah Purdum, Reed Riddle, Jesper Sollerman, Anastasia Wei, Mansi M. Kasliwal

Abstract: Modern time-domain surveys such as the Zwicky Transient Facility (ZTF) generate hundreds of thousands of alerts each night, making real-time decisions for follow-up observations a central challenge in time-domain astronomy. Robust early classification is crucial for making informed decisions, but is hindered by sparse light curves and degeneracies between classes. In this work, we leverage multimo… ▽ More Modern time-domain surveys such as the Zwicky Transient Facility (ZTF) generate hundreds of thousands of alerts each night, making real-time decisions for follow-up observations a central challenge in time-domain astronomy. Robust early classification is crucial for making informed decisions, but is hindered by sparse light curves and degeneracies between classes. In this work, we leverage multimodality to substantially improve real-time classification and demonstrate the practicality of our approach by deploying our model on the ZTF alert stream. Building on the Online Ranked Astrophysical CLass Estimator (ORACLE), we introduce the ORACLE-2 models, which combine light curves, metadata, and images for real-time hierarchical classification. Using both real and simulated datasets, we show that incorporating additional modalities consistently improves classification performance. On observations from ZTF's Bright Transient Survey, our best-performing model, ORACLE-2 Omni, achieves a macro F1 score of 0.73 -- an improvement of up to 11% over models using light curves and metadata alone, and up to 40% over light-curve-only models, with the strongest gains realized at early times. To demonstrate applicability to the Legacy Survey of Space and Time, which will increase alert volume by more than an order of magnitude, we train a light curve + metadata variant on the simulated ELAsTiCC dataset. This model achieves a macro F1 score of 0.88, an improvement of up to 13% over the light-curve-only variant, matching the performance of other state-of-the-art models. Finally, we quantify the trade-offs between performance and throughput, identifying regimes where multimodal approaches offer the greatest benefit. These results show that combining multiple modalities improves early-time classification, enabling more effective triage of high-volume alert streams for current and future time-domain surveys. △ Less

Submitted 30 June, 2026; originally announced July 2026.

Comments: 29 Pages, 15 Figures, 8 Tables. Comments welcome

Constraints on Late-Time Flaring from Luminous Fast Blue Optical Transients using the Transiting Exoplanet Survey Satellite and the Zwicky Transient Facility

Authors: Rahul Jayaraman, Anna Y. Q. Ho, Michael M. Fausnaugh, Eran Ofek, Daniel A. Perley, Ruslan Konno, Martti Kristiansen, Tracy X. Chen, Michael W. Coughlin, Steven L. Groom, George Helou, K. -Ryan Hinds, Mansi M. Kasliwal, Zoë McGrath, Josiah N. Purdum, Argyro Sasli, Jesper Sollerman

Abstract: The Luminous Fast Blue Optical Transient (LFBOT) AT2022tsd exhibited minutes-timescale optical flares in the tens of days following the initial transient event, likely due to a central engine -- either an accreting black hole or a magnetar. In this paper, we use data from the Transiting Exoplanet Survey Satellite (TESS) and the Zwicky Transient Facility (ZTF) to constrain the occurrence of similar… ▽ More The Luminous Fast Blue Optical Transient (LFBOT) AT2022tsd exhibited minutes-timescale optical flares in the tens of days following the initial transient event, likely due to a central engine -- either an accreting black hole or a magnetar. In this paper, we use data from the Transiting Exoplanet Survey Satellite (TESS) and the Zwicky Transient Facility (ZTF) to constrain the occurrence of similar flares in the 12 (of 14) known LFBOTs that had observational coverage with TESS from tens of days to thousands of days after the transient's initial emission. We find seven flare-like signals at the locations of four unique LFBOTs; all seven can likely be attributed to a solar system object (SSO) moving through the TESS aperture. Assuming all seven flares arise from SSOs, for the LFBOT AT2024qfm we rule out flaring with a similar timescale (40--65 d) and luminosity ($νL_ν\sim10^{43}$ erg s$^{-1}$) as in AT2022tsd, while for AT2022tsd itself we rule out flares between 380--430 d after the initial transient that were as luminous as the earlier flares. This observation suggests that the engine power in AT2022tsd declined or shut off on a timescale of hundreds of days. We also find that there is no late-time activity detectable in TESS thousands of days after the prototype LFBOT, AT2018cow. We discuss our constraints on the duty cycle of such flaring and then present estimates for the number of minutes-duration flares detectable with ongoing and upcoming high-cadence ($\ll1$ d) wide-field surveys. △ Less

Submitted 15 June, 2026; originally announced June 2026.

Comments: 23 pages (including references), 8 figures, 6 tables. Submitted to ApJ

Method to get Better Sky Maps in a GstLAL Low-Latency Analysis

Authors: Prathamesh Joshi, Becca Ewing, Chad Hanna, Zach Yarbrough, Jolien D. E. Creighton, Shomik Adhicary, Pratyusava Baral, Amanda Baylor, Kipp Cannon, Sarah Caudill, Michael W. Coughlin, Bryce Cousins, Heather Fong, Richard N. George, Shaon Ghosh, Patrick Godwin, Reiko Harada, Yun-Jing Huang, Rachael Huxford, James Kennington, Soichiro Kuwahara, Alvin K. Y. Li, Ryan Magee, Duncan Meacher, Cody Messick , et al. (18 additional authors not shown)

Abstract: Modeled gravitational wave searches correlate the strain data with a bank of gravitational wave template waveforms to make detections of gravitational wave candidates, and these results are processed by downstream tools to calculate the likely sky location and distance of the source of the candidates. This is crucial for multi-messenger efforts, since it informs astronomers where to point their te… ▽ More Modeled gravitational wave searches correlate the strain data with a bank of gravitational wave template waveforms to make detections of gravitational wave candidates, and these results are processed by downstream tools to calculate the likely sky location and distance of the source of the candidates. This is crucial for multi-messenger efforts, since it informs astronomers where to point their telescopes to facilitate electromagnetic follow-up of the gravitational wave candidates. We present a novel method to improve the low-latency results of the GstLAL gravitational wave search pipeline, and thus improving sky location estimates of low-latency candidates. This method involves ingesting the GstLAL low-latency results, and performing a small targeted hierarchical search to recover the candidates with more accurate parameters, in a medium-latency timescale (few seconds to five minutes). To test our method, we perform a GstLAL low-latency analysis on forty days of data from the third observing run of LIGO, Virgo, and KAGRA, and show that our method improves the GstLAL results by 5.38% and the subsequent sky location results by 16.75% on average. In addition to this increase in precision, we also show that these results are more accurate as compared to the GstLAL results. This method has been adopted by GstLAL for the fourth observing run. △ Less

Submitted 8 June, 2026; originally announced June 2026.

Comments: 15 pages, 12 figures

A Multi-Wavelength View of the First Type Ic-BL Supernova with an Einstein Probe X-ray Shock Breakout

Authors: Jillian C. Rastinejad, Gokul Srinivasaragavan, Nikhil Sarin, Tanner O'Dwyer, S. Bradley Cenko, James K. Leung, Anya E. Nugent, Daniel A. Perley, Genevieve Schroeder, Shreya Anand, Tomas Ahumada, Igor Andreoni, Aleksandra Bochenek, Alessandra Corsi, Christoffer Fremling, Anna Y. Q. Ho, Mansi M. Kasliwal, Geoffrey Mo, Anirudh Salgundi, Kendall I. Sippy, J. Sollerman, Eric C. Bellm, Tracy X. Chen, Michael W. Coughlin, Michael C. Davis , et al. (14 additional authors not shown)

Abstract: In March 2026, the Einstein Probe (EP) discovered its most nearby (z = 0.0343) Fast X-ray Transient (FXT), EP260321a, the first EP FXT to provide a strong match to expectations for X-ray "shock breakout'" (SBO) emission. Here, we present our multi-wavelength follow-up campaign of EP260321a and its broad-line Type Ic (Ic-BL) supernova (SN) counterpart, SN2026gzf. We show that our radio follow-up ex… ▽ More In March 2026, the Einstein Probe (EP) discovered its most nearby (z = 0.0343) Fast X-ray Transient (FXT), EP260321a, the first EP FXT to provide a strong match to expectations for X-ray "shock breakout'" (SBO) emission. Here, we present our multi-wavelength follow-up campaign of EP260321a and its broad-line Type Ic (Ic-BL) supernova (SN) counterpart, SN2026gzf. We show that our radio follow-up extending over 5.8 - 54.5 days post-FXT rules out an on-axis jet counterpart of isotropic-equivalent kinetic energy $E_{K} \gtrsim 10^{49}$ erg for circumburst densities $n > 10^{-2}~{\rm cm}^{-3}$ and assuming microphysical parameters $ε_e = ε_B = 0.1$. Our radio data also constrains a median mass-loss rate of $\dot{M} \lesssim 1.2 \times 10^{-5} M_{\odot}~{\rm yr}^{-1}$ for a Wolf-Rayet progenitor. In addition, we derive SN2026gzf's properties, including $^{56}$Ni mass, diffusion timescale, and expansion velocities, from our $\sim$nightly-cadence optical data and compare them with those of optically discovered Type Ic-BL SNe, finding that SN 2026gzf is well within the 90\% confidence interval across all properties. We further fit SN2026gzf's light curve and determine that combined emission from both interaction with CSM and $^{56}$Ni radioactive decay provides the best fit with plausible model parameters. Finally, using the rate of Ic-BL SNe from the ZTF Bright Transient Survey and assuming all Type Ic-BL SNe produce EP260321a-like FXTs, we infer an expected rate of EP-detected SBOs of 4.4 - 16 year$^{-1}$. This is inconsistent at the 90% confidence level with current EP detection rates, potentially indicating that most Type Ic-BL SNe produce less luminous X-ray SBO signals compared to EP260321a. △ Less

Submitted 15 July, 2026; v1 submitted 8 June, 2026; originally announced June 2026.

Comments: Version accepted to ApJL with minor changes

Narrow spectral artifact investigation and mitigation in LIGO data from the fourth LIGO-Virgo-KAGRA observing run

Authors: E. Goetz, A. Neunzert, A. M. Knee, A. Calafat, X. Fan, J-R. Mérou, K. A. Pham, T. Starkman, N. Aggarwal, Z. Bhalla, P. Baxi, J. Bayley, Y. Bu, J. B. Carlin, P. Charlton, X. Chen, G. Cheng, T. Cheunchitra, N. Christensen, A. Claveus, C. M. Compton, M. W. Coughlin, F. De Lillo, L. Dunn, S. E. Dwyer , et al. (237 additional authors not shown)

Abstract: We present efforts to identify, characterize, and mitigate narrow spectral artifacts in LIGO detector data during the fourth LIGO-Virgo-KAGRA observing run. Narrow spectral artifacts in gravitational-wave detectors are non-astrophysical noise sources that can degrade searches for narrowband persistent gravitational waves. Identifying and, where possible, mitigating these noise sources is one of th… ▽ More We present efforts to identify, characterize, and mitigate narrow spectral artifacts in LIGO detector data during the fourth LIGO-Virgo-KAGRA observing run. Narrow spectral artifacts in gravitational-wave detectors are non-astrophysical noise sources that can degrade searches for narrowband persistent gravitational waves. Identifying and, where possible, mitigating these noise sources is one of the core efforts of the LIGO Detector Characterization group. Key software tools have been updated and new tools deployed for the fourth LIGO-Virgo-KAGRA observing run to facilitate investigations and data monitoring. We discuss these tool upgrades, and present several identified narrowband artifacts that have been successfully investigated and mitigated in LIGO data. Regardless of whether artifacts are mitigated or not, narrowband persistent gravitational-wave searches require information on which frequency bands contain non-astrophysical artifacts. Minimizing the number of bands containing non-astrophysical artifacts is essential to maximize the potential for discovery of a new class of gravitational-wave signals. △ Less

Submitted 23 July, 2026; v1 submitted 2 June, 2026; originally announced June 2026.

Comments: 40 pages, 15 figures

Toward decision-aware AI for LSST-scale time-domain astronomy

Authors: C. R. Bom, A. Mahabal, F. Bianco, P. Darc, B. Fraga, R. Bonito, S. Chaini, M. W. Coughlin, S. Dillmann, F. Fontinele Nunes, A. Gomboc, N. Hernitschek, X. Li, F. Z. Majidi, A. I. Malz, A. Melandri, V. Petrecca, S. Piranomonte, M. Rabus, F. Ragosta, O. Razim, M. C. Romão, N. Sarin, A. Sasli, V. A. Srećković , et al. (5 additional authors not shown)

Abstract: The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will generate approximately (10^7) alerts per night, pushing time-domain astronomy beyond pipelines that treat discovery as a static labeling problem. We argue that LSST is better understood as a partially observed dynamical environment, in which scientific return depends on the quality of follow-up decisions made under uncerta… ▽ More The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will generate approximately (10^7) alerts per night, pushing time-domain astronomy beyond pipelines that treat discovery as a static labeling problem. We argue that LSST is better understood as a partially observed dynamical environment, in which scientific return depends on the quality of follow-up decisions made under uncertainty and finite observational resources. The central challenge is therefore to maintain evolving, uncertainty-aware representations of astrophysical sources and to select actions that maximize long-term scientific value. We propose that foundation models trained on heterogeneous time-domain data can learn survey-scale representations of source state, while decision-theoretic policies support principled, auditable allocation of follow-up resources. Embedded within human-supervised agentic systems, these components position AI as part of the operational inference loop rather than as a downstream predictive tool. The way such systems represent belief, optimize utility, and expose their reasoning will shape observational efficiency, the distribution of scientific agency, including who participates in discovery and the scientific questions that receive priority. △ Less

Submitted 3 June, 2026; originally announced June 2026.

Comments: 12 pages, 2 figures, submitted as a Perspective article

Early Multiwavelength Observations of AT 2026fgk: The Luminous Afterglow to Sub-luminous GRB 260310A, Identified Independently of a Gamma-ray Trigger

Authors: K. -R. Hinds, A. Y. Q. Ho, Y. Wagh, R. Jayaraman, D. A. Perley, G. Waratkar, A. Bochenek, B. P. Gompertz, C. Fremling, J. Rastinejad, N. Sarin, G. Schroeder, R. A. Perley, G. P. Srinivasaragavan, K. Ackley, T. Ahumada, M. F. Aller, I. Andreoni, A. Aryan, S. Belkin, E. C. Bellm, S. Ben-Ami, T. de Boer, M. Bremer, R. P. Breton , et al. (81 additional authors not shown)

Abstract: The origins of sub-luminous ($L_\mathrm{γ,\mathrm{iso}} < 10^{49.5}$\,erg\,s$^{-1}$) gamma-ray bursts (GRBs) associated with broad-lined Type~Ic supernovae (Ic-BL SNe) are poorly understood, in part due to the low discovery rate and faint afterglows. Here we present the identification of the optical afterglow of Fermi-GBM-detected GRB\,260310A (AT\,2026fgk) as a rapidly rising ($>1\,$mag\,d… ▽ More The origins of sub-luminous ($L_\mathrm{γ,\mathrm{iso}} < 10^{49.5}$\,erg\,s$^{-1}$) gamma-ray bursts (GRBs) associated with broad-lined Type~Ic supernovae (Ic-BL SNe) are poorly understood, in part due to the low discovery rate and faint afterglows. Here we present the identification of the optical afterglow of Fermi-GBM-detected GRB\,260310A (AT\,2026fgk) as a rapidly rising ($>1\,$mag\,d$^{-1}$), red ($g-r=0.4$\,mag) transient using the Gravitational-wave Optical Transient Observatory, Large Array Survey Telescope, and Zwicky Transient Facility (ZTF) data streams. We present multiwavelength follow-up observations from the first 50\,days, which reveal that GRB 260310A/AT\,2026fgk was sub-luminous ($L_\mathrm{γ,iso}=10^{48.8}\,$erg\,s$^{-1}$); it was the most nearby ($z=0.153$) afterglow identified blindly by an optical survey; and that it is one of the brightest afterglows ever observed at X-ray, optical, and radio (cm to mm) wavelengths. We spectroscopically confirm an underlying Ic-BL SN with properties typical of GRB-SNe ($M_\mathrm{ej}\approx3\,M_\odot$, $E_{\rm K}\approx 10^{52}\,$erg). With basic modeling of the afterglow, including the long optical rise ($\approx10^{3}\,$s), we infer either a low initial Lorentz factor ($Γ_0\approx40$) or a slightly off-axis viewing angle ($\lesssim3^\circ$). The host galaxy's mass and star formation rate are similar to the hosts of other sub-luminous GRBs. ZTF's flux-limited survey gives a volumetric rate of AT\,2026fgk-like events of $0.30^{+1.37}_{-0.29}\,$Gpc\,$^{-3}$\,yr$^{-1}$, which is consistent with the on-axis, high luminosity ($L_{\rm γ,iso}>10^{49.5}$\,erg\,s$^{-1}$) long-GRB rate. The similarity in the rates strongly constrains the prevalence of low-$Γ_0$ bursts and the beaming of the initial relativistic material in GRBs. △ Less

Submitted 3 June, 2026; originally announced June 2026.

GWTC-5.0: Constraints on the Cosmic Expansion Rate and Modified Gravitational-wave Propagation

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1788 additional authors not shown)

Abstract: We employ 236 gravitational-wave (GW) sources in the fifth LIGO--Virgo--KAGRA Collaboration (LVK) Gravitational-Wave Transient Catalog (GWTC-5.0) to estimate the Hubble constant $H_0$. We compare the luminosity distance measured from GWs to the redshift inferred i) using features in the mass spectrum, and ii) using statistical host galaxy association. Probing the relationship between source lumino… ▽ More We employ 236 gravitational-wave (GW) sources in the fifth LIGO--Virgo--KAGRA Collaboration (LVK) Gravitational-Wave Transient Catalog (GWTC-5.0) to estimate the Hubble constant $H_0$. We compare the luminosity distance measured from GWs to the redshift inferred i) using features in the mass spectrum, and ii) using statistical host galaxy association. Probing the relationship between source luminosity distances and redshifts obtained in this way yields constraints on cosmological parameters. We estimate $H_0 = {71.7}_{-7.5}^{+9.4}\,{\text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}}$ (median with $68\%$ symmetric credible interval). This combines information from the source-frame mass distribution with the $H_0$ measurement from GW170817 and its electromagnetic counterpart as well as galaxy catalog information from Dark Energy Survey Year 6 (DES-Y6). We improve over the GWTC-4.0 measurement by using more GW sources, some with significantly smaller sky localization volumes, which leads to a reduction by $22.0\%$ of the $H_0$ uncertainty and a reconstructed mass distribution with lower uncertainties. We also constrain deviations from general relativity (GR) which affect GW propagation, specifically that modify the luminosity distance inferred from the GW signal. We find no departures from GR in parameterized tests of GW propagation. △ Less

Submitted 4 August, 2026; v1 submitted 26 May, 2026; originally announced May 2026.

Comments: main paper: 20 pages and 7 figures; total with appendices: 51 pages and 13 figures. This article draws heavily from the corresponding GWTC-4 article, arXiv:2509.04348

Report number: LIGO-P2600018

GWTC-5.0: Population Properties of Merging Compact Binaries

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1791 additional authors not shown)

Abstract: We present the population properties of merging compact binaries inferred using 267 mergers from the cumulative Gravitational-Wave Transient Catalog 5.0. As this data set contains no new sources with a neutron star, we primarily focus on the properties of the binary black hole mergers. We infer the merger rate of binary black holes with component masses between $2.5\,\mathrm{M}_\odot $ and… ▽ More We present the population properties of merging compact binaries inferred using 267 mergers from the cumulative Gravitational-Wave Transient Catalog 5.0. As this data set contains no new sources with a neutron star, we primarily focus on the properties of the binary black hole mergers. We infer the merger rate of binary black holes with component masses between $2.5\,\mathrm{M}_\odot $ and $200\,\mathrm{M}_\odot $ to be $27.5\text{--} 49.4 \, \mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$ (all intervals at $90\%$ credible levels) at redshift $z = 0.2$. We find evidence for a subpopulation of binary black hole mergers that host a rapidly spinning black hole (dimensionless spins $χ\sim 0.7$), consistent with signatures of hierarchical mergers. We find that these occur at two mass scales, the first at primary masses $\sim 10$--$20\,\mathrm{M}_\odot $ and the second above $\sim 45\,\mathrm{M}_\odot $, and we estimate their total rate at $z=0.2$ to be $0.2\text{--} 3.11 \, {\rm Gpc}^{-3} {\rm yr}^{-1}$. We infer that, above $40\,\mathrm{M}_\odot $, the mass distribution of the less massive (secondary) black hole declines more steeply than that of the more massive (primary) one. This is consistent with a flatter mass-ratio distribution and indicates the prevalence of unequal-mass binaries with large primary masses. We find evidence for two features in the black hole mass spectrum: a peak around $10\,\mathrm{M}_\odot $ and a change of slope at around $35\,\mathrm{M}_\odot $. Black holes of $\sim 35\,\mathrm{M}_\odot $ pair preferentially with companions of similar mass. Additionally, we find that the effective inspiral spin distribution of binary black holes is asymmetric about zero, based on which we infer that at least $9 \%$ of mergers occur in channels with some preference for spin-orbit alignment. We find evidence that... △ Less

Submitted 1 July, 2026; v1 submitted 26 May, 2026; originally announced May 2026.

Comments: Abstract truncated in Arxiv metadata. The paper appendices draw heavily from the corresponding GWTC-4 article, arXiv:2508.18083

Report number: LIGO-P2600045

GWTC-5.0: Observations from the Second Part of the Fourth LIGO-Virgo-KAGRA Observing Run and Updates to the Gravitational-Wave Transient Catalog

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1805 additional authors not shown)

Abstract: Version 5.0 of the Gravitational-Wave Transient Catalog (GWTC-5.0) adds new candidates detected by the LIGO Virgo KAGRA network of observatories through the second part of the fourth observing run (O4b: 2024 April 10 15:00:00 to 2025 January 28 17:00:00 UTC) and four days of the preceding engineering run (2024 April 6 to 2024 April 10). We find 161 compact binary coalescence candidates that are id… ▽ More Version 5.0 of the Gravitational-Wave Transient Catalog (GWTC-5.0) adds new candidates detected by the LIGO Virgo KAGRA network of observatories through the second part of the fourth observing run (O4b: 2024 April 10 15:00:00 to 2025 January 28 17:00:00 UTC) and four days of the preceding engineering run (2024 April 6 to 2024 April 10). We find 161 compact binary coalescence candidates that are identified by at least one of our search algorithms with a probability of astrophysical origin $p_\mathrm{astro} \geq 0.5$ and that are not vetoed during event validation. We also provide detailed source property measurements for 104 candidates that have a false-alarm rate < 1yr$^{-1}$. Based on the inferred component masses, all these candidates are consistent with signals from binary black holes. Median inferred component masses in the new candidates range from 5.14$M_\odot$ (GW241109_115924) to 70$M_\odot$ (GW241116_151753). Improvements in detector sensitivity allow us to observe compact binary coalescences with increasing clarity: 5 binary-black-hole signals have network signal-to-noise ratio exceeding 30, with a maximum to date of 76.9 for GW250114_082203. Such loud signals enable more precise studies of properties of their astrophysical sources and tests of general relativity. We also present updated results up to the first part of the fourth observing run, identifying 229 candidates. This brings the total number of transients in the cumulative GWTC having $p_\mathrm{astro} \geq 0.5$ to 390, further expanding the size of the catalog and our view of the gravitational-wave universe. △ Less

Submitted 23 June, 2026; v1 submitted 26 May, 2026; originally announced May 2026.

Comments: main paper: 30 pages, 8 figures; total with appendices: 43 pages, 9 figures. This article draws heavily from the corresponding GWTC-4 article, arXiv:2508.18082

Report number: LIGO-P2600152

GWTC-5.0: Methods for Identifying and Characterizing Gravitational-wave Transients

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1800 additional authors not shown)

Abstract: The Gravitational-Wave Transient Catalog (GWTC) is a collection of candidate gravitational-wave transient signals identified and characterized by the LIGO-Virgo-KAGRA Collaboration. Producing the contents of the GWTC from detector data requires complex analysis methods. These comprise techniques to model the signal; identify the transients in the data; evaluate the quality of the data and mitigate… ▽ More The Gravitational-Wave Transient Catalog (GWTC) is a collection of candidate gravitational-wave transient signals identified and characterized by the LIGO-Virgo-KAGRA Collaboration. Producing the contents of the GWTC from detector data requires complex analysis methods. These comprise techniques to model the signal; identify the transients in the data; evaluate the quality of the data and mitigate possible instrumental issues; infer the parameters of each transient; compare the data with the waveform models for compact binary coalescences, and handle the large amount of results associated with all these different analyses. In this paper, we describe the methods employed to produce the catalog's fifth release, GWTC-5.0, focusing on the analysis of the second part of the fourth observing run of LIGO, Virgo and KAGRA. △ Less

Submitted 23 June, 2026; v1 submitted 26 May, 2026; originally announced May 2026.

Comments: This article supersedes the corresponding GWTC-4 article, arXiv:2508.18081

Report number: LIGO-P2600166

GWTC-5.0: An Introduction to Version 5.0 of the Gravitational-Wave Transient Catalog

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1800 additional authors not shown)

Abstract: The Gravitational-Wave Transient Catalog (GWTC) is a collection of short-duration (transient) gravitational-wave signals identified by the LIGO-Virgo-KAGRA Collaboration in gravitational-wave data produced by the eponymous detectors. The catalog provides information about the identified candidates, such as the arrival time and amplitude of the signal and properties of the signal's source as inferr… ▽ More The Gravitational-Wave Transient Catalog (GWTC) is a collection of short-duration (transient) gravitational-wave signals identified by the LIGO-Virgo-KAGRA Collaboration in gravitational-wave data produced by the eponymous detectors. The catalog provides information about the identified candidates, such as the arrival time and amplitude of the signal and properties of the signal's source as inferred from the observational data. GWTC is the release of this dataset and version 5.0 extends the catalog to include observations made during the second part of the fourth LIGO-Virgo-KAGRA observing run up until 2025 January 28. This paper marks an introduction to a collection of articles related to this version of the catalog, GWTC-5.0. This update significantly increases the number of detected merging binary systems of black holes and neutron stars to over 300, enabling many follow-up studies toward understanding the gravitational-wave universe. The collection of articles accompanying the catalog provides documentation of the methods used to analyze the data, summaries of the catalog of events, observational measurements drawn from the population, and detailed discussions of selected candidates. △ Less

Submitted 23 June, 2026; v1 submitted 26 May, 2026; originally announced May 2026.

Comments: This article supersedes the corresponding GWTC-4 article, arXiv:2508.18080

Report number: LIGO-P2500701

Open Data from LIGO, Virgo, and KAGRA through the Second Part of the Fourth Observing Run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1787 additional authors not shown)

Abstract: LIGO, Virgo, KAGRA, and GEO 600 form a network of gravitational-wave observatories. Data and analysis results from this network are made publicly available through the Gravitational Wave Open Science Center (GWOSC). This paper describes open data from this network, including the addition of data from the second part of the fourth observing run (O4b) and selected periods from the preceding engineer… ▽ More LIGO, Virgo, KAGRA, and GEO 600 form a network of gravitational-wave observatories. Data and analysis results from this network are made publicly available through the Gravitational Wave Open Science Center (GWOSC). This paper describes open data from this network, including the addition of data from the second part of the fourth observing run (O4b) and selected periods from the preceding engineering run (ER16), which were collected from times spanning April 6th, 2024 to January 28th, 2025. The public data set includes calibrated strain time series for each instrument, data from additional channels used for noise subtraction and detector characterization, and new analysis data products in the online GWOSC release associated with version 5.0 of the Gravitational-Wave Transient Catalog. △ Less

Submitted 17 June, 2026; v1 submitted 26 May, 2026; originally announced May 2026.

Comments: This version has an updated author list, updated references and some text has been modified to improve clarity. This article draws heavily from the corresponding O4a article, arXiv:2508.18079

Report number: LIGO-P2600085

skysurvey: a pure python package to simulate the transient sky

Authors: M. Rigault, M. Ginolin, L. Dellazzeri, B. Popovic, J. O. Hjortlund, A. Gilles-Lordet, S. Conseil, M. Coughlin, F. Ruppin, M. Smith, A. Townsend, A. Trigui, C. Barjou-Delayre, R. Kebadian, J. Nordin

Abstract: Accurate simulation of astronomical observations is a critical element for any modern analyses, be it to measure event rates, analyses population properties, validate or train pipelines, account for selection effects, or correct biases. We present a novel pure python package named skysurvey made to enable the user to quickly simulate astrophysical transients as observed by a survey. The package is… ▽ More Accurate simulation of astronomical observations is a critical element for any modern analyses, be it to measure event rates, analyses population properties, validate or train pipelines, account for selection effects, or correct biases. We present a novel pure python package named skysurvey made to enable the user to quickly simulate astrophysical transients as observed by a survey. The package is structured to make the implementation of any complex population modeling fast and easy. The skysurvey package relies on three core objects: a Target, that models how an astrophysical target exists in nature, a Survey, that specifies how the sky has been observed and, a DataSet that combine these two to generate data as they would have been acquired. In addition, we present a side stand-alone package named modeldag that contains the core structure that simplifies the parameter modeling. We present in this paper how skysurvey is structured and we clearly illustrate how the code can straightforwardly be used to simulate complex populations, such as Type Ia Supernovae with varying color-brightness $β$ term. We also illustrate how the package can be made to replicate the rate and redshift distribution of the ZTF SNe Ia DR2 dataset. The skysurvey package, already used in recent scientific publications, is now ready for general usage and paves the way for future use of simulations such as simulation based inference. △ Less

Submitted 25 May, 2026; originally announced May 2026.

Comments: Submitted to Astronomy & Astrophysics

SN~2018erx: A fast-evolving, dust-reddened Type Icn supernova with broad C II emission lines

Authors: Kaustav K. Das, Anjasha Gangopadhyay, Mansi M. Kasliwal, Jesper Sollerman, Daniel A. Perley, Steve Schulze, Christoffer Fremling, Michael W. Coughlin, Kishalay De, Avishay Gal-Yam, Ariel Goobar, Matthew Graham, Frank J. Masci, Takashi J. Moriya, Josiah Purdum, Sam Rose, Ben Rusholme, Nikhil Sarin, Roger Smith, Daichi Tsuna

Abstract: We present the discovery and characterization of SN~2018erx (ZTF18abkmbpy), a fast-evolving, unusually red, interacting stripped-envelope supernova. Spectroscopically, SN~2018erx shows broad \ion{C}{2} emission with characteristic widths of $\sim\!3800$~km~s$^{-1}$, consistent with interaction with carbon-rich circumstellar material and a Type~Icn core-collapse SN classification. Photometrically,… ▽ More We present the discovery and characterization of SN~2018erx (ZTF18abkmbpy), a fast-evolving, unusually red, interacting stripped-envelope supernova. Spectroscopically, SN~2018erx shows broad \ion{C}{2} emission with characteristic widths of $\sim\!3800$~km~s$^{-1}$, consistent with interaction with carbon-rich circumstellar material and a Type~Icn core-collapse SN classification. Photometrically, it evolves rapidly, rising from half-maximum to peak in 2.1~d and declining back in 3.1~d. Semi-analytical CSM-interaction modeling favors a compact, shell-like CSM with $M_{\rm CSM}\approx0.3\,M_\odot$, $R_0\approx0.7$~AU, and a low ejecta mass of $M_{\rm ej}\approx0.11\,M_\odot$. The radioactive yield is also small, with $M_{\rm Ni}\lesssim(3$--$5)\times10^{-3}\,M_\odot$, placing SN~2018erx at the low end of the H-poor distribution. At +29~d after peak, we detect a near-infrared excess consistent with pre-existing local circumstellar dust, with $M_{\rm d}\sim10^{-6}$--$10^{-5}\,M_\odot$. Together, the rapid evolution, strong local reddening, carbon-rich emission, and dust point to a multi-component circumstellar environment: a dense inner interaction region from enhanced pre-SN mass loss and an outer dusty layer from an earlier mass-loss episode roughly $10$--$200$~yr before core collapse. These properties favor an ultra-stripped core-collapse explosion of a low-mass He star in a binary system, with fallback-modified Wolf--Rayet collapse or merger-driven mass loss remaining possible alternatives. SN~2018erx provides rare insight into the mass-loss history of stripped-envelope SNe and suggests that dust-enshrouded explosions of this kind may be underrepresented in optical surveys. △ Less

Submitted 23 May, 2026; originally announced May 2026.

Comments: Submitted to PASP. Comments are welcome

Low-Luminosity Type IIP Supernovae from the Zwicky Transient Facility Census of the Local Universe. III: Hunting for electron-capture supernovae using nebular spectroscopy

Authors: Kaustav K. Das, Anders Jerkstrand, Mansi M. Kasliwal, Jesper Sollerman, Christoffer Fremling, Steve Schulze, Avishay Gal-Yam, Tomas Ahumada, Shreya Anand, Bart van Baal, Michael W. Coughlin, Sofia Covarrubias, Richard Dekany, Nicholas Earley, W. V. Jacobson-Galán, Joahan Castaneda Jaimes, Frank J. Masci, Yu-Jing Qin, Reed Riddle, Sam Rose, Yashvi Sharma

Abstract: Electron-capture supernovae (ECSNe) may arise from ONeMg-core collapse in super-asymptotic giant branch (sAGB) stars near the low-mass core-collapse limit ($\approx\!8$--$10$\,\Msun). At early times, models predict that ECSNe resemble low-mass red supergiant iron-core-collapse SNe (FeCCSNe), making the two channels difficult to distinguish. Nebular spectroscopy, however, can reveal differences in… ▽ More Electron-capture supernovae (ECSNe) may arise from ONeMg-core collapse in super-asymptotic giant branch (sAGB) stars near the low-mass core-collapse limit ($\approx\!8$--$10$\,\Msun). At early times, models predict that ECSNe resemble low-mass red supergiant iron-core-collapse SNe (FeCCSNe), making the two channels difficult to distinguish. Nebular spectroscopy, however, can reveal differences in ejecta composition. We present a systematic sample of nebular spectra of 19 low-luminosity Type IIP (LLIIP) SNe from the ZTF CLU survey, obtained 115$-$450\,d after explosion. Their low velocities expose narrow lines blended in brighter SNe, which we identify and model to constrain progenitor properties. We find a strong correlation between the FWHM of H\,\textsc{i}\,$λ$6563 and peak luminosity, showing that LLIIP SNe occupy the low-energy end of the core-collapse population, but no correlation with plateau duration, suggesting that envelope and core properties are not tightly linked. Only one SN reaches the extremely low H\,\textsc{i}\,$λ$6563 widths predicted for the weakest $\sim$9\,M$_\odot$ explosion models, implying that such low-energy events are intrinsically rare. Combining our sample with 118 literature nebular spectra of Type II SNe, we infer an IMF slope of $2.1\pm1.2$. We also introduce an `ECSN score'' based on the absence of He- and O-shell emission lines, and identify two plausible ECSN candidates, SN~2023bvj and SN~2024btj. However, neither shows the extremely narrow nebular lines predicted by current ECSN models. If ECSNe arise predominantly through the LLIIP channel, we infer an upper limit on the ECSN rate of $\lesssim (5$--$8)\times10^{2}\,\mathrm{Gpc^{-3}\,yr^{-1}}$, corresponding to a narrow sAGB progenitor mass window of $ΔM_{\rm sAGB} \lesssim 0.02$--$0.06\,\mathrm{M_\odot}$. △ Less

Submitted 18 May, 2026; originally announced May 2026.

Comments: Submitted to PASP. Comments are welcome

Hyrax: An Extensible Framework for Rapid ML Experimentation and Unsupervised Discovery in the Era of Rubin, Roman, and Euclid

Authors: Aritra Ghosh, Drew Oldag, Michael Tauraso, Andrew J. Connolly, Peter Ferguson, Derek Jones, Gourav Khullar, Argyro Sasli, Samarth Venkatesh, Gracia Wang, Maxine West, Dylan Berry, Neven Caplar, Colin Orion Chandler, Tanawan Chatchadanoraset, Michael W. Coughlin, Melissa DeLucchi, Alexandra Junell, Diego Miura, Felipe Fontinele Nunes, Wilson Beebe, Doug Branton, Sandro Campos, Liam Cunningham, Mi Dai , et al. (7 additional authors not shown)

Abstract: The NSF-DOE Vera C. Rubin Observatory, Roman Space Telescope, Euclid, and other next-generation surveys will deliver imaging, spectroscopic, and time-domain data at scales that increasingly shift the bottleneck in astronomical machine learning (ML) projects from model design to infrastructure. We present Hyrax, an open-source, modular, GPU-enabled Python framework that supports the full ML lifecyc… ▽ More The NSF-DOE Vera C. Rubin Observatory, Roman Space Telescope, Euclid, and other next-generation surveys will deliver imaging, spectroscopic, and time-domain data at scales that increasingly shift the bottleneck in astronomical machine learning (ML) projects from model design to infrastructure. We present Hyrax, an open-source, modular, GPU-enabled Python framework that supports the full ML lifecycle in astronomy: from data acquisition and training to inference and experiment comparison, with capabilities including multimodal dataset support, integrated vector databases for similarity search, and interactive two- and three-dimensional latent-space exploration for unsupervised discovery. We demonstrate Hyrax's versatility through five representative applications on real survey data: (i) unsupervised representation learning on $\sim 4\times10^5$ Rubin Legacy Survey of Space and Time (LSST) Data Preview 1 (DP1) galaxies, surfacing new merger and low-surface-brightness candidates missing from reference Euclid and Dark Energy Survey catalogs, while also isolating imaging artifacts -- all without labeled training data; (ii) hybrid density-based clustering for identifying cluster-scale gravitational lens candidates in DP1 data; (iii) multimodal early-time transient classification in the Zwicky Transient Facility leveraging light curves, spectra, images, and metadata; (iv) supervised false-positive filtering in shift-and-stack searches for distant solar system objects in the Dark Energy Camera Ecliptic Exploration Project survey; and (v) supervised detection of semi-resolved dwarf galaxies in Hyper Suprime-Cam and LSST-like imaging using synthetic source injection. Together, these results demonstrate that Hyrax provides astronomy-specific ML infrastructure that enables systematic discovery and rapid methodological iteration across next-generation astronomical surveys. △ Less

Submitted 18 May, 2026; originally announced May 2026.

Comments: 28 pages, 20 figures, submitted to AJ

GW240925 and GW250207: Astrophysical Calibration of Gravitational-wave Detectors

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, C. Adamcewicz, S. Adhicary, D. Adhikari, N. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith , et al. (1817 additional authors not shown)

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

Submitted 17 August, 2026; v1 submitted 12 May, 2026; originally announced May 2026.

Comments: 38 pages (11 pages author list, 7 pages main paper, 8 pages references, 12 pages supplemental material), 15 figures (4 in main text, 11 in supplemental material); data products available from https://doi.org/10.5281/zenodo.18600070

Report number: LIGO-P2500565

Journal ref: Phys. Rev. Lett. 137, 071401 (2026)

Electromagnetic Follow-up of the Sub-Solar Mass Gravitational Wave Candidate S251112cm: Kilonova Constraints and a Coincident IIb Supernova

Authors: Xander J. Hall, Tomas Ahumada, Julius Gassert, Antonella Palmese, Brian D. Metzger, Mansi M. Kasliwal, Mattia Bulla, Daniel Gruen, Robert Stein, Christoffer Fremling, Shreya Anand, Igor Andreoni, Malte Busmann, Tomás Cabrera, Ryan Christinzio, James Freeburn, Ignacio Magaña Hernandez, Lei Hu, Brendan O'Connor, Ji-an Jiang, Zhengyan Liu, Wen Zhao, Eric C. Bellm, David Cook, Michael W. Coughlin , et al. (3 additional authors not shown)

Abstract: On November 12th, 2025 the LIGO--Virgo--KAGRA (LVK) collaboration reported gravitational waves (GWs) from a compact object merger candidate (S251112cm) with at least one sub-solar mass component. Using the Dark Energy Camera (DECam), the Fraunhofer Telescope at Wendelstein Observatory (FTW), and the Zwicky Transient Facility (ZTF), we surveyed $56\%$ of the GW localization region beginning $2.4$~h… ▽ More On November 12th, 2025 the LIGO--Virgo--KAGRA (LVK) collaboration reported gravitational waves (GWs) from a compact object merger candidate (S251112cm) with at least one sub-solar mass component. Using the Dark Energy Camera (DECam), the Fraunhofer Telescope at Wendelstein Observatory (FTW), and the Zwicky Transient Facility (ZTF), we surveyed $56\%$ of the GW localization region beginning $2.4$~hours after the GW alert. We find no kilonova (KN) counterpart, and use radiative-transfer models to rule out $42\%$ (ZTF), $68\%$ (DECam), and $92\%$ (FTW) of the KN models as possible emission from this GW candidate. Within the recently proposed disk-fragmentation (``superkilonova'') model for generating sub-solar mass neutron star mergers from stellar core-collapse, the delay between the supernova explosion time and the GW merger time is estimated to be less than a few days. Searching this time window prior to the GW event, we identify and spectroscopically classify a IIb supernova (SN~2025adtq), with a spatial association odds ratio of $\log_{10}\mathcal{I} \approx 4.8$, a chance coincidence probability of ${\sim}2$--$9\%$, and an estimated explosion time ${\sim}2$ days prior to S251112cm. SN~2025adtq is the second Type~IIb supernova found in spatial and temporal coincidence with a sub-solar mass GW candidate, following the previously reported S250818k/SN~2025ulz association; jointly, we measure an odds ratio that favors the association hypothesis over the null, however, when conditioned on finding a coincident supernova by chance, the odds ratio disfavors association. Together, these results provide suggestive but inconclusive evidence for the superkilonova formation channel. △ Less

Submitted 11 May, 2026; originally announced May 2026.

Searches for Binary Mergers with Sub-solar Mass Components in Data from the First Part of LIGO--Virgo--KAGRA's Fourth Observing Run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, C. Adamcewicz, S. Adhicary, D. Adhikari, N. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith , et al. (1810 additional authors not shown)

Abstract: We report on a gravitational wave search for compact binary coalescences involving at least one component with mass between $0.2\,M_\odot$ to $1\,M_\odot$, and ratio of component masses between 0.1 and 1. The analysis uses data collected by the LIGO detectors between May 24 2023 15:00 UTC and January 16 2024 16:00 UTC. No statistically significant sub-solar mass candidates were identified by the p… ▽ More We report on a gravitational wave search for compact binary coalescences involving at least one component with mass between $0.2\,M_\odot$ to $1\,M_\odot$, and ratio of component masses between 0.1 and 1. The analysis uses data collected by the LIGO detectors between May 24 2023 15:00 UTC and January 16 2024 16:00 UTC. No statistically significant sub-solar mass candidates were identified by the participating search algorithms. We report the detection sensitivity of the current searches to the target sub-solar mass black hole population. With the absence of detections, we place upper limits on the merger rate of sub-solar mass black holes, ranging from 110 ${\rm Gpc^{-3}\,yr^{-1}}$ to 10000 ${\rm Gpc^{-3}\,yr^{-1}}$ at 90\% confidence. We constrain two illustrative dark matter scenarios that can form sub-solar mass compact objects with these searches: primordial black holes, and dark black holes forming in a dissipative dark matter model. For late-forming primordial black hole binaries, our search excludes the fraction of dark matter in primordial black holes to be $\leq 1$ only for masses above $0.9\,M_\odot$. In the early-formation scenario, we limit this fraction to be $\leq 7\%$ at $1\,M_\odot$, and $\leq 40\%a$ at $0.35\,M_\odot$. For the dissipative model, the excluded region in the parameter space of dark matter fraction in dark black holes and their minimum possible mass extends down to (0.9 to 1.2) $\times 10^{-5}$ at $1\,M_\odot$ with no constraints below $0.02\,M_\odot$. For the first time, we report the detection sensitivity of our searches to binaries with sub-solar mass neutron stars, and place the 90\% confidence merger rate limit at (570 to 710) ${\rm Gpc^{-3}\,yr^{-1}}$ for a population with component masses distributed uniformly down to $0.5\,M_\odot$. △ Less

Submitted 23 July, 2026; v1 submitted 6 May, 2026; originally announced May 2026.

Report number: LIGO-P2400387

EMBER: Machine-Learning Detection of Modulated Ion Acoustic Waves and Associated Core-Electron Heating in the Solar Wind with Parker Solar Probe

Authors: Argyro Sasli, Karish Seebaluck, Chris Colpitts, Michael Coughlin

Abstract: Modulated ion acoustic waves (IAWs) -- including triggered ion acoustic waves (TIAWs) and frequency-dispersed ion acoustic waves (FDIAWs) -- are increasingly recognized as efficient drivers of electron heating in the solar wind through nonlinear wave-particle interactions. Identification of these events in the Parker Solar Probe (PSP) FIELDS burst-mode archive has so far relied on expert visual in… ▽ More Modulated ion acoustic waves (IAWs) -- including triggered ion acoustic waves (TIAWs) and frequency-dispersed ion acoustic waves (FDIAWs) -- are increasingly recognized as efficient drivers of electron heating in the solar wind through nonlinear wave-particle interactions. Identification of these events in the Parker Solar Probe (PSP) FIELDS burst-mode archive has so far relied on expert visual inspection and does not scale to the full mission. We present EMBER (Electron heating from Modulated Burst-mode Event Recognition), an open-source pipeline that converts PSP FIELDS Digital Burst Memory (DBM) voltage bursts into log-scaled Fourier spectrograms and applies a multi-detector, background-only anomaly detection suite. The suite combines physics-motivated detectors, classical outlier detectors, and deep learning detectors. The EMBER ensemble recovers 93% of the anomalous events at 1% FAR (1 false positive per 100 held-out backgrounds). Coincident SWEAP/SPAN diagnostics show that flagged intervals exhibit core perpendicular electron temperatures above the adiabatic cooling expectation and elevated Te/Ti, reproducing the preferential-heating phenomenology established by prior manual studies without any use of electron temperatures in the detection step. △ Less

Submitted 30 April, 2026; originally announced May 2026.

Comments: 16 pages, 5 figures

Prospects for GRB Afterglow Discovery with the Eric and Wendy Schmidt Observatory System

Authors: James Freeburn, Dougal Dobie, Akash Anumarlapudi, Igor Andreoni, Brendan O'Connor, Eric Burns, Jonathan Carney, Hank Corbett, Michael W. Coughlin, Adam Goldstein, Anna Tartaglia, V. Ashley Villar

Abstract: Two time domain surveys, recently funded as part of the Eric and Wendy Schmidt Observatory System; the Argus Array, in the optical, and the Deep Synoptic Array (DSA), in the radio, will transform gamma-ray burst (GRB) science via the serendipitous discovery of hundreds of GRB afterglows per year. In this work, we simulate DSA and Argus observations of GRB afterglows. We find that, of the long-dura… ▽ More Two time domain surveys, recently funded as part of the Eric and Wendy Schmidt Observatory System; the Argus Array, in the optical, and the Deep Synoptic Array (DSA), in the radio, will transform gamma-ray burst (GRB) science via the serendipitous discovery of hundreds of GRB afterglows per year. In this work, we simulate DSA and Argus observations of GRB afterglows. We find that, of the long-duration GRBs (LGRBs) detected by the Fermi Gamma-ray Burst Monitor, (26$\pm$3)% will yield afterglow detections with Argus and (44$\pm$3)% with DSA, corresponding to rates of 51$\pm$5 and 86$\pm$8 per year respectively. We also compute rates for both upcoming and proposed GRB monitors; the forthcoming StarBurst Multi-messenger Pioneer, with 84$\pm$6 detections per year in Argus and 158$\pm$11 detections per year in DSA, and the Moon Burst Energetics All-sky Monitor (MoonBEAM) concept, with 90$\pm$9 per year in Argus and 160$\pm$14 per year in DSA. The observatory system will also detect 118$\pm$10 optical and 199$\pm$17 radio afterglows per year, independent of GRB triggers, exceeding the current annual rate with global follow-up. Afterglow counterparts to short-duration GRBs, originating from neutron star mergers, will be detected at 5-10% of the LGRB afterglow rate, which is promising for multi-messenger detections of gravitational wave sources and constraining the neutron star merger rate. The Argus Array, with its second-minute cadence, will detect afterglows before they peak $\sim$18% of the time, which will dramatically increase the sample of observed reverse shock and prompt optical emission. △ Less

Submitted 22 July, 2026; v1 submitted 15 April, 2026; originally announced April 2026.

Comments: 22 pages, 9 figures, 6 tables, accepted for publication in ApJ

A Natural $\gtrsim 100\times$ Telescope: Discovery of the Strongly Lensed Type II SN 2025mkn at $z=1.37$

Authors: Cameron Lemon, Ariel Goobar, Joel Johansson, Edvard Mörtsell, Steve Schulze, Igor Andreoni, Aleksandra Bochenek, Seán J. Brennan, Malte Busmann, Michael Coughlin, Kaustav K. Das, Suhail Dhawan, Christoffer Fremling, Anjasha Gangopadhyay, Daniel Gruen, Xander J. Hall, Anna Y. Q. Ho, Mansi M. Kasliwal, Daniel A. Perley, Mickael Rigault, Genevieve Schroeder, Mathew Smith, Jesper Sollerman, Jean J. Somalwar, Robert Stein , et al. (68 additional authors not shown)

Abstract: We present the discovery of SN 2025mkn, a gravitationally lensed Type II supernova. First detected as a blue transient in ZTF, 0.83$^{\prime\prime}$ from a $z=0.42$ elliptical galaxy, follow-up SNIFS/UH2.2m and LRIS/Keck spectra revealed absorption lines at $z=1.371$. Later JWST NIRCam imaging shows that the bright transient is a close pair of point sources separated by $\sim 0.07^{\prime\prime}$,… ▽ More We present the discovery of SN 2025mkn, a gravitationally lensed Type II supernova. First detected as a blue transient in ZTF, 0.83$^{\prime\prime}$ from a $z=0.42$ elliptical galaxy, follow-up SNIFS/UH2.2m and LRIS/Keck spectra revealed absorption lines at $z=1.371$. Later JWST NIRCam imaging shows that the bright transient is a close pair of point sources separated by $\sim 0.07^{\prime\prime}$, and a 30 times fainter counterimage opposite the lens, for which NIRSpec reveals strong H$α$ emission also at $z=1.371$. The light curves and spectra are consistent with the Type II supernova source being magnified $\gtrsim 100$ times, with $\sim 250$ required to reconcile its luminosity with that of nearby events such as SN 2023ixf. Lens models are consistent with such high magnifications, and always show that the faint image arrived first (undetected in earlier ZTF imaging), consistent with the later spectral phase of this fainter image. A fourth image is also predicted and possibly detected in the NIRSpec data. Light-curve-based time-delay measurements are not possible due to the first image being the faintest; however, the resolved NIRSpec spectra offer a future opportunity for time-delay cosmography through supernova phase measurements. △ Less

Submitted 9 April, 2026; originally announced April 2026.

Comments: 14 pages, 4 figures, 1 table, submitted to ApJL

The ZTF-ULTRASAT experiment: Characterizing the non-transients in ULTRASAT's high cadence survey

Authors: Daniel Warshofsky, Michael W. Coughlin, Theophile Jegou Du Laz, Anna Y. Q. Ho, S. Bradley Cenko, Andrew Drake, Jesper Sollerman, Argyro Sasli, Ben Rusholme, Frank J. Masci, Roger Smith, A. M. Krassilchtchikov, David Berge, Eran O. Ofek, Yossi Shvartzvald, Reed L. Riddle, Mansi M. Kasliwal, Matthew J. Graham, Eric C. Bellm

Abstract: The forthcoming launch of the Ultraviolet Transient Astronomy Satellite (ULTRASAT) will transform our understanding of the transient ultraviolet sky by increasing our ability to identify transients due to its unprecedented 204 deg2 field of view. While rapid (extragalactic) transients are a priority science area for the mission, flaring stars and AGN can often contaminate searches for such objects… ▽ More The forthcoming launch of the Ultraviolet Transient Astronomy Satellite (ULTRASAT) will transform our understanding of the transient ultraviolet sky by increasing our ability to identify transients due to its unprecedented 204 deg2 field of view. While rapid (extragalactic) transients are a priority science area for the mission, flaring stars and AGN can often contaminate searches for such objects. To prepare for these challenges, the Zwicky Transient Facility (ZTF)-ULTRASAT experiment observed five fields at high cadence over three nights, in close proximity to ULTRASAT's three northern high-cadence fields. A real-time filter identified seven transient candidates, of which five were persistent variable sources and two were spurious. Periods and amplitudes derived from the ZTF Source Classification Project (SCoPe) showed that three candidates were RR Lyrae stars with short periods and high amplitudes, while the remaining two displayed flaring behavior. We demonstrate that short-timescale, high-amplitude variables can systematically mimic transient alerts in high-cadence UV surveys, and we provide a concrete strategy to this contamination using pre-existing machine learning catalogs. △ Less

Submitted 28 April, 2026; v1 submitted 8 April, 2026; originally announced April 2026.

Comments: 12 pages,13 figures

Narrowband searches for continuous gravitational waves from known pulsars in the first two parts of the fourth LIGO--Virgo--KAGRA observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, C. Adamcewicz, S. Adhicary, D. Adhikari, N. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith , et al. (1831 additional authors not shown)

Abstract: Rotating non-axisymmetric neutron stars (NSs) are promising sources for continuous gravitational waves (CWs). Such CWs can, if detected, inform us about the internal structure and equation of state of NSs. Here, we present a narrowband search for CWs from known pulsars, for which an efficient and sensitive matched-filter search can be applied. Narrowband searches are designed to be robust to misma… ▽ More Rotating non-axisymmetric neutron stars (NSs) are promising sources for continuous gravitational waves (CWs). Such CWs can, if detected, inform us about the internal structure and equation of state of NSs. Here, we present a narrowband search for CWs from known pulsars, for which an efficient and sensitive matched-filter search can be applied. Narrowband searches are designed to be robust to mismatches between the electromagnetic (EM) and gravitational emissions, in contrast to fully targeted searches where the CW emission is assumed to be phase-locked to the EM one. In this work, we search for the CW counterparts emitted by 34 pulsars using data from the first and second parts of the fourth LIGO--Virgo--KAGRA observing run. This is the largest number of pulsars so far targeted for narrowband searches in the advanced detector era. We use the 5n-vector narrowband pipeline, which applies frequency-domain matched filtering. In previous searches, it covered a narrow range in the frequency -- frequency time derivative ($f$ -- $\dot{f}$) space. Here, we also explore a range in the second time derivative of the frequency $\ddot{f}$ around the value indicated by EM observations. Additionally, for the first time, we target sources in a binary system with this kind of search. We find no evidence for CWs and therefore set upper limits on the strain amplitude emitted by each pulsar, using simulated signals added in real data. For 20 analyses, we report an upper limit below the theoretical spin-down limit. The tightest constraint is for pulsar PSR J0534+2200 (the Crab pulsar), for which our strain upper limit on the CW amplitude is $\lesssim 2\%$ of its spin-down limit, corresponding to less than $0.04\%$ of the spin-down power being radiated in the CW channel. △ Less

Submitted 8 July, 2026; v1 submitted 26 March, 2026; originally announced March 2026.

Comments: 30 pages, 6 figures, submitted to ApJ

Report number: LIGO-P2500612

Journal ref: The Astrophysical Journal, 1005:221, 2026 July 10

Searches for Continuous Gravitational Waves from Supernova Remnants in the first part of the LIGO-Virgo-KAGRA Fourth Observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, D. Adhikari, N. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1742 additional authors not shown)

Abstract: We present results from directed searches for continuous gravitational waves from a sample of 15 nearby supernova remnants, likely hosting young neutron star candidates, using data from the first eight months of the fourth observing run (O4) of the LIGO-Virgo-KAGRA Collaboration. The analysis employs five pipelines: four semi-coherent methods -- the Band-Sampled-Data directed pipeline, Weave and t… ▽ More We present results from directed searches for continuous gravitational waves from a sample of 15 nearby supernova remnants, likely hosting young neutron star candidates, using data from the first eight months of the fourth observing run (O4) of the LIGO-Virgo-KAGRA Collaboration. The analysis employs five pipelines: four semi-coherent methods -- the Band-Sampled-Data directed pipeline, Weave and two Viterbi pipelines (single- and dual-harmonic) -- and PyStoch, a cross-correlation-based pipeline. These searches cover wide frequency bands and do not assume prior knowledge of the targets' ephemerides. No evidence of a signal is found from any of the 15 sources. We set 95\% confidence-level upper limits on the intrinsic strain amplitude, with the most stringent constraints reaching $\sim 4 \times 10^{-26}$ near 300 Hz for the nearby source G266.2$-$1.2 (Vela Jr.). We also derive limits on neutron star ellipticity and $r$-mode amplitudes for the same source, with the best constraints reaching $\lesssim 10^{-7}$ and $\lesssim 10^{-5}$, respectively, at frequencies above 400 Hz. These results represent the most sensitive wide-band directed searches for continuous gravitational waves from supernova remnants to date. △ Less

Submitted 2 April, 2026; v1 submitted 26 March, 2026; originally announced March 2026.

Journal ref: The Astrophysical Journal, 1004 204 (2026)

SN 2024iss: A Multi-Wavelength Exposé of a Type IIb Supernova with an Early-Time Ultraviolet Spectrum and Shock Breakout Constraints

Authors: Rujula Yete, Wynn Jacobson-Galan, Ferdinand Ferdinand, Luc Dessart, Mansi Kasliwal, Kyle Davis, Anthony Piro, Victoria Villar, Igor Andreoni, Katie Auchettl, Kenneth Chambers, Ryan Chornock, Michael Coughlin, David Coulter, Richard Dekany, Ryan Foley, Ori Fox, Lluís Galbany, Christa Gall, Anjasha Gangopadhyay, Jared Goldberg, Gowri Govindaraj, Xinze Guo, Mark Huber, Joahan Jaimes , et al. (27 additional authors not shown)

Abstract: We present multi-wavelength observations and a comprehensive analysis of the nearby (D$\sim$14 Mpc) Type IIb supernova (SN IIb) 2024iss. Observations of SN2024iss include an early ZTF detection at $\sim$40 minutes after first light and the earliest Hubble Space Telescope UV spectrum for a SN IIb to date at 7 days after first light. With the bolometric light curve and He-star models, we estimate an… ▽ More We present multi-wavelength observations and a comprehensive analysis of the nearby (D$\sim$14 Mpc) Type IIb supernova (SN IIb) 2024iss. Observations of SN2024iss include an early ZTF detection at $\sim$40 minutes after first light and the earliest Hubble Space Telescope UV spectrum for a SN IIb to date at 7 days after first light. With the bolometric light curve and He-star models, we estimate an ejecta mass range of $\sim 1.1-3.3~M_{\odot}$ and a $^{56}\textrm{Ni}$ mass of $0.11 \pm 0.01~M_{\odot}$. We fit shock-cooling emission models to the first peak in the light curve and estimate a progenitor radius of $100-320~R_{\odot}$ and a H-rich envelope mass of $0.07-0.46~M_{\odot}$. We also compared optical/UV spectra to binary progenitor model spectra, which indicate a stripped H-rich envelope mass of $0.19-0.28~M_{\odot}$. We use early-time X-ray detections to calculate CSM densities that are consistent with a progenitor mass-loss rate of $5\times10^{-4}~M_{\odot}$ ($v_w = 100~$km/s), corresponding to a period of significant mass ejection in the final ~2-5 years before core collapse. In the UV spectrum, we observe strong Mg II emission extending to $\sim15,000 ~$km/s as well as weak P-Cygni profiles of iron-group elements (e.g., Fe, Ti, Al, Ni) present in the outer SN ejecta during the end of shock cooling phase. We find that the overall spectroscopic evolution of SN2024iss is comparable to other SNe IIb, but that the increased brightness following the initial light curve peak is likely influenced by SN ejecta-CSM interaction. Finally, optical/NIR nebular spectroscopy of SN2024iss at $\sim 260-412~$ days reveals multi-peaked forbidden line profiles of O I and Mg I] indicative of inner ejecta asymmetry and/or clumping. We demonstrate the utility of a rich, multi-wavelength dataset for constraining the progenitor systems and explosion dynamics of SNe IIb. △ Less

Submitted 30 March, 2026; v1 submitted 20 March, 2026; originally announced March 2026.

Comments: 30 pages, 18 figures, submitted to ApJ

GWTC-4.0: Tests of General Relativity. III. Tests of the Remnants

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, D. Adhikari, N. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1757 additional authors not shown)

Abstract: This is the third paper of the set recording the results of the suite of tests of general relativity (GR) performed on the signals from the fourth Gravitational-Wave Transient Catalog (GWTC-4.0), where we focus on the remnants of the binary mergers. We examine for the first time 42 events from the first part of the fourth observing run of the LIGO-Virgo-KAGRA detectors, alongside events from the p… ▽ More This is the third paper of the set recording the results of the suite of tests of general relativity (GR) performed on the signals from the fourth Gravitational-Wave Transient Catalog (GWTC-4.0), where we focus on the remnants of the binary mergers. We examine for the first time 42 events from the first part of the fourth observing run of the LIGO-Virgo-KAGRA detectors, alongside events from the previous observation runs, restricting our analysis to the confident signals, which were measured in at least two detectors and that have false alarm rates $\le 10^{-3} \mathrm{yr}^{-1}$. This paper focuses on seven tests of the coalescence remnants. Three of these are tests of the ringdown and its consistency with the expected quasinormal mode spectrum of a Kerr black hole. Specifically, two tests analyze just the ringdown in the time domain, and the third test analyzes the entire signal in the frequency domain. Four tests allow for the existence of possible echoes arriving after the end of the ringdown, which are not expected in GR. We find overall consistency of the remnants with GR. When combining events by multiplying likelihoods (hierarchically), one analysis finds that the GR prediction lies at the boundary of the $98.6^{+1.4}_{-9.4}\%$ ($99.3^{+0.7}_{-4.5}\%$) credible region, an increase from $93.8^{+6.1}_{-20.0}\%$ ($94.9^{+4.4}_{-18.2}\%$) for GWTC-3.0. Here the ranges of values comes from bootstrapping to account for the finite number of events analyzed and suggest that some of the apparently significant deviation could be attributed to variance due to the finite catalog. Since the significance also decreases to 92.2% (96.2%) when including the more recent very loud event GW250114, there is no strong evidence for a GR deviation. We find no evidence for post-merger echoes in the events that were analyzed. (Abridged) △ Less

Submitted 19 March, 2026; originally announced March 2026.

Comments: As part of the Astrophysical Journal Letters Focus Issue on the Gravitational Wave Transient Catalog

Report number: LIGO-P2500067

GWTC-4.0: Tests of General Relativity. II. Parameterized Tests

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, D. Adhikari, N. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1763 additional authors not shown)

Abstract: In this second of three papers on tests of general relativity (GR) applied to the compact binary coalescence signals in the 4th Gravitational-Wave Transient Catalog (GWTC-4.0), we present the results of the parameterized tests of GR and constraints on line-of-sight acceleration (LOSA). We include events up to and including the 1st part of the 4th observing run (O4a) of the LIGO-Virgo-KAGRA detecto… ▽ More In this second of three papers on tests of general relativity (GR) applied to the compact binary coalescence signals in the 4th Gravitational-Wave Transient Catalog (GWTC-4.0), we present the results of the parameterized tests of GR and constraints on line-of-sight acceleration (LOSA). We include events up to and including the 1st part of the 4th observing run (O4a) of the LIGO-Virgo-KAGRA detectors. As in the other two papers in this series, we restrict our analysis to the 42 confident signals, measured by at least two detectors, that have FAR < 10^{-3}/yr from O4a, in addition to the 49 such events from previous observing runs. This paper focuses on the 8 tests that constrain parameterized deviations from the expected GR (or unaccelerated) values. These include modifications of post-Newtonian (PN) parameters, spin-induced quadrupole moments different from those of a binary black hole (BH), and possible dispersive or birefringent propagation effects. Overall, we find no evidence for physics beyond GR, for spin-induced quadrupole moments different from those of a Kerr BH in GR, or for LOSA, with more than 90% of the events including the null result (no deviation) within their 90% credible intervals. We discuss possible systematics affecting the other events and tests, even though they are statistically not surprising, given noise. The increased number of events analyzed allow us to improve the constraints on deviations from GR. For instance, for the PN coefficients, we improve the constraints by factors of 1.2-5.5, though some of this improvement is due to allowing the PN coefficient deviations to affect more of the waveform. We also provide illustrative translations to some modified theories. We update the bound on the graviton mass, at 90% credibility, to $m_g\leq1.92\times10^{-23}\mathrm{eV}/c^2$. Many of the bounds on possible deviations derived from our events are the best to date. △ Less

Submitted 20 July, 2026; v1 submitted 19 March, 2026; originally announced March 2026.

Comments: As part of the Astrophysical Journal Letters Focus Issue on the Gravitational Wave Transient Catalog; updated with ApJL resubmission and updated data release

Report number: LIGO-P2500066

GWTC-4.0: Tests of General Relativity. I. Overview and General Tests

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, D. Adhikari, N. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1759 additional authors not shown)

Abstract: The worldwide LIGO-Virgo-KAGRA network of gravitational-wave (GW) detectors continues to increase in sensitivity, thus increasing the quantity and quality of the detected GW signals from compact binary coalescences. These signals allow us to perform ever-more sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. This paper is the first of three, where we p… ▽ More The worldwide LIGO-Virgo-KAGRA network of gravitational-wave (GW) detectors continues to increase in sensitivity, thus increasing the quantity and quality of the detected GW signals from compact binary coalescences. These signals allow us to perform ever-more sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. This paper is the first of three, where we present the results of a suite of tests of GR using the binary signals included in the fourth GW Transient Catalog (GWTC-4.0), i.e., up to and including the first part of the fourth observing run of the detectors (O4a). We restrict our analysis to the 91 confident signals, henceforth called events, that were measured by at least two detectors, and have false alarm rates $\le 10^{-3} \mathrm{yr}^{-1}$. These include 42 events from O4a. This first paper presents an overview of the methods, selection of events and GR tests, and serves as a guidemap for all three papers. Here we focus on the four general tests of consistency, where we find no evidence for deviations from our models. Specifically, for all the events considered, we find consistency of the residuals with noise. The final mass and final spin as inferred from the low- and high-frequency parts of the waveform are consistent with each other. We also find no evidence for deviations from the GR predictions for the amplitudes of subdominant GW multipole moments, or for non-GR modes of polarization. We thus find that GR, without new physics beyond it, is still consistent with these GW events. The results of the two additional papers in this trio also find overall consistency with vacuum GR, with more than 90% of the events being consistent with GR at the 90% credible level. While one of the ringdown analyses finds the GR value in the tails for its combined results, this may be due in part to catalog variance. △ Less

Submitted 19 March, 2026; originally announced March 2026.

Comments: As part of the Astrophysical Journal Letters Focus Issue on the Gravitational Wave Transient Catalog

Report number: LIGO-P2500065

GRB 241030A: a bright afterglow challenging forward shock emission

Authors: J. -G. Ducoin, C. Pellouin, V. Aivazyan, D. Akl, F. Alvarez, C. Andrade, C. Angulo, S. Antier, J. -L. Atteia, S. Basa, R. L. Becerra, Z. Benkhaldoun, E. Bissaldi, A. Breeveld, E. de. Bruin, E. Burns, N. R. Butler, M. W. Coughlin, F. Daigne, T. Dietrich, D. Dornic, C. Douzet, T. du Laz, P. -A. Duverne, H. B. Eggenstein , et al. (75 additional authors not shown)

Abstract: Gamma-Ray Burst GRB 241030A (z = 1.411) exhibited a bright afterglow (similar to GRB 221009A), detected across gamma-ray, X-ray, UV, and optical bands, providing a probe of GRB afterglow physics. We compiled multi-wavelength observations spanning from a minute to a week after the prompt emission, processing the data through a unified photometry pipeline. We analysed the observations both analytica… ▽ More Gamma-Ray Burst GRB 241030A (z = 1.411) exhibited a bright afterglow (similar to GRB 221009A), detected across gamma-ray, X-ray, UV, and optical bands, providing a probe of GRB afterglow physics. We compiled multi-wavelength observations spanning from a minute to a week after the prompt emission, processing the data through a unified photometry pipeline. We analysed the observations both analytically and using Bayesian inference with two independent models. Our models assume that the afterglow emission arises from the strong forward shock of a laterally structured jet, with possible contributions from synchrotron self-Compton (SSC) scatterings. Our models reproduce X-ray to optical data, favouring a jet propagating into a constant-density interstellar medium, with a viewing angle within the jet core. However, both analyses require parameter values that are extreme compared to expectations from standard theory. In particular, our results imply extremely energetic jets despite regular prompt energy, leading to a very inefficient prompt emission. Furthermore, the jets are inefficient at accelerating particles, with low electron and magnetic energy fractions, leading to significant SSC emission. Our analyses indicate that the jets have large opening angles and propagate in high-density media. If the afterglow is indeed powered by radiation emitted behind a strong forward shock, our results place GRB 241030A within a sub-class of GRBs characterised by extreme kinetic energies, large jet opening angles, and very low prompt emission efficiencies, with strong SSC radiation. These predictions are difficult to reconcile with typical expectations from other GRBs. We therefore suggest that the afterglow of GRB 241030A is not solely powered by forward shock emission. △ Less

Submitted 19 March, 2026; originally announced March 2026.

All-sky Searches for Continuous Gravitational Waves from Isolated Neutron Stars in the Data from the First Part of the Fourth LIGO-Virgo-KAGRA Observing Run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, C. Adamcewicz, S. Adhicary, D. Adhikari, N. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith , et al. (1804 additional authors not shown)

Abstract: We present results from an all-sky search for continuous gravitational waves, using three different methods applied to the first eight months of LIGO data from the fourth LIGO-Virgo-KAGRA Collaboration s observing run. We aim at signals potentially emitted by rotating, non-axisymmetric isolated neutron star in the Milky Way. The analysis spans a frequency range from 20 Hz to 2000 Hz and accommodat… ▽ More We present results from an all-sky search for continuous gravitational waves, using three different methods applied to the first eight months of LIGO data from the fourth LIGO-Virgo-KAGRA Collaboration s observing run. We aim at signals potentially emitted by rotating, non-axisymmetric isolated neutron star in the Milky Way. The analysis spans a frequency range from 20 Hz to 2000 Hz and accommodates frequency derivative magnitudes up to $10^{-8}$ Hz/s. No statistically significant periodic gravitational wave signals were detected. We establish 95% confidence-level (CL) frequentist upper limits on the dimensionless strain amplitudes. The most stringent population-averaged strain upper limits reach 9.7 $\times$ $10^{-26}$ near 290 Hz, matching the best previous constraints from 250 to $\sim$1700 Hz while extending coverage to a much broader spin-down range. At higher frequencies, the new limits improve upon previous results by factors of approximately $\sim$1.6. These constraints are applied to three astrophysical scenarios: 1) the distribution of galactic neutron stars as a function of spin frequency and ellipticity; 2) the contribution of millisecond pulsars to the GeV excess near the galactic center; and 3) the possible dark matter fraction composed of nearby inspiraling primordial binary black holes with asteroid-scale masses. △ Less

Submitted 14 March, 2026; originally announced March 2026.

Comments: 45 pages, 17 figures

Report number: LIGO-P2500416

Nebular Phase Evolution of SN 2023ixf (I): From Circumstellar Infrared Echo to the onset of in-situ Dust Formation in a Type II Supernova

Authors: Avinash Singh, S. Goto, A. Sarangi, J. Johansson, C. Fransson, S. Barmentloo, J. Sollerman, R. S. Teja, K. Maeda, T. Hamada, N. Sarin, M. Yamanaka, T. Nakaoka, K. S. Kawabata, S. Schulze, A. Jerkstrand, S. Rose, D. K. Sahu, A. Gangopadhyay, G. C. Anupama, T. Ahumada, S. Anand, A. Bochenek, S. J. Brennan, X. Chen , et al. (37 additional authors not shown)

Abstract: We present optical and near-infrared (NIR) photometric and spectroscopic observations of the Type II supernova SN 2023ixf spanning 150 to 750 days, combined with published early-time optical and infrared photometry, and JWST NIRSpec and MIRI spectroscopy, to disentangle circumstellar echo emission from newly formed internal dust. The combined dataset reveals an early infrared excess by 1.8 days, a… ▽ More We present optical and near-infrared (NIR) photometric and spectroscopic observations of the Type II supernova SN 2023ixf spanning 150 to 750 days, combined with published early-time optical and infrared photometry, and JWST NIRSpec and MIRI spectroscopy, to disentangle circumstellar echo emission from newly formed internal dust. The combined dataset reveals an early infrared excess by 1.8 days, a broad secondary NIR rebrightening over about 89 to 175 days, progressive attenuation of the red wing of H-alpha from about 132 days, and CO emission detected by about 217 days. We identify the onset of H-alpha asymmetry as the first direct signature for internal dust formation, and modeling of the H-alpha profile over 140 to 418 days yields an internal silicate-equivalent dust mass of about 1.5e-6 to 6e-5 solar masses. By contrast, the early infrared evolution is best interpreted as echo-dominated: the 1.8 to 33.6 day excess is consistent with a radiative-flash infrared echo from pre-existing circumstellar dust, while the 89 to 175 day rebrightening is more naturally explained by a more extended echo arising from structured wind material. JWST spectral energy distribution modeling further reveals a multi-component infrared continuum in which a cold graphite component traces lingering echo emission, while a colder silicate-bearing component grows to about 2e-3 solar masses, providing the strongest late-time spectral energy distribution evidence that internal CDS/ejecta dust becomes substantial. SN 2023ixf therefore provides one of the clearest time-resolved case studies of dust signatures in a Type II supernova, linking early circumstellar reprocessing with increasingly important in situ dust formation. △ Less

Submitted 14 March, 2026; originally announced March 2026.

Comments: 63 Pages, 26 Figures, 7 Tables, Submitted to AAS Journals, Comments are welcome!

Long GRB 250916A: an Off-axis Powerlaw Jet with Thermal Cocoon

Authors: Utkarsh Pathak, Sameer K. Patil, Hitesh Tanenia, Tanishk Mohan, Xander J. Hall, Yogesh Wagh, Viswajeet Swain, Aditya Pawan Saikia, Varun Bhalerao, Tomas Ahumada, G. C. Anupama, Sudhanshu Barway, Malte Busmann, Michael W. Coughlin, Matthew J. Graham, Daniel Gruen, Assaf Horesh, Mansi M. Kasliwal, Russ R. Laher, Frank J. Masci, Antonella Palmese, Josiah Purdum, Argyro Sasli, Roger Smith, Xiaoxiong Zuo

Abstract: Some gamma-ray bursts (GRBs) exhibit precursor emission episodes preceding the main emission, with a quiescent period in between. The properties of the precursor emission and the duration of the quiescent interval are related to the central engine activity and jet formation processes, thus providing insights into the physics of GRBs. We present a comprehensive analysis of the prompt emission and… ▽ More Some gamma-ray bursts (GRBs) exhibit precursor emission episodes preceding the main emission, with a quiescent period in between. The properties of the precursor emission and the duration of the quiescent interval are related to the central engine activity and jet formation processes, thus providing insights into the physics of GRBs. We present a comprehensive analysis of the prompt emission and multi-wavelength afterglow of GRB 250916A. Using detailed afterglow modeling, we find that the broadband data are best described by a powerlaw structured jet with a relatively narrow core ($θ_c \approx 0.8^\circ$), viewed moderately off-axis at a viewing angle $θ_v \approx 2.7^\circ$. The isotropic-equivalent kinetic energy of the jet ($E_{k,iso} \approx 2.4 \times 10^{54}$ erg) is on the higher side for typical GRBs. The precursor emission is well described by a blackbody spectrum with a temperature of kT $\approx$ 13.2 keV and is separated from the main emission by a long quiescent interval of 150 s. Put together, our results indicate that the precursor is likely to be a shock breakout from a cocoon formed by the interaction of the relativistic jet with the progenitor star. The resulting cocoon pressure and shock collimation naturally lead to the launch of a narrowly collimated jet, consistent with the jet geometry inferred from afterglow observations. The long quiescent interval may imply the central engine turn-off in addition to the effect of the off-axis geometry. △ Less

Submitted 11 March, 2026; originally announced March 2026.

Comments: 16 pages, 6 figures, 4 Tables. Submitted to ApJ

Multi-epoch afterglow rebrightenings in GRB 250129A: Evidence for successive shock interactions

Authors: D. Akl, S. Antier, H. Koehn, P. T. H. Pang, J. J. Geng, R. Gill, E. Abdikamalov, C. Adami, V. Aivazyan, L. Almeida, S. Alshamsi, C. Andrade, Q. André, C. Angulo-Valdez, J. -L. Atteia, K. Barkaoui, S. Basa, R. L. Becerra, P. Bendjoya, D. Berdikhan, E. Bernaud, S. Boissier, S. Brunier, A. Y. Burdanov, N. R. Butler , et al. (83 additional authors not shown)

Abstract: Most long gamma-ray bursts (GRBs) exhibit afterglows broadly consistent with external forward-shock emission, typically described by smooth broken power-law decays in the multiband light curve. However, a minority of well-sampled GRBs deviate from this behavior, including GRB 250129A. This object shows multiple late-time rebrightenings at X-ray and optical wavelengths. Rebrightenings are often att… ▽ More Most long gamma-ray bursts (GRBs) exhibit afterglows broadly consistent with external forward-shock emission, typically described by smooth broken power-law decays in the multiband light curve. However, a minority of well-sampled GRBs deviate from this behavior, including GRB 250129A. This object shows multiple late-time rebrightenings at X-ray and optical wavelengths. Rebrightenings are often attributed to energy injection from prolonged central engine activity, refreshed shocks from delayed shell collisions, density jumps in the ambient medium, or angular jet structure and viewing-angle effects. After analysing the prompt emission observed in gamma-rays and the near-infrared, we investigate the origin of X-ray and optical flaring episodes in GRB 250129A. Physical processes in the afterglow light curves were investigated using methods ranging from empirical fitting to Bayesian inference. The well-sampled flares and the connection between the prompt and afterglow emission allow us to test the consistency of the fireball model and alternative scenarios. Conducting the prompt and time-resolved analyses, we obtained an isotropic-equivalent energy of E_iso,gamma = (1.35 +/- 0.12) x 10^53 erg. By modeling the afterglow using an agnostic Bayesian framework (NMMA), we rule out both a single external-shock evolution and a one-time energy-injection scenario. Numerical calculations show that the rebrightening episodes are consistent with refreshed shocks from delayed collisions between relativistic shells. Based on the consistency between our analyses of the prompt and afterglow GRB 250129A data, we find that two statistically significant rebrightening episodes occur within 1.1 days post trigger and can be explained by a sequence of refreshed shocks. Temporally and spectrally rich GRB datasets such as the one presented in this work, provide a powerful means to test current modeling frameworks. △ Less

Submitted 9 March, 2026; originally announced March 2026.

Journal ref: A&A 710, A211 (2026)

Beyond Gaussian Assumptions: A new robust statistical framework for gravitational-wave data analysis

Authors: Argyro Sasli, Minas Karamanis, Nikolaos Karnesis, Michael W. Coughlin, Vuk Mandic, Uroš Seljak, Nikolaos Stergioulas

Abstract: Many traditional algorithms applied in gravitational-wave astronomy rely on the assumption of Gaussian noise, a condition not always met. To meet this need, this study extends a robust statistical framework, advancing previous work on heavy-tailed likelihoods, that adapts the hyperbolic likelihood method for full frequency domain applications. The framework is designed to maintain high performance… ▽ More Many traditional algorithms applied in gravitational-wave astronomy rely on the assumption of Gaussian noise, a condition not always met. To meet this need, this study extends a robust statistical framework, advancing previous work on heavy-tailed likelihoods, that adapts the hyperbolic likelihood method for full frequency domain applications. The framework is designed to maintain high performance under ideal conditions while improving robustness against non-Gaussian noise and outliers in real-world data. We demonstrate the efficacy of this approach through two key case studies. The first case study analyzes a massive black hole binary merger in simulated Laser Interferometer Space Antenna (LISA) data with Gaussian noise, showing that the extended hyperbolic likelihood method performs comparably to the more commonly used Whittle likelihood. The second case study examines a stellar-mass black hole binary merger using real ground-based gravitational-wave data containing non-Gaussian noise or overlapping signals, where our framework exhibits increased robustness and yields more accurate parameter estimations. Our results show that the hyperbolic likelihood better captures the true noise distribution, providing a flexible and physically motivated alternative for GW data analysis across current and future detectors. △ Less

Submitted 25 February, 2026; originally announced February 2026.

Comments: 20 pages, 13 figures

AT2024lhc and AT2024kmq in the landscape of featureless tidal disruption events

Authors: Yuhan Yao, Ryan Chornock, Andrew Mummery, Raffaella Margutti, Marat Gilfanov, Muryel Guolo, Eric R. Coughlin, Wenbin Lu, Joheen Chakraborty, Dheeraj R. Pasham, Kate D. Alexander, Olivia Aspegren, Charlotte R. Angus, Xinze Guo, Xander J. Hall, Erica Hammerstein, K. -Ryan Hinds, Anna Y. Q. Ho, Xiaoshan Huang, Elias Kammoun, Natalie LeBaron, Matteo Lucchini, Zoë McGrath, Matt Nicholl, Daniel A. Perley , et al. (17 additional authors not shown)

Abstract: We study AT2024kmq and AT2024lhc, two tidal disruption events (TDEs) with blue featureless spectra associated with high-mass black holes ($M_{\rm BH}\sim 10^8\,M_\odot$). Both events show optical precursors consistent with shock dissipation from stream self-intersection. Their X-ray emission is luminous ($L_{\rm X}\sim 10^{44}\,{\rm erg\,s^{-1}}$), highly variable (with minimum observed variabilit… ▽ More We study AT2024kmq and AT2024lhc, two tidal disruption events (TDEs) with blue featureless spectra associated with high-mass black holes ($M_{\rm BH}\sim 10^8\,M_\odot$). Both events show optical precursors consistent with shock dissipation from stream self-intersection. Their X-ray emission is luminous ($L_{\rm X}\sim 10^{44}\,{\rm erg\,s^{-1}}$), highly variable (with minimum observed variability timescales of 1.3\,hr and 4.8\,hr for factor of $\sim3$ flux changes), long-lasting ($>1\,\rm yr$), emerging no later than the optical peak, and well characterized by power-laws with $1.7<Γ<3$ (where $f_ν\propto ν^{1-Γ}$). The X-ray properties and radio non-detections support a compact corona ($\lesssim 10 r_{\rm g}$) producing Comptonized X-ray emission. Using all published featureless TDEs, we find statistically significant bimodality in the distribution of their peak UV/optical blackbody luminosities and radii. We assemble a comparison TDE sample with early-time X-ray observations with eROSITA, in which we find different $M_{\rm BH}$ distributions in TDEs with different X-ray spectral evolution properties: low-mass black holes ($M_{\rm BH} \sim 10^6 M_\odot$) remain soft ($Γ>4$) within $t\lesssim 2$\,yr, intermediate masses ($\sim 10^7 M_\odot$) transition from soft to hard at $\sim$1 yr, while high masses ($\sim 10^8 M_\odot$) are hard ($1.5<Γ\lesssim 3$) from the outset. We interpret this result as evidence that the soft-to-hard state transition in TDEs occurs at the critical threshold of $\dot{M}_{\rm acc} \sim 0.03 \dot M_{\rm Edd}$ (similar to X-ray binaries), using the fact that the transition timescale predicted by simple disk theory scales with black hole mass as $t_{\rm tr}\propto M_{\rm BH}^{-3/4}$. △ Less

Submitted 25 February, 2026; originally announced February 2026.

Comments: 24 pages, 16 figures, submitted to MNRAS

Journal ref: Mon Not R Astron Soc (2026)

Multiwavelength Analysis of Six Luminous Fast Blue Optical Transients

Authors: Cassie Sevilla, Anna Y. Q. Ho, Nayana A. J., Steve Schulze, Daniel A. Perley, Michael Bremer, Igor Andreoni, Ivan Altunin, Thomas G. Brink, Michael Camilo, Poonam Chandra, Ping Chen, Ashley A. Chrimes, Michael W. Coughlin, Kaustav K. Das, Andrew Drake, Alexei V. Filippenko, Christoffer Fremling, James Freeburn, Avishay Gal Yam, Mary Gerhart, Matthew J. Graham, George Helou, K-Ryan Hinds, Natalya Johnson , et al. (21 additional authors not shown)

Abstract: We present multiwavelength observations and analysis of six luminous fast blue optical transients (LFBOTs) discovered in Zwicky Transient Facility (ZTF) survey data. We identified these LFBOTs from their fast light-curve evolution ($t_{1/2}\leq 12 $d), blue colors at peak brightness ($g-r\leq-0.5 $mag), a visible host galaxy, high optical luminosity ($M_g<-20$), and an X-ray or radio detection. Wi… ▽ More We present multiwavelength observations and analysis of six luminous fast blue optical transients (LFBOTs) discovered in Zwicky Transient Facility (ZTF) survey data. We identified these LFBOTs from their fast light-curve evolution ($t_{1/2}\leq 12 $d), blue colors at peak brightness ($g-r\leq-0.5 $mag), a visible host galaxy, high optical luminosity ($M_g<-20$), and an X-ray or radio detection. With the exception of AT2024aehp (ZTF24abygbss), these transients exhibit peaks in their $10\,$GHz radio light curves at $t_{\text{rest}} \approx 50-100$ d, with peak radio luminosities ranging from $10^{38}-10^{40}$ erg s$^{-1}$. Modeling the radio emission as synchrotron radiation indicates a fast ($v=0.1-0.3c$) shock in a dense ($n_e\approx10^{3}-10^{4}$ cm$^{-3}$) medium. The X-ray emission varies by $\approx2$ orders of magnitude in luminosity ($10^{42}-10^{44}$ erg s$^{-1}$) at $t_{\text{rest}}\sim20 $d. Analysis of the host-galaxy photometry and spectroscopy for each transient shows that they are predominantly nonnuclear (a few kpc offset) with star-forming host galaxies of stellar masses $10^{9}-10^{11} ,M_\odot$. Unlike all other LFBOTs to date, AT2024aehp exhibited a luminous ($M<-19 $mag) plateau in the optical light curve; spectra during this plateau phase showed a featureless blue continuum. The $6-15$ GHz radio emission of AT2024aehp brightened by over an order of magnitude from $t_{\text{rest}} \approx70 $d to $t_{\mathrm{rest}} \approx130 $d. The mostly consistent radio behavior between optically selected LFBOTs implies a similar circumburst medium, leading us to prefer a progenitor scenario in which mass is lost in a consistent way shortly prior to the terminal event, such as a massive star merging with a compact object. △ Less

Submitted 10 July, 2026; v1 submitted 26 January, 2026; originally announced January 2026.

Comments: 42 pages, 12 figures

Opportunities in AI/ML for the Rubin LSST Dark Energy Science Collaboration

Authors: LSST Dark Energy Science Collaboration, Eric Aubourg, Camille Avestruz, Matthew R. Becker, Biswajit Biswas, Rahul Biswas, Boris Bolliet, Adam S. Bolton, Clecio R. Bom, Raphaël Bonnet-Guerrini, Alexandre Boucaud, Jean-Eric Campagne, Chihway Chang, Aleksandra Ćiprijanović, Johann Cohen-Tanugi, Michael W. Coughlin, John Franklin Crenshaw, Juan C. Cuevas-Tello, Juan de Vicente, Seth W. Digel, Steven Dillmann, Mariano Javier de León Dominguez Romero, Alex Drlica-Wagner, Sydney Erickson, Alexander T. Gagliano , et al. (41 additional authors not shown)

Abstract: The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will produce unprecedented volumes of heterogeneous astronomical data (images, catalogs, and alerts) that challenge traditional analysis pipelines. The LSST Dark Energy Science Collaboration (DESC) aims to derive robust constraints on dark energy and dark matter from these data, requiring methods that are statistically powerful… ▽ More The Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) will produce unprecedented volumes of heterogeneous astronomical data (images, catalogs, and alerts) that challenge traditional analysis pipelines. The LSST Dark Energy Science Collaboration (DESC) aims to derive robust constraints on dark energy and dark matter from these data, requiring methods that are statistically powerful, scalable, and operationally reliable. Artificial intelligence and machine learning (AI/ML) are already embedded across DESC science workflows, from photometric redshifts and transient classification to weak lensing inference and cosmological simulations. Yet their utility for precision cosmology hinges on trustworthy uncertainty quantification, robustness to covariate shift and model misspecification, and reproducible integration within scientific pipelines. This white paper surveys the current landscape of AI/ML across DESC's primary cosmological probes and cross-cutting analyses, revealing that the same core methodologies and fundamental challenges recur across disparate science cases. Since progress on these cross-cutting challenges would benefit multiple probes simultaneously, we identify key methodological research priorities, including Bayesian inference at scale, physics-informed methods, validation frameworks, and active learning for discovery. With an eye on emerging techniques, we also explore the potential of the latest foundation model methodologies and LLM-driven agentic AI systems to reshape DESC workflows, provided their deployment is coupled with rigorous evaluation and governance. Finally, we discuss critical software, computing, data infrastructure, and human capital requirements for the successful deployment of these new methodologies, and consider associated risks and opportunities for broader coordination with external actors. △ Less

Submitted 20 January, 2026; originally announced January 2026.

Comments: 84 pages. This is v1.0 of the DESC's white paper on AI/ML, a collaboration document that is being made public but which is not planned for submission to a journal

Report number: FERMILAB-PUB-25-0886-CSAID-PPD