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)

Tuesday 8 September

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

Detection prospects for heavy WIMP dark matter around M31* in microwave band

This work analyzes the detection prospects for weakly interacting massive particles (WIMPs) in dark matter (DM) density spike around the supermassive black hole (SMBH) in Andromeda galaxy M31 in microwave band. WIMP annihilation produces relativistic electrons and positrons, which generate the synchrotron emission in strong magnetic field of the accretion flow medium at microwave frequencies. Properties of this emission, as well as their dependence on DM and medium parameters, were modeled. Atacama Large Millimeter Array (ALMA) sensitivity to this emission and underlying WIMP parameters was estimated. It was obtained that M31* in microwave band represents the unique target, where DM spike signal can be detected in case of realistic spike and halo density profiles. In such case ALMA may reach the thermal s-wave annihilating WIMPs with masses up to around 1 TeV with just 10–100 hours exposure. Optimistic parameter configurations may allow probing all possible WIMP masses up to 100 TeV. Related systematic uncertainties are discussed.
Read the full abstract

Cosmology and Nongalactic Astrophysics

$\texttt{LensFactory.jl}$: A general-purpose strong lens modeling package

We present LensFactory.jl, an efficient, open-source, general-purpose strong lens modeling package written in Julia. It supports single- and multi-plane lensing with a wide range of analytic lens profiles, and its modular design, built on Julia's multiple dispatch, makes it straightforward to add new profiles or inference methods, while native multi-threading parallelizes the computationally expensive steps, enabling full cluster-scale lens reconstructions within hours on a desktop-class machine. We validate the mass-reconstruction pipeline of LensFactory.jl on two simulated galaxy clusters (Ares and Hera), recovering the convergence and magnification maps well where the lensing constraints lie, and on the JWST-observed cluster SMACS J0723.3-7327, obtaining a strong lens model (image-plane RMS of 0.46”) consistent with previously published results. In addition, we introduce singularity maps as a complementary diagnostic for comparing lens models, probing differences in the higher-order derivatives of the lensing potential, which remain largely unconstrained by standard observables such as image positions, time delays, and flux ratios. LensFactory.jl is accompanied by a set of worked examples covering all core functionality.
Read the full abstract

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

A Missing Tool for Calculating Auto/Cross-correlation Function under Nonuniform Sampling Observations

Nonuniform sampling presents a long-standing challenge in astrophysical time-domain analysis, invalidating the standard autocorrelation and cross-correlation functions and forcing researchers to adopt ad-hoc methods like interpolation or binning, which introduce unquantified biases and lack rigorous error estimation. Here we introduce a new method for calculating the nonuniform autocorrelation function (NUACF) and nonuniform cross-correlation function (NUCCF) for irregularly sampled time series. Instead of relying on interpolation, it naturally evaluates the correlation function by incorporating time-interval weights and misalignment penalties. Monte Carlo simulations provide confidence bands for significance assessment and a complete error budget for the time delays that accounts for both flux uncertainties and sampling irregularity (essential but generally absent from existing methods). Through extensive simulations, we demonstrate our method outperforms traditional methods across various conditions, from strictly periodic to complex repeating variability patterns (e.g., intermittent but aperiodic). Its effectiveness is demonstrated via various real astrophysical data sets, revealing repetitive variability in stellar light curves, measuring time delays for multi-band disc reverberation in the AGN Fairall 9, and providing model-independent validation of time delays for the gravitationally lensed quasar HE 0435-1223. The method provides a rigorous and general solution to the ubiquitous problem of nonuniform sampling, positioning it as a useful tool for large-scale time-domain survey data analysis. The framework is also directly applicable to emerging time-domain phenomena such as fast radio bursts (FRBs), enabling, e.g., the study of correlations between persistent radio source luminosity and repeating FRB activity, or among the multi-parameter variability curves of FRB emission itself.
Read the full abstract

All authors Chen-Ran Hu, Yong-Feng Huang, Jin-Jun Geng, Orkash Amat, Ze-Cheng Zou, Chen Deng, Fan Xu, Xiao-Fei Dong, Chen Du, Nurimangul Nurmamat, Pei Wang, Lang Cui, Cheng-Ming Li

Cosmology and Nongalactic Astrophysics

Harnessing stellar kinematics to constrain dark energy with the double-source-plane gravitational lens SDSS J0946+1006

SDSS J0946+1006 is an attractive target for measuring cosmological parameters through lens modelling. It is the best-studied galaxy-scale strong gravitational lens with multiple sources whose redshift separations are well-suited to constraining the dark energy equation of state. However, multi-plane lens models with free cosmological parameters risk a multi-plane mass-sheet degeneracy, although this can be lifted by a non-lensing deflector density profile tracer. We simultaneously reconstruct near-infrared and near-ultraviolet HST imaging whilst including a velocity dispersion measurement from VLT-MUSE to constrain the foreground deflector. Imaging is reconstructed faithfully regardless of whether the inferred kinematics are realistic, though we find the kinematic constraint essential for shifting the preferred cosmology into a region not in significant tension with other dark energy probes. Deflector density profile perturbations, via substructure and multipolar halo shape deformations, have only a modest effect on inferred cosmology. Combining our fiducial model with Planck CMB data yields $w=-1.01^{+0.08}_{-0.13}$; or, combined with Pantheon SNe Ia, $w=-0.99^{+0.13}_{-0.15}$. We further show that this system proves a remarkably valuable complementary probe in the $w_{0}$-$w_{a}$ plane of an evolving dark energy model, and that the DESI BAO tension with $\Lambda$CDM seen when combined with other datasets is not reproduced when combined with this lens, yielding $(w_{0}, w_{a})=(-0.87^{+0.10}_{-0.11}, -0.22^{+0.28}_{-0.25})$. The system's third source, visible with MUSE, only weakly constrains $w$CDM but may strengthen $w_{0}w_{a}$CDM constraints, though further mass-model complexity along its line of sight is required. Overall, kinematics-informed multi-plane lens modelling is a robust route to competitive dark energy constraints, even with a single system.
Read the full abstract

All authors Daniel J. Ballard, Wolfgang J. R. Enzi, Thomas E. Collett, Coleman M. Krawczyk, Giovanni Granata, Tian Li, Geraint F. Lewis, Karl Glazebrook, Huimin Qu

Cosmology and Nongalactic Astrophysics

Vacuum Decay and Baryogenesis Associated with Primordial Black Holes at Finite Temperature

Within finite-temperature effective potential framework, we study the effects of PBHs on vacuum decay, the associated gravitational wave spectrum, and the baryon asymmetry. For the PBHs considered, the Hawking temperatures are below the ambient temperature, and their dominant effect is taken to be the gravitational distortion of the Higgs field configuration. Relative to flat spacetime, the PBH background reduces the vacuum decay exponent, enhances the gravitational wave peak amplitude, and shifts the peak frequency downward. For PBH masses of roughly $10^{11}~{\rm g}-10^{13}~{\rm g}$, the sphaleron energy is reduced, while the small PBH abundance renders the associated washout negligible. Additionally, over roughly the same mass range, PBH-emitted BSM particle decays generate a baryon asymmetry too small to spoil the agreement between the baryon abundances inferred from BBN and CMB.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Finding the distribution of matter using lenses - II: deconvolution-based reconstruction with 3x2pt measurements

We present a deconvolution-based framework for testing scale-dependent departures of the late-time matter power spectrum from a fiducial cosmological model using $3\times2$pt measurements. We introduce a free-form scale-dependent modulation $A(k)$ of the fiducial nonlinear matter power spectrum and construct the linear response of binned galaxy-clustering, galaxy-galaxy-lensing, and cosmic-shear spectra to the discretized modulation $A(k)$. The response is evaluated with full-sky, beyond-Limber kernels including density, redshift-space-distortion, gravitational-shear, and intrinsic-alignment contributions. We reconstruct $A(k)$ using a regularized modified Richardson-Lucy algorithm, with weak diffusion in $\ln k$ and selection of the minimum-$\chi^2$ solution along the iteration history. Using Rubin/LSST Year 10-like synthetic data, we find that oscillatory modulations with amplitudes $\gtrsim1\%$ can be recovered over $0.1\lesssim k\lesssim0.5\,{\rm Mpc}^{-1}$, provided the oscillation frequency $f\lesssim10$ on $\log_{10}[k/(0.2\,{\rm Mpc}^{-1})]$. We further introduce a posterior-weighted consistency statistic calibrated with posterior-predictive null mocks, thereby accounting for cosmological and nuisance-parameter uncertainties without relying on Wilks' theorem. The null case is consistent with $A(k)=1$, while a $1\%$ oscillatory modulation is detected at $\sim 2.6\sigma$. These results demonstrate the potential of regularized deconvolution as a model-independent consistency test of the matter power spectrum in future $3\times2$pt surveys.
Read the full abstract

All authors Jun-Qian Jiang, Ajoy Dawn, Dhiraj Kumar Hazra, Benjamin L'Huillier, Arman Shafieloo

Cosmology and Nongalactic Astrophysics

Finding the distribution of matter using lenses - I: deconvolution-based reconstruction with CMB lensing

The matter power spectrum is one of the primary statistical descriptors of the large-scale distribution of matter in the Universe and provides a powerful probe of cosmic structure formation. Measurements of cosmic microwave background (CMB) lensing offer an integrated view of the matter distribution over a wide range of redshifts, enabling the reconstruction of the underlying matter power spectrum. In this work, we reconstruct the reference linear matter power spectrum $ P_\text{lin}(k,0)$ from the baseline joint CMB lensing measurements of Planck PR4, ACT DR6, and SPT-3G using a covariance-weighted modified Richardson-Lucy(MRL) deconvolution algorithm. The reconstructed spectrum is found to be consistent with the fiducial linear prediction on large scales, while exhibiting a systematic enhancement for $k \gtrsim 0.1\,{\rm Mpc}^{-1}$, where nonlinear gravitational evolution becomes important. To investigate this behavior, we introduce a scale-dependent correction factor, $A(k)$, defined through $P(k)=A(k)\,P_{\rm nl}(k),$ where $P_{\rm nl}(k)$ is the fiducial nonlinear matter power spectrum obtained from 2LPT simulations. The reconstructed correction factor remains consistent with unity within $2\sigma$ confidence over the reconstructed range, indicating that the observed enhancement is well explained by the standard nonlinear evolution of the matter power spectrum. In addition, the reconstruction shows agreement with the fiducial BAO template around the BAO feature at $k\sim(0.04-0.06)\ {\rm Mpc}^{-1}$, indicating that some BAO-scale information survives the lensing projection.
Read the full abstract

All authors Ajoy Dawn, Jun-Qian Jiang, Dhiraj Kumar Hazra, Benjamin L'Huillier, Arman Shafieloo

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

Revisiting EBL Constraints from Gamma-Ray Observations: A Critical Assessment and Methodological Improvements

The extragalactic background light (EBL), ranging from the infrared to the ultraviolet bands, is the second most intense photon field in the universe, surpassed only by the cosmic microwave background (CMB). It is primarily generated by starlight in galaxies, either directly or through absorption by dust and re-emission at longer wavelengths. Very high energy (VHE, $E$ > $100\,$GeV) photons can be absorbed via $\gamma \gamma$ interactions with the EBL during their propagation across cosmological distances, providing an indirect method to probe the EBL by studying its impact on the energy spectra of distant VHE sources. This paper examines the robustness of EBL constraints derived from gamma-ray data, critically assessing the assumptions made in previous studies about the intrinsic source spectra, the uncertainties of the observations and the validity of the analysis tools. We find that earlier studies likely underestimated the uncertainties in the EBL intensity constraints, underscoring the need to account for systematic uncertainties comprehensively. By employing a Monte Carlo (MC) simulation and a plausible model for unknown systematic errors, we compute more realistic uncertainties. Additionally, we discuss possible alternatives to set EBL constraints, relaxing the assumptions on the intrinsic spectra of gamma-ray sources, with the goal of achieving more robust constraints.
Read the full abstract

All authors K. Abe, S. Abe, J. Abhir, A. Abhishek, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, I. Albanese, T. Aniello, S. Ansoldi, A. Arbet Engels, C. Arcaro, T. T. H. Arnesen, A. Babić, C. Bakshi, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak, E. Bronzini, I. Burelli, A. Campoy-Ordaz, A. Carosi, R. Carosi, M. Carretero-Castrillo, A. J. Castro-Tirado, D. Cerasole, G. Ceribella, A. Cerviño, Y. Chai, G. Chon, A. Cifuentes Santos, J. L. Contreras, J. Cortina, S. Covino, G. D'Amico, P. Da Vela, F. Dazzi, A. De Angelis, B. De Lotto, M. Delfino, J. Delgado, C. Delgado Mendez, F. Di Pierro, R. Di Tria, L. Di Venere, A. Dinesh, D. Dominis Prester, A. Donini, D. Dorner, M. Doro, L. Eisenberger, D. Elsaesser, L. Foffano, L. Font, F. Frías García-Lago, S. Fröse, Y. Fukazawa, S. García Soto, M. Gaug, J. G. Giesbrecht Paiva, N. Giglietto, F. Giordano, P. Gliwny, T. Gradetzke, R. Grau, J. G. Green, P. Günther, D. Hadasch, A. Hahn, G. Harutyunyan, T. Hassan, J. Herrera Llorente, D. Hrupec, D. Israyelyan, J. Jahanvi, I. Jiménez Martínez, J. Jiménez Quiles, S. Kankkunen, T. Kayanoki, G. W. Kluge, J. Konrad, P. M. Kouch, H. Kubo, J. Kushida, M. Láinez, A. Lamastra, E. Lindfors, S. Lombardi, F. Longo, R. López-Coto, M. López-Moya, A. López-Oramas, S. Loporchio, L. Lulić, E. Lyard, P. Majumdar, M. Makariev, M. Mallamaci, G. Maneva, M. Manganaro, S. Mangano, K. Mannheim, M. Mariotti, M. Martínez, P. Maruševec, D. Mazin, S. Menchiari, J. Méndez Gallego, S. Menon, D. Miceli, J. M. Miranda, R. Mirzoyan, M. Molero González, E. Molina, H. A. Mondal, A. Moralejo, C. Nanci, A. Negro, V. Neustroev, M. Nievas Rosillo, C. Nigro, L. Nikolić, K. Nilsson, S. Nozaki, A. Okumura, J. Otero-Santos, S. Paiano, D. Paneque, R. Paoletti, J. M. Paredes, M. Peresano, M. Persic, M. Pihet, F. Podobnik, P. G. Prada Moroni, E. Prandini, W. Rhode, M. Ribó, J. Rico, A. Roy, N. Sahakyan, F. G. Saturni, F. Schiavone, K. Schmitz, F. Schmuckermaier, A. Sciaccaluga, G. Silvestri, A. Simongini, J. Sitarek, V. Sliusar, D. Sobczynska, A. Stamerra, J. Strišković, D. Strom, M. Strzys, Y. Suda, M. Takahashi, R. Takeishi, J. Tartera Barberà, P. Temnikov, T. Terzić, A. Tutone, S. Ubach, M. Vazquez Acosta, S. Ventura, G. Verna, I. Viale, A. Vigliano, C. F. Vigorito, E. Visentin, V. Vitale, M. Vorbrugg, I. Vovk, R. Walter, C. Walther, F. Wersig

Cosmology and Nongalactic Astrophysics

Probing Projection Effects in Optically-Selected Clusters with Velocity Dispersions Outliers

The evolution of galaxy cluster abundance is a powerful probe of cosmology, but its efficacy depends on our ability to understand and control for systematics in the selection and mass estimation of clusters. We study redMaPPer-selected clusters in the Cardinal simulation, which simulates the DES Y3 observational conditions. Specifically, we investigate the line-of-sight velocity distributions of cluster member galaxies in 12,272 simulated clusters. We find that a significant fraction of clusters, 35%, have velocity dispersions that are high compared to their richness with an offset of $\sim 400$ km s$^{-1}$ between the outliers and the main population. A similar population of velocity dispersion outliers was found in the DES Y3 data. Nearly all of the outliers (98%) in a low-richness subsample have non-Gaussian velocity distributions, and they are more likely than the main population to have significant velocity substructure. These clusters also tend to have higher redshifts than clusters with typical velocity dispersions, and they show a larger scatter between halo mass and richness. Overall, our results point to a significant contamination of redMaPPer-selected samples from line of sight structure, particularly at high redshifts and low richness. Comparison to the Buzzard simulations shows similar results.
Read the full abstract

All authors K. Graham, T. E. Jeltema, C.-H. To, V. Wetzell, Kyle W. Davis, Conghao Zhou, H.-Y. Wu, E.S. Rykoff, J. DeRose, Michel Aguena, R.H. Wechsler

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

CTree: Branch-wise compatible merger tree builder

We introduce a merger tree builder, CTree. Its key feature is the identification of progenitors using the evolutionary information encoded in object branches, rather than relying solely on the instantaneous particle membership in a single snapshot. In this sense, CTree can be described as a branch-wise merger tree builder. Specifically, CTree extracts particles associated with objects over a defined timescale and uses this information to compute the merit functions employed in the progenitor search across multiple snapshots. During this process, CTree requires only the identifiers of particles and objects. This ensures high compatibility with different halo and galaxy catalogues and also with pre-existing merger trees. To identify `core' particles that remain bound to the host halo over extended physical timescales, CTree employs a simple and efficient algorithm that does not rely on explicit radial or binding-energy information. The timescale adopted for core particle extraction and the number of snapshots used in the progenitor search can be user-configured. Basic validation tests using suites of dark matter-only simulations demonstrate that CTree produces stable and physically consistent halo merger trees, particularly in high-density regions or in simulations with high output cadence. The code is written in C++ and employs a hybrid message passing interface and OpenMP parallelization, making it well-suited for high-resolution cosmological simulations. The code is publicly available and actively maintained to improve performance, reliability, and compatibility.
Read the full abstract

Monday 7 September

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

Black Holes, Gravitational Waves and Space-Time Singularities Lemaitre Conference 2024

This editorial introduces the topical collection arising from the Lemaitre Conference 2024 - the second in a series of meetings dedicated to the scientific legacy of Georges Lemaitre - held at the Vatican Observatory in Castel Gandolfo, 17-20 June 2024, fifty-seven years after Lemaitre's death and on the eve of the centenary of his foundational 1927 paper on the expanding universe. The nineteen contributions collected here address the most pressing open problems at the interface of cosmology, gravitation, and quantum theory: the Hubble tension and the future of ${\Lambda}$CDM in light of recent DESI data; inflationary cosmology and dark-energy model-building in string theory and supergravity, together with the Swampland constraints that bound them; the observer-dependence of the quantum-cosmological wave function; the removal, or consistent crossing, of spacetime singularities - from the canonical quantization of Lemaitre's own 1933 dust model to string-theoretic pre-Big-Bang bounces and many-body constructions of singularity-free black-hole cores; primordial black holes and their gravitational-wave signatures; the search for a background-independent formulation of quantum gravity and a first-principles account of horizon thermodynamics; and the foundational status of the quantum-to-classical transition in gravitational contexts. A proposal for lunar-based cosmological observations, and a historical reconstruction of the genesis of Lemaitre's primeval-atom hypothesis, complete the collection. A recurring unifying thread was the legacy of Lemaitre himself - who first distinguished coordinate from physical singularities and coined the term "horizon" - whose pioneering vision continues to guide research at the frontier of cosmology and fundamental physics.
Read the full abstract

All authors Massimo Bianchi, Sergio Luigi Cacciatori, Matteo Galaverni, Gabriele Gionti, Fabio Scardigli

Cosmology and Nongalactic Astrophysics

Interacting fluid cosmologies from the metric-affine framework

Metric-affine gravity provides a natural geometric framework in which spacetime curvature, torsion, and non-metricity are treated as independent degrees of freedom, leading to novel cosmological dynamics beyond General Relativity. A generic consequence of such theories is the emergence of effective interactions among the cosmological fluids, even in the absence of explicit phenomenological couplings. Interestingly, the above happens while the Standard Model of Particle physics remains unmodified. In this work, we investigate interacting cosmological models that arise from metric-affine gravity and analyze their implications for the background evolution of the Universe. We derive the modified continuity equations governing the matter, radiation, and dark-energy components, highlighting the geometric origin of energy exchange in the dark sector. We discuss how these interactions modify the expansion history and can effectively mimic evolving dark-energy behavior. The resulting cosmological scenarios are confronted with the most recent observational data from Type Ia Supernovae (Pantheon+), Cosmic Chronometers, and Baryon Acoustic Oscillations (DESI DR2). Moreover, using model selection criteria, e.g., Akaike Information Criterion, Bayesian Information Criterion and the Bayesian Evidence, we compare the aforementioned models against the standard $\Lambda$CDM model. We find that the proposed framework performs slightly better than $\Lambda CDM$, in the light of the data sets considered here. Our analysis demonstrates that metric-affine-induced interactions constitute a viable and theoretically motivated alternative for the description of the Universe.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Non-linear Pairwise Velocities as a Cosmological Probe

Peculiar velocities trace gravitational dynamics directly, complementing density-based clustering as wide-field spectroscopic surveys (DESI, Euclid, PFS) and kinetic Sunyaev-Zel'dovich (kSZ) measurements enter the precision era. Observational analyses of pairwise velocities and the pairwise kSZ signal remain largely restricted to linear and quasi-linear scales ($\gtrsim20-30$ $h^{-1}Mpc$), despite substantial cosmological information available at smaller separations. We present a simulation-calibrated likelihood for $v_{12}(r,a)$, built on the exact pair conservation equation and PyCAMB HMcode-2020 non-linear clustering, and validate it against the QUIJOTE and TNG300-3 suites across resolution, particle sampling, box size, and redshift. Using QUIJOTE, we show that extending the fit to non-linear scales tightens the $1\sigma$ uncertainty on $\Omega_\mathrm{m}$ and $\sigma_8$ by $74\%$ and $81\%$ at $z=0.5$ ($r_{\rm min}=4$ $h^{-1}Mpc$), and by $70\%$ and $74\%$ at $z=0$ ($r_{\rm min}=8$ $h^{-1}Mpc$), relative to the linear-regime baseline ($r_{\rm min}=50$ $h^{-1}Mpc$) (all at fixed $r_{max}=140$ $h^{-1}Mpc$); $\sigma(f\sigma_8)$ tightens by up to $86\%$, and the joint $\Omega_\mathrm{m}$-$f\sigma_8$ Figure-of-Merit improves by up to $9.1\times$ relative to the same baseline. These constraints are conditional on $h$ fixed at its QUIJOTE fiducial value, due to a near-degenerate response with $\sigma_8$; jointly sampling $h$ shifts and widens them substantially. This framework provides a validated route to exploiting that information in upcoming direct peculiar-velocity, redshift-space, and kSZ analyses.
Read the full abstract

All authors Mariana Jaber, Antonela Taverna, Rasmus Strid, Wojciech A. Hellwing

Cosmology and Nongalactic Astrophysics

Observation-driven simulations of strong lensing galaxy clusters

Galaxy clusters are the most powerful strong lenses: they greatly magnify the flux of distant and faint sources. Strong lensing also allows for the reconstruction of their mass distribution with-1% level accuracy and enables investigating cosmological parameters. The number of known systems of this type is bound to increase in the next years thanks to wide imaging surveys. Using simulations in this context is crucial to validate the analysis methods before they are applied to real data, and to train machine learning algorithms that can handle large volumes of images. In this work, we present a simulated set of one hundred images of galaxy clusters that we have produced with a novel code for simulating cluster-scale strong lenses. One of the main novelties of our approach, distinguishing it from other existing codes, is the use of empirical relations, derived from state-of-the-art observations, for modelling the characteristics, such as morphology, color, and spatial distribution of the cluster member population. This allows us to reliably reproduce the complexity of real observations. The simulations are partly carried out with the latest version of SkyLens, a code that creates mock observations of strong lensing events in different systems and observational setups. The main improvements we introduce are the use of the message passage interface (MPI) paradigm, a standard for parallel programming that leverages the use of several processors to perform some given task, and the implementation of score-based diffusion models to augment the images of the background sources. Together, they lead to more efficient and realistic image simulations. We also present the validation of the code and simulations by comparing the properties of the mock clusters to those of real ones. [Abridged]
Read the full abstract

All authors L. Leuzzi, M. Meneghetti, A. Adam, L. Moscardini, C. Giocoli, P. Bergamini, Y. Hezaveh, M. Maturi, A. Mercurio, A. Moretti, Andrés A. Plazas Malagón, P. Rosati

Cosmology and Nongalactic Astrophysics

The SRG/eROSITA All-Sky Survey: Early dark energy and Hubble constant from the cluster mass function

The evolution of the mass function of massive galaxy clusters is a well-established probe of cosmology. Using the eROSITA X-ray instrument on board the Spectrum Roentgen Gamma (SRG) mission, the western Galactic hemisphere of eROSITA's first All-Sky Survey (eRASS1) delivers a uniformly selected and securely confirmed sample of 5259 galaxy clusters spanning $0.1 < z < 0.8$. When combined with overlapping weak lensing data from DES Year 3, KiDS, and HSC for mass calibration, this dataset enables precise tests of the standard $\Lambda$CDM framework and beyond. In this work, we use the $0.1 < z < 0.45$ subsample of the eRASS1 cluster catalog to place constraints on an axion-like Early Dark Energy (EDE) scenario. Such models have been proposed as a possible mechanism to alleviate the tension between early- and late-Universe measurements of the Hubble constant. In this framework, a scalar field temporarily enhances the cosmic expansion rate around the epoch of recombination before rapidly diluting at later times, leaving distinct imprints on structure growth. We leverage these signatures to present the first constraints on EDE derived from galaxy cluster number counts. Using eRASS1 number counts alone, we obtain an upper limit on the maximum EDE fraction of $f_{\rm EDE} < 0.3$, consistent with primary CMB analyses. Combining our cluster analysis with primary CMB measurements, baryon acoustic oscillations, and CMB lensing, while notably excluding distance-ladder calibration data, yields results consistent with a nonzero EDE contribution and $H_0$ values compatible with late-Universe measurements. The most significant detection, at $\sim 3.2\sigma$, arises from the joint analysis of the eRASS1 sample with the full set of external datasets, yielding $f_{\rm EDE} = 0.10^{+0.04}_{-0.03}$ and $H_0 = 71.4 \pm 1.4$ km/s/Mpc. [ABRIDGED]
Read the full abstract

All authors J. Strunk, E. Artis, E. Bulbul, S. Grandis, N. Clerc, V. Ghirardini, M. Kluge, A. Liu, F. Balzer, J. Comparat, I. Chiu, Z. Ding, N. Malavasi, A. Merloni, T. Mistele, H. Miyatake, S. Miyazaki, K. Nandra, F. Kleinebreil, N. Okabe, M. E. Ramos-Ceja, J. S. Sanders, T. Schrabback, R. Seppi, S. Zelmer, X. Zhang

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

Neural Posterior Estimation for Tomographic Weak Lensing Mass Mapping

Weak gravitational lensing shear and convergence trace the distribution of baryonic and dark matter across space, making them a powerful probe of cosmic structure. Inferring shear and convergence from images is a challenging inverse problem. The prevailing approach to this task estimates shear from weighted averages of galaxy ellipticities, calibrates these estimates to account for systematic biases, and transforms them to reconstruct convergence, a multistage procedure that requires substantial computational resources and meticulous handling of statistical uncertainties. As an alternative, we propose a probabilistic approach to field-level weak lensing inference in which we train a deep neural network to directly map a multiband image to a variational distribution over the underlying tomographic shear and convergence fields. This neural posterior estimation (NPE) procedure implicitly marginalizes over nuisance variables in the cosmological forward model and does not require evaluating the likelihood function. It is also amortized, so it enables rapid posterior inference for astronomical surveys once the neural network is trained. When evaluated on synthetic images from the LSST-DESC DC2 Simulated Sky Survey, NPE produces well-calibrated variational distributions for shear and convergence that are consistent with the ground truth. We describe how maps sampled from these variational distributions could be used in a subsequent simulation-based inference procedure to approximate the posterior distribution over cosmological parameters.
Read the full abstract

All authors Tim White, Shreyas Chandrashekaran, Camille Avestruz, Jeffrey Regier, the LSST Dark Energy Science Collaboration

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

Inelastic Singlet-Doublet Fermion Dark Matter in light of the 248 keV LZ event

Recently, the LUX-ZEPLIN (LZ) collaboration reported the observation of a single dark matter (DM)-nucleus scattering event at a nuclear recoil energy of $248\pm23_{\rm stat}\pm23_{\rm sys}$ keV, corresponding to an exposure of 2.84 tonne-year. The absence of events at lower nuclear recoil energies in the predicted spectrum is naturally explained if the underlying process is inelastic DM-nucleus scattering. Motivated by this, we investigate the singlet-doublet fermion DM model, in which the DM consists of two pseudo-Dirac states: the Majorana nature of the lighter state forbids tree level $Z$-mediated elastic scattering identically, while the same states enable inelastic DM-nucleus scattering via $Z$ exchange, with any residual elastic scattering proceeding only through a suppressed Higgs-mediated channel. We further extend the model with a $Z_2$-even scalar triplet, which is responsible both for generating the pseudo-Dirac splitting and for realizing Majorana neutrino masses via the Type-II seesaw mechanism.
Read the full abstract

All authors Debasish Borah, Sujit Kumar Sahoo, Narendra Sahu, Shashwat Sharma

Cosmology and Nongalactic Astrophysics

Filament-Arm Node Systems (FANS): A new probe of the winding and chirality of the cosmic web

We define Filament-Arm Node Systems (FANS), a new class of cosmic-web objects consisting of multiple filament arms connected to a common node, and use their geometry to construct an object-based test of cosmic chirality. With the cosmic-web node providing a natural reference point, FANS span characteristic scales of ~60Mpc/h, allowing the coherent winding of their arms around the node to be measured directly. We define cosmic spirality as the amplitude of this winding and handedness as its signed clockwise (CW) or counter-clockwise (CCW) sense. We further introduce an axis-conditioned extrinsic estimator that searches for the projected signature of a coherent three-dimensional axial pattern of the handedness field, corresponding to a preferred axis of winding. Applying these measurements to two distinct SDSS filament reconstructions, and assessing their significance with sign-flip nulls and 27 controlled variations of the FANS construction parameters, we find no significant preference for CW or CCW winding, either globally or around a preferred axis. The data are therefore consistent with parity symmetry, as expected in standard LCDM cosmology. We additionally find that cosmic spirality decreases with cosmic time, a trend that may reflect the gravitational straightening of filament arms. These first constraints establish FANS as a new probe of cosmic-web symmetry and large-scale morphology and provide a baseline for future parity tests with the much larger filament samples expected from upcoming surveys.
Read the full abstract

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

LZ nuclear recoil event from inelastic singlet-doublet scalar dark matter

We study the possibility of explaining the recently reported high-energy nuclear recoil event by the LUX-ZEPLIN (LZ) collaboration within the framework of singlet-doublet scalar dark matter (DM). Considering DM to be the CP even mass eigenstate formed out of a scalar doublet and a real scalar singlet, both being odd under an unbroken $Z_2$ symmetry, we find the parameter space of the model consistent with correct relic abundance and direct-detection limits on elastic scattering rates. A large part of the parameter space can also lead to inelastic up-scattering of DM with a rate and recoil energy consistent with the recent LZ event. Depending upon the singlet-doublet mixing, the model allows a much wider range of currently allowed parameter space compared to the pure scalar doublet DM limit, which is disfavored due to non-observations of high-energy neutrinos from solar-captured DM annihilation by the IceCube experiment. Presence of $Z_2$-odd right-handed neutrinos also leads to other interesting phenomenology related to the origin of light neutrino masses and leptogenesis.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Ultra-slow-roll Inflation with Non-perturbative Non-Gaussianity and Scalar Induced Gravitational Waves

An ultra-slow-roll phase during inflation could potentially produce the right abundance of primordial black holes (PBHs) for them to form all of dark matter (DM). We consider such a scenario, carefully treating the transitions from slow to ultra slow roll inflation and back, using a parametrisation of the inflaton potential that can describe inflation from the time when CMB scales crossed out of the horizon until its end. A $\delta N$ analysis shows that this model can possess $O(1)$ non-Gaussianity. When computing the primordial black hole abundance, we keep the full non-linear relation between curvature perturbations and the density contrast and consider the non-Gaussianity to all orders. We find that $O(1)$ non-Gaussianity is sufficient to enhance the PBH abundance by orders of magnitude, leading to reduced gravitational wave (GW) signals for fixed PBH abundance. The signal to noise ratio is computed for future gravitational wave observatories. We find that, in spite of the reduced strength, the GW signal is well above the sensitivity of LISA if PBHs form a significant component of DM.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Recovering Large-Scale Clustering of Missing Galaxies for Galaxy–Gravitational-Wave Cross-correlations

We present a framework for completing galaxy catalogs from flux-limited surveys which reconstructs the missing galaxy population while preserving its clustering properties. This is particularly important for cosmological analyses using gravitational-wave (GW) dark sirens without electromagnetic counterparts, such as binary black hole mergers, which constitute the majority of GW detections. As GW detectors probe increasingly larger distances, galaxy catalogs become progressively incomplete owing to survey flux limits. Current state-of-the-art statistical host identification methods typically account for this incompleteness by assuming that the missing galaxies are uniformly distributed in comoving volume. While this approximation can mitigate biases due to incompleteness, it neglects the clustering of galaxies within the cosmic web. Since galaxies and GW sources are both expected to trace the underlying large-scale structure, their spatial distributions are correlated, making a homogeneous reconstruction physically unrealistic, a shortcoming which our framework addresses. As a proof of concept, we apply the method to simulated galaxy catalogs with different flux limits and observational selection functions. Our approach provides a significantly more realistic completion of incomplete galaxy catalogs than the commonly adopted homogeneous-in-comoving-volume approximation, enabling more robust cosmological inference from gravitational-wave dark sirens.
Read the full abstract

Cosmology and Nongalactic Astrophysics

Conformally Interacting Dark Energy with Early and Late-Time Measurements

The conformal interacting dark energy (CIDE) model is investigated by introducing a scalar field representing dark energy (DE) coupled to dark matter (DM) via a conformal transformation, thereby yielding a class of scalar-tensor theories. The interaction term ${Q}$, that indicts the flow of energy between the dark sectors is proportional to the trace of the energy momentum tensor of the DM fluid, $T^{\mathrm{DM}}_{\mu\nu}$, as ${Q}=\frac{C'(\phi)}{2C(\phi)}g^{\mu\nu}T^{\mathrm{DM}}_{\mu\nu}$, where $C(\phi)$ is the conformal function of a scalar potential (coupling) field. We assume the power-law parametrization $C(\phi)\propto (1+\phi)^m$, with $m$ a coupling parameter, and determine the direction and magnitude of the energy flow between the DM and DE components. The dynamical nature of DE is modeled by the two conformal interacting scenarios, CIDE and $w$CIDE, with $w$ the DE equation of state parameter. To test the viability of each model, we constrain them using a combination of early- and late-time cosmological data, namely: CMB measurements from the South Pole Telescope, Planck 2018, and the Atacama Cosmology Telescope (DR6) (CMB SPA); BAO data from the Dark Energy Survey (DESI DR2 BAO); and Supernova Type Ia distance compilations (PantheonPlus, PantheonPlus + SH0ES, Union3, and DES Dovekie). Parameter inference is performed with Monte Carlo Markov Chain (MCMC) simulations using COBAYA and a modified CLASS. Further statistical analysis using the Akaike information criterion (AIC) and the Bayesian information criterion (BIC) is summarized to assess the viability of the model in comparison with the standard cosmological model, investigating the model's potential to alleviate the $H_0$ and $S_8$ tensions.
Read the full abstract

All authors Shambel Sahlu, Abunie Gezahegn, Amare Abebe, Gonzalo J. Olmo, Diego Rubiera-Garcia

Cosmology and Nongalactic Astrophysics

Primordial Black Hole Abundances and Scalar Induced Gravitational Waves from Finite-Width Power Spectra in a Stiff Thermal History

We study the formation of primordial black holes (PBHs) from large primordial perturbations that re-enter the horizon during an epoch with equation of state ${\rm w}\geq1/3$. We consider a log-normal curvature power spectrum of finite width $\Delta$ and determine the collapse amplitude by numerical-relativity simulations of a self-gravitating perfect fluid. Threshold scans are performed for five values of $\rm w$ and five spectral widths, and the resulting numerical thresholds are used in the PBH abundance and scalar-induced gravitational-wave (SIGW) calculations. For comparison, we also evaluate the semi-analytical $q$-function prescription. It reproduces the numerical trend close to radiation domination and for nearly monochromatic profiles, but it is not a reliable threshold estimator for generic finite-width profiles in stiffer backgrounds. We show how collapse thresholds increase with both $\Delta$ and $\rm w$, changing the curvature amplitude required for PBHs to constitute all of the dark matter and, consequently, the normalization of the accompanying SIGW signal.
Read the full abstract

All authors Abolhassan Mohammadi, Yogesh, Cristian Joana, Hongwei Tan, Ying-li Zhang

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

Recovering the characteristic halo mass of rare populations from clustering: a benchmark using the Millennium Simulation

We measure how accurately the characteristic halo mass of a rare, massive population can be recovered from clustering alone, using dark matter haloes in the Millennium simulation, where the true masses are known. Massive samples of 25 to 10000 haloes are cross-correlated against field reference samples of 5000 to 50000 haloes at six redshifts over $0 \leq z \leq 5.29$. The bias is measured at $r = 8\,h^{-1}\,\mathrm{Mpc}$ and inverted to return a characteristic mass, which we compare with the mean logarithmic halo mass measured directly from the simulation. Across 192 configurations the characteristic mass is recovered with a median absolute offset of $0.103$ dex. We show that (i) accuracy depends only weakly on the number of rare objects while the uncertainty falls from $\pm 0.71$ dex at 25 objects to $\pm 0.08$ dex at 10000; (ii) enlarging the smallest samples gives the greatest return, with a doubling from 25 to 50 objects reducing the uncertainty by about a third, against a tenth from 5000 to 10000; and (iii) the reference population carries an error of its own, which propagates coherently into every cross-correlation measured against it and is not reduced by observing more rare objects. Working with haloes rather than galaxies isolates the clustering measurement from the galaxy–halo connection, providing a benchmark before population-specific complications are introduced. For the best configurations, the uncertainty contributed by the calibration of the bias–mass relation, at $\sim 0.07$ dex, already matches that of the measurement itself.
Read the full abstract

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

NEXUS: Transient Searches and First Results from Year One Observations

We describe ongoing efforts of high-redshift transient searches using multi-epoch NIRCam imaging (F200W+F444W) and NIRSpec MSA/PRISM spectroscopy from the NEXUS JWST multi-cycle Treasury program targeting the north ecliptic pole region. The transient search area covers $\sim 61.5\,{\rm arcmin^2}$ between the reference epoch and each subsequent NEXUS-Deep epoch at a cadence of $\sim2$~months. In the first year of observations from NEXUS, we detect 68 robust transients, with the host photometric redshift distribution declining rapidly at $z>2$ but extending to $z_{\rm phot}\approx 6$. In addition, we obtained secure spectroscopic redshifts for 37 transients ($\sim54\%$) from NIRSpec/PRISM and NIRCam/WFSS, with seventeen at $1 < z < 2$, eight at $2 < z < 3$, two at $3 < z < 4$, and one tentative host association at z = 6.151. Overall, the NEXUS program recovers observed supernova (SN) rates broadly consistent with other SN search programs with JWST. While NEXUS achieves the highest transient detection efficiency, 6.6 SNe per imaging hr ($2.6\times$ COSMOS-SN and $26\times$ JADES-SN), the limited filter coverage (F200W+F444W only) limits robust identification and classification of high-$z$ SNe for deep spectroscopic follow-up. We describe the details of the data reduction and transient detection pipeline, and report the Year-1 transient sample along with their light curves, available MSA spectroscopy, host association, and the raw detection rate. We also describe and release a new PSF photometry package that properly accounts for correlated pixel noise from combining drizzled images.
Read the full abstract

All authors Ming-Yang Zhuang, Lei Hu, Justin D. R. Pierel, Zhiwei Pan, Yue Shen, Adam J. Burgasser, Xiaohui Fan, Luis C. Ho, Junyao Li, Alice E. Shapley, Zachary Stone, Feige Wang, Lifan Wang

Cosmology and Nongalactic Astrophysics

Little Red Dot Cosmology: A Matter-Era Baryon Acoustic Oscillations Probe of $\Lambda$CDM

The discovery of Little Red Dots (LRDs) with JWST provides a new source population for cosmology in a largely unmapped period of cosmic history. LRDs are most numerous at $4\lesssim z\lesssim9$, thereby bridging existing low-redshift large-scale-structure measurements and the cosmic microwave background. In this epoch, the Universe is deep in the matter-dominated era, where dark energy is dynamically negligible and the expansion history is tightly predicted in $\Lambda$CDM, making these redshifts a clean test of the standard cosmological model. Recent measurements indicate that LRDs have number densities of $\bar n\sim10^{-4}$ $h^3$ Mpc$^{-3}$, a bias rising from $b\sim 3$ to $\sim8$ over this interval, and distinctive spectrophotometric signatures – a remarkably favorable combination for baryon acoustic oscillation (BAO) measurements. Using a Fisher forecast with these fiducial inputs, we find that a wide-field spectroscopic survey of LRDs over a DESI-like (14,000 deg$^2$) footprint delivers a percent-level isotropic distance scale measurement of $D_V/r_d$ across four bins spanning $4\lesssim z\lesssim9$. LRDs are therefore a compelling target for future wide-field near- to mid-infrared spectroscopy aimed at mapping large-scale structure deep into the matter-dominated era.
Read the full abstract

Sunday 6 September

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

Inelastic Self-interacting Dark Matter and LUX-ZEPLIN 248 keV Event in a Dirac Modular Inverse Seesaw

We propose a novel framework that simultaneously addresses the origin of Dirac neutrino masses and the nature of self-interacting dark matter (SIDM). The model is based on an $A_{4}$ modular symmetry to ensure the Diracness of neutrinos as well as the stability of the DM. The neutrino sector realizes a Dirac Inverse Seesaw mechanism where the smallness of the neutrino mass is governed by the vacuum expectation value (VEV) of a singlet scalar $\phi$. This same scalar couples to a vector-like fermion DM candidate, inducing a tiny Majorana mass splitting that renders the DM pseudo-Dirac and inelastic. Crucially, the scalar also acts as a light mediator for DM self-interactions, potentially solving the small-scale structure problems of Cold DM. In light of the recent 248 keV nuclear-recoil event, LZ230616, observed by LUX-ZEPLIN (LZ) DM direct detection experiment, we demonstrate that our inelastic SIDM parameter space naturally accommodates this signal via endothermic scattering kinematics. Furthermore, the spontaneous breaking of the dark parity required for this inelasticity produces a network of cosmological domain walls. We show that the explicit symmetry breaking needed to safely annihilate these walls generates a stochastic gravitational wave background. The non-holomorphic modular symmetry reduces the free parameters, correlating neutrino observables, addressing DM phenomenology, $\Delta N_{\rm eff}$, and gravitational-wave signatures.
Read the full abstract

All authors Pritam Das, Biswajit Karmakar, Satyabrata Mahapatra, Partha Kumar Paul

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

ADF22-WEB: Massive galaxy formation shaped by cosmic web filaments in the z = 3.1 proto-cluster core

We present a census of dust, molecular gas, and galaxy structure in 18 dusty star-forming galaxies (DSFGs) at $z=3.09$ embedded in Ly$\alpha$-traced cosmic web filaments in the core of the SSA22 proto-cluster, using multi-band ALMA, JWST imaging, and JVLA CO(1–0). Using up to six ALMA bands and Herschel/SPIRE data, we construct rest-frame far-infrared spectral energy distributions. The DSFGs span two orders of magnitude in far-infrared luminosity ($L_{\rm FIR}\sim10^{11}$–$10^{13}L_\odot$), with median values of $\log(L_{\rm FIR}/L_\odot)=11.56^{+0.32}_{-0.37}$ and $T_{\rm dust}=25^{+6}_{-3}{\rm K}$. For eight DSFGs with high-resolution ($\sim1$ kpc) JWST and ALMA imaging, we measure stellar and dust surface densities and find a positive correlation, suggesting that structural compaction of the stellar component proceeds together with that of the interstellar medium. We report 12 detections of CO(1–0), 18 of CO(3–2), four each of CO(8–7) and CO(9–8), and one of CO(12–11). The median brightness temperature ratio between CO(3–2) and CO(1–0) is $r_{31}=0.66^{+0.05}_{-0.04}$, consistent with field galaxies at similar redshifts. The high-$J$ CO lines show relatively low excitation, and the CO(8–7)/CO(3–2) ratio correlates with star-formation rate surface density, suggesting that molecular gas excitation is primarily regulated by star formation. Lower-mass galaxies show excess molecular gas fractions and depletion times relative to field scaling relations. By combining quiescent galaxy samples in the same field, we identify a sequence from gas-rich systems with low stellar surface density to gas-poor, compact systems with high stellar surface density. These results suggest that massive galaxy evolution in dense environments is governed by a baryon cycle linking gas supply, star formation, and structural transformation within the cosmic web environment.
Read the full abstract

All authors Hideki Umehata, Ian Smail, Charles C. Steidel, Bret D. Lehmer, Kouichiro Nakanishi, Bunyo Hatsukade, Mariko Kubo, Norika Okauchi, Yoichi Tamura, Masato Hagimoto

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

A Common $\sim$55 GeV Inelastic Dark Matter Origin of the LZ, Fermi-LAT, and AMS-02 Excesses and Paleo-Detector Prospects

We investigate the 248~keV nuclear recoil excess reported by LUX-ZEPLIN (LZ) in a leptophobic inelastic dark matter framework, in which the mass splitting naturally suppresses low-energy recoils. We find two viable branches, including the endothermic solution with a high dark matter mass $m_\chi \gtrsim400$~GeV and the exothermic one with a mass of tens GeV. Remarkably, the exothermic dark matter scenario with $m_\chi\sim55$~GeV can simultaneously account for the LZ excess, the Galactic Center GeV gamma-ray excess, and the AMS-02 antiproton excess. The favored annihilation cross-section is slightly below the one required for the observed relic abundance, resulting in an overabundant relic in the minimal setup. An additional neutrino-portal mediator can provide the required depletion while leaving the gamma-ray and antiproton signals essentially unchanged. We further show that paleo-detectors provide a powerful probe of the inelastic interpretation, requiring exposures of $\sim10$ and $\sim100$~mg Gyr for the endothermic and exothermic scenarios, respectively.
Read the full abstract

All authors Meiwen Yang, Quan-feng Wu, Kai-Kai Duan, Zi-Qing Xia, Yue-Lin Sming Tsai, Yi-Zhong Fan

Cosmology and Nongalactic Astrophysics · Cross-listed from stat.ML

Recovering Weak Signals with Normalizing Flows

In many scientific disciplines, weak signals of interest are obscured by dominant nuisance signals that are several orders of magnitude stronger. Recovering these weak signals requires subtracting the dominant ones; however, this calibration process inherently distorts or partially suppresses the underlying signal of interest. To address this problem, we propose the use of normalizing flow models to reconstruct calibration-affected weak signals. By leveraging the statistical invariance of the target signals and assuming minimal initial suppression, our framework effectively recovers the lost signal components. We provide a comprehensive theoretical overview of this normalizing flow-based recovery method and demonstrate its efficacy using simulated data.
Read the full abstract