---
title: Onsager–Machlup Principle for Active Systems
url: https://www.emergentmind.com/topics/onsager-machlup-principle-for-active-systems
type: topic
---

# Onsager–Machlup Principle for Active Systems

Searching arXiv for the core papers and closely related work on Onsager–Machlup formulations for active and nonreciprocal systems.
Search result 1: "Covariant Onsager and Onsager-Machlup principles for active and inertial dynamics" [2604.22371]  
Search result 2: "Onsager's variational principle for nonreciprocal systems with odd elasticity" [2209.15363]  
Search result 3: "Neural optimization of the most probable paths of 3D active Brownian particles" [2511.16178]  
Search result 4: "Statistical formulation of Onsager-Machlup variational principle" [2404.18611]  
Search result 5: "Non-equilibrium Onsager-Machlup theory" [2311.00120]
The Onsager–Machlup principle for active systems denotes a path-variational formulation in which actively driven mesoscopic trajectories are weighted by an action built from an **active Rayleighian**, rather than from purely passive relaxational dynamics. In its most direct current formulation, activity enters through the mechanical power term \(\dot W=X_i(\mathbf q)\dot q^i\) in
\[
R(\mathbf q,\dot{\mathbf q})=\Phi(\mathbf q,\dot{\mathbf q})+\dot F(\mathbf q,\dot{\mathbf q})-\dot W(\mathbf q,\dot{\mathbf q}),
\]
with deterministic drift obtained by minimizing \(R\) and stochastic trajectory weights assigned through the excess Rayleighian \(L=R-R_{\min}\) and the path probability
\[
\rho[\mathbf q(t)\mid \mathbf q_0]=N\exp\!\bigl(-O[\mathbf q(t)]\bigr),\qquad
O[\mathbf q(t)]=\frac{1}{2k_{\rm B}T}\int_0^\tau dt\,L(\mathbf q(t),\dot{\mathbf q}(t)).
\]
This construction extends classical Onsager–Machlup theory to systems with nonconservative active driving and, in the covariant formulation, keeps the theory geometrically consistent under smooth coordinate transformations [2604.22371].

## 1. Formal setting and departure from passive Onsager theory

Classical Onsager theory treats irreversible dynamics as the balance between a free-energy gradient and linear dissipation. That structure is insufficient for active matter because active systems are driven out of equilibrium by forces or propulsion mechanisms that cannot be reduced to a free-energy gradient. In the active formulation, the mesoscopic state is described by generalized coordinates
\[
\mathbf q=(q^1,q^2,\dots),
\]
equipped with a metric \(g_{ij}(\mathbf q)\), and thermodynamic quantities are required to be scalars on the reduced state space. A key definition is the mesoscopic entropy
\[
S_{\mathrm{mes}}(\mathbf q)=-k_{\rm B}\ln \rho(\mathbf q)
=-k_{\rm B}\ln\!\left[\frac{P(\mathbf q)}{\sqrt{g(\mathbf q)}}\right],
\]
where \(\rho(\mathbf q)\) is the scalar probability density and \(P(\mathbf q)\) the coordinate PDF. This yields a free energy \(F(\mathbf q)\) that is itself a scalar under point transformations \(Q^i=Q^i(\mathbf q)\) [2604.22371].

The requirement of covariance is not ornamental. In reduced descriptions of soft and active matter, the relevant variables are often curvilinear or collective coordinates, such as angles, shape modes, surface coordinates, or defect coordinates. A noncovariant formulation can depend spuriously on the chosen parametrization. The covariant active Onsager–Machlup construction therefore insists that the dissipation function, free-energy rate, active power, and resulting equations all transform consistently under smooth one-to-one coordinate changes [2604.22371].

A complementary non-equilibrium extension shifts emphasis away from equilibrium and toward stationarity. In that formulation, the essential generalization is from a globally stationary Ornstein–Uhlenbeck process to a globally non-stationary but locally stationary, piecewise-stationary stochastic process. The authors explicitly argue that stationarity, not equilibrium, is the essential mathematical requirement behind the generalized Langevin structure. This does not itself yield an active-matter theory, but it clarifies that equilibrium is not the only admissible reference structure for Onsager–Machlup-type reasoning [2311.00120].

## 2. Active Rayleighian and deterministic phenomenological equations

The deterministic active Onsager principle is formulated through the Rayleighian
\[
R(\mathbf q,\dot{\mathbf q})=\Phi(\mathbf q,\dot{\mathbf q})+\dot F(\mathbf q,\dot{\mathbf q})-\dot W(\mathbf q,\dot{\mathbf q}),
\]
where the dissipation function is quadratic,
\[
\Phi(\mathbf q,\dot{\mathbf q})=\frac{1}{2}\dot q^i\zeta_{ij}(\mathbf q)\dot q^j,
\]
with \(\zeta_{ij}=\zeta_{ji}\) symmetric and positive-definite, the free-energy rate is
\[
\dot F(\mathbf q,\dot{\mathbf q})=\dot q^i\partial_iF(\mathbf q),
\]
and the active or external power is
\[
\dot W(\mathbf q,\dot{\mathbf q})=X_i(\mathbf q)\dot q^i.
\]
The physical velocity is selected by minimizing \(R\) with respect to \(\dot{\mathbf q}\) at fixed \(\mathbf q\),
\[
\left(\frac{\partial R(\mathbf q,\dot{\mathbf q})}{\partial \dot q^i}\right)_{\mathbf q}=0,
\]
which gives the active phenomenological equation
\[
\zeta_{ij}(\mathbf q)\dot q^j+\partial_iF(\mathbf q)-X_i(\mathbf q)=0.
\]
In the passive limit \(X_i=0\), this reduces to the usual relaxational balance between dissipation and thermodynamic force [2604.22371].

This extension is conceptually minimal but structurally decisive. Activity enters linearly as power injection into the Rayleighian, rather than as an effective modification of the free energy. The deterministic active drift is therefore not a gradient flow in general. A plausible implication is that the active Onsager equation preserves the variational architecture of passive soft-matter theory while abandoning the equilibrium identification of drift with free-energy descent.

The same logic can be enlarged to inertial dynamics. The generalized Rayleighian becomes
\[
\mathcal R(\mathbf q,\dot{\mathbf q},\mathbf a)=\Phi+\dot F-\dot W+\dot K,
\]
with kinetic energy
\[
K=\frac{1}{2}\dot q^i\mu_{ij}(\mathbf q)\dot q^j
\]
and covariant acceleration
\[
a^i=\ddot q^i+\Gamma^i_{jk}\dot q^j\dot q^k.
\]
Minimization must then be carried out at fixed \(\mathbf q\) and fixed covariant acceleration \(\mathbf a\), which yields
\[
\zeta_{ij}(\mathbf q)\dot q^j+\partial_iF(\mathbf q)-X_i(\mathbf q)+\mu_{ij}(\mathbf q)a^j=0.
\]
The emphasis on \(a^i\), rather than \(\ddot q^i\), is the paper’s criterion for covariance in the inertial sector [2604.22371].

## 3. Stochastic path weight, most probable paths, and fluctuation theorems

Thermal fluctuations are incorporated by defining the Onsager–Machlup Lagrangian as the excess Rayleighian above its local minimum,
\[
L(\mathbf q,\dot{\mathbf q})=R(\mathbf q,\dot{\mathbf q})-R_{\min}(\mathbf q),\qquad
R_{\min}(\mathbf q)=\min_{\dot{\mathbf q}\,|\,\mathbf q}R(\mathbf q,\dot{\mathbf q}),
\]
so that \(L=0\) on the deterministic active drift. The Onsager–Machlup functional is then
\[
O[\mathbf q(t)] = \frac{1}{2k_{\rm B}T}\int_0^\tau dt\,L(\mathbf q(t),\dot{\mathbf q}(t)),
\]
and the conditioned path probability is postulated as
\[
\rho[\mathbf q(t)\mid \mathbf q_0]=N e^{-O[\mathbf q(t)]}.
\]
The most probable trajectory is the minimizer of \(O\), equivalently the extremal of the Euler–Lagrange equation associated with the scalar Lagrangian \(L\) [2604.22371].

A central consistency condition is the detailed fluctuation theorem. Defining entropy production by
\[
\sigma = k_{\rm B}\ln\frac{\rho[\mathbf q(t)]}{\rho[\mathbf q_{\rm rev}(t)]},
\qquad
\mathbf q_{\rm rev}(t)=\mathbf q(\tau-t),
\]
the active Onsager–Machlup construction yields
\[
\sigma=\frac{W}{T}-\frac{\Delta F}{T}+\Delta S_{\mathrm{mes}}.
\]
The same expression is obtained from thermodynamic bookkeeping, so the path-probability formulation is explicitly consistent with stochastic thermodynamics. In this sense, active dynamics break equilibrium detailed balance without invalidating the fluctuation-theorem structure of forward and backward path weights [2604.22371].

For non-reciprocal active systems, a thermodynamically consistent Markov framework based on microscopic local detailed balance provides closely related ingredients even though it does not print a final Onsager–Machlup functional in the standard diffusion form. At the coarse-grained level, the deterministic evolution is identified as the “most likelihood path” of the Doi–Peliti action, namely the “Instanton,” and the mean entropy production decomposes into reciprocal, non-reciprocal, chemical, and self-propulsion sectors. In particular, the non-reciprocal contribution takes the form
\[
\langle \dot{\Sigma}^{nr} \rangle = \int_{\mathcal{V}} \sum_{i,j} V_{ij}^{nr}\langle \omega_{ij}\rangle,
\qquad
\omega_{ij}=\rho_i\partial_t\rho_j-\rho_j\partial_t\rho_i,
\]
so dissipation is directly tied to vorticity currents in macrostate space [2504.10515]. This suggests that, in active non-reciprocal matter, the time-antisymmetric part of an Onsager–Machlup action is naturally organized by circulating currents rather than by free-energy relaxation alone.

## 4. Hidden active variables, nonreciprocity, and odd elasticity

A distinct but closely related line of work adapts Onsager’s variational principle to active systems with odd elasticity. That formulation is deterministic, overdamped, and coarse-grained; it does **not** formulate stochastic path weights or an Onsager–Machlup action. Its relevance lies in showing how a hidden actively driven coordinate can generate nonreciprocal effective dynamics after elimination [2209.15363].

In the finite-dimensional toy model, the potential energy is
\[
U=\frac{k}{2}\left(u_x^2+u_y^2\right)-fs,
\]
with \(\mathbf u=(u_x,u_y)\) the particle displacement and \(s\) the angular coordinate of a driven disk. The dissipation function is
\[
\Phi = \frac{\zeta_u}{2}\left(\dot u_x^2+\dot u_y^2\right) +\frac{\zeta_s}{2}\dot s^2 +\frac{\xi}{2}\left(u_x\dot s-\dot u_y\right)^2 +\frac{\xi}{2}\left(u_y\dot s+\dot u_x\right)^2,
\]
and the Rayleighian is \(R=\dot U+\Phi\). Minimization with respect to \(\dot u_x,\dot u_y,\dot s\) yields coupled deterministic equations, and eliminating \(\dot s\) produces an effective friction matrix \(\boldsymbol{\Gamma}\) that remains symmetric together with an elastic matrix \(\mathbf E\) whose off-diagonal entries are antisymmetric and proportional to the nonequilibrium drive \(f\). In the small-displacement limit,
\[
(\zeta_u+\xi)
\begin{pmatrix}\dot u_x\\ \dot u_y\end{pmatrix}
=
\begin{pmatrix}
-k & -f\xi/\zeta_s\\
f\xi/\zeta_s & -k
\end{pmatrix}
\begin{pmatrix}u_x\\ u_y\end{pmatrix},
\]
so the odd elastic part is
\[
\mathbf E^{\rm odd}= \frac{f\xi}{\zeta_s}
\begin{pmatrix}0 & -1\\ 1 & 0\end{pmatrix}.
\]
The odd coupling is therefore proportional to the nonequilibrium force \(f\) and the ratio of friction coefficients \(\xi/\zeta_s\) [2209.15363].

The continuum construction makes the same point at field level. With displacement field \(\mathbf u(\mathbf r)\), hidden active variable \(s(\mathbf r)\), and deformation basis
\[
\tau_0 = u_{xx} + u_{yy},\qquad \tau_1 = u_{yx} - u_{xy},\qquad \tau_2 = u_{yy} - u_{xx},\qquad \tau_3 = u_{xy}+u_{yx},
\]
the free energy includes the active term \(-fs\), the dissipation function contains quadratic couplings between \(\dot s\) and strain-rate sectors, and in the small-strain limit one obtains odd elastic moduli
\[
A= - \frac{f \mu}{\zeta_s},\qquad K^{\rm o}= \frac{f \nu}{\zeta_s}.
\]
The equilibrium limit \(f=0\) removes these antisymmetric moduli and restores reciprocal passive dynamics [2209.15363].

The paper explicitly states that it does not discuss Onsager–Machlup path probabilities, stochastic action functionals, large deviations, or path integrals. Conceptually, it sits closer to the deterministic drift structure that a stochastic Onsager–Machlup theory would use. This distinction is important because “Onsager–Machlup principle for active systems” is sometimes used loosely to refer to any active Onsagerian formulation, whereas the odd-elasticity construction remains a deterministic Onsager variational framework rather than a path-probability theory.

## 5. Active Brownian particles and explicit most-probable-path computations

A direct application of the active Onsager–Machlup principle is the 3D active Brownian particle. The particle position is \(\mathbf r(t)=[x(t),y(t),z(t)]\), the self-propulsion velocity is \(\mathbf U(t)=U\mathbf e(t)\), and the orientation unit vector is parametrized by spherical angles,
\[
\mathbf e(t)= [\sin\theta(t)\cos\phi(t),\,  \sin\theta(t)\sin\phi(t),\,  \cos\theta(t)].
\]
The angular velocity satisfies \(\dot{\mathbf e}(t)=\bm\omega(t)\times \mathbf e(t)\) with \(\bm\omega\cdot\mathbf e=0\), hence
\[
\bm\omega^2=\dot{\theta}^2+\dot{\phi}^2\sin^2\theta.
\]
For a free active Brownian particle the Rayleighian is
\[
R=\Phi-\dot W_{\rm a},
\]
with
\[
\Phi=\frac{\zeta_{\rm t}}{2}\mathbf v^2+\frac{\zeta_{\rm r}}{2}\bm\omega^2,
\qquad
\dot W_{\rm a}=\zeta_{\rm t}\mathbf U\cdot \mathbf v.
\]
Minimization gives the deterministic active path
\[
\mathbf v=\mathbf U,\qquad \bm\omega=\mathbf 0,
\]
and
\[
R_{\min}=-\frac{\zeta_{\rm t}U^2}{2}.
\]
The Onsager–Machlup integral therefore becomes
\[
O[\mathbf r(t),\mathbf e(t)] = \int_0^{t_{\rm f}}dt\,
\left[
\frac{\zeta_{\rm t}}{2}(\mathbf v-\mathbf U)^2
+
\frac{\zeta_{\rm r}}{2}\bm\omega^2
\right],
\]
with path probability
\[
P[\mathbf r(t),\mathbf e(t)\,|\,\mathbf r_{\rm i},\mathbf e_{\rm i}]
=
C\exp\!\left[-\frac{O[\mathbf r(t),\mathbf e(t)]}{2k_{\rm B}T}\right].
\]
The most probable path is the minimizer of \(O\) under the Dirichlet boundary conditions \(\mathbf r(0)=\mathbf r_{\rm i}\), \(\mathbf e(0)=\mathbf e_{\rm i}\), \(\mathbf r(t_{\rm f})=\mathbf r_{\rm f}\), and \(\mathbf e(t_{\rm f})=\mathbf e_{\rm f}\) [2511.16178].

In dimensionless variables,
\[
L=\sqrt{\frac{\zeta_{\rm r}}{\zeta_{\rm t}}},
\qquad
\overline{\mathbf r}=\frac{\mathbf r}{L},
\qquad
\overline t=\frac{tU}{L},
\qquad
{\rm Pe}=\frac{U}{D_{\rm r}L},
\qquad
D_{\rm r}=\frac{k_{\rm B}T}{\zeta_{\rm r}},
\]
the action is
\[
\frac{O}{k_{\rm B}T}
=
\frac{\rm Pe}{2}
\int_0^{\overline t_{\rm f}} d\overline t\,
\Big[
(\dot{\overline x}-\sin\theta\cos\phi)^2
+
(\dot{\overline y}-\sin\theta\sin\phi)^2
+
(\dot{\overline z}-\cos\theta)^2
+
\dot\theta^2
+
\dot\phi^2\sin^2\theta
\Big].
\]
The Euler–Lagrange equations are nonlinear coupled boundary-value equations, so the paper minimizes the action directly with a neural-network ansatz. The resulting most probable paths exhibit geometric transitions between in-plane **I-paths**, planar **U-paths**, fully 3D **H-paths**, and, for winding boundary conditions, planar \(\ell\)-paths. The transition structure is diagnosed by the curvature and torsion,
\[
\kappa(t)=\frac{|\dot{\mathbf r}\times \ddot{\mathbf r}|}{|\dot{\mathbf r}|^3},
\qquad
\tau(t)=\frac{(\dot{\mathbf r}\times \ddot{\mathbf r})\cdot \dddot{\mathbf r}}{|\dot{\mathbf r}\times \ddot{\mathbf r}|^2},
\]
and by direct comparison of the Onsager–Machlup integral [2511.16178].

The same work also gives the thermal-bath entropy change through local detailed balance,
\[
e^{\Delta S_{\rm b}[\mathbf r(t),\mathbf e(t)]/k_{\rm B}}
=
\frac{P[\mathbf r(t),\mathbf e(t)|\mathbf r_{\rm i},\mathbf e_{\rm i}]}
{P[\mathbf r^{\rm rev}(t),\mathbf e^{\rm rev}(t)|\mathbf r_{\rm f},\mathbf e_{\rm f}]},
\]
with
\[
\mathbf r^{\rm rev}(t)=\mathbf r(t_{\rm f}-t),\qquad
\mathbf e^{\rm rev}(t)=\mathbf e(t_{\rm f}-t).
\]
For the 3D active Brownian particle,
\[
\frac{\Delta S_{\rm b}}{k_{\rm B}}
=
{\rm Pe}
\int_0^{\overline t_{\rm f}} d\overline t\,
\left[
\dot{\overline x}\sin\theta\cos\phi
+
\dot{\overline y}\sin\theta\sin\phi
+
\dot{\overline z}\cos\theta
\right].
\]
This makes the path action, the most-probable-path problem, and the irreversibility measure part of the same active Onsager–Machlup structure [2511.16178].

## 6. Broader generalizations, open directions, and current limitations

A broader non-equilibrium Onsager–Machlup theory has been proposed in which the classical globally stationary Ornstein–Uhlenbeck backbone is replaced by a globally non-stationary but locally stationary, piecewise-stationary stochastic process. The mean state evolves according to
\[
\frac{d\overline{\mathbf a}(t)}{dt}=\mathcal R[\overline{\mathbf a}(t)]
=\mathcal L[\overline{\mathbf a}(t)]\cdot \Delta \mathbf F[\overline{\mathbf a}(t)],
\]
while fluctuations around the evolving mean obey, at fixed \(t\), a stationary non-Markov equation in the fluctuation time \(\tau\),
\[
\frac{d \delta {\textbf a}(t+\tau)}{d\tau}
=
\int_0^\tau {\textbf H}(\tau-\tau';t) \cdot \delta {\textbf a}(t+\tau') d\tau'
+
{\textbf f}(t+\tau).
\]
The framework allows memory, additive noise, local fluctuation-dissipation structure, and a Lyapunov function \(\Lambda[\mathbf a]\), identified in thermodynamic applications with \(-S[\mathbf a]/k_B\). It does **not** provide an explicit Onsager–Machlup path action, and it does **not** explicitly treat active matter, but it clarifies how a locally stationary stochastic description can extend Onsager–Machlup reasoning beyond equilibrium [2311.00120].

The limitations of that extension are also clear and particularly relevant for active systems. The paper itself notes that active matter often lacks a thermodynamic entropy \(S[\mathbf a]\) whose negative is a Lyapunov function; active noise is usually non-thermal; self-propulsion variables are missing; non-reciprocal couplings and irreversible circulating currents are not developed explicitly; multiplicative noise is absent; and time reversal and entropy production are not formulated pathwise [2311.00120]. These are not peripheral issues. They identify the precise points at which a generic non-equilibrium Onsager–Machlup scaffold ceases to be a ready-made active-matter theory.

A different statistical extension reformulates the Onsager–Machlup variational principle as a tool for cumulant generating functions. Starting from
\[
O[\mathbf v(\mathbf r,t),\mathbf x(\mathbf r,t)]
=
\int_0^\tau dt\,\bigl(R[\mathbf v(\mathbf r,t),\mathbf x(\mathbf r,t)]-R_{\min}[\mathbf x(\mathbf r,t)]\bigr),
\]
it introduces the modified functional
\[
\Omega[\mathbf v(\mathbf r,t),\mathbf x(\mathbf r,t)]
=
qX-\frac{O[\mathbf v(\mathbf r,t),\mathbf x(\mathbf r,t)]}{2k_{\rm B}T}+\Gamma,
\]
whose maximum approximates the cumulant generating function \(K_X(q)\). That construction is demonstrated for a Brownian particle in a viscous fluid and in steady shear flow, and the authors explicitly propose the active extension
\[
R=\Phi+\dot F+\dot W+C,
\]
where \(\dot W\) is the power generated by active elements due to energy input or chemical reactions. The same paper also states that the present statistical formulation “cannot deal with active non-thermal fluctuations” [2404.18611]. This is therefore a formal route toward active Onsager–Machlup theory, but not yet a complete treatment of canonical active-noise models.

Taken together, the current literature supports a relatively precise statement. The active Onsager–Machlup principle is now explicit for thermally fluctuating active systems whose mesoscopic dynamics can be written through an active Rayleighian with direct power injection, and it has been worked out in detail for active Brownian particles [2604.22371]. Deterministic Onsagerian treatments with hidden active variables show how nonreciprocity and odd elasticity can emerge by elimination, but they are not themselves stochastic Onsager–Machlup theories [2209.15363]. More general non-equilibrium formulations based on local stationarity, statistical variational principles, Doi–Peliti coarse-graining, or non-reciprocal stochastic thermodynamics provide adjoining structures rather than a single universal active Onsager–Machlup formalism [2311.00120]. A plausible implication is that the subject is converging on a family of active Onsager–Machlup principles rather than on one unique canonical action, with the decisive distinctions set by covariance, hidden active coordinates, noise statistics, and the thermodynamic status of the driving.

Source: https://www.emergentmind.com/topics/onsager-machlup-principle-for-active-systems