Consider a d-dimensional Riemannian manifold M and a diffusion process X = {Xt : 0 ≤ t ≤ T} on M with infinitesimal generator 1/2ΔM + b, where ΔM is the Laplace–Beltrami operator and b is a vector field. For any two smooth curves φ1, φ2 : [0, T] → M,
where ρ is the Riemannian distance, φ ˙ 1 , φ ˙ 2 {\displaystyle \scriptstyle {\dot {\varphi }}_{1},{\dot {\varphi }}_{2}} denote the first derivatives of φ1, φ2, and L is called the Onsager–Machlup function.
The Onsager–Machlup function is given by345
where || ⋅ ||x is the Riemannian norm in the tangent space Tx(M) at x, div b(x) is the divergence of b at x, and R(x) is the scalar curvature at x.
The following examples give explicit expressions for the Onsager–Machlup function of a continuous stochastic processes.
The Onsager–Machlup function of a Wiener process on the real line R is given by6
Proof: Let X = {Xt : 0 ≤ t ≤ T} be a Wiener process on R and let φ : [0, T] → R be a twice differentiable curve such that φ(0) = X0. Define another process Xφ = {Xtφ : 0 ≤ t ≤ T} by Xtφ = Xt − φ(t) and a measure Pφ by
For every ε > 0, the probability that |Xt − φ(t)| ≤ ε for every t ∈ [0, T] satisfies
By Girsanov's theorem, the distribution of Xφ under Pφ equals the distribution of X under P, hence the latter can be substituted by the former:
By Itō's lemma it holds that
where φ ¨ {\displaystyle \scriptstyle {\ddot {\varphi }}} is the second derivative of φ, and so this term is of order ε on the event where |Xt| ≤ ε for every t ∈ [0, T] and will disappear in the limit ε → 0, hence
The Onsager–Machlup function in the one-dimensional case with constant diffusion coefficient σ is given by7
In the d-dimensional case, with σ equal to the unit matrix, it is given by8
where || ⋅ || is the Euclidean norm and
Generalizations have been obtained by weakening the differentiability condition on the curve φ.9 Rather than taking the maximum distance between the stochastic process and the curve over a time interval, other conditions have been considered such as distances based on completely convex norms10 and Hölder, Besov and Sobolev type norms.11
The Onsager–Machlup function can be used for purposes of reweighting and sampling trajectories,12 as well as for determining the most probable trajectory of a diffusion process.1314
Onsager, L. and Machlup, S. (1953) ↩
Stratonovich, R. (1971) ↩
Takahashi, Y. and Watanabe, S. (1980) ↩
Fujita, T. and Kotani, S. (1982) ↩
Wittich, Olaf ↩
Ikeda, N. and Watanabe, S. (1980), Chapter VI, Section 9 ↩
Dürr, D. and Bach, A. (1978) ↩
Zeitouni, O. (1989) ↩
Shepp, L. and Zeitouni, O. (1993) ↩
Capitaine, M. (1995) ↩
Adib, A.B. (2008). ↩
Dürr, D. and Bach, A. (1978). ↩