skip to content

Kavli Institute for Cosmology, Cambridge

 

Updated dark pixel fraction constraints on reionization's end from the Lyman-series forests of XQR-30

Fri, 31/10/2025 - 08:17
arXiv:2510.25829v1 Announce Type: new Abstract: The fraction of "dark pixels" in the Ly$\alpha$ and other Lyman-series forests at $z\sim 5-6$ provides a powerful constraint on the end of the reionization process. Any spectral region showing transmission must be highly ionized, while dark regions could be ionized or neutral, thus the dark pixel fraction provides a (nearly) model independent upper limit to the volume-filling fraction of the neutral intergalactic medium, modulo choices in binning scale and dark pixel definition. Here we provide updated measurements of the 3.3 comoving Mpc dark pixel fraction at $z=4.85-6.25$ in the Ly$\alpha$, Ly$\beta$, and Ly$\gamma$ forests of 34 deep $5.8\lesssim z\lesssim 6.6$ quasar spectra from the (enlarged) XQR-30 sample. Using the negative pixel method to measure the dark pixel fraction, we derive fiducial $1\sigma$ upper limits on the volume-average neutral hydrogen fraction of $\langle x_{\rm{HI}}\rangle \leq \{0.030 + 0.048, 0.095 + 0.037, 0.191 + 0.056, 0.199 + 0.087\}$ at $\bar{z} = \{5.481, 5.654, 5.831, 6.043\}$ from the optimally sensitive combination of the Ly$\beta$ and Ly$\gamma$ forests. We further demonstrate an alternative method that treats the forest flux as a mixture of dark and transparent regions, where the latter are modeled using a physically-motivated parametric form for the intrinsic opacity distribution. The resulting model-dependent upper limits on $\langle x_{\rm{HI}}\rangle$ are similar to those derived from our fiducial model-independent analysis. We confirm that the bulk of reionization must be finished at $z>6$, while leaving room for an extended "soft landing" to the reionization history down to $z\sim 5.4$ suggested by Ly$\alpha$ forest opacity fluctuations.

Learning the Universe: Cosmological and Astrophysical Parameter Inference with Galaxy Luminosity Functions and Colours

Thu, 30/10/2025 - 10:29
arXiv:2411.13960v2 Announce Type: replace Abstract: We perform the first direct cosmological and astrophysical parameter inference from the combination of galaxy luminosity functions and colours using a simulation based inference approach. Using the Synthesizer code we simulate the dust attenuated ultraviolet-near infrared stellar emission from galaxies in thousands of cosmological hydrodynamic simulations from the CAMELS suite, including the Swift-EAGLE, IllustrisTNG, Simba & Astrid galaxy formation models. For each galaxy we calculate the rest-frame luminosity in a number of photometric bands, including the SDSS $\textit{ugriz}$ and GALEX FUV & NUV filters; this dataset represents the largest catalogue of synthetic photometry based on hydrodynamic galaxy formation simulations produced to date, totalling >200 million sources. From these we compile luminosity functions and colour distributions, and find clear dependencies on both cosmology and feedback. We then perform simulation based (likelihood-free) inference using these distributions to constrain $\Omega_m$, $\sigma_8$, and four parameters controlling the strength of stellar and AGN feedback. Both colour distributions and luminosity functions provide complementary information on certain parameters when performing inference. We achieve constraints on the stellar feedback parameters, as well as $\Omega_m$ and $\sigma_8$. The latter is attributable to the fact that the photometry encodes the star formation-metal enrichment history of each galaxy; galaxies in a universe with a higher $\sigma_8$ tend to form earlier and have higher metallicities, which leads to redder colours. We find that a model trained on one galaxy formation simulation generalises poorly when applied to another, and attribute this to differences in the subgrid prescriptions, and lack of flexibility in our emission modelling. The photometric catalogues are publicly available at: https://camels.readthedocs.io/

A great diversity of spectral shapes in the ionising spectra of z ~ 0.6 - 1 galaxies revealed by HST/COS and possible detection of nebular LyC emission

Tue, 28/10/2025 - 09:44
arXiv:2510.22152v1 Announce Type: new Abstract: We present observations of eleven compact star-forming galaxies in the redshift range z = 0.6145 - 1.0053, with the Cosmic Origins Spectrograph (COS) on board the Hubble Space Telescope (HST). We aim to spectroscopically measure for the first time the Lyman continuum (LyC) over a wider rest-frame wavelength range of ~ 600 - 900A compared to ~ 850 - 900A in previous studies of galaxies at z ~ 0.3 - 0.4. The HST data are supplemented by SDSS spectra of all galaxies and by a VLT/Xshooter spectrum of one galaxy, J0232+0025. These data are used to derive the spectral energy distribution in the entire UV and optical range, the stellar mass, and the chemical composition from the nebular emission lines. We detect stellar LyC emission in seven out of eleven galaxies with escape fractions, f_esc(LyC), in the range of ~ 2 - 60%, and establish upper limits for f_esc(LyC) in the remaining galaxies. We discover for the first time nebular LyC emission as a bump just bluewards of the LyC limit at 912A in two galaxies, J0232+0025 and J1021+0436. We find a similar bump among our earlier studies in a less distant galaxy J1243+4646 with z = 0.4317. We conclude that the use of the LyC continuum in the wavelength range close to the LyC limit, which contains both the stellar and nebular continua, requires special consideration to not overestimate the observed f_esc(LyC).

$\texttt{SBi3PCF:}$ Simulation-based inference with the integrated 3PCF

Mon, 27/10/2025 - 11:10
arXiv:2510.13805v2 Announce Type: replace Abstract: We present $\texttt{SBi3PCF}$, a simulation-based inference (SBI) framework for analysing a higher-order weak lensing statistic, the integrated 3-point correlation function (i3PCF). Our approach forward-models the cosmic shear field using the $\texttt{CosmoGridV1}$ suite of N-body simulations, including a comprehensive set of systematic effects such as intrinsic alignment, baryonic feedback, photometric redshift uncertainty, shear calibration bias, and shape noise. Using this, we have produced a set of DES Y3-like synthetic measurements for 2-point shear correlation functions $\xi_{\pm}$ (2PCFs) and i3PCFs $\zeta_{\pm}$ across 6 cosmological and 11 systematic parameters. Having validated these measurements against theoretical predictions and thoroughly examined for potential systematic biases, we have found that the impact of source galaxy clustering and reduced shear on the i3PCF is negligible for Stage-III surveys. Furthermore, we have tested the Gaussianity assumption for the likelihood of our data vector and found that while the sampling distribution of the 2PCF can be well approximated by a Gaussian function, the likelihood of the combined 2PCF + i3PCF data vector including filter sizes of $90'$ and larger can deviate from this assumption. Our SBI pipeline employs masked autoregressive flows to perform neural likelihood estimation and is validated to give statistically accurate posterior estimates. On mock data, we find that including the i3PCF yields a substantial $63.8\%$ median improvement in the figure of merit for $\Omega_m - \sigma_8 - w_0$. These findings are consistent with previous works on the i3PCF and demonstrate that our SBI framework can achieve the accuracy and realism needed to analyse the i3PCF in wide-area weak lensing surveys.

Tentative rotation in a galaxy at z$\sim$14 with ALMA

Mon, 27/10/2025 - 11:05
arXiv:2503.10751v3 Announce Type: replace Abstract: We re-analysed ALMA observations of the [OIII]$\lambda$88$\mu$m emission line in JADES-GS-z14.0, so far the most distant spectroscopically confirmed galaxy at z=14.18. Our analysis shows a tentative detection of a velocity gradient of [OIII]$\lambda$88$\mu$m using three independent tests: 1) construction of moment maps; 2) extraction of integrated spectra from a grid of apertures; and 3) spectro-astrometry in both the image and uv planes. We performed kinematical fitting using the KinMS code and estimated a dynamical mass of log$_{10}$(M$_{\rm dyn}$/$\rm M_\odot$)= 9.4$^{+0.8}_{-0.4}$, with the bulk of the uncertainties due to the degeneracy between dynamical mass and inclination. We measure an upper limit on the velocity dispersion ($\sigma_{v}$) of $<40~$ km/s~which results in an estimate of V$_{\rm rot}/\sigma>$ 2.5. This result, if confirmed with higher-resolution observations, would imply that kinematically cold discs are already in place at $z\sim14$. Comparison with mock observations from the SERRA cosmological simulations confirms that even low-resolution observations are capable of detecting a velocity gradient in $z>10$ galaxies as compact as JADES-GS-z14.0. This work shows that deeper ALMA or JWST/NIRSpec IFS observations with high spatial resolution will be able to estimate an accurate dynamical mass for JADES-GS-z14.0, providing an upper limit to the stellar mass of this over-luminous galaxy.

Net-zero gas inflow: Measurement of the gas consumption history of a massive quiescent galaxy

Mon, 27/10/2025 - 11:04
arXiv:2405.19401v2 Announce Type: replace Abstract: JWST is discovering increasing numbers of quiescent galaxies 1--2 billion years after the Big Bang, whose redshift, high mass, and old stellar ages indicate that their formation and quenching were surprisingly rapid. This fast-paced evolution seems to require that feedback from AGN (active galactic nuclei) be faster and/or more efficient than previously expected \citep{Xie24}. We present deep ALMA observations of cold molecular gas (the fuel for star formation) in a massive, fast-rotating, post-starburst galaxy at $z=3.064$. This galaxy hosts an AGN, driving neutral-gas outflows with a mass-outflow rate of $60\pm20$ M$_{\odot}$ yr$^{-1}$, and has a star-formation rate of $<5.6$ M$_{\odot}$ yr$^{-1}$. Our data reveal this system to be the most distant gas-poor galaxy confirmed with direct CO observations (molecular-gas mass $< 10^{9.1}$ M$_{\odot}$; <0.8 % of its stellar mass). Combining ALMA and JWST observations, we estimate the gas-consumption history of this galaxy, showing that it evolved with net zero gas inflow, i.e., gas consumption by star formation matches the amount of gas this galaxy is missing relative to star-forming galaxies. This could arise both from preventive feedback stopping further gas inflow, which would otherwise refuel star formation or, alternatively, from fine-tuned ejective feedback matching precisely gas inflows. Our methods, applied to a larger sample, promise to disentangle ejective vs preventive feedback.

JADES: Carbon-enhanced, Nitrogen-normal compact galaxy at z=11.2

Mon, 27/10/2025 - 10:47
arXiv:2507.17809v3 Announce Type: replace Abstract: Over the past few years \textit{JWST} has been a major workhorse in detecting and constraining the metal enrichment of the first galaxies in the early Universe and finding the source of the ionisation of their interstellar medium. In this work, we present new deep JWST/NIRSpec spectroscopy of GS-z11-1, a galaxy at z = 11.28, in which we report the detection of multiple rest-frame UV and optical emission lines: CIII]$\lambda\lambda$1907,09, CIV]$\lambda\lambda$1548,51, [OII]$\lambda\lambda$3726,29, [NeIII]$\lambda$3869, H$\gamma$ and tentative evidence for HeII$\lambda$1640. The ionisation properties of GS-z11-1 are consistent with star formation, with potential contribution from an active galactic nucleus (AGN). We estimate a galaxy stellar mass of log(M$_{*}$/M$_{\odot}$) = 7.8$\pm$0.2 and log(SFR/(M$_{\odot}$ yr$^{-1}$))= 0.32$\pm$0.11 for the fiducial SF-only models. We measured C/O from the SED modelling of C/O = 1.20$\pm0.15 \times$ solar. This is one of the highest C/O abundances at z$>$10, and it is consistent with either PopII and PopIII enrichment paths. Despite this source being extremely compact, with a half-light radius of 73$\pm$10 pc, we see no increased equivalent width of NIV] and NIII] emission lines as seen in some other compact sources at similar redshifts, a potential signature of second-generation stars in GCs. Overall, this galaxy exhibits low metallicity and high ionisation parameter consistent with intense star-formation or AGN activity in the early Universe, possibly observed before the enrichment by the second generation of stars in proto-globular clusters in the core of the galaxy.

Characterizing CMB noise anisotropies from CMB delensing

Fri, 24/10/2025 - 14:54
arXiv:2504.20919v2 Announce Type: replace Abstract: Un-doing the effect of gravitational lensing on the Cosmic Microwave Background (`de-lensing') is essential in shaping constraints on weak signals limited by lensing effects on the CMB, for example on a background of primordial gravitational waves. Removing these anisotropies induced by large-scale structures from the CMB maps also generally helps our view of the primordial Universe by sharpening the acoustic peaks and the damping tail. However, practical implementations of delensing transfer parts of these anisotropies to the noise maps. This will induce a new large scale `mean-field' bias to any anisotropy estimator applied to the delensed CMB, and this bias directly traces large-scale structures. This paper analytically quantifies this delensed noise mean-field and its impact on quadratic (QE) and likelihood-based lensing estimators. We show that for Simons-Observatory-like surveys, this mean-field bias can reach 15\% in cross-correlation with large-scale structures if unaccounted for. We further demonstrate that this delensed noise mean-field can be safely neglected in likelihood-based estimators without compromising the quality of lensing reconstruction or $B$-mode delensing, provided the resulting lensing map is properly renormalized.

Classification of Radio Backgrounds at Cosmic Dawn and 21 cm Signal Confirmation Using Neural Networks

Fri, 24/10/2025 - 14:52
arXiv:2410.04792v3 Announce Type: replace Abstract: Several ongoing and upcoming radio telescopes aim to detect either the global 21 cm signal or the 21 cm power spectrum. The extragalactic radio background, as detected by ARCADE-2 and LWA-1, suggests a strong radio background from cosmic dawn, which can significantly alter the cosmological 21 cm signal, enhancing both the global signal amplitude and the 21 cm power spectrum. In this paper, we employ an artificial neural network (ANN) to check if there is a radio excess over the cosmic microwave background in mock data, and if present, we classify its type into one of two categories, a background from high-redshift radio galaxies or a uniform exotic background from the early Universe. Based on clean data (without observational noise), the ANN can predict the background radiation type with $96\%$ accuracy for the power spectrum and $90\%$ for the global signal. Although observational noise reduces the accuracy, the results remain quite useful. We also apply ANNs to map the relation between the 21 cm power spectrum and the global signal. By reconstructing the global signal using the 21 cm power spectrum, an ANN can estimate the global signal range consistent with an observed power spectrum from SKA-like experiments. Conversely, we show that an ANN can reconstruct the 21 cm power spectrum over a wide range of redshifts and wavenumbers given the global signal over the same redshifts. Such trained networks can potentially serve as a valuable tool for cross-confirmation of the 21 cm signal.

Rapid, out of equilibrium metal enrichment indicated by a flat mass-metallicity relation at z~6 from NIRCam grism spectroscopy

Fri, 24/10/2025 - 14:05
arXiv:2510.19959v1 Announce Type: new Abstract: We aim to characterise the mass-metallicity relation (MZR) and the 3D correlation between stellar mass, metallicity and star-formation rate (SFR) known as the fundamental metallicity relation (FMR) for galaxies at $5

JWST Observations of SN 2023ixf I: Completing the Early Multi-Wavelength Picture with Plateau-phase Spectroscopy

Thu, 23/10/2025 - 09:31
arXiv:2507.18785v2 Announce Type: replace Abstract: We present and analyze panchromatic (0.35--14 $\mu$m) spectroscopy of the Type II supernova 2023ixf, including near- and mid-infrared spectra obtained 33.6 days after explosion during the plateau-phase, with the James Webb Space Telescope (JWST). This is the first in a series of papers examining the evolution of SN 2023ixf with JWST spanning the initial 1000 days after explosion, monitoring the formation and growth of molecules and dust in ejecta and surrounding environment. The JWST infrared spectra are overwhelmingly dominated by H lines, whose profiles reveal ejecta structures, including flat tops, blue notches, and red shoulders, unseen in the optical spectra. We characterize the nature of these structures, concluding that they likely result from a combination of ejecta geometry, viewing angle, and opacity effects. We find no evidence for the formation of dust precursor molecules such as carbon-monoxide (CO), nor do we observe an infrared excess attributable to dust. These observations imply that the detections of molecules and dust in SN 2023ixf at later epochs arise either from freshly synthesized material within the ejecta or circumstellar material at radii not yet heated by the supernova at this epoch.

The impact of superradiance on the spin evolution of variably accreting massive black holes

Thu, 23/10/2025 - 09:25
arXiv:2510.19443v1 Announce Type: new Abstract: This paper explores how time-varying increases in mass accretion onto rapidly spinning black holes influence their long-term spin evolution when affected by superradiance - a process where energy is extracted from the black hole by a surrounding axion field. Using simulations the study tracks how sudden accretion boosts affect a critical spin-down phase (the superradiance drop) during which the black hole's spin rapidly decreases while its mass remains nearly constant. The black hole spin evolution is controlled by the competition between two processes: how fast angular momentum is added through accretion, and how fast it is removed by the axion cloud. One major conclusion is that boosts to the accretion rate before the superradiance drop have the strongest effect, as they can delay or reshape the drop and significantly shrink the region of the mass-spin plane depopulated due to the superradiance. In particular, a super-Eddington accretion rate of 5 times Eddington accretion, lasting for 4 Myr and occurring 30 Myr before the superradiance drop can reduce the superradiance exclusion region in the mass-spin plane by 40 percent. In contrast, boosts to the accretion rate after the superradiance drop only cause temporary changes in the black hole spin. The study also shows that black holes with lighter axion clouds are more sensitive to these early boosts and can show observable spin changes lasting tens to hundreds of millions of years. Heavier axion clouds, however, require much stronger or longer-lasting boosts to produce similar effects, making them more stable under variable accretion.

The slope and scatter of the star forming main sequence at z~5 : reconciling observations with simulations

Thu, 23/10/2025 - 09:12
arXiv:2510.19044v1 Announce Type: new Abstract: Galaxies exhibit a tight correlation between their star-formation rate and stellar mass over a wide redshift range known as the star-forming main sequence (SFMS). With JWST, we can now investigate the SFMS at high redshifts down to masses of $\sim10^6$ M$_{\odot}$, using sensitive star-formation rate tracers such as H$\alpha$ emission -- which allow us to probe the variability in star formation histories. We present inferences of the SFMS based on 316 H$\alpha$-selected galaxies at $z\sim4$-$5$ with $\log(\rm M_\star/M_\odot) = 6.4$ -$10.6$. These galaxies were identified behind the Abell 2744 lensing cluster with NIRCam grism spectroscopy from the ``All the Little Things'' (ALT) survey. At face value, our data suggest a shallow slope of the SFMS (SFR $\propto \mathrm{M}_\star^\alpha$, with $\alpha=0.45$). After correcting for the H$\alpha$-flux limited nature of our survey using a Bayesian framework, the slope steepens to $\alpha = 0.59^{+0.10}_{-0.09}$, whereas current data on their own are inconclusive on the mass dependence of the scatter. These slopes differ significantly from the slope of $\approx1$ expected from the observed evolution of the galaxy stellar mass function and from simulations. When fixing the slope to $\alpha=1$, we find evidence for a decreasing intrinsic scatter with stellar mass (from $\approx 0.5$ dex at M$_\star=10^8$ M$_\odot$ to $0.4$ dex at M$_\star=10^{10}$ M$_\odot$). This tension might be explained by a (combination of) luminosity-dependent SFR(H$\alpha$) calibration, a population of (mini)-quenched low-mass galaxies, or underestimated dust attenuation in high-mass galaxies. Future deep observations across facilities can quantify these processes, enabling better insights into the variability of star formation histories.

On the consistency of jet feedback modelling across different astrophysics hydrodynamical codes

Wed, 22/10/2025 - 10:57
arXiv:2510.18293v1 Announce Type: new Abstract: Active Galactic Nuclei (AGN) feedback is essential in cosmological simulations of galaxy formation, yet its implementation has to rely on subgrid models due to limited resolution. We present a novel subgrid jet-launching method for galaxy formation simulations and implement it in three hydrodynamical codes: the smoothed particle hydrodynamics (SPH) code SWIFT, the moving-mesh code AREPO, and the Eulerian grid code PLUTO. To isolate the impact of hydrodynamical solvers on jet evolution, we compare idealised jets and their remnants in uniform and stratified media across resolutions and jet parameters. In uniform media, all jets drive bow shocks, inflate hot lobes, exhibit backflows, and evolve self-similarly. For the parameters explored, SWIFT lobes are shorter, wider, and hotter; AREPO lobes are longer, thinner, and cooler; while PLUTO lobes display complex flows with intermediate characteristics. In stratified media, jets deviate from self-similar evolution, inflating longer and thinner lobes due to lower external ram pressure. After switch-off, SWIFT jets evolve into smooth cylindrical bubbles, AREPO jets produce long filamentary remnants, and PLUTO jets yield intermediate-length remnants with varying degrees of mixing. Despite such differences, all jets and remnants have a similar impact on the ambient medium. We conclude that variations in lobe properties between codes emerge even for identical subgrid prescriptions, since the coupling of jet feedback to resolvable scales and the effective resolution depend on the hydrodynamical method. In structure formation simulations, these solver differences are likely subdominant to uncertainties in subgrid modelling and calibration, while averaging over galaxy populations may lessen their impact.

The ALPINE-CRISTAL-JWST Survey: Stellar and nebular dust attenuation of main-sequence galaxies at z~4-6

Wed, 22/10/2025 - 10:57
arXiv:2510.18248v1 Announce Type: new Abstract: Characterizing dust attenuation is crucial for revealing the intrinsic physical properties of galaxies. We present an analysis of dust attenuation in 18 spectroscopically confirmed star-forming main-sequence galaxies at $z = 4.4-5.7$ observed with JWST/NIRSpec IFU and NIRCam, selected from the ALPINE and CRISTAL ALMA large programs. We fit the emission line fluxes from NIRSpec and the broad-band photometry from NIRCam with Prospector, using both spatially integrated emission and $\sim0.6$ kpc pixel-by-pixel measurements. We derive the stellar-to-nebular dust attenuation ratio ($f=E(B-V)_{\mathrm{star}}/E(B-V)_{\mathrm{neb}}$) from the SED fits and the Balmer decrement with H$\alpha$ and H$\beta$. Although individual galaxies show large scatter, the best-fit value is $f = 0.51^{+0.04}_{-0.03}$, slightly higher than that measured for local starburst galaxies. We find weak correlations of $f$ with galaxy properties, increasing with higher specific star-formation rates, younger stellar ages, and more recent star-formation. For the range of $E(B-V)_{\mathrm{star}} = 0.009-0.15$ mag for in our sample, assuming $f = 1$ (often adopted in high-redshift studies) instead of $f = 0.51$ underestimate line luminosities and ionizing photon production efficiency $\xi_\text{ion}$ by $\sim3-36\%$ and $\sim4-46\%$, respectively. We also find that the total stellar masses estimated from spatially-integrated SED fits with a delayed-$\tau$ star-formation histories are systematically smaller than the sum of pixel-by-pixel SED fits, with a median offset of $\sim 0.26$ dex, likely because the integrated fits are biased toward luminous young stellar populations.

The thermal and kinematic Sunyaev-Zeldovich effect in galaxy clusters and filaments using multifrequency temperature maps of the cosmic microwave background: Abell 399--Abell 401 cluster pair case study

Wed, 22/10/2025 - 10:54
arXiv:2510.18153v1 Announce Type: new Abstract: We present a multifrequency and multi-instrument methodology to study the physical properties of galaxy clusters and cosmic filaments using cosmic microwave background observations. Our approach enables simultaneous measurement of both the thermal (tSZ) and kinematic Sunyaev-Zeldovich (kSZ) effects, incorporates relativistic corrections, and models astrophysical foregrounds such as thermal dust emission. We do this by jointly fitting a single physical model across multiple maps from multiple instruments at different frequencies, rather than fitting a model to a single Compton-$y$ map. We demonstrate the success of this method by fitting the Abell 399-Abell 401 galaxy cluster pair and filament system using archival data from the Planck satellite and new, targeted deep data from the Atacama Cosmology Telescope, covering 11 different frequencies over 14 maps from 30 GHz to 545 GHz. Our tSZ results are consistent with previous work using Compton-$y$ maps. We measure the line-of-sight peculiar velocities of the cluster-filament system using the kSZ effect and find statistical uncertainties on individual cluster peculiar velocities of $\lesssim $600 km s$^{-1}$, which are competitive with current state-of-the-art measurements. Additionally, we measure the optical depth of the filament component with a signal-to-noise of 8.5$\sigma$ and reveal hints of its morphology. This modular approach is well-suited for application to future instruments across a wide range of millimeter and sub-millimeter wavebands.

MISTRAL: a model for AGN winds from radiatively efficient accretion in cosmological simulations

Wed, 22/10/2025 - 10:18
arXiv:2504.08041v2 Announce Type: replace Abstract: Feedback from active galactic nuclei (AGN) is crucial for regulating galaxy evolution. Motivated by observations of broad absorption line winds from rapidly accreting supermassive black holes (SMBHs), we introduce the Mistral AGN feedback model, implemented in the Arepo code. Mistral comes in two versions: continuous radial (Mistral-continuous) and stochastic bipolar momentum deposition (Mistral-stochastic). Using the framework of the IllustrisTNG simulations, we explore the effect of Mistral on BH and galaxy properties, through an idealized Milky Way-mass galaxy and cosmological zoom simulations run down to $z=2$. Unlike standard thermal AGN feedback prescriptions, Mistral generates galaxy-scale winds that mimic outflows driven by BH accretion. Mistral-continuous produces short-lived galactic fountains, and is inefficient at regulating the growth of massive galaxies at $z=2$. In contrast, Mistral-stochastic efficiently suppresses star formation in massive galaxies, reproduces the empirical stellar-to-halo mass relation, and yields a consistent trend of BH-stellar mass evolution. By supporting large-scale outflows while simultaneously preventing gas inflows, Mistral-stochastic additionally regulates the cold and hot gas fractions at both galaxy and halo scales. Mistral-stochastic therefore works self-consistently across the halo mass range explored $\left(10^{12}-3\times10^{13}\,\rm M_\odot\right)$, without adopting a SMBH-mass dependent AGN feedback scheme such as the one used in IllustrisTNG. Our model is a promising tool for predicting the impact of AGN winds on galaxy evolution, and interpreting the growing population of high-redshift galaxies and quasars observed by JWST. This work is part of the "Learning the Universe" collaboration, which aims to infer the physical processes governing the evolution of the Universe.

The impact of cosmic ray feedback during the epoch of reionisation

Tue, 21/10/2025 - 12:14
arXiv:2501.17239v2 Announce Type: replace Abstract: Galaxies form and evolve via a multitude of complex physics. In this work, we investigate the role of cosmic ray (CR) feedback in galaxy evolution and reionisation, by examining its impact on the escape of ionising radiation from galaxies. For this purpose, we present two Sphinx cosmological radiation-magneto-hydrodynamics simulations, allowing for the first time a study of the impact of CR feedback on thousands of resolved galaxies during the Epoch of Reionisation (EoR). The simulations differ in their feedback prescriptions: one adopts a calibrated strong supernova (SN) feedback, while the other simulation reduces the strength of SN feedback and includes CR feedback instead. We show that both comparably regulate star formation, reasonably match observations of high-redshift UV luminosity functions, and produce a similar amount of hydrogen ionising photons. In contrast to the model with strong SN feedback, the model with CRs lead to incomplete reionisation, which is in strong disagreement with observational estimates of the reionisation history. This is due to CR feedback shaping the ISM differently, filling with gas the low density cavities otherwise carved by SN explosions. As a result, this reduces the escape of ionising photons, at any halo mass, and primarily in the close vicinity of the stars. Our study indicates that CR feedback regulates galaxy growth during the EoR, but negatively affects reionisation, a tension which paves the way for further exploration and refinement of existing galaxy formation and feedback models. Such improvements are crucial in order to capture and understand the process of reionisation and the underlying evolution of galaxies through cosmic time.

QUIJOTE scientific results XIX. New constraints on the synchrotron spectral index using a semi-blind component separation method

Tue, 21/10/2025 - 12:02
arXiv:2510.17761v1 Announce Type: new Abstract: We introduce a novel approach to estimate the spectral index, $\beta_s$, of polarised synchrotron emission, combining the moment expansion of CMB foregrounds and the constrained-ILC method. We reconstruct the maps of the first two synchrotron moments, combining multi-frequency data, and apply the `T-T plot' technique between two moment maps to estimate the synchrotron spectral index. This approach offers a new technique for mapping the foreground spectral parameters, complementing the model-based parametric component separation methods. Applying this technique, we derive a new constraint on the spectral index of polarised synchrotron emission using QUIJOTE MFI wide-survey 11 and 13 GHz data, Wilkinson Microwave Anisotropy Probe (WMAP) data at K and Ka bands, and Planck LFI 30 GHz data. In the Galactic plane and North Polar Spur regions, we obtain an inverse-variance-weighted mean synchrotron index of $\beta_s = -3.11$ with a standard deviation of $0.21$ due to intrinsic scatter, consistent with previous results based on parametric methods using the same dataset. We find that the inverse-variance-weighted mean spectral index, including both statistical and systematic uncertainties, is $\beta_s^{\rm plane} = -3.05 \pm 0.01$ in the Galactic plane and $\beta_s^{\rm high\text{-}lat} = -3.13 \pm 0.02$ at high latitudes, indicating a moderate steepening of the spectral index from low to high Galactic latitudes. Our analysis indicates that, within the current upper limit on the AME polarisation fraction, our results are not subject to any appreciable bias. Furthermore, we infer the spectral index over the entire QUIJOTE survey region, partitioning the sky into 21 patches. This technique can be further extended to constrain the synchrotron spectral curvature by reconstructing higher-order moments when better-quality data become available.

The ALPINE-CRISTAL-JWST Survey: JWST/IFU Optical Observations for 18 Main-Sequence Galaxies at z=4-6

Tue, 21/10/2025 - 11:55
arXiv:2510.16111v1 Announce Type: new Abstract: To fully characterize the formation and evolution of galaxies, we need to observe their stars, gas, and dust on resolved spatial scales. We present the ALPINE-CRISTAL-JWST survey, which combines kpc-resolved imaging and spectroscopy from HST, JWST, and ALMA for 18 representative main-sequence galaxies at z=4-6 and log(M/$M_\odot$) > 9.5 to study their star formation, chemical properties, and extended gas reservoirs. The co-spatial measurements resolving the ionized gas, molecular gas, stars, and dust on 1-2 kpc scales make this a unique benchmark sample for the study of galaxy formation and evolution at z~5, connecting the Epoch of Reionization with the cosmic noon. In this paper, we outline the survey goals and sample selection, and present a summary of the available data for the 18 galaxies. In addition, we measure spatially integrated quantities (such as global gas metallicity), test different star formation rate indicators, and quantify the presence of H$\alpha$ halos. Our targeted galaxies are relatively metal rich (10-70% solar), complementary to JWST samples at lower stellar mass, and there is broad agreement between different star formation indicators. One galaxy has the signature of an active galactic nuclei (AGN) based on its emission line ratios. Six show broad H$\alpha$ emission suggesting type 1 AGN candidates. We conclude with an outlook on the exciting science that will be pursued with this unique sample in forthcoming papers.