We have conducted extensive lattice simulations to study the post-inflation dynamics of multifield models involving nonminimal couplings. We explore the parameter dependence of preheating in these models and describe the various time-scales that control such nonlinear processes as energy transfer, re-scattering, and the approach to radiation-domination and thermalization. In the limit of large nonminimal couplings ($\xi_I \sim 100$), we find that efficient transfer of energy from the inflaton condensate to radiative degrees of freedom, emergence of a radiation-dominated equation of state, and the onset of thermalization each consistently occur within $N_{\rm reh} \lesssim 3$ $e$-folds after the end of inflation, largely independent of the values of the other couplings in the models. The exception is the case of negative ellipticity, in which there is a misalignment between the dominant direction in field-space along which the system evolves and the larger of the nonminimal couplings $\xi_I$. In those cases, the field-space-driven parametric resonance is effectively shut off. More generally, the competition between the scalar fields' potential and the field-space manifold structure can yield interesting phenomena such as two-stage resonances. Despite the explosive particle production, which can lead to a quick depletion of the background energy density, the nonlinear processes do not induce any super-horizon correlations after the end of inflation in these models, which keeps predictions for CMB observables unaffected by the late-time amplification of isocurvature fluctuations. Hence the excellent agreement between primordial observables and recent observations is preserved for this class of models, even when we consider post-inflation dynamics.

We study the post-inflation dynamics of multifield models involving nonminimal couplings using lattice simulations to capture significant nonlinear effects like backreaction and rescattering. We measure the effective equation of state and typical time-scales for the onset of thermalization, which could affect the usual mapping between predictions for primordial perturbation spectra and measurements of anisotropies in the cosmic microwave background radiation. For large values of the nonminimal coupling constants, we find efficient particle production that gives rise to nearly instantaneous preheating. Moreover, the strong single-field attractor behavior that was previously identified persists until the end of preheating, thereby suppressing typical signatures of multifield models. We therefore find that predictions for primordial observables in this class of models retain a close match to the latest observations.

Neutrino-matter interactions play an important role in the post-merger evolution of neutron star-neutron star and black hole-neutron star mergers. Most notably, they determine the properties of the bright optical/infrared transients observable after a merger. Unfortunately, Boltzmann's equations of radiation transport remain too costly to be evolved directly in merger simulations. Simulations rely instead on approximate transport algorithms with unquantified modeling errors. In this paper, we use for the first time a time-dependent general relativistic Monte-Carlo (MC) algorithm to solve Boltzmann's equations and estimate important properties of the neutrino distribution function ~10ms after a neutron star merger. We do not fully couple the MC algorithm to the fluid evolution, but use a short evolution of the merger remnant to critically assess errors in our approximate gray two-moment transport scheme. We demonstrate that the analytical closure used by the moment scheme is highly inaccurate in the polar regions, but performs well elsewhere. While the average energy of polar neutrinos is reasonably well captured by the two-moment scheme, estimates for the neutrino energy become less accurate at lower latitudes. The two-moment formalism also overestimates the density of neutrinos in the polar regions by ~50%, and underestimates the neutrino pair-annihilation rate at the poles by factors of 2-3. Although the latter is significantly more accurate than one might have expected before this study, our results indicate that predictions for the properties of polar outflows and for the creation of a baryon-free region at the poles are likely to be affected by errors in the two-moment scheme, thus limiting our ability to reliably model kilonovae and gamma-ray bursts.