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.

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.