We explore the impact of elastic, anisotropic, velocity-dependent dark matter (DM) self-interactions on the host halo and subhalos of Milky Way (MW)--mass systems. We consider a generic self-interacting dark matter (SIDM) model parameterized by the masses of a light mediator and the DM particle. The ratio of these masses, $w$, sets the velocity scale above which momentum transfer due to DM self-interactions becomes inefficient. We perform high-resolution zoom-in simulations of an MW-mass halo for values of $w$ that span scenarios in which self-interactions either between the host and its subhalos or only within subhalos efficiently transfer momentum, and we study the effects of self-interactions on the host halo and on the abundance, radial distribution, orbital dynamics, and density profiles of subhalos in each case. The abundance and properties of surviving subhalos are consistent with being determined primarily by subhalo--host halo interactions. In particular, subhalos on radial orbits in models with larger values of the cross section at the host halo velocity scale are more susceptible to tidal disruption owing to mass loss from ram pressure stripping caused by self-interactions with the host. This mechanism suppresses the abundance of surviving subhalos relative to collisionless DM simulations, with stronger suppression for larger values of $w$. Thus, probes of subhalo abundance around MW-mass hosts can be used to place upper limits on the self-interaction cross section at velocity scales of $\sim 200\ \rm{km\ s}^{-1}$, and combining these measurements with the orbital properties and internal dynamics of subhalos may break degeneracies among velocity-dependent SIDM models.

We study the covariance in the angular power spectrum estimates of CMB fluctuations when the primordial fluctuations are non-Gaussian. The non-Gaussian covariance comes from a nonzero connected four-point correlation function -- or the trispectrum in Fourier space -- and can be large when long-wavelength (super-CMB) modes are strongly coupled to short-wavelength modes. The effect of such non-Gaussian covariance can be modeled through additional freedom in the theoretical CMB angular power spectrum and can lead to different inferred values of the standard cosmological parameters relative to those in $\Lambda$CDM. Taking the collapsed limit of the primordial trispectrum in the quasi-single field inflation model as an example, we study how the six standard $\Lambda$CDM parameters shift when two additional parameters describing the trispectrum are allowed. The reduced statistical significance of the Hubble tension in the extended model allows us to combine the {\it Planck} temperature data and the type Ia supernovae data from Panstarrs with the distance-ladder measurement of the Hubble constant. This combination of data shows strong evidence for a primordial trispectrum-induced non-Gaussian covariance, with a likelihood improvement of $\Delta \chi^2 \approx -15$ (with two additional parameters) relative to $\Lambda$CDM.

Gravitational waves produced from the merger of binary neutron stars (BNSs) are accompanied by electromagnetic counterparts, making it possible to identify the associated host galaxy. We explore how properties of the host galaxies relate to the astrophysical processes leading to the mergers. It is thought that the BNS merger rate within a galaxy at a given epoch depends primarily on the galaxy's star-formation history as well as the underlying merger time-delay distribution of the binary systems. The stellar history of a galaxy, meanwhile, depends on the cosmological evolution of the galaxy through time, and is tied to the growth of structure in the Universe. We study the hosts of BNS mergers in the context of structure formation by populating the Universe Machine simulations with gravitational-wave events~ according to a simple time-delay model. We find that different time-delay distributions predict different properties of the associated host galaxies, including the distributions of stellar mass, star-formation rate, halo mass, and local and large-scale clustering of hosts. BNSs that merge today with short delay times prefer to be in hosts that have high star-formation rates, while those with long delay times live in dense regions within massive halos that have low star formation. We show that with ${\mathcal O}(10)$ events from current gravitational-wave detector networks, it is possible to make preliminary distinctions between formation channels which trace stellar mass, halo mass, or star-formation rate. We also find that strategies to follow up gravitational-wave events with electromagnetic telescopes can be significantly optimized using the clustering properties of their hosts.

The splashback radius ($R_{\rm sp}$) of dark matter halos has recently been detected using weak gravitational lensing and cross-correlations with galaxies. However, different methods have been used to measure $R_{\rm sp}$ and to assess the significance of its detection. In this paper, we use simulations to study the precision and accuracy to which we can detect the splashback radius with 3D density, 3D subhalo, and weak lensing profiles. We study how well various methods and tracers recover $R_{\rm sp}$ by comparing it with the value measured directly from particle dynamics. We show that estimates of $R_{\rm sp}$ from density and subhalo profiles correspond to different percentiles of the underlying $R_{\rm sp}$ distribution of particle orbits. At low accretion rates, a second caustic appears and can bias results. Finally, we show that upcoming lensing surveys may be able to constrain the splashback-accretion rate relation directly.

Non-gravitational interactions between dark matter particles with strong scattering, but relatively small annihilation and dissipation, has been proposed to match various observables on cluster and group scales. In this paper, we present the results from large cosmological simulations which include the effects of different self-interaction scenarios. In particular we explore a model with the differential cross section that can depend on both the relative velocity of the interacting particles and the angle of scattering. We focus on how quantities, such as the stacked density profiles, subhalo counts and the splashback radius change as a function of different forms of self-interaction. We find that self-interactions not only affect the central region of the cluster, the effect well known from previous studies, but also significantly alter the distribution of subhalos and the density of particles out to the splashback radius. Our results suggest that current weak lensing data can already put constraints on the self-interaction cross-section that are only slightly weaker than the Bullet Cluster constraints ($\sigma/m \lesssim 2$ cm$^2/$g), and future lensing surveys should be able to tighten them even further making halo profiles on cluster scales a competitive probe for DM physics.

We present a detection of the splashback feature around galaxy clusters selected using their Sunyaev-Zel'dovich (SZ) signal. Recent measurements of the splashback feature around optically selected galaxy clusters have found that the splashback radius, $r_{\rm sp}$, is smaller than predicted by N-body simulations. A possible explanation for this discrepancy is that $r_{\rm sp}$ inferred from the observed radial distribution of galaxies is affected by selection effects related to the optical cluster-finding algorithms. We test this possibility by measuring the splashback feature in clusters selected via the SZ effect in data from the South Pole Telescope SZ survey and the Atacama Cosmology Telescope Polarimeter survey. The measurement is accomplished by correlating these clusters with galaxies detected in the Dark Energy Survey Year 3 data. The SZ observable used to select clusters in this analysis is expected to have a tighter correlation with halo mass and to be more immune to projection effects and aperture-induced biases than optically selected clusters. We find that the measured $r_{\rm sp}$ for SZ-selected clusters is consistent with the expectations from simulations, although the small number of SZ-selected clusters makes a precise comparison difficult. In agreement with previous work, when using optically selected redMaPPer clusters, $r_{\rm sp}$ is $\sim$ $2\sigma$ smaller than in the simulations. These results motivate detailed investigations of selection biases in optically selected cluster catalogs and exploration of the splashback feature around larger samples of SZ-selected clusters. Additionally, we investigate trends in the galaxy profile and splashback feature as a function of galaxy color, finding that blue galaxies have profiles close to a power law with no discernible splashback feature, which is consistent with them being on their first infall into the cluster.

Astrophysical and cosmological observations currently provide the only robust, empirical measurements of dark matter. Future observations with Large Synoptic Survey Telescope (LSST) will provide necessary guidance for the experimental dark matter program. This white paper represents a community effort to summarize the science case for studying the fundamental physics of dark matter with LSST. We discuss how LSST will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational couplings to the Standard Model, and compact object abundances. Additionally, we discuss the ways that LSST will complement other experiments to strengthen our understanding of the fundamental characteristics of dark matter. More information on the LSST dark matter effort can be found at https://lsstdarkmatter.github.io/ .

Astrophysical observations currently provide the only robust, empirical measurements of dark matter. In the coming decade, astrophysical observations will guide other experimental efforts, while simultaneously probing unique regions of dark matter parameter space. This white paper summarizes astrophysical observations that can constrain the fundamental physics of dark matter in the era of LSST. We describe how astrophysical observations will inform our understanding of the fundamental properties of dark matter, such as particle mass, self-interaction strength, non-gravitational interactions with the Standard Model, and compact object abundances. Additionally, we highlight theoretical work and experimental/observational facilities that will complement LSST to strengthen our understanding of the fundamental characteristics of dark matter.

Splashback refers to the process of matter that is accreting onto a dark matter halo reaching its first orbital apocenter and turning around in its orbit. The cluster-centric radius at which this process occurs, r_sp, defines a halo boundary that is connected to the dynamics of the cluster. A rapid decline in the halo profile is expected near r_sp. We measure the galaxy number density and weak lensing mass profiles around redMaPPer galaxy clusters in the first year Dark Energy Survey (DES) data. For a cluster sample with mean M_200m mass ~2.5 x 10^14 M_sun, we find strong evidence of a splashback-like steepening of the galaxy density profile and measure r_sp=1.13 +/- 0.07 Mpc/h, consistent with earlier SDSS measurements of More et al. (2016) and Baxter et al. (2017). Moreover, our weak lensing measurement demonstrates for the first time the existence of a splashback-like steepening of the matter profile of galaxy clusters. We measure r_sp=1.34 +/- 0.21 Mpc/h from the weak lensing data, in good agreement with our galaxy density measurements. For different cluster and galaxy samples, we find that consistent with LCDM simulations, r_sp scales with R_200m and does not evolve with redshift over the redshift range of 0.3--0.6. We also find that potential systematic effects associated with the redMaPPer algorithm may impact the location of r_sp. We discuss progress needed to understand the systematic uncertainties and fully exploit forthcoming data from DES and future surveys, emphasizing the importance of more realistic mock catalogs and independent cluster samples.

Several unexpected features are observed at large angular scales in the CMB temperature anisotropy measurements by both WMAP and Planck. These include the lack of both variance and correlation, alignment of the lowest multipole moments with one another, hemispherical power asymmetry, and an odd-to-even parity excess. In this work, we study the statistics of eight representative large-angle CMB features in order to evaluate their covariance in LCDM. We do so using two sets of simulated CMB temperature maps; an ensemble of 100,000 simple Gaussian simulations, and 1000 Full Focal Plane (FFP) simulations provided by the Planck collaboration. In measuring feature probabilities, we pay particular attention to analysis choices, making sure that we can reproduce previous results in the literature, and explain differences where appropriate. The covariance structure we find is consistent with expectations given that many of the features studied are functions of the angular power spectrum. Notably, we find significant differences in the covariance entries associated with the quadrupole-octopole alignments derived from the Gaussian and FFP simulations. We additionally perform a principal component analysis to quantitatively gauge what combinations of features capture the most information about how simulation measurements vary, and to provide an alternative assessment of the ways in which the real sky is anomalous. The first four principal components explain about 90\% of the simulations' variance, with the first two roughly quantifying the lack of large-angle correlations, and the next two quantifying the phase-dependent anomalies (multipole alignments and power asymmetry). Though the results of this analysis are fairly unsurprising, its comprehensive approach serves to tie together a number of previous results, and will therefore provide context for future studies of large-angle anomalies.

We study statistical anisotropies generated in the observed two-point function of the cosmic microwave background (CMB) fluctuations if the primordial statistics are non-Gaussian. Focusing on the dipole modulations of the anisotropies, we find that the hemispherical power asymmetry observed in the CMB temperature fluctuations can be modeled by a local-type trispectrum with amplitude $\tau_{\rm NL}(k_p=0.05~{\rm Mpc}^{-1}) \approx 2 \times 10^4$ and a large red tilt $n\approx -0.68$. We numerically evaluate the non-Gaussian covariance of the modulation estimators for both temperature and E-mode polarization fluctuations and discuss the prospects of constraining the model using Planck satellite data. We then discuss other effects of the scale-dependent trispectrum that could be used to distinguish this scenario from other explanations of the power asymmetry: higher-order modulations of the two-point function and the non-Gaussian angular power spectrum covariance. As an important consequence of the non-Gaussian power spectrum covariance, we discuss how the CMB-inferred spectral index of primordial scalar fluctuations can be significantly biased in the presence of a scale-dependent local-type trispectrum.

Tensions between cosmological measurements by different surveys or probes have always been important --- and are presently much discussed --- as they may lead to evidence of new physics. Several tests have been devised to probe the consistency of datasets given a cosmological model, but they often have undesired features such as dependence on the prior volume, or burdensome requirements such as that of near-Gaussian posterior distributions. We propose a new quantity, defined in a similar way as the Bayesian evidence ratio, in which these undesired properties are absent. We test the quantity on simple models with Gaussian and non-Gaussian likelihoods. We then apply it to data from the Planck satellite: we investigate the consistency of $\Lambda$CDM model parameters obtained from TT and EE angular power spectrum measurements, as well as the mutual consistency of cosmological parameters obtained from large scale (multipoles, $\ell<1000$) and small scale ($\ell \geq 1000$) portions of each measurement and find no significant discrepancy in the six-dimensional $\Lambda$CDM parameter space.

We study the kinematics of galaxies within massive clusters, as a probe of the physics of star-formation quenching within clusters. Using N-body simulations, we argue that satellite kinematics provide information about galaxy infall that is complementary to the (instantaneous) spatial distribution of satellites. Comparing the simulation results with measurements of real cluster galaxies, we find evidence that the kinematics of red (quiescent) satellite galaxies are consistent with earlier infall times than that of blue (star-forming) satellites.

Self-Interacting Dark Matter (SIDM) has long been proposed as a solution to small scale problems posed by standard Cold Dark Matter (CDM). We use numerical simulations to study the effect of dark matter interactions on the morphology of disk galaxies falling into galaxy clusters. The effective drag force on dark matter leads to offsets of the stellar disk with respect to the surrounding halo, causing distortions in the disk. For anisotropic scattering cross-sections of 0.5 and 1.0$\,\textrm{cm}^{2}\textrm{g}^{-1}$, we show that potentially observable warps, asymmetries, and thickening of the disk occur in simulations. We discuss observational tests of SIDM with galaxy surveys and more realistic simulations needed to obtain detailed predictions.

The splashback radius is a physical scale in dark matter halos that is set by the gravitational dynamics of recently accreted shells. We use analytical models and N-body simulations to study the dependence of splashback on dark energy and screened modified gravity theories. In modified gravity models, the transition from screened to unscreened regions typically occurs in the cluster outskirts, suggesting potentially observable signatures in the splashback feature. We investigate the location of splashback in both chameleon and Vainshtein screened models and find significant differences compared with $\Lambda$CDM predictions. We also find an interesting interplay between dynamical friction and modified gravity, providing a distinctive signature for modified gravity models in the behavior of the splashback feature as a function of galaxy luminosity.

Cosmological inflation predicts that the scalar spectral index "runs" with scale. Constraints on the values of the spectral runnings, $\alpha_s\equiv \textrm{d} n_s/\textrm{d}\ln k$ and $\beta_s\equiv \textrm{d}\alpha_s/\textrm{d}\ln k$, therefore provide a fundamental test of the physics of inflation. Here we study the feasibility of measuring the runnings when information from upcoming large-volume galaxy surveys is used to supplement the information provided by a CMB-S4 experiment, particularly focusing on the effect of including high-$k$, nonlinear scales. Since these measurements will be sensitive to modeling uncertainties for the nonlinear power spectrum, we examine how three different ways of parameterizing those systematics---introducing zero, two, or several hundred nuisance parameters---affect constraints and protect against parameter biases. Considering statistical errors alone, we find that including strongly nonlinear scales can substantially tighten constraints. However, these constraints weaken to levels not much better than those from a CMB-S4 experiment alone when we limit our analysis to scales where estimates are not strongly affected by systematic biases. Given these considerations, near-future large-scale structure surveys are unlikely to add much information to the CMB-S4 measurement of the first running $\alpha_s$. There is more potential for improvement for the second running, $\beta_s$, for which large-scale structure information will allow constraints to be improved by a factor of 3--4 relative to using the CMB alone. Though these constraints are still above the value predicted by slow roll inflation, they do probe regions of parameter space relevant to nonstandard inflationary models with large runnings, for example those that can generate an appreciable abundance of primordial black holes.

Mass around dark matter halos can be divided into "infalling" material and "collapsed" material that has passed through at least one pericenter. Analytical models and simulations predict a rapid drop in the halo density profile associated with the transition between these two regimes. Using data from SDSS, we explore the evidence for such a feature in the density profiles of galaxy clusters and investigate the connection between this feature and a possible phase space boundary. We first estimate the steepening of the outer galaxy density profile around clusters: the profiles show an abrupt steepening, providing evidence for truncation of the halo profile. Next, we measure the galaxy density profile around clusters using two sets of galaxies selected based on color. We find evidence of an abrupt change in the galaxy colors that coincides with the location of the steepening of the density profile. Since galaxies are likely to be quenched of star formation and turn red inside of clusters, this change in the galaxy color distribution can be interpreted as the transition from an infalling regime to a collapsed regime. We also measure this transition using a model comparison approach which has been used recently in studies of the "splashback" phenomenon, but find that this approach is not a robust way to quantify the significance of detecting a splashback-like feature. Finally, we perform measurements using an independent cluster catalog to test for potential systematic errors associated with cluster selection. We identify several avenues for future work: improved understanding of the small-scale galaxy profile, lensing measurements, identification of proxies for the halo accretion rate, and other tests. With upcoming data from the DES, KiDS and HSC surveys, we can expect significant improvements in the study of halo boundaries.

We develop and study the position-dependent bispectrum. It is a generalization of the recently proposed position-dependent power spectrum method of measuring the squeezed-limit bispectrum. The position-dependent bispectrum can similarly be used to measure the squeezed-limit trispectrum in which one of the wavelengths is much longer than the other three. In this work, we will mainly consider the case in which the three smaller wavelengths are nearly the same (the equilateral configuration). We use the Fisher information matrix to forecast constraints on bias parameters and the amplitude of primordial trispectra from the position-dependent bispectrum method. We find that the method can constrain the local-type $g_{\rm NL}$ at a level of $\sigma(g_{\rm NL}^{\rm local}) \approx 3 \times 10^5$ for a large volume SPHEREx-like survey; improvements can be expected by including all the triangular configurations of the bispectra rather than just the equilateral configuration. However, the same measurement would also constrain a much larger family of trispectra than local $g_{\rm NL}$ model. We discuss the implications of the forecasted reach of future surveys in terms of super cosmic variance uncertainties from primordial non-Gaussianities.

We present a novel method to detect the effects of dynamical friction in observed galaxy clusters. Following accretion into clusters, massive satellite galaxies will backsplash to systematically smaller radii than less massive satellites, an effect that may be detected by stacking the number density profiles of galaxies around clusters. We show that this effect may be understood using a simple toy model which reproduces the trends with halo properties observed in simulations. We search for this effect using SDSS redMaPPer clusters with richness 10<lambda<20, and find that bright (M_i<-21.5) satellites have smaller splashback radii than fainter (M_i>-20) satellites at 99% confidence.

We derive a general expression for the probability of observing deviations from statistical isotropy in the cosmic microwave background (CMB) if the primordial fluctuations are non-Gaussian and extend to superhorizon scales. The primary motivation is to properly characterize the monopole and dipole modulations of the primordial power spectrum that are generated by the coupling between superhorizon and subhorizon perturbations. Unlike previous proposals for generating the hemispherical power asymmetry, we do not assume that the power asymmetry results from a single large superhorizon mode. Instead, we extrapolate the observed power spectrum to superhorizon scales and compute the power asymmetry that would result from a specific realization of non-Gaussian perturbations on scales larger than the observable universe. Our study encompasses many of the scenarios that have been put forward as possible explanations for the CMB hemispherical power asymmetry. We confirm our analytic predictions for the probability of a given power asymmetry by comparing them to numerical realizations of CMB maps. We find that non-local models of non-Gaussianity and scale-dependent local non-Gaussianity produce scale-dependent modulations of the power spectrum, thereby potentially producing both a monopolar and a dipolar power modulation on large scales. We then provide simple examples of finding the posterior distributions for the parameters of the bispectrum from the observed monopole and dipole modulations.

We employ robust weak gravitational lensing measurements to improve cosmological constraints from measurements of the galaxy cluster mass function and its evolution, using X-ray selected clusters detected in the ROSAT All-Sky Survey. Our lensing analysis constrains the absolute mass scale of such clusters at the 8 per cent level, including both statistical and systematic uncertainties. Combining it with the survey data and X-ray follow-up observations, we find a tight constraint on a combination of the mean matter density and late-time normalization of the matter power spectrum, $\sigma_8(\Omega_m/0.3)^{0.17}=0.81\pm0.03$, with marginalized, one-dimensional constraints of $\Omega_m=0.26\pm0.03$ and $\sigma_8=0.83\pm0.04$. For these two parameters, this represents a factor of two improvement in precision with respect to previous work, primarily due to the reduced systematic uncertainty in the absolute mass calibration provided by the lensing analysis. Our new results are in good agreement with constraints from cosmic microwave background (CMB) data, both WMAP and Planck (plus WMAP polarization), under the assumption of a flat $\Lambda$CDM cosmology with minimal neutrino mass. Consequently, we find no evidence for non-minimal neutrino mass from the combination of cluster data with CMB, supernova and baryon acoustic oscillation measurements, regardless of which all-sky CMB data set is used (and independent of the recent claimed detection of B-modes on degree scales). We also present improved constraints on models of dark energy (both constant and evolving), modifications of gravity, and primordial non-Gaussianity. Assuming flatness, the constraints for a constant dark energy equation of state from the cluster data alone are at the 15 per cent level, improving to $\sim 6$ per cent when the cluster data are combined with other leading probes.

Recently, it was shown that local variance maps of temperature anisotropy are simple and useful tools for the study of large scale hemispherical power asymmetry. This was done by studying the distribution of dipoles of the local variance maps. In this work, we extend the study of the dipolar asymmetry in local variance maps using foreground cleaned Planck 143 GHz and 217 GHz data to smaller scales. In doing so, we include the effect of the CMB Doppler dipole. Further, we show that it is possible to use local variance maps to measure the Doppler dipole in these Planck channel maps, after removing large scale features (up to $l=600$), at a significance of about $3 \sigma$. At these small scales, we do not find any power asymmetry in the direction of the anomalous large scale power asymmetry beyond that expected from cosmic variance. At large scales, we verify previous results i.e. the presence of hemispherical power asymmetry at a significance of at least $3.3 \sigma$.

Recent work has shown that density profiles in the outskirts of dark matter halos can become extremely steep over a narrow range of radius. This behavior is produced by splashback material on its first apocentric passage after accretion. We show that the location of this splashback feature may be understood quite simply, from first principles. We present a simple model, based on spherical collapse, that accurately predicts the location of splashback without any free parameters. The important quantities that determine the splashback radius are accretion rate and redshift.

We propose a method to measure the ellipticities of dark matter halos using the lens-shear-shear 3-point correlation function. This method is immune to effects of galaxy-halo misalignments that can potentially limit 2-point galaxy-galaxy lensing measurements of halo anisotropy. Using a simple model for the projected mass distributions of dark matter halos, we construct an ellipticity estimator that sums over all possible triangular configurations of the 3-point function. By applying our estimator to halos from N-body simulations, we find that systematic errors in the recovered ellipticity will be at the <5% fractional level. We estimate that future imaging surveys like LSST will have sufficient statistics to detect halo ellipticities using 3-point lensing.

We perform cosmological N-body simulations with non-Gaussian initial conditions generated from two independent fields. The dominant contribution to the perturbations comes from a purely Gaussian field, but we allow the second field to have local non-Gaussianity that need not be weak. This scenario allows us to adjust the relative importance of non-Gaussian contributions beyond the skewness, producing a scaling of the higher moments different from (and stronger than) the scaling in the usual single field local ansatz. We compare semi-analytic prescriptions for the non-Gaussian mass function, large scale halo bias, and stochastic bias against the simulation results. We discuss applications of this work to large scale structure measurements that can test a wider range of models for the primordial fluctuations than is usually explored.