We explore a class of primordial power spectra that can fit the observed anisotropies in the cosmic microwave background well and that predicts a value for the Hubble parameter consistent with the local measurement of $H_0 = 74$ km/s/Mpc. This class of primordial power spectrum consists of a continuous deformation between the best-fit power law primordial power spectrum and the primordial power spectrum derived from the modified Richardson-Lucy deconvolution algorithm applied to the $C_\ell$s of best-fit power law primordial power spectrum. We find that linear interpolation half-way between the power law and modified Richardson-Lucy power spectra fits the Planck data better than the best-fit $\Lambda$CDM by $\Delta\log\mathcal{L} = 2.5$. In effect, this class of deformations of the primordial power spectra offer a new dimension which is correlated with the Hubble parameter. This correlation causes the best-fit value for $H_0$ to shift and the uncertainty to expand to $H_0 = 70.2 \pm 1.2$ km/s/Mpc. When considering the Planck dataset combined with the Cepheid $H_0$ measurement, the best-fit $H_0$ becomes $H_0 = 71.8 \pm 0.9$ km/s/Mpc. We also compute a Bayes factor of $\log K = 5.7$ in favor of the deformation model.

We determine the Hubble constant $H_0$ precisely ($2.3\%$ uncertainty) in a manner independent of cosmological model through Gaussian process regression, using strong lensing and supernova data. Strong gravitational lensing of a variable source can provide a time-delay distance $D_{\Delta t}$ and angular diameter distance to the lens $D_{\rm{d}}$. These absolute distances can anchor Type Ia supernovae, which give an excellent constraint on the shape of the distance-redshift relation. Updating our previous results to use the H0LiCOW program's milestone dataset consisting of six lenses, four of which have both $D_{\Delta t}$ and $D_{\rm{d}}$ measurements, we obtain $H_0=72.8_{-1.7}^{+1.6}\rm{\ km/s/Mpc}$ for a flat universe and $H_0=77.3_{-3.0}^{+2.2}\rm{\ km/s/Mpc}$ for a non-flat universe. We carry out several consistency checks on the data and find no statistically significant tensions, though a noticeable redshift dependence persists in a particular systematic manner that we investigate. Speculating on the possibility that this trend of derived Hubble constant with lens distance is physical, we show how this can arise through modified gravity light propagation, which would also impact the weak lensing $\sigma_8$ tension.

The extended excess towards the Galactic Center (GC) in gamma rays inferred from Fermi-LAT observations has been interpreted as being due to dark matter (DM) annihilation. Here, we perform new likelihood analyses of the GC and show that when including templates for the stellar galactic and nuclear bulges, the GC shows no significant detection of a DM annihilation template, even after generous variations in the Galactic diffuse emission models and a wide range of DM halo profiles. We include Galactic diffuse emission models with combinations of 3D inverse Compton maps, variations of interstellar gas maps, and a central source of electrons. For the DM profile, we include both spherical and ellipsoidal DM morphologies and a range of radial profiles from steep cusps to kiloparsec-sized cores, motivated in part by hydrodynamical simulations. Our derived upper limits on the dark matter annihilation flux place strong constraints on DM properties. In the case of the pure $b$-quark annihilation channel, our limits on the annihilation cross section are more stringent than those from the Milky Way dwarfs up to DM masses of $\sim$TeV, and rule out the thermal relic cross section up to $\sim$300 GeV. Better understanding of the DM profile, as well as the Fermi-LAT data at its highest energies, would further improve the sensitivity to DM properties.

We reconstruct the expansion history of the Universe using type Ia supernovae (SN Ia) in a manner independent of any cosmological model assumptions. To do so, we implement a non-parametric iterative smoothing method on the Joint Light-curve Analysis (JLA) data while exploring the SN Ia light-curve hyper-parameter space by Monte Carlo Markov Chain sampling. We test to see how the posteriors of these hyper-parameters depend on cosmology, whether using different dark energy models or reconstructions shift these posteriors. Our constraints on the SN Ia light-curve hyper-parameters from our model-independent analysis are very consistent with the constraints from using different parametrizations of the equation of state of dark energy, namely the flat $\Lambda\mathrm{CDM}$ cosmology, the Chevallier-Polarski-Linder (CPL) model, and the Phenomenologically Emergent Dark Energy (PEDE) model. This implies that the distance moduli constructed from the JLA data are mostly independent of the cosmological models. We also studied that the possibility the light-curve parameters evolve with redshift and our results show consistency with no evolution. The reconstructed expansion history of the Universe and dark energy properties also seem to be in good agreement with the expectations of the standard $\Lambda\mathrm{CDM}$ model. However, our results also indicate that the data still allow for considerable flexibility in the expansion history of the Universe.

Accurate estimation of the Hubble constant, and other cosmological parameters, from distances measured by cosmic gravitational wave sirens requires sufficient allowance for the dark energy evolution. We demonstrate how model independent statistical methods, specifically Gaussian process regression, can remove bias in the reconstruction of $H(z)$, and can be combined model independently with supernova distances. This allows stringent tests of both $H_0$ and $\Lambda$CDM, and can detect unrecognized systematics. We also quantify the redshift systematic control necessary for the use of dark sirens, showing that it must approach spectroscopic precision to avoid significant bias.

We explore the implications of a rapid appearance of dark energy between the redshifts ($z$) of one and two on the expansion rate and growth of perturbations. Using both Gaussian process regression and a parameteric model, we show that this is the preferred solution to the current set of low-redshift ($z<3$) distance measurements if $H_0=73~\rm km\,s^{-1}\,Mpc^{-1}$ to within 1\% and the high-redshift expansion history is unchanged from the $\Lambda$CDM inference by the Planck satellite. Dark energy was effectively non-existent around $z=2$, but its density is close to the $\Lambda$CDM model value today, with an equation of state greater than $-1$ at $z<0.5$. If sources of clustering other than matter are negligible, we show that this expansion history leads to slower growth of perturbations at $z<1$, compared to $\Lambda$CDM, that is measurable by upcoming surveys and can alleviate the $\sigma_8$ tension between the Planck CMB temperature and low-redshift probes of the large-scale structure.

Strongly lensed quasar systems with time delay measurements provide "time delay distances", which are a combination of three angular diameter distances and serve as powerful tools to determine the Hubble constant $H_0$. However, current results often rely on the assumption of the $\Lambda$CDM model. Here we use a model-independent method based on Gaussian process to directly constrain the value of $H_0$. By using Gaussian process regression, we can generate posterior samples of unanchored supernova distances independent of any cosmological model and anchor them with strong lens systems. The combination of a supernova sample with large statistics but no sensitivity to $H_0$ with a strong lens sample with small statistics but $H_0$ sensitivity gives a precise $H_0$ measurement without the assumption of any cosmological model. We use four well-analyzed lensing systems from the state-of-art lensing program H0LiCOW and the Pantheon supernova compilation in our analysis. Assuming the Universe is flat, we derive the constraint $H_0=72.2 \pm 2.1\,$km/s/Mpc, a precision of $2.9\%$. Allowing for cosmic curvature with a prior of $\Omega_{k}=[-0.2,0.2]$, the constraint becomes $H_0=73.0_{-3.0}^{+2.8}\,$km/s/Mpc.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

While Type Ia Supernovae (SNe Ia) are one of the most mature cosmological probes, the next era promises to be extremely exciting in the number of different ways SNe Ia are used to measure various cosmological parameters. Here we review the experiments in the 2020s that will yield orders of magnitudes more SNe Ia, and the new understandings and capabilities to constrain systematic uncertainties at a level to match these statistics. We then discuss five different cosmological probes with SNe Ia: the conventional Hubble diagram for measuring dark energy properties, the distance ladder for measuring the Hubble constant, peculiar velocities and weak lensing for measuring sigma8 and strong-lens measurements of H0 and other cosmological parameters. For each of these probes, we discuss the experiments that will provide the best measurements and also the SN Ia-related systematics that affect each one.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

During the next decade, gravitational waves will be observed from hundreds of binary inspiral events. When the redshifts of the host galaxies are known, these events can be used as `standard sirens', sensitive to the expansion rate of the Universe. Measurements of the Hubble constant $H_0$ from standard sirens can be done independently of other cosmological probes, and events occurring at $z<0.1$ will allow us to infer $H_0$ independently of cosmological models. The next generation of spectroscopic galaxy surveys will play a crucial role in reducing systematic uncertainties in $H_0$ from standard sirens, particularly for the numerous `dark sirens' which do not have an electromagnetic counterpart. In combination with large spectroscopic data sets, standard sirens with an EM counterpart are expected to constrain $H_0$ to $\sim 1-2\%$ precision within the next decade. This is competitive with the best estimates of $H_0$ obtained to date and will help illuminate the current tension between existing measurements. Information on the galaxies that host the gravitational wave events will also shed light on the origin and evolution of compact object binaries.

The hot dense environment of the early universe is known to have produced large numbers of baryons, photons, and neutrinos. These extreme conditions may have also produced other long-lived species, including new light particles (such as axions or sterile neutrinos) or gravitational waves. The gravitational effects of any such light relics can be observed through their unique imprint in the cosmic microwave background (CMB), the large-scale structure, and the primordial light element abundances, and are important in determining the initial conditions of the universe. We argue that future cosmological observations, in particular improved maps of the CMB on small angular scales, can be orders of magnitude more sensitive for probing the thermal history of the early universe than current experiments. These observations offer a unique and broad discovery space for new physics in the dark sector and beyond, even when its effects would not be visible in terrestrial experiments or in astrophysical environments. A detection of an excess light relic abundance would be a clear indication of new physics and would provide the first direct information about the universe between the times of reheating and neutrino decoupling one second later.

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.

Lack of knowledge about the background expansion history of the Universe from independent observations makes it problematic to obtain a precise and accurate estimation of the Hubble constant $H_0$ from gravitational wave standard sirens, even with electromagnetic counterpart redshifts. Simply fitting simultaneously for the matter density in a flat \lcdm\ model can reduce the precision on $H_0$ from 1\% to 5\%, while not knowing the actual background expansion model of the universe (e.g.\ form of dark energy) can introduce substantial bias in estimation of the Hubble constant. When the statistical precision is at the level of 1\% uncertainty on $H_0$, biases in non-\lcdm\ cosmologies that are consistent with current data could reach the 3$\sigma$ level. To avoid model-dependent biases, statistical techniques that are appropriately agnostic about model assumptions need to be employed.

We model the expansion history of the Universe as a Gaussian Process and find constraints on the dark energy density and its low-redshift evolution using distances inferred from the Luminous Red Galaxy (LRG) and Lyman-alpha (Ly$\alpha$) datasets of the Baryon Oscillation Spectroscopic Survey, supernova data from the Joint Light-curve Analysis (JLA) sample, Cosmic Microwave Background (CMB) data from the Planck satellite, and local measurement of the Hubble parameter from the Hubble Space Telescope ($\mathsf H0$). Our analysis shows that the CMB, LRG, Ly$\alpha$, and JLA data are consistent with each other and with a $\Lambda$CDM cosmology, but the ${\mathsf H0}$ data is inconsistent at moderate significance. Including the presence of dark radiation does not alleviate the ${\mathsf H0}$ tension in our analysis. While some of these results have been noted previously, the strength here lies in that we do not assume a particular cosmological model. We calculate the growth of the gravitational potential in General Relativity corresponding to these general expansion histories and show that they are well-approximated by $\Omega_{\rm m}^{0.55}$ given the current precision. We assess the prospects for upcoming surveys to measure deviations from $\Lambda$CDM using this model-independent approach.

The Milky Way's Galactic Center harbors a gamma-ray excess that is a candidate signal of annihilating dark matter. Dwarf galaxies remain predominantly dark in their expected commensurate emission. In this work we quantify the degree of consistency between these two observations through a joint likelihood analysis. In doing so we incorporate Milky Way dark matter halo profile uncertainties, as well as an accounting of diffuse gamma-ray emission uncertainties in dark matter annihilation models for the Galactic Center Extended gamma-ray excess (GCE) detected by the Fermi Gamma-Ray Space Telescope. The preferred range of annihilation rates and masses expands when including these unknowns. Even so, using two recent determinations of the Milky Way halo's local density leave the GCE preferred region of single-channel dark matter annihilation models to be in strong tension with annihilation searches in combined dwarf galaxy analyses. A third, higher Milky Way density determination, alleviates this tension. Our joint likelihood analysis allows us to quantify this inconsistency. We provide a set of tools for testing dark matter annihilation models' consistency within this combined dataset. As an example, we test a representative inverse Compton sourced self-interacting dark matter model, which is consistent with both the GCE and dwarfs.

We incorporate Milky Way dark matter halo profile uncertainties, as well as an accounting of diffuse gamma-ray emission uncertainties in dark matter annihilation models for the Galactic Center Extended gamma-ray excess (GCE) detected by the Fermi Gamma Ray Space Telescope. The range of particle annihilation rate and masses expand when including these unknowns. However, two of the most precise empirical determinations of the Milky Way halo's local density and density profile leave the signal region to be in considerable tension with dark matter annihilation searches from combined dwarf galaxy analyses for single-channel dark matter annihilation models. The GCE and dwarf tension can be alleviated if: one, the halo is very highly concentrated or strongly contracted; two, the dark matter annihilation signal differentiates between dwarfs and the GC; or, three, local stellar density measures are found to be significantly lower, like that from recent stellar counts, increasing the local dark matter density.

Self-interacting dark matter (SIDM) models have been proposed to solve the small-scale issues with the collisionless cold dark matter (CDM) paradigm. We derive equilibrium solutions in these SIDM models for the dark matter halo density profile including the gravitational potential of both baryons and dark matter. Self-interactions drive dark matter to be isothermal and this ties the core sizes and shapes of dark matter halos to the spatial distribution of the stars, a radical departure from previous expectations and from CDM predictions. Compared to predictions of SIDM-only simulations, the core sizes are smaller and the core densities are higher, with the largest effects in baryon-dominated galaxies. As an example, we find a core size around 0.5 kpc for dark matter in the Milky Way, more than an order of magnitude smaller than the core size from SIDM-only simulations, which has important implications for indirect searches of SIDM candidates.