arxiv: 2604.25882 · v1 · submitted 2026-04-28 · ❄️ cond-mat.dis-nn

Recognition: unknown

Excluded volume and molecular field in the Lennard-Jones fluid: a modified first-order perturbation theory

A. Trokhymchuk , V. Hordiichuk , R. Melnyk , I. Nezbeda

Authors on Pith no claims yet

Pith reviewed 2026-05-07 13:41 UTC · model grok-4.3

classification ❄️ cond-mat.dis-nn

keywords Lennard-Jones fluidfirst-order perturbation theoryequation of statereference systemexcluded volumemolecular fieldsupercritical thermodynamicsperturbation expansion

0 comments

The pith

A range-based split of the Lennard-Jones potential lets first-order perturbation theory match simulation data when state derivatives are retained.

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

The paper examines a modified first-order perturbation theory for the Lennard-Jones fluid in which the reference system takes the entire short-range part of the interaction and the perturbation is limited to the long-range tail. This decomposition turns the perturbation into a small, smoothly varying, near-mean-field term over a wide supercritical domain. When the density and temperature derivatives of this term are kept consistently, the resulting equation of state reproduces high-accuracy reference data with excellent fidelity. The results indicate that the success of first-order perturbation theory depends primarily on the physical content of the reference system and on consistent treatment of its state dependence, rather than on the formal order of the expansion.

Core claim

In this modified first-order perturbation theory the reference system incorporates the full short-range part of the Lennard-Jones interaction while the perturbation is confined to the remaining long-range tail. The range-based decomposition transforms the perturbation contribution into a small, smoothly varying, near-mean-field quantity. Retaining its density and temperature derivatives produces an equation of state that reproduces high-accuracy reference data with excellent fidelity across a broad supercritical thermodynamic domain. The work establishes that the success of first-order perturbation theory is governed by the physical content of the reference system and by the consistent state

What carries the argument

The range-based decomposition that assigns the entire short-range interaction to the reference system and the long-range tail to the perturbation, together with consistent retention of its density and temperature derivatives.

Load-bearing premise

The range-based split of the potential turns the perturbation into a small, smoothly varying near-mean-field quantity over a broad supercritical domain without extra parameters or adjustments.

What would settle it

High-accuracy molecular simulations at multiple supercritical temperatures and densities that show large systematic deviations in pressure or internal energy from the predicted equation of state would falsify the claim of excellent fidelity.

Figures

Figures reproduced from arXiv: 2604.25882 by A. Trokhymchuk, I. Nezbeda, R. Melnyk, V. Hordiichuk.

**Figure 1.** Figure 1: Lennard-Jones potential u(r) (the thin solid line) and two choices of its decomposition. The left panel: The CA decomposition suggested within the WCA method – into the repulsive interaction (the thick solid line) and the attractive interaction (the thin dashed line). The right panel: The NCA decomposition suggested in Ref. [15] and used in present study – into the short-range interaction uSR(r) (the thick… view at source ↗

**Figure 2.** Figure 2: Vapor-liquid phase diagram of the system with short-range attraction uSR(r) according to Eq. (7) (the empty circles) from Ref. [15] in comparison against the same for the Lennard-Jones fluid (the solid triangles) from Ref. [16]. The empty small circles in supercritical region indicate thermodynamic states for which MD simulations were performed. Also shown in view at source ↗

**Figure 3.** Figure 3: Dependence of the perturbation term a(T, ρ) on density for the set of fixed temperatures T ⋆ shown in figure. The results of calculations according to Eq. (11) are marked by symbols. The thin solid curves correspond to the results of polynomial parametrization, Eq. (14), while the thick solid horizontal line presents the result of mean-field approximation by assuming gSR(r) = 1 in Eq. (11). The origin of t… view at source ↗

**Figure 4.** Figure 4: Dependence of the perturbation term a(T, ρ) on temperature for the set of fixed densities ρ ⋆ shown in figure. The results of calculations according to Eq. (11) are marked by symbols. The thin solid curves correspond to the results of polynomial parametrization, Eq. (14), while the thick solid horizontal line presents the result of mean-field approximation by assuming gSR(r) = 1 in Eq. (11). corresponding … view at source ↗

**Figure 5.** Figure 5: Dependence of the pressure βp = p/kT on density for the set of fixed temperatures T ⋆ = kT /ϵLJ. The solid lines are the results of of computer-based equation of state for LJ fluid [1], the filled symbols show results of Eq. (17). Nezbeda [1]. This equation of state, derived from theoretical considerations combined with extensive simulation data, is widely regarded as one of the most reliable benchmarks fo… view at source ↗

**Figure 6.** Figure 6: The same as in view at source ↗

**Figure 7.** Figure 7: Dependence of the internal energy βE/N = E/N kT on temperature for the set of fixed densities ρ ⋆ = ρσ3 LJ shown in figure. The solid lines are the results of computer-based approach due to Kolafa and Nezbeda [1], the filled symbols show results of Eq. (18). and strongly state-dependent, the present range-based decomposition transforms the perturbation into a weak and regular contribution. A central outcom… view at source ↗

read the original abstract

The equation of state and, more generally, the thermodynamics of the Lennard-Jones fluid have long served as a benchmark problem in the statistical theory of fluids. Among available theoretical approaches, first-order perturbation theory occupies a special position: only at this level does the correction to the Helmholtz free energy admit an exact statistical-mechanical expression. In this work, we present a systematic, simulation-based assessment of a non-classical first-order perturbation theory in which the reference system incorporates the entire short-range part of the interaction, while the perturbation is confined to the remaining long-range tail. We show that this range-based decomposition transforms the perturbation contribution into a small, smoothly varying, near-mean-field quantity over a broad supercritical thermodynamic domain. When its density and temperature derivatives are consistently retained, the resulting equation of state reproduces high-accuracy reference data with excellent fidelity. The results demonstrate that the success of first-order perturbation theory is governed primarily by the physical content of the reference system and by the consistent treatment of its state dependence, rather than by the formal truncation order of the expansion.

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

The paper refines first-order perturbation theory for the LJ fluid by splitting the potential at its minimum and retaining full state dependence in the reference, yielding close matches to simulation data without extra parameters.

read the letter

The main thing to know is that this approach gets solid results for the Lennard-Jones equation of state by letting the reference system take the full short-range part of the potential up to the minimum, then treating the long-range tail as a small perturbation whose density and temperature derivatives are kept when building the free energy and its derivatives. That choice turns the perturbation into something smooth and nearly mean-field over a wide supercritical range, and the numerics line up well with high-accuracy reference data. No extra fitting parameters are added along the way. What is actually new here is the systematic simulation check of this particular range-based split inside a strictly first-order framework, plus the explicit demonstration that keeping the state dependence matters more than pushing to higher orders. They compute the perturbation integral directly from reference runs and differentiate numerically, which avoids some of the usual approximations in older perturbation schemes. The paper does this cleanly and shows the improvement over more conventional splits. The comparisons look credible on the evidence given, and the construction avoids obvious circularity since the split is fixed by the potential shape rather than tuned to the target results. One minor soft spot is that the full details on how the short-range reference is simulated (cutoffs, system size, or convergence checks) would help for full reproducibility, though nothing in the reported logic suggests hidden inconsistencies. The central claim holds up without load-bearing flaws. This is for people who work on liquid-state theory or practical equations of state for simple fluids. Readers who care about perturbation methods or supercritical modeling will find the explicit results and the argument about reference content useful. It is worth sending to peer review because the work is grounded, the validation is independent, and the refinement is testable.

Referee Report

0 major / 3 minor

Summary. The manuscript develops a modified first-order perturbation theory for the Lennard-Jones fluid in which the potential is split at its minimum: the reference system incorporates the entire short-range portion (repulsive core plus the attractive well), while the perturbation is restricted to the long-range tail. The first-order correction to the Helmholtz free energy is computed directly from simulations of the reference system, and the equation of state is obtained by consistent numerical differentiation of the total free energy with respect to density and temperature. The resulting theory is shown to reproduce high-accuracy reference data with excellent fidelity over a broad supercritical domain without adjustable parameters.

Significance. If the reported agreement holds under detailed scrutiny, the work would establish that the physical content of the reference system and the retention of its full state dependence are the dominant factors controlling the accuracy of first-order perturbation theory, rather than the formal truncation order. This parameter-free, simulation-informed approach could simplify thermodynamic calculations for LJ-type fluids and provide a template for range-decomposed treatments of other potentials, bridging analytic theory and direct simulation.

minor comments (3)

The precise location of the range split (relative to the potential minimum) and its relation to the standard WCA division should be stated explicitly in the methods section, including any rationale for retaining the full short-range well in the reference.
Figures comparing the EOS to reference data would benefit from explicit uncertainty bands on the simulation-derived perturbation integrals and a statement of the thermodynamic range (densities and temperatures) over which the comparisons were performed.
A brief discussion of how the numerical differentiation for the pressure and internal energy is implemented (finite-difference step size, smoothing, or analytic alternatives) would improve reproducibility.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of our manuscript and for the accurate summary of its central results. The referee's comments correctly highlight the role of the range-based split and the consistent differentiation of the free energy. No specific major comments were raised in the report.

Circularity Check

0 steps flagged

No circularity; derivation uses independent simulations and external benchmarks

full rationale

The paper decomposes the LJ potential at its minimum (WCA-style reference plus long-range tail), evaluates the first-order perturbation integral directly from reference-system simulations, and obtains the EOS via numerical differentiation of the total free energy while retaining full state dependence. All reported fidelity is measured against independent high-accuracy simulation data; the range split is an external physical choice, not fitted to the target EOS. No load-bearing self-citation, no fitted parameter renamed as prediction, and no self-definitional reduction appear in the construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard statistical-mechanical framework for fluids plus the domain assumption that a clean range-based split of the Lennard-Jones potential is physically meaningful and produces a near-mean-field perturbation. No free parameters or invented entities are indicated in the abstract.

axioms (1)

domain assumption The Lennard-Jones interaction can be partitioned into a short-range reference part and a long-range perturbation tail on the basis of distance.
This partition is the defining step of the modified theory and is invoked to make the perturbation small and mean-field.

pith-pipeline@v0.9.0 · 5501 in / 1337 out tokens · 56613 ms · 2026-05-07T13:41:31.952183+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

16 extracted references · 4 canonical work pages

[1]

Kolafa, I

J. Kolafa, I. Nezbeda, Fluid Phase Equil.100, 1-34 (1994) DOI: 10.1016/0378-3812(94)80 0 01-4

work page doi:10.1016/0378-3812(94)80 1994
[2]

Stephan, J

S. Stephan, J. Staubach, H. Hasse, Fluid Phase Equil.523, 112772 (2020) DOI: 10.1016/j.fluid.2020.112772

work page doi:10.1016/j.fluid.2020.112772 2020
[3]

Barker, D

J.A. Barker, D. Henderson, Rev. Mod. Phys.,48, 587 (1976)

1976
[4]

Weeks, D

J.D. Weeks, D. Chandler, H. Andersen, Science220, 787 (1983)

1983
[5]

Hansen, I.R

J.-P. Hansen, I.R. McDonald,Theory of Simple Liquids, Academic Press, London (1986)

1986
[6]

Hansen, I.R

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

2013
[7]

van der Waals,On the continuity of the gaseous and liquid states, Ph.D

J.D. van der Waals,On the continuity of the gaseous and liquid states, Ph.D. thesis, Leiden University (1873)
[8]

Barker, D

J.A. Barker, D. Henderson, J. Chem. Phys.,47, 4714 (1967)

1967
[9]

Andersen, J.D

H.C. Andersen, J.D. Weeks, D. Chandler, Phys. Rev. A4, 1597 (1971)

1971
[10]

Zwanzig, J

R.W. Zwanzig, J. Chem. Phys.22, 1420 (1954)

1954
[11]

van Westen, J

T. van Westen, J. Gross, J. Chem. Phys.147, 014503 (2017)

2017
[12]

Melnyk, I

R. Melnyk, I. Nezbeda, D. Henderson, A. Trokhymchuk, Fluid Phase Equil.279, 1 (2009)

2009
[13]

Nezbeda, R

I. Nezbeda, R. Melnyk, A. Trokhymchuk, J. Supercritical Fluids55, 448 (2010)

2010
[14]

Toxvaerd, J.C

S. Toxvaerd, J.C. Dyre, J. Chem. Phys.134, 081102 (2011)

2011
[15]

Hordiichuk, J

V. Hordiichuk, J. Skvara, A. Trokhymchuk, I. Nezbeda, J. Mol. Liquids386, 122483 (2023) DOI: 10.1016/j.molliq.2023.122483

work page doi:10.1016/j.molliq.2023.122483 2023
[16]

Trokhymchuk, J

A. Trokhymchuk, J. Alejandre, J. Chem. Phys.111, 8510 (1999) DOI: 10.1063/1.480192

work page doi:10.1063/1.480192 1999