2026/03/17 by Michael K. Gilson, Tom Kurtzman · 1 voice
Physics and Astronomy · #Advanced Thermodynamics and Statistical Mechanics #Quantum Mechanics and Non-Hermitian Physics #Statistical Mechanics and Entropy
paper · pdf · doi:10.26434/chemrxiv.15000996/v1
openalex publication_date 2026/03/17 · openalex created_date 2026/03/18 · openalex updated_date 2026/07/14
Systems in which the free energy density is nonuniform in space are familiar: the surface tension of a water droplet and the surface energy of a solid are good examples. Some such cases can be treated with prior theory, notably inhomogeneous solvation theory (IST), but IST is applicable only to liquids. Despite this limitation, IST has proven useful as a guide to the design of ligands to bind a targeted protein, based on the idea that ligands which displace water at high free energy will tend to bind more tightly, \em ceteris paribus. Here, we present Generalized Thermodynamic Mapping Theory (GTMT), a more general theory that is applicable to the entirety of a biomolecular or other chemical system, and which thus may provide additional guidance for molecular design. For example, it might highlight parts of a ligand whose local free energy rises on binding, thus suggesting where modifications could improve affinity. Starting from classical statistical thermodynamics, we derive structural decompositions that assign energy and entropy to individual atoms or to larger chemical components such as amino acid residues, and spatial decompositions that define continuously varying thermodynamic densities throughout the system. The potential energy is decomposed via the multibody expansion, and the entropy via the mutual information expansion. The resulting thermodynamic densities satisfy key desiderata: their spatial integrals yield the correct total thermodynamic quantities, the densities vanish where the atomic number density is zero, and all entropy terms above first order go to zero in the absence of correlation. The entropy decomposition in GTMT is closely related to that of IST, but it is simpler, largely because it expresses the entropy in terms of normalized probability density functions instead of non-normalized correlation functions. We show that, while GTMT is formally exact for any number of particles, N, IST is formally inexact for finite N, although it becomes exact in the thermodynamic limit. The unified framework presented here enables thermodynamic mapping of solute and solvent alike and is expected to support a range of applications, including in structure-based drug design, protein design, the analysis of allostery, and materials science.