This is the first of a three-part series of papers, in which we study the preheating phase for multifield models of inflation involving nonminimal couplings. In this paper, we study the single-field attractor behavior that these models exhibit during inflation and quantify its strength and parameter dependence. We further demonstrate that the strong single-field attractor behavior persists after the end of inflation. Preheating in such models therefore generically avoids the "de-phasing" that typically affects multifield models with minimally coupled fields, allowing efficient transfer of energy from the oscillating inflaton condensate(s) to coupled perturbations across large portions of parameter space. We develop a doubly-covariant formalism for studying the preheating phase in such models and identify several features specific to multifield models with nonminimal couplings, including effects that arise from the nontrivial field-space manifold. In papers II and III, we apply this formalism to study how the amplification of adiabatic and isocurvature perturbations varies with parameters, highlighting several distinct regimes depending on the magnitude of the nonminimal couplings $\xi_I$.

This is the second in a series of papers on preheating in inflationary models comprised of multiple scalar fields coupled nonminimally to gravity. In this paper, we work in the rigid-spacetime approximation and consider field trajectories within the single-field attractor, which is a generic feature of these models. We construct the Floquet charts to find regions of parameter space in which particle production is efficient for both the adiabatic and isocurvature modes, and analyze the resonance structure using analytic and semi-analytic techniques. Particle production in the adiabatic direction is characterized by the existence of an asymptotic scaling solution at large values of the nonminimal couplings, $\xi_I \gg 1$, in which the dominant instability band arises in the long-wavelength limit, for comoving wavenumbers $k \rightarrow 0$. However, the large-$\xi_I$ regime is not reached until $\xi_I \geq {\cal O} (100)$. In the intermediate regime, with $\xi_I \sim {\cal O}(1 - 10)$, the resonance structure depends strongly on wavenumber and couplings. The resonance structure for isocurvature perturbations is distinct and more complicated than its adiabatic counterpart. An intermediate regime, for $\xi_I \sim {\cal O} (1 - 10)$, is again evident. For large values of $\xi_I$, the Floquet chart consists of densely spaced, nearly parallel instability bands, suggesting a very efficient preheating behavior. The increased efficiency arises from features of the nontrivial field-space manifold in the Einstein frame, which itself arises from the fields' nonminimal couplings in the Jordan frame, and has no analogue in models with minimal couplings. Quantitatively, the approach to the large-$\xi_I$ asymptotic solution for isocurvature modes is slower than in the case of the adiabatic modes.

This paper concludes our semi-analytic study of preheating in inflationary models comprised of multiple scalar fields coupled nonminimally to gravity. Using the covariant framework of paper I in this series, we extend the rigid-spacetime results of paper II by considering both the expansion of the universe during preheating, as well as the effect of the coupled metric perturbations on particle production. The adiabatic and isocurvature perturbations are governed by different effective masses that scale differently with the nonminimal couplings and evolve differently in time. The effective mass for the adiabatic modes is dominated by contributions from the coupled metric perturbations immediately after inflation. The metric perturbations contribute an oscillating tachyonic term that enhances an early period of significant particle production for the adiabatic modes, which ceases on a time-scale governed by the nonminimal couplings $\xi_I$. The effective mass of the isocurvature perturbations, on the other hand, is dominated by contributions from the fields' potential and from the curvature of the field-space manifold (in the Einstein frame), the balance between which shifts on a time-scale governed by $\xi_I$. As in papers I and II, we identify distinct behavior depending on whether the nonminimal couplings are small ($\xi_I \lesssim {\cal O} (1)$), intermediate ($\xi_I \sim {\cal O} (1 - 10)$), or large ($\xi_I \geq 100$).