The apparent properties of distant objects encode information about the way the light they emit propagates to an observer, and therefore about the curvature of the underlying spacetime. Measuring the relationship between the redshift $z$ and the luminosity distance $D_{\rm L}$ of a standard candle, for example, yields information on the Universe's matter content. In practice, however, in order to decode this information the observer needs to make an assumption about the functional form of the $D_{\rm L}(z)$ relation; in other words, a cosmological model needs to be assumed. In this work, we use numerical-relativity simulations, equipped with a new ray-tracing module, to numerically obtain this relation for a few black-hole--lattice cosmologies and compare it to the well-known Friedmann-Lema\^itre-Robertson-Walker case, as well as to other relevant cosmologies and to the Empty-Beam Approximation. We find that the latter provides the best estimate of the luminosity distance and formulate a simple argument to account for this agreement. We also find that a Friedmann-Lema\^itre-Robertson-Walker model can reproduce this observable exactly, as long as a time-dependent cosmological constant is included in the fit. Finally, the dependence of these results on the lattice mass-to-spacing ratio $\mu$ is discussed: we discover that, unlike the expansion rate, the $D_{\rm L}(z)$ relation in a black-hole lattice does not tend to that measured in the corresponding continuum spacetime as $\mu \to 0$.

The applications of numerical relativity to cosmology are on the rise, contributing insight into such cosmological problems as structure formation, primordial phase transitions, gravitational-wave generation, and inflation. In this paper, I present the infrastructure for the computation of inhomogeneous dust cosmologies which was used recently to measure the effect of nonlinear inhomogeneity on the cosmic expansion rate. I illustrate the code's architecture, provide evidence for its correctness in a number of familiar cosmological settings, and evaluate its parallel performance for grids of up to several billion points. The code, which is available as free software, is based on the Einstein Toolkit infrastructure, and in particular leverages the automated-code-generation capabilities provided by its component Kranc.

We construct a three-dimensional, fully relativistic numerical model of a universe filled with an inhomogeneous pressureless fluid, starting from initial data that represent a perturbation of the Einstein-de Sitter model. We then measure the departure of the average expansion rate with respect to this homogeneous and isotropic reference model, comparing local quantities to the predictions of linear perturbation theory. We find that collapsing perturbations reach the turnaround point much earlier than expected from the reference spherical top-hat collapse model and that the local deviation of the expansion rate from the homogeneous one can be as high as $28\%$ at an underdensity, for an initial density contrast of $10^{-2}$. We then study, for the first time, the exact behavior of the backreaction term ${\cal Q}_{\cal D}$. We find that, for small values of the initial perturbations, this term exhibits a $1/a$ scaling, and that it is negative with a linearly growing absolute value for larger perturbation amplitudes, thereby contributing to an overall deceleration of the expansion. Its magnitude, on the other hand, remains very small even for relatively large perturbations.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

We describe the Einstein Toolkit, a community-driven, freely accessible computational infrastructure intended for use in numerical relativity, relativistic astrophysics, and other applications. The Toolkit, developed by a collaboration involving researchers from multiple institutions around the world, combines a core set of components needed to simulate astrophysical objects such as black holes, compact objects, and collapsing stars, as well as a full suite of analysis tools. The Einstein Toolkit is currently based on the Cactus Framework for high-performance computing and the Carpet adaptive mesh refinement driver. It implements spacetime evolution via the BSSN evolution system and general-relativistic hydrodynamics in a finite-volume discretization. The toolkit is under continuous development and contains many new code components that have been publicly released for the first time and are described in this article. We discuss the motivation behind the release of the toolkit, the philosophy underlying its development, and the goals of the project. A summary of the implemented numerical techniques is included, as are results of numerical test covering a variety of sample astrophysical problems.

We use high-redshift type Ia supernova and compact radio source data in order to test the infrared (IR) fixed point model of the late Universe which was proposed recently. It describes a cosmology with a time dependent cosmological constant and Newton constant whose dynamics arises from an underlying renormalization group flow near an IR-attractive fixed point. Without any finetuning or quintessence field it yields $\Omega_{\rm M}=\Omega_{\Lambda}=1/2$. Its characteristic $t^{4/3}$-dependence of the scale factor leads to a distance-redshift relation whose predictions are compared both to the supernova and to the radio source data. According to the $\chi^2$ test, the fixed point model reproduces the data at least as well as the best-fit (Friedmann-Robertson-Walker) standard cosmology.Furthermore, we extend the original fixed point model by assuming that the fixed point epoch is preceded by an era with constant $G$ and $\Lambda$. By means of a Monte Carlo simulation we show that the data expected from the forthcoming SNAP satellite mission could detect the transition to the fixed point regime provided it took place at a redshift of less than about 0.5.