Eric J. Baxter

Gravitational lensing by massive galaxy clusters distorts the observed cosmic microwave background (CMB) on arcminute scales, and these distortions carry information about cluster masses. Standard approaches to extracting cluster mass constraints from the CMB cluster lensing signal are either sub-optimal, ignore important physical or observational effects, are computationally intractable, or require additional work to turn the lensing measurements into constraints on cluster masses. We apply simulation based inference (SBI) using neural likelihood models to the problem. We show that in circumstances where the exact likelihood can be computed, the SBI constraints on cluster masses are in agreement with the exact likelihood, demonstrating that the SBI constraints are close to optimal. In scenarios where the exact likelihood cannot be feasibly computed, SBI still recovers unbiased estimates of individual cluster masses and combined constraints from multiple clusters. SBI will be a powerful tool for constraining the masses of galaxy clusters detected by future cosmic surveys. Code to run the analyses presented here will be made publicly available.

The two most favored explanations of the Fermi Galactic Center gamma-ray excess (GCE) are millisecond pulsars and self annihilation of the smooth dark matter halo of the galaxy. In order to distinguish between these possibilities, we would like to optimally use all information in the available data, including photon direction and energy information. To date, analyses of the GCE have generally treated directional and energy information separately, or have ignored one or the other completely. Here, we develop a method for analyzing the GCE that relies on simulation-based inference with neural posterior models to jointly analyze photon directional and spectral information while correctly accounting for the spatial and energy resolution of the telescope, here assumed to be the Fermi Large Area Telescope (LAT). Our results also have implications for analyses of the diffuse gamma-ray background, which we discuss.

Gravitational lensing by massive galaxy clusters distorts the observed cosmic microwave background (CMB) on arcminute scales, and these distortions carry information about cluster masses. Standard approaches to extracting cluster mass constraints from the CMB cluster lensing signal are either sub-optimal, ignore important physical or observational effects, are computationally intractable, or require additional work to turn the lensing measurements into constraints on cluster masses. We apply simulation based inference (SBI) using neural likelihood models to the problem. We show that in circumstances where the exact likelihood can be computed, the SBI constraints on cluster masses are in agreement with the exact likelihood, demonstrating that the SBI constraints are close to optimal. In scenarios where the exact likelihood cannot be feasibly computed, SBI still recovers unbiased estimates of individual cluster masses and combined constraints from multiple clusters. SBI will be a powerful tool for constraining the masses of galaxy clusters detected by future cosmic surveys. Code to run the analyses presented here will be made publicly available.

We present the Texas Euclid Survey for Lyα (TESLA), a spectroscopic survey in the 10 deg² of the Euclid North Ecliptic Pole (NEP) field. Using TESLA, we study how the physical properties of Lyα emitters (LAEs) correlate with Lyα emission to understand the escape of Lyα emission from galaxies at redshifts of 2–3.5. We present an analysis of 43 LAEs performed in the NEP field using early data from the TESLA survey. We use Subaru Hyper Suprime-Cam imaging in the grizy bands, Spitzer/IRAC channels 1 and 2 from the Hawaii 20 deg² (H20) survey, and spectra acquired by the Visible Integral-Field Replicable Unit Spectrograph (VIRUS) on the Hobby–Eberly Telescope. We perform spectral energy distribution (SED) fitting to compute the galaxy properties of 43 LAEs, and study correlations between stellar mass, star formation rate (SFR), and dust to the Lyα rest-frame equivalent width (W _Lyα ). We uncover marginal (1σ significance) correlations between stellar mass and W _Lyα , and SFR and W _Lyα , with a Spearman correlation coefficient of −0. 34 − .14 + .17 and −0. 37 − .14 + .16 , respectively. We show that the W _Lyα distribution of the 43 LAEs is consistent with being drawn from an exponential distribution with an e-folding scale of W ₀ = 150 Å. Once complete the TESLA survey will enable the study of ≳50,000 LAEs to explore more correlations between galaxy properties and W _Lyα . The large sample size will allow the construction of a predictive model for W _Lyα as a function of SED-derived galaxy properties, which could be used to improve Lyα-based constraints on reionization.

We consider gamma-ray signals of dark matter annihilation in extragalactic halos in the case where dark matter annihilates from ap-wave ord-wave state. In these scenarios, signals from extragalactic halos are enhanced relative to other targets, such as the Galactic Center or dwarf spheroidal galaxies, because the typical relative speed of the dark matter is larger in extragalactic halos. We perform a mock data analysis of gamma rays produced by dark matter annihilation in halos detected by the Sloan Digital Sky Survey. We include a model for uncorrelated galactic and extragalactic gamma ray backgrounds, as well as a simple model for backgrounds due to astrophysical processes in the extragalactic halos detected by the survey. We find that, for models which are still allowed by other gamma ray searches, searches of extragalactic halos with the current Fermi exposure can produce evidence for dark matter annihilation, though it is difficult to distinguish thep-wave andd-wave scenarios. With a factor 10× larger exposure, though, discrimination of the velocity-dependence is possible.

The cosmic web contains filamentary structure on a wide range of scales. On the largest scales, superclustering aligns multiple galaxy clusters along intercluster bridges, visible through their thermal Sunyaev–Zel’dovich signal in the cosmic microwave background. We demonstrate a new, flexible method to analyze the hot gas signal from multiscale extended structures. We use a Compton y-map from the Atacama Cosmology Telescope (ACT) stacked on redMaPPer cluster positions from the optical Dark Energy Survey (DES). Cutout images from the y-map are oriented with large-scale structure information from DES galaxy data such that the superclustering signal is aligned before being overlaid. We find evidence of an extended quadrupole moment of the stacked y signal at the 3.5σ level, demonstrating that the large-scale thermal energy surrounding galaxy clusters is anisotropically distributed. We compare our ACT × DES results with the Buzzard simulations, finding broad agreement. Using simulations, we highlight the promise of this novel technique for constraining the evolution of anisotropic, non-Gaussian structure using future combinations of microwave and optical surveys.

In the cold dark matter scenario, galactic dark matter halos are populated with a large number of smaller subhalos. Previous work has shown that dark matter annihilations in subhalos can generate a distinctive, non-Poisson signal in the gamma-ray photon counts probability distribution function (PDF). Here we show that the gamma-ray PDF also carries information about the velocity dependence of the dark matter annihilation cross section. After calculating the PDF assuming s-wave and Sommerfeld-enhanced annihilation, we perform a mock data analysis to illustrate how current and future observations can constrain the microphysics of the dark matter annihilation. We find that, with current Fermi data, and assuming a dark matter annihilation cross section roughly at the limit of current bounds from annihilation in dwarf spheroidal galaxies, one can potentially distinguish the non-Poissonian fluctuations expected from dark matter annihilation in subhalos from Poisson sources, as well as from dark matter models with an incorrect velocity-dependence. We explore how robust these results are to assumptions about the modeling of the galactic gamma-ray background, but further work is needed to determine the impact of realistic astrophysical source populations on our results. We also point out a four-parameter degeneracy between the velocity dependence of the dark matter annihilation, the minimum subhalo mass, the power law index of the subhalo mass function, and the normalization of the dark matter signal. This degeneracy can be broken with priors from N-body simulations or from observational constraints on the subhalo mass function.

We show the improvement to cosmological constraints from galaxy cluster surveys with the addition of cosmic microwave background (CMB)-cluster lensing data. We explore the cosmological implications of adding mass information from the 3.1σ detection of gravitational lensing of the CMB by galaxy clusters to the Sunyaev–Zel’dovich (SZ) selected galaxy cluster sample from the 2500 deg² SPT-SZ survey and targeted optical and X-ray follow-up data. In the ΛCDM model, the combination of the cluster sample with the Planck power spectrum measurements prefers σ 8 Ω m / 0.3 0.5 = 0.831 ± 0.020 . Adding the cluster data reduces the uncertainty on this quantity by a factor of 1.4, which is unchanged whether the 3.1σ CMB-cluster lensing measurement is included or not. We then forecast the impact of CMB-cluster lensing measurements with future cluster catalogs. Adding CMB-cluster lensing measurements to the SZ cluster catalog of the ongoing SPT-3G survey is expected to improve the expected constraint on the dark energy equation of state w by a factor of 1.3 to σ(w) = 0.19. We find the largest improvements from CMB-cluster lensing measurements to be for σ ₈, where adding CMB-cluster lensing data to the cluster number counts reduces the expected uncertainty on σ ₈ by respective factors of 2.4 and 3.6 for SPT-3G and CMB-S4.

We search for the signature of cosmological shocks in stacked gas pressure profiles of galaxy clusters using data from the South Pole Telescope (SPT). Specifically, we stack the latest Compton-y maps from the 2500 deg² SPT-SZ survey on the locations of clusters identified in that same data set. The sample contains 516 clusters with mean mass |$\langle M_{\rm 200m}\rangle = 10^{14.9} \, {\rm M}_\odot$| and redshift 〈z〉 = 0.55. We analyse in parallel a set of zoom-in hydrodynamical simulations from the three hundred project. The SPT-SZ data show two features: (i) a pressure deficit at R/R_200m = 1.08 ± 0.09, measured at 3.1σ significance and not observed in the simulations, and; (ii) a sharp decrease in pressure at R/R_200m = 4.58 ± 1.24 at 2.0σ significance. The pressure deficit is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions, and the second feature is consistent with accretion shocks seen in previous studies. We split the cluster sample by redshift and mass, and find both features exist in all cases. There are also no significant differences in features along and across the cluster major axis, whose orientation roughly points towards filamentary structure. As a consistency test, we also analyse clusters from the Planck and Atacama Cosmology Telescope Polarimeter surveys and find quantitatively similar features in the pressure profiles. Finally, we compare the accretion shock radius (⁠|$R_{\rm sh,\, acc}$|⁠) with existing measurements of the splashback radius (R_sp) for SPT-SZ and constrain the lower limit of the ratio, |$R_{\rm sh,\, acc}/R_{\rm sp}\gt 2.16 \pm 0.59$|⁠.