Adriaan J. Duivenvoorden

Simulation-based inference (SBI) is a powerful inference technique for cases where the exact functional form of the likelihood is not known. A prime example is the likelihood of cross-correlation power spectra of the cosmic microwave background (CMB) fields at low multipoles, ℓ ≲ 10. In this paper, we investigate a parity-violating cross-correlation between E - and B - mode polarization fields using SBI. The EB correlation at low ℓ is essential to distinguish between possible axion dark energy and dark matter interpretations of 'cosmic birefringence', a rotation of the plane of linear polarization of the CMB, recently reported from WMAP , Planck , and Atacama Cosmology Telescope data. We use neural likelihood estimation to infer the likelihood of the EB correlation at low ℓ and show that it is highly non-Gaussian. We then employ neural posterior estimation to constrain the scalar field mass ( m _ϕ ), the cosmic birefringence amplitude ( gϕ _in /2), and the instrumental miscalibration angle ( α ), from simulated datasets. We find that the posterior on m _ϕ shows two regimes, with a transition marked by 10 ^-32 eV, highlighting a strong sensitivity to the scale dependence of cosmic birefringence. To quantify this behavior, we compute the probability p ( m _ϕ < 10 ^-32 eV) for various fiducial values of m _ϕ . We find that α and the contribution of lensed B modes ultimately limit our ability to exclude the dark energy scenario fully.

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 A399–A401 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 ≲600 km s ⁻¹ , 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 σ and reveal hints of its morphology. This modular approach is well-suited for application to future instruments across a wide range of millimeter and submillimeter wavebands.

We present updated forecasts for the scientific performance of the degree-scale (0.5 deg FWHM at 93 GHz), deep-field survey to be conducted by the Simons Observatory (SO). By 2027, the SO Small Aperture Telescope (SAT) complement will be doubled from three to six telescopes, including a doubling of the detector count in the 93 GHz and 145 GHz channels to 48,160 detectors. Combined with a planned extension of the survey duration to 2035, this expansion will significantly enhance SO's search for a B -mode signal in the polarisation of the cosmic microwave background, a potential signature of gravitational waves produced in the very early Universe. Assuming a 1/ f noise model with knee multipole ℓ _knee = 50 and a moderately complex model for Galactic foregrounds, we forecast a 1 σ (or 68% confidence level) constraint on the tensor-to-scalar ratio r of σ _r = 1.2 × 10 ^-3 , assuming no primordial B-modes are present. This forecast assumes that 70% of the B -mode lensing signal can ultimately be removed using high resolution observations from the SO Large Aperture Telescope (LAT) and overlapping large-scale structure surveys. For more optimistic assumptions regarding foregrounds and noise, and assuming the same level of delensing, this forecast constraint improves to σ _r = 7 × 10 ^-4 . These forecasts represent a major improvement in SO's constraining power, being a factor of around 2.5 times better than what could be achieved with the originally planned campaign, which assumed the existing three SATs would conduct a five-year survey.

We present power spectra of the cosmic microwave background (CMB) anisotropy in temperature and polarization, measured from the Data Release 6 maps made from Atacama Cosmology Telescope (ACT) data. These cover 19,000 deg ² of sky in bands centered at 98, 150 and 220 GHz, with white noise levels three times lower than Planck in polarization. We find that the ACT angular power spectra estimated over 10,000 deg ² , and measured to arcminute scales in TT, TE and EE, are well fit by the sum of CMB and foregrounds, where the CMB spectra are described by the ΛCDM model. Combining ACT with larger-scale Planck data, the joint P-ACT dataset provides tight limits on the ingredients, expansion rate, and initial conditions of the universe. We find similar constraining power, and consistent results, from either the Planck power spectra or from ACT combined with WMAP data, as well as from either temperature or polarization in the joint P-ACT dataset. When combined with CMB lensing from ACT and Planck , and baryon acoustic oscillation data from the Dark Energy Spectroscopic Instrument (DESI DR1), we measure a baryon density of Ω _b h ² = 0.0226 ± 0.0001, a cold dark matter density of Ω _c h ² = 0.118 ± 0.001, a Hubble constant of H ₀ = 68.22 ± 0.36 km/s/Mpc, a spectral index of n _s = 0.974 ± 0.003, and an amplitude of density fluctuations of σ ₈ = 0.813 ± 0.005. Including the DESI DR2 data tightens the Hubble constant to H ₀ = 68.43 ± 0.27 km/s/Mpc; ΛCDM parameters agree between the P-ACT and DESI DR2 data at the 1.6 σ level. We find no evidence for excess lensing in the power spectrum, and no departure from spatial flatness. The contribution from Sunyaev-Zel'dovich (SZ) anisotropy is detected at high significance; we find evidence for a tilt with suppressed small-scale power compared to our baseline SZ template spectrum, consistent with hydrodynamical simulations with feedback.

We present Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) maps of the Cosmic Microwave Background temperature and polarization anisotropy at arcminute resolution over three frequency bands centered on 98, 150 and 220 GHz. The maps are based on data collected with the AdvancedACT camera over the period 2017–2022 and cover 19,000 square degrees with a median combined depth of 10 μK arcmin. We describe the instrument, mapmaking and map properties and illustrate them with a number of figures and tables. The ACT DR6 maps and derived products are available on LAMBDA at https://lambda.gsfc.nasa.gov/product/act/actadv_prod_table.html . We also provide an interactive web atlas at https://phy-act1.princeton.edu/public/snaess/actpol/dr6/atlas and HiPS data sets in Aladin (e.g. https://alasky.cds.unistra.fr/ACT/DR4DR6/color_CMB ).

We use new cosmic microwave background (CMB) primary temperature and polarization anisotropy measurements from the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6) to test foundational assumptions of the standard cosmological model, ΛCDM, and set constraints on extensions to it. We derive constraints from the ACT DR6 power spectra alone, as well as in combination with legacy data from the Planck mission. To break geometric degeneracies, we include ACT and Planck CMB lensing data and baryon acoustic oscillation data from DESI Year-1. To test the dependence of our results on non-ACT data, we also explore combinations replacing Planck with WMAP and DESI with BOSS, and further add supernovae measurements from Pantheon+ for models that affect the late-time expansion history. We verify the near-scale-invariance (running of the spectral index dn _s / d ln k = 0.0062 ± 0.0052) and adiabaticity of the primordial perturbations. Neutrino properties are consistent with Standard Model predictions: we find no evidence for new light, relativistic species that are free-streaming ( N _eff = 2.86 ± 0.13, which combined with astrophysical measurements of primordial helium and deuterium abundances becomes N _eff = 2.89 ± 0.11), for non-zero neutrino masses (∑ m _ν < 0.089 eV at 95% CL), or for neutrino self-interactions. We also find no evidence for self-interacting dark radiation ( N _idr < 0.134), or for early-universe variation of fundamental constants, including the fine-structure constant ( α _EM / α _EM,0 = 1.0043 ± 0.0017) and the electron mass ( m _e / m _e,0 = 1.0063 ± 0.0056). Our data are consistent with standard big bang nucleosynthesis (we find Y _p = 0.2312 ± 0.0092), the COBE/FIRAS -inferred CMB temperature (we find T _CMB = 2.698 ± 0.016 K), a dark matter component that is collisionless and with only a small fraction allowed as axion-like particles, a cosmological constant ( w = -0.986 ± 0.025), and the late-time growth rate predicted by general relativity ( γ = 0.663 ± 0.052). We find no statistically significant preference for a departure from the baseline ΛCDM model. In fits to models invoking early dark energy, primordial magnetic fields, or an arbitrary modified recombination history, we find H ₀ = 69.9 ^+0.8 _-1.5 , 69.1 ± 0.5, or 69.6 ± 1.0 km/s/Mpc, respectively; using BOSS instead of DESI BAO data reduces the central values of these constraints by 1–1.5 km/s/Mpc while only slightly increasing the error bars. In general, models introduced to increase the Hubble constant or to decrease the amplitude of density fluctuations inferred from the primary CMB are not favored over ΛCDM by our data.

We discuss the model of astrophysical emission at millimeter wavelengths used to characterize foregrounds in the multi-frequency power spectra of the Atacama Cosmology Telescope (ACT) Data Release 6 (DR6), expanding on Louis et al. (2025) ( 2503.14452 ). We detail several tests to validate the capability of the DR6 parametric foreground model to describe current observations and complex simulations, and show that cosmological parameter constraints are robust against model extensions and variations. We demonstrate consistency of the model with pre-DR6 ACT data and observations from Planck and the South Pole Telescope. We evaluate the implications of using different foreground templates and extending the model with new components and/or free parameters. In all scenarios, the DR6 ΛCDM and ΛCDM+ N _eff cosmological parameters shift by less than 0.5 σ relative to the baseline constraints. Some foreground parameters shift more; we estimate their systematic uncertainties associated with modeling choices. From our constraint on the kinematic Sunyaev-Zel'dovich power, we obtain a conservative limit on the duration of reionization of Δ z _rei < 4.4, assuming a reionization midpoint consistent with optical depth measurements and a minimal low-redshift contribution, with varying assumptions for this component leading to tighter limits. Finally, we analyze realistic non-Gaussian, correlated microwave sky simulations containing Galactic and extragalactic foreground fields, built independently of the DR6 parametric foreground model. Processing these simulations through the DR6 power spectrum and likelihood pipeline, we recover the input cosmological parameters of the underlying cosmic microwave background field, a new demonstration for small-scale CMB analysis. These tests validate the robustness of the ACT DR6 foreground model and cosmological parameter constraints.

The Simons Observatory (SO) Small Aperture Telescopes (SATs) will observe the Cosmic Microwave Background (CMB) temperature and polarization at six frequency bands. Within these bands, the angular response of the telescope (beam) is convolved with the instrument's spectral response (commonly called bandpass) and the signal from the sky, which leads to the band-averaged telescope beam response, which is sampled and digitized. The spectral properties of the band-averaged beam depend on the natural variation of the beam within the band, referred to as beam chromaticity. In this paper, we quantify the impact of the interplay of beam chromaticity and intrinsic frequency scaling from the various components that dominate the polarized sky emission on the tensor-to-scalar ratio, r, and foreground parameters. We do so by employing a parametric power-spectrum-based foreground component separation algorithm, namely BBPower, to which we provide beam-convolved time domain simulations performed with the beamconv software while assuming an idealized version of the SO SAT optics. We find a small, 0.02σ, bias on r, due to beam chromaticity, which seems to mostly impact the dust spatial parameters, causing a maximum 0.77σ bias on the dust B-mode spectra amplitude, A_d , when employing Gaussian foreground simulations. However, we find all parameter biases to be smaller than 1σ at all times, independently of the foreground model. This includes the case where we introduce additional uncertainty on the bandpass shape, which accounts for approximately half of the total allowed gain uncertainty, as estimated in previous work for the SO SATs.