Dynamics for spherical spin glasses: Gibbs distributed initial conditions

We derive the coupled non-linear integro-differential equations for the thermodynamic limit of the empirical correlation and response functions in the Langevin dynamics at temperature $T$, for spherical mixed $p$-spin disordered mean-field models, initialized according to a Gibbs measure for temperature $T_0$, in the replica-symmetric (RS) or $1$-replica-symmetry-breaking (RSB) phase. For any $T_0=T$ above the dynamical phase transition point $T_c^{\rm dyn}$ the resulting stationary relaxation dynamics coincide with the FDT solution for these equations, while for lower $T_0=T$ in the $1$-RSB phase, the relaxation dynamics coincides with the FDT solution, now concentrated on the single spherical band within the Gibbs measure's support on which the initial point lies.