Perturbation Theory of the Free Energy via the Mesoscopic Combined Partition Function

Bob Osano

arxiv: 2605.18121 · v1 · pith:BL2J7RL2new · submitted 2026-05-18 · ❄️ cond-mat.stat-mech

Perturbation Theory of the Free Energy via the Mesoscopic Combined Partition Function

Bob Osano This is my paper

Pith reviewed 2026-05-20 00:50 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech

keywords free energyperturbation theorymesoscopic partition functionmutual informationcoarse-grainingclassical statistical mechanicsextensivitymany-body systems

0 comments

The pith

The free energy of a classical N-body system equals a mesoscopic free energy minus inter-cell mutual information corrections plus a small error.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a systematic perturbation theory for the Helmholtz free energy of a classical many-body system by introducing a mesoscopic partition function obtained through coarse-graining of single-particle phase space into product cells. The reference mesoscopic partition function factorizes exactly via the multinomial theorem, which enables a controlled expansion in the inter-cell interaction terms. This yields an exact relation expressing the full free energy in terms of the mesoscopic free energy and corrections from pairwise mutual informations between cells, along with a remainder that vanishes exponentially with cell size relative to the correlation length. The first-order result recovers the van der Waals equation, while higher orders connect to known structure-factor expressions in appropriate limits. For long-range forces the same mutual-information terms directly measure the breakdown of extensivity.

Core claim

The full free energy satisfies F(λ)=ℱ_meso(λ)−k_BT ∑_{i<j} I(i,j;λ)+O(|Λ|ℓ^{-d}e^{-2ℓ/ξ}), where the inter-cell mutual informations I(i,j;λ) are the corrections identified in the extensivity analysis. This equality follows from perturbation theory applied to the mesoscopic partition function ℤ_meso(λ) that is generated by the combined coarse-graining operator acting on single-particle phase space and whose reference level factorizes exactly as (Z_1^{(0)})^N.

What carries the argument

The combined coarse-graining operator C=C_x∘C_p that partitions single-particle phase space into product cells C_{i,α}=V_i×Π_α and thereby produces a mesoscopic partition function whose unperturbed level factorizes exactly by the multinomial theorem.

If this is right

The first-order perturbation recovers the van der Waals equation and the Barker-Henderson result.
The second-order term converges to the structure-factor formula for the free energy in the fine-cell limit.
The perturbation expansion produces the mesoscopic Gibbs-Bogoliubov inequality.
An exact coupling-parameter integration formula for the free energy is obtained directly from the mesoscopic partition function.
For long-range interactions the mutual-information corrections quantify the resulting non-extensivity.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Choosing cell size ℓ ≫ ξ makes the exponential error negligible, suggesting a practical route to accurate free-energy estimates from modest simulations.
The mutual-information corrections supply an information-theoretic diagnostic for extensivity violations that could be computed in Monte Carlo sampling.
The same factorization-plus-correction structure may generalize to time-dependent or driven systems by promoting the coarse-graining operator to a dynamical map.

Load-bearing premise

The reference mesoscopic partition function factorizes exactly as Z_meso^{(0)}=(Z_1^{(0)})^N by the multinomial theorem under the combined coarse-graining operator.

What would settle it

Compute the exact free energy for a finite system with a short-range potential at a given cell size ℓ much larger than the correlation length ξ, subtract the mesoscopic free energy and the sum of inter-cell mutual informations, and check whether the residual difference lies inside the stated O(|Λ|ℓ^{-d}e^{-2ℓ/ξ}) bound.

read the original abstract

We develop a systematic perturbation theory for the Helmholtz free energy of a classical $N$-body system within the mesoscopic framework of~\cite{OsanoMeso,OsanoExtensivity}. The combined coarse-graining operator $\mathcal{C}=\mathcal{C}_x\circ\mathcal{C}_p$ acting on single-particle phase space partitions it into product cells $C_{i,\alpha}=V_i\times\Pi_\alpha$ and generates a mesoscopic partition function $\mathcal{Z}_{\rm meso}(\lambda)$ whose reference level factorises by the multinomial theorem: $\mathcal{Z}_{\rm meso}^{(0)}=(Z_1^{(0)})^N$. Perturbation theory for $\mathcal{F}_{\rm meso}(\lambda)=-k_BT\ln\mathcal{Z}_{\rm meso}(\lambda)$ in the inter-cell perturbation $\mathcal{V}_{\rm meso}$ yields the mesoscopic Gibbs--Bogoliubov inequality and an exact coupling-parameter integration formula. The full free energy satisfies \begin{equation*} F(\lambda)=\mathcal{F}_{\rm meso}(\lambda)-k_BT\!\sum_{i<j}I(i,j;\lambda)+O\!\left(|\Lambda|\ell^{-d}e^{-2\ell/\xi}\right), \end{equation*} where the inter-cell mutual informations $I(i,j;\lambda)$ are the corrections identified in the extensivity analysis. The first-order theory recovers the van der Waals equation and the Barker--Henderson result; the second-order term converges to the structure-factor formula in the fine-cell limit. For long-range interactions, factorisation fails, and the mutual-information corrections quantify the resulting non-extensivity.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This extends the mesoscopic approach to a free-energy perturbation series that recovers known limits, but the exact factorization of the reference partition function needs explicit verification.

read the letter

The main point here is a perturbation expansion for the free energy built on the mesoscopic partition function. The claimed relation is F(λ) equals the mesoscopic free energy minus kT times the sum of inter-cell mutual informations, plus a small error that decays exponentially with cell size. What is new is applying the coupling-parameter method inside this coarse-grained picture to generate the expansion and to identify the mutual-information terms as the non-extensivity corrections. It recovers the van der Waals equation at first order and the structure-factor result at second order in the appropriate limit. For long-range interactions the same terms quantify how extensivity breaks down. The construction looks consistent on its own terms and gives a systematic way to add those corrections without starting from scratch. The soft spot is the asserted exact factorization of the reference mesoscopic partition function as (Z_1^(0))^N via the multinomial theorem. The stress-test note raises a fair question about whether the 1/N! factor for indistinguishable particles is properly absorbed by the cell definitions. If it is only approximate, the downstream steps that rely on exact factorization could pick up corrections. The abstract does not walk through the algebra, so this needs to be shown clearly in the full paper. The heavy dependence on the two prior self-citations for the mesoscopic setup and the mutual informations is another limitation; the current work stands or falls with those earlier claims. This paper is for statistical mechanicians interested in mesoscopic methods or corrections to extensivity in classical systems. A reader who has followed the author's previous papers on the topic will find the perturbation construction useful. It is not aimed at a broad audience. The work shows honest engagement with the problem of non-extensivity and organizes the expansion logically. It deserves a serious referee to examine the derivations and test the factorization assumption. I recommend sending it to peer review.

Referee Report

3 major / 2 minor

Summary. The manuscript develops a perturbation theory for the Helmholtz free energy of a classical N-body system using a mesoscopic combined partition function generated by the coarse-graining operator C = C_x ∘ C_p. The central claim is that the full free energy satisfies F(λ) = ℱ_meso(λ) − k_B T ∑_{i<j} I(i,j;λ) + O(|Λ| ℓ^{-d} e^{-2ℓ/ξ}), where the inter-cell mutual informations I(i,j;λ) correct for non-extensivity. The reference mesoscopic partition function is asserted to factorize exactly as Z_meso^{(0)} = (Z_1^{(0)})^N by the multinomial theorem, enabling a mesoscopic Gibbs–Bogoliubov inequality and exact coupling-parameter integration. First-order results recover the van der Waals equation and Barker–Henderson perturbation; second order converges to the structure-factor formula in the fine-cell limit.

Significance. If the central relation and its supporting derivations hold, the work supplies a systematic route from mesoscopic coarse-graining to the full free energy while quantifying extensivity violations through mutual information. Recovery of standard perturbative limits is a concrete strength. The approach may prove useful for long-range interactions where conventional extensivity assumptions break down, provided the foundational factorization and error bounds are placed on a rigorous footing.

major comments (3)

[Abstract] Abstract, paragraph beginning 'The combined coarse-graining operator': the exact factorization Z_meso^{(0)} = (Z_1^{(0)})^N is stated to follow from the multinomial theorem once single-particle phase space is partitioned into product cells C_{i,α} = V_i × Π_α. For indistinguishable classical particles the full partition function carries an explicit 1/N! factor; it is not shown whether the mesoscopic sum over occupation numbers n_{iα} with multinomial coefficients N!/∏n_{iα}! exactly cancels this factor or leaves a residual Stirling correction. Because the subsequent Gibbs–Bogoliubov inequality, coupling-parameter integration, and identification of the mutual-information corrections all presuppose an exactly factorized reference, even a small violation would propagate directly into the claimed relation for F(λ).
[Central equation] Central equation (the displayed formula for F(λ)): the error term O(|Λ| ℓ^{-d} e^{-2ℓ/ξ}) is presented without an explicit derivation or bound on the correlation length ξ. It is unclear whether this estimate follows from the same mesoscopic framework or is imported from the self-cited works; a self-contained sketch of how the exponential decay arises from the coarse-graining length ℓ would be required to assess its validity.
[Perturbation theory and recovery statements] Recovery of known results (van der Waals, Barker–Henderson, structure-factor formula): the manuscript asserts that first- and second-order terms reproduce these limits, yet supplies no intermediate derivation steps, explicit expansion of the mesoscopic perturbation V_meso, or comparison of coefficients. Without these steps it is impossible to verify that the mapping is direct rather than formal.

minor comments (2)

[Notation] Notation: the distinction between script F_meso and ordinary F is introduced only in the abstract; repeating the definition at the start of the main text would improve readability.
[Introduction / references] Self-containedness: the mesoscopic partition function and mutual-information corrections are defined via the self-cited papers OsanoMeso and OsanoExtensivity. A short appendix recapitulating the essential definitions would make the present manuscript more independent.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address each major point below with clarifications and indicate revisions that will be incorporated to strengthen the presentation and derivations.

read point-by-point responses

Referee: [Abstract] Abstract, paragraph beginning 'The combined coarse-graining operator': the exact factorization Z_meso^{(0)} = (Z_1^{(0)})^N is stated to follow from the multinomial theorem once single-particle phase space is partitioned into product cells C_{i,α} = V_i × Π_α. For indistinguishable classical particles the full partition function carries an explicit 1/N! factor; it is not shown whether the mesoscopic sum over occupation numbers n_{iα} with multinomial coefficients N!/∏n_{iα}! exactly cancels this factor or leaves a residual Stirling correction. Because the subsequent Gibbs–Bogoliubov inequality, coupling-parameter integration, and identification of the mutual-information corrections all presuppose an exactly factorized reference, even a small violation would propagate directly into the claimed relation for F(λ).

Authors: We thank the referee for highlighting this important point on indistinguishability. In the mesoscopic construction, the reference partition function sums over occupation-number distributions {n_{iα}} with each term weighted by the multinomial coefficient N! / ∏ n_{iα}! times the product of single-cell factors. This multinomial prefactor exactly cancels the 1/N! from the classical phase-space measure, yielding Z_meso^{(0)} = (Z_1^{(0)})^N with no residual Stirling correction in the thermodynamic limit. The cancellation is implicit in the framework of the cited works but we will add an explicit one-paragraph derivation in the revised manuscript to make the factorization self-contained. revision: yes
Referee: [Central equation] Central equation (the displayed formula for F(λ)): the error term O(|Λ| ℓ^{-d} e^{-2ℓ/ξ}) is presented without an explicit derivation or bound on the correlation length ξ. It is unclear whether this estimate follows from the same mesoscopic framework or is imported from the self-cited works; a self-contained sketch of how the exponential decay arises from the coarse-graining length ℓ would be required to assess its validity.

Authors: The error bound follows directly from the mesoscopic coarse-graining: the probability that particles in cells separated by more than ℓ remain correlated decays as e^{-ℓ/ξ} (with ξ the correlation length of the underlying microscopic system). Summing the pairwise contributions over the volume |Λ| at cell density ℓ^{-d} produces the stated O(|Λ| ℓ^{-d} e^{-2ℓ/ξ}) remainder. While the full derivation appears in our prior extensivity analysis, we will insert a concise, self-contained sketch immediately after the central equation in the revised manuscript, showing the exponential suppression step by step from the definition of the coarse-graining operator. revision: yes
Referee: [Perturbation theory and recovery statements] Recovery of known results (van der Waals, Barker–Henderson, structure-factor formula): the manuscript asserts that first- and second-order terms reproduce these limits, yet supplies no intermediate derivation steps, explicit expansion of the mesoscopic perturbation V_meso, or comparison of coefficients. Without these steps it is impossible to verify that the mapping is direct rather than formal.

Authors: We agree that the intermediate steps were omitted for brevity and that their inclusion will improve verifiability. In the revised manuscript we will expand the perturbation-theory section to provide: (i) the explicit first-order average of V_meso and its reduction to the van der Waals mean-field term for uniform density; (ii) the coupling-parameter integration that recovers the Barker–Henderson perturbation; and (iii) the second-order term in the fine-cell limit ℓ → 0, demonstrating convergence to the structure-factor formula via Fourier transform of the pair-correlation function, with direct coefficient comparison. These additions will occupy approximately one additional page. revision: yes

Circularity Check

3 steps flagged

Free-energy relation imports mutual-information corrections and mesoscopic factorization via self-citations

specific steps

self citation load bearing [Abstract, first sentence]
"We develop a systematic perturbation theory for the Helmholtz free energy of a classical N-body system within the mesoscopic framework of~[OsanoMeso,OsanoExtensivity]."

The entire perturbation construction for F_meso and the subsequent relation to the full F(λ) is declared to occur inside this framework; the framework's definitions of Z_meso, the reference factorization, and the cell structure are therefore presupposed rather than re-derived.
self citation load bearing [Abstract, equation and following clause]
"The full free energy satisfies F(λ)=F_meso(λ)−k_BT ∑_{i<j}I(i,j;λ)+O(|Λ|ℓ^{-d}e^{-2ℓ/ξ}), where the inter-cell mutual informations I(i,j;λ) are the corrections identified in the extensivity analysis."

The mutual-information corrections I(i,j;λ) that appear in the central identity are not derived in the present manuscript; they are imported by reference to the extensivity analysis (OsanoExtensivity), so the claimed relation between F and F_meso reduces to a prior self-cited result.
self definitional [Abstract, sentence on reference level]
"generates a mesoscopic partition function Z_meso(λ) whose reference level factorises by the multinomial theorem: Z_meso^{(0)}=(Z_1^{(0)})^N"

The exact factorization is asserted to follow directly from the multinomial theorem once the combined operator C=C_x∘C_p is applied; this factorization is the starting point for all subsequent Gibbs-Bogoliubov and coupling-parameter steps, yet the operator and its action on phase space are defined inside the same self-cited mesoscopic framework.

full rationale

The derivation chain begins by placing the perturbation theory inside the mesoscopic framework of the author's prior papers, asserts exact factorization of the reference partition function by multinomial theorem under the combined coarse-graining operator, and then states the central relation for F(λ) by subtracting the inter-cell mutual informations that are explicitly identified as coming from the extensivity analysis. These steps reduce the claimed result to quantities and assumptions introduced in the self-cited works rather than deriving them independently here. The perturbation steps themselves (Gibbs-Bogoliubov, coupling integration) are not shown to be circular, but the load-bearing reference level and correction terms are.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claim rests on the mesoscopic framework and extensivity analysis developed in the author's self-cited prior works together with the assumption that the coarse-graining operator produces non-overlapping product cells whose inter-cell interactions admit a controlled perturbation expansion.

axioms (1)

domain assumption The combined coarse-graining operator C = C_x ∘ C_p partitions single-particle phase space into product cells C_{i,α} = V_i × Π_α such that the reference partition function factorises exactly by the multinomial theorem.
Invoked immediately after the definition of Z_meso(λ) to obtain Z_meso^{(0)} = (Z_1^{(0)})^N.

invented entities (1)

Mesoscopic partition function Z_meso(λ) no independent evidence
purpose: To enable exact factorisation of the reference state and a subsequent perturbation expansion in inter-cell interactions.
Introduced as the central object of the mesoscopic framework; no independent falsifiable prediction is supplied in the abstract.

pith-pipeline@v0.9.0 · 5825 in / 1624 out tokens · 80256 ms · 2026-05-20T00:50:29.884784+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages · 1 internal anchor

[1]

The Mesoscopic Partition Function:A Combined Spatial and Phase-Space Cell Structure

B. Osano, The Mesoscopic Partition Function: A combined spatial and Phase–Space Cell Structure.PrepreintarXiv:2605.00958

work page internal anchor Pith review Pith/arXiv arXiv
[2]

Osano, Entropy Additivity from exponential decay of correlations: a coarse-grained operator approachsubmit/7606660.(2026)

B. Osano, Entropy Additivity from exponential decay of correlations: a coarse-grained operator approachsubmit/7606660.(2026)

work page arXiv 2026
[3]

R. W. Zwanzig, High-temperature equation of state by a perturbation method. I. Nonpolar gases,J. Chem. Phys.22, 1420–1426 (1954)

work page 1954
[4]

J. A. Barker and D. Henderson, Perturbation theory and equation of state for fluids. II. A successful theory of liquids,J. Chem. Phys.47, 4714–4721 (1967). 17

work page 1967
[5]

J. G. Kirkwood, Statistical mechanics of fluid mixtures,J. Chem. Phys.3, 300–313 (1935)

work page 1935
[6]

J. E. Mayer and M. G. Mayer,Statistical Mechanics(Wiley, New York, 1940)

work page 1940
[7]

Ruelle,Statistical Mechanics: Rigorous Results(W

D. Ruelle,Statistical Mechanics: Rigorous Results(W. A. Benjamin, New York, 1969)

work page 1969
[8]

J. D. Weeks, D. Chandler, and H. C. Andersen, Role of repulsive forces in determining the equilibrium structure of simple liquids,J. Chem. Phys.54, 5237–5247 (1971)

work page 1971
[9]

N. F. Carnahan and K. E. Starling, Equation of state for nonattracting rigid spheres, J. Chem. Phys.51, 635–636 (1969)

work page 1969
[10]

Hansen and I

J.-P. Hansen and I. R. McDonald,Theory of Simple Liquids, 4th ed. (Academic Press, Oxford, 2013)

work page 2013
[11]

Peierls, On a minimum property of the free energy,Phys

R. Peierls, On a minimum property of the free energy,Phys. Rev.54, 918–919 (1938)

work page 1938
[12]

Frenkel and B

D. Frenkel and B. Smit,Understanding Molecular Simulation: From Algorithms to Applications, 2nd ed. (Academic Press, San Diego, 2002)

work page 2002
[13]

J. L. Lebowitz and O. Penrose, Rigorous treatment of the van der Waals–Maxwell theory of the liquid-vapour transition,J. Math. Phys.7, 98–113 (1966)

work page 1966
[14]

Georgii,Gibbs Measures and Phase Transitions(de Gruyter, Berlin, 1988)

H.-O. Georgii,Gibbs Measures and Phase Transitions(de Gruyter, Berlin, 1988)

work page 1988
[15]

Buchert, On average properties of inhomogeneous fluids in general relativity: dust cosmologies,Gen

T. Buchert, On average properties of inhomogeneous fluids in general relativity: dust cosmologies,Gen. Rel. Grav.32, 105–125 (2000)

work page 2000
[16]

Osano, Second-order perturbation theory: a covariant approach involving a barotropic equation of stateGen

B. Osano, Second-order perturbation theory: a covariant approach involving a barotropic equation of stateGen. Rel. Grav.34(12), 125004 (2017) 18

work page 2017

[1] [1]

The Mesoscopic Partition Function:A Combined Spatial and Phase-Space Cell Structure

B. Osano, The Mesoscopic Partition Function: A combined spatial and Phase–Space Cell Structure.PrepreintarXiv:2605.00958

work page internal anchor Pith review Pith/arXiv arXiv

[2] [2]

Osano, Entropy Additivity from exponential decay of correlations: a coarse-grained operator approachsubmit/7606660.(2026)

B. Osano, Entropy Additivity from exponential decay of correlations: a coarse-grained operator approachsubmit/7606660.(2026)

work page arXiv 2026

[3] [3]

R. W. Zwanzig, High-temperature equation of state by a perturbation method. I. Nonpolar gases,J. Chem. Phys.22, 1420–1426 (1954)

work page 1954

[4] [4]

J. A. Barker and D. Henderson, Perturbation theory and equation of state for fluids. II. A successful theory of liquids,J. Chem. Phys.47, 4714–4721 (1967). 17

work page 1967

[5] [5]

J. G. Kirkwood, Statistical mechanics of fluid mixtures,J. Chem. Phys.3, 300–313 (1935)

work page 1935

[6] [6]

J. E. Mayer and M. G. Mayer,Statistical Mechanics(Wiley, New York, 1940)

work page 1940

[7] [7]

Ruelle,Statistical Mechanics: Rigorous Results(W

D. Ruelle,Statistical Mechanics: Rigorous Results(W. A. Benjamin, New York, 1969)

work page 1969

[8] [8]

J. D. Weeks, D. Chandler, and H. C. Andersen, Role of repulsive forces in determining the equilibrium structure of simple liquids,J. Chem. Phys.54, 5237–5247 (1971)

work page 1971

[9] [9]

N. F. Carnahan and K. E. Starling, Equation of state for nonattracting rigid spheres, J. Chem. Phys.51, 635–636 (1969)

work page 1969

[10] [10]

Hansen and I

J.-P. Hansen and I. R. McDonald,Theory of Simple Liquids, 4th ed. (Academic Press, Oxford, 2013)

work page 2013

[11] [11]

Peierls, On a minimum property of the free energy,Phys

R. Peierls, On a minimum property of the free energy,Phys. Rev.54, 918–919 (1938)

work page 1938

[12] [12]

Frenkel and B

D. Frenkel and B. Smit,Understanding Molecular Simulation: From Algorithms to Applications, 2nd ed. (Academic Press, San Diego, 2002)

work page 2002

[13] [13]

J. L. Lebowitz and O. Penrose, Rigorous treatment of the van der Waals–Maxwell theory of the liquid-vapour transition,J. Math. Phys.7, 98–113 (1966)

work page 1966

[14] [14]

Georgii,Gibbs Measures and Phase Transitions(de Gruyter, Berlin, 1988)

H.-O. Georgii,Gibbs Measures and Phase Transitions(de Gruyter, Berlin, 1988)

work page 1988

[15] [15]

Buchert, On average properties of inhomogeneous fluids in general relativity: dust cosmologies,Gen

T. Buchert, On average properties of inhomogeneous fluids in general relativity: dust cosmologies,Gen. Rel. Grav.32, 105–125 (2000)

work page 2000

[16] [16]

Osano, Second-order perturbation theory: a covariant approach involving a barotropic equation of stateGen

B. Osano, Second-order perturbation theory: a covariant approach involving a barotropic equation of stateGen. Rel. Grav.34(12), 125004 (2017) 18

work page 2017