Researcher profile

David N. Spergel

David N. Spergel contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

22 published item(s)

preprint2026arXiv

Cosmological back-reaction of baryons on dark matter in the CAMELS simulations

Baryonic processes such as radiative cooling and feedback from massive stars and active galactic nuclei (AGN) directly redistribute baryons in the Universe but also indirectly redistribute dark matter due to changes in the gravitational potential. In this work, we investigate this &#34;back-reaction&#34; of baryons on dark matter using thousands of cosmological hydrodynamic simulations from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project, including parameter variations in the SIMBA, IllustrisTNG, ASTRID, and Swift-EAGLE galaxy formation models. Matching haloes to corresponding N-body (dark matter-only) simulations, we find that virial masses decrease owing to the ejection of baryons by feedback. Relative to N-body simulations, halo profiles show an increased dark matter density in the center (due to radiative cooling) and a decrease in density farther out (due to feedback), with both effects being strongest in SIMBA (> 450% increase at r < 0.01 Rvir). The clustering of dark matter strongly responds to changes in baryonic physics, with dark matter power spectra in some simulations from each model showing as much as 20% suppression or increase in power at k ~ 10 h/Mpc relative to N-body simulations. We find that the dark matter back-reaction depends intrinsically on cosmology (Omega_m and sigma_8) at fixed baryonic physics, and varies strongly with the details of the feedback implementation. These results emphasize the need for marginalizing over uncertainties in baryonic physics to extract cosmological information from weak lensing surveys as well as their potential to constrain feedback models in galaxy evolution.

preprint2026arXiv

Unified and consistent structure growth measurements from joint ACT, SPT and \textit{Planck} CMB lensing

We present the tightest cosmic microwave background (CMB) lensing constraints to date on the growth of structure by combining CMB lensing measurements from the Atacama Cosmology Telescope (ACT), the South Pole Telescope (SPT) and \textit{Planck}. Each of these surveys individually provides lensing measurements with similarly high statistical power, achieving signal-to-noise ratios of approximately 40. The combined lensing bandpowers represent the most precise CMB lensing power spectrum measurement to date with a signal-to-noise ratio of 61 and an amplitude of $A_\mathrm{lens}^\mathrm{recon} = 1.025 \pm 0.017$ with respect to the theory prediction from the best-fit CMB \textit{Planck}-ACT cosmology. The bandpowers from all three lensing datasets, analyzed jointly, yield a $1.6\%$ measurement of the parameter combination $S_8^\mathrm{CMBL} \equiv σ_8\,(Ω_m/0.3)^{0.25} = 0.825^{+0.015}_{-0.013}$. Including Dark Energy Spectroscopic Instrument (DESI) Baryon Acoustic Oscillation (BAO) data improves the constraint on the amplitude of matter fluctuations to $σ_8 = 0.829 \pm 0.009$ (a $1.1\%$ determination). When combining with uncalibrated supernovae from \texttt{Pantheon+}, we present a $4\%$ sound-horizon-independent estimate of $H_0=66.4\pm2.5\,\mathrm{km\,s^{-1}\,Mpc^{-1}} $. The joint lensing constraints on structure growth and present-day Hubble rate are fully consistent with a $Λ$CDM model fit to the primary CMB data from \textit{Planck} and ACT. While the precise upper limit is sensitive to the choice of data and underlying model assumptions, when varying the neutrino mass sum within the $Λ\mathrm{CDM}$ cosmological model, the combination of primary CMB, BAO and CMB lensing drives the probable upper limit for the mass sum towards lower values, comparable to the minimum mass prior required by neutrino oscillation experiments.

preprint2022arXiv

Breaking baryon-cosmology degeneracy with the electron density power spectrum

Uncertain feedback processes in galaxies affect the distribution of matter, currently limiting the power of weak lensing surveys. If we can identify cosmological statistics that are robust against these uncertainties, or constrain these effects by other means, then we can enhance the power of current and upcoming observations from weak lensing surveys such as DES, Euclid, the Rubin Observatory, and the Roman Space Telescope. In this work, we investigate the potential of the electron density auto-power spectrum as a robust probe of cosmology and baryonic feedback. We use a suite of (magneto-)hydrodynamic simulations from the CAMELS project and perform an idealized analysis to forecast statistical uncertainties on a limited set of cosmological and physically-motivated astrophysical parameters. We find that the electron number density auto-correlation, measurable through either kinematic Sunyaev-Zel&#39;dovich observations or through Fast Radio Burst dispersion measures, provides tight constraints on $Ω_{m}$ and the mean baryon fraction in intermediate-mass halos, $\bar{f}_{\mathrm{bar}}$. By obtaining an empirical measure for the associated systematic uncertainties, we find these constraints to be largely robust to differences in baryonic feedback models implemented in hydrodynamic simulations. We further discuss the main caveats associated with our analysis, and point out possible directions for future work.

preprint2022arXiv

Cosmology with one galaxy?

Galaxies can be characterized by many internal properties such as stellar mass, gas metallicity, and star-formation rate. We quantify the amount of cosmological and astrophysical information that the internal properties of individual galaxies and their host dark matter halos contain. We train neural networks using hundreds of thousands of galaxies from 2,000 state-of-the-art hydrodynamic simulations with different cosmologies and astrophysical models of the CAMELS project to perform likelihood-free inference on the value of the cosmological and astrophysical parameters. We find that knowing the internal properties of a single galaxy allow our models to infer the value of $Ω_{\rm m}$, at fixed $Ω_{\rm b}$, with a $\sim10\%$ precision, while no constraint can be placed on $σ_8$. Our results hold for any type of galaxy, central or satellite, massive or dwarf, at all considered redshifts, $z\leq3$, and they incorporate uncertainties in astrophysics as modeled in CAMELS. However, our models are not robust to changes in subgrid physics due to the large intrinsic differences the two considered models imprint on galaxy properties. We find that the stellar mass, stellar metallicity, and maximum circular velocity are among the most important galaxy properties to determine the value of $Ω_{\rm m}$. We believe that our results can be explained taking into account that changes in the value of $Ω_{\rm m}$, or potentially $Ω_{\rm b}/Ω_{\rm m}$, affect the dark matter content of galaxies. That effect leaves a distinct signature in galaxy properties to the one induced by galactic processes. Our results suggest that the low-dimensional manifold hosting galaxy properties provides a tight direct link between cosmology and astrophysics.

preprint2022arXiv

Predicting the Thermal Sunyaev-Zel&#39;dovich Field using Modular and Equivariant Set-Based Neural Networks

Theoretical uncertainty limits our ability to extract cosmological information from baryonic fields such as the thermal Sunyaev-Zel&#39;dovich (tSZ) effect. Being sourced by the electron pressure field, the tSZ effect depends on baryonic physics that is usually modeled by expensive hydrodynamic simulations. We train neural networks on the IllustrisTNG-300 cosmological simulation to predict the continuous electron pressure field in galaxy clusters from gravity-only simulations. Modeling clusters is challenging for neural networks as most of the gas pressure is concentrated in a handful of voxels and even the largest hydrodynamical simulations contain only a few hundred clusters that can be used for training. Instead of conventional convolutional neural net (CNN) architectures, we choose to employ a rotationally equivariant DeepSets architecture to operate directly on the set of dark matter particles. We argue that set-based architectures provide distinct advantages over CNNs. For example, we can enforce exact rotational and permutation equivariance, incorporate existing knowledge on the tSZ field, and work with sparse fields as are standard in cosmology. We compose our architecture with separate, physically meaningful modules, making it amenable to interpretation. For example, we can separately study the influence of local and cluster-scale environment, determine that cluster triaxiality has negligible impact, and train a module that corrects for mis-centering. Our model improves by 70 % on analytic profiles fit to the same simulation data. We argue that the electron pressure field, viewed as a function of a gravity-only simulation, has inherent stochasticity, and model this property through a conditional-VAE extension to the network. This modification yields further improvement by 7 %, it is limited by our small training set however. (abridged)

preprint2022arXiv

Reionization with Simba: How much does astrophysics matter in modeling cosmic reionization?

Traditional large-scale models of reionization usually employ simple deterministic relations between halo mass and luminosity to predict how reionization proceeds. We here examine the impact on modelling reionization of using more detailed models for the ionizing sources as identified within the $100~{\rm Mpc/h}$ cosmological hydrodynamic simulation {\sc Simba}, coupled with post-processed radiative transfer. Comparing with simple (one-to-one) models, the main difference of using {\sc Simba} sources is the scatter in the relation between dark matter halos and star formation, and hence ionizing emissivity. We find that, at the power spectrum level, the ionization morphology remains mostly unchanged, regardless of the variability in the number of sources or escape fraction. In particular, the power spectrum shape remains unaffected and its amplitude changes slightly by less than 5-10\%, throughout reionization, depending on the scale and neutral fraction. Our results show that simplified models of ionizing sources remain viable to efficiently model the structure of reionization on cosmological scales, although the precise progress of reionization requires accounting for the scatter induced by astrophysical effects.

preprint2022arXiv

The Atacama Cosmology Telescope: Constraints on Pre-Recombination Early Dark Energy

The early dark energy (EDE) scenario aims to increase the value of the Hubble constant ($H_0$) inferred from cosmic microwave background (CMB) data over that found in $Λ$CDM, via the introduction of a new form of energy density in the early universe. The EDE component briefly accelerates cosmic expansion just prior to recombination, which reduces the physical size of the sound horizon imprinted in the CMB. Previous work has found that non-zero EDE is not preferred by Planck CMB power spectrum data alone, which yield a 95% confidence level (CL) upper limit $f_{\rm EDE} < 0.087$ on the maximal fractional contribution of the EDE field to the cosmic energy budget. In this paper, we fit the EDE model to CMB data from the Atacama Cosmology Telescope (ACT) Data Release 4. We find that a combination of ACT, large-scale Planck TT (similar to WMAP), Planck CMB lensing, and BAO data prefers the existence of EDE at $>99.7$% CL: $f_{\rm EDE} = 0.091^{+0.020}_{-0.036}$, with $H_0 = 70.9^{+1.0}_{-2.0}$ km/s/Mpc (both 68% CL). From a model-selection standpoint, we find that EDE is favored over $Λ$CDM by these data at roughly $3σ$ significance. In contrast, a joint analysis of the full Planck and ACT data yields no evidence for EDE, as previously found for Planck alone. We show that the preference for EDE in ACT alone is driven by its TE and EE power spectrum data. The tight constraint on EDE from Planck alone is driven by its high-$\ell$ TT power spectrum data. Understanding whether these differing constraints are physical in nature, due to systematics, or simply a rare statistical fluctuation is of high priority. The best-fit EDE models to ACT and Planck exhibit coherent differences across a wide range of multipoles in TE and EE, indicating that a powerful test of this scenario is anticipated with near-future data from ACT and other ground-based experiments.

preprint2022arXiv

The Atacama Cosmology Telescope: The Persistence of Neutrino Self-Interaction in Cosmological Measurements

We use data from the Atacama Cosmology Telescope (ACT) DR4 to search for the presence of neutrino self-interaction in the cosmic microwave background. Consistent with prior works, the posterior distributions we find are bimodal, with one mode consistent with $Λ$CDM and one where neutrinos strongly self-interact. By combining ACT data with large-scale information from WMAP, we find that a delayed onset of neutrino free streaming caused by significantly strong neutrino self-interaction is compatible with these data at the $2-3σ$ level. As seen in the past, the preference shifts to $Λ$CDM with the inclusion of Planck data. We determine that the preference for strong neutrino self-interaction is largely driven by angular scales corresponding to $700 \lesssim \ell \lesssim 1000$ in the ACT E-mode polarization data. This region is expected to be key to discriminate between neutrino self-interacting modes and will soon be probed with more sensitive data.

preprint2021arXiv

Combining information from multiple cosmological surveys: inference and modeling challenges

The tightest and most robust cosmological results of the next decade will be achieved by bringing together multiple surveys of the Universe. This endeavor has to happen across multiple layers of the data processing and analysis, e.g., enhancements are expected from combining Euclid, Rubin, and Roman (as well as other surveys) not only at the level of joint processing and catalog combination, but also during the post-catalog parts of the analysis such as the cosmological inference process. While every experiment builds their own analysis and inference framework and creates their own set of simulations, cross-survey work that homogenizes these efforts, exchanges information from numerical simulations, and coordinates details in the modeling of astrophysical and observational systematics of the corresponding datasets is crucial.

preprint2021arXiv

Finding universal relations in subhalo properties with artificial intelligence

We use a generic formalism designed to search for relations in high-dimensional spaces to determine if the total mass of a subhalo can be predicted from other internal properties such as velocity dispersion, radius, or star-formation rate. We train neural networks using data from the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) project and show that the model can predict the total mass of a subhalo with high accuracy: more than 99% of the subhalos have a predicted mass within 0.2 dex of their true value. The networks exhibit surprising extrapolation properties, being able to accurately predict the total mass of any type of subhalo containing any kind of galaxy at any redshift from simulations with different cosmologies, astrophysics models, subgrid physics, volumes, and resolutions, indicating that the network may have found a universal relation. We then use different methods to find equations that approximate the relation found by the networks and derive new analytic expressions that predict the total mass of a subhalo from its radius, velocity dispersion, and maximum circular velocity. We show that in some regimes, the analytic expressions are more accurate than the neural networks. We interpret the relation found by the neural network and approximated by the analytic equation as being connected to the virial theorem.

preprint2021arXiv

The CAMELS Multifield Dataset: Learning the Universe&#39;s Fundamental Parameters with Artificial Intelligence

We present the Cosmology and Astrophysics with MachinE Learning Simulations (CAMELS) Multifield Dataset, CMD, a collection of hundreds of thousands of 2D maps and 3D grids containing many different properties of cosmic gas, dark matter, and stars from 2,000 distinct simulated universes at several cosmic times. The 2D maps and 3D grids represent cosmic regions that span $\sim$100 million light years and have been generated from thousands of state-of-the-art hydrodynamic and gravity-only N-body simulations from the CAMELS project. Designed to train machine learning models, CMD is the largest dataset of its kind containing more than 70 Terabytes of data. In this paper we describe CMD in detail and outline a few of its applications. We focus our attention on one such task, parameter inference, formulating the problems we face as a challenge to the community. We release all data and provide further technical details at https://camels-multifield-dataset.readthedocs.io.

preprint2021arXiv

The GIGANTES dataset: precision cosmology from voids in the machine learning era

We present GIGANTES, the most extensive and realistic void catalog suite ever released -- containing over 1 billion cosmic voids covering a volume larger than the observable Universe, more than 20 TB of data, and created by running the void finder VIDE on QUIJOTE&#39;s halo simulations. The expansive and detailed GIGANTES suite, spanning thousands of cosmological models, opens up the study of voids, answering compelling questions: Do voids carry unique cosmological information? How is this information correlated with galaxy information? Leveraging the large number of voids in the GIGANTES suite, our Fisher constraints demonstrate voids contain additional information, critically tightening constraints on cosmological parameters. We use traditional void summary statistics (void size function, void density profile) and the void auto-correlation function, which independently yields an error of $0.13\,\mathrm{eV}$ on $\sum\,m_ν$ for a 1 $h^{-3}\mathrm{Gpc}^3$ simulation, without CMB priors. Combining halos and voids we forecast an error of $0.09\,\mathrm{eV}$ from the same volume. Extrapolating to next generation multi-Gpc$^3$ surveys such as DESI, Euclid, SPHEREx, and the Roman Space Telescope, we expect voids should yield an independent determination of neutrino mass. Crucially, GIGANTES is the first void catalog suite expressly built for intensive machine learning exploration. We illustrate this by training a neural network to perform likelihood-free inference on the void size function. Cosmology problems provide an impetus to develop novel deep learning techniques, leveraging the symmetries embedded throughout the universe from physical laws, interpreting models, and accurately predicting errors. With GIGANTES, machine learning gains an impressive dataset, offering unique problems that will stimulate new techniques.

preprint2020arXiv

A new probe of Axion-Like Particles: CMB polarization distortions due to cluster magnetic fields

We propose using the upcoming Cosmic Microwave Background (CMB) ground based experiments to detect the signal of ALPs (Axion like particles) interacting with magnetic fields in galaxy clusters. The conversion between CMB photons and ALPs in the presence of the cluster magnetic field can cause a polarized spectral distortion in the CMB around a galaxy cluster. The strength of the signal depends upon the redshift of the galaxy cluster and will exhibit a distinctive spatial profile around it depending upon the structure of electron density and magnetic field. This distortion produces a different shape from the other known spectral distortions like $y$-type and $μ$-type and hence are separable from the multi-frequency CMB observation. The spectrum is close to kinematic Sunyaev-Zeldovich (kSZ) signal but can be separated from it using the polarization information. For the future ground-based CMB experiments such as Simons Observatory and CMB-S4, we estimate the measurability of this signal in the presence of foreground contamination, instrument noise and CMB anisotropies. This new avenue can probe the photon-ALP coupling over the ALP mass range from $10^{-13}$ eV to $10^{-12}$ eV with two orders of magnitude better accuracy from CMB-S4 than the current existing bounds.

preprint2020arXiv

Limitations to the &#34;basic&#39;&#39; HOD model and beyond

We make use of the IllustrisTNG cosmological, hydrodynamical simulations to test fundamental assumptions of the mass-based Halo Occupation Distribution (HOD) approach to modelling the galaxy-halo connection. By comparing the clustering of galaxies measured in the 300 Mpc TNG box (TNG300) with that predicted by the standard (``basic&#39;&#39;) HOD model, we find that, on average, the ``basic&#39;&#39; HOD model underpredicts the real-space correlation function in the TNG300 box by $\sim$ 15\% on scales of $1 \ {\rm Mpc}/h < r < 20 \ {\rm Mpc}/h$, which is well beyond the target precision demanded of next-generation galaxy redshift surveys. We perform several tests to establish the robustness of our findings to systematic effects, including the effect of finite box size and the choice of halo finder. In our exploration of ``secondary&#39;&#39; parameters with which to augment the ``basic&#39;&#39; HOD, we find that the local environment of the halo, the velocity dispersion anisotropy, $β$, and the product of the half-mass radius and the velocity dispersion, $σ^2 R_{\rm halfmass}$, are the three most effective measures of assembly bias that help reconcile the ``basic&#39;&#39; HOD-predicted clustering with that in TNG300. In addition, we test other halo properties such as halo spin, formation epoch and halo concentration. We also find that at fixed halo mass, galaxies in one type of environment cluster differently from galaxies in another. We demonstrate that a more complete model of the galaxy-halo connection can be constructed if we combine both mass and local environment information about the halo.

preprint2020arXiv

The Atacama Cosmology Telescope: Component-separated maps of CMB temperature and the thermal Sunyaev-Zel&#39;dovich effect

Optimal analyses of many signals in the cosmic microwave background (CMB) require map-level extraction of individual components in the microwave sky, rather than measurements at the power spectrum level alone. To date, nearly all map-level component separation in CMB analyses has been performed exclusively using satellite data. In this paper, we implement a component separation method based on the internal linear combination (ILC) approach which we have designed to optimally account for the anisotropic noise (in the 2D Fourier domain) often found in ground-based CMB experiments. Using this method, we combine multi-frequency data from the Planck satellite and the Atacama Cosmology Telescope Polarimeter (ACTPol) to construct the first wide-area, arcminute-resolution component-separated maps (covering approximately 2100 sq. deg.) of the CMB temperature anisotropy and the thermal Sunyaev-Zel&#39;dovich (tSZ) effect sourced by the inverse-Compton scattering of CMB photons off hot, ionized gas. Our ILC pipeline allows for explicit deprojection of various contaminating signals, including a modified blackbody approximation of the cosmic infrared background (CIB) spectral energy distribution. The cleaned CMB maps will be a useful resource for CMB lensing reconstruction, kinematic SZ cross-correlations, and primordial non-Gaussianity studies. The tSZ maps will be used to study the pressure profiles of galaxies, groups, and clusters through cross-correlations with halo catalogs, with dust contamination controlled via CIB deprojection. The data products described in this paper are available on LAMBDA.

preprint2020arXiv

The Atacama Cosmology Telescope: Constraints on Cosmic Birefringence

We present new constraints on anisotropic birefringence of the cosmic microwave background polarization using two seasons of data from the Atacama Cosmology Telescope covering $456$ square degrees of sky. The birefringence power spectrum, measured using a curved-sky quadratic estimator, is consistent with zero. Our results provide the tightest current constraint on birefringence over a range of angular scales between $5$ arcminutes and $9$ degrees. We improve previous upper limits on the amplitude of a scale-invariant birefringence power spectrum by a factor of between $2$ and $3$. Assuming a nearly-massless axion field during inflation, our result is equivalent to a $2\,σ$ upper limit on the Chern-Simons coupling constant between axions and photons of $g_{αγ}<4.0\times 10^{-2}/H_I$ where $H_I$ is the inflationary Hubble scale.

preprint2020arXiv

The Effective Halo Model: Creating a Physical and Accurate Model of the Matter Power Spectrum and Cluster Counts

We introduce a physically-motivated model of the matter power spectrum, based on the halo model and perturbation theory. This model achieves 1\% accuracy on all $k-$scales between $k=0.02h\,\mathrm{Mpc}^{-1}$ to $k=1h\,\mathrm{Mpc}^{-1}$. Our key ansatz is that the number density of halos depends on the non-linear density contrast filtered on some unknown scale $R$. Using the Effective Field Theory of Large Scale Structure to evaluate the two-halo term, we obtain a model for the power spectrum with only two fitting parameters: $R$ and the effective `sound speed&#39;, which encapsulates small-scale physics. This is tested with two suites of cosmological simulations across a broad range of cosmologies and found to be highly accurate. Due to its physical motivation, the statistics can be easily extended beyond the power spectrum; we additionally derive the one-loop covariance matrices of cluster counts and their combination with the matter power spectrum. This yields a significantly better fit to simulations than previous models, and includes a new model for super-sample effects, which is rigorously tested with separate universe simulations. At low redshift, we find a significant ($\sim 10\%$) exclusion covariance from accounting for the finite size of halos which has not previously been modeled. Such power spectrum and covariance models will enable joint analysis of upcoming large-scale structure surveys, gravitational lensing surveys and cosmic microwave background maps on scales down to the non-linear scale. We provide a publicly released Python code.

preprint2020arXiv

Ultra-light Dark Matter is Incompatible with the Milky Way&#39;s Dwarf Satellites

The density profiles of dwarf galaxies are a highly varied set. If the dark matter is an ultra-light particle such as axions, then simulations predict a distinctive and unique profile. If the axion mass is large enough to fit the ultra-faint dwarf (UFD) satellites($m\gtrapprox 10^{-21}$ eV), then the models do not fit the density profile of Fornax and Sculptor and are ruled out by more than $3-σ$ confidence. If the axion mass is in the mass range that can fit mass profiles of Fornax and Sculptor dwarf spheroidals, then its extended profile implies enormous masses ($\approx10^{11}-10^{12}M_{\odot}$) for the UFDs. These large masses for the UFDS are ruled out by more than $3-σ$ confidence by dynamical friction arguments. The tension would increase further considering star formation histories and stellar masses of the UFDs. Unless future ultra-light dark matter (ULDM) simulations with baryonic feedback show a significant change in the density structure of the halos, the current data is incompatible with the ULDM scenario. Relaxing the slope constraint from classical dwarf galaxies would lead to excluding ULDM with mass less than $6\times10^{-22}$ eV.

preprint2020arXiv

Using the Marked Power Spectrum to Detect the Signature of Neutrinos in Large-Scale Structure

Cosmological neutrinos have their greatest influence in voids: these are the regions with the highest neutrino to dark matter density ratios. The marked power spectrum can be used to emphasize low density regions over high density regions, and therefore is potentially much more sensitive than the power spectrum to the effects of neutrino masses. Using 22,000 N-body simulations from the Quijote suite, we quantify the information content in the marked power spectrum of the matter field, and show that it outperforms the standard power spectrum by setting constraints improved by a factor larger than 2 on all cosmological parameters. The combination of marked and standard power spectrum allows to place a 4.3σ constraint on the minimum sum of the neutrino masses with a volume equal to 1 (Gpc/h)^3 and without CMB priors. Combinations of different marked power spectra yield a 6σ constraint within the same conditions.

preprint2020arXiv

What does the Marked Power Spectrum Measure? Insights from Perturbation Theory

The marked power spectrum is capable of placing far tighter constraints on cosmological parameters (particularly the neutrino mass) than the conventional power spectrum. What new information does it contain beyond conventional statistics? Through the development of a perturbative model, we find that the mark induces a significant coupling between non-Gaussianities, which are usually found on small scales, and large scales, leading to the additional information content. The model is derived in the context of one-loop perturbation theory and validated by comparison to $N$-body simulations across a variety of mark parameters. At moderate redshifts, including for massive neutrino cosmologies, the theory is in good agreement with the simulations. The importance of non-linear gravitational effects on the large-scale spectra complicates the modeling as there is no well-defined convergence radius of the theory at low $z$. Extension to higher perturbative order and biased tracers is possible via a similar approach, and a simple model of the latter is shown to yield promising results. The theory becomes non-perturbative at redshift zero for small smoothing scales, with important contributions from higher-order terms: these will need to be studied before the full power of this tool can be realized.

preprint2019arXiv

Dynamical Friction in a Fuzzy Dark Matter Universe

We present an in-depth exploration of the phenomenon of dynamical friction in a universe where the dark matter is composed entirely of so-called Fuzzy Dark Matter (FDM), ultralight bosons of mass $m\sim\mathcal{O}(10^{-22})\,$eV. We review the classical treatment of dynamical friction before presenting analytic results in the case of FDM for point masses, extended mass distributions, and FDM backgrounds with finite velocity dispersion. We then test these results against a large suite of fully non-linear simulations that allow us to assess the regime of applicability of the analytic results. We apply these results to a variety of astrophysical problems of interest, including infalling satellites in a galactic dark matter background, and determine that \emph{(1)}~for FDM masses $m\gtrsim 10^{-21}\, {\rm eV}\, c^{-2}$, the timing problem of the Fornax dwarf spheroidal&#39;s globular clusters is no longer solved and \emph{(2)}~the effects of FDM on the process of dynamical friction for satellites of total mass $M$ and relative velocity $v_{\rm rel}$ should require detailed numerical simulations for $\left(M/10^9~M_{\odot}\right) \left(m/10^{-22}~{\rm eV}\right)\left(100~{\rm km}~{\rm s}^{-1}/v_{\rm rel}\right) \sim 1$, parameters which would lie outside the validated range of applicability of any currently developed analytic theory, due to transient wave structures in the time-dependent regime.

preprint2018arXiv

Measuring the Duration of Last Scattering

The cosmic microwave background (CMB) fluctuations effectively measure the basic properties of the universe during the recombination epoch. CMB measurements fix the distance to the surface of last scatter, the sound horizon of the baryon-photon fluid and the fraction of the energy density in relativistic species. We show that the microwave background observations can also very effectively constrain the thickness of the last scattering surface, which is directly related to the ratio of the small-scale E-mode polarization signal to the small-scale temperature signal. The current cosmological data enables a 0.1\% measurement of the thickness of the surface of last scatter: $19 \pm 0.065$ Mpc. This constraint is relatively model-independent, so it can provide a new metric for systematic errors and an independent test of the $Λ{\rm CDM}$ model. On the other hand, it is sensitive to models which affect the reionization history of the universe such as models with annihilating dark matter and varying fundamental constants (e.g., the fine-structure constant, $α_{\rm EM}$, and electron rest mass, $m_{\rm e}$) and as such can be used as a viable tool to constrain them.