We introduce a new method for evaluating the oscillatory integrals which describe natural interference patterns. As an illustrative example of contemporary interest, we consider astrophysical plasma lensing of coherent sources like pulsars and fast radio bursts in radioastronomy. Plasma lenses are known to occur near the source, in the interstellar medium, as well as in the solar wind and the earth's ionosphere. Such lensing is strongest at long wavelengths hence it is generally important to go beyond geometric optics and into the full wave optics regime. Our computational method is a spinoff of new techniques two of us, and our collaborators, have developed for defining and performing Lorentzian path integrals. Cauchy's theorem allows one to transform a computationally fragile and expensive, highly oscillatory integral into an exactly equivalent sum of absolutely and rapidly convergent integrals which can be evaluated in polynomial time. We require only that it is possible to analytically continue the lensing phase, expressed in the integrated coordinates, into the complex domain. We give a first-principles derivation of the Fresnel-Kirchhoff integral, starting from Feynman's path integral for a massless particle in a refractive medium. We then demonstrate the effectiveness of our method by computing the interference patterns of Thom's caustic catastrophes, both in their "normal forms" and within a variety of more realistic, local lens models, over all wavelengths. Our numerical method, implemented in a freely downloadable code, provides a fast, accurate tool for modeling interference patterns in radioastronomy and other fields of physics.

Inflation is most often described using quantum field theory (QFT) on a fixed, curved spacetime background. Such a description is valid only if the spatial volume of the region considered is so large that its size and shape moduli behave classically. However, if we trace an inflating universe back to early times, the volume of any comoving region of interest -- for example the present Hubble volume -- becomes exponentially small. Hence, quantum fluctuations in the trajectory of the background cannot be neglected at early times. In this paper, we develop a path integral description of a flat, inflating patch (approximated as de Sitter spacetime), treating both the background scale factor and the gravitational wave perturbations quantum mechanically. We find this description fails at small values of the initial scale factor, because \emph{two} background saddle point solutions contribute to the path integral. This leads to a breakdown of QFT in curved spacetime, causing the fluctuations to be unstable and out of control. We show the problem may be alleviated by a careful choice of quantum initial conditions, for the background and the fluctuations, provided that the volume of the initial, inflating patch is larger than $\gg H^{-1}$ in Planck units with $H$ the Hubble constant at the start of inflation. The price of the remedy is high: not only the inflating background, but also the stable, Bunch-Davies fluctuations must be input by hand. Our discussion emphasizes that, even if the inflationary scale is far below the Planck mass, new physics is required to explain the initial quantum state of the universe.

The topology and geometry of random fields - in terms of the Euler characteristic and the Minkowski functionals - has received a lot of attention in the context of the Cosmic Microwave Background (CMB), as the detection of primordial non-Gaussianities would form a valuable clue on the physics of the early Universe. The virtue of both the Euler characteristic and the Minkowski functionals in general, lies in the fact that there exist closed form expressions for their expectation values for Gaussian random fields. However, the Euler characteristic and Minkowski functionals are summarizing characteristics of topology and geometry. Considerably more topological information is contained in the homology of the random field, as it completely describes the creation, merging and disappearance of topological features in superlevel set filtrations. In the present study we extend the topological analysis of the superlevel set filtrations of two-dimensional Gaussian random fields by analysing the statistical properties of the Betti numbers - counting the number of connected components and loops - and the persistence diagrams - describing the creation and mergers of homological features. Using the link between homology and the critical points of a function - as illustrated by the Morse-Smale complex - we derive a one-parameter fitting formula for the expectation value of the Betti numbers and forward this formalism to the persistent diagrams. We, moreover, numerically demonstrate the sensitivity of the Betti numbers and persistence diagrams to the presence of non-Gaussianities.

This study presents a numerical analysis of the topology of a set of cosmologically interesting three-dimensional Gaussian random fields in terms of their Betti numbers $\beta_0$, $\beta_1$ and $\beta_2$. We show that Betti numbers entail a considerably richer characterization of the topology of the primordial density field. Of particular interest is that Betti numbers specify which topological features - islands, cavities or tunnels - define its spatial structure. A principal characteristic of Gaussian fields is that the three Betti numbers dominate the topology at different density ranges. At extreme density levels, the topology is dominated by a single class of features. At low levels this is a \emph{Swiss-cheeselike} topology, dominated by isolated cavities, at high levels a predominantly \emph{Meatball-like} topology of isolated objects. At moderate density levels, two Betti number define a more \emph{Sponge-like} topology. At mean density, the topology even needs three Betti numbers, quantifying a field consisting of several disconnected complexes, not of one connected and percolating overdensity. A {\it second} important aspect of Betti number statistics is that they are sensitive to the power spectrum. It reveals a monotonic trend in which at a moderate density range a lower spectral index corresponds to a considerably higher (relative) population of cavities and islands. We also assess the level of complementary information that Betti numbers represent, in addition to conventional measures such as Minkowski functionals. To this end, we include an extensive description of the Gaussian Kinematic Formula (GKF), which represents a major theoretical underpinning for this discussion.

We present a general formalism for identifying the caustic structure of an evolving mass distribution in an arbitrary dimensional space. For the class of Hamiltonian fluids the identification corresponds to the classification of singularities in Lagrangian catastrophe theory. Based on this we develop a theoretical framework for the formation of the cosmic web, and specifically those aspects that characterize its unique nature: its complex topological connectivity and multiscale spinal structure of sheetlike membranes, elongated filaments and compact cluster nodes. The present work represents an extension of the work by Arnol'd et al., who classified the caustics for the 1- and 2-dimensional Zel'dovich approximation. His seminal work established the role of emerging singularities in the formation of nonlinear structures in the universe. At the transition from the linear to nonlinear structure evolution, the first complex features emerge at locations where different fluid elements cross to establish multistream regions. The classification and characterization of these mass element foldings can be encapsulated in caustic conditions on the eigenvalue and eigenvector fields of the deformation tensor field. We introduce an alternative and transparent proof for Lagrangian catastrophe theory, and derive the caustic conditions for general Lagrangian fluids, with arbitrary dynamics, including dissipative terms and vorticity. The new proof allows us to describe the full 3-dimensional complexity of the gravitationally evolving cosmic matter field. One of our key findings is the significance of the eigenvector field of the deformation field for outlining the spatial structure of the caustic skeleton. We consider the caustic conditions for the 3-dimensional Zel'dovich approximation, extending earlier work on those for 1- and 2-dimensional fluids towards the full spatial richness of the cosmic web.

We present a general formalism for identifying the caustic structure of an evolving mass distribution in an arbitrary dimensional space. For the class of Hamiltonian fluids the identification corresponds to the classification of singularities in Lagrangian catastrophe theory. Based on this we develop a theoretical framework for the formation of the cosmic web, and specifically those aspects that characterize its unique nature: its complex topological connectivity and multiscale spinal structure of sheetlike membranes, elongated filaments and compact cluster nodes. The present work represents an extension of the work by Arnol'd et al., who classified the caustics for the 1- and 2-dimensional Zel'dovich approximation. His seminal work established the role of emerging singularities in the formation of nonlinear structures in the universe. At the transition from the linear to nonlinear structure evolution, the first complex features emerge at locations where different fluid elements cross to establish multistream regions. The classification and characterization of these mass element foldings can be encapsulated in caustic conditions on the eigenvalue and eigenvector fields of the deformation tensor field. We introduce an alternative and transparent proof for Lagrangian catastrophe theory, and derive the caustic conditions for general Lagrangian fluids, with arbitrary dynamics, including dissipative terms and vorticity. The new proof allows us to describe the full 3-dimensional complexity of the gravitationally evolving cosmic matter field. One of our key findings is the significance of the eigenvector field of the deformation field for outlining the spatial structure of the caustic skeleton. We consider the caustic conditions for the 3-dimensional Zel'dovich approximation, extending earlier work on those for 1- and 2-dimensional fluids towards the full spatial richness of the cosmic web.

The cosmic web is a complex spatial pattern of walls, filaments, cluster nodes and underdense void regions. It emerged through gravitational amplification from the Gaussian primordial density field. Here we infer analytical expressions for the spatial statistics of caustics in the evolving large-scale mass distribution. In our analysis, following the quasi-linear Zeldovich formalism and confined to the 1D and 2D situation, we compute number density and correlation properties of caustics in cosmic density fields that evolve from Gaussian primordial conditions. The analysis can be straightforwardly extended to the 3D situation. We moreover, are currently extending the approach to the non-linear regime of structure formation by including higher order Lagrangian approximations and Lagrangian effective field theory.

We study the topology of the Megaparsec Cosmic Web in terms of the scale-dependent Betti numbers, which formalize the topological information content of the cosmic mass distribution. While the Betti numbers do not fully quantify topology, they extend the information beyond conventional cosmological studies of topology in terms of genus and Euler characteristic. The richer information content of Betti numbers goes along the availability of fast algorithms to compute them. For continuous density fields, we determine the scale-dependence of Betti numbers by invoking the cosmologically familiar filtration of sublevel or superlevel sets defined by density thresholds. For the discrete galaxy distribution, however, the analysis is based on the alpha shapes of the particles. These simplicial complexes constitute an ordered sequence of nested subsets of the Delaunay tessellation, a filtration defined by the scale parameter, $\alpha$. As they are homotopy equivalent to the sublevel sets of the distance field, they are an excellent tool for assessing the topological structure of a discrete point distribution. In order to develop an intuitive understanding for the behavior of Betti numbers as a function of $\alpha$, and their relation to the morphological patterns in the Cosmic Web, we first study them within the context of simple heuristic Voronoi clustering models. Subsequently, we address the topology of structures emerging in the standard LCDM scenario and in cosmological scenarios with alternative dark energy content. The evolution and scale-dependence of the Betti numbers is shown to reflect the hierarchical evolution of the Cosmic Web and yields a promising measure of cosmological parameters. We also discuss the expected Betti numbers as a function of the density threshold for superlevel sets of a Gaussian random field.