On the generalized Langevin equation and the Mori projection operator technique

Abstract In statistical physics, the Nakajima-Mori-Zwanzig projection operator formalism is used to derive an integro-differential equation for observables in a Hilbert space, the generalized Langevin equation (GLE). This technique relies on the splitting of the dynamics into a projected and an orthogonal part. However, the well-posedness of the abstract Cauchy problem for the orthogonal dynamics remains an open problem. Moreover, it is rarely discussed under which assumptions the Dyson identity, which is used to derive the GLE, holds. In this article, we address this issue for rank-one projections (Mori’s projection). For the Mori projection, the orthogonal dynamics is obtained from the bounded perturbation theorem. The variation of constants formula for strongly continuous semigroups then yields the GLE and the second fluctuation dissipation theorem (2FDT). We show that the variation of constants can be replaced by a limiting process in order to give a general proof of the GLE and 2FDT that does not require the differentiability of the fluctuating forces. In addition, we offer an alternative approach that does not require the bounded perturbation theorem. Our starting point is the observation that the GLE and 2FDT uniquely determine the fluctuating forces as well as the memory kernel. Furthermore, the orbit maps for the orthogonal dynamics can be directly defined via solutions of linear Volterra equations. All desired properties of the orthogonal dynamics are then proven directly from this definition. In particular, the orthogonal dynamics is a strongly continuous semigroup generated by <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:mover> <mml:mrow> <mml:mi class="MJX-tex-calligraphic">QL</mml:mi> </mml:mrow> <mml:mo accent="true">―</mml:mo> </mml:mover> <mml:mrow> <mml:mi class="MJX-tex-calligraphic">Q</mml:mi> </mml:mrow> <mml:mo>=</mml:mo> <mml:mrow> <mml:mi class="MJX-tex-calligraphic">QLQ</mml:mi> </mml:mrow> </mml:mrow> </mml:math> , where <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:mrow> <mml:mi class="MJX-tex-calligraphic">L</mml:mi> </mml:mrow> </mml:mrow> </mml:math> is the generator of the time evolution operator, and <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:mrow> <mml:mi class="MJX-tex-calligraphic">P</mml:mi> </mml:mrow> <mml:mo>=</mml:mo> <mml:mn>1</mml:mn> <mml:mo>−</mml:mo> <mml:mrow> <mml:mi class="MJX-tex-calligraphic">Q</mml:mi> </mml:mrow> </mml:mrow> </mml:math> is the Mori projection operator. If <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" overflow="scroll"> <mml:mrow> <mml:mrow> <mml:mi class="MJX-tex-calligraphic">L</mml:mi> </mml:mrow> </mml:mrow> </mml:math> is skew-adjoint, the orthogonal dynamics is a unitary group and the fluctuating forces are stationary. Our results apply to general autonomous dynamical systems whose time evolution is given by a strongly continuous semigroup.