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 small-scale non-Gaussianities 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 mixing of large and small scales 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. This suggests that marked spectrum may have important contributions from baryonic effects even at low $k$: these will need to be studied before the full power of this tool can be realized.

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', 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.

We present cosmological constraints from a joint analysis of the pre- and post-reconstruction galaxy power spectrum multipoles from the final data release of the Baryon Oscillation Spectroscopic Survey (BOSS). Geometric constraints are obtained from the positions of BAO peaks in reconstructed spectra, analyzed in combination with the unreconstructed spectra in a full-shape (FS) likelihood using a joint covariance matrix, giving stronger parameter constraints than FS-only or BAO-only analyses. We introduce a new method for obtaining constraints from reconstructed spectra based on a correlated theoretical error, which is shown to be simple, robust, and applicable to any flavor of density-field reconstruction. Assuming $\Lambda$CDM with massive neutrinos, we analyze data from two redshift bins $z_\mathrm{eff}=0.38,0.61$ and obtain $1.6\%$ constraints on the Hubble constant $H_0$, using only a single prior on the current baryon density $\omega_b$ from Big Bang Nucleosynthesis (BBN) and no knowledge of the power spectrum slope $n_s$. This gives $H_0 = 68.6\pm1.1\,\mathrm{km\,s}^{-1}\mathrm{Mpc}^{-1}$, with the inclusion of BAO data sharpening the measurement by $40\%$, representing one of the strongest current constraints on $H_0$ independent of cosmic microwave background data. Restricting to the best-fit slope $n_s$ from Planck (but without additional priors on the spectral shape), we obtain a $1\%$ $H_0$ measurement of $67.8\pm 0.7\,\mathrm{km\,s}^{-1}\mathrm{Mpc}^{-1}$. We find strong constraints on the cosmological parameters from a joint analysis of the FS, BAO, and Planck data. This sets new bounds on the sum of neutrino masses $\sum m_\nu < 0.14\,\mathrm{eV}$ (at $95\%$ confidence) and the effective number of relativistic degrees of freedom $N_\mathrm{eff} = 2.90^{+0.15}_{-0.16}$, though contours are not appreciably narrowed by the inclusion of BAO data.

We present $\mathcal{O}(N^2)$ estimators for the small-scale power spectrum and bispectrum in cosmological simulations. In combination with traditional methods, these allow spectra to be efficiently computed across a vast range of scales, requiring orders of magnitude less computation time than Fast Fourier Transform based approaches alone. These methods are applicable to any tracer; simulation particles, halos or galaxies, and take advantage of the simple geometry of the box and periodicity to remove almost all dependence on large random particle catalogs. By working in configuration-space, both power spectra and bispectra can be computed via a weighted sum of particle pairs up to some radius, which can be reduced at larger $k$, leading to algorithms with decreasing complexity on small scales. These do not suffer from aliasing or shot-noise, allowing spectra to be computed to arbitrarily large wavenumbers. The estimators are rigorously derived and tested against simulations, and their covariances discussed. The accompanying code, HIPSTER, has been publicly released, incorporating these algorithms. Such estimators will be of great use in the analysis of large sets of high-resolution simulations.

We present a new class of estimators for computing small-scale power spectra and bispectra in configuration-space via weighted pair- and triple-counts, with no explicit use of Fourier transforms. Particle counts are truncated at $R_0\sim 100h^{-1}\,\mathrm{Mpc}$ via a continuous window function, which has negligible effect on the measured power spectrum multipoles at small scales. This gives a power spectrum algorithm with complexity $\mathcal{O}(NnR_0^3)$ (or $\mathcal{O}(Nn^2R_0^6)$ for the bispectrum), measuring $N$ galaxies with number density $n$. Our estimators are corrected for the survey geometry and have neither self-count contributions nor discretization artifacts, making them ideal for high-$k$ analysis. Unlike conventional Fourier transform based approaches, our algorithm becomes more efficient on small scales (since a smaller $R_0$ may be used), thus we may efficiently estimate spectra across $k$-space by coupling this method with standard techniques. We demonstrate the utility of the publicly available power spectrum algorithm by applying it to BOSS DR12 simulations to compute the high-$k$ power spectrum and its covariance. In addition, we derive a theoretical rescaled-Gaussian covariance matrix, which incorporates the survey geometry and is found to be in good agreement with that from mocks. Computing configuration- and Fourier-space statistics in the same manner allows us to consider joint analyses, which can place stronger bounds on cosmological parameters; to this end we also discuss the cross-covariance between the two-point correlation function and the small-scale power spectrum.

To make use of clustering statistics from large cosmological surveys, accurate and precise covariance matrices are needed. We present a new code to estimate large scale galaxy two-point correlation function (2PCF) covariances in arbitrary survey geometries that, due to new sampling techniques, runs $\sim 10^4$ times faster than previous codes, computing finely-binned covariance matrices with negligible noise in less than 100 CPU-hours. As in previous works, non-Gaussianity is approximated via a small rescaling of shot-noise in the theoretical model, calibrated by comparing jackknife survey covariances to an associated jackknife model. The flexible code, RascalC, has been publicly released, and automatically takes care of all necessary pre- and post-processing, requiring only a single input dataset (without a prior 2PCF model). Deviations between large scale model covariances from a mock survey and those from a large suite of mocks are found to be be indistinguishable from noise. In addition, the choice of input mock are shown to be irrelevant for desired noise levels below $\sim 10^5$ mocks. Coupled with its generalization to multi-tracer data-sets, this shows the algorithm to be an excellent tool for analysis, reducing the need for large numbers of mock simulations to be computed.

Constraining parameters such as the initial mass function high-mass slope and the frequency of type Ia supernovae is of critical importance in the ongoing quest to understand galactic physics and create realistic hydrodynamical simulations. In this paper, we demonstrate a method to precisely determine these using individual chemical abundances from a large set of stars, coupled with some estimate of their ages. Inference is performed via the simple chemical evolution model Chempy in a Bayesian framework, marginalizing over each star's specific interstellar medium parameters, including an element-specific `model error' parameter to account for inadequacies in our model. Hamiltonian Monte Carlo (HMC) methods are used to sample the posterior function, made possible by replacing Chempy with a trained neural network at negligible error. The approach is tested using data from both Chempy and the IllustrisTNG simulation, showing sub-percent agreement between inferred and true parameters using data from up to 1600 individual stellar abundances. For IllustrisTNG, strongest constraints are obtained from metal ratios, competitive with those from other methods including star counts. Analysis using a different set of nucleosynthetic yields shows that incorrectly assumed yield models can give non-negligible bias in the derived parameters; this is reduced by our model errors, which further show how well the yield tables match data. We also find a significant bias from analyzing only a small set of stars, as is often done in current analyses. The method can be easily applied to observational data, giving tight bounds on key galactic parameters from chemical abundances alone.

We present configuration-space estimators for the auto- and cross-covariance of two- and three-point correlation functions (2PCF and 3PCF) in general survey geometries. These are derived in the Gaussian limit (setting higher-order correlation functions to zero), but for arbitrary non-linear 2PCFs (which may be estimated from the survey itself), with a shot-noise rescaling parameter included to capture non-Gaussianity. We generalize previous approaches to include Legendre moments via a geometry-correction function calibrated from measured pair and triple counts. Making use of importance sampling and random particle catalogs, we can estimate model covariances in fractions of the time required to do so with mocks, obtaining estimates with negligible sampling noise in $\sim 10$ ($\sim 100$) CPU-hours for the 2PCF (3PCF) auto-covariance. We compare results to sample covariances from a suite of BOSS DR12 mocks and find the matrices to be in good agreement, assuming a shot-noise rescaling parameter of $1.03$ ($1.20$) for the 2PCF (3PCF). To obtain strongest constraints on cosmological parameters we must use multiple statistics in concert; having robust methods to measure their covariances at low computational cost is thus of great relevance to upcoming surveys.

Searches for inflationary gravitational wave signals in the CMB B-mode polarisation are expected to reach unprecedented power over the next decade. A major difficulty in these ongoing searches is that galactic foregrounds such as dust can easily mimic inflationary signals. Though typically foregrounds are separated from primordial signals using the foregrounds' different frequency dependence, in this paper we investigate instead the extent to which the galactic dust B-modes' statistical anisotropy can be used to distinguish them from inflationary B-modes, building on the work of Kamionkowski and Kovetz (2014). In our work, we extend existing anisotropy estimators and apply them to simulations of polarised dust to forecast their performance for future experiments. Considering the application of this method as a null-test for dust contamination to CMB-S4, we find that we can detect residual dust levels corresponding to $r\sim0.001$ at $2\sigma$, which implies that statistical anisotropy estimators will be a powerful diagnostic for foreground residuals (though our results show some dependence on the dust simulation used). Finally, considering applications beyond a simple null test, we demonstrate how anisotropy statistics can be used to construct an estimate of the dust B-mode map, which could potentially be used to clean the B-mode sky.

To fully harvest the rich library of stellar elemental abundance data available, we require reliable models that facilitate our interpretation of them. Galactic chemical evolution (GCE) models are one such set, and a key part of which are the selection of chemical yields from different nucleosynthetic enrichment channels, predominantly asymptotic giant branch (AGB) stars, Type Ia supernovae (SNe Ia), and core-collapse supernovae (CC-SNe). Here, we present a scoring system for yield tables based on their ability to reproduce proto-solar abundances within a simple parametrisation of the GCE modelling software Chempy, which marginalises over galactic parameters describing simple stellar populations (SSPs) and interstellar medium physics. Two statistical scoring methods are presented, based on Bayesian evidence and leave-one-out cross-validation and are applied to five CC-SN tables, (a) for all mutually available elements and (b) for a subset of the 9 most abundant elements. We find that the yields of Prantzos et al. (P18, including stellar rotation) and Chieffi & Limongi (C04) best reproduce proto-solar abundances for the two cases, respectively. The inferred best-fit SSP parameters for (b) are $\alpha_\mathrm{IMF}=-2.45^{+0.15}_{-0.11}$ for the initial mass function high-mass slope and $\mathrm{N}_\mathrm{Ia}=1.29^{+0.45}_{-0.31}\times10^{-3}$ M$_\odot^{-1}$ for the SN Ia normalisation, which are broadly consistent across tested yield tables. Additionally, we demonstrate how Chempy can be used to dramatically improve elemental abundance predictions of hydrodynamical simulations by plugging tailored best-fit SSP parameters into a Milky Way analogue from Gutcke & Springel. Our code, including a comprehensive tutorial, is freely available and can additionally provide SSP enrichment tables for any combination of parameters and yield tables.