 Solution

For $0<\omega <ϵ$, let $x=\theta -{\theta }_s$ and $\kappa =\sqrt{{ϵ}^2-{\omega }^2}$. In a single well, sufficiently small fluctuations obey the linear Overdamped Langevin dynamicsThis is an Ornstein-Uhlenbeck process. Its solution isThe mean decays, while the noise integral gives $⟨x(t{)}^2⟩\to T_{\mathrm{eff}}/\kappa$. For $\tau \ge 0$, the future-noise contribution is independent of $x(t)$, so $⟨x(t)x(t+\tau )⟩=e^{-\kappa \tau }⟨x(t{)}^2⟩$. The intrawell phase autocorrelation is thereforeThe connected unnormalized time autocorrelation is just the exponentially decaying second term. The constant first term is necessary because the requested raw correlation is not centered at the equilibrium point. This is the harmonic-well approximation. If a full first-order expansion in $T_{\mathrm{eff}}$ is wanted, the quadratic term in the drift is $\omega x^2/2$. Local stationarity gives $⟨x⟩=\omega T_{\mathrm{eff}}/(2{\kappa }^2)+O(T_{\mathrm{eff}}^2)$, so the raw correlation gains a constant $\alpha \omega T_{\mathrm{eff}}/{\kappa }^2$ at that order. The connected correlation is unchanged to first order; the boxed formula is the conventional linearized result rather than a complete nonlinear first-order raw-correlation expansion.

The required time window is long compared with ${\kappa }^{-1}$ but short compared with the escape time, with the lag also short compared with escape. Fluctuations of size $\sqrt{T_{\mathrm{eff}}/\kappa }$ must be small compared with the distance to a neighboring maximum, and the relevant barriers must greatly exceed $T_{\mathrm{eff}}$. For a fixed nonzero noise strength, the unwrapped phase eventually makes phase slips; its raw correlation does not have this stationary infinite-time limit. Thus “large time” here refers to local relaxation within the occupied well, not to the limit after arbitrarily many barrier crossings.

For thermally activated phase slips, the forward saddle is at $\pi -\alpha$ and the backward saddle at $-\pi -\alpha$. Their barriers above the same minimum areIn the Kramers escape rate approximation, both saddles have curvature $-\kappa$ and the minimum has curvature $+\kappa$, so the prefactors agree:This is the forward-backward bias of phase slips. It favors forward motion because each forward step lowers the tilted potential by $2\pi \omega$. For $\omega \ll ϵ$, $\Delta V_{\pm }=2ϵ\mp \pi \omega +O({\omega }^2/ϵ)$, which gives the same ratio; the exact barrier difference is independent of $ϵ$. These probabilities may be read as rates in a common short observation interval, or as normalized competing exit probabilities. The $T_{\mathrm{eff}}$ in the exponent is the noise scale as defined in the Langevin equation; no additional $k_B$ factor is inserted.