We prove joint functional convergence of the stationary and quench autocorrelations of zero-field Sherrington–Kirkpatrick heat-bath dynamics at inverse temperature β = 1, with the same random limit for Gaussian and Rademacher couplings. Each site has a rate-one clock, mean spin autocorrelations are multiplied by n^1/3, and waiting times and lags are measured in units n^2/3. The quench starts from independent fair spins, and convergence is uniform on compact sets of positive waiting times and lags. The quench limit is selected by these initial states and relaxes to the stationary limiting autocorrelation as the waiting time tends to infinity.

We prove that the stationary spin autocorrelation of the zero-field Sherrington–Kirkpatrick model at inverse temperature β = 1 has a common functional scaling limit for Gaussian and Rademacher couplings. With rate-one heat-bath clocks at every site, time is scaled by n^2/3 and the mean spin autocorrelation is multiplied by n^1/3. The functions converge in law uniformly on compact positive-time intervals. Their limiting law is not a point mass: it retains sample-to-sample randomness, and its functions decay to zero at large times.

For every fixed inverse temperature $0\lt \beta\lt 1$, we prove that the unscaled spectral gap of single-site heat-bath dynamics for the zero-field Gaussian Sherrington–Kirkpatrick model is bounded away from zero with probability tending to one over the disorder. Equivalently, the Gibbs law satisfies a dimension-free Poincaré inequality for all functions.

We prove worst-case total-variation cutoff for the zero-field Gaussian Sherrington–Kirkpatrick heat-bath dynamics at every fixed inverse temperature $0\leq\beta\lt 1$. With rate-one refresh at each spin, the cutoff location is $\log n/(2\lambda(\beta))$ for a positive deterministic rate $\lambda(\beta)$. The location for uniformly chosen single-site update attempts is n times as large. Convergence is in probability over the disorder.

At the critical inverse temperature β = 1, we prove slow mixing for zero-field Gaussian Sherrington–Kirkpatrick heat-bath dynamics from typical equilibrium configurations held fixed as initial states. For every deterministic sequence $t_n=o(n^{2/3})$ in rate-one-per-site time, the Gibbs mass of initial states whose time-t_n total-variation distance from equilibrium exceeds 1/4 tends to one in probability over the disorder. The same statement holds for every deterministic integer sequence $k_n=o(n^{5/3})$ of uniform-site update attempts.

At every fixed inverse temperature β > 1, we prove a stretched-exponential obstruction to mixing for the zero-field Sherrington–Kirkpatrick model from typical equilibrium configurations held fixed as initial states. The Gibbs mass of states whose total-variation distance from equilibrium at time $\exp(n^{1/10000})$ exceeds 1/4 tends to one in probability over the disorder. This holds for rate-one-per-site heat-bath dynamics and after the same stated number of discrete update attempts.

For every fixed inverse temperature β > 1, we prove subexponential mixing for rate-one-per-site heat-bath dynamics in the zero-field Gaussian Sherrington–Kirkpatrick model from a typical equilibrium configuration. If one Gibbs-sampled configuration is held fixed as the initial state, its total-variation distance from equilibrium is at most 1/4 by time $\exp(n^{1-1/40000000})$, with joint probability tending to one over the disorder and the sampled state.

At the critical inverse temperature β = 1, we determine the worst-start total-variation mixing exponent for single-site heat-bath dynamics in the zero-field Gaussian Sherrington–Kirkpatrick model. The mixing time is $n^{2/3+o(1)}$ when every site has rate one, or $n^{5/3+o(1)}$ attempted uniform-site updates, in probability over the disorder.