Pith. sign in

REVIEW 4 major objections 6 minor 1 cited by

Quantifying many-body contributions to depletion forces

T0 review · 4 major / 6 minor · reviewed 2026-08-15 · deepseek-v4-flash

Pith's one-line read Below a size ratio of about 0.15, the depletion attraction between large colloids is unaffected by how many large colloids are present, confirming a purely geometric prediction by direct simulation of quasi-hard-sphere mixtures.

desk verdict Useful new data on density-dependent depletion forces for q=0.45 and 0.60, but the q=0.15 null result is overinterpreted as confirmation of the hard-sphere geometric threshold when the simulated particles are visibly soft. read the letter →

arxiv 2509.03342 v1 pith:FIXHF3GK submitted 2025-09-03 cond-mat.soft

classification cond-mat.soft
keywords depletionforcesmany-bodycontributionseffectivepairpotentialsbinarycolloidmixturessemi-grandcanonicalensemblecontractionofbaresizeratiothresholdintegralequationtheory
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper establishes when many-body effects can be ignored in depletion interactions, the entropic attraction between large colloidal particles immersed in a sea of much smaller ones. Using a force-contraction simulation scheme that can separate the contributions of pairs, triples, and larger groups, the authors show that for size ratios $q = \sigma_s/\sigma_l \le 0.15$ the density of large particles has no effect on their mutual depletion force, provided the chemical potential of the small depletants is held fixed; the force curves for large-particle packing fractions $\phi_l = 0.1$ to $0.4$ superpose within numerical error. This confirms a decades-old geometric argument: below $q_3 = 2/\sqrt{3} - 1 \approx 0.1547$ a small sphere can slip through the gap between three touching large spheres, so it can never couple three large particles at once. For less asymmetric mixtures ($q = 0.45$ and $0.60$) the large-particle density strongly modulates the force, with roughly a 20% spread at $q = 0.6$, and that density dependence is captured by low-degree polynomials. If correct, the work gives a directly testable criterion for when pair-additive depletion potentials are good enough, plus a method to quantify how much of any effective force comes from three or more particles.

What carries the argument

The central mechanism is a geometric threshold for hard spheres: the size ratio $q_3 = 2/\sqrt{3} - 1 \approx 0.1547$ at which a small sphere just fits through the triangular cavity formed by three mutually touching large spheres; below it a depletant can never be in triple contact with large particles, so three-body and higher effective interactions are geometrically forbidden, and a second ratio $q_4 = \sqrt{3/2} - 1 \approx 0.2247$ marks where four-body terms switch on. The simulation engine is the contraction of bare forces (CBF), a least-squares inversion that extracts effective pair forces $G(r)$ from the noisy instantaneous forces on large particles averaged over configurations of the small ones; by resampling only configurations free of large-particle triangles, the same engine produces a many-body-free version of the force. The semi-grand ensemble is the enabling thermodynamic frame: the small-particle chemical potential $\mu_s$ is measured by a particle-insertion method in a reservoir and held constant while $\phi_l$ varies, with the small-particle packing fraction $\phi_s$ adjusted accordingly, which is what makes the predicted independence from $\phi_l$ a well-posed, testable statement. The Ornstein-Zernike integral-equation approach with the modified-Verlet closure supplies an independent, much cheaper route to the same effective forces.

What would settle it

Repeat the fixed-$\mu_s$ protocol at $q = 0.15$ with systematically steeper and softer repulsive potentials, or with an event-driven true-hard-sphere algorithm, and test whether the superposition of depletion-force curves across $\phi_l = 0.1$ to $0.4$ survives in every case; if it breaks down at some softness, the geometric explanation of the null result fails even if the empirical rule remains useful. Independently, hold two large particles near contact with a third large particle present and measure the three-body depletion force directly: at $q$ just below $q_3$ it should vanish if the no-triple-contact argument is right.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is that the many-body content of a depletion force is a quantitative, size-ratio-controlled quantity, and that it can be measured order by order. In the semi-grand ensemble, with the chemical potential of the small particles fixed by a particle-insertion method, the CBF simulations show that the effective two-body depletion force between large colloids is independent of the large-particle packing fraction $\phi_l$ whenever $q \le q_3 = 2/\sqrt{3} - 1 \approx 0.1547$: the curves for $\phi_l = 0.1, 0.2, 0.3, 0.4$ superpose within numerical error, directly verifying the geometric prediction that a small particle can never simultaneously touch three large ones below this threshold. The same scheme, by discarding configurations in which large particles form triangles, isolates the purely two-body contribution and shows three-body terms are negligible at $q = 0.15$ yet responsible for the density effect at $q = 0.60$. For $q = 0.45$ and $0.60$ the depletion force weakens as $\phi_l$ grows, about a 20% spread between $\phi_l = 0.1$ and $0.4$ at $q = 0.6$, and this dependence is well described by third-order polynomials in $\phi_l$ whose coefficients can extrapolate the force to the infinite-dilution limit $\phi_l \to 0$. The independent Ornstein-Zernike integral-equation route, closed with the modified-Verlet approximation, quantitatively reproduces the CBF results, and the authors present the explicit separation of two-, three-, and higher-order contributions as the work's most significant outcome.

Load-bearing premise

The load-bearing premise, stated in Section II as 'the previous geometrical arguments extend approximately to semi-hard spheres,' is that the slightly soft repulsive potential used in the simulations behaves enough like hard spheres for the geometric threshold $q_3 = 2/\sqrt{3} - 1 \approx 0.1547$, derived for perfect hard spheres, to apply unchanged to the simulated particles; the finite slope of the repulsive force in the measured curves shows the softness is real, and the paper provides no systematic check of how it shifts the effective threshold.

Editorial extensions

If this is right

  • For any binary mixture with $q \le 0.15$, effective potentials are pair-additive up to large-particle packing fractions of at least 0.4, so models of phase behavior in that regime can omit three-body and higher terms.
  • At fixed $\mu_s$, a depletion-force measurement made once at low large-particle density transfers unchanged to large-particle concentrations across the investigated range (up to $\phi_l = 0.4$), shortening the simulation or experimental effort needed to characterize such mixtures.
  • The published third-order polynomial coefficients for $q = 0.45$ and $0.60$ let other researchers interpolate depletion forces at un-simulated packing fractions and extrapolate them to $\phi_l \to 0$ without running new simulations.
  • Because the effective two-body force is shown to carry many-body contributions at the bare-interaction level, the CBF decomposition yields an explicit budget of how much of any measured attraction is pairwise versus three-body and higher at each state point.
  • For $q > q_3$, the weakening of the attraction as dense aggregates form could frustrate depletion-driven phase transitions, with condensation weakening the very attraction that drives it, a possibility the authors flag as worth exploring.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Applying the same fixed-$\mu_s$ protocol at $q \approx 0.2$, between $q_3$ and $q_4$, should reveal a staircase of many-body orders: two- and three-body terms present, four-body terms still absent. The paper does not run this test.
  • The superposition test itself could serve as a calibrated probe of how faithfully a soft potential mimics hard-sphere geometry, since softness visibly rounds the force curves and may shift the effective threshold.
  • The screening of the depletion attraction by large-particle density suggests a density-dependent effective pair potential, which coarse-grained dynamical simulations could absorb as a local rescaling of the attraction.
  • The close agreement between the integral-equation route and CBF in the many-body regime suggests the cheap theory could screen parameter space for simulations, while the residual discrepancies when bridge functions are removed locate precisely where the closure relation needs improvement.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. This manuscript studies depletion forces in binary mixtures of quasi-hard (mWCA) spheres, with size ratios q=0.15, 0.45, and 0.60. The authors use the contraction-of-bare-forces (CBF) method in the semi-grand ensemble, fixing the chemical potential of the depletant by Widom insertion, and compare with Ornstein–Zernike integral-equation theory with a modified Verlet closure. The central claims are: (i) for q≤0.15 (below the hard-sphere geometric threshold q3=2/√3−1), the depletion force is independent of the large-particle packing fraction φ_l at fixed μ_s, confirming a purely geometric prediction; (ii) for q=0.45 and 0.60 the force depends strongly on φ_l and can be fitted/extrapolated with low-degree polynomials; and (iii) the CBF scheme distinguishes two-body from three- and higher-body contributions by comparing runs with and without triangular large-particle configurations. The paper also reports agreement between CBF and OZ theory.

Significance. The significance of a verified criterion for pair-additive depletion potentials is considerable for coarse-grained modeling of colloid–polymer mixtures. The paper's strengths are its direct force-based methodology (avoiding integration drift), the semi-grand canonical protocol, the systematic scan over φ_l and μ_s, and the independent OZ cross-check. The promise of an explicit order-by-order separation of many-body contributions, if realized, would be a substantial methodological advance. However, in its present form the evidence supports the qualitative trends but not the two strongest statements: the hard-sphere interpretation of the q=0.15 null result and the claim of an explicit two/three/higher-order separation. With the additional quantitative checks and claim revisions requested below, the paper would be a solid contribution to the soft-matter literature.

major comments (4)
  1. [Section II and IV.B] The paper's central confirmation claim rests on the assertion that the geometric threshold q3=2/√3−1 applies to the simulated mWCA particles at kBT/ε=1. Section II states that for semi-hard spheres "the eventual touching of three larger spheres with a smaller one will be extremely unlikely," but no quantitative estimate is given. At q=0.15 the small–large contact distance is σ_sl=0.575σ_l, while three large spheres compressed to r_ll≈0.95σ_l create a triangular interstice whose centroid is about 0.549σ_l from each large center; the small–large repulsion at that point is already nonzero, so such three-body configurations are not exponentially suppressed at kBT/ε=1. To support the interpretation of Fig. 3(a) as a confirmation of the hard-sphere prediction, the authors need to characterize the effective softness of the mWCA potential (e.g., an effective hard-sphere diameter or Barker–Henderson mapping) and show that the threshold remains below q=0.15, or, alternatively, state the independence as an empirical observation independent of the geometric argument.
  2. [Fig. 3(a) and Section V.B] The null result at q=0.15—that the depletion force is independent of φ_l under fixed μ_s—is asserted from visual superposition without error bars. The text says the curves "become indistinguishable within the level of numerical error" (Fig. 3 caption), but no measure of the numerical error or a bound on a possible φ_l drift is provided. The paper should show statistical error bars (the ten independent realizations mentioned in §IV.B make this straightforward) and, ideally, a quantitative comparison such as the maximum difference between φ_l curves divided by the estimated uncertainty. This is load-bearing because the abstract and conclusions claim verification "with numerical precision."
  3. [Section I (final paragraph), Section III, V.C] The paper claims that the CBF scheme "is able to explicitly distinguish between the two-, three-, and higher-order many-body contributions" and calls this the most significant result. The actual analysis, however, compares the full force (including all many-body effects) with a force measured after excluding configurations in which large particles form triangles; it does not separately extract the f^(2), F^(3), F^(4), ... terms of Eq. (9). The comparison can at most separate "with" from "without" triangular configurations, and even that separation is not an order-by-order decomposition. Either provide a density-resolved fit of Eq. (9) that yields individual F^(n), or revise the claim to say that the method separates the total many-body contribution from the pair contribution.
  4. [Section V.C and Fig. 7] The operational definition of the noMBC sampling is not given. Fig. 7 is a sketch; the text says "configurations where large particles do not form triangles" are discarded, but no threshold distance or graph criterion is specified for what counts as a triangle. Because the exclusion conditions the sample on the absence of triangular large-particle arrangements, it can bias the measured pair force (e.g., through correlations between pair distance and third-particle proximity). The authors should define the criterion precisely (e.g., three large particles with mutual separations below a cutoff r_c) and test the sensitivity of the noMBC results to r_c.
minor comments (6)
  1. [Throughout] Throughout the text there are typos ("a priory" in Section I, "interations" in Ref. [8], "Webwer" in Ref. [2], "unified" in the abstract); the manuscript would benefit from a careful proofread.
  2. [Figs. 8 and 9] In Figs. 8 and 9 the legends use "BMBC" and "WMBC" but these abbreviations are never defined; the captions only define MBC and noMBC. These should be replaced by "with MBC" and "without MBC" or clarified.
  3. [Table II] The caption of Table II states that error bars are "not shown" and can be large for C and D; since these coefficients are used for extrapolation in Fig. 6, the fitted coefficients should be reported with standard errors.
  4. [Figs. 8(a) and 8(b)] The discrepancy between the CBF and OZ noMBC results in Figs. 8(a) and 8(b) (no visible difference in CBF, clear reduction in OZ) is attributed to the inadequacy of removing the bridge function, but the text does not explain which of the two procedures is a more faithful estimate of the true noMBC force; a few sentences of quantitative discussion would help.
  5. [Figures 3–5] The figures use color gradients to represent φ_l values, which is difficult to read in print; additional line styles or symbols would improve clarity.
  6. [Reference [34]] Reference [34] provides a virial-based mapping between hard-sphere and continuous potentials and could be used to justify or calibrate the mWCA potential as a hard-sphere surrogate; the authors only cite it, without applying it.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the central null result is a direct simulation observation benchmarked against an external geometric prediction, and self-cited methods are tools rather than the target.

full rationale

The paper's central claim, that the depletion force between large colloids is independent of the large-particle packing fraction for q <= 0.15 at fixed chemical potential of the small species, is established by direct simulation using the CBF procedure (Appendix A.1). The CBF method is a least-squares reduction of measured instantaneous forces, not a fit of the conclusion; the observed superposition of force curves in Fig. 3(a) is raw simulation output. The q3 threshold is an external geometric result from Dijkstra et al. (refs. 31-33), not derived or fitted in this work. Self-citations to the authors' own CBF and integral-equation closures are used as methodological tools and are cross-checked against each other, but the physics claim does not reduce to the validity of any fitted parameter or to a self-citation chain. The softness of the mWCA potential relative to hard spheres is a modeling approximation and a potential correctness risk, but it is not a circular reduction: the empirical superposition would stand even if the geometric interpretation were imperfect.

Assumptions & free parameters 2 free parameters · 5 assumptions · 0 invented entities

The central claims rest on standard statistical mechanics expansions, an empirical closure relation, and two domain approximations: the soft-sphere stand-in for hard spheres, and the triangle-exclusion rule for isolating many-body effects. No new physical entities are introduced.

free parameters (2)
  • Cubic polynomial fit coefficients A, B, C, D for depletion forces = Reported in Table II for combinations of phi_s^res = 0.025, 0.050, 0.075, 0.100 and phi_l = 0.1, 0.2, 0.3, 0.4
    Coefficients of third-order polynomial fits to the attractive region of depletion forces (1 < r/sigma_l < 1.4). Used to extrapolate to phi_l -> 0. Errors can be large for C and D, as acknowledged in Table II caption.
  • mWCA potential parameters epsilon and sigma = epsilon = 1, sigma_l = 1, sigma_s = q
    Model parameters for the modified Weeks-Chandler-Andersen potential. Treated as inputs, not fitted to the target result, but the mapping to hard-sphere behavior is assumed.
assumptions (5)
  • standard math The effective many-body potential can be expanded in cluster functions U(n) and the effective force in a power series in the large-particle density rho_0 (Eqs. 1-9).
    Standard statistical mechanics of effective interactions; the series truncation is an approximation tested by the polynomial fits.
  • domain assumption The contracted effective force between two large particles depends only on their scalar separation r_ij (Eq. A2).
    Assumes isotropy and pairwise-additive contracted forces; standard in effective interaction approaches.
  • domain assumption The mWCA potential with k_BT/epsilon = 1 approximates hard-sphere behavior well enough for the geometric q3 and q4 thresholds to apply.
    Invoked in Section II for semi-hard spheres; softness is visible in the finite slopes of the repulsive force, so the approximation is not exact.
  • domain assumption The modified-Verlet closure (Eq. A14) is accurate for short-range repulsive potentials.
    Empirical closure relation used in the OZ theory; previously tested by the authors in Refs. [14,15,27,29,52].
  • ad hoc to paper Excluding configurations where large particles form triangles removes three-body and higher-order many-body contributions from the CBF force measurement without other bias (Fig. 7).
    This is the new heuristic procedure used to define noMBC; its validity is not independently established.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Quantifying many-body contributions to depletion forces." pith.science (2026). https://pith.science/paper/FIXHF3GK

@misc{pith2026250903342,
  author       = {Pith},
  title        = {Pith review of: Quantifying many-body contributions to depletion forces},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/FIXHF3GK}},
  note         = {Machine review of arXiv:2509.03342}
}
abstract

Effective interactions inherently encompass many-body effects that appear unified. Analyzing these in reverse, that is, separating them into contributions from pairs, triples, or larger groups, is typically intricate and seldom pursued. However, this could offer essential insights into the structural architecture of complex systems, such as soft materials. This contribution tackles this issue employing a new simulation-based approach to accurately determine effective interactions. The approach is sufficiently sensitive to assess the effects of higher-order terms in the depletion potentials between large particles. Previous research has primarily focused on the role of small particles, and this work expands on these findings. However, understanding the contributions of large particles on the effective forces, particularly beyond the dilute limit, remains challenging and is not yet fully grasped, leading us to primarily focus on exploring this topic. Within the range of particle concentrations examined here, while maintaining a constant chemical potential for smaller particles, we have observed that the concentration of larger particles has no impact on the entropic potential between large colloids, as long as the size ratio remains below $q=0.15$. This confirms a long-established prediction founded on purely geometric considerations, which we have confirmed by means of direct observation. In contrast, for mixtures with less size asymmetry, such effects become considerably influential. Specifically, we have conducted an in-depth examination of scenarios with size asymmetries of $q=0.45$ and $0.60$. Our results for the depletion forces were directly compared with those obtained from the integral equation theory, which has also enabled us to improve and refine the approaches involved.

Figures

Figures reproduced from arXiv: 2509.03342 by the authors.

Figure 1
Figure 1. FIG. 1. Curves of [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Curves of [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. (b) displays the depletion forces obtained from the integral equations formalism for the cases reported in [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figures from the paper (3 more)
Figure 5
Figure 5. Figure 5: (a). Insets show an enlargement of the negative-forces region for ϕ res s = 0.075. ϕl . A second fitting, now of the coefficients themselves as functions of ϕl , can be used to extrapolate to the limit ϕl → 0, allowing us to predict in this way the shape of the depleti…
Figure 7
Figure 7. Figure 7: FIG. 7. Representative diagram of how we proceed in compu [PITH_FULL_IMAGE:figures/full_fig_p010_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Depletion forces, [PITH_FULL_IMAGE:figures/full_fig_p011_8.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. A General Theory for Phenotypic Association in Biological Systems

    physics.bio-ph 2026-08 conditional novelty 6.0 of 10

    Multivalent binding maps exactly onto the monomer-dimer problem, so recognition thresholds are necessarily smooth crossovers with log-concave bond distributions, and a single free-energy functional with many-body coup...

Reference graph

Works this paper leans on

56 extracted references · 55 canonical work pages · cited by 1 Pith paper

  1. [54]

    Perera-Burgos, J

    J. Perera-Burgos, J. M. M´ endez-Alcaraz, G. P´ erez-´Angel, and R. Casta˜ neda-Priego, Assessment of the micro struc- ture and depletion potentials in two dimensional binary mixtures of additive hard disks, J. Chem. Phys. 145, 104905 (2016)

  2. [1]

    best fit

    Contraction of the bare forces Effective interactions, and in particular depletion in- teractions, can be determined at the level of forces using an average of the effects of the small particle-big parti- cle forces in the dynamics of the large ones. This is an approach that we have named contraction of the forces (CBF) in previous works [14, 15, 29, 50]....

  3. [2]

    Integral equations theory for effective interactions In the following, we refer to the integral equations the- ory for effective interactions in complex fluids [25, 26, 29], which takes advantage of the covariance of the Ornstein- Zernike equation under contractions of the description. For homogeneous mixtures that contain p species of spherical particles...

  4. [3]

    B. M. Axilrod and E. Teller, Interaction of the van der Waals Type Between Three Atoms , J. Chem. Phys. 11, 299 (1943)

  5. [4]

    F. H. Stillinger and T. A. Webwer, Computer simulation of local order in condensed phases of silicon , Phys. Rev. B 31, 5262 (1985)

  6. [5]

    Bukowsky and K

    R. Bukowsky and K. Szalewicz, Complete ab initio three- body nonadditive potential in Monte Carlo simulations of vapor-liquid equilibria and pure phases of argon, J. Chem. Phys. 114, 9518 (2001)

  7. [6]

    Yan, L.-Y

    P.-G. Yan, L.-Y. Tang, Z.-C. Yan, and J. F. Babb, Long- range additive and nonadditive potentials in a hybrid sys- tem: Ground-state atom, excited-state atom, and ion , Phys. Rev. A 104, 022807 (2021)

  8. [7]

    GROMACS developement team, GROMACS documen- tation, Release 2023.3 (2023), available at https:// manual.gromacs.org/2023.3/manual-2023.3.pdf

Show all 56 references
  1. [8]

    LAMMPS Documentation, (2024), available at https: //docs.lammps.org/Manual.html

  2. [9]

    Goulding and S

    D. Goulding and S. Melchionna, Accurate calculation of three-body depletion interactions , Phys. Rev. E 64, 011403 (2001)

  3. [10]

    Frischknecht and A

    A. Frischknecht and A. Yethiraj, Two- and three-body interations among nanoparticles in a polymer melt , J. Chem. Phys. 134, 174901 (2011)

  4. [11]

    D. Zhu, W. Li, and H. R. Ma,On pair additivity of the de- pletion force, J. Phys.: Condens. Matter 15, 8281 (2003)

  5. [12]

    Lowen and E

    H. Lowen and E. Allahyarov, The role of effective triplet interactions in charged colloidal suspensions , J. Phys.: Condens. Matter 10, 4147 (1998)

  6. [13]

    D. J. Ashton and N. B. Wilding, Three-body interactions in complex fluids: Viral coefficients from finite size ef- fects, J. Chem. Phys. 140, 244118 (2014)

  7. [14]

    Santos, M

    A. Santos, M. L´ opez de Haro, G. Fiumara, and F. Saija, The effective colloid interaction in the Asakura-Oosawa model. Assessement of non-pairwise terms from the virial expansion, J. Chem. Phys. 142, 224903 (2015)

  8. [15]

    Dobnikar, R

    J. Dobnikar, R. Casta˜ neda-Priego, H. H. von Gr¨ unberg, and E. Trizac, Testing the relevance of effective interac- tion potentials between highly-charged colloids in suspen- sion, New Journal of Physics 8, 277 (2006)

  9. [16]

    N. M. de los Santos-L´ opez, G. P´ erez- ´Angel, R. Casta˜ neda-Priego, and J. M. M´ endez-Alcaraz, Determining depletion interactions by contracting forces , J. Chem. Phys. 157, 074903 (2022)

  10. [17]

    N. M. de los Santos-L´ opez, M. A. Ram´ ırez- Gu´ ızar, G. P´ erez-´Angel, J. M. M´ endez-Alcaraz, and R. Casta˜ neda-Priego,Depletion forces beyond the diluted limit, Physica A 655, 130180 (2024)

  11. [18]

    Asakura and F

    S. Asakura and F. Oosawa, On Interaction between Two Bodies Immersed in a Solution of Macromolecules , J. Chem. Phys. 22, 1255 (1954)

  12. [19]

    Asakura and F

    S. Asakura and F. Oosawa, Interaction between Particles Suspended in Solutions of Macromolecules, J. Polym. Sci. 33, 183 (1958)

  13. [20]

    Vrij, Polymers at interfases and the interactions in colloidal dispersions, Pure Appl

    A. Vrij, Polymers at interfases and the interactions in colloidal dispersions, Pure Appl. Chem. 48, 471 (1976)

  14. [21]

    Oosawa, The history of the birth of the Asakura- Oosawa theory, J

    F. Oosawa, The history of the birth of the Asakura- Oosawa theory, J. Chem. Phys. 155, 084104 (2021)

  15. [22]

    A. P. Gast, C. K. Hall, and W. B. Russell, Polymer- Induced Phase Separations in Nonaqueous Colloidal Sus- pensions, J. Colloid. Interphase Sci. 96, 251 (1983)

  16. [23]

    H. N. W. Lekkerkerker and R. Tuinier, Colloids and the depletion interaction (Springer, 2011)

  17. [24]

    Dijkstra and R

    M. Dijkstra and R. van Roij, Entropic Wetting and Many-Body Induced Layering in a Model Colloid- Polymer Mixture, Phys. Rev. Lett. 89, 208303 (2002)

  18. [25]

    Dijkstra, R

    M. Dijkstra, R. van Roij, R. Roth, and A. Fortini, Ef- fect of many-body interactions on the bulk and interfacial phase behavior of a model colloid polymer mixture , Phys. Rev. E 73, 041404 (2006)

  19. [26]

    R. Roth, R. Evans, and S. Dietrich, Depletion potential in hard-sphere mixtures: Theory and applications, Phys. Rev. E 62, 5360 (2000)

  20. [27]

    J. M. M´ endez-Alcaraz and R. Klein, Depletion forces in colloidal mixtures, Phys. Rev. E 61, 4095 (2000)

  21. [28]

    Casta˜ neda-Priego, A

    R. Casta˜ neda-Priego, A. Rodr´ ıguez-L´ opez, and J. M. M´ endez-Alcaraz, Entropic forces in dilute colloidal sys- tems, Phys. Rev. E 73, 051404 (2006)

  22. [29]

    L´ opez-S´ anchez, C

    E. L´ opez-S´ anchez, C. D. Estrada- ´Alvarez, G. P´ erez- ´Angel, J. M. M´ endez-Alcaraz, P. Gonz´ alez-Mozuelos, and R. Casta˜ neda-Priego, Demixing transition, struc- ture, and depletion forces in binary mixtures of hard spheres: The role of bridge functions, J. Chem. Phys....

  23. [30]

    Heinen, Calculating particle pair potentials from fluid-state pair correlations: Iterative ornstein–zernike inversion, Journal of Computational Chemistry 39, 1531 (2018)

    M. Heinen, Calculating particle pair potentials from fluid-state pair correlations: Iterative ornstein–zernike inversion, Journal of Computational Chemistry 39, 1531 (2018)

  24. [31]

    N. M. de los Santos-L´ opez, G. P´ erez- ´Angel, J. M. M´ endez-Alcaraz, and R. Casta˜ neda-Priego,Competing interations in the depletion forces of ternary colloidal mixtures, J. Chem. Phys. 155, 024901 (2021)

  25. [32]

    C. N. Likos, Effective interations in condensed soft matter physics, Phys. Rep. 348, 267 (2001)

  26. [33]

    Dijkstra, R

    M. Dijkstra, R. van Roij, and R. Evans, Phase Behav- ior and Structure of Binary Hard Sphere Mixtures , Phys. Rev. Lett. 81, 2268 (1998)

  27. [34]

    Dijkstra, R

    M. Dijkstra, R. van Roij, and R. Evans, Phase diagram of highly asymmetric binary hard sphere mixtures , Phys. Rev. E 59, 5744 (1999)

  28. [35]

    Dijkstra, J

    M. Dijkstra, J. M. Brader, and R. Evans,Phase behaviour and structure of model colloid-polymer mixtures, J. Phys.: Condens. Matter 11, 10079 (1999)

  29. [36]

    C. A. B´ aez, A. Torres-Carbajal, R. Casta˜ neda-Priego, A. Villada-Balbuena, J. M. M´ endez-Alcaraz, and S. Herrera-Velarde, Using the second virial coefficient as physical criterion to map the hard-sphere potential onto a continuous potential, The Journal of Chemical Physics...

  30. [37]

    Dhont, An Introduction to Dynamics of Colloids, Stud- ies in Interface Science (Elsevier Science, 1996)

    J. Dhont, An Introduction to Dynamics of Colloids, Stud- ies in Interface Science (Elsevier Science, 1996)

  31. [38]

    Gardiner, Stochastic Methods, 4th ed., Springer Series in Synergetics, Vol

    C. Gardiner, Stochastic Methods, 4th ed., Springer Series in Synergetics, Vol. 4 (Springer Berlin, Heidelberg, 2009) pp. XVIII + 447, hardcover ISBN: 978-3-540-70712-7, Softcover ISBN: 978-3-642-08962-6. Series ISSN: 0172- 7389, E-ISSN: 2198-333X

  32. [39]

    Kolafa and S

    J. Kolafa and S. Lab´ ık, Density expansion of the radial distribution and bridge functions of the hard sphere fluid, Molecular Physics 104, 1915–1924 (2006)

  33. [40]

    B. R. A. Nijboer and L. V. Hove, Radial distribution function of a gas of hard spheres and the superposi- tion approximation, Phys. Rev. Journals Archive 85, 777 (1952)

  34. [41]

    Santos, A Concise Course on the Theory of Classi- cal Liquids , 1st ed

    A. Santos, A Concise Course on the Theory of Classi- cal Liquids , 1st ed. (Springer International Publishing, 2016)

  35. [42]

    M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids, 2nd ed. (Oxford U. Press, 2017)

  36. [43]

    Boinepalli and P

    S. Boinepalli and P. Attard, Grand Canonical Molecular Dynamics, J. Chem. Phys. 119, 12769 (2003)

  37. [44]

    Eslami and F

    H. Eslami and F. M¨ uller-Plathe, Molecular Dynamics Simulation in the Grand Canonical Ensemble , J. Com- put. Chem. 28, 1763 (2007)

  38. [45]

    Belloni, Non-equilibrium hybrid insertion/extraction through the 4th dimension in grand-canonical simulation, J

    L. Belloni, Non-equilibrium hybrid insertion/extraction through the 4th dimension in grand-canonical simulation, J. Chem. Phys. 151, 021101 (2019)

  39. [46]

    Widom, Some Topics in the Theory of Fluids , J

    B. Widom, Some Topics in the Theory of Fluids , J. Chem. Phys. 39, 2808 (1963)

  40. [47]

    Frenkel, Simulations: The dark side , Eur

    D. Frenkel, Simulations: The dark side , Eur. Phys. J. Plus 128, 10 (2013)

  41. [48]

    Bussi, D

    G. Bussi, D. Donadio, and M. Parrinello, Canonical sam- pling through velocity rescaling , J. Chem. Phys. 126, 014101 (2007)

  42. [49]

    J. P. Hansen and I. R. McDonald, Theory of Simple Liq- uids, 2nd ed. (Academic Press, 1986)

  43. [50]

    P. E. Kloeden and E. Platen, Numerical Solution of Stochastic Differential Equations (Springer Berlin Hei- delberg, 1992)

  44. [51]

    H. S. Nguyen, J. Forsman, and C. E. Woodward, Many- body depletion forces of colloids in a polydisperse polymer dispersant in the long-chain limit , Soft Matter 14, 6921 (2018)

  45. [52]

    J. J. Pedrozo-Romero and G. P´ erez- ´Angel, Depletion forces in dense mixtures of spheres and rods , J. Chem. Phys. 155, 024901 (2022)

  46. [53]

    W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, Numerical Recipes in C, 2nd ed. (Cambridge U. Press, 1992)

  47. [55]

    Ng, Hypernetted chain solutions for the classical one-component plasma up to γ = 7000, J

    K.-C. Ng, Hypernetted chain solutions for the classical one-component plasma up to γ = 7000, J. Chem. Phys. 61, 2680 (1974)

  48. [56]

    E. A. Bedolla-Montiel, R. A. Casta˜ neda-Cerd´ an, and R. Casta˜ neda-Priego, Phase diagram of the hard-sphere potential model in three and four dimensions using a pseudo-hard-sphere potential, The Journal of Chemical Physics 162, 1234567 (2025). 16 TABLE II. Coefficients for ...

Pith tools

Reviewed August 15, 2026 · model on record in the stance chip above.