A number of codes for general-relativistic simulations of cosmological structure formation have been developed in recent years. Here we demonstrate that a sample of these codes produce consistent results beyond the Newtonian regime. We simulate solutions to Einstein's equations dominated by gravitomagnetism -- a vector-type gravitational field that doesn't exist in Newtonian gravity and produces frame-dragging, the leading-order post-Newtonian effect. We calculate the coordinate-invariant effect on intersecting null geodesics by performing ray tracing in each independent code. With this observable quantity, we assess and compare each code's ability to compute relativistic effects.

We examine how cosmological expansion arises in a universe containing a lattice of spinning black holes. We study averaged expansion properties as a function of fundamental properties of the black holes, including the bare mass of the black holes and black hole spin. We then explore how closely the expansion properties correspond to properties of a corresponding matter-dominated FLRW universe. As residual radiation present in the initial data decays, we find good agreement with a matter-dominated FLRW solution, and the effective density in the volume is well-described by the horizon mass of the black hole.

Compensated isocurvature perturbations (CIPs) are modulations of the relative baryon and dark matter density that leave the total matter density constant. The best current constraints from the primary cosmic microwave background (CMB) are consistent with CIPs some two orders of magnitude larger in amplitude than adiabatic perturbations, suggesting that there may be a huge gap in our knowledge of the early Universe. However, it was recently suggested by Barreira~et.~al. that CIPs which are correlated with the primordial curvature perturbation, as arises in some versions of the curvaton model, lead to a new observable: scale dependent galaxy bias. Combining a galaxy survey with an unbiased tracer of the density field facilitates a measurement of the amplitude of correlated CIPs that is free from cosmic variance, the main limitation on constraints from the primary CMB. Among the most promising tracers to use for this purpose is the remote dipole field, reconstructed using the technique of kinetic Sunyaev Zel'dovich (kSZ) tomography. In this paper, we evaluate the detection significance on the amplitude of correlated CIPs possible with next-generation CMB and galaxy surveys using kSZ tomography. Our analysis includes all relativistic contributions to the observed galaxy number counts and allows for both CIPs and primordial non-Gaussianity, which also gives rise to a scale dependent galaxy bias. We find that kSZ tomography can probe CIPs of comparable amplitude to the adiabatic fluctuations, representing an improvement of over two orders of magnitude upon current constraints, and an order of magnitude over what will be possible using future CMB or galaxy surveys alone.

The Simons Observatory (SO) is a ground-based cosmic microwave background (CMB) experiment sited on Cerro Toco in the Atacama Desert in Chile that promises to provide breakthrough discoveries in fundamental physics, cosmology, and astrophysics. Supported by the Simons Foundation, the Heising-Simons Foundation, and with contributions from collaborating institutions, SO will see first light in 2021 and start a five year survey in 2022. SO has 287 collaborators from 12 countries and 53 institutions, including 85 students and 90 postdocs. The SO experiment in its currently funded form ('SO-Nominal') consists of three 0.4 m Small Aperture Telescopes (SATs) and one 6 m Large Aperture Telescope (LAT). Optimized for minimizing systematic errors in polarization measurements at large angular scales, the SATs will perform a deep, degree-scale survey of 10% of the sky to search for the signature of primordial gravitational waves. The LAT will survey 40% of the sky with arc-minute resolution. These observations will measure (or limit) the sum of neutrino masses, search for light relics, measure the early behavior of Dark Energy, and refine our understanding of the intergalactic medium, clusters and the role of feedback in galaxy formation. With up to ten times the sensitivity and five times the angular resolution of the Planck satellite, and roughly an order of magnitude increase in mapping speed over currently operating ("Stage 3") experiments, SO will measure the CMB temperature and polarization fluctuations to exquisite precision in six frequency bands from 27 to 280 GHz. SO will rapidly advance CMB science while informing the design of future observatories such as CMB-S4.

High-resolution, low-noise observations of the cosmic microwave background (CMB) planned for the near-future will enable new cosmological probes based on re-scattered CMB photons -- the secondary CMB. At the same time, enormous galaxy surveys will map out huge volumes of the observable Universe. Using the technique of kinetic Sunyaev Zel'dovich (kSZ) tomography these new probes can be combined to reconstruct the remote dipole field, the CMB dipole as observed from different vantage points in our Universe. The volume accessible to future galaxy surveys is large enough that general relativistic corrections to the observed distribution of galaxies must be taken into account. These corrections are interesting probes of gravity in their own right, but can also obscure potential signatures of primordial non-Gaussianity. In this paper, we demonstrate that correlations between the reconstructed remote dipole field and the observed galaxy number counts can in principle be used to detect general relativistic corrections. We show that neglecting general relativistic corrections leads to an $\mathcal{O}(1)$ bias on the inferred amplitude of primordial non-Gaussianity, $f_{\rm NL}$. In addition, we demonstrate that the reconstructed remote dipole field can provide useful constraining power on various bias parameters appearing in the galaxy number counts, and can significantly mitigate the effects of alignment bias.

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.

Gravitational potentials which change in time induce fluctuations in the observed cosmic microwave background (CMB) temperature. Cosmological structure moving transverse to our line of sight provides a specific example known as the moving lens effect. Here we explore how the observed CMB temperature fluctuations combined with the observed matter over-density can be used to infer the transverse velocity of cosmological structure on large scales. We show that near-future CMB surveys and galaxy surveys will have the statistical power to make a first detection of the moving lens effect, and we discuss applications for the reconstructed transverse velocity.

Standard cosmological models rely on an approximate treatment of gravity, utilizing solutions of the linearized Einstein equations as well as physical approximations. In an era of precision cosmology, we should ask: are these approximate predictions sufficiently accurate for comparison to observations, and can we draw meaningful conclusions about properties of our Universe from them? In this work we examine the accuracy of linearized gravity in the presence of collisionless matter and a cosmological constant utilizing fully general relativistic simulations. We observe the gauge-dependence of corrections to linear theory, and note the amplitude of these corrections. For perturbations whose amplitudes are in line with expectations from the standard $\Lambda$CDM model, we find that the full, general relativistic metric is well-described by linear theory in Newtonian and harmonic gauges, while the metric in comoving-synchronous gauge is not. For the largest observed structures in our Universe, our results suggest that corrections to linear gravitational theory can reach or surpass the percent-level.

The kinetic Sunyaev Zel'dovich (kSZ) effect, cosmic microwave background (CMB) anisotropies induced by scattering from free electrons in bulk motion, is a primary target of future CMB experiments. Measurements of the kSZ effect have the potential to address fundamental questions about the structure and evolution of our Universe on the largest scales and at the earliest times. This potential is unlocked by combining measurements of small-scale CMB anisotropies with large-scale structure surveys, a technique known as kSZ tomography. Previous work established a quadratic estimator for the remote dipole field, the CMB dipole observed at different locations in the Universe. This previous work did not include gravitational lensing, redshift space distortions, or non-linear evolution of structure. In this paper, we investigate how well the remote dipole field can be reconstructed in the presence of such effects by using mock data from a suite of simulations. To properly model both large and small scales, we develop a novel box-in-box simulation pipeline, where small-scale information is obtained from N-body simulations, and large-scale information obtained by evolving fields using linear theory and adding the resulting corrections to the N-body particle data. This pipeline allows us to create properly correlated maps of the primary CMB including lensing, as well as the kSZ effect and density maps on the past light cone of an observer. Analyzing an ensemble of mocks, we find that the dipole field can be reconstructed with high fidelity over a range of angular scales and redshift bins, although there is evidence of excess power from nonlinear structure. We also analyze correlations with the primary CMB, investigating the ability of kSZ tomography to reconstruct the intrinsic CMB dipole. Our results constitute a proof-of-principle that kSZ tomography is a promising technique for future datasets.

We compute the weak lensing convergence power spectrum, $C^{\kappa\kappa}(\theta)$, in a dust-filled universe using fully non-linear general relativistic simulations. The spectrum is then compared to more standard, approximate calculations by computing the Bardeen (Newtonian) potentials in linearized gravity and utilizing the Born approximation. We find corrections to the angular power spectrum amplitude of order ten percent at very large angular scales, $\ell ~ 2-3$, and percent-level corrections at intermediate angular scales of $\ell ~ 20-30$.

How does inhomogeneity affect our interpretation of cosmological observations? It has long been wondered to what extent the observable properties of an inhomogeneous universe differ from those of a corresponding Friedman-Lemaitre-Robertson-Walker (FLRW) model, and how the inhomogeneities affect that correspondence. Here, we use numerical relativity to study the behavior of light beams traversing an inhomogeneous universe and construct the resulting Hubble diagrams. The universe that emerges exhibits an average FLRW behavior, but inhomogeneous structures contribute to deviations in observables across the observer's sky. We also investigate the relationship between angular diameter distance and the angular extent of a source, finding deviations that grow with source redshift. These departures from FLRW are important path-dependent effects with implications for using real observables in an inhomogeneous universe such as our own.

We present cosmological-scale numerical simulations of an evolving universe in full general relativity (GR) and introduce a new numerical tool, {\sc CosmoGRaPH}, which employs the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formalism on a 3-dimensional grid. Using {\sc CosmoGRaPH}, we calculate the effect of an inhomogeneous matter distribution on the evolution of a spacetime. We also present the results of a set of standard stability tests to demonstrate the robustness of our simulations.

While the use of numerical general relativity for modeling astrophysical phenomena and compact objects is commonplace, the application to cosmological scenarios is only just beginning. Here, we examine the expansion of a spacetime using the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formalism of numerical relativity in synchronous gauge. This work represents the first numerical cosmological study that is fully relativistic, non-linear and without symmetry. The universe that emerges exhibits an average Friedmann-Lema\"itre-Robertston-Walker (FLRW) behavior, however this universe also exhibits locally inhomogeneous expansion beyond that expected in linear perturbation theory around a FLRW background.

First-order phase transitions are a source of stochastic gravitational radiation. Precision calculations of the gravitational waves emitted during these processes, sourced by both the degrees of freedom undergoing the transition and anisotropic stress of the coupled, ambient constituents, have reached an age of maturity. Here we present high-resolution numerical simulations of a scalar field coupled to a fluid and parameterize the final gravitational wave spectrum as a function of the ratio of the energies of the two sectors and the coupling between the two sectors for a set of models that represent different types of first-order phase transitions. In most cases, the field sector is the dominant source of gravitational radiation, but it is possible in certain scenarios for the fluid to have the most important contribution.

In this work we investigate the effect on weak-lensing shear and convergence measurements due to distortions from the Lorentz boost induced by our Galaxy's motion. While no ellipticity is induced in an image from the Lorentz boost to first order in beta = v/c, the image is magnified. This affects the inferred convergence at a 10 per cent level, and is most notable for low multipoles in the convergence power spectrum C {\kappa}{\kappa} and for surveys with large sky coverage like LSST and DES. Experiments which image only small fractions of the sky and convergence power spectrum determinations at l > 5 can safely neglect the boost effect to first order in beta.