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") HOD model (bHOD), we find that, on average, the bHOD 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" parameters with which to augment the bHOD, we find that the local environment of the halo, the velocity dispersion anisotropy, $\beta$, and the product of the half-mass radius and the velocity dispersion, $\sigma^2 R_{\rm halfmass}$, are the three most effective measures of assembly bias that help reconcile the bHOD-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 both the mass and information regarding the local environment in which a halo is embedded are combined.

Over the past decade, the gravitational lensing of the Cosmic Microwave Background (CMB) has become a powerful tool for probing the matter distribution in the Universe. The standard technique used to reconstruct the CMB lensing signal employs the quadratic estimator (QE) method, which has recently been shown to be suboptimal for lensing measurements on very small scales in temperature and polarization data. We implement a simple, more optimal method for the small-scale regime, which involves taking the direct inverse of the background gradient. We derive new techniques to make continuous maps of lensing using this "Gradient-Inversion" (GI) method and validate our method with simulated data, finding good agreement with predictions. For idealized simulations of lensing cross- and autospectra that neglect foregrounds, we demonstrate that our method performs significantly better than previous quadratic estimator methods in temperature; at $L=5000-9000$, it reduces errors on the lensing auto-power spectrum by a factor of $\sim 4$ for both idealized CMB-S4 and Simons Observatory-like experiments and by a factor of $\sim 2.6$ for cross-correlations of CMB-S4-like lensing reconstruction and the true lensing field. We caution that the level of the neglected small-scale foreground power, while low in polarization, is very high in temperature; though we briefly outline foreground mitigation methods, further work on this topic is required. Nevertheless, our results show the future potential for improved small-scale CMB lensing measurements, which could provide stronger constraints on cosmological parameters and astrophysics at high redshifts.

Of the many proposed extensions to the $\Lambda$CDM paradigm, a model in which neutrinos self-interact until close to the epoch of matter-radiation equality has been shown to provide a good fit to current cosmic microwave background (CMB) data, while at the same time alleviating tensions with late-time measurements of the expansion rate and matter fluctuation amplitude. Interestingly, CMB fits to this model either pick out a specific large value of the neutrino interaction strength, or are consistent with the extremely weak neutrino interaction found in $\Lambda$CDM, resulting in a bimodal posterior distribution for the neutrino self-interaction cross section. In this paper, we explore why current cosmological data select this particular large neutrino self-interaction strength, and by consequence, disfavor intermediate values of the self-interaction cross section. We show how it is the $\ell \gtrsim 1000$ CMB temperature anisotropies, most recently measured by the Planck satellite, that produce this bimodality. We also establish that smaller scale temperature data, and improved polarization data measuring the temperature-polarization cross-correlation, will best constrain the neutrino self-interaction strength. We forecast that the upcoming Simons Observatory should be capable of distinguishing between the models.

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 $\Lambda {\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, $\alpha_{\rm EM}$, and electron rest mass, $m_{\rm e}$) and as such can be used as a viable tool to constrain them.

Calculations of the Cosmic Microwave Background lensing power implemented into the standard cosmological codes such as CAMB and CLASS usually treat the surface of last scatter as an infinitely thin screen. However, since the CMB anisotropies are smoothed out on scales smaller than the diffusion length due to the effect of Silk damping, the photons which carry information about the small-scale density distribution come from slightly earlier times than the standard recombination time. The dominant effect is the scale dependence of the mean redshift associated with the fluctuations during recombination. We find that fluctuations at $k = 0.01 {\rm \ Mpc^{-1}}$ come from a characteristic redshift of $z \approx 1090$, while fluctuations at $k = 0.3 {\rm \ Mpc^{-1}}$ come from a characteristic redshift of $z \approx 1130$. We then estimate the corrections to the lensing kernel and the related power spectra due to this effect. We conclude that neglecting it would result in a deviation from the true value of the lensing kernel at the half percent level at small CMB scales. For an all-sky, noise-free experiment, this corresponds to a $\sim 0.1 \sigma$ shift in the observed temperature power spectrum on small scales ($ 2500 \lesssim l \lesssim 4000$).

We present a method to recover and study the projected gravitational tidal forces from a galaxy survey containing little or no redshift information. The method and the physical interpretation of the recovered tidal maps as a tracer of the cosmic web are described in detail. We first apply the method to a simulated galaxy survey and study the accuracy with which the cosmic web can be recovered in the presence of different observational effects, showing that the projected tidal field can be estimated with reasonable precision over large regions of the sky. We then apply our method to the 2MASS survey and present a publicly available full-sky map of the projected tidal forces in the local Universe. As an example of an application of these data we further study the distribution of galaxy luminosities across the different elements of the cosmic web, finding that, while more luminous objects are found preferentially in the most dense environments, there is no further segregation by tidal environment.

We present the first attempt to use a combination of CMB, LIGO, and PPTA data to constrain both the tilt and the running of primordial tensor power spectrum through constraints on the gravitational wave energy density generated in the early universe. Combining measurements at different cosmological scales highlights how complementary data can be used to test the predictions of early universe models including the inflationary consistency relation. Current data prefers a slightly positive tilt ($n_t = 0.06^{+0.63}_{-0.89}$) and a negative running ($n_{t, {\rm run}} < -0.22$) for the tensor power spectrum spectrum. Interestingly, the addition of direct gravitational wave detector data alone puts strong bounds on the tensor-to-scalar ratio $r < 0.2$ since the large positive tensor tilt preferred by the \textit{Planck} temperature power spectrum is no longer allowed. Adding the recently released BICEP2/KECK and Planck 353 GHz Polarization cross correlation data puts an even stronger bound $r<0.1$ . We comment on possible effects of a large positive tilt on the background expansion and show that depending on the assumptions regarding the UV cutoff ($k_{\rm UV}/k_* = 10^{24}$) of the primordial spectrum of gravitational waves, the strongest bounds on $n_t = 0.07^{+0.52}_{-0.80}$ are derived from this effect.