Thermodynamic Geometry Insights

• Thermodynamic geometry is a mathematical framework that endows thermodynamic state spaces with a Riemannian metric, enabling coordinate-invariant analysis of fluctuations and dissipation.
• It quantifies minimal-dissipation protocols via optimal control paths (geodesics) that minimize work, offering practical guidelines for traversing phase transitions.
• The approach bridges microscopic behaviors and macroscopic phenomena, using scalar curvature and metric properties to signal stability and critical fluctuations.

Thermodynamic geometry is the mathematical framework that endows the space of thermodynamic states or control parameters with a Riemannian structure, allowing physical properties such as fluctuations, dissipation, stability, and phase transitions to be rigorously analyzed through geometric invariants. This approach generalizes classical thermodynamic and statistical characterizations, offering coordinate-invariant, quantitative criteria for microscopic and macroscopic behavior—including phase transitions, optimal control protocols, stability boundaries, and correlations of fluctuations.

1. Foundational Principles: Metric Structure and Thermodynamic Length

The hallmark of thermodynamic geometry is the construction of a Riemannian metric $g_{\mu\nu}$ on the space of control parameters $\lambda^\mu$, typically defined by equilibrium correlation and linear-response functions. In the linear-response regime, the excess work (“availability”) dissipated under a protocol $\lambda(t)$ over time $T$ is given by the quadratic form: $$A = \int_0^T g_{\mu\nu}(\lambda(t))\,\dot\lambda^\mu(t)\,\dot\lambda^\nu(t)\,dt$$ where

$$g_{\mu\nu}(\lambda) = \beta \int_0^\infty \langle \delta X_\mu(t)\,\delta X_\nu(0)\rangle_{\mathrm{eq}}\,dt = \beta\,T_{\mu\nu}(\lambda)\,\partial_\mu\partial_\nu\ln Z(\lambda)$$

with $\beta = 1/T$ and $X_\mu$ the thermodynamic variables conjugate to $\lambda^\mu$ (Basri et al., 1 Dec 2025).

The associated “thermodynamic length” of a path $\lambda(s)$, $\lambda^\mu$0, is defined as: $\lambda^\mu$1 Through a Cauchy–Schwarz inequality, the dissipated work is bounded below by $\lambda^\mu$2, attaining equality for constant-speed geodesics. The minimization of $\lambda^\mu$3 for fixed $\lambda^\mu$4 yields optimal, minimal-dissipation protocols as geodesics of $\lambda^\mu$5.

2. Divergence Structure near Second-Order Phase Transitions

Approaching critical points (e.g., $\lambda^\mu$6), the equilibrium response functions and relaxation times diverge, leading to divergence of $\lambda^\mu$7. In statistical models described by Widom scaling, the singular part of the free-energy density $\lambda^\mu$8 scales as $\lambda^\mu$9, where $\lambda(t)$0, $\lambda(t)$1, and $\lambda(t)$2 is a scaling function (Basri et al., 1 Dec 2025).

Key scaling relations:

• Static susceptibilities:

$\lambda(t)$3

• Relaxation time: $\lambda(t)$4, for $\lambda(t)$5, with $\lambda(t)$6 and $\lambda(t)$7 the dynamical and spatial exponents.

The metric divergence leads to possible breakdowns of the geometric picture, including non-existence or incompleteness of geodesics.

3. Criteria for Finite-Length Crossing of Critical Points

Despite the divergence of $\lambda(t)$8 at criticality, the integrated thermodynamic length $\lambda(t)$9 may remain finite if the divergence is sufficiently mild. For trajectories approaching the critical point as $T$0, the line element scales as $T$1, with

$T$2

and the length integral converges near $T$3 if $T$4, i.e., $T$5. Similarly, along $T$6, the exponent is $T$7, requiring $T$8. These conditions delimit universality classes in which a finite-length crossing is possible.

Concrete implications:

• $T$9 infinite thermodynamic length (e.g., 2D Ising, Potts)
• $$A = \int_0^T g_{\mu\nu}(\lambda(t))\,\dot\lambda^\mu(t)\,\dot\lambda^\nu(t)\,dt$$0 finite-length crossing possible (e.g., 3D Ising with $$A = \int_0^T g_{\mu\nu}(\lambda(t))\,\dot\lambda^\mu(t)\,\dot\lambda^\nu(t)\,dt$$1) However, even with finite $$A = \int_0^T g_{\mu\nu}(\lambda(t))\,\dot\lambda^\mu(t)\,\dot\lambda^\nu(t)\,dt$$2, the curvature (Riemann tensor) of the thermodynamic manifold generically diverges at the critical point.

4. Numerical Computation of Minimal-Dissipation Paths

When $$A = \int_0^T g_{\mu\nu}(\lambda(t))\,\dot\lambda^\mu(t)\,\dot\lambda^\nu(t)\,dt$$3 is only available numerically (e.g., from mean-field or simulation data), geodesic distances and optimal protocols can be constructed via discrete gridding and fast-marching or level-set methods. In practice, for the antiferromagnetic Ising system:

• Discretize the $$A = \int_0^T g_{\mu\nu}(\lambda(t))\,\dot\lambda^\mu(t)\,\dot\lambda^\nu(t)\,dt$$4 plane.
• Solve the mean-field equations at each grid point for magnetizations.
• Compute $$A = \int_0^T g_{\mu\nu}(\lambda(t))\,\dot\lambda^\mu(t)\,\dot\lambda^\nu(t)\,dt$$5 from the linear-response integrals.
• Use the Fast Marching Method to obtain the minimal distance field $$A = \int_0^T g_{\mu\nu}(\lambda(t))\,\dot\lambda^\mu(t)\,\dot\lambda^\nu(t)\,dt$$6.
• Recover the geodesic by steepest descent from target to source in $$A = \int_0^T g_{\mu\nu}(\lambda(t))\,\dot\lambda^\mu(t)\,\dot\lambda^\nu(t)\,dt$$7.

Four regime types appear: (i) Both endpoints disordered: optimal path runs along infinite temperature direction (metric rank-deficient). (ii) Endpoints in different phases: cross the critical line rapidly, then minimize via disordered-phase directions. (iii) Both endpoints in ordered phase with direct shorter geodesic: stay in ordered manifold. (iv) Shortest route crosses the transition twice: excursions into disordered regions are “cheap” in metric length.

These computations confirm that, when $$A = \int_0^T g_{\mu\nu}(\lambda(t))\,\dot\lambda^\mu(t)\,\dot\lambda^\nu(t)\,dt$$8 is finite, the true minimum-dissipation protocol can exploit phase transitions even if start and end points share a phase (Basri et al., 1 Dec 2025).

5. Geometric Interpretation of Fluctuations, Stability, and Correlation Lengths

Thermodynamic geometry generalizes the equilibrium fluctuation theory by relating metric components to variances and covariances of thermodynamic variables; in the entropy representation: $$A = \int_0^T g_{\mu\nu}(\lambda(t))\,\dot\lambda^\mu(t)\,\dot\lambda^\nu(t)\,dt$$9 On the state manifold, local heat capacities correspond to $$g_{\mu\nu}(\lambda) = \beta \int_0^\infty \langle \delta X_\mu(t)\,\delta X_\nu(0)\rangle_{\mathrm{eq}}\,dt = \beta\,T_{\mu\nu}(\lambda)\,\partial_\mu\partial_\nu\ln Z(\lambda)$$0, with global stability requiring positive-definiteness of $$g_{\mu\nu}(\lambda) = \beta \int_0^\infty \langle \delta X_\mu(t)\,\delta X_\nu(0)\rangle_{\mathrm{eq}}\,dt = \beta\,T_{\mu\nu}(\lambda)\,\partial_\mu\partial_\nu\ln Z(\lambda)$$1 and $$g_{\mu\nu}(\lambda) = \beta \int_0^\infty \langle \delta X_\mu(t)\,\delta X_\nu(0)\rangle_{\mathrm{eq}}\,dt = \beta\,T_{\mu\nu}(\lambda)\,\partial_\mu\partial_\nu\ln Z(\lambda)$$2 (Bellucci et al., 2010).

The scalar curvature $$g_{\mu\nu}(\lambda) = \beta \int_0^\infty \langle \delta X_\mu(t)\,\delta X_\nu(0)\rangle_{\mathrm{eq}}\,dt = \beta\,T_{\mu\nu}(\lambda)\,\partial_\mu\partial_\nu\ln Z(\lambda)$$3 encodes the strength and nature of correlations:

• $$g_{\mu\nu}(\lambda) = \beta \int_0^\infty \langle \delta X_\mu(t)\,\delta X_\nu(0)\rangle_{\mathrm{eq}}\,dt = \beta\,T_{\mu\nu}(\lambda)\,\partial_\mu\partial_\nu\ln Z(\lambda)$$4 (2D): divergence of $$g_{\mu\nu}(\lambda) = \beta \int_0^\infty \langle \delta X_\mu(t)\,\delta X_\nu(0)\rangle_{\mathrm{eq}}\,dt = \beta\,T_{\mu\nu}(\lambda)\,\partial_\mu\partial_\nu\ln Z(\lambda)$$5 marks phase transitions.
• Sign of $$g_{\mu\nu}(\lambda) = \beta \int_0^\infty \langle \delta X_\mu(t)\,\delta X_\nu(0)\rangle_{\mathrm{eq}}\,dt = \beta\,T_{\mu\nu}(\lambda)\,\partial_\mu\partial_\nu\ln Z(\lambda)$$6 distinguishes attractive ($$g_{\mu\nu}(\lambda) = \beta \int_0^\infty \langle \delta X_\mu(t)\,\delta X_\nu(0)\rangle_{\mathrm{eq}}\,dt = \beta\,T_{\mu\nu}(\lambda)\,\partial_\mu\partial_\nu\ln Z(\lambda)$$7) vs. repulsive ($$g_{\mu\nu}(\lambda) = \beta \int_0^\infty \langle \delta X_\mu(t)\,\delta X_\nu(0)\rangle_{\mathrm{eq}}\,dt = \beta\,T_{\mu\nu}(\lambda)\,\partial_\mu\partial_\nu\ln Z(\lambda)$$8) regimes.

At spinodal or critical points, divergence of $$g_{\mu\nu}(\lambda) = \beta \int_0^\infty \langle \delta X_\mu(t)\,\delta X_\nu(0)\rangle_{\mathrm{eq}}\,dt = \beta\,T_{\mu\nu}(\lambda)\,\partial_\mu\partial_\nu\ln Z(\lambda)$$9 signals instability and onset of critical fluctuations.

6. Riemannian Structure, Geodesic Equations, and Broader Implications

The optimal control problem is reduced to solving the geodesic equations: $\beta = 1/T$0 where the Christoffel symbols are constructed from $\beta = 1/T$1. This formalism is applicable across classical fluids, spin systems, quantum phase transitions, and computational protocols in finite-time thermodynamics (Bellucci et al., 2010, Basri et al., 1 Dec 2025).

Thermodynamic geometry thereby provides:

• Quantitative criteria for stability and criticality via curvature and determinant conditions.
• Unification of fluctuation theory, optimal dissipation, and linear response.
• Connection between phase transition theory and transport-based optimal control.

Extensions to quantum systems utilize contact and fiber-bundle structures, where the divergence of geodesic length toward rank-deficient quantum states geometrizes the third law of thermodynamics (Tejero et al., 13 Nov 2025).

7. Summary Table: Key Features of Thermodynamic Geometry

Feature	Mathematical Formulation	Physical Interpretation
Metric tensor $\beta = 1/T$2	Equilibrium covariances, linear response	Fluctuations, response, dissipation
Thermodynamic length $\beta = 1/T$3	Integral along protocol geodesic	Minimal work, protocol optimality
Scalar curvature $\beta = 1/T$4	Determinant/Hessian of $\beta = 1/T$5	Correlation strength, phase transition
Geodesics	Euler–Lagrange eqs, Christoffel symbols	Minimum-dissipation protocol trajectories
Divergence at criticality	$\beta = 1/T$6, $\beta = 1/T$7 blow up	Breakdown of stability, fluctuation scaling

In total, thermodynamic geometry rigorously connects equilibrium fluctuation phenomena, optimal nonequilibrium control, and phase-transition theory in a unified, coordinate-invariant framework. Geometric criteria characterize when dissipation can be minimized across phase transitions and elucidate the role of critical exponents and universality. Numerical and analytical methods enable explicit protocol construction and stability diagnostics even in regimes with divergent fluctuation behavior (Basri et al., 1 Dec 2025, Bellucci et al., 2010).