Metriplectic Bracket Formalism

• Metriplectic bracket formalism is an algebraic framework combining a Poisson bracket for reversible dynamics with a symmetric bracket for irreversible dissipation.
• It conserves the Hamiltonian while ensuring monotonic entropy production or free-energy decrease under appropriate degeneracy conditions.
• The framework extends to Lie–Poisson systems, field theories, and higher-bracket geometries, supporting structure-preserving computation and data-driven modeling.

Searching arXiv for recent and foundational papers on metriplectic bracket formalism. Metriplectic bracket formalism is an algebraic framework for dissipative dynamics that combines a reversible Hamiltonian component with an irreversible dissipative component in a single Leibniz-type evolution law. In its classical form, the reversible part is generated by a Poisson bracket and the dissipative part by a symmetric positive semidefinite bracket, with compatibility conditions chosen so that the Hamiltonian is conserved and an entropy-like functional is produced monotonically. Across later developments, this basic split has been extended to Lie–Poisson systems, field theories, 4-bracket and triple-bracket constructions, contact and jet-bundle geometries, reduced-order modeling, and data-driven identification, while retaining the central idea that reversible and irreversible effects should be encoded by distinct but coupled bracket structures (Teng et al., 2024, Morrison et al., 2023).

1. Classical algebraic structure

Classical metriplectic dynamics on a smooth manifold $P$ uses four ingredients: a Poisson bracket $\{\cdot,\cdot\}$, a symmetric positive semidefinite bracket $(\cdot,\cdot)$, a Hamiltonian $H$, and an entropy $S$. The Poisson bracket is bilinear, skew-symmetric, satisfies the Leibniz rule and the Jacobi identity, and can be represented by a Poisson bivector $\Pi$. The metric bracket is induced by a symmetric positive semidefinite covector inner product $\mathcal K$, or in coordinates by a symmetric matrix $K$. In this setting,

$$\dot x=\sharp_\Pi(dH)+\sharp_{\mathcal K}(dS)=\Pi\nabla H+K\nabla S.$$

Equivalent formulations write the evolution of an observable $\psi$ as

$\{\cdot,\cdot\}$0

or define a combined Leibniz bracket $\{\cdot,\cdot\}$1 (Teng et al., 2024, Materassi, 2014, Marcucci et al., 2018).

The classical compatibility conditions are the two degeneracies

$\{\cdot,\cdot\}$2

and

$\{\cdot,\cdot\}$3

Thus $\{\cdot,\cdot\}$4 lies in the kernel of the metric bracket, while $\{\cdot,\cdot\}$5 is a Casimir of the Poisson bracket. In the language of complete metriplectic systems, the dissipative bracket redistributes energy without destroying the total Hamiltonian, and the entropy-like generator is invisible to the Hamiltonian part (Teng et al., 2024, Marcucci et al., 2018).

This decomposition is often written in operator form as

$\{\cdot,\cdot\}$6

with $\{\cdot,\cdot\}$7, $\{\cdot,\cdot\}$8, and the algebraic core expressed as $\{\cdot,\cdot\}$9 and $(\cdot,\cdot)$0. That formulation is particularly common in reduced-order modeling and structure-preserving numerics, where skew-symmetry, symmetry, and degeneracy become the primary invariants to preserve (Gruber et al., 2022).

2. Thermodynamic content and free-energy principles

Under the classical degeneracy conditions, the formalism encodes the standard thermodynamic identities

$(\cdot,\cdot)$1

Accordingly, the Hamiltonian is conserved and the entropy is nondecreasing. Several formulations package these generators into a free-energy functional, written for example as $(\cdot,\cdot)$2, $(\cdot,\cdot)$3, or $(\cdot,\cdot)$4, depending on sign conventions and the choice of entropy variable (Materassi et al., 2011, Marcucci et al., 2018, Bosboom et al., 2023).

In this setting, the entropy-like Casimir functions as a Lyapunov functional. The MHD formulation states the corresponding $(\cdot,\cdot)$5-theorem as

$(\cdot,\cdot)$6

while the radiative-transfer formulation obtains

$(\cdot,\cdot)$7

for its chosen sign convention and negative semidefinite metric bracket. Equilibria are characterized by stationarity of the free energy, typically through $(\cdot,\cdot)$8, and the resulting asymptotic states are static thermodynamic equilibria or absorbed states depending on the model (Materassi et al., 2011, Bosboom et al., 2023, Marcucci et al., 2018).

The classical interpretation is especially transparent in examples where the dissipative channel transfers energy among subsystems while preserving the total Hamiltonian. In the radiation–matter interaction toy model, the field variables $(\cdot,\cdot)$9 form the conservative subsystem and the excited-state population $H$0 acts as an environmental variable that absorbs energy irreversibly; the metric bracket is constructed so that the total Hamiltonian is preserved while the entropy-like functional increases monotonically until the absorbed state is reached (Marcucci et al., 2018).

3. Generalizations beyond the classical Casimir picture

The classical dictum “energy conserved, entropy produced” is not universal within metriplectic theory. A generalized metriplectic system can be obtained by relaxing the requirement that the entropy be a Casimir of the Poisson bracket and relaxing the condition that the metric annihilate $H$1 in the classical way. In that case one still writes

$H$2

but $H$3 is allowed to participate in the Poisson term. The central consequence is the free-energy law

$H$4

Here $H$5 is interpreted as a free-energy analogue, so the monotone object is no longer mechanical energy or entropy alone, but the combined quantity $H$6 (Teng et al., 2024).

This generalization permits nonconservative dissipation that does not preserve mechanical energy, and the 2D and $H$7 examples in that work show trajectories converging to zero level sets of a prescribed free energy. In the 2D example, the zero level set is described as $H$8-shaped; in the $H$9 example it is a sphere. The stated thermodynamic claim is that monotonic change of free energy can serve as a more precise criterion than mechanical energy or entropy alone (Teng et al., 2024).

Open-system extensions push the departure from the classical scheme further. For gravitational subsystems, the symmetric bracket is introduced as a boundary term representing charge loss due to radiation. On the resulting metriplectic space, Komar charges are Hamiltonian, yet they are not conserved under their own flow: $S$0 This construction is explicitly motivated by the statement that open gravitational dynamics are not symplectic in the strict sense, because a bounded region exchanges charge with its environment through the boundary (Kabel et al., 2022).

A related distinction appears in contact geometry. Standard contact Hamiltonian systems on $S$1 generally satisfy

$S$2

so thermodynamic consistency is not automatic. By contrast, the metriplectic system constructed on the same one-jet bundle enforces

$S$3

with the coordinate $S$4 serving as entropy. This suggests that contact dissipation and metriplectic thermodynamic consistency are distinct structures rather than interchangeable descriptions (Morrison et al., 10 May 2026).

4. Higher-bracket and geometric formulations

A major line of development replaces the symmetric bilinear dissipative bracket by higher-order brackets. In metriplectic 4-bracket dynamics, one defines

$S$5

with symmetries modeled on the Riemann curvature tensor: antisymmetry in the first pair, antisymmetry in the second pair, symmetry under exchange of the pairs, and a cyclic/Bianchi-type identity. Freezing the Hamiltonian twice produces the induced dissipative 2-bracket

$S$6

which is automatically symmetric and satisfies $S$7, while entropy production is given by

$S$8

under the positivity condition on the sectional-curvature-like quantity (Morrison et al., 2023).

This curvature-like viewpoint is implemented through Kulkarni–Nomizu constructions and related geometric mechanisms. The 4-bracket formalism is explicitly described as more inclusive than earlier binary-bracket approaches and is said to include all known previous binary bracket theories for dissipation or relaxation as special cases. It also recovers Kaufman–Morrison brackets, double-bracket dynamics, and a symmetrized or linearized version of GENERIC (Morrison et al., 2023).

In continuum models, the Morrison–Updike 4-bracket is used as an algorithmic device for dissipative transport. For Cahn–Hilliard–Navier–Stokes systems, the 4-bracket is multilinear, antisymmetric within each pair, symmetric under pair exchange, and built so that $S$9. The dissipative evolution

$\Pi$0

then yields thermodynamically consistent viscosity, heat conduction, and diffusion, and naturally supports anisotropic surface energy effects (Zaidni et al., 2024).

Variational thermodynamics provides another route to the same structure. Constrained Lagrangian variational principles for non-equilibrium thermodynamics can be rewritten as

$\Pi$1

with the dissipative 4-bracket constructed directly from friction laws and, in several cases, given a Kulkarni–Nomizu interpretation. The same program covers simple systems, discrete systems, Euler–Poincaré reduced systems, and systems with no symplectic part (Carlier, 2024).

An alternative 4-bracket construction defines irreversibility as the product of two skew-symmetric brackets. The induced 2-bracket then acquires symmetry, degeneracy, and positive semidefiniteness automatically: $\Pi$2 That formulation is applied to generalized 2D quasigeostrophic ocean models with advected quantities, where internal energy is conserved and entropy is generated (Beron-Vera et al., 24 Jun 2026).

A related but distinct algebraic route derives metriplectic systems from completely antisymmetric triple brackets on quadratic Lie algebras. In that construction the symmetric dissipative bracket is generated by freezing one slot to the Hamiltonian, yielding generalized double-bracket dissipation and recovering examples such as relaxing rigid-body dynamics and dissipative Toda flows (Bloch et al., 2012).

5. Representative applications

The formalism has been applied to a wide range of finite- and infinite-dimensional systems in which one wants explicit separation of reversible transport from irreversible relaxation.

Domain	Variables or structure	Encoded dissipative content
Visco-resistive MHD	$\Pi$3	Viscosity, resistivity, heat transport, Galilean invariance, asymptotic equilibria (Materassi et al., 2011)
Dissipative XMHD	$\Pi$4	Thermophoresis, current viscosity, viscous cross couplings, nonlocal Lagrangian dissipative brackets (Coquinot et al., 2019)
CHNS systems	$\Pi$5 or $\Pi$6	Viscosity, heat conduction, diffusion, anisotropic surface energy (Zaidni et al., 2024)
Polarized radiative transfer	Coherence matrix $\Pi$7	Scattering dissipation with Jacobi-satisfying conservative bracket after basis transformation (Bosboom et al., 2023)
Guiding-center Vlasov–Maxwell–Landau	$\Pi$8	Collisional Landau bracket with energy-momentum, angular-momentum conservation and $\Pi$9-theorem (Brizard et al., 27 Jun 2025)
Rigid body control	Lie–Poisson algebra of $\mathcal K$0	Energy-preserving servo torque driving alignment to a stable principal axis (Materassi et al., 2018)

In fluid mechanics, the formalism is often used to clarify the status of entropy. The Lagrangian formulation of viscous fluids rewrites the metriplectic algebra in parcel variables $\mathcal K$1, where $\mathcal K$2 and $\mathcal K$3 describe macroscopic parcel motion and $\mathcal K$4 represents microscopic or internal degrees of freedom. In that representation, entropy is described as a C2-type Casimir because it does not appear in the Poisson bracket at all (Materassi, 2014).

In optics and radiation theory, metriplectic brackets encode irreversible transfer and scattering without abandoning Hamiltonian transport. The two-photon absorption toy model is a complete metriplectic system in which radiation variables decay while atomic excitation grows toward an asymptotically stable absorbed state (Marcucci et al., 2018). The polarized radiative-transfer formulation uses a matrix-valued Poisson-like bracket that satisfies the Jacobi identity after transforming to a polarization basis in which the optical rotation term vanishes, and the metric bracket is negative semidefinite so that the first two laws of thermodynamics hold automatically (Bosboom et al., 2023).

Kinetic and plasma applications emphasize conservation of collision invariants under dissipation. In guiding-center Vlasov–Maxwell–Landau theory, the symmetric Landau bracket reproduces the collisional guiding-center Fokker–Planck operator while conserving guiding-center energy-momentum and angular momentum, and satisfying a guiding-center $\mathcal K$5-theorem (Brizard et al., 27 Jun 2025). In resistive MHD, the symmetric bracket extends the Navier–Stokes dissipative bracket to magnetic-field dynamics and produces Joule heating, viscous heating, and thermal conduction while preserving the total energy of the closed plasma system (Materassi et al., 2011).

6. Structure-preserving computation, reduction, and identification

Because the formalism is defined by algebraic identities rather than only by differential equations, numerical and reduced models are often designed to preserve those identities explicitly. For metriplectic Euler–Poincaré equations, discrete gradient methods are used to construct integrators that preserve energy exactly and produce the correct entropy production rate. In the rigid-body example, the midpoint scheme preserves the Hamiltonian exactly and satisfies the expected monotonicity of the quadratic entropy/Casimir (Bloch et al., 2024).

Reduced-order modeling poses a specific challenge because naive projection generally breaks the degeneracy conditions. A POD-based reduced model can nevertheless be made metriplectic-preserving by projecting tensor representations of the skew and symmetric operators in a way that preserves the algebraic structure. The resulting reduced system satisfies

$\mathcal K$6

and the method is stated to be provably convergent. A caveat is also stated explicitly: if the Poisson operator $\mathcal K$7 depends on the state, the reduced Poisson bracket may fail to satisfy the Jacobi identity even though the thermodynamic structure remains intact (Gruber et al., 2022).

Data-driven identification has also been formulated directly in bracket language. In the generalized metriplectic setting, a bilevel convex optimization method identifies $\mathcal K$8 and $\mathcal K$9 from trajectory data under the assumption that the Poisson structure $K$0 is known. The optimization enforces $K$1 and bounded polynomial parameterizations, and is built around the generalized form

$K$2

The stated aim is to recover the bracket-based decomposition from measurements, not to redefine the dynamics (Teng et al., 2024).

A broader algorithmic program appears in the unified thermodynamic algorithm, which turns the 4-bracket formalism into a construction method. Its basic pattern is: choose variables, choose Hamiltonian and entropy functionals, construct a noncanonical Poisson bracket with entropy as a Casimir, and then build a metriplectic 4-bracket from force–flux relations. This program is applied to Navier–Stokes–Fourier, Cahn–Hilliard–Navier–Stokes, and Brenner–Navier–Stokes–Fourier systems, yielding thermodynamically consistent generalizations (Zaidni et al., 2024).

7. Conceptual scope, limitations, and neighboring frameworks

Metriplectic formalism is sometimes treated as synonymous with “Hamiltonian plus friction,” but the literature gives a narrower and more structured definition. The dissipative bracket is not arbitrary damping: it is constrained by symmetry or higher-bracket symmetries, semidefiniteness, and degeneracy relative to the Hamiltonian or free-energy generators. This is why the formalism is repeatedly used to encode exact energy conservation together with monotone entropy production in closed systems (Materassi et al., 2017, Materassi et al., 2011).

At the same time, the formalism is not confined to one universal recipe. Classical complete metriplectic systems use entropy as a Poisson Casimir and Hamiltonian degeneracy of the metric bracket. Generalized metriplectic systems replace that with monotone free-energy dynamics, while gravitational subsystem models treat the symmetric term as genuine flux through a boundary and therefore allow the generator itself to change under its own flow (Teng et al., 2024, Kabel et al., 2022). A plausible implication is that the term “metriplectic” now denotes a family of algebraically related thermodynamic structures rather than a single fixed bracket ansatz.

Relations to neighboring frameworks are likewise nontrivial. The 4-bracket literature presents standard metriplectic binary brackets as special cases, and states that a symmetrized or linearized version of GENERIC can be reconstructed within the 4-bracket framework, while raw GENERIC is not automatically in the metriplectic class (Morrison et al., 2023). Contact dynamics, although dissipative, does not in general guarantee $K$3 and $K$4, which is precisely the distinction emphasized by the metriplectic construction on $K$5 (Morrison et al., 10 May 2026).

The recurrent conceptual theme is therefore not merely the presence of dissipation, but the algebraic organization of dissipation. Poisson brackets encode reversible transport, symmetric or curvature-like brackets encode irreversible relaxation, and the chosen degeneracies determine whether the natural monotone object is entropy, free energy, or a boundary-sensitive quasi-local charge.