skip to content

Kavli Institute for Cosmology, Cambridge

 

COBRA: Optimal Factorization of Cosmological Observables

Wed, 16/04/2025 - 10:50
arXiv:2407.04660v2 Announce Type: replace Abstract: We introduce COBRA (Cosmology with Optimally factorized Bases of Radial Approximants), a novel framework for rapid computation of large-scale structure observables. COBRA separates scale dependence from cosmological parameters in the linear matter power spectrum while also minimising the number of necessary basis terms $N_b$, thus enabling direct and efficient computation of derived and nonlinear observables. Moreover, the dependence on cosmological parameters is efficiently approximated using radial basis function interpolation. We apply our framework to decompose the linear matter power spectrum in the standard $\Lambda$CDM scenario, as well as by adding curvature, dynamical dark energy and massive neutrinos, covering all redshifts relevant for Stage IV surveys. With only a dozen basis terms $N_b$, COBRA reproduces exact Boltzmann solver calculations to $\sim 0.1\%$ precision, which improves further to $0.02\%$ in the pure $\Lambda$CDM scenario. Using our decomposition, we recast the one-loop redshift space galaxy power spectrum in a separable minimal-basis form, enabling $\sim 4000$ model evaluations per second at $0.02\%$ precision on a single thread. This constitutes a considerable improvement over previously existing methods (e.g., FFTLog) opening a window for efficient computations of higher loop and higher order correlators involving multiple powers of the linear matter power spectra. The resulting factorisation can also be utilised in clustering, weak lensing and CMB analyses. Our implementation will be made public upon publication.

Flat-sky Angular Power Spectra Revisited

Wed, 16/04/2025 - 10:24
arXiv:2307.13768v2 Announce Type: replace Abstract: We revisit the flat-sky approximation for evaluating the angular power spectra of projected random fields by retaining information about the correlations along the line of sight. With broad, overlapping radial window functions, these line-of-sight correlations are suppressed and are ignored in the Limber approximation. However, retaining the correlations is important for narrow window functions or unequal-time spectra but introduces significant computational difficulties due to the highly oscillatory nature of the integrands involved. We deal with the integral over line-of-sight wave-modes in the flat-sky approximation analytically, using the FFTlog expansion of the 3D power spectrum. This results in an efficient computational method, which is a substantial improvement compared to any full-sky approaches. We apply our results to galaxy clustering (with and without redshift-space distortions), CMB lensing and galaxy lensing observables. For clustering, we find excellent agreement with the full-sky results on large (percent-level agreement) and intermediate or small (subpercent agreement) scales, dramatically out-performing the Limber approximation for both wide and narrow window functions, and in equal- and unequal-time cases. In the case of lensing, we show on the full sky that the angular power spectrum of the convergence can be very well approximated by projecting the 3D Laplacian (rather than the correct angular Laplacian) of the gravitational potential, even on large scales. Combining this approximation with our flat-sky techniques provides an efficient and accurate evaluation of the CMB lensing angular power spectrum on all scales.

Near-field imaging of local interference in radio interferometric data: Impact on the redshifted 21 cm power spectrum

Mon, 14/04/2025 - 11:45
arXiv:2503.21728v2 Announce Type: replace Abstract: Radio-frequency interference (RFI) is a major systematic limitation in radio astronomy, particularly for science cases requiring high sensitivity, such as 21 cm cosmology. Traditionally, RFI is dealt with by identifying its signature in the dynamic spectra of visibility data and flagging strongly affected regions. However, for RFI sources that do not occupy narrow regions in the time-frequency space, such as persistent local RFI, modeling these sources could be essential to mitigating their impact. This paper introduces two methods for detecting and characterizing local RFI sources from radio interferometric visibilities: matched filtering and maximum a posteriori (MAP) imaging. These algorithms use the spherical wave equation to construct three-dimensional near-field image cubes of RFI intensity from the visibilities. The matched filter algorithm can generate normalized maps by cross-correlating the expected contributions from RFI sources with the observed visibilities, while the MAP method performs a regularized inversion of the visibility equation in the near field. We developed a full polarization simulation framework for RFI and demonstrated the methods on simulated observations of local RFI sources. The stability, speed, and errors introduced by these algorithms were investigated, and, as a demonstration, the algorithms were applied to a subset of NenuFAR observations to perform spatial, spectral, and temporal characterization of two local RFI sources. We used simulations to assess the impact of local RFI on images, the uv plane, and cylindrical power spectra, and to quantify the level of bias introduced by the algorithms in order to understand their implications for the estimated 21 cm power spectrum with radio interferometers. The near-field imaging and simulation codes are publicly available in the Python library nfis.

JAX-bandflux: differentiable supernovae SALT modelling for cosmological analysis on GPUs

Mon, 14/04/2025 - 11:35
arXiv:2504.08081v1 Announce Type: new Abstract: JAX-bandflux is a JAX implementation of critical supernova modelling functionality for cosmological analysis. The codebase implements key components of the established library SNCosmo in a differentiable framework, offering efficient parallelisation and gradient-based optimisation capabilities through GPU acceleration. The package facilitates differentiable computation of supernova light curve measurements, supporting the inference of SALT parameters necessary for cosmological analysis.

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

Mon, 14/04/2025 - 11:31
arXiv:2504.08041v1 Announce Type: new 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, and reproduces the empirical stellar-to-halo mass and ($z=0$) BH-stellar mass relations. By supporting large-scale ($>50\,\rm kpc$) 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 radiatively efficient 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 Luminosity Function and Clustering of H$\alpha$ Emitting Galaxies at $z\approx4-6$ from a Complete NIRCam Grism Redshift Survey

Mon, 14/04/2025 - 11:28
arXiv:2504.08028v1 Announce Type: new Abstract: We study the luminosity function (LF) and clustering properties of 888 H$\alpha$ emitters (HAEs) at $3.75 < z < 6$ in the GOODS-N field. The sample, built from JWST CONGRESS and FRESCO NIRCam grism surveys using a novel redshift assignment algorithm, spans $\sim$62 arcmin$^2$ and reaches $L_{\rm H\alpha} \sim 10^{41.2} {\rm erg s^{-1}}$. We identify two prominent filamentary protoclusters at $z \approx 4.41$ and $z \approx 5.19$, hosting 98 and 144 HAEs, respectively. The observed H$\alpha$ LFs show similar shallow faint-end slopes for both protocluster and field galaxies at $z=3.75-5$, and for the protocluster at $z=5-6$ ($\alpha\approx 1.2$ to $-1.3$). In contrast, the field LF at $z=5-6$ is much steeper ($\alpha=-1.87_{-0.23}^{+0.30}$), suggesting that protocluster galaxies at $z > 5$ are more evolved, resembling those at $z=3.75-5$. The observed star formation rate density from H$\alpha$, integrated down to 0.45 ${\rm M_\odot yr^{-1}}$, is $0.050^{+0.002}_{-0.003}$ and $0.046^{+0.006}_{-0.004} M_\odot {\rm yr}^{-1} {\rm Mpc}^{-3}$ at $z=3.75-5$ and $z=5-6$, with protoclusters contributing $\sim$25% and 55%, respectively. This implies that a large fraction of star formation at $z > 4$ occurs in protoclusters. We conduct the first star-formation-rate-limited 3D clustering analysis at $z > 4$. We find the filamentary protocluster geometry flattens the power-law shape of the HAE auto-correlation functions, with slopes much shallower than typically assumed. The auto-correlation function of field HAEs have correlation lengths of $r_0 = 4.61^{+1.00}_{-0.68} h^{-1}{\rm Mpc}$ at $z \approx 4-5$ and $r_0 = 6.23^{+1.68}_{-1.13} h^{-1}{\rm Mpc}$ at $z=5-6$. Comparing the observed correlation functions with the UniverseMachine simulation, we infer the dark matter (sub-)halo masses of HAEs to be $\log (M_h/M_\odot)=11.0-11.2$ at $z\approx 4-6$, with a scatter of 0.4 dex.

Constraining Mixed Dark Matter models with high redshift Lyman-alpha forest data

Thu, 10/04/2025 - 10:03
arXiv:2504.06367v1 Announce Type: new Abstract: This study sets new constraints on Cold+Warm Dark Matter (CWDM) models by leveraging the small-scale suppression of structure formation imprinted in the Lyman-$\alpha$ forest. Using the Sherwood-Relics suite, we extract high-fidelity flux power spectra from simulated Lyman-$\alpha$ forest data, spanning a broad range of cosmologies and thermal histories. This enables precise constraints on the warm dark matter (WDM) fraction, $f_{\mathrm{WDM}}$, and the mass of the WDM particle, $m_{\mathrm{WDM}}$. A key advancement of our analysis is the integration of a neural network emulator directly at the likelihood level, significantly accelerating Bayesian parameter inference. With new observations of high-redshift ($z$ = 4.2$-$5.0) quasar spectra from UVES and HIRES, we establish stringent upper limits: for $m_{\mathrm{WDM}}$ = 1 keV, we find $f_{\mathrm{WDM}} < 0.16$ (2$\sigma$), with constraints loosening to 35\%, 50\%, and 67\% for $m_{\mathrm{WDM}}$ = 2, 3, and 4 keV, respectively. Our results for pure WDM reaffirm the lower bounds of previous work. Crucially, we account for the fixed resolution of simulations and the impact of patchy reionization, demonstrating their minimal influence on mixed dark matter constraints. This robustness paves the way for tighter bounds with improved statistical samples in the future. Our findings suggest that CWDM models can naturally accommodate mild suppression of matter clustering in the high redshift Lyman-$\alpha$ forest 1D flux power, potentially offering a resolution to some of the ongoing cosmological tensions at low redshifts, namely the $S_{8}$ tension.

The Ultraviolet Spectra of 2003fg-like Type Ia Supernovae

Wed, 09/04/2025 - 10:57
arXiv:2504.05869v1 Announce Type: new Abstract: 2003fg-like Type Ia supernovae (03fg-like SNe Ia) are a rare subtype of SNe Ia, photometrically characterized by broader optical light curves and bluer ultraviolet (UV) colors compared to normal SNe Ia. In this work, we study four 03fg-like SNe Ia using Swift UltraViolet and Optical Telescope (UVOT) grism observations to understand their unique UV properties and progenitor scenario(s). We report 03fg-like SNe Ia to have similar UV features and elemental compositions as normal SNe Ia, but with higher UV flux relative to optical. Previous studies have suggested that the UV flux levels of normal SNe Ia could be influenced by their progenitor properties, such as metallicity, with metal-poor progenitors producing higher UV flux levels. While 03fg-like SNe were previously reported to occur in low-mass and metal-poor host environments, our analysis indicates that their UV excess cannot be explained by their host-galaxy parameters. Instead, we demonstrate that the addition of a hot blackbody component, likely arising from the interaction with the circumstellar material (CSM), to the normal SN Ia spectrum, can reproduce their distinctive UV excess. This supports the hypothesis that 03fg-like SNe Ia could explode in a CSM-rich environment.

Decoding the variability in the star-formation histories of z ~ 0.8 galaxies

Tue, 08/04/2025 - 11:51
arXiv:2504.05281v1 Announce Type: new Abstract: The scatter of the star-forming main sequence (SFMS) holds a wealth of information about how galaxies evolve. The timescales encoded in this scatter can provide valuable insight into the relative importance of the physical processes regulating star formation. In this paper, we present a detailed observational analysis of the timescales imprinted in galaxy star-formation history (SFH) fluctuations by using the stochastic SFH model to fit 1928 massive, z ~ 0.8 galaxies in the LEGA-C survey. We find that the total intrinsic scatter of the SFMS is ~0.3 dex in galaxies with stellar masses $\gtrsim 10^{10}~\mathrm{M}_\odot$. This scatter decreases as the timescale over which SFRs are averaged increases, declining to a non-negligible ~0.15 - 0.25 dex at 2 Gyr, underscoring the importance of long-timescale SFH diversity to the SFMS scatter. Furthermore, galaxies currently above (below) the SFMS tend to have been above (below) the SFMS for at least ~1 Gyr, providing evidence that individual galaxies may follow different median tracks through SFR$-\mathrm{M}_*$ space. On shorter timescales (~30 - 100 Myr), galaxies' SFRs also vary on the order of ~0.1 - 0.2 dex. Our work supports the idea that the SFMS emerges from a population average of the pathways that individual galaxies trace through the SFR$-\mathrm{M}_*$ plane. The scatter reflects the long-term heterogeneity of these paths likely set by the evolutionary timescales of halo growth and cooling, accentuated by short-term variations reflecting the dynamical timescale of the galaxy and its interstellar medium. Our results emphasize the dynamic nature of the SFMS and the importance of understanding the diverse processes governing star formation.

MARTA: Temperature-temperature relationships and strong-line metallicity calibrations from multiple auroral lines detections at cosmic noon

Tue, 08/04/2025 - 11:42
arXiv:2504.03839v1 Announce Type: new Abstract: We present the first results from MARTA (Measuring Abundances at high Redshift with the T$_e$ Approach), a programme leveraging ultra-deep, medium-resolution JWST/NIRSpec spectroscopy to probe the interstellar medium (ISM) of star-forming galaxies at $z \sim 2 - 3$. We report detections of one or more auroral lines, including [O III]$\lambda4363$, [O II]$\lambda\lambda7320,7330$, [S II] $\lambda4068$, and [S III] $\lambda6312$ for 16 galaxies in the sample, providing measurements of multiple ionic temperatures. We tested the validity of the T[O II]-T[O III] relation at high redshift considering a total sample of 21 objects including literature data, and obtained a shallower slope than in the low-$z$ literature. However, such a slope is consistent with low-redshift data when ultra-low metallicity objects are considered. We assessed the correlation of the T[O II]-T[O III] relationship and its scatter on different physical parameters, finding a mild correlation with the ionisation parameter and radiation field hardness, while no significant correlation with gas density. The location of high-redshift data is also consistent with the low-$z$ literature in the T[O II]-T[S II], and T[S III]-T[O III] relations, although this conclusion is limited with low-number statistics. Finally, we leveraged our sample together with a comprehensive compilation of galaxies with [O III]$\lambda4363$ detections from the literature to recalibrate classical strong-line diagnostics at high redshift. MARTA represents a key addition in this space because it provides direct metallicities at moderately high oxygen abundances (12+log(O/H) $\sim8.0-8.4$).

ALMA Observations of [OI]145um and [NII]205um Emission lines from Star-Forming Galaxies at $z\sim7$

Tue, 08/04/2025 - 11:40
arXiv:2504.03831v1 Announce Type: new Abstract: We present results of new observations of [OI]145um and [NII]205um emission lines from four star-forming galaxies at redshifts between $z=6.58$ and $7.68$ that have previous detections of \Ciium\ and dust continua. Using ALMA, we successfully detect [OI]145um emission from all targets at $>4\,\sigma$ significance. However, [NII]205um emission is undetected in all galaxies (SNR $<3.5\,\sigma$) except for a tentative detection from A1689-zD1. From the observed high [CII]/[NII] emission line ratios ($\gtrsim20 - 80$), we find that most of the [CII]158um emission arise from neutral gas regions ($3\,\sigma$ lower limits of $\gtrsim 74 - 96\%$). From [OI]145um, [CII]158um lines, and infrared luminosities, we estimate the neutral gas densities of $n_{\rm H}=10^{3.5}$ - $10^6\,{\rm cm^{-3}}$ and the far-ultraviolet (FUV) radiation strengths of $G_0\sim10^{2.5}$-$10^{3}$. While the neutral gas densities are similar to those of high-redshift starburst galaxies, the FUV strengths are lower compared to both local and high-redshift starbursts. Finally, we estimate atomic hydrogen masses using [OI]145um emission lines and the oxygen abundances measured from recent JWST observations. We find gas mass ratios of $f_{\rm gas}\sim0.3$ - $0.8$, which are similar to earlier studies using [CII]158um. Starting from this pilot observation, future large [OI]145um emission line surveys will provide us with currently little-known neutral gas properties of star-forming galaxies in the early Universe.

The Dark Ages Explorer (DEX): a filled-aperture ultra-long wavelength radio interferometer on the lunar far side

Mon, 07/04/2025 - 11:52
arXiv:2504.03418v1 Announce Type: new Abstract: The measurement of the spatial fluctuations of the neutral hydrogen 21 cm signal arising during the Dark Ages and Cosmic Dawn periods of our Universe (z from 200 to 10) holds the potential to resolve these still-unexplored earliest phases of the evolution of matter structures. As these cosmological signals are very weak, large distributed telescopes are required at locations free from terrestrial radio interference and ionospheric disturbances. This paper presents a description of the scientific aims, the instrumental concept, and technological developments of an experiment - dubbed the Dark-ages EXplorer (DEX) - which would allow us to (a) measure the Global Signal and (b) measure the angular density fluctuations and conduct line-of-sight tomography in the Dark Ages and Cosmic Dawn epochs. Additional scientific goals are also briefly described. The experiment consists of a low-frequency radio interferometer, which should ideally be located on the far side of the Moon. The paper presents findings from an ESA Concurrent Design Facility (CDF) study, which was conducted to assess the feasibility of such a system using present-day technologies with a high TRL (Technology Readiness Level). Although the study finds that the number of antennas needed to achieve the primary scientific goals is not yet feasible at the moment, it points to a path of technological development that can lead to a realistic and valuable experiment in the medium-term future (i.e., the next decade(s)), as well as development of multi-purpose use technology that can be applied on Earth, and towards other lunar operations.

Superclustering with the Atacama Cosmology Telescope and Dark Energy Survey: II. Anisotropic large-scale coherence in hot gas, galaxies, and dark matter

Mon, 07/04/2025 - 11:49
arXiv:2409.04535v2 Announce Type: replace Abstract: Statistics that capture the directional dependence of the baryon distribution in the cosmic web enable unique tests of cosmology and astrophysical feedback. We use constrained oriented stacking of thermal Sunyaev-Zel'dovich (tSZ) maps to measure the anisotropic distribution of hot gas $2.5-40$ Mpc away from galaxy clusters embedded in massive filaments and superclusters. The cluster selection and orientation (at a scale of $\sim15$ Mpc) use Dark Energy Survey (DES) Year 3 data, while expanded tSZ maps from the Atacama Cosmology Telescope Data Release 6 enable a $\sim3\times$ more significant measurement of the extended gas compared to the technique's proof-of-concept. Decomposing stacks into cosine multipoles of order $m$, we detect a dipole ($m=1$) and quadrupole ($m=2$) at $8-10\sigma$, as well as evidence for $m=4$ signal at up to $6\sigma$, indicating sensitivity to late-time non-Gaussianity. We compare to the Cardinal simulations with spherical gas models pasted onto dark matter halos. The fiducial tSZ data can discriminate between two models that deplete pressure differently in low-mass halos (mimicking astrophysical feedback), preferring higher average pressure in extended structures. However, uncertainty in the amount of cosmic infrared background contamination reduces the constraining power. Additionally, we apply the technique to DES galaxy density and weak lensing to study for the first time their oriented relationships with tSZ. In the tSZ-to-lensing relation, averaged on 7.5 Mpc (transverse) scales, we observe dependence on redshift but not shape or radial distance. Thus, on large scales, the superclustering of gas pressure, galaxies, and total matter is coherent in shape and extent.

The DESI Y1 RR Lyrae catalog I: Empirical modeling of the cyclic variation of spectroscopic properties and a chemodynamical analysis of the outer halo

Mon, 07/04/2025 - 11:19
arXiv:2504.02924v1 Announce Type: new Abstract: We present the catalog of RR Lyrae stars observed in the first year of operations of the Dark Energy Spectroscopic Instrument (DESI) survey. This catalog contains 6,240 RR Lyrae stars out to $\sim100$\,kpc from the Galactic center and over 12,000 individual epochs with homogeneously-derived stellar atmospheric parameters. We introduce a novel methodology to model the cyclical variation of the spectroscopic properties of RR Lyrae from single-epoch measurements. We employ this method to infer the radial velocity and effective temperature variation of fundamental mode and first-overtone RR Lyrae stars and to determine their systemic velocities and mean temperatures. For fundamental mode pulsators, we obtain radial velocity curves with amplitudes of $\sim$30--50\,km\,s$^{-1}$ and effective temperature curves with 300--1,000\,K variations, whereas for first-overtone pulsators these amplitudes are $\sim20$\,km\,s$^{-1}$ and $\sim 600$\,K, respectively. We use our sample to study the metallicity distribution of the halo and its dependence on Galactocentric distance ($R_{\rm GC}$). Using a radius-dependent mixture model, we split the data into chemodynamically distinct components and find that our inner halo sample ($R_{\rm GC}\lesssim50$\,kpc) is predominantly composed of stars with [Fe/H] $\sim-1.5$\,dex and largely radial orbits (with an anisotropy parameter $\beta\sim0.94$), that we associate with the Gaia-Sausage-Enceladus merger event. Stars in the outer halo field exhibit a broader and more metal-poor [Fe/H] distribution with more circular orbits ($\beta\sim0.39$). The metallicity gradient of the metal-rich and the metal-poor components is found to be $0.005$ and $0.010$\,dex\,kpc$^{-1}$, respectively. Our catalog highlights DESI's tantalizing potential for studying the Milky Way and the pulsation properties of RR Lyrae stars in the era of large spectroscopic surveys.

JADES: comprehensive census of broad-line AGN from Reionization to Cosmic Noon revealed by JWST

Mon, 07/04/2025 - 11:09
arXiv:2504.03551v1 Announce Type: new Abstract: The depth and coverage of the first years of JWST observations have revealed low luminosity active galactic nuclei (AGN) across a wide redshift range, shedding light on black hole (BH) assembly and feedback. We present our spectroscopic sample of 34 Type 1 AGN obtained from JADES survey data and spanning $1.5 < z < 9$. Our sample of AGN probes a BH mass range of $10^{6-9}$~M$_{\odot}$ at bolometric luminosities down to $10^{43}$~erg~s$^{-1}$, implying generally sub-Eddington ratios of $<0.5L_{\rm Edd}$. Most of these AGN are hosted in low mass ($M_{\star}\sim10^8$~M$_{\odot}$) galaxies and are overmassive relative to the local $M_{BH}-M_{\star}$ relation, while remaining consistent with the local $M_{BH}$-$\sigma_*$ relation. The wide redshift range provided by our sample allows us to trace the emergence of local $M_{BH}$-$M_*$ scaling relation across the cosmic epoch. Additionally, we explore the capability of narrow-line diagnostics in identifying Type 2 AGN and find that a significant fraction of our AGN would be missed by them due to low metallicity or lack of high energy ionizing photons (potentially due to dust absorption, dense gas blanketing the broad and narrow line regions, or intrinsically soft ionizing spectra). We explore the UV luminosity function of AGN and their hosts and find that it is subject to significant cosmic variance and is also dependent on the AGN bolometric luminosity. Finally, we show that the electron and Balmer scattering scenarios recently proposed to explain the broad components of the Balmer lines are untenable on multiple grounds. There is no evidence that the black hole masses have been overestimated by orders of magnitude as proposed in those scenarios.

TDCOSMO XVII. New time delays in 22 lensed quasars from optical monitoring with the ESO-VST 2.6m and MPG 2.2m telescopes

Mon, 07/04/2025 - 11:03
arXiv:2504.02932v1 Announce Type: new Abstract: We present new time delays, the main ingredient of time delay cosmography, for 22 lensed quasars resulting from high-cadence r-band monitoring on the 2.6 m ESO VLT Survey Telescope and Max-Planck-Gesellschaft 2.2 m telescope. Each lensed quasar was typically monitored for one to four seasons, often shared between the two telescopes to mitigate the interruptions forced by the COVID-19 pandemic. The sample of targets consists of 19 quadruply and 3 doubly imaged quasars, which received a total of 1 918 hours of on-sky time split into 21 581 wide-field frames, each 320 seconds long. In a given field, the 5-{\sigma} depth of the combined exposures typically reaches the 27th magnitude, while that of single visits is 24.5 mag - similar to the expected depth of the upcoming Vera-Rubin LSST. The fluxes of the different lensed images of the targets were reliably de-blended, providing not only light curves with photometric precision down to the photon noise limit, but also high-resolution models of the targets whose features and astrometry were systematically confirmed in Hubble Space Telescope imaging. This was made possible thanks to a new photometric pipeline, lightcurver, and the forward modelling method STARRED. Finally, the time delays between pairs of curves and their uncertainties were estimated, taking into account the degeneracy due to microlensing, and for the first time the full covariance matrices of the delay pairs are provided. Of note, this survey, with 13 square degrees, has applications beyond that of time delays, such as the study of the structure function of the multiple high-redshift quasars present in the footprint at a new high in terms of both depth and frequency. The reduced images will be available through the European Southern Observatory Science Portal.

In-situ vs. ex-situ drivers of galaxy quenching: ubiquity of main sequence and critical black hole mass from the FLAMINGO simulation

Fri, 04/04/2025 - 10:26
arXiv:2504.02027v1 Announce Type: new Abstract: Exploiting a large sample of 5.3 million galaxies with $M_\ast\,{=}\,10^{10-11}\,{\rm M}_\odot$ from the highest-resolution FLAMINGO simulation, we carry out a statistical analysis of quiescent and star-forming galaxies to explore quenching mechanisms. From redshift $z\,{\simeq}\,7$ to 0, we find that the median star-formation rate of main-sequence galaxies is independent of the environment and of whether a galaxy is a central or satellite, whereas the fraction of quiescent galaxies is highly sensitive to both. By employing Random Forest (RF) classifiers, we demonstrate that black hole (BH) feedback is the most responsible quenching mechanism for both centrals and satellites, while halo mass is the second most significant. For satellites, a notable importance given by RF to stellar mass implies in-situ pre-quenching rather than ex-situ preprocessing prior to infall to the current host halo. In the cosmic afternoon of $z\,{=}\,$0--1, we identify two distinct regimes of evolution: at $M_{\rm BH}\,{\gtrsim}\,10^7\,{\rm M}_\odot$, essentially all galaxies are quenched regardless of their environment; at $M_{\rm BH}\,{\lesssim}\,10^7\,{\rm M}_\odot$, quenching is determined mainly by halo mass. Galaxies undergo a sharp transition from the main sequence to quiescence once their BH mass reaches $M_{\rm BH}\,{\simeq}\,10^7\,{\rm M}_\odot$ (typically when $M_\ast\,{\simeq}\,10^{10.5}\,{\rm M}_\odot$ and $M_{\rm h}\,{\simeq}\,10^{12}\,{\rm M}_\odot$) with a short quenching timescale of ${<}$1 Gyr. This transition is driven by a sudden change in the gas mass in the inner circum-galactic medium. Our results indicate that galaxy quenching arises from a combination of in-situ and ex-situ physical processes.

JWST MIRI reveals the diversity of nuclear mid-infrared spectra of nearby type-2 quasars

Thu, 03/04/2025 - 10:43
arXiv:2504.01595v1 Announce Type: new Abstract: Type-2 quasars (QSO2s) are active galactic nuclei (AGN) seen through a significant amount of dust and gas that obscures the central supermassive black hole and the broad line region. Here we present new mid-infrared spectra of the central kiloparsec of five optically-selected QSO2s at redshift z~0.1 obtained with JWST/MIRI/MRS. These QSO2s belong to the QSOFEED sample and they have log Lbol=45.5-46.0 erg/s, global SFRs that place them above the main sequence, and practically identical optical spectral shape and [OIII] luminosity, but their nuclear mid-infrared spectra exhibit an unexpected diversity of both continua and features. They show: 1) 9.7 micron silicate features going from emission (strength of S9.7=0.5) to relatively strong absorption (S9.7=-1.0) and 18 and 23 micron silicates either in emission or flat. In addition, two of the QSO2s show absorption bands of CO, H2O, and aliphatic grains, indicating different levels of nuclear obscuration across the sample. 2) [NeV]/[NeII] ratios ranging from 0.1 to 2.1 and [NeIII]/[NeII] from 1.0 to 3.5, indicating different coronal line and ionizing continuum strengths. 3) Warm molecular gas masses of 1-4x10^7 Msun and warm-to-cold gas mass ratios of 1-2%, with molecular gas excitation likely due to jet-induced shocks in J1430+1339, and to UV heating and/or turbulence in J1509+0434. 4) PAH emission features with equivalent widths ranging from <0.002 to 0.075 micron, from which we measure a larger contribution from neutral molecules (PAH 11.3/6.2=1.3-3.4) and SFRs<3-7 Msun/yr. This unprecedented dataset allowed us to start exploring the role of various AGN and galaxy properties including ionizing continuum, obscuration, electron density, and jet-ISM interactions on some of the spectral differences listed above, but larger samples are now required to fully understand the diversity of QSO2s' nuclear mid-infrared spectra.

Deciphering the Nature of Virgil: An Obscured AGN Lurking Within an Apparently Normal Lyman-{\alpha} Emitter During Cosmic Reionization

Thu, 03/04/2025 - 10:40
arXiv:2504.01852v1 Announce Type: new Abstract: We present a comprehensive analysis of the MIRI Extremely Red Object Virgil, a Lyman-$\alpha$ emitter at $z_{spec} = 6.6379 \pm 0.0035$ with the photometric properties of a Little Red Dot. Leveraging new JWST/MIRI imaging from the MIDIS and PAHSPECS programs, we confirm Virgil's extraordinary nature among galaxies in JADES/GOODS-South, exhibiting a strikingly red NIRCam-to-MIRI color (F444W $-$ F1500W = $2.84\pm0.04$~mag). Deep NIRSpec/PRISM spectroscopy from the OASIS program offers key insights into the host galaxy, revealing properties of an average star-forming galaxy during Cosmic Reionization, such as a subsolar metallicity, low-to-moderate dust content, and a relatively high ionization parameter and electron temperature. By estimating the star formation rate of Virgil from UV and H$\alpha$, we find evidence that the galaxy is either entering or fading out of a bursty episode. Although line-ratio diagnostics employed at high-$z$ would classify Virgil as an Active Galactic Nucleus (AGN), this classification becomes ambiguous once redshift evolution is considered. Nonetheless, Virgil occupies the same parameter space as recently confirmed AGNs at similar redshifts. The new deep MIRI data at 15~$\mu$m reinforce the AGN nature of Virgil, as inferred from multiple spectral energy distribution (SED) fitting codes. Virgil's rising infrared SED and UV excess resemble those of Dust-Obscured Galaxies (DOGs) studied with Spitzer at Cosmic Noon, particularly blue-excess HotDOGs. Our results highlight the need for a multi-wavelength approach incorporating MIRI to uncover such extreme sources at $z\gtrsim6$ and to shed light on the interplay between galaxy evolution and early black hole growth during Cosmic Reionization.

The Importance of Dust Distribution in Ionizing-photon Escape: NIRCam and MIRI Imaging of a Lyman Continuum-emitting Galaxy at z ~ 3.8

Thu, 03/04/2025 - 10:16
arXiv:2504.01067v1 Announce Type: new 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.