Elia Stefano Battistelli

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.

Galaxy clusters form through hierarchical assembly, where smaller substructures merge to build the largest gravitationally bound objects in the Universe. These mergers, combined with feedback from active galactic nucleus (AGN), filamentary accretion, and other energy injection processes, generate turbulence and perturbations within the intracluster medium (ICM). X-ray and Sunyaev–Zel’dovich (SZ) observations can be utilized to measure these ICM density and pressure inhomogeneities, in turn providing constraints on the effective equation of state of the perturbations and ICM velocities. In this work, we analyze deep SRG-eROSITA (hereafter eROSITA) and Planck/South Pole Telescope (SPT) observations of A3266, a dynamically complex merging cluster with elongated morphology and significant substructure. We measure pressure and density fluctuations, and compute the power spectra and deprojected 3D amplitudes of these perturbations. We estimate the ratio of pressure-to-density fluctuation amplitudes as 1.07 ± 0.55 and nonthermal pressure support 0.067 ± 0.049. Density fluctuations are found to be stronger in the northern sector of the cluster compared to the south, consistent with ongoing accretion along a filamentary structure revealed by eROSITA. Further, we find the amplitude of density fluctuations increases with radius, qualitatively consistent with the trend found in cosmological simulations. Uncertainties in our results are dominated by the relatively low sensitivity of current Planck/SPT data, suggesting that improvements in SZ data quality could substantially improve our understanding of ICM energy injection, transport, and dissipation from this technique.

We present a multi-frequency, multi-epoch catalog of extragalactic sources. The catalog is based on 150, 220, and 280 GHz observations carried out in 2008, 2009, and 2010 using the Millimeter Bolometric Array Camera on the Atacama Cosmology Telescope. We also present and release 280 GHz maps from 2008 and 2010. The catalog contains 483 sources found in a sky area of ∼600 square degrees. It was obtained by cross-matching sources found in 11 sub-catalogs, one for each season and frequency band. We also include coadded data from ∼150 and ∼160 square degrees using two and three years of overlapping observations. We divide the sources into two populations, synchrotron and dusty emitters, based on their spectral behavior in the 150-280 GHz frequency range. We find 284 synchrotron sources and 183 dusty source candidates. Our cross-matching with catalogs from radio to X-ray results in 251 synchrotron sources (88%) and 92 dusty sources (51%) with counterparts and suggests that 91 dusty candidates are not in existing catalogs. We study the variability and number counts of each population. In the case of synchrotron sources, we find year-to-year variability, with a mean value around 35%. As expected, we find no evidence of dusty source variability. Our number counts generally agree with previous measurements and models, except for dusty sources at 280 GHz, where some models overestimate our results.

The search for relics from the inflation era in the form of B-mode polarization of the CMB is a major challenge in cosmology. The main obstacle appears to come from the complexity of Galactic foregrounds that need to be removed. Multi-frequency observations are key to mitigating their contamination and mapping primordial fluctuations. We present “Spectral-Imaging”, a method to reconstruct sub-frequency maps of the CMB polarization within the instrument's physical bandwidth, a unique feature of Bolometric Interferometry that could be crucial for foreground mitigation as it provides an increased spectral resolution. Our technique uses the frequency evolution of the shape of the Bolometric Interferometer's synthesized beam to reconstruct frequency information from the time domain data. We reconstruct sub-frequency maps using an inverse problem approach based on detailed modeling of the instrument acquisition. We use external data to regularize the convergence of the estimator and account for bandpass mismatch and varying angular resolution. The reconstructed maps are unbiased and allow exploiting the spectral-imaging capacity of QUBIC. Using end-to-end simulations of the QUBIC instrument, we perform a cross-spectra analysis to extract a forecast on the tensor-to-scalar ratio constraint of σ ( r )= 0.0225 after component separation.

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.

The detection of B-modes in the CMB polarization pattern is a major issue in modern cosmology and must therefore be handled with analytical methods that produce reliable results. We describe a method that uses the frequency dependence of the QUBIC synthesized beam to perform component separation at the map-making stage, to obtain more precise results. We aim to demonstrate the feasibility of component separation during the map-making stage in time domain space. This new technique leads to a more accurate description of the data and reduces the biases in cosmological analyses. The method uses a library for highly parallel computation which facilitates the programming and allows the description of experiments as easily manipulated operators. These operators can be combined to obtain a joint analysis using several experiments to maximize precision. The results show that the method works well and permits end-to-end analysis for the CMB experiments, and in particular, for QUBIC. The method includes astrophysical foregrounds, and also systematic effects like gain variation in the detectors. We developed a software pipeline that produces uncertainties on the tensor-to-scalar ratio at the level of σ(r) ∼ 0.023 using only QUBIC simulated data.

The powerful jets of blazars have been historically considered as likely sites of high-energy cosmic-ray acceleration. However, the particulars of the launched jet and the locations of leptonic and hadronic jet loading remain unclear. In the case when leptonic and hadronic particle injection occur jointly, a temporal correlation between synchrotron radiation and neutrino production is expected. We use a first catalog of millimeter wavelength (95–225 GHz) blazar light curves from the Atacama Cosmology Telescope for a time-dependent correlation with 12 yr of muon neutrino events from the IceCube South Pole Neutrino Observatory. Such millimeter emission traces activity of the bright jet base, which is often self-absorbed at lower frequencies and potentially gamma-ray opaque. We perform an analysis of the population, as well as analyses of individual, selected sources. We do not observe a significant signal from the stacked population. TXS 0506+056 is found as the most significant, individual source, though this detection is not globally significant in our analysis of selected active galactic nuclei. Our results suggest that the majority of millimeter-bright blazars are neutrino dim. In general, it is possible that many blazars have lighter, leptonic jets, or that only selected blazars provide exceptional conditions for neutrino production.

The high electron temperature in galaxy clusters ( > 1 keV or > 10 7 K ) leads to corrections at the level of a few percent in their thermal Sunyaev-Zeldovich effect signatures. Both the size and frequency dependence of these corrections, which are known as relativistic temperature corrections, depend upon the temperature of the objects. In this work we exploit this effect to measure the average temperature of a stack of Compton- y selected clusters. Specifically, we apply the “spectroscopic method” and search for the temperature that best fits the clusters’ signal measured at frequencies from 30 to 545 GHz by the Atacama Cosmology Telescope and satellite. We measure the average temperature of clusters detected in the Atacama Cosmology Telescope maps to be 8.5 ± 2.4 keV , with an additional systematic error of comparable amplitude dominated by passband uncertainty. Upcoming surveys, such as the Simons Observatory and CMB-S4, have the potential to dramatically improve upon these measurements and thereby enable precision studies of cluster temperatures with millimeter observations. The key challenge for future observations will be mitigating instrumental systematic effects, which already limit this analysis.

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.