2026-08-19 | | Total: 18
Variations in dynamical states of galaxy clusters can introduce biases and scatter in observable-mass relations. The dynamical state of a cluster is an emergent feature of its mass accretion history (MAH), it is therefore useful to constrain the MAH of the cluster. In this work, we characterize 305 massive clusters from The300 project by connecting features from their projected stellar distributions to their mass accretion histories (MAH). As a baseline, we first correlate host dark matter halo dynamical state indicators at $z=0$ with their MAH via the Spearman rank correlation coefficient $ρ_{\mathrm{sp}}$. Both substructure mass fraction and center-of-mass offset measurements correlate strongly with the MAH measured between $0.1\lesssim z\lesssim 1$. We repeat this exercise with morphological measurements of projected stellar density maps, many of which exhibit moderate correlation strength with different times in the MAH. Broadly, core morphological measurements ($r \leq 30\,\mathrm{kpc}$) correlate better with early-time MAH. Core-excised ($50\,\mathrm{kpc} \leq r \leq 1\,\mathrm{Mpc}$) morphological measurements correlate better with late-time MAH. We further quantify the MAH prediction power of both traditional dynamical state indicators and morphological parameters using Multivariable Conditional Abundance Matching (MultiCAM). MultiCAM employs simple rank-ordering operations, making it straightforward to translate to observed datasets. We find reasonable ($ρ_{\mathrm{sp}} \geq 0.6$) performance for predictions of the mass fraction between $1\lesssim z\lesssim 0.1$, though with notable information loss when using projected quantities. In one example application of our methodology, we use the coefficients of the MultiCAM models to select subsamples of galaxy clusters that have accreted more (or less) of their $z = 0$ mass budget over a given time frame.
We aim to exploit the large area coverage, good sensitivity, and low instrumental background of eROSITA to detect the faint surface brightness emission of galaxy groups from the Two Micron All Sky Survey Redshift Survey (2MRS). Using the data from eROSITA-DE Data Release 1, including images, exposure maps, and local background maps, we performed a wavelet decomposition of image mosaics in the 0.6--2.3 keV band at angular scales of 1/8-16'. We adopted 8-16' scales for source detection and 2-4' scales to improve catalog purity. A novel identification method based on the ranked partial Hausdorff distance fully exploits the X-ray image and group membership information. Random catalogs were used to control match purity, and the identification threshold was chosen to maximize the catalog size at a fixed purity. {We present a catalog of 619 X-ray galaxy groups with 80% purity, and define subsamples with 90% and 97% purity. Bright sources closely match the AXES-2MRS catalog (which is based on ROSAT All Sky Survey data analysis on spatial scales of 12-24'). The X-ray luminosity function of our groups agrees with previous studies down to 5.e41 erg/s. Using dynamical mass estimates, we find that the X-ray counterpart completeness for groups with >=4 members exceeds 60% for masses >2e13 Msun. We modeled the 2MRS group catalog and justify the inclusion of two-member groups in the identification. This study demonstrates that large X-ray sources on spatial scales relevant for cosmological studies of baryonic distributions can be reliably detected and identified using nearby galaxy group catalogs.
We investigate infrared line signatures from radiatively decaying sterile-neutrino dark matter using publicly available JWST/NIRSpec IFU blank-sky observations. The main signal considered is the sterile-to-sterile transition $N_1\to N_2γ$, induced by the transition magnetic dipole coefficient $d_{NNγ}$, with $m_1>m_2$. In contrast to ordinary two-photon decays of axion-like or Majoron-like particles, the observed photon energy is not fixed by the full dark matter mass, but by the small mass splitting $Δm=m_1-m_2$. Thus, a keV-scale sterile-neutrino dark matter state can generate an eV-scale infrared photon line in the JWST band. We construct a $χ^2$-based line-search analysis using the NIRSpec IFU $\rm F170LP$-$\rm G235M$ blank-sky data toward $\rm GN\text{-}z11$, modelling the smooth continuum with a cubic spline and including the Milky Way halo decay flux. In the absence of a significant excess, we derive projected limits on $d_{NNγ}$ and on the decay width $Γ_{N_1\to N_2γ}$ for $0.1~{\rm eV}\lesssimΔm\lesssim1~{\rm eV}$. For a sterile component saturating the dark matter abundance, the strongest sensitivity reaches $d_{NNγ}\lesssim7\times10^{-14}~{\rm GeV}^{-1}$ and $Γ_{N_1\to N_2γ}\lesssim10^{-25}~{\rm s}^{-1}$. We also include a minimal anomalous-Majoron benchmark, $ω\toγγ$, obtaining JWST sensitivity to $λ_{ωγγ}\sim10^{-11}$-$10^{-9}~{\rm GeV}^{-1}$ for $m_ω\sim0.6$-$1~{\rm eV}$.
Fast Radio Bursts probe ionised baryons through their observed dispersion measures. We present \turbofrb, a semi-analytic stochastic forward model for the intergalactic dispersion-measure distribution, $P({\rm DM}_{\rm IGM}\mid z)$, that resolves the diffuse IGM, halo, and filament contributions as explicit physical channels, with the halo and filament encounter rates coupled by a latent line-of-sight environmental variable. Only four effective parameters are calibrated against hydrodynamical ray-traced IllustrisTNG benchmark. The model matches the benchmark mean DM to the percent level and yields a per-redshift Jensen-Shannon divergence of at most $5\times10^{-3}$ across $z = 0.5$-$2.5$. The per-sightline channel decomposition makes explicit what closed-form parametric descriptions cannot show: the diffuse IGM sets the body of the distribution, while halos and filaments populate the high-DM tail. Applied to representative localised FRBs, the forward likelihood quantifies host-excess events independently of their astrophysical signatures and recovers the injected $H_0$ within $1σ$ in a closed-loop consistency test. The \turbofrb package is available at \href{https://github.com/jefersonfortunato/turbofrb}{github.com/jefersonfortunato/turbofrb}.
We introduce the Dianoga set of cosmological simulations of galaxy clusters and groups, specifically aimed at studying the impact of the implementation of AGN feedback and star formation. Using the OpenGadget3 code, we carry out simulations of 28 regions centred on massive galaxy clusters, and of a cosmological box. This generates a sample of 293 halos with M_{200}> 1.5 x 10^{13} M_{\odot}. Parameters of AGN feedback in the reference implementation were minimally calibrated exclusively to match the local relation between SMBH masses and stellar masses of host galaxies. Simulations are compared to observed galaxy stellar mass function (GSMF), stellar mass fraction in clusters and groups, BCG masses, scaling relations between ICM/IGM properties and profiles of their thermodynamical properties. In the appendix, we show how results vary as we modify the reference feedback model in six alternative configurations. Our reference model predicts a GSMF in general agreement with observations, albeit overestimated in the high end. BCG stellar masses and mass fractions are higher than observed in massive clusters, while being closer to observations for groups. Predicted properties of the ICM/IGM are in general agreement with observations, with the core regions of simulated clusters having entropy and temperature profiles that are slightly less "cool-cored" than observed. A comparison with other implementations of AGN feedback highlights that models including thermal evaporation of the sub-resolution interstellar medium succeed to bring BCG masses and stellar mass fractions closer to observation, and to increase the cool-coreness of simulated clusters. Our results demonstrate that the details of the interface between AGN energy injection and the sub-resolution interstellar medium model are at least as critical as the total feedback efficiency itself.
Early Dark Energy (EDE) models have been suggested as a possible solution to the so-called Hubble Tension, a discrepancy of the measurement of the Hubble constant at early and late times. In this paper, we investigate the effects of EDE on large-scale structure probes of cosmology via the EDENS (Early Dark Energy N-body Simulations) suite. The EDENS suite extends the reach of available EDE simulations considerably by adding volume and resolution. We derive several key metrics, such as the halo mass function, the nonlinear power spectrum, and the concentration-mass relation. Furthermore, we implement a halo occupation distribution model to populate the simulations with synthetic galaxies. This allows us to measure the galaxy-galaxy correlation function and galaxy bias. We choose an EDE model that is consistent with observations across several probes and allows for the increase of the present day value of the Hubble constant to resolve the Hubble tension. This model adds three more parameters to the standard $Λ$CDM model. Additionally, it requires small shifts in the best-fit $Λ$CDM cosmological parameters to accommodate existing observational constraints. We find significant differences between the standard $Λ$CDM model and the EDE model, suggesting that some of our chosen metrics may allow us to distinguish EDE from $Λ$CDM, and future observations should help further constrain possible cosmological models. We release outputs from our simulations via the OpenCosmo data portal.
We introduce a unified framework for dark energy diagnostics based on the joint phase space of the equation of state $w$ and the sound speed $c_s^2$. The resulting $w-c_s^2$ plane provides a minimal extension beyond background cosmology, capturing both the expansion history and perturbative properties within a single representation. Building on this framework, we define the microphysical flow parameter $F=dc_s^2/dw$, which encodes the dynamical relation between background evolution and perturbative response. We derive a direct connection between the present-day value $F_0$, $H_0$, and $σ_8$, and show that the microphysical flow parameter enables a hierarchy of increasingly stringent consistency tests that substantially reduce the viable dark energy theory space. We further demonstrate how trajectories in the $w-c_s^2$ plane distinguish models that are nearly degenerate at the level of $w(a)$, including canonical quintessence, Chaplygin gas models, and noncanonical scalar field realizations. This framework provides a compact phenomenological bridge between dark energy microphysics and future perturbation-sensitive observations, establishing $F_0$ as a useful discriminator of the kinetic structure underlying cosmic acceleration.
Finch et al (2002) derived relations for the summed absolute magnifications of quadruply lensed images that vary inversely with flattening in two alternative isothermal gravitational potential models. However, they did not elaborate on the selection effects this "flatter-fainter'' relation induces in actual lensed systems. We test the relation against the Luhtaru et al (2021) sample of 39 quadruply lensed quasars (38 of which we model successfully), using predicted rather than observed magnifications to avoid the complication of microlensing. We found that the summed predicted magnification decreases by a factor of ten over the observed range of flattening.
Galaxy shapes are unique tensor tracers of large-scale structure, providing a promising avenue to both enhance current cosmological programs and detect new physics beyond the scalar sector. We develop a general formalism to describe the full 3D structure of galaxy shapes and their statistics, including the breaking of isotropy by the line of sight and redshift space distortions. We constructively show that the redshift-space mapping generates a kinematic basis whose form factors are strictly polynomial in the line-of-sight angle $μ= \hat{k} \cdot \hat{n}$, and that parity selection rules restrict scalar-tensor and tensor-tensor correlators to 3 and 13 independent form factors, respectively, with the latter further reduced to 9 by exchange symmetry. We show that this polynomiality is preserved transforming into a total helicity basis denoted by total angular momentum $M$---sourced to be nonzero by powers of the line-of-sight $\hat{n}$---and that this is equivalent to the form factors having spin weights $(1 - μ^2)^{|M|/2}$ and described by associated Legendre polynomials $P^{|M|}_\ell$. We construct estimators for form factors in the normalized total-helicity basis that provide the optimal angular weighting to extract shape information, and connect the full tensor basis to projected shape statistics. To validate our formalism, we study the above tensor form factors within a toy model, showing that all allowed channels are generated even within simplified assumptions, and apply our estimators to halo shape statistics in N-body simulations wherein all channels up to total angular momentum $|M| \leq 2$ are detected. We anticipate that the methods developed here will have applications ranging from optimal extraction of intrinsic-alignments in hydrodynamical simulations to identifying new physics in tensor channels forbidden by selection rules in the standard model of cosmology.
Stochastic gravitational-wave backgrounds from different primordial mechanisms can produce identical peak frequencies and amplitudes, obscuring their physical origin. We show that the broadband spectral shape breaks this degeneracy. Comparing scalar-induced and gauge-field-induced gravitational waves for identical primordial spectral profiles and cosmological evolution, we find that their spectra retain distinct source-dependent signatures even after matching both the peak frequency and amplitude. While the two mechanisms share a universal deep-infrared behavior, their tensor-source kernels generate different peak morphologies and ultraviolet asymptotic scalings. In particular, the ultraviolet spectrum provides a direct diagnostic of the underlying source mechanism that is insensitive to the overall normalization. Our analysis shows that broadband spectral information contains additional source-dependent signatures that are inaccessible from peak observables alone.
Merging galaxy clusters represent one of the most powerful macroscopic laboratories in the Universe for searching for new physics within the dark sector. High-velocity cosmic collisions inherently separate the dark matter and stellar components from the highly collisional, X-ray-emitting intracluster gas. These massive systems provide an ideal environment to probe the fundamental nature of dark matter, specifically testing whether it behaves as a strictly collisionless particle or exhibits non-zero self-interactions. While pioneering systems like the Bullet Cluster historically demonstrated the macroscopic decoupling of dark and ordinary matter, the field has evolved into a sophisticated discipline driven by multi-disciplinary methodologies. This review synthesizes recent theoretical and empirical advances in interpreting post-collision dynamics. It examines how the synergy of combined approaches---integrating multi-wavelength observations from gravitational lensing and X-ray mapping with high-fidelity N-body hydrodynamical simulations---allows the translation of macroscopic spatial observables into stringent constraints on microscopic particle properties. Through this synthesis, the work evaluates how leveraging heterogeneous merger ensembles can reliably advance the ongoing search for physics beyond the standard cosmological model.
We derive a compact analytic formula for the one-loop triangle correlator of conformally coupled scalars in de Sitter space. The result is organised as six leading-singularity prefactors multiplying pure weight-two functions of the six energy variables. It contains forty-two dilogarithms, compared with approximately one hundred and twenty in the previously known closed-form representation, and requires no auxiliary regulator. The dilogarithms occur in Galois-conjugate pairs, making each contribution separately real throughout the physical region. We validate the result numerically. Moreover, we show that its symbol can be derived directly from the dressed integral representation or, independently, bootstrapped from Landau singularities and general consistency conditions. Finally, we show that the correlator (in a suitable normalization) is a Stieltjes function of each squared energy separately and is jointly completely monotone in all six squared energies. These structures suggest a route towards higher-point one-loop cosmological correlators.
Non-gravitational processes that expel hot gas beyond the virial regions of groups and clusters of galaxies, known collectively as baryonic feedback, play a key role in reshaping the matter distribution of the Universe on Mpc scales. We use eROSITA observations of a complete sample of 25 galaxy groups selected from the first public release of the eROSITA-DE data (eRASS1) and identified with the Two Micron Redshift Survey optical group catalogue (2MRS). We extract and fit surface brightness (SBx) profiles and present hot gas mass and hot gas fraction profiles out to $R_{200}$. We perform a Bayesian analysis of $M_{\mathrm{gas}}-M_{\mathrm{tot}}$, $L_{\mathrm{X}}-M_{\mathrm{tot}}$, and $L_{\mathrm{X}}-M_{\mathrm{gas}}$ relations, taking into account the aperture covariance effects. At $R_{500}$, we report uniformly flat SBx profiles with a mean $β$ parameter of $0.38 \pm0.04$, steepening to $β= 0.76\pm0.19$ beyond $R_{500}$. We measure a sub-cosmic hot gas fraction at the median mass of our sample $M_{500} = 2.54\times10^{13}M_{\odot}$ of $ f_{\mathrm{gas,500}} = 4.32\pm0.42\%$. Similarly, at $R_{200}$ and the median mass $M_{\mathrm{ 200}} = 3.69\times10^{13}M_{\odot}$, we obtain $f_{\mathrm{gas,200}}=5.78\pm0.69\%$. Our $f_{\mathrm{gas}}-M_{\mathrm{tot}}$ and $L_{\mathrm{X}}-M_{\mathrm{tot}}$ relations show significant deviations from the predictions of the strong feedback variants of the FLAMINGO simulation ($2.5σ$ to $8.0σ$ tension), while fiducial FLAMINGO and BAHAMAS provide the closest match to our measurements. Using our measured baryon fractions and the \texttt{SP(k)} model, we infer a $10\%-15\%$ reduction in the matter power spectrum at $k = 5\ h\ \mathrm{Mpc}^{-1}$ relative to a dark matter-only universe, in agreement with fiducial FLAMINGO and BAHAMAS, while revealing a growing tension on smaller scales with the strong feedback variants.
Modern machine learning is leading to substantial gains in precision, flexibility, and computational efficiency in fundamental physics. Statistical validation, uncertainty quantification, and robustness assessment are less systematically addressed. The VERaiPHY initiative (Validation & Evaluation for Robust AI in PHYsics) is a series of articles developed within the PHYSTAT programme, aimed at establishing statistical standards for the development, evaluation, and deployment of ML techniques. Each article focuses on a specific methodological domain from a statistics perspective and clarifies statistical questions, tests, and the interpretation of results. This opening article establishes the probabilistic, statistical, and machine learning foundations that the later contributions assume, together with the notation used throughout.
It has been proposed that primordial quark pellets or PQPs --ultra-dense quark-matter mini- stars-- formed at $T\sim 1$ GeV maybe a good candidate accounting for the dark matter of the universe. If correct, dark matter would consist of very compact objects with a maximum mass of $10^{-2} \ M_\odot$ and radii of approximately 100 m, although smaller objects would be much more abundant and encompass the bulk of the dark matter content. Here we describe a viable formation mechanism based on the enhancement of the local baryon density when supra-horizon Peccei-Quinn domain walls formed at an earlier epoch sweep and accumulate quarks and gluons before entering the horizon. Assuming an efficient baryon concentration by contracting Peccei-Quinn domain walls, we derive the resulting primordial quark pellet mass spectrum, minimum mass and cosmological abundance The results confirm that PQPs potentially constitute a conservative, Standard-Model-based, and observationally viable solution to the dark-matter puzzle.
We report the discovery of AT 2020afjz (TSS2020a): the first hour-scale extragalactic transient discovered in optical wavelengths whose complete evolution -- from explosion onset to decay -- is temporally resolved, and the first such transient discovered by TESS. AT 2020afjz was identified as a $>10σ$ detection in the pilot HiLaTS program run within the TESSELLATE Sky Survey, which blindly searches for transient phenomena in TESS data with the TESSELLATE pipeline. Through cross-matching with legacy imaging, we associate it with DES J042144.37$-$383311.3, a member of an interacting galaxy pair at $z_{\rm phot}=0.67^{+0.07}_{-0.10}$. While AT 2020afjz is similar in duration and brightness to GRB afterglows, it exhibits a slow rise time of $\sim1$ hr and lasts for only 2.4 hr above the half-max brightness; modeling the TESS light curve with VegasAfterglow finds that it is best described as either an on-axis "dirty-fireball" or off-axis orphan afterglow. Each of these rare classifications hinge upon a non-detection at gamma-ray energies, but as Fermi-GBM was Earth-occulted at the time of explosion, AT 2020afjz's gamma-quiet nature cannot be definitively confirmed. Regardless, AT 2020afjz demonstrates TESS's power to discover fast extragalactic transients, and heralds a new population awaiting discovery with TESSELLATE.
The detection of the global 21-cm signal from the Cosmic Dawn and Epoch of Reionisation remains a challenge to experiments worldwide. Emitted at a rest-frame frequency of 1420.405~MHz, this signal from the early Universe is redshifted to 40-200~MHz with a maximum brightness temperature of a few 100~mK. Efforts to detect this sky-averaged signal include experiments such as the Shaped Antenna measurement of the background RAdio Spectrum (SARAS) and Probing ReionizATion of the Universe using Signal from Hydrogen (PRATUSH). Detecting this faint signal requires precise calibration of the antenna, which includes a high-precision measurement of its reflection coefficient. This measurement must be performed \textit{in situ} at the observation site, as the antenna characteristics vary significantly with the environment. PRATUSH, a space-based radiometer, faces the additional challenge of structural distortions induced by thermal cycling, necessitating multiple measurements of the reflection coefficient. This work highlights the development of an \textit{in situ} Vector Network Analyser, which utilises a novel noise source-based calibration scheme and a cross-correlation spectrometer to perform magnitude and phase measurements of the complex reflection coefficient of the antenna. Further, we demonstrate the performance of the designed network analyser using independent measurements from a precision network analyser and reflection coefficient modelling of the device under test. We find the level of non-smooth calibration systematics, which need critical control for 21-cm signal detection, to be about $10^{-5}$. Finally, we study the impact of reflection coefficient correction on sky measurements, highlighting its usability for precision 21-cm observations.
The standard cosmological model is challenged by an ever-growing collection of observations, which invites (and stimulates) inquiry into possible additions and/or alterations. One such alteration comes from letting cosmic isotropy -- as demanded by the cosmological principle -- go, whilst maintaining only homogeneity. This study concerns Bianchi models, a class of anisotropic, homogeneous spacetimes, and in particular their perturbations. Knowledge of their properties under perturbations (such as allowed wavemodes) aids in understanding cosmological signatures of such universes, e.g. CMBs, and thus allows for comparsion to observation and the theory of the standard model. This study develops linear perturbation theory of general Bianchi models, by working in a frame such that metric components depend solely on (cosmic) time. Perturbation equations in the Newtonian gauge, but for arbitrary metric, are derived for energy density $ρ$, (relativistic) pressure $p$, momentum density $q$, and anisotropic stress $π$, for the case of scalar and pure tensor perturbations. For the former, the equations for density and pressure are combined to yield the equivalent of the Mukhanov-Sasaki equation for Bianchi models. For a specific choice of metric and fluid flow $u$, the Friedmann equations for Bianchi models are also formulated, as this knowledge is necessary to fully formulate the perturbation equations. Finally, the obtained results are applied to the formulations of density contrasts in an Einstein-de Sitter universe and a Bianchi I universe.