Peer Revenue Publish or profit

Papers

Fresh from arXiv, typeset properly. Highlighted words trade on the market.

Physics Astrophysics Cosmology and Nongalactic Astrophysics

astro-ph.CO Watch Cosmology and Nongalactic Astrophysics (sign in to use)

Thursday 24 September

Cosmology and Nongalactic Astrophysics

Joint measurements of thermal Sunyaev-Zel'dovich and cosmic infrared background cross-correlations with cosmic shear

The cross-correlation between weak lensing (WL) and the thermal Sunyaev-Zel'dovich (tSZ) effect probes the connection between the matter distribution and the thermal state of baryons at low redshift ($z\lesssim1$). A key systematic in such measurements is contamination by extragalactic foregrounds, particularly the cosmic infrared background (CIB). Here we describe the CIB spectral energy distribution (SED) with a small number of physically motivated spectral parameters, calibrated against external measurements of the radiation field intensity in galaxies. Their uncertainty is then propagated into the recovered WL$\times$tSZ signal in component separation at the level of the multi-frequency cross-power spectra rather than the map level. Applying this method to cosmic shear from the Dark Energy Survey and frequency maps from Planck, we show that the resulting measurements of the WL$\times$tSZ power spectrum are insensitive to variations in the effective CIB spectral index, radio source contamination, and the scale and redshift dependence of the SED. The measurements are robust to realistic foreground uncertainties, supporting their use as probes of baryonic feedback and cosmology, and providing a framework for future high-precision surveys. We compare our measurements with the Flamingo suite of hydrodynamical simulations and find that they are compatible with the Flamingo predictions for a low amplitude of matter fluctuations, parametrised by $S_8$, strongly disfavouring the larger $S_8$ preferred by Planck. This result is independent of the impact of baryonic feedback as modelled in Flamingo.
Read the full abstract

All authors Amy Wayland, Adrien La Posta, David Alonso, Baptiste Jego, Matthieu Béthermin

Cosmology and Nongalactic Astrophysics · Cross-listed from astro-ph.GA

Breaking the AGB Color Degeneracy with Variability and H$_2$O: A Consistent Picture of AGB Chemistry, Dust, and the J-Region Across Environments

Resolved asymptotic giant branch (AGB) populations present many observational puzzles. We consider three often treated separately: photometric C/O classification boundaries must be redrawn empirically in new environments; hydrostatic model grids are too blue in bands sampling H2O absorption; and J-region (JAGB) luminosities can differ between fields of a single galaxy. Pulsation connects them into a coherent picture. Using SPHEREx spectroscopy of AGB stars in the Large Magellanic Cloud (LMC) with JWST/NIRCam photometry and HST time-series of the SN Ia host M101, placed on a common scale by synthetic photometry, we investigate how AGB variability, surface chemistry, and circumstellar dust depend on environment. We develop an F182M index sensitive to 1.9$\mu$m H2O absorption that chemically classifies variable AGB stars and reveals oxygen-rich (O-rich) stars degenerate with carbon-rich (C-rich) stars in broadband colors. Because these features arise in extended atmospheric layers, hydrostatic grids struggle to reproduce them. Pulsation and circumstellar dust are asymmetrically linked: nearly all stars with substantial infrared excess are variable, while many variables show little dust. This supports mass loss initiated by pulsation and enhanced by dust formation. Variable O-rich AGB stars in M101 are up to ~0.7 mag redder than their LMC counterparts in dust-sensitive colors, consistent with more efficient silicate-dust production at higher metallicity, while C-rich populations depend less on environment. Color windows that select C-rich stars cleanly in the LMC are therefore 33-68% O-rich in the metal-rich M101. This challenges the J-region's universality as a standard candle while evading standard self-consistency checks. We recommend combining variability, or the infrared excess that implies it, with photometry of molecular features to improve AGB classification at 0.75-5.0$\mu$m.
Read the full abstract

All authors Caroline D. Huang, Massimo Marengo, Antonio Capodagli, Adam G. Riess, Louise Breuval, Stefano Casertano, Patricia A. Whitelock, Gagandeep S. Anand, Randall K. Smith, Wenlong Yuan

Cosmology and Nongalactic Astrophysics · Cross-listed from astro-ph.GA

SN Helios: A Multiply Imaged Type II Supernova Opening Time-Delay Cosmography beyond redshift of 3

We report the discovery and spectroscopic confirmation of SN Helios, a strongly lensed, multiply imaged Type II supernova (SN) at $z=3.34$ behind the galaxy cluster RXC J0018.5+1626. The transient was identified as a red F150W-dropout source (F277W = 26.1 mag) in JWST/NIRCam imaging taken on 2026 July 3 for the VENUS survey, with no counterpart in 2024 December imaging. Director's Discretionary Time observations on 2026 August 17 added second-epoch NIRCam imaging and a deep NIRSpec prism spectrum. The spectrum shows a strong, spectrally resolved H$\alpha$ emission line (FWHM $\approx 7900$ km/s) with a P-Cygni absorption blueshifted by 8800 km/s, together with corresponding H$\beta$ and He I $\lambda$10830 features, establishing a hydrogen-rich photosphere at $z=3.34$. Template fits classify SN Helios as a Type II SN, and light-curve comparisons with well-studied local SNe II yield a phase of $27.5\pm2.8$ rest-frame days past B-band maximum. Independent cluster lens models magnify the observed image by $\mu\approx 10$-20 and predict another image within $\simeq 2$-6 yr, offering a rare opportunity for a time-delay measurement of $H_{0}$ with a source at $z>3$. The host galaxy is undetected in deep rest-frame ultraviolet and optical imaging, in emission lines, and in ALMA 2 mm dust continuum, leaving a heavily obscured counterpart unlikely. Its de-lensed $M_{\rm UV}\gtrsim -13$ places it far below blank-field limits and adds support to an elevated core-collapse rate per unit star formation in ultra-faint galaxies. Without the multi-epoch data, SN Helios would masquerade as a convincing $z\approx 16$ Lyman-break galaxy, with galaxy-template fits strongly disfavoring any low-redshift solution even with multiple medium bands. Its spectrum thus provides an empirical SED template with which high-redshift searches can vet similar dropout candidates.
Read the full abstract

All authors Seiji Fujimoto, Conor Larison, Justin D. R. Pierel, Lukas J. Furtak, Masamune Oguri, Adi Zitrin, Jose M. Diego, Keren Sharon, David A. Coulter, Gabriel Brammer, Vasily Kokorev, Chelsea Nash, Pedram Abedi, Christa DeCoursey, Armin Rest, Andreas L. Faisst, Brenda L. Frye, Tiger Yu-Yang Hsiao, Kohei Inayoshi, Kotaro Kohno, Jorryt Matthee, Matteo Messa, Yoshiaki Ono, Fengwu Sun, Eros Vanzella, Rogier A. Windhorst, Abdurro'uf, Aadya Agrawal, Joseph F. V. Allingham, Yoshihisa Asada, Hakim Atek, Maruša Bradač, Larry D. Bradley, Mateusz Bronikowski, Dan Coe, Christopher J. Conselice, Pratika Dayal, Eiichi Egami, Michael Engesser, Gavin Scott Farley, Qinyue Fei, Ori D. Fox, Yoshinobu Fudamoto, Miriam Golubchik, Yuichi Harikane, Gourav Khullar, Tomokazu Kiyota, Paulo A. A. Lopes, Ray A. Lucas, Guillaume Mahler, Ashish Kumar Meena, Rohan P. Naidu, Gautham Narayan, Gaël Noirot, Pascal A. Oesch, Masami Ouchi, José María Palencia, Massimo Pascale, Luke Robbins, M. R. Siebert, Roberta Tripodi, Francesco Valentino, Darach Watson, Hayley Williams, Hiroto Yanagisawa, Erik Zackrisson, Anita Zanella

Cosmology and Nongalactic Astrophysics

Analyzing $\mu$ Parameter Range and the Duration of Reheating in the Polynomial $\alpha$-attractor Models using CMB Observations

Cosmological $\alpha$-attractor models are renowned for their universal predictions for inflationary observables, largely independent of the detailed form of the inflaton potential. In this work, we investigate the extent to which this universality persists for polynomial $\alpha$-attractor models, with particular emphasis on the scalar spectral index ($n_s$) and the post-inflationary reheating history. We derive an extended analytical expression for the scalar spectral index ($n_s$) and the tensor-to-scalar ratio ($r$), retaining sub-leading contributions beyond the commonly used universal approximation. Interestingly, we find that these corrections introduce an explicit dependence of $n_s$ on the parameter characterizing the polynomial attractor potential $\mu$, a dependence that is absent in the leading universal expression. We systematically investigate how this model dependence modifies the predictions in the $(n_s-r)$ plane and affects the model parameter ($\mu$). We further incorporate the Duration of reheating ($\Delta N_{rh}$) and effective equation of state parameter of the reheating epoch ($w_{rh}$), accounting for their impact on the number of e-folds between horizon crossing and the end of inflation ($\Delta N_{\mathrm{CMB}}$). By confronting the resulting predictions with current CMB observations, we identify the viable parameter space and examine the interplay between the model parameter ($\mu$) and reheating. Our results demonstrate that sub-leading corrections can provide a window into the underlying structure of polynomial attractors, offering a possible way to distinguish different members of the attractor class observationally.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from astro-ph.GA

TRICEPS. I. Tracing cold gas, dust, and stars at sub-kpc scales in massive star-forming galaxies at z = 4-5

We introduce the TRICEPS (Tracing Resolved Ionized Carbon Emission in Primeval Systems) survey. TRICEPS focuses on a sample of 16 massive galaxies at z = 4 - 5 with top-quality data: (i) ALMA band-7 observations of the [CII]-158 micron line, tracing the cold gas distribution and kinematics at spatial resolutions of about 0.1" (~0.8 kpc), (ii) ALMA maps of the rest-frame 160 micron continuum emission, tracing the dust distribution at similar spatial resolutions, (iii) multi-band JWST images from NIRCam and MIRI, tracing the stellar mass distribution and star formation activity at spatial resolutions of 0.1"-0.2", and (iv) multi-line CO observations from ALMA and the VLA, including the CO(1-0) line, tracing the molecular gas content. In this first paper, we describe the TRICEPS scientific goals, the basic properties of the galaxy sample, and the ALMA band-7 observations. TRICEPS galaxies were selected based only on bright [CII] fluxes (>2 Jy km/s). They lie on the top-end of the star-forming main sequence at z = 4-5, having stellar masses between about 10^10.5-10^12 Msun and star-formation rates between about 100 - 3000 Msun/yr. To date, TRICEPS is the largest sample of galaxies at z = 4-5 for which the gas distribution and kinematics are exceptionally well resolved: on average, 46 resolution elements per galaxy with signal-to-noise ratio higher than 3. By survey design, half of the TRICEPS targets are in galaxy pairs; clear signs of tidal interactions are present in three of the four pairs. In all galaxies, the [CII] velocity fields display regular rotation, albeit the interacting galaxies tend to be somewhat more disturbed. The TRICEPS observations suggest that regular rotation is ubiquitous among massive star-forming galaxies when the Universe was less than ~1.5 Gyr old.
Read the full abstract

All authors Federico Lelli, Lingrui Lin, Zhi-Yu Zhang, Cecilia Bacchini, Davide Bevacqua, Joe Bhangal, Marco Castellano, Carlos De Breuck, Giovanni Gandolfi, Allison W. S. Man, Antonino Marasco, Diego Paris, Borja Pérez-Díaz, Paola Santini

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-th

Tracing out massive fields in cosmology

Cosmological correlators are calculated through the in-in formalism in (quasi-)de Sitter background. Its effective description is best represented by a reduced density matrix. In this work, we study how integrating out a massive field in de Sitter space determines the reduced density matrix of a light field at a finite observation time. Starting from a quasi-single-field model with non-linear couplings, we first evaluate the massive field wavefunctional at fixed final value and then trace over this value. The trace generates branch-mixing terms which are generally present in the Schwinger–Keldysh description. For a conformally coupled massive field, we show explicitly that the light field bispectrum from the effective description agrees with a full theory calculation. This example also shows that decay of a field at the future boundary does not by itself justify discarding its contribution to the trace. We discuss the differences between an exact nonlocal effective description and the local effective field theory obtained through a large mass expansion. We then identify which local cubic terms affect the late-time probability and which contribute only a wavefunctional phase. Finally, we organize the resulting kernels by response, noise, and Schwinger–Keldysh consistency conditions, and discuss the symmetry constraints on the reduced density matrix.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

Hybrid Partial Dressing: Correct Effective Potentials at All Temperatures

We develop Hybrid Partial Dressing (HPD), a simple diagrammatic resummation scheme for finite-temperature effective potentials in scalar field theory which is valid at all temperatures. HPD improves on earlier resummation schemes such as Daisy Resummation and Partial Dressing by being two-loop exact at all temperatures and free of the overlapping momentum problem. Renormalization group improvement is straightforward to implement. We illustrate the scheme in multi-scalar theories with a spontaneously broken $\mathbb{Z}_2$ symmetry, and compare it against existing approaches, showing in particular that it reproduces the two-loop Dimensional Reduction (DR) potential exactly in the high-temperature regime where DR applies, while remaining valid at all temperatures. This makes HPD suitable for the general study of strong first order phase transitions beyond the Standard Model, as we demonstrate through gravitational wave predictions from HPD and its RG-improved version, RGHPD.
Read the full abstract

All authors Raphaël Berthiaume, David Curtin, Michael Luke, Andrija Rasovic, Jyotirmoy Roy

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

Metastable Neutron Stars as Transient-Rate Detectors of Heavy Dark Matter

Metastable neutron stars can undergo catastrophic conversion to quark-containing stars when a localized energy deposition nucleates a critical bubble. We show that dark-matter-induced conversions can be constrained by the present transient rate of the neutron star population, rather than by the survival of any individual hadronic star as often assumed in previous literature. Population depletion makes the present transition rate nonmonotonic: weak triggering produces few transitions, whereas sufficiently rapid triggering predominantly converts stars formed in the past. Consequently, transient observations can exclude a finite interval of interaction rates, qualitatively different from conventional one-sided survival bounds. Using the beaming-corrected short gamma ray burst rate as a benchmark, we derive constraints on dark matter annihilation, decay, inelastic particle scattering, and elastic scattering of macroscopic dark matter, including parameter space not covered by existing probes. Transient demographics thus provide a new population level probe of dark sector interactions and the phase structure of ultradense matter.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Statistics of large-scale structure from the transition probability of Lagrangian trajectories

We describe a probabilistic approach to gravitational clustering. Working purely in configuration space, we show that the two-point function can be formulated in terms of a transition probability describing the likely final position of a particle given its initial position. The transition probability is initially sharply peaked and then over time broadens out in a diffusion-like manner. The effect is more pronounced for widely-separated particles (e.g.\ at the baryon acoustic scale), leading to an effective suppression of the two-point function in what is normally referred to as infrared resummation.
Read the full abstract

Cosmology and Nongalactic Astrophysics

The Optical Expansion Scalar in Cosmological Perturbation Theory

Within the geometric optics approximation, the optical expansion scalar describes the local rate of change in the amplitude of a wave or the area of a light bundle propagating through the Universe. As such, it is fundamental for describing almost all the observations we make in cosmology. It was recently pointed out that this quantity can, in principle, be reconstructed from measurements of the angular diameter distance $D_A$ and the effective observed expansion rate $\mathfrak{H}$. Motivated by this possibility, we introduce the observable $\vartheta \equiv \theta/\omega_o$, where $\theta$ is Sachs' optical expansion scalar, and $\omega_o$ the observed photon frequency. We then derive the linear perturbative expression and angular power spectrum for $\vartheta$ in a flat FLRW universe with scalar perturbations. Comparing the power spectra with those of the fluctuations in $D_A$ and with power spectra from relativistic N-body simulations reveals no immediate striking advantage of studying the fluctuations of $\vartheta$ to those of $D_A$. However, the anisotropies of $\vartheta$ do exhibit several unique qualities such as an enhanced sensitivity to the peculiar velocity and relativistic potentials. We suggest that our results may hint at $\vartheta$ being useful for developing consistency tests of cosmological observations and possibly probing departures from FLRW cosmology and modified gravity theories. The formalism developed and presented here provides the foundation for future studies into these possibilities.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Two-rung ladder: $H_0$ from Tip of the Red Giant Branch and geometric anchors alone

We infer the Hubble constant $H_0$ from the ~400 Extragalactic Distance Database hosts with a tip of the red giant branch (TRGB) distance, calibrated on the geometric anchors of the Large Magellanic Cloud and NGC 4258. We forward-model the tip magnitudes and host redshifts including magnitude-limited selection, sky exposure, and density-dependent galaxy bias, conditioning on each of the 80 Manticore-Local realisations of the local density and peculiar-velocity fields in turn. As most hosts lie at $cz_{\rm CMB} \lesssim 1000$ km/s, where peculiar velocities are important, the inference is limited by the variance of these fields rather than by the TRGB calibration. For a fixed realisation, the redshift likelihood and density smoothing each move $H_0$ by less than 1 km/s/Mpc and a velocity monopole by $5.4 \pm 3.4$ km/s/Mpc, whereas the field-to-field scatter in $H_0$ is 5.2 km/s/Mpc. Weighting the posteriors from each realisation equally gives $H_0 = 68.4 \pm 5.9$ km/s/Mpc, while weighting each by its Bayesian evidence gives $72.3 \pm 2.6$ km/s/Mpc. The latter value is set by a single particularly high-evidence realisation, making it sensitive to the finite sampling of the reconstruction posterior. This ladder therefore cannot yet arbitrate the Hubble tension. Extending the TRGB data beyond the present ~1000 km/s, together with improved reconstructions of the local Universe, would reduce the field-to-field scatter well below the 5.2 km/s/Mpc we find here.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-th

Infrared Screening of the Cosmological Constant

The observed cosmological constant is set by the deepest infrared scale of the present universe rather than by the microscopic scales that contribute to the vacuum energy. We show that this emerges as the quantum fluctuations of the gravitational field are progressively taken into account. Considering the Einstein-Hilbert+$R^2$ truncation of gravity, we find that in a wide region of the parameter space the theory develops instabilities that give rise to large amplitude quantum fluctuations which govern the renormalization group flow. Once the flow enters this region, it remains there all the way towards the infrared and progressively loses memory of its boundary conditions. The running cosmological constant $\Lambda_k$ ($k$ is the running scale) is driven towards the universal infrared scaling $\Lambda_k\sim \lambda_{_{\rm IR}}\,k^2$, with $\lambda_{_{\rm IR}}= \frac{1}{2(\pi-2)}$. At the deepest infrared scale $k_{_{\rm IR}}\sim H_0$ (the current Hubble scale), the cosmological constant is then screened to its observed value $\Lambda_{_{\rm IR}}\sim H_0^2$.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

Plasma Effects Suppress Mixing-Induced Collisional Freeze-In

Dark matter that is weakly mixed with particles in a thermal plasma can be produced through both, oscillatory and collisional freeze-in. While plasma effects have long been known to suppress the former, their role in collisional freeze-in remains less clear quantitatively, leading to conflicting results in sterile neutrino dark matter production via collisional freeze-in. Using first-principles nonequilibrium Kadanoff-Baym equations in the flavour-covariant formulation, we establish the general result that plasma effects also suppress collisional freeze-in: the vacuum mixing angle in the collision rate is replaced by its plasma-corrected counterpart. Applying this result, we conclude that sterile neutrinos cannot be dark matter in the keV-MeV mass range when produced from the Standard Model plasma.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from astro-ph.GA

Discovery of an isolated RELHIC Candidate: J1351+0039

We report the discovery of J1351+0039, a new isolated dark galaxy candidate identified through a blind HI survey with the Five-hundred-meter Aperture Spherical radio Telescope (FAST), combined with deep optical imaging from the DESI Legacy Survey. The object displays a compact HI structure with a systemic velocity of V_{LSR} = 429.0 \kms and an inferred HI mass of 4.8*10^5 Msun. No optical counterpart is detected in the DESI images, with a g-band surface brightness limit of 28.1 mag arcsec^{-2}, corresponding to a stellar mass limit of M_{star} of 1.5*10^{5} Msun. The dynamical mass is estimated to be M_{dyn} of 1.3*10^{8} Msun, yielding a dark-to-baryonic mass ratio of 200, implying that J1351+0039 is overwhelmingly dominated by dark matter. While week rotation is detected in J1351+0039, the kinematics are dominated by internal velocity dispersion, confirming pressure support and reminiscent of Reionization-Limited HI Clouds (RELHICs). J1351+0039 stands out as one of the most promising isolated RELHIC candidates, offering a unique laboratory for studying baryonic physics and dark matter in the absence of star formation.
Read the full abstract

All authors Jin-Long Xu, Alejandro Benitez-Llambay, Ming Zhu, Chuan-Peng Zhang, Nai-Ping Yu, Xiao-Lan Liu, Mei Ai, Ruilei Zhou, Bin Liu

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

Impact of Heavy Modes on Primordial Black Hole Formation

Primordial black holes (PBHs) form from rare fluctuations in the high-density tail of the primordial distribution, making their abundance highly sensitive to primordial non-Gaussianity. Heavy fields during inflation can source such non-Gaussianity, producing characteristic signatures widely explored in cosmological collider physics, yet its impact on PBH formation remains largely unexamined. We investigate this question within quasi-single-field inflation (QSFI), the minimal realization that couples a heavy field to the inflaton. Unlike most PBH studies that consider only the bispectrum, the QSFI can generate a sizable trispectrum alongside the bispectrum, allowing us to investigate their combined effects on PBH formation. In our leading-order calculation, heavy-field non-Gaussianities can amplify the high-density tail and substantially enhance PBH production. Within the parameter scan, the largest enhancement occurs at the smallest heavy-field mass and the largest sound-speed ratio considered. The trispectrum contribution can be comparable to that of the bispectrum and may even dominate the total non-Gaussian correction. Our results establish PBH abundance as a novel observational channel for cosmological collider physics, demonstrating how heavy fields during inflation can significantly shape PBH formation.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from astro-ph.GA

Spatial segregation as the origin of anisotropic satellite quenching in haloes and beyond

Satellite galaxies in massive halos are seen to be preferentially quenched along the major axes of their central galaxies. In contrast, the AGN jet tends to point along the minor axis, which imposes interesting constraints on the role of AGN feedback in satellite galaxy quenching. Here we use the SIMBA, TNG100, EAGLE, and feedback-free SIMBA cosmological simulations at $z=0$ to show that this satellite anisotropy primarily originates from spatial segregation between quenched and unquenched populations. Quenched satellites are more centrally concentrated and more strongly aligned with the central galaxy's major axis, whereas unquenched satellites occupy larger radii and are more nearly isotropic. The persistence of this behaviour without stellar or active galactic nucleus feedback suggests that feedback is not important for satellite quenching. Beyond the halo boundary, the signal varies strongly among simulations and correlates with the excess number of galaxies along the major axis. This links large-scale anisotropic quenching to the direction-dependent distribution of groups, clusters and filamentary environments. Our results provide a unified geometrical interpretation of anisotropic satellite quenching across halo and extra-halo scales.
Read the full abstract

All authors Zhuoming Zhang, Weiguang Cui, Yun Chen, Romeel Davé

Cosmology and Nongalactic Astrophysics

Limits on Primordial Black Hole Evaporation from LUX-ZEPLIN

The LUX-ZEPLIN (LZ) collaboration recently reported a singular anomalous nuclear recoil event at $248 \pm 23 \, \rm keV$, accompanied by a strict null result for unexplained excesses in their primary $5.4-50 \, \rm keV$ search window. The macroscopic momentum transfer required to generate this event is kinematically inaccessible to standard halo cold dark matter (CDM), naturally motivating models featuring light, relativistic relics. In this work, we investigate whether a contemporary flux of Hawking-boosted dark matter (DM) emitted by evaporating Primordial Black Holes (PBHs) in the $10^{11}-10^{13} \, \rm g $ mass range can source the anomaly. By parameterizing the nuclear scattering with heavy non-relativistic effective field theory (NREFT) operators and endothermic inelastic mass transitions, we formulate a strict kinematic exclusion. We demonstrate that the inherent thermal nature of the Hawking emission, coupled with the relativistic kinematics of the incident flux, inevitably overproduces low-energy recoils, aggressively violating the LZ background bounds. Utilizing this failure, we map the LZ low-energy null results into novel, highly stringent upper limits on the PBH abundance fraction, excluding $f_{\rm PBH} \gtrsim 10^{-6}$ for $M_{\rm PBH} \sim 2 \times 10^{11} \, \rm g $ at $ 90 \% $ confidence level.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from gr-qc

Joint Model-independent Constraints on the Post-Newtonian Parameter and Cosmic Curvature from Galaxy-scale Strong Lensing, SNe Ia, BAO, and Cosmic Chronometers

This investigation aims to test the validity of general relativity (GR) on kiloparsec scales via galaxy-scale strong gravitational lensing (SGL) combined with multiple cosmological probes. To circumvent the circularity problem induced by the presumption of a cosmological model based on GR, we apply the distance sum rule in the Friedmann-Lemaitre-Robertson-Walker (FLRW) metric to constrain the post-Newtonian parameter gamma PPN and spatial curvature Omega k independently of any cosmological model. The study introduces two methodological improvements: (i) parameterizing the dark-energy evolution factor with second-order Chebyshev orthogonal polynomials, which ensures uniformly high convergence across the full redshift baseline while retaining Omega k as an explicit free parameter; (ii) extending the original two-probe configuration (SGL + Type Ia supernovae) to a four-probe joint analysis that includes baryon acoustic oscillations and cosmic chronometers. Our results show that gamma PPN = 1.119 (+0.072/-0.072) and Omega k = 0.101 (+0.077/-0.076) at 68% confidence. The gamma PPN value is consistent with the GR prediction of unity at 1.7 sigma, while Omega k favors a mildly open universe, with flatness consistent at 2 sigma. Taking model-independent constraints as a benchmark, we quantify the systematic shift in gamma PPN induced by eight representative dark energy models, and find that all offsets are negligible, well below the baseline 1 sigma uncertainty, with no significant bias at current precision. These results confirm that conventional strong-lensing gravity tests introduce no measurable systematic bias, validating our model-independent framework.
Read the full abstract

All authors Fu-you Zheng, Xiao-jun Gao, Hou-bing Tang, Wei Xiang

Cosmology and Nongalactic Astrophysics

Testing cosmic acceleration from thermogravity without vacuum energy

We present a first background-level observational test of a thermogravity theory in which the Einstein equations are clipped to a trace-free version, the cosmological constant $\Lambda$ does not gravitate, and controlled violations of energy conservation obstruct the usual reinstatement of $\Lambda$ as an integration constant. We therefore set $\Lambda=0$ and ask whether late-time acceleration can instead be generated solely by energy non-conservation, with no extra background parameter relative to flat $\Lambda\mathrm{CDM}$. We contrast a minimal universal implementation with a model in which only CDM partakes in non-conservation. The universal model is strongly disfavored when combining Supernovae and BAO distances, with the inclusion of DESI BAO worsening the fit by $\Delta\chi^2\simeq 51$ relative to $\Lambda\mathrm{CDM}$, because matter creation ties the intermediate-redshift normalization of $H(z)$ too rigidly to the present acceleration. This conclusion, however, should be interpreted with caution, since universal non-conservation would modify the observational dictionary itself. By contrast, restricting non-conservation to CDM leaves baryonic and photon observables unaffected at the background level and therefore allows a self-contained analysis. The resulting $\xi\mathrm{CDM}$ model, which can be thought of as a one parameter extension of $\Lambda\mathrm{CDM}$, provides an improved fit relative to $\Lambda\mathrm{CDM}$ for the background dataset combinations considered, with improvements reaching $\Delta\chi^2=-4.85$ and a maximum Bayesian preference of $\ln\mathcal{B}=2.77$. However, when the BAO ruler is calibrated using BBN or CMB information, the reduced effective early-time CDM density increases the sound horizon and drives $H_0$ towards lower values, thereby increasing rather than alleviating the tension with the distance-ladder calibration.
Read the full abstract

All authors Dong Ha Lee, João Magueijo, Carsten van de Bruck, Eleonora Di Valentino

Cosmology and Nongalactic Astrophysics

QCD Crossover Transfer Functions for Scalar-Induced Gravitational Waves in the PTA Band

Pulsar timing array (PTA) collaborations have reported evidence for a stochastic gravitational wave (GW) background in the nHz band. Should scalar-induced GWs, sourced at second order by enhanced primordial curvature perturbations, contribute to this signal, a coincidence of scales makes them directly sensitive to the softening of the equation of state around the QCD crossover. We solve the tensor and scalar equations of motion across the Standard Model (SM) thermal history and tabulate the transfer functions for direct use in present and future PTA analyses. We show that this SM effect modifies the height of the induced spectrum by up to $\approx55\%$ across the PTA band relative to the radiation-domination expectation, with either sign depending on whether the source modes cross the horizon before or after the crossover. Fitting the NANOGrav 15-year data with a broken-power-law curvature power spectrum, we find that including the crossover shifts the inferred peak amplitude and scale by an amount that could already be relevant for the comparison with primordial black hole overproduction bounds. The importance of this SM effect will grow as the statistical uncertainty on the amplitude and scale shrinks with future, more sensitive PTA datasets, at which point neglecting it could significantly bias the inference.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from gr-qc

Matter density perturbations in Quadratic Gravity

We investigate linear matter-density perturbations in full quadratic gravity. Working in the longitudinal gauge and adopting the sub-horizon and quasi-static approximations, we obtain an algebraic expression for the effective gravitational coupling in the theory and the equation for evolution of matter density perturbations. We identify the positive-real poles of this coupling and map their dependence on the quadratic couplings. In the formal infrared continuation, the reduced coupling recovers its general-relativistic value, whereas at large physical wavenumber it is suppressed when both quadratic couplings are nonzero; the single-coupling limits reproduce distinct behaviors. Numerical integrations of the growing mode are performed up to the first pole and show that pole-free regions of the reduced system can support regular matter-growth solutions, while poles obstruct a continuous late-time evolution. Our results show that quadratic gravity cannot generally reproduce the standard cosmological evolution of perturbations and in the regions where it can, it closely matches the behavior of $f(R)$ theories.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from astro-ph.HE

DETECT: Real-Time Identification of Transients with DESI Spectroscopic Redshift

Wide-field surveys now report thousands of transients per month, while spectra are obtained for only a few per cent of them. We present DETECT (DESI–Transient Event Cross-matching Tool), a pipeline that turns each Transient Name Server (TNS) alert into a distance-informed candidate within an hour by placing it in the context of the Dark Energy Spectroscopic Instrument (DESI) archive. DETECT cross-matches every new report against the 22 million spectra of DESI EDR and DR1 through a HEALPix index, associates the transient with a host galaxy by a directional-light-radius rule whose morphology-dependent thresholds are calibrated on 4,474 supernovae with known redshifts (92.5% completeness, 0.4% wrong hosts), and passes ambiguous cases to a reviewer through a web interface. With a validated host redshift it derives the peak absolute magnitude, the projected offset and basic host properties of each event, and ranks it for spectroscopic follow-up. Applied retrospectively to ~1.2x10^5 TNS reports from 2020–2024, DETECT recovers a spectroscopic host for 15% of them, from these we build a Gold Sample of ~5,400 well-sampled light curves whose luminosity distributions by class reproduce those of untargeted surveys. Because the absolute magnitude is known at discovery, intrinsically overluminous events stand out immediately: in the prospective 2025 run this is how the lensed superluminous supernova SN~2025wny was flagged within hours of its report, and how faint, fast-declining kilonova candidates were screened against model grids and archival photometry. DETECT shows how an archival spectroscopic survey can make host redshifts a routine part of transient triage ahead of the Rubin Observatory's LSST.
Read the full abstract

All authors Yu-Hsing Lee, Ting-Wan Chen, Ze-Ning Wang, Sheng Yang, Limeng Deng, Chuan-Jui Li, K. C. Chambers, Thomas de Boer, Chien-Cheng Lin, Thomas B. Lowe, Paloma Mínguez, Gregory S. H. Paek, Richard Wainscoat

Wednesday 23 September

Cosmology and Nongalactic Astrophysics · Cross-listed from astro-ph.IM

FPCA-Enhanced Simulation-Based Inference for Robust Type Ia Supernova Cosmology

Precision cosmology is a cornerstone of modern astronomy, and upcoming wide-field photometric surveys make it crucial to develop robust, data-driven inference methods that can handle complex survey systematics without relying on closed-form likelihoods. Simulation-Based Inference (SBI) meets this need by enabling forward simulations that encode complex survey characteristics. Previous SBI analyses in supernova cosmology have used SALT2 light-curve parameters as summary statistics; however, SALT2's two-basis spectral template imposes rigid modeling assumptions, motivating a more flexible representation. In this work, we present the first application of Functional Principal Component Analysis (FPCA) light-curve parameters as summary statistics for SBI-based cosmological inference under a non-flat {\Lambda}CDM model. On identically generated simulations, FPCA+SBI yields constraints comparable to both SALT2-based SBI and explicit-likelihood analyses, while providing more robust constraints on out-of-domain simulations than SALT2. Applying the model trained on LSST-like light curves to a spectroscopically confirmed DES Year 5 supernova sample, we recover constraints consistent with those of the DES collaboration to within 0.12 {\sigma} and 0.23 {\sigma} in {\Omega}m and {\Omega}{\Lambda} respectively, establishing the method's generalizability to real survey data. By introducing host-dependent systematics through a mass-dependent dust extinction law, we find FPCA-based summary statistics implicitly encode host-dependent systematics, suggesting dedicated host modeling may be unnecessary within this framework. Combined with its previously demonstrated effectiveness for photometric classification, FPCA offers a compelling foundation for unified, data-driven pipelines capable of jointly performing supernova classification and cosmological inference for upcoming wide-field photometric surveys.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

A Stepping Stone Solution to the QCD Axion Isocurvature Problem

We propose a ”stepping stone” realization of the varying Peccei Quinn (PQ) scale solution to the QCD axion isocurvature problem. In conventional implementations, suppressing isocurvature during high scale inflation requires a large inflationary PQ scale, while the subsequent evolution toward the much smaller late time vacuum can induce parametric resonance and nonthermal restoration of the PQ symmetry. We ameliorate this tension by introducing an intermediate PQ-breaking vacuum that is reached before the end of inflation, so that the later evolution to the present day PQ scale spans only a moderate hierarchy. The vacuum structure arises radiatively: a Yukawa coupling between the PQ scalar and PQ-charged fermions drives the quartic coupling negative at high scales while keeping it positive at lower scales. A positive dimension six operator then stabilizes a large field minimum, while the positive low energy quartic together with the quadratic terms generates the intermediate vacuum. We study two mechanisms in which the rolling inflaton removes the large field minimum, either by scanning a hidden non-Abelian gauge coupling or by directly scanning the Yukawa coupling. In both cases the PQ field is released toward the intermediate vacuum during inflation, after CMB modes have exited the horizon. We exhibit benchmark regions with subdominant backreaction on the inflationary sector and suppressed vacuum tunneling, and show that the construction permits substantially less hierarchical quartic and higher-dimensional couplings than conventional realizations that depend on a varying PQ scale. The underlying logic may have analogues for extra-dimensional or string axions, where the effective axion decay constant is controlled by geometric moduli, although realizing the required intermediate vacuum in a controlled modulus potential is nontrivial.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from gr-qc

Unpolarized Low-Frequency Tail of Causal Chiral Gravitational Waves

Parity violation in the early Universe may lead to a chiral gravitational wave background. We show that, for finite-time parity-violating sources, the low-frequency tail of the generated chiral gravitational wave background is universally unpolarized, and that the chirality fraction decreases linearly with frequency. Our result generally applies to any parity-violating gravitational-wave source that appears quadratically in the wave equation. We consider both Gaussian sources and explicit trispectrum templates. We attribute the lack of chirality of the low-frequency tail to the fact that, on large scales, gravitational wave generation is almost equally inefficient for both helicities.
Read the full abstract

All authors Tikinas Chabour, Guillem Domènech, Alexander Ganz, Jan Tränkle

Cosmology and Nongalactic Astrophysics · Cross-listed from astro-ph.IM

PSFSim: PhySics First Simulations for the Point Spread Function of the Roman Space Telescope

We present a new tool for generating simulated point spread functions (PSFs) in the Wide Field Instrument (WFI) of the Nancy Grace Roman Space Telescope. The simulation tool is physics forward, and performs a raytrace of the optical system of the instrument to compute the PSF, while also accounting for various detector effects. The tool also has the capability to add perturbations to the optical system and detectors. We explore several applications of the tool that may arise during the commissioning of the telescope, including optical ghosts, measuring the optical aberrations from diffraction spikes, and studying the effect of mirror polarization.
Read the full abstract

All authors Nihar Dalal, Chris Hirata, Charuhas Shiveshwarkar, Anthony Harbo Torres, Chun-Hao To, Elise Moore, Katherine Laliotis, David Kuhtenia, Emily Macbeth

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-ph

Heavy Higgsino Dark Matter at the Scale of the Vanishing Higgs Quartic

The Standard Model Higgs quartic coupling runs to zero around $10^9-10^{12}$ GeV at $2\sigma$, suggestive of an ultraviolet boundary condition $\lambda\simeq0$ at this scale. In a well-motivated class of high-scale supersymmetric (SUSY) models, the quartic is automatically set to zero at tree level at the superpartner scale. The Higgsino can naturally emerge as the lightest SUSY particle, with its two neutral Weyl components sufficiently degenerate to form an effectively Dirac fermion with an unsuppressed vector coupling to the $Z$. For such Higgsino dark matter, direct detection requires a mass above $\sim10^{11}$ GeV, leading to a coincidence with the scale independently selected by the vanishing quartic. Such a heavy Higgsino cannot freeze out with the observed relic density. However, if its mass exceeds the reheating temperature, the observed dark matter abundance can be reproduced through Boltzmann suppressed freeze-in. This scenario connects the dark matter mass to the measured Higgs mass, with a parameter space already constrained by LUX-ZEPLIN (LZ) and which will be decisively tested at DARWIN/XLZD. Finally, we consider PeV Higgsinos with Majorana gauginos of order $10^7$ GeV, permitting an inelastic Higgsino interpretation of the tentative LZ event while retaining the connection between the scalar superpartner scale and the vanishing Higgs quartic scale.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from hep-th

Breaking de Sitter Symmetry with Cosmic Topology: Cosmological Role of Infrared Vacuum Freedom

When spatial sections of spacetime are topologically nontrivial, the symmetry arguments usually used to choose the Bunch-Davies vacuum are no longer guaranteed to apply. We study this problem for a scalar field on de Sitter spacetime with toroidal spatial sections. The Hadamard conditions, which ensure the same short-distance structure as in flat space, are local and naturally insensitive to the compactification. By contrast, we find that the infrared sector is affected by the toroidal identifications that break global de Sitter invariance, thereby removing the usual argument for discarding zero-mode contributions. We treat the remaining zero-mode freedom as an effective topology-induced source and analyze its effects on an inflationary de Sitter expansion within a semiclassical setting. We show that, for a massless scalar, the zero-mode sector can contain a contribution to the stress-energy tensor with equation of state $P=-\rho$, behaving as a cosmological constant. Thus, the source driving an early-Universe de Sitter phase need not be determined solely by local terms in the action, but can also depend on the infrared structure of the quantum state. Nontrivial topology, as a concrete realization of de Sitter symmetry violation, turns the otherwise-unfixed infrared freedom in vacuum selection to a contribution directly relevant to the semiclassical expansion. This contribution could potentially be applicable to the current accelerated expansion of the Universe.
Read the full abstract

All authors Anna Negro, Kurt Hinterbichler, Glenn D. Starkman, Yashar Akrami, Stefano Anselmi, Javier Carrón Duque, Mikel Martin Barandiaran, Andrius Tamosiunas

Cosmology and Nongalactic Astrophysics

Wave-Optics Imprints of Warm Dark Matter Subhalos with Prompt Cusps on Strongly Lensed Gravitational Waves

Dark matter halos are expected to form with a prompt $\rho \propto r^{-3/2}$ density cusp at their centres, and in warm dark matter (WDM) cosmologies these cusps can dominate the inner structure of the low-mass subhalos that survive free-streaming suppression. We investigate whether that inner structure is visible to gravitational waves lensed in the wave-optics (WO) regime. Extending the strong lensing WO diffraction-integral framework to WDM, we generate subhalo populations with the sashimi-w semi-analytic model, assign each subhalo a prompt cusp through a cusp-halo relation, and compute the frequency-dependent amplification factor $F(f)$ across the LISA band over 500 independent realisations for $m_\chi \in \{6, 10, 20, 40, 80, 100\}$~keV. The signal is governed by the overlap between the WDM-suppressed subhalo mass function and the $m_{\rm sub} \sim 10^{4}$ to $10^{7}\,M_\odot$ range to which the LISA band is sensitive. Modulations are negligible for $m_\chi \lesssim 10$~keV, reach the percent level by 40~keV, and then saturate to a CDM-like plateau, so the observable reliably probes $m_\chi$ only over the intermediate range $\sim 20$ to 40~keV. Repeating the 40~keV ensemble on identical subhalo catalogs with the cusp amplitude set to zero, we find that prompt cusps enhance the amplitude modulation only mildly. Where cusps make up a large fraction of the subhalo mass, free-streaming has already depleted the signal-carrying mass range, and where that range is populated, the cusp lies well inside the Fresnel radius, to which the diffraction integral responds through projected enclosed mass rather than central density. A percent-level detection in a strongly lensed massive black-hole binary would therefore establish the presence of substructure on Fresnel scales without diagnosing its inner profile.
Read the full abstract

Cosmology and Nongalactic Astrophysics · Cross-listed from astro-ph.IM

Differentiable astrophysics at scale: solving and differentiating ODE ensembles on the GPU

Astronomers increasingly fit their models with gradient-based methods, such as Hamiltonian Monte Carlo, which need the derivatives of the model with respect to its parameters. In many analyses a prediction requires solving a small system of ordinary differential equations (ODEs) for thousands to millions of parameter sets, and this ensemble of integrations often sets the cost of the analysis. We introduce to the astronomical community GRADSOLVE, a JAX library that moves this computation to graphics processing units (GPUs): it integrates each member of the ensemble in its own GPU thread with its own step size and returns the derivatives with respect to the parameters in the same pass. We measure the speed-up it provides in three examples from different fields of astronomy, stellar orbits in the Galactic potential, the expansion history of a dark-energy cosmology and the spin precession of binary black holes, with every code held to the same accuracy requirement. For a million trajectories GRADSOLVE runs 8.5 to 1500 times faster than the serial CPU code of each example running on all 128 cores of a CPU, and several orders of magnitude faster than on one core. On the same GPU it is also 11 to 15 times faster than DIFFRAX, the state-of-the-art ODE library in JAX, in all three examples. With GRADSOLVE a million such integrations take seconds or less on one GPU, so gradient-based analyses that need ensembles of this size become routine. The code is publicly available at https://github.com/ECLIPSE-AI4Science/gradsolve.
Read the full abstract