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.