We analyze the onset of inflation for a simple single-field model when the system begins with significant inhomogeneities on length-scales shorter than the initial Hubble radius. We incorporate certain nonlinear interactions among the coupled degrees of freedom by using the nonperturbative Hartree approximation. Consistent with recent, more computationally intensive numerical-relativity studies, we find inflation to be robust for large-field models, even when the system begins with significant structure on sub-Hubble scales. We consider the space of initial conditions $(\varphi (t_0), \dot{\varphi} (t_0))$, where $\varphi$ is the vacuum expectation value of the quantized field that drives inflation. Although some regions of $(\varphi (t_0), \dot{\varphi} (t_0))$ that would have yielded sufficient inflation in the absence of inhomogeneities fail to do so when backreaction from inhomogeneities is incorporated, an equal volume of such regions succeeds in producing sufficient inflation which did not do so in the absence of inhomogeneities. For large-field models, in other words, the total volume of the space of initial conditions $(\varphi (t_0), \dot{\varphi} (t_0))$ that yields sufficient inflation is {\it conserved} when we incorporate nonlinear backreaction from inhomogeneities, compared to the case in which inhomogeneities are neglected.

Chi-squared random fields arise naturally from the study of fluctuations in field theories with SO(n) symmetry. The extrema of chi-squared fields are of particular physical interest. In this paper, we undertake a statistical analysis of the stationary points of chi-squared fields, with particular emphasis on extrema. We begin by describing the neighborhood of a stationary point in terms of a biased chi-squared random field, and then compute the expected profile of this field, as well as a variety of associated statistics. We are interested in understanding how spherically symmetric the neighborhood of a stationary point is, on average. To this end, we decompose the biased field into its spherical harmonic modes about this point, and explore their statistics. Using these mode statistics, we construct a metric to gauge the degree of spherical symmetry of the field in this neighborhood. Finally, we show how to leverage the harmonic decomposition to efficiently sample both Gaussian and chi-squared fields about a stationary point.

We investigate the statistics of stationary points in the sum of squares of $N$ Gaussian random fields, which we call a "chi-squared" field. The behavior of such a field at a point is investigated, with particular attention paid to the formation of topological defects. An integral to compute the number density of stationary points at a given field amplitude is constructed. We compute exact expressions for the integral in various limits and provide code to evaluate it numerically in the general case. We investigate the dependence of the number density of stationary points on the field amplitude, number of fields, and power spectrum of the individual Gaussian random fields. This work parallels the work of Bardeen, Bond, Kaiser and Szalay, who investigated the statistics of peaks of Gaussian random fields. A number of results for integrating over matrices are presented in appendices.

Scalar field theories that possess a Vainshtein mechanism are able to dynamically suppress the associated fifth forces in the presence of massive sources through derivative non-linearities. The resulting equations of motion for the scalar are highly non-linear and therefore very few analytic solutions are known. Here we present a brief investigation of the structure of Vainshtein screening in symmetrical configurations, focusing in particular on the spherical, cylindrical and planar solutions that are relevant for observations of the cosmic web. We consider Vainshtein screening in both the Galileon model, where the non-linear terms involve second derivatives of the scalar, and a k-essence theory, where the non-linear terms involve only first derivatives of the scalar. We find that screening, and consequently the suppression of the scalar force, is most efficient around spherical sources, weaker around cylindrical sources and can be absent altogether around planar sources.

We present a comprehensive formalism for the description of primordial black hole formation in spherical symmetry based on the formalisms of Misner, Sharp, and Hernandez, which can be used to predict whether or not a black hole will form, and extract the resulting black hole mass when formation does occur. Rigorous derivations of all aspects of the formalism are provided, including a thorough investigation of appropriate initial and boundary conditions. We connect our formalism with numerous other approaches in the literature. Some implementation details for numerical code are provided. We include animations of simulated primordial black hole formation as supplemental material.

The most general covariant action describing gravity coupled to a scalar field with only second order equations of motion, Horndeski's theory (also known as "Generalized Galileons"), provides an all-encompassing model in which single scalar dark energy models may be constrained. However, the generality of the model makes it cumbersome to manipulate. In this paper, we demonstrate that when considering linear perturbations about a Friedmann-Robertson-Walker background, the theory is completely specified by only six functions of time, two of which are constrained by the background evolution. We utilise the ideas of the Effective Field Theory of Inflation/Dark Energy to explicitly construct these six functions of time in terms of the free functions appearing in Horndeski's theory. These results are used to investigate the behavior of the theory in the quasistatic approximation. We find that only four functions of time are required to completely specify the linear behavior of the theory in this limit, which can further be reduced if the background evolution is fixed. This presents a significantly reduced parameter space from the original presentation of Horndeski's theory, giving hope to the possibility of constraining the parameter space. This work provides a cross-check for previous work on linear perturbations in this theory, and also generalizes it to include spatial curvature.

Cosmic superstring loops within the galaxy microlens background point sources lying close to the observer-string line of sight. For suitable alignments, multiple paths coexist and the (achromatic) flux enhancement is a factor of two. We explore this unique type of lensing by numerically solving for geodesics that extend from source to observer as they pass near an oscillating string. We characterize the duration of the flux doubling and the scale of the image splitting. We probe and confirm the existence of a variety of fundamental effects predicted from previous analyses of the static infinite straight string: the deficit angle, the Kaiser-Stebbins effect, and the scale of the impact parameter required to produce microlensing. Our quantitative results for dynamical loops vary by O(1) factors with respect to estimates based on infinite straight strings for a given impact parameter. A number of new features are identified in the computed microlensing solutions. Our results suggest that optical microlensing can offer a new and potentially powerful methodology for searches for superstring loop relics of the inflationary era.

We present a formalism for the numerical implementation of general theories of dark energy, combining the computational simplicity of the equation of state for perturbations approach with the generality of the effective field theory approach. An effective fluid description is employed, based on a general action describing single-scalar field models. The formalism is developed from first principles, and constructed keeping the goal of a simple implementation into CAMB in mind. Benefits of this approach include its straightforward implementation, the generality of the underlying theory, the fact that the evolved variables are physical quantities, and that model-independent phenomenological descriptions may be straightforwardly investigated. We hope this formulation will provide a powerful tool for the comparison of theoretical models of dark energy with observational data.

We take an Effective Field Theory (EFT) approach to unifying existing proposals for the origin of cosmic acceleration and its connection to cosmological observations. Building on earlier work where EFT methods were used with observations to constrain the background evolution, we extend this program to the level of the EFT of the cosmological perturbations - following the example from the EFT of Inflation. Within this framework, we construct the general theory around an assumed background which will typically be chosen to mimic Lambda-CDM, and identify the parameters of interest for constraining dark energy and modified gravity models with observations. We discuss the similarities to the EFT of Inflation, but we also identify a number of subtleties including the relationship between the scalar perturbations and the Goldstone boson of the spontaneously broken time translations. We present formulae that relate the parameters of the fundamental Lagrangian to the speed of sound, anisotropic shear stress, effective Newtonian constant, and Caldwell's varpi parameter emphasizing the connection to observations. It is anticipated that this framework will be of use in constraining individual models, as well as for placing model-independent constraints on dark energy and modified gravity model building.

We explore a class of effective field theory models of cosmic acceleration involving a metric and a single scalar field. These models can be obtained by starting with a set of ultralight pseudo-Nambu-Goldstone bosons whose couplings to matter satisfy the weak equivalence principle, assuming that one boson is lighter than all the others, and integrating out the heavier fields. The result is a quintessence model with matter coupling, together with a series of correction terms in the action in a covariant derivative expansion, with specific scalings for the coefficients. After eliminating higher derivative terms and exploiting the field redefinition freedom, we show that the resulting theory contains nine independent free functions of the scalar field when truncated at four derivatives. This is in contrast to the four free functions found in similar theories of single-field inflation, where matter is not present. We discuss several different representations of the theory that can be obtained using the field redefinition freedom. For perturbations to the quintessence field today on subhorizon lengthscales larger than the Compton wavelength of the heavy fields, the theory is weakly coupled and natural in the sense of t'Hooft. The theory admits a regime where the perturbations become modestly nonlinear, but very strong nonlinearities lie outside its domain of validity.