pith. sign in

arxiv: 1906.12062 · v1 · pith:I4GV2MDLnew · submitted 2019-06-28 · ❄️ cond-mat.stat-mech · cond-mat.soft

Effect of aggregation on adsorption phenomena

Marek Litniewski , Alina Ciach This is my paper

Pith reviewed 2026-05-25 13:47 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech cond-mat.soft
keywords adsorptionaggregationclustersmolecular dynamicsattractive surfacedepletion regionpower-law dependenceYukawa interaction
0
0 comments X

The pith

The temperature marking the onset of cluster adsorption near an attractive wall increases as a power law with wall-particle attraction strength.

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

This paper uses molecular dynamics simulations to study how particles that self-assemble into small clusters adsorb onto an attractive surface at low densities. The key finding is that the critical temperature for the onset of significant cluster adsorption near the surface increases relative to the bulk value according to a power law in the strength of the wall attraction. Below this temperature, the surface layer consists of intact clusters for moderate attractions, but stronger attractions cause flattening into monolayers or stripes. The growing layer creates a repulsive barrier that dramatically slows further adsorption, and a depletion zone forms adjacent to the wall. These results demonstrate the interplay between bulk aggregation and surface phenomena in the chosen interaction model.

Core claim

In systems with particles forming small clusters via Lennard-Jones plus repulsive Yukawa interactions at low densities, the relative increase of the temperature at the critical cluster concentration near the attractive surface follows a power-law dependence on the wall-particle attraction strength. Below the CCCS temperature the adsorbed layer has undeformed clusters if attraction is not too strong; above it or for strong attraction, a monolayer of flattened clusters or stripes forms. Accumulated repulsion slows adsorption with accession time growing rapidly with attraction strength, while a depletion region of thickness comparable to the Yukawa range appears with density decreasing as wall,

What carries the argument

The critical cluster concentration near the attractive surface (CCCS), defined via the temperature where significant cluster adsorption begins, whose relative shift obeys a power law in wall attraction strength and is measured in MD simulations of the Lennard-Jones plus repulsive Yukawa model.

Load-bearing premise

The Lennard-Jones plus repulsive Yukawa interaction model at small densities produces the intended small clusters and molecular dynamics trajectories adequately sample equilibrium adsorption states without significant finite-size or equilibration artifacts.

What would settle it

A set of simulations or measurements at several wall attraction strengths showing that the relative CCCS temperature increase does not follow a power law would falsify the central claim.

Figures

Figures reproduced from arXiv: 1906.12062 by Alina Ciach, Marek Litniewski.

Figure 1
Figure 1. Figure 1: Fig.1 [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 1
Figure 1. Figure 1: FIG. 1: The interaction potential [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Histograms for the probability of finding a particle i [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Cartoon showing the clusters composed of 5 and 7 parti [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Pair distribution function [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5 [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6: The ratio of the critical cluster concentration temp [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7: The density profile for [PITH_FULL_IMAGE:figures/full_fig_p014_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8: The density profile for [PITH_FULL_IMAGE:figures/full_fig_p015_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9: The inverse accessibility time [PITH_FULL_IMAGE:figures/full_fig_p015_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10: The adsorption (9) for [PITH_FULL_IMAGE:figures/full_fig_p016_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11: A projection of a representative configuration of th [PITH_FULL_IMAGE:figures/full_fig_p017_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12: 2D pair distribution function for centers of mass fo [PITH_FULL_IMAGE:figures/full_fig_p018_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: FIG. 13: A projection of the representative configuration of [PITH_FULL_IMAGE:figures/full_fig_p018_13.png] view at source ↗
read the original abstract

Adsorption at an attractive surface in a system with particles self-assembling into small clusters is studied by Molecular dynamics (MD) simulation. We assume Lennard-Jones plus repulsive Yukawa tail interactions, and focus on small densities. The relative increase of the temperature at the critical cluster concentration near the attractive surface (CCCS) shows a power-law dependence on the strength of the wall-particle attraction. At temperatures below the CCCS, the adsorbed layer consists of undeformed clusters if the wall-particle attraction is not too strong. Above the CCCS, or for strong attraction leading to flattening of the adsorbed aggregates, we obtain a monolayer that for strong or very strong attraction consists of flattened clusters or stripes respectively. The accumulated repulsion from the particles adsorbed at the wall leads to a repulsive barrier that slows down the adsorption process, and the accession time grows rapidly with the strength of the wall-particle attraction. Beyond the adsorbed layer of particles, a depletion region of a thickness comparable with the range of the repulsive tail of interactions occurs, and the density in this region decreases with increasing strength of the wall-particle attraction. At larger separations, the exponentially damped oscillations of density agree with theoretical predictions for self-assembling systems. Structural and thermal properties of the bulk are also determined. In particular, a new structural crossover associated with the maximum of the specific heat, and a double-peaked histogram of the cluster size distribution are observed.

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

2 major / 2 minor

Summary. The manuscript reports molecular dynamics simulations of low-density particles with Lennard-Jones plus repulsive Yukawa interactions adsorbing onto an attractive wall. The central claim is that the relative increase in the temperature at the critical cluster concentration near the surface (CCCS) follows a power-law dependence on wall-particle attraction strength. Below CCCS the adsorbed layer contains undeformed clusters (for moderate attraction); above CCCS or for strong attraction the layer forms a monolayer of flattened clusters or stripes. A repulsive barrier from adsorbed particles slows adsorption, with accession time growing rapidly with attraction strength; a depletion zone and exponentially damped density oscillations beyond it are reported and compared to theory. Bulk structural and thermal properties, including a specific-heat maximum and double-peaked cluster-size histograms, are also presented.

Significance. If equilibration is demonstrated, the work supplies a concrete, falsifiable power-law relation between wall attraction and the shift in CCCS temperature together with direct comparison of simulated density oscillations to existing theoretical predictions for self-assembling fluids. These elements would constitute a useful benchmark for theories of adsorption in aggregating systems.

major comments (2)
  1. [Abstract] Abstract: the power-law claim for the relative CCCS temperature increase is load-bearing for the central result, yet the same paragraph states that accession time grows rapidly with wall attraction because of the accumulated repulsion barrier. This directly raises the possibility that MD trajectories for stronger attractions have not reached equilibrium adsorption states, undermining reliable location of the specific-heat maximum or cluster-size histogram used to define CCCS.
  2. [Abstract] Abstract (and implied Methods/Results sections): no information is given on trajectory lengths, number of independent runs, finite-size checks, or explicit equilibration diagnostics (e.g., time series of cluster-size histograms or energy). At the low densities studied, such controls are essential to rule out the sampling artifacts flagged by the weakest assumption.
minor comments (2)
  1. [Abstract] Abstract: the acronym CCCS is introduced without an explicit parenthetical expansion on first use.
  2. [Abstract] Abstract: the exponent of the reported power-law is not stated, nor is the range of attraction strengths over which the fit was performed.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and for highlighting important issues regarding the demonstration of equilibration in our molecular dynamics simulations. The concerns are well-taken and point to areas where the manuscript can be strengthened by providing additional methodological details. We address each major comment below and will incorporate the necessary revisions.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the power-law claim for the relative CCCS temperature increase is load-bearing for the central result, yet the same paragraph states that accession time grows rapidly with wall attraction because of the accumulated repulsion barrier. This directly raises the possibility that MD trajectories for stronger attractions have not reached equilibrium adsorption states, undermining reliable location of the specific-heat maximum or cluster-size histogram used to define CCCS.

    Authors: We agree that the rapid growth of accession time with attraction strength raises a legitimate question about whether equilibrium was reached for the strongest attractions. In our simulations, run lengths were chosen to substantially exceed the observed accession times for each wall attraction value, allowing the system to reach a regime where the specific-heat maximum and the double-peaked cluster-size histograms became stationary. The power-law dependence was extracted only from those equilibrated portions of the trajectories. Nevertheless, to remove any ambiguity we will add explicit statements in the revised manuscript confirming that the reported CCCS temperatures correspond to post-equilibration data. revision: yes

  2. Referee: [Abstract] Abstract (and implied Methods/Results sections): no information is given on trajectory lengths, number of independent runs, finite-size checks, or explicit equilibration diagnostics (e.g., time series of cluster-size histograms or energy). At the low densities studied, such controls are essential to rule out the sampling artifacts flagged by the weakest assumption.

    Authors: The referee is correct that the current manuscript does not report trajectory lengths, number of independent runs, finite-size checks, or explicit equilibration diagnostics. We will revise the Methods section to include these details: typical production runs of 5×10^7–2×10^8 MD steps (after equilibration), at least three independent realizations per state point, system-size checks up to N=8000 particles, and time-series monitoring of total energy, wall adsorption, and cluster-size histograms to verify convergence. These additions will directly address the sampling concerns at the low densities examined. revision: yes

Circularity Check

0 steps flagged

No significant circularity; simulation observations compared to external theory

full rationale

The paper is a molecular dynamics simulation study reporting numerical observations of adsorption and cluster behavior in an LJ+Yukawa system. The central claim (power-law dependence of relative CCCS temperature increase on wall attraction) is presented as a direct simulation result, not a fitted parameter renamed as prediction or a self-definitional relation. Density profile oscillations are stated to agree with external theoretical predictions for self-assembling systems, providing an independent benchmark. Bulk properties such as specific-heat maxima and cluster-size histograms are reported as observed features without reduction to prior self-citations or ansatzes. No equations, self-citations, or derivation steps in the abstract or described content reduce the outputs to the inputs by construction. The work is self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The study rests on a standard pair-potential model chosen to produce clustering and on the assumption that MD yields representative equilibrium statistics; no new entities are postulated.

free parameters (2)
  • Lennard-Jones well depth and range
    Chosen to set the attraction scale for clustering at low density.
  • Yukawa repulsion amplitude and decay length
    Selected to produce small stable clusters rather than macroscopic phase separation.
axioms (1)
  • domain assumption Particles interact via Lennard-Jones attraction plus repulsive Yukawa tail at small densities
    Core model assumption enabling cluster formation and adsorption study.

pith-pipeline@v0.9.0 · 5775 in / 1134 out tokens · 30082 ms · 2026-05-25T13:47:45.119676+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

28 extracted references · 28 canonical work pages

  1. [1]

    Imperio and L

    A. Imperio and L. Reatto, J. Phys.: Condens. Matter 16, S3769 (2004)

    work page 2004

  2. [2]

    A. J. Archer and N. B. Wilding, Phys. Rev. E 76, 031501 (2007)

    work page 2007

  3. [3]

    A. J. Archer, Phys. Rev. E 78, 031402 (2008)

    work page 2008

  4. [4]

    Ciach, Phys

    A. Ciach, Phys. Rev. E 78, 061505 (2008)

    work page 2008

  5. [5]

    Ciach and W

    A. Ciach and W. T. G´ o´ zd´ z,Condens. Matter Phys. 13, 23603 (2010)

    work page 2010

  6. [6]

    Kowalczyk, A

    P. Kowalczyk, A. Ciach, P. A. Gauden, and A. P. Terzyk, Int. J. Colloid Interface Sci. 363, 579 (2011). 20

    work page 2011

  7. [7]

    Zhuang, K

    Y. Zhuang, K. Zhang, and P. Charbonneau, Phys. Rev. Lett. 116, 098301 (2016)

    work page 2016

  8. [8]

    Zhuang and P

    Y. Zhuang and P. Charbonneau, J. Phys. Chem. B 120, 6178 (2016)

    work page 2016

  9. [9]

    Edelmann and R

    M. Edelmann and R. Roth, Phys. Rev. E 93, 062146 (2016)

    work page 2016

  10. [10]

    Pini and A

    D. Pini and A. Parola, Soft Matter 13, 9259 (2017)

    work page 2017

  11. [11]

    A. I. Campbell, V. J.Anderson, J. S. van Duijneveldt, an d P. Bartlett, Phys. Rev. Lett. 94, 208301 (2005)

    work page 2005

  12. [12]

    Stradner et al

    A. Stradner et al. , Nature 432, 492 (2004)

    work page 2004

  13. [13]

    C. P. Royall, Soft Matter 14, 4020 (2018)

    work page 2018

  14. [14]

    Bergman, T

    M. Bergman, T. Garting, P. Schurtenberger, and A. Strad ner, J. Phys. Chem. B 123, 2432 (2019)

    work page 2019

  15. [15]

    A. P. Santos, J. P¸ ekalski, and A. Z. Panagiotopoulos, Soft Matter 13, 8055 (2017)

    work page 2017

  16. [16]

    Hu and P

    Y. Hu and P. Charbonneau, Soft Matter 14, 4101 (2018)

    work page 2018

  17. [17]

    D. Amos, J. Markels, S. Lynn, and C. Radke, J. Phys. Chem. B 102, 2739 (1998)

    work page 1998

  18. [18]

    M. A. Floriano, E. Caponetti, and A. Z. Panagiotopoulos , Langmuir 15, 3143 (1999)

    work page 1999

  19. [19]

    N. G. Almarza, J. P¸ ekalski, and A. Ciach, J. Chem. Phys. 140, 164708 (2014)

    work page 2014

  20. [20]

    P¸ ekalski, A

    J. P¸ ekalski, A. Ciach, and N. G. Almarza, J. Chem. Phys. 140, 114701 (2014)

    work page 2014

  21. [21]

    Litniewski, J

    M. Litniewski, J. Phys. Chem. 97, 3842 (1993)

    work page 1993

  22. [22]

    Morishita, J

    T. Morishita, J. Chem. Phys. 119, 7075 (2003)

    work page 2003

  23. [23]

    Collins, G

    P. Collins, G. S. Ezra and S. Wiggins, J. Chem. Phys. 133, 014105 (2010)

    work page 2010

  24. [24]

    M. P. Allen and D. J. Tildesley, Computer Simulations of Liquids (Clarendon Press, Oxford, 1990)

    work page 1990

  25. [25]

    P. D. Godfrin, N. E. Valadez-Perez, R. Castaneda-Prieg o, N.J. Wagner and Y. Liu, Soft Mater 10, 5061 (2014)

    work page 2014

  26. [26]

    J. B. Pedersen and P. Sibani, J. Chem. Phys. 75, 5368 (1981)

    work page 1981

  27. [27]

    Szabo, J.Phys

    A. Szabo, J.Phys. Chem. 93, 6929 (1989)

    work page 1989

  28. [28]

    Zhuang and P

    Y. Zhuang and P. Charbonneau, J. Chem. Phys. 147, 091102 (2017). 21

    work page 2017