KICC papers

arXiv:2507.13905v1 Announce Type: new Abstract: We systematically extend the framework of galaxy bias renormalization to two-loop order. For the minimal complete basis of 29 deterministic bias operators up to fifth order in the density field and at leading order in gradient expansion we explicitly work out one- and two-loop renormalization. The latter is provided in terms of double-hard limits of bias kernels, which we find to depend on only one function of the ratio of the loop momenta. After including stochasticity in terms of composite operator renormalization, we apply the framework to the two-loop power spectrum of biased tracers and provide a simple result suitable for numerical evaluation. In addition, we work out one- and two-loop renormalization group equations (RGE) for deterministic bias coefficients related to bias operators constructed from a smoothed density field, generalizing previous works. We identify a linear combination of bias operators with enhanced UV sensitivity, related to a positive eigenvalue of the RGE. Finally, we present an analogy with the RGE as used in quantum field theory, suggesting that a resummation of large logarithms as employed in the latter may also yield useful applications in the study of large-scale galaxy bias.

arXiv:2504.01067v2 Announce Type: replace Abstract: We present deep JWST/NIRCam and MIRI imaging of Ion1, a previously confirmed Lyman Continuum (LyC)-emitting galaxy at $z_{spec}=3.794$. Together with existing HST imaging, these new observations from the JADES program enable a joint analysis of Ion1's LyC, rest-frame UV, stellar, and dust emission with unprecedented detail. We report the first detection of dust emission at rest-frame $\sim3 \mu$m in a high-redshift LyC-emitting galaxy using MIRI/F1500W. Our analysis suggests a porous distribution of dust in Ion1, with regions exhibiting evidence of dust deficit coinciding both with LyC-emitting regions and with the peak of H$\alpha$ emission. Furthermore, multi-band NIRCam imaging reveals a strong FUV-to-optical color gradient, where LyC-emitting regions appear significantly bluer than the rest of Ion1. Spatially resolved SED fitting confirms that this color gradient is primarily driven by spatially varying dust attenuation. Together, these findings suggest that Ion1's LyC emission originates from a compact star-forming complex near its stellar-light centroid, where stellar feedback carves out low HI column density channels, facilitating LyC escape. However, only a fraction of these LyC photons - specifically those along sightlines with minimal HI obscuration - ultimately escape and reach observers. This work underscores the critical role of dust and neutral gas geometry in shaping LyC escape in galaxies at high redshifts. Anisotropic LyC escape may be a common feature in the early Universe, which must be properly incorporated to constrain the Epoch of Reionization.

arXiv:2507.14136v1 Announce Type: new Abstract: We present new constraints on the halo masses and matter density profiles of DESI galaxy groups by cross-correlating samples of Luminous Red Galaxies (LRGs) and Bright Galaxy Survey (BGS) galaxies with the publicly available CMB lensing convergence map from ACT DR6. This provides an independent, lensing-based calibration of halo masses, complementary to methods relying on clustering or dynamics. We derive constraints on the mean halo mass for three DESI-selected samples, finding $\log(M_{\rm halo}/(M_\odot/h)) \approx 13.18$, 13.03 and 13.02 for the Main LRG, Extended LRG, and BGS samples, respectively. Using a halo model approach, we also compare the projected galaxy-matter density profiles with previously reported gas profiles inferred from measurements of the kinematic Sunyaev-Zel'dovich (kSZ) effect. This work addresses one of the key uncertainties in interpreting kSZ signals -- the unknown host halo mass distribution -- by providing an independent and consistent mass calibration. The agreement between the gas and total mass profiles at large aperture suggests that sufficiently far from the group center (2--3 virial radii), we recover all the baryons, offering a resolution to the 'missing baryon' problem. We further study the cumulative gas fractions for all galaxies as well as for the most massive galaxy groups in the sample ($\log(M_{\rm halo}/(M_\odot/h)) \approx 13.5$), finding values that are physically sensible and in agreement with previous findings using kSZ and X-ray data: compared to the TNG300 simulation, the observed gas fractions are systematically lower at fixed radius by $\gtrsim$4$\sigma$, providing compelling, independent evidence for stronger baryonic feedback in the real Universe. These findings highlight the power of combining CMB lensing with galaxy surveys to probe the interplay between baryons and dark matter in group-sized halos.

arXiv:2507.13495v1 Announce Type: new Abstract: Simulation-Based Inference (SBI) offers a principled and flexible framework for conducting Bayesian inference in any situation where forward simulations are feasible. However, validating the accuracy and reliability of the inferred posteriors remains a persistent challenge. In this work, we point out a simple diagnostic approach rooted in ensemble learning methods to assess the internal consistency of SBI outputs that does not require access to the true posterior. By training multiple neural estimators under identical conditions and evaluating their pairwise Kullback-Leibler (KL) divergences, we define a consistency criterion that quantifies agreement across the ensemble. We highlight two core use cases for this framework: a) for generating a robust estimate of the systematic uncertainty in parameter reconstruction associated with the training procedure, and b) for detecting possible model misspecification when using trained estimators on real data. We also demonstrate the relationship between significant KL divergences and issues such as insufficient convergence due to, e.g., too low a simulation budget, or intrinsic variance in the training process. Overall, this ensemble-based diagnostic framework provides a lightweight, scalable, and model-agnostic tool for enhancing the trustworthiness of SBI in scientific applications.

arXiv:2507.11530v2 Announce Type: new Abstract: We present direct measurements of the intrinsic alignments (IA) of over 2 million spectroscopic galaxies using DESI Data Release 1 and imaging from four lensing surveys: DES, HSC, KiDS, and SDSS. In this uniquely data-rich regime, we take initial steps towards a more tailored IA modelling approach by building a library of IA measurements across colour, luminosity, stellar mass, and redshift. We map the dependence between galaxy type -- in terms of rest-frame colour, strength of the 4000 Angstrom break, and specific star formation rate -- and IA amplitude; the bluest galaxies have an alignment consistent with zero, across low ($0.05

arXiv:2507.11530v1 Announce Type: new Abstract: We present direct measurements of the intrinsic alignments (IA) of over 2 million spectroscopic galaxies using DESI Data Release 1 and imaging from four lensing surveys: DES, HSC, KiDS, and SDSS. In this uniquely data-rich regime, we take initial steps towards a more tailored IA modelling approach by building a library of IA measurements across colour, luminosity, stellar mass, and redshift. We map the dependence between galaxy type -- in terms of rest-frame colour, strength of the 4000 Angstrom break, and specific star formation rate -- and IA amplitude; the bluest galaxies have an alignment consistent with zero, across low (0.05

arXiv:2507.09228v1 Announce Type: new Abstract: Constraints on the cosmological parameters of Torsion Condensation (TorC) are investigated using Planck 2018 Cosmic Microwave Background data. TorC is a case of Poincar\'e gauge theory -- a formulation of gravity motivated by the gauge field theories underlying fundamental forces in the standard model of particle physics. Unlike general relativity, TorC incorporates intrinsic torsion degrees of freedom while maintaining second-order field equations. At specific parameter values, it reduces to the $\Lambda$CDM model, providing a natural extension to standard cosmology. The base model of TorC introduces two parameters beyond those in $\Lambda$CDM: the initial value of the torsion scalar field and its time derivative -- one can absorb the latter by allowing the dark energy density to float. To constrain these parameters, `PolyChord` nested sampling algorithm is employed, interfaced via `Cobaya` with a modified version of `CAMB`. Our results indicate that TorC allows for a larger inferred Hubble constant, offering a potential resolution to the Hubble tension. Tension analysis using the $R$-statistic shows that TorC alleviates the statistical tension between the Planck 2018 and SH0Es 2020 datasets, though this improvement is not sufficient to decisively favour TorC over $\Lambda$CDM in a Bayesian model comparison. This study highlights TorC as a compelling theory of gravity, demonstrating its potential to address cosmological tensions and motivating further investigations of extended theories of gravity within a cosmological context. As current and upcoming surveys -- including Euclid, Roman Space Telescope, Vera C. Rubin Observatory, LISA, and Simons Observatory -- deliver data on gravity across all scales, they will offer critical tests of gravity models like TorC, making the present a pivotal moment for exploring extended theories of gravity.

arXiv:2507.08929v1 Announce Type: new Abstract: The majority of Little Red Dots (LRDs) hosting Active Galactic Nuclei (AGN) exhibits broad H$\alpha$ emission, which recent studies propose originates from scattering off free electrons within an ionized and dense medium embedding the Broad Line Region (BLR), rather than directly from the BLR itself. This model suggests that the observed broad lines may be intrinsically narrower than observed, which would lead to black hole masses that are up to two orders of magnitude smaller than what inferred when assuming that the whole broad line comes from the BLR. To test this model, we present a joint analysis of multiple hydrogen recombination lines in the ''Rosetta Stone''AGN, the brightest known LRD at $z$=2.26. We show that H$\alpha$, H$\beta$ and Pa$\beta$ have different spectral profiles, which is inconsistent with the predictions of the simple electron scattering scenario. Additionally, we test a variety of exponential models and show that none of them can simultaneously reproduce all three line profiles with physically plausible parameters. The inadequacy of these models for the Rosetta Stone implies that the scenario of electron scattering by an ionized medium surrounding the BLR is not universally applicable to LRDs and AGN, and therefore provides a counterexample to the claim of a universal and systematic overestimation of black hole masses.

arXiv:2507.08787v1 Announce Type: new Abstract: We present an analysis of metallicities and chemical abundances at $3-0.6$) without the need for exotic yields in our chemical network. Instead, bursty star formation naturally generates order-of-magnitude excursions in N/O on $\lesssim$100 Myr timescales due to temporally differential galactic winds; after a starburst, stellar feedback expels gas, leaving a large population of asymptotic-giant-branch stars to dominate the enrichment of the relatively low-mass interstellar medium. NRGs lie below the main sequence and typically exhibit $\mathrm{EW}[H$\beta$]\lesssim40$ \r{A}, in apparent tension with observed high-EW NRGs. This tension is reconciled if observed NRGs are in the initial stages of a subsequent starburst, illuminating previously enriched gas, which is supported by the finding of high SFR surface density nitrogen-rich giant molecular clouds.

arXiv:2507.07346v2 Announce Type: replace Abstract: These notes are very much work-in-progress and simply intended to showcase, in various degrees of details (and rigour), some of the cosmology calculations that class_sz can do. We describe the class_sz code in C, Python and Jax. Based on the Boltzmann code class, it can compute a wide range of observables relevant to current and forthcoming CMB and Large Scale Structure surveys. This includes galaxy shear and clustering, CMB lensing, thermal and kinetic Sunyaev and Zeldovich observables, Cosmic Infrared Background, cross-correlations and three-point statistics. Calculations can be done either within the halo model or the linear bias model. For standard $\Lambda$CDM cosmology and extensions, class_sz uses high-accuracy cosmopower emulators of the CMB and matter power spectrum to accelerate calculations. With this, along with efficient numerical integration routines, most class_sz output can be obtained in less than 500 ms (CMB $C_\ell$'s or matter $P(k)$ take $\mathcal{O}(1\mathrm{ms})$), allowing for fast or ultra-fast parameter inference analyses. Parts of the calculations are "jaxified", so the software can be integrated into differentiable pipelines.

arXiv:2507.07346v1 Announce Type: new Abstract: These notes are very much work-in-progress and simply intended to showcase, in various degrees of details (and rigour), some of the cosmology calculations that class_sz can do. We describe the class_sz code in C, Python and Jax. Based on the Boltzmann code class, it can compute a wide range of observables relevant to current and forthcoming CMB and Large Scale Structure surveys. This includes galaxy shear and clustering, CMB lensing, thermal and kinetic Sunyaev and Zeldovich observables, Cosmic Infrared Background, cross-correlations and three-point statistics. Calculations can be done either within the halo model or the linear bias model. For standard $\Lambda$CDM cosmology and extensions, class_sz uses high-accuracy cosmopower emulators of the CMB and matter power spectrum to accelerate calculations. With this, along with efficient numerical integration routines, most class_sz output can be obtained in less than 500 ms (CMB $C_\ell$'s or matter $P(k)$ take $\mathcal{O}(1\mathrm{ms})$), allowing for fast or ultra-fast parameter inference analyses. Parts of the calculations are "jaxified", so the software can be integrated into differentiable pipelines.

arXiv:2507.07155v1 Announce Type: new Abstract: We evaluate 9 Retrieval Augmented Generation (RAG) agent configurations on 105 Cosmology Question-Answer (QA) pairs that we built specifically for this purpose.The RAG configurations are manually evaluated by a human expert, that is, a total of 945 generated answers were assessed. We find that currently the best RAG agent configuration is with OpenAI embedding and generative model, yielding 91.4\% accuracy. Using our human evaluation results we calibrate LLM-as-a-Judge (LLMaaJ) system which can be used as a robust proxy for human evaluation. These results allow us to systematically select the best RAG agent configuration for multi-agent system for autonomous scientific discovery in astrophysics (e.g., cmbagent presented in a companion paper) and provide us with an LLMaaJ system that can be scaled to thousands of cosmology QA pairs. We make our QA dataset, human evaluation results, RAG pipelines, and LLMaaJ system publicly available for further use by the astrophysics community.

arXiv:2405.05772v3 Announce Type: replace Abstract: We report the detection of continuum excess in the rest-frame UV between 3000 {\AA} and 3550 {\AA} in the JWST/NIRSpec spectrum of GN-z11, a luminous galaxy $z=10.603$. The shape of the continuum excess resembles a Balmer continuum but has a break around 3546 {\AA}. The fitting result of this excess depends on the assumed origin of the continuum. If the continuum of GN-z11 is dominated by a stellar population with a small Balmer break, the apparent blueshift of the Balmer continuum is not significant and the best-fit Balmer continuum model indicates a temperature of $T_e = 1.78^{+0.25}_{-0.21}\times 10^4$ K. In contrast, if the continuum is dominated by AGN emission, a nebular continuum model cannot fit the spectrum properly. The absence of the Balmer jump indicates an electron temperature of $\sim 3\times 10^4$ K, significantly higher than the temperature of $T_{e}({\rm O^{2+}}) = 1.36\pm 0.13\times 10^{4}$ K inferred from [OIII]$\lambda 4363$ and [OIII]$\lambda 5007$. The temperature difference can result from mixing of different ionized regions: the Balmer emission mainly arises from dense and hot clouds in the Broad Line Region, whereas the forbidden lines originate from less dense and colder gas. An alternative explanation for the observed continuum excess is the FeII emission, which shows a characteristic jump blueward of the Balmer limit as previously seen in the spectra of many lower-redshift quasars. Through comparisons with Cloudy models, we show an Fe abundance above $\sim 1/3$ solar is likely needed, which could be achieved via enrichment from Type-Ia supernovae, hypernovae, or pair-instability supernovae.

arXiv:2507.00351v1 Announce Type: new Abstract: We present Dark from Light (DfL) - a novel method to infer the dark sector in wide-field galaxy surveys, leveraging a machine learning approach trained on contemporary cosmological simulations. The aim of this algorithm is to provide a fast, straightforward, and accurate route to estimating dark matter halo masses and group membership in wide-field spectroscopic galaxy surveys. This approach requires a highly limited number of input parameters and yields full probability distribution functions for the output halo masses. To achieve this, we train a series of Random Forest (RF) regression models on the IllustrisTNG and EAGLE simulations at z=0-3, which provide model-dependent mappings from luminous tracers to dark matter halo properties. We incorporate the individual regression models into a virial group-finding algorithm (DfL), which outputs halo properties for observational-like input data. We test the method at z=0-2 for both the EAGLE and IllustrisTNG models, as well as in a cross-validation mode. We demonstrate that known halo masses can be recovered with a mean systematic bias of $\langle b \rangle = \pm 0.10\,$dex (resulting from simulation choice), a mean statistical uncertainty of $\langle \sigma \rangle = 0.12 \,$dex across epochs, and a central - (core) satellite classification accuracy of 96%. We establish that this approach yields superior halo mass recovery to standard abundance matching applied to groups identified through a friends-of-friends algorithm. Additionally, we compare the outputs of DfL to observational constraints on the $M_* - M_{\rm Halo}$ relation from strong gravitational lensing at $z \sim 0$, demonstrating the promise of this novel approach. Finally, we systematically quantify how DfL performs on observational-like input data with varying stellar mass uncertainty and spectroscopic incompleteness, enabling robust error calibration.

arXiv:2506.22999v1 Announce Type: new Abstract: The recent measurement of baryonic acoustic oscillations (BAO) by the Dark Energy Spectroscopic Instrument (DESI) reveals a mild tension with observations of the cosmic microwave background (CMB) within the standard $\Lambda$CDM cosmological model. This discrepancy leads to a preference for a total neutrino mass that is lower than the minimum value inferred from neutrino oscillation experiments. Alternatively, this tension can be eased within $\Lambda$CDM by assuming a higher optical depth ($\tau \simeq 0.09$), but such a value conflicts with large-scale CMB polarization data. We point out that cosmic birefringence, as suggested by recent Planck reanalyses, resolves this discrepancy if the birefringence angle varies significantly during reionization. Specifically, we consider the fact that the measured cosmic birefringence angle $\beta_0=0.34\pm0.09\,(1\,\sigma)\,$deg has the phase ambiguity, i.e., the measured rotation angle is described by $\beta=\beta_0+180n\,$deg ($n\in \mathbb{Z}$). We show that cosmic birefringence induced by axion-like particles with nonzero $n$ suppresses the reionization bump, allowing a higher $\tau$ consistent with data. We provide a viable parameter space where the birefringence effect simultaneously accounts for the low-$\ell$ polarization spectra, the Planck $EB$ correlations, and the elevated value of $\tau$, suggesting a key role for cosmic birefringence in current cosmological tensions.

arXiv:2401.10328v2 Announce Type: replace Abstract: The spectra of high-redshift ($z\gtrsim 6$) quasars contain valuable information on the progression of the Epoch of Reionization (EoR). At redshifts $z<6$, the observed Lyman-series forest shows that the intergalactic medium (IGM) is nearly ionized, while at $z>7$ the observed quasar damping wings indicate high neutral gas fractions. However, there remains a gap in neutral gas fraction constraints at $6\lesssim z \lesssim 7$ where the Lyman series forest becomes saturated but damping wings have yet to fully emerge. In this work, we use a sample of 18 quasar spectra at redshifts $6.0

arXiv:2506.22147v1 Announce Type: new Abstract: JWST has revealed a large population of active galactic nuclei (AGN) in the distant universe, which are challenging our understanding of early massive black hole seeding and growth. We expand the exploration of this population to lower luminosities by stacking $\sim 600$ NIRSpec grating spectra from the JWST Advanced Deep Extragalactic Survey (JADES) at $3

arXiv:2410.14404v3 Announce Type: replace Abstract: According to CMB measurements, baryonic matter constitutes about $5\%$ of the mass-energy density of the universe. A significant population of these baryons, for a long time referred to as `missing', resides in a low density, warm-hot intergalactic medium (WHIM) outside galaxy clusters, tracing the ``cosmic web'', a network of large scale dark matter filaments. Various studies have detected this inter-cluster gas, both by stacking and by observing individual filaments in compact, massive systems. In this paper, we study short filaments (< 10 Mpc) connecting massive clusters ($M_{500} \approx 3\times 10^{14} M_{\odot}$) detected by the Atacama Cosmology Telescope (ACT) using the scattering of CMB light off the ionised gas, a phenomenon known as the thermal Sunyaev-Zeldovich (tSZ) effect. The first part of this work is a search for suitable candidates for high resolution follow-up tSZ observations. We identify four cluster pairs with an intercluster signal above the noise floor (S/N $>$ 2), including two with a tentative $>2\sigma$ statistical significance for an intercluster bridge from the ACT data alone. In the second part of this work, starting from the same cluster sample, we directly stack on ${\sim}100$ cluster pairs and observe an excess SZ signal between the stacked clusters of $y=(7.2^{+2.3}_{-2.5})\times 10^{-7}$ with a significance of $3.3\sigma$. It is the first tSZ measurement of hot gas between clusters in this range of masses at moderate redshift ($\langle z\rangle\approx 0.5$). We compare this to the signal from simulated cluster pairs with similar redshifts and separations in the THE300 and MAGNETICUM Pathfinder cosmological simulations and find broad consistency. Additionally, we show that our measurement is consistent with scaling relations between filament parameters and mass of the embedded halos identified in simulations.

arXiv:2503.02927v2 Announce Type: replace Abstract: We investigate the impact of ionizing external ultraviolet (UV) radiation on low-mass haloes ($M_{h}<10^{10}M_\odot$) at high redshift using $1140M_\odot$ baryonic resolution zoom-in simulations of seven regions from the THESAN-ZOOM project. We compare three simulation sets that differ in the treatment of external UV radiation: one employing a uniform UV background initiated at z=10.6 in addition to radiation transport for local sources, another with the same background starting at z=5.5, and the default configuration in which the large-scale radiation field from the parent THESAN-1 simulation box acts as a boundary condition. The multi-phase interstellar medium (ISM) model, combined with its high mass resolution, allows us to resolve all star-forming haloes and capture the back-reaction of ionizing radiation on galaxy properties during the epoch of reionization. When present, external UV radiation efficiently unbinds gas in haloes with masses below $10^9M_\odot$ and suppresses subsequent star formation. As a result, in simulations with early reionization, minihaloes fail to form stars from pristine gas, leading to reduced metal enrichment of gas later accreted by more massive haloes. Consequently, haloes with masses below $10^{10}M_\odot$ at all simulated epochs (z>3) exhibit lower metallicities and altered metallicity distributions. The more accurate and realistic shielding from external UV radiation, achieved through self-consistent radiative transfer, permits the existence of a cold but low-density gas phase down to z=3. These findings highlight the importance of capturing a patchy reionization history in high-resolution simulations targeting high-redshift galaxy formation. We conclude that at minimum, a semi-numerical model that incorporates spatially inhomogeneous reionization and a non-uniform metallicity floor is necessary to accurately emulate metal enrichment in minihaloes.

arXiv:2506.20824v1 Announce Type: new Abstract: We present constraints on isotropic cosmic birefringence induced by axion-like particles (ALPs), derived from the analysis of cosmic microwave background (CMB) polarization measurements obtained with the high-frequency channels of Planck. Recent measurements report a hint of isotropic cosmic birefringence, though its origin remains uncertain. The detailed dynamics of ALPs can leave characteristic imprints on the shape of the $EB$ angular power spectrum, which can be exploited to constrain specific models of cosmic birefringence. We first construct a multi-frequency likelihood that incorporates an intrinsic nonzero $EB$ power spectrum. We also show that the likelihood used in previous studies can be further simplified without loss of generality. Using this framework, we simultaneously constrain the ALP model parameters, the instrumental miscalibration angle, and the amplitudes of the $EB$ power spectrum of a Galactic dust foreground model. We find that, if ALPs are responsible for the observed cosmic birefringence, ALP masses at $\log_{10}m_{\phi}[{\rm eV}]\simeq-27.8$, $-27.5$, $-27.3$, $-27.2$, $-27.1$, as well as $\log_{10}m_{\phi}[{\rm eV}]\in[-27.0,-26.5]$, are excluded at more than $2\,\sigma$ statistical significance.