pith. sign in

arxiv: 2606.12275 · v1 · pith:JAXKIULJnew · submitted 2026-06-10 · ❄️ cond-mat.soft

Approximate additivity in the solvent-mediated potential of mean force for ultrasoft particle systems

Joshua F. Robinson , Gary Yu , Patrick B. Warren This is my paper

Pith reviewed 2026-06-27 07:50 UTC · model grok-4.3

classification ❄️ cond-mat.soft
keywords solvent-mediated potential of mean forcehypernetted-chain closureOrnstein-Zernike equationsultrasoft particlesdissipative particle dynamicsAsakura-Oosawa potentialgeneralised excluded volumedepletion interaction
0
0 comments X

The pith

The solvent-mediated potential of mean force between solutes equals a convolution of their generalised excluded volume functions in the infinite-dilution limit.

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

The paper establishes that the solvent-mediated PMF extracted via the HNC closure of the Ornstein-Zernike equations takes the form of a convolution between solute-specific generalised excluded volume functions. This identity holds exactly in the infinite-dilution limit and recovers the classic Asakura-Oosawa depletion potential when the solvent consists of point particles and the solutes have hard cores. The same construction remains accurate for ultrasoft repulsive potentials of the kind used in dissipative particle dynamics, where it reveals that effective pair interactions depend on the overlap of soft excluded volumes and therefore on the choice of intramolecular bond lengths.

Core claim

In the infinite dilution limit, the solvent-mediated potential of mean force (PMF) between solutes, extracted from the hypernetted-chain (HNC) closure of the Ornstein-Zernike equations, can be expressed as a convolution between solute-specific generalised excluded volume functions. In the limit of a structureless solvent of point particles and hard core solutes, this recovers the exact Asakura-Oosawa depletion potential as the overlap between excluded volume spheres. The methodology recovers the solvent-mediated PMF with considerable accuracy for ultrasoft particle systems such as those in dissipative particle dynamics.

What carries the argument

Convolution of solute-specific generalised excluded volume functions obtained from the hypernetted-chain closure of the Ornstein-Zernike equations.

If this is right

  • The solvent-mediated PMF can be computed without solving the full many-body problem once the generalised excluded volume functions are known.
  • In DPD simulations the non-bonded repulsion parameters become sensitive to intramolecular bond lengths whenever those lengths fall inside the range of the soft potential.
  • The same convolution structure supplies an exact reference for hard-sphere depletion forces and an approximate but accurate one for soft particles.
  • Parametrisation of coarse-grained models must account for the overlap of soft excluded volumes rather than treating them as additive.

Where Pith is reading between the lines

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

  • Effective pair potentials derived this way could be used to pre-screen bond-length choices before running expensive DPD simulations.
  • The approach may extend to other integral-equation closures or to mixtures with multiple solute species.
  • It offers a route to analytic expressions for depletion forces in systems where direct simulation is costly.

Load-bearing premise

The hypernetted-chain closure remains accurate enough to extract the solvent-mediated PMF when the particles interact through ultrasoft repulsive potentials.

What would settle it

Direct numerical comparison of the convolution expression against the PMF obtained from explicit molecular dynamics of the full DPD system at low but finite solute concentration, for a chosen soft potential and bond length.

Figures

Figures reproduced from arXiv: 2606.12275 by Gary Yu, Joshua F. Robinson, Patrick B. Warren.

Figure 2
Figure 2. Figure 2: FIG. 2. Effective hard-sphere radius calculated from Eq. ( [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 1
Figure 1. Figure 1: FIG. 1. Generalised indicator functions extracted by the two methods [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Excess chemical potential of monomers inserted into a [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Solvent-mediated potential between a pair of identical solutes [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Partition coefficients of ultrasoft homodimers in a monoatomic [PITH_FULL_IMAGE:figures/full_fig_p007_5.png] view at source ↗
read the original abstract

In the infinite dilution limit, we show that the solvent-mediated potential of mean force (PMF) between solutes, extracted from the hypernetted-chain (HNC) closure of the Ornstein-Zernike equations, can expressed as a convolution between solute-specific generalised excluded volume functions. In the limit of a structureless solvent of point particles and hard core solutes, this recovers the exact Asakura-Oosawa depletion potential as the overlap between excluded volume spheres. The methodology can be deployed for ultrasoft particle systems such as those encountered in dissipative particle dynamics (DPD), where the solvent-mediated PMF can be recovered with considerable accuracy. These results confirm that in coarse-grained molecular DPD simulations the parametrisation of the non-bonded repulsions is sensitive to the assumed intramolecular bond lengths if they are smaller than the range of the DPD potential, due to the overlap of the soft excluded volume functions.

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.

Referee Report

1 major / 0 minor

Summary. In the infinite dilution limit, the solvent-mediated potential of mean force (PMF) between solutes extracted from the hypernetted-chain (HNC) closure of the Ornstein-Zernike equations can be expressed as a convolution between solute-specific generalised excluded volume functions. This recovers the exact Asakura-Oosawa depletion potential for hard-core solutes in a structureless point-particle solvent. The methodology is applied to ultrasoft particle systems such as those in dissipative particle dynamics (DPD), where the solvent-mediated PMF is recovered with considerable accuracy, implying that parametrisation of non-bonded repulsions in coarse-grained DPD simulations is sensitive to intramolecular bond lengths smaller than the DPD potential range.

Significance. If the accuracy for ultrasoft potentials holds, the result supplies an efficient analytical route to solvent-mediated PMFs in soft-matter systems and clarifies effective interactions in DPD coarse-graining. The algebraic derivation within the HNC framework and the exact recovery of the known Asakura-Oosawa limit are strengths.

major comments (1)
  1. [Abstract] Abstract: the claim that the solvent-mediated PMF 'can be recovered with considerable accuracy' for ultrasoft DPD potentials lacks any independent numerical benchmark (e.g., explicit MD or Monte Carlo inversion of g(r)) that would quantify the HNC closure error for the same soft potentials and state points; this is load-bearing for the applicability statement.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading and positive assessment of the work's significance. We address the single major comment below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that the solvent-mediated PMF 'can be recovered with considerable accuracy' for ultrasoft DPD potentials lacks any independent numerical benchmark (e.g., explicit MD or Monte Carlo inversion of g(r)) that would quantify the HNC closure error for the same soft potentials and state points; this is load-bearing for the applicability statement.

    Authors: The algebraic derivation establishes that the solvent-mediated PMF is exactly the stated convolution within the HNC closure; this is not an approximation. The phrase 'recovered with considerable accuracy' in the abstract refers to numerical verification that the convolution expression reproduces the PMF obtained by direct iterative solution of the HNC-OZ equations for the DPD potentials and state points examined. We did not perform, and therefore do not claim, an independent validation against explicit MD or Monte Carlo data that would quantify the intrinsic error of the HNC closure itself. We agree that such a benchmark would be required to make stronger statements about applicability to DPD coarse-graining. We will revise the abstract to remove any ambiguity and add a short paragraph in the discussion noting the known limitations of HNC for ultrasoft repulsions. revision: yes

Circularity Check

0 steps flagged

No circularity: HNC-OZ convolution identity is algebraic consequence of closure, independent of inputs

full rationale

The central result is obtained by direct substitution of the HNC closure c(r) = −βu(r) + h(r) − ln(1 + h(r)) into the Ornstein–Zernike equation at infinite dilution, yielding an exact convolution expression for the solvent-mediated PMF within that closure. This is a mathematical identity, not a fit or self-definition. The paper explicitly recovers the known Asakura–Oosawa depletion potential as a special case, confirming the derivation is self-contained and externally verifiable. No self-citations, fitted parameters renamed as predictions, or ansatz smuggling appear in the load-bearing steps. The claim of “considerable accuracy” for DPD is presented as a numerical observation, not as a derivation that reduces to itself.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The result rests on the HNC closure and introduces new functions whose definition is internal to the derivation.

axioms (1)
  • domain assumption Hypernetted-chain (HNC) closure to the Ornstein-Zernike equations accurately captures the relevant pair correlations for the solvent-mediated PMF
    Invoked to extract the PMF from the integral equations in the infinite-dilution limit.
invented entities (1)
  • solute-specific generalised excluded volume functions no independent evidence
    purpose: To enable the convolution representation of the solvent-mediated PMF for ultrasoft particles
    New functions introduced to generalise the hard-sphere excluded volume to soft interactions.

pith-pipeline@v0.9.1-grok · 5689 in / 1412 out tokens · 22697 ms · 2026-06-27T07:50:39.771808+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

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

  1. [1]

    and Del Regno, Annalaura and Anderson, Richard L

    Bray, David J. and Del Regno, Annalaura and Anderson, Richard L. , year = 2020, journal =

    2020

  2. [2]

    and Anderson, Richard L

    Seaton, Michael A. and Anderson, Richard L. and Metz, Sebastian and Smith, William , year = 2013, journal =

    2013

  3. [3]

    Morphological Thermodynamics of Fluids: Shape Dependence of Free Energies , author =. Phys. Rev. Lett. , volume =

  4. [4]

    Density-functional theory for structure and freezing of star polymer solutions , author =. J. Chem. Phys. , volume = 114, pages =

  5. [5]

    Density-functional theory for soft interactions by dimensional crossover , author =. Phys. Rev. E , volume = 60, number = 6, pages =

  6. [6]

    Fundamental measure theory for non-spherical hard particles: predicting liquid crystal properties from the particle shape , author =. J. Phys.: Condens. Matter , volume = 28, number = 24, pages = 244003, doi =

  7. [7]

    Europhys

    Fundamental mixed measure theory for non-spherical colloids , author =. Europhys. Lett. , volume = 109, pages = 26003, doi =

  8. [8]

    The fluid structures for soft-sphere potentials via the zero-separation theorems on molecular distribution functions , author =. J. Chem. Phys. , volume = 104, pages =

  9. [9]

    Role of repulsive forces in determining the equilibrium structure of simple liquids , author =. J. Chem. Phys. , volume = 54, pages = 5237, doi =

  10. [10]

    Introduction to geometric probability , author =

  11. [11]

    New and proper integral equations for site-site equilibrium correlations in molecular fluids , author =. Mol. Phys. , volume = 46, pages =

  12. [12]

    Statistical mechanics of rigid spheres , author =. J. Chem. Phys. , volume = 31, pages =

  13. [13]

    Scaled field particle theory of the structure and the thermodynamics of isotropic hard particle fluids , author =. J. Chem. Phys. , volume = 89, pages =

  14. [14]

    Free-energy model for the inhomogeneous hard-sphere fluid mixture and density-functional theory of freezing , author =. Phys. Rev. Lett. , volume = 63, pages =

  15. [15]

    Fundamental measure theory for hard-sphere mixtures: a review , author =. J. Phys.: Condens. Matter , volume = 22, pages =

  16. [16]

    On interaction between two bodies immersed in a solution of macromolecules , author =. J. Chem. Phys. , volume = 22, pages =

  17. [17]

    Interaction between particles suspended in solutions of macromolecules , author =. J. Polym. Sci. , volume = 33, pages =

  18. [18]

    Metadensity functional theory for classical fluids: extracting the pair potential , author =. Phys. Rev. Lett. , volume = 134, pages = 107301, doi =

  19. [19]

    Metadensity functional learning for classical fluids: Regularizing with pair correlations

    Metadensity functional learning for classical fluids: regularizing with pair correlations , author =. doi:10.48550/arXiv.2603.11973 , eprint =

    work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2603.11973

  20. [20]

    Comparison of integral equation theories of the liquid state , author =. Phys. Rev. E , volume = 110, pages =

  21. [21]

    Coarse-grained models for automated fragmentation and parametrization of molecular databases , author =. J. Chem. Inf. Model. , volume = 56, pages =

  22. [22]

    Morphometric approach to many-body correlations in hard spheres , author =. Phys. Rev. Lett. , volume = 122, pages =

  23. [23]

    Some topics in the theory of fluids , author =. J. Chem. Phys. , volume = 39, pages =

  24. [24]

    Potential-distribution theory and the statistical mechanics of fluids , author =. J. Phys. Chem. , volume = 86, pages =

  25. [25]

    Solvent-mediated interactions and solvation close to fluid--fluid phase separation: a density functional treatment , author =. J. Chem. Phys. , volume = 118, pages =

  26. [26]

    Dissipative particle dynamics parametrisation using infinite dilution activity coefficients: the impact of bonding , author =. Phys. Chem. Chem. Phys. , volume = 27, pages =

  27. [27]

    Dissipative particle dynamics: systematic parametrization using water-octanol partition coefficients , author =. J. Chem. Phys. , volume = 147, pages =

  28. [28]

    Dissipative particle dynamics: bridging the gap between atomistic and mesoscopic simulation , author =. J. Chem. Phys. , volume = 107, pages =

  29. [29]

    and McDonald, I

    Hansen, J.-P. and McDonald, I. R. , year = 2006, publisher =

    2006

  30. [30]

    Louis, A. A. and Bolhuis, P. G. and Hansen, J. P. , year = 2000, journal =. Mean-field fluid behavior of the

    2000