GENERIC SDE: Geometric Formulation of Stochastic Dynamics

• GENERIC SDE is a coordinate-invariant formulation uniting reversible Hamiltonian dynamics with irreversible dissipative processes, ensuring almost-sure energy conservation.
• It employs a degenerate Poisson structure, a dissipative co-metric, and a unimodular volume form to preserve Boltzmann-type invariant measures.
• The framework supports structure-preserving numerical methods and bridges Itô and Stratonovich forms via geometric insights.

The GENERIC stochastic differential equation is a coordinate-invariant geometric formulation that unifies reversible Hamiltonian and irreversible dissipative dynamics on a smooth finite-dimensional manifold. Its defining ingredients are a degenerate Poisson structure, a degenerate co-metric, and a smooth positive volume form satisfying a unimodularity condition. In this formulation, the stochastic dynamics preserves a Boltzmann-type measure, conserves energy almost surely, and reduces in the zero-noise limit to the deterministic metriplectic formulation (Peletier et al., 11 Sep 2025).

1. Geometric data and non-interaction structure

The geometric formulation starts with a smooth finite-dimensional manifold $M$, together with smooth energy and entropy functions

$E:M\to \mathbb{R}, \qquad S:M\to \mathbb{R}.$

The reversible part of the dynamics is encoded by a possibly degenerate Poisson bivector

$\Pi \in \Gamma(\Lambda^2 TM),$

while the irreversible part is encoded by a possibly degenerate dissipative co-metric

$K \in \Gamma(S^2 TM),$

which is symmetric and positive semidefinite (Peletier et al., 11 Sep 2025).

These tensors define two brackets on $C^\infty(M)$: $$\{f,g\}:=\langle df,\Pi\,dg\rangle, \qquad (f,g):=\langle df,K\,dg\rangle.$$ The first is the Poisson bracket associated with reversible Hamiltonian motion, and the second is the symmetric dissipative bracket associated with irreversible relaxation. The two structures are coupled by the degeneracy, or non-interaction, conditions

$\Pi\,dS=0, \qquad K\,dE=0.$

These express that entropy is a Casimir of the Poisson structure and energy is a dissipation Casimir of the metric structure (Peletier et al., 11 Sep 2025).

A further structural ingredient is a smooth positive volume form $\mu$ on $M$. Divergence with respect to $\mu$ is defined by

$E:M\to \mathbb{R}, \qquad S:M\to \mathbb{R}.$0

and the unimodularity condition requires that every Hamiltonian vector field $E:M\to \mathbb{R}, \qquad S:M\to \mathbb{R}.$1 be $E:M\to \mathbb{R}, \qquad S:M\to \mathbb{R}.$2-divergence-free: $E:M\to \mathbb{R}, \qquad S:M\to \mathbb{R}.$3 This is equivalent to vanishing of the modular class with respect to $E:M\to \mathbb{R}, \qquad S:M\to \mathbb{R}.$4, and ensures that the reversible Hamiltonian part preserves the reference volume form (Peletier et al., 11 Sep 2025).

2. Coordinate-free stochastic dynamics

The stochastic part of the construction is built from the horizontal distribution

$E:M\to \mathbb{R}, \qquad S:M\to \mathbb{R}.$5

The co-metric $E:M\to \mathbb{R}, \qquad S:M\to \mathbb{R}.$6 induces a sub-Riemannian metric $E:M\to \mathbb{R}, \qquad S:M\to \mathbb{R}.$7 on $E:M\to \mathbb{R}, \qquad S:M\to \mathbb{R}.$8 through

$E:M\to \mathbb{R}, \qquad S:M\to \mathbb{R}.$9

The associated horizontal gradient is

$\Pi \in \Gamma(\Lambda^2 TM),$0

and the sub-Laplacian relative to $\Pi \in \Gamma(\Lambda^2 TM),$1 is

$\Pi \in \Gamma(\Lambda^2 TM),$2

Locally, if $\Pi \in \Gamma(\Lambda^2 TM),$3 is an orthonormal horizontal frame for $\Pi \in \Gamma(\Lambda^2 TM),$4, then

$\Pi \in \Gamma(\Lambda^2 TM),$5

with $\Pi \in \Gamma(\Lambda^2 TM),$6, and

$\Pi \in \Gamma(\Lambda^2 TM),$7

(Peletier et al., 11 Sep 2025).

The noise vector fields are chosen as

$\Pi \in \Gamma(\Lambda^2 TM),$8

and the drift is

$\Pi \in \Gamma(\Lambda^2 TM),$9

The manifold GENERIC SDE is then written in Stratonovich form as

$K \in \Gamma(S^2 TM),$0

Equivalently, if $K \in \Gamma(S^2 TM),$1 is a noise tensor satisfying

$K \in \Gamma(S^2 TM),$2

then

$K \in \Gamma(S^2 TM),$3

(Peletier et al., 11 Sep 2025).

The use of Stratonovich form is geometrically natural because it is coordinate-invariant. More generally, the distinction between Itô and Stratonovich formulations is essential in stochastic analysis, and numerical methods depend on that choice (Gevorkyan et al., 2016). In the GENERIC setting, the equivalent Itô form is

$K \in \Gamma(S^2 TM),$4

where the additional drift is precisely the Stratonovich-to-Itô correction (Peletier et al., 11 Sep 2025).

A crucial structural property is that the stochastic perturbations preserve energy. Since $K \in \Gamma(S^2 TM),$5, one has $K \in \Gamma(S^2 TM),$6, hence

$K \in \Gamma(S^2 TM),$7

Thus the noise acts only along directions compatible with the energy constraint (Peletier et al., 11 Sep 2025).

3. Generator and forward equation

The infinitesimal generator of the GENERIC SDE acts on $K \in \Gamma(S^2 TM),$8 as

$K \in \Gamma(S^2 TM),$9

(Peletier et al., 11 Sep 2025).

If $C^\infty(M)$0 denotes the density with respect to $C^\infty(M)$1, the forward Kolmogorov, or Fokker–Planck, equation is

$C^\infty(M)$2

This can be rewritten as

$C^\infty(M)$3

The reversible term is purely transport by the Hamiltonian vector field, while the irreversible and stochastic terms appear through the co-metric $C^\infty(M)$4 and the diffusion contribution $C^\infty(M)$5 (Peletier et al., 11 Sep 2025).

This structure is a specialized geometric instance of the standard SDE-to-PDE correspondence. For a generic Itô SDE

$C^\infty(M)$6

the forward density evolves by the usual Fokker–Planck equation with first-order advection and second-order diffusion terms (Gevorkyan et al., 2016). The GENERIC formulation imposes additional algebraic and geometric constraints on the drift, diffusion, and invariant volume, rather than treating them as unrelated coefficients.

Unimodularity plays a decisive role in the forward equation. Because

$C^\infty(M)$7

the reversible Hamiltonian part does not create or destroy $C^\infty(M)$8-volume; it induces $C^\infty(M)$9-measure-preserving transport along the Hamiltonian flow (Peletier et al., 11 Sep 2025).

4. Invariant measure, energy conservation, and zero-noise limit

A central result of the geometric formulation is the existence of a family of invariant measures of Boltzmann type: $$\{f,g\}:=\langle df,\Pi\,dg\rangle, \qquad (f,g):=\langle df,K\,dg\rangle.$$0 where $$\{f,g\}:=\langle df,\Pi\,dg\rangle, \qquad (f,g):=\langle df,K\,dg\rangle.$$1 is any smooth nonnegative function (Peletier et al., 11 Sep 2025). The dependence on $$\{f,g\}:=\langle df,\Pi\,dg\rangle, \qquad (f,g):=\langle df,K\,dg\rangle.$$2 reflects the fact that the dynamics conserves energy almost surely, so invariant measures are supported on energy level sets.

Almost-sure energy conservation follows directly from the structural conditions. For solutions $$\{f,g\}:=\langle df,\Pi\,dg\rangle, \qquad (f,g):=\langle df,K\,dg\rangle.$$3,

$$\{f,g\}:=\langle df,\Pi\,dg\rangle, \qquad (f,g):=\langle df,K\,dg\rangle.$$4

Indeed, the reversible drift $$\{f,g\}:=\langle df,\Pi\,dg\rangle, \qquad (f,g):=\langle df,K\,dg\rangle.$$5, the dissipative drift $$\{f,g\}:=\langle df,\Pi\,dg\rangle, \qquad (f,g):=\langle df,K\,dg\rangle.$$6, and all noise vector fields are annihilated by $$\{f,g\}:=\langle df,\Pi\,dg\rangle, \qquad (f,g):=\langle df,K\,dg\rangle.$$7 (Peletier et al., 11 Sep 2025).

The deterministic limit $$\{f,g\}:=\langle df,\Pi\,dg\rangle, \qquad (f,g):=\langle df,K\,dg\rangle.$$8 eliminates both the stochastic forcing and the Itô correction, leaving the metriplectic flow

$$\{f,g\}:=\langle df,\Pi\,dg\rangle, \qquad (f,g):=\langle df,K\,dg\rangle.$$9

Along smooth solutions of this limit system,

$\Pi\,dS=0, \qquad K\,dE=0.$0

while

$\Pi\,dS=0, \qquad K\,dE=0.$1

Thus the deterministic dynamics conserves energy and has non-decreasing entropy, matching the thermodynamic interpretation of the formalism (Peletier et al., 11 Sep 2025).

The fluctuation–dissipation relation is encoded geometrically by

$\Pi\,dS=0, \qquad K\,dE=0.$2

This relation ties the diffusion tensor directly to the dissipative co-metric, so the stochastic forcing is not an independent addition but is structurally matched to the irreversible part of the dynamics (Peletier et al., 11 Sep 2025).

5. Euclidean form and a representative Lie–Poisson example

On $\Pi\,dS=0, \qquad K\,dE=0.$3, the geometric data reduce to matrix fields. The Poisson structure becomes an antisymmetric matrix field $\Pi\,dS=0, \qquad K\,dE=0.$4, the co-metric becomes a symmetric positive semidefinite matrix field $\Pi\,dS=0, \qquad K\,dE=0.$5, and the noise amplitude is a matrix $\Pi\,dS=0, \qquad K\,dE=0.$6 satisfying

$\Pi\,dS=0, \qquad K\,dE=0.$7

The Stratonovich GENERIC SDE takes the form

$\Pi\,dS=0, \qquad K\,dE=0.$8

subject to the degeneracy conditions

$\Pi\,dS=0, \qquad K\,dE=0.$9

(Peletier et al., 11 Sep 2025).

The corresponding Itô equation acquires the correction $\mu$0, so the Itô drift is

$\mu$1

The forward equation with respect to Lebesgue measure is

$\mu$2

and the invariant density family is

$\mu$3

(Peletier et al., 11 Sep 2025).

A representative example is the Lie–Poisson rigid-body-type manifold on $\mu$4, with bracket

$\mu$5

Choosing

$\mu$6

one has $\mu$7 whenever $\mu$8 depends only on the Casimir $\mu$9. If $M$0 is chosen symmetric positive semidefinite with $M$1, and $M$2 satisfies $M$3, then the noise preserves energy. For the linear Poisson structure on $M$4, Hamiltonian vector fields are divergence-free with respect to Lebesgue volume, so unimodularity holds. The invariant density is then

$M$5

(Peletier et al., 11 Sep 2025).

This example exhibits the intended division of labor in the formalism: the Poisson structure governs reversible coadjoint motion, the co-metric governs dissipative relaxation subject to the energy constraint, and the noise acts compatibly with both.

6. Numerical implications, scope, and extensions

The geometric formulation is designed to support structure-preserving numerical methods. The paper proposes splitting the dynamics into a reversible Hamiltonian step and an irreversible metric/noise step. For the reversible part, one may use symplectic or Poisson-preserving integrators that preserve $M$6, or a discrete analogue, together with the energy constraint. For the irreversible stochastic part, one may use stochastic integrators consistent with

$M$7

and respecting the degeneracy condition $M$8, for example by projecting noise increments onto the horizontal distribution $M$9 (Peletier et al., 11 Sep 2025).

To preserve the invariant measure

$\mu$0

at the discrete level, the paper further suggests combining such splitting with exact preservation of energy shells and a stochastic step that is reversible with respect to $\mu$1, for example through Metropolis–Hastings corrections on energy shells (Peletier et al., 11 Sep 2025).

The scope of the formulation is explicitly geometric and coordinate-free. It separates system-specific inputs,

$\mu$2

from ambient geometric data,

$\mu$3

and thereby clarifies which features belong to the physical model and which belong to the surrounding differential-geometric structure (Peletier et al., 11 Sep 2025). A plausible implication is that this separation is particularly useful when transporting the formalism across coordinate systems or comparing models on different manifolds.

The main assumptions are smoothness of $\mu$4; symmetry and positive semidefiniteness of $\mu$5; the Jacobi identity for $\mu$6; the degeneracy conditions $\mu$7 and $\mu$8; and unimodularity of the Hamiltonian part with respect to $\mu$9 (Peletier et al., 11 Sep 2025). Under these conditions, the resulting equation is well-posed in the Stratonovich sense and yields almost-sure energy conservation together with the Boltzmann-type invariant measure.

The framework is presented as a foundation for analytic and numerical methods and for future extensions to coarse-grained and quantum-inspired systems (Peletier et al., 11 Sep 2025). Within the broader SDE literature, this places the GENERIC SDE as a highly structured subclass of stochastic differential equations: unlike generic Itô models, where drift and diffusion are specified directly, the GENERIC construction derives them from coupled reversible, dissipative, and measure-theoretic principles.