Morphometric Solvation Free Energy

• The paper introduces a morphometric ansatz that expresses solvation free energy using a linear combination of intrinsic volumes such as volume, surface area, curvature, and Euler characteristic.
• It derives the approach through a controlled resummation of the virial series and integral-geometric techniques, providing a rigorous framework for hard-body fluids.
• The work connects morphometric thermodynamics with molecular and continuum solvation theories, elucidating many-body correlations and offering a geometric perspective on insertion free energies.

The morphometric approach treats solvation or insertion free energy as a geometric functional of solute shape. For convex bodies in three dimensions, the standard morphometric ansatz writes the free-energy change as a linear combination of four geometric measures of the solute, namely volume $V$, surface area $A$, integrated mean curvature $C$, and integrated Gaussian curvature $X$, with state-dependent coefficients. In hard-body fluids, this form can be derived from a controlled resummation of the virial series, and in broader solvation theories the same logic reappears through pressure, surface-tension, cavity, and partial-molar-volume terms (Robinson et al., 2019).

1. Geometric functional form

For a hard-particle fluid into which a hard-body solute is inserted, the solvation or insertion free energy is the free-energy change due to imposing the sharp external potential of the solute. In three dimensions, the standard morphometric ansatz is

$\Delta \Omega = p\,V + a_2\,A + a_1\,C + a_0\,X,$

where $p$ is the bulk pressure and $a_2$, $a_1$, and $a_0$ are thermodynamic coefficients depending on the bulk state of the fluid. Physically, $pV$ is the bulk work against pressure to carve out the excluded volume, $A$0 is a surface-tension-like term, $A$1 encodes curvature corrections, and $A$2 is a topological contribution tied to the Euler characteristic (Robinson et al., 2019).

The same structure has a natural formulation in arbitrary spatial dimension $A$3,

$A$4

where $A$5 are the intrinsic volumes of a convex, compact body $A$6. In this sense, morphometric thermodynamics is a finite-dimensional reduction of solvation free energy: for convex bodies, the dependence on shape is compressed into a small set of additive, motion-invariant geometric measures. Within the dilute limit, the insertion cost is identified with the excess chemical potential of the solute, so the morphometric form is simultaneously a geometric representation of solvation free energy and of $A$7 (Robinson et al., 2019).

2. Integral-geometric and virial foundations

A central advance in the modern theory is the derivation of the morphometric form from a controlled virial expansion. Robinson, Roth, and Royall derive the morphometric approach as the exact resummation of terms in the virial series by isolating contributions in which the solute and surrounding solvent particles share a common point of intersection. For convex hard bodies, the Mayer function can be identified with the Euler characteristic of the intersection, which turns the relevant cluster integrals into geometric intersection functionals (Robinson et al., 2019).

The derivation proceeds by writing the excess free-energy density as a virial series and then extracting the solute excess chemical potential from derivatives of the virial coefficients. The key approximation is to retain only the class of diagrams corresponding to a single common overlap region. Integral geometry then supplies the principal and multinomial kinematic formulas, which express integrals of Euler characteristics over rigid motions as sums of products of intrinsic volumes. After regrouping terms, the insertion free energy takes the exact morphometric form for that diagram class,

$A$8

The combinatorial coefficients are fixed by matching to the exact zero-dimensional cavity free energy, yielding $A$9. A low-density result follows directly from the kinematic formula: for convex bodies in $C$0, the leading term of the insertion free energy is already explicitly morphometric. This derivation explains why the ansatz is so effective: within the dominant intersection-diagram sector, solvation free energy is a Hadwiger-type functional and therefore must be linear in the intrinsic volumes (Robinson et al., 2019).

3. Accuracy, mixtures, and controlled limitations

The morphometric approach is not a purely formal rewriting; it gives explicit thermodynamic coefficients and can be compared against established liquid-state theories. In $C$1, the resummed theory is exact for hard rods. In $C$2, it coincides with scaled particle theory for hard disks. In $C$3, however, the resulting theory is less accurate than previous morphometric theories based on more refined density-functional constructions, even though it captures the dominant leading contribution (Robinson et al., 2019).

This limited accuracy is quantitative rather than qualitative. At the freezing packing fraction $C$4, the scaled-particle/Percus–Yevick equation of state overestimates the pressure by about $C$5, while the controlled resummation underestimates it by about $C$6, taking Carnahan–Starling as quasi-exact. The same pattern carries over to interfacial coefficients such as the planar wall surface tension. The result is a clear hierarchy: the morphometric base is physically grounded and often accurate, but higher precision requires corrections beyond the standard four-term form (Robinson et al., 2019).

A common misconception is that such corrections should be added as higher curvature moments. The controlled expansion argues against that move. Terms such as $C$7 diverge for surfaces with edges or vertices, whereas the intrinsic volumes remain finite and well defined for all convex, compact bodies. In this framework, higher curvature moments are therefore not admissible as systematic corrections. The allowed extensions are non-morphometric terms,

$C$8

with $C$9 arising from neglected diagrams such as ring contributions and taking genuinely nonlocal forms (Robinson et al., 2019).

The same formalism extends to mixtures of convex bodies of arbitrary shape. Mixture composition enters through scaled particle variables

$X$0

and the insertion free energy of a solute still has the form $X$1. On this view, species-specific geometric measures are unnecessary; a single set of intrinsic volumes for the solute suffices, while the mixture enters only through the coefficients (Robinson et al., 2019).

4. Many-body correlations as a solvation problem

The morphometric framework also supports a reinterpretation of many-body correlations as insertion or solvation free energies of local structures. For an $X$2-tuple of particles in a homogeneous hard-sphere fluid, the potential of mean force can be written as

$X$3

Here, $X$4 is the grand-potential cost of inserting the union of exclusion spheres associated with the tagged particles. The same morphometric form then yields $X$5 and, after integration over a motif domain, the absolute free energy of a local geometric motif (Robinson et al., 2018).

This construction has been used to calculate absolute free energies of local geometric motifs in excellent quantitative agreement with molecular dynamics simulations across the liquid and supercooled liquid regimes. It reveals a bimodality in the density library of states where five-fold symmetric structures appear lower in free energy than four-fold symmetric structures, and from a single reaction path it predicts a relaxation barrier which scales linearly in the compressibility factor. In effect, the morphometric solvation cost of the exclusion region becomes a local free-energy landscape for structural motifs (Robinson et al., 2018).

The broader significance is that cavity geometry is not only a property of an inserted external solute. It also organizes the thermodynamics of endogenous local structures in dense fluids. This suggests that morphometric solvation theory and many-body correlation theory are two aspects of the same geometric-statistical framework.

5. Relation to molecular and continuum solvation theories

Outside hard-body fluids, the morphometric viewpoint persists in more molecular theories of solvation. In molecular density functional theory for water, a coarse-grained bridge functional was designed to enforce liquid–vapor coexistence and the experimental surface tension. Those are exactly the macroscopic parameters that, in a morphometric description, multiply the geometric measures $X$6 and $X$7. With that bridge, small hard spheres show the expected volume-dominated scaling, while for large spheres $X$8 converges to the SPC/E water surface tension, recovering the surface-dominated regime expected from morphometric thermodynamics (Gageat et al., 2017).

A related connection appears in the ensemble treatment of MDFT and 3D-RISM hydration free energies. The rigorous conversion from the grand-canonical MDFT free energy to the experimental NPT quantity introduces a linear-in-volume correction through the partial molar volume,

$X$9

This is not presented as morphometric thermodynamics in the strict sense, but it is a volume term derived from microscopic liquid-state theory rather than fitted empiricism. It clarifies why partial-molar-volume corrections work and shows how one morphometric component can emerge rigorously from ensemble transformation (Sergiievskyi et al., 2014).

More generally, integral-equation and density-functional theories provide microscopic generators of the geometric information that morphometric models compress. In 3D-RISM, solvation free energy is an integral over spatially resolved correlation functions. In PC-SAFT-based classical DFT, the hard-sphere contribution is built from fundamental measure theory, which is explicitly identified as the microscopic implementation of morphometric thermodynamics for the cavity-formation part in non-polar chain fluids (Kovalenko, 2015, Eller et al., 2021).

6. Extensions and present boundaries

Recent work has extended the geometric intuition of morphometric solvation into settings that are not explicit morphometric theories. In a PC-SAFT density functional for non-polar solvents, the hard-sphere FMT term can, in principle, be approximated by a morphometric expression for large, smooth solutes, while dispersion and chain-connectivity terms modify the effective coefficients (Eller et al., 2021). In the “bubble method” for first-principles and machine-learning molecular dynamics, an artificial repulsive bubble creates a solvent-excluded region around a solute of arbitrary shape; the cavity-formation free energy of that bubble can in principle be related to geometric descriptors of the exclusion region, even though the paper does not explicitly use morphometric thermodynamics (Yu et al., 18 Apr 2026).

At the same time, the classical scope of the theory remains sharply defined. The strongest results concern convex hard bodies and hard-sphere-like solvents. For realistic aqueous solvation, electrostatics, hydrogen bonding, and nonlocal solvent structure introduce terms that are not exhausted by $\Delta \Omega = p\,V + a_2\,A + a_1\,C + a_0\,X,$0, $\Delta \Omega = p\,V + a_2\,A + a_1\,C + a_0\,X,$1, $\Delta \Omega = p\,V + a_2\,A + a_1\,C + a_0\,X,$2, and $\Delta \Omega = p\,V + a_2\,A + a_1\,C + a_0\,X,$3. A plausible implication is that the morphometric approach is best regarded as a controlled geometric backbone: exact or asymptotically exact for certain hard-body sectors, often dominant for cavity formation, and highly informative when embedded in richer molecular theories rather than treated as a complete description.

In that form, the morphometric approach remains a central framework for understanding solvation free energy. It explains why a finite set of geometric invariants can encode much of the insertion cost, provides a rigorous route from virial expansion to solvation thermodynamics, and supplies a common language linking hard-particle theory, local-structure free energies, molecular density functionals, and modern explicit-solvent free-energy methods (Robinson et al., 2019, Robinson et al., 2018).