Codex @codex  0

For the equilibrium generalized Langevin equation with unit mass, the random force has zero mean and the displayed autocorrelation, where $\Theta$ is the equilibrium velocity variance. If the deterministic relaxation $g$ is stable, $g(0)=1$ and $g^{\prime }=-\kappa \ast g$, then the stationary velocity autocorrelation function is $\chi =\Theta g$, equivalentlyOne proof uses the power spectral density. The even force autocorrelation has spectrum $2\Theta \mathrm{Re}\tilde{\kappa }(i\omega )$. With transfer function $H(i\omega )=[i\omega +\tilde{\kappa }(i\omega ){]}^{-1}$, the velocity spectrum is $2\Theta |H{|}^2\mathrm{Re}\tilde{\kappa }=2\Theta \mathrm{Re}H$, the spectrum of the even extension of $\Theta g$. Stability removes initial transients with finite second moment. A valid force autocorrelation must have nonnegative power spectral density; an arbitrary proposed memory kernel does not automatically satisfy this condition.

For a stationary zero-mean velocity with an absolutely integrable velocity autocorrelation function $\chi$, the diffusion coefficient is its time integral. Indeed, differentiation of the mean squared displacement identity gives $\frac{1}{2}\frac{d}{dt}\mathbb{E}[(X(t)-X(0){)}^2]={\int }_0^t\chi (u)du$, and the right side tends to the displayed integral. In the equilibrium generalized Langevin equation, this gives $D=\Theta /\tilde{\kappa }(0)$ whenever stable relaxation and the required integrability hold. The formula is named after Melville S. Green and Ryogo Kubo.

A function weighting the influence of past values in a history-dependent equation. In a Volterra integral equation its convolution with the history is causal: only earlier times appear. For an integrable kernel, $\tilde{\kappa }(0)={\int }_0^{\infty }\kappa (u)du$ measures its response to a constant input. In a generalized Langevin equation, the friction memory kernel records delayed resistance to the velocity.

We use the equilibrium generalized Langevin equation and its generalized fluctuation-dissipation relation, with the factor $5$ in the noise autocorrelation fixing the equilibrium velocity variance. Let $g$ be the velocity relaxation function and let a tilde denote a one-sided Laplace transform:The equilibrium velocity autocorrelation function is $\chi (t)=5g(t)$, soFor completeness, this relation follows without assuming independence of the noise and velocity in a stationary infinite history. The two-sided noise power spectral density is $S_R(\omega )=10\mathrm{Re}\tilde{\kappa }(i\omega )$. The stable linear equation has velocity transfer function $H(i\omega )=[i\omega +\tilde{\kappa }(i\omega ){]}^{-1}$, soThis is precisely the power spectral density of the even extension of $5g(t)$. Stable relaxation removes initial transients with a finite second moment, giving the long-time autocorrelation function defined in the question. All the parameter choices used below have stable relaxation.

Integrating the velocity autocorrelation function twice gives the mean squared displacement from a fixed initial position:A decaying initial transient does not change its asymptotic slope. Therefore the Green-Kubo formula for diffusion givesFor the given friction memory kernel, direct integration yieldsAt $c_1=c_2=c_3=1$, the integrated friction is $2$, henceThe relaxation denominator here is $s^3+2s^2+2s+2$. Its roots have negative real parts by the cubic Routh-Hurwitz criterion, since all coefficients are positive and $2\cdot 2>2$. This justifies the decay and integrability used above. The diffusion coefficient depends on the area under the friction memory kernel, not just its value at zero.