Pith. sign in

REVIEW 4 major objections 5 minor 53 references

Unveiling the Role of Friction in Coarse-Grained Clay: A Hybrid Framework Integrating Long-Range Interactions and Granular Contact Mechanics

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

Pith's one-line read A hybrid coarse-grained model shows clay shear strength comes from inter-particle friction, not geometric interlocking alone.

desk verdict Useful and honest sensitivity study of damping and friction in CG clay, but the headline friction claim is largely baked into the force law and the repulsion-only potential weakens the broader conclusions. read the letter →

arxiv 2608.04359 v1 pith:HTOUQBJR submitted 2026-08-05 cond-mat.mtrl-sci

classification cond-mat.mtrl-sci
keywords coarse-grainedmoleculardynamicsmontmorilloniteinter-particlefrictionHertziancontactmodelBuckinghampotentialclayfabricviscoelasticdampinguniaxialcompression
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

This paper tries to establish that coarse-grained molecular dynamics of clay must include inter-particle friction and viscous damping, not just conservative interaction forces. It builds a hybrid framework in which a Buckingham potential handles long-range interactions and Hertzian contact mechanics handles short-range hydration contacts with explicit friction and damping. Using this model, the authors show that frictionless clay assemblies undergo fluid-like yielding under compression, while frictional assemblies sustain much higher loads, and that excessive damping or missing thermal fluctuations traps platelets in loose, disordered states. If right, this identifies friction as a primary driver of shear strength in clay models and provides a practical way to bridge atomistic potentials with granular contact physics.

What carries the argument

The central object is a hybrid inter-particle force law: a Buckingham potential $E = A e^{-r/\rho} - C/r^6$ for long-range van der Waals and electrostatic interactions, combined with a Hertzian contact force whose normal component is $k_n R_{\text{eff}}^{1/2}\delta_{ij}^{3/2} - \eta_n v_n$ and whose tangential component is capped by $\mu\|\mathbf{F}_n\|$. The Hertz term is integrated with the Buckingham term and fitted to reference energy curves, allowing particle overlap to represent hydration-shell compression while the tangential term provides explicit frictional resistance. This coupling is what lets friction and damping be tuned rather than smoothed away, and it carries the paper's argument that these dissipative forces, not conservative potentials alone, determine strength and fabric.

What would settle it

Run the same uniaxial compression test at a strain rate orders of magnitude slower with the friction coefficient set to zero; if the frictionless assembly then strain-hardens and reaches stresses close to the frictional case, the claim that friction is the primary strength source fails. An experimental counterpart would be shearing clay with deliberately lubricated grain contacts and checking whether the shear strength drops to near zero.

Watch

Extended reading notes

Core claim

The central claim is that inter-particle friction is the primary driver of shear strength and structural integrity in coarse-grained clay models: in uniaxial compression, a frictionless montmorillonite assembly yields at roughly 2 MPa with strain softening, whereas the same assembly with friction coefficient 0.1 locks sliding interfaces, strain-hardens, and sustains significantly higher loads. The paper also finds that viscoelastic damping controls the fabric produced by compression, with high normal damping nearly doubling the void ratio and suppressing platelet ordering, and that removing the thermostat causes damping forces to quench the assembly into metastable, low-density states. These results come from a hybrid coarse-grained model that couples the Buckingham potential with a Hertzian granular contact law, calibrated against reference energy curves for face-to-face, edge-to-edge, and face-to-edge platelet configurations.

Load-bearing premise

The load-bearing premise is that the reference energy curves used for calibration correctly represent sodium montmorillonite pairwise interactions and that zeroing the Buckingham C term does not remove an interaction that matters; if either fails, the validated baseline and all subsequent friction results inherit the error.

Editorial extensions

If this is right

  • Coarse-grained clay models that omit friction will systematically underestimate shear strength and may show fluid-like yielding instead of realistic jamming.
  • Damping coefficients must be chosen carefully: overdamped contacts freeze the clay fabric into high-void-ratio, disordered configurations.
  • Thermostat control is necessary during compression so thermal fluctuations can drive the assembly toward equilibrium; damping-only runs quench platelets into metastable states.
  • The calibrated hybrid model reproduces macroscopic bentonite compressibility, with a compression index of 2.17 at low-to-mid pressures and 0.39 at high pressures, supporting its use for geotechnical-scale predictions.
  • Explicit friction combined with long-range potentials offers a route to multi-scale clay simulation that avoids purely macroscopic empirical calibration.

Reading between the lines

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

  • Beyond the paper, the same hybrid coupling could be refitted to potentials for other clay minerals, which would test whether friction's dominance holds across different clay chemistries.
  • The stochastic domain orientation observed at low tangential damping hints that 1000 platelets may be too few for quantitative fabric predictions; larger assemblies could average local domains into a more isotropic fabric.
  • If frictionless coarse-grained clay models are widely used, their published strengths and consolidation behaviors may need to be rechecked against frictional runs.
  • The friction and damping coefficients are chosen from literature ranges rather than derived from the calibration data; fitting them to macroscopic strength measurements would make the framework predictive rather than demonstrative.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 5 minor

Summary. The paper proposes a coarse-grained molecular dynamics (CGMD) framework for montmorillonite assemblies that couples a long-range Buckingham potential with Hertzian granular contact mechanics, including normal and tangential damping and a Coulomb-type friction coefficient. The model is calibrated to the reference energy profiles of a previous coarse-grained model of Zhang et al., validated by isotropic compression against experimental compressibility data, and then interrogated via parametric studies of damping, temperature control, and friction. The central conclusion is that inter-particle friction is the primary driver of shear strength and structural integrity, without which clay assemblies exhibit unphysical fluid-like yielding.

Significance. If the central claims are correct, the framework would be a useful step toward incorporating non-conservative contact physics into clay CGMD, and the qualitative trends (damping suppresses densification, thermostat prevents kinetic trapping) are plausible and of practical value. Strengths include the stepwise fitting to reference MD data, the use of a common initial configuration across the damping study, and the attempt to relate simulated compressibility to experimental measurements. The paper also explicitly acknowledges some limitations, such as the small magnitude of edge-to-edge attraction and the stochastic orientation of domains in finite systems. However, the main friction result is not an independent finding: it is built into the force law and is tested in a cohesionless assembly. The validation lacks error bars or replicate configurations and rests on a broad experimental window, and the temperature-control comparison is confounded by a different friction coefficient. These issues weaken the quantitative and causal conclusions as currently stated.

major comments (4)
  1. [§2.2 and §3.4] Setting C=0 in Eq. (1) removes all attractive interactions, so the frictionless control in Section 3.4 is a purely repulsive, cohesionless assembly. The tangential force in Eq. (3) vanishes identically when μ=0, meaning the observed fluid-like yielding of the frictionless case is a direct consequence of the constitutive model rather than an emergent physical prediction. The concluding claim that 'explicit friction is the primary driver of shear strength and structural integrity' (Section 4, conclusion 3) is therefore not established, because the comparison cannot separate the role of friction from the role of omitted cohesion. To support the claim, the authors would need to repeat the control with a nonzero attractive term (or otherwise include cohesion) while keeping μ=0, and they should show whether the reference energy curves of Zhang et al. are actually well reproduced without the C/r^6 term.
  2. [§3.1, Figure 6] The validation of the baseline model is based on a single simulation trajectory with no replicated initial configurations or error bars. The compression index Cc=2.17 is compared with a wide experimental range (1.0–2.6 for bentonites), and the 'bimodal' response is inferred from only four pressure points (1, 3, 10, and 100 atm). With this scatter and a single sample, the agreement with experiments is qualitative at best, and the claim that the baseline is 'validated' is overstated. The compressibility check is also not independent in the sense claimed, because it tests a repulsion-only potential: if the fitted Buckingham A and ρ parameters compensate for the omitted C term, the same validation would not detect that distortion.
  3. [§2.2, Table 1 and Figure 4] The calibration procedure fits face-to-face and edge-to-edge interactions separately, but the paper acknowledges that 'some attraction exists in edge-to-edge cases'. Suppressing this attraction by setting C=0 may bias the effective potential for edge contacts, which are known to be important for clay fabric and strength. The paper does not show the fitted energy curves against the reference data, only force-distance profiles, so the reader cannot assess how well the repulsion-only form reproduces the reference energy landscape, in particular whether an attractive well is present. This is a load-bearing assumption for all subsequent mechanical results and needs to be justified with quantitative fitting errors or a sensitivity analysis.
  4. [§3.3] The thermostatted versus unthermostatted comparison uses a friction coefficient of μ=1.0, whereas the baseline and all other parameter studies use μ=0.1. Since Section 3.4 shows that the friction coefficient strongly affects the mechanical response, the differences in final void ratio and order parameter between the thermostatted and unthermostatted cases cannot be attributed solely to the presence or absence of a thermostat. The unthermostatted case should be repeated with the baseline μ=0.1 to cleanly isolate the role of thermal fluctuations.
minor comments (5)
  1. [Abstract] The abstract contains typos and spacing errors, such as 'propose s', 'methodological limitations ,', and 'compression index ( Cc)', which should be corrected.
  2. [§3.2.1] The text says the stability criteria are defined in Section 2.2, but the steady-state criterion (slope below 0.01% per nanosecond) is actually introduced in Section 2.3; the cross-reference should be fixed.
  3. [§2.2 (Table 1)] The Hertzian stiffness kn is listed with units kcal/(mol·Å^3/2) in Table 1, but Eq. (2) and the LAMMPS granular package convention may imply different units; clarifying the unit system and conversion would improve reproducibility.
  4. [§2.2] The description of the 'hybrid atom style (sphere and molecular)' is vague; it would be helpful to specify how this is implemented in LAMMPS (e.g., rigid bodies composed of granular spheres and the associated fix rigid command).
  5. [§4] The claim of a 'first-of-its-kind solution' is stronger than warranted given reference [20], which already studies the absence of friction in clay CGMD; the authors should frame the contribution as a hybrid framework that explicitly unifies long-range potentials with contact mechanics rather than as the first demonstration that friction matters.

Circularity Check

1 steps flagged · score 6.0 of 10

Friction conclusion reduces to the Coulomb law in Eq. (3): setting μ=0 removes tangential resistance by construction, so the central claim is partially circular despite an external compressibility validation.

  1. self definitional [Section 2.1, Eq. (3); Section 3.4, Fig. 10]
    "F_t = − min(μ‖F_n‖, |x_{γ,t} η_n v_t|) t ... To isolate the specific contribution of friction, the response of the baseline was compared against a frictionless control mode (μ = 0), as shown in Figure 10. ... the absence of shear resistance allows platelets to slide past one another with minimal energy cost. This leads to a much lower peak stress (~ 2 MPa) followed by a noticeable strain-softening response ... This divergence provides strong evidence of the necessity for explicitly incorporating friction in CG models."

    Because Eq. (3) defines the tangential contact force as proportional to μ, setting μ=0 removes all tangential resistance by construction. The uniaxial comparison therefore demonstrates the constitutive input rather than an emergent prediction: any model using this Coulomb-style law will show that frictionless contacts slide more easily and sustain lower loads. The quantitative stress-strain curve, strain at transition, and peak-stress values are genuine simulation outputs, but the headline conclusion that friction is the primary driver of shear strength is a restatement of the force law, not an independent derivation. The conclusion is further conditional on the model's zero-attraction assumption, since no competing cohesive term is present to test against friction.

full rationale

The baseline calibration is not circular: the hybrid potential is fitted to the external Zhang et al. [13] reference energy curves, and the resulting compressibility is compared against independent experimental data on bentonite, smectite, and clayey soils. That external check supports the model's macroscopic reasonableness and lowers the overall circularity burden. The central friction claim, however, is partially circular: Section 3.4 compares μ=0.1 against μ=0, and Eq. (3) makes the tangential force vanish identically when μ=0. The observed fluid-like yielding and low peak stress are therefore the expected output of the constitutive model, not a discovered role for friction. The paper's assertion in Section 2.2 that 'the C parameter in the Buckingham potential is set to zero' is flagged as a limitation: with all long-range attraction removed, the model cannot separate the contribution of friction from the omitted cohesive mechanism, and the fitted A and ρ parameters may partially absorb the missing attractive well. This weakens the external validity of the 'friction is the primary driver' conclusion, though it is a modeling assumption rather than a circular derivation. No load-bearing self-citation chain was found: the authors' own prior works are used for context or parameter ranges, not to force the central result. Overall, one key prediction reduces by construction, giving a partial circularity score of 6.

Assumptions & free parameters 13 free parameters · 5 assumptions · 1 invented entities

The model's conservative parameters are fitted to the reference CG energy curves of Zhang et al. [13], while the dissipative parameters (mu, eta_n, x_gamma,t) are chosen by hand and varied. The compressibility comparison with experiments provides external validation, which lowers circularity, but the headline finding that friction controls shear strength is largely an expected consequence of the Coulomb friction law because setting mu=0 removes tangential resistance by construction.

free parameters (13)
  • Friction coefficient mu = 0.1 baseline; 0 in frictionless control
    Chosen from the reported range for saturated clay minerals [33-35]; central to the shear-strength claim.
  • Normal damping coefficient eta_n = 0.2 baseline; 1.0 and 2.0 in sensitivity
    Selected as sub-critical relative to estimated critical damping range 0.18-0.73; controls densification and fabric order in Section 3.2.1.
  • Tangential damping scaling factor x_gamma,t = 5.0 baseline; 1.0 and 10.0 in sensitivity
    No independent source is given; strongly affects final orientation and void ratio in Section 3.2.2.
  • Buckingham repulsion coefficient A, inner-inner = 1200 kcal/mol
    Fitted to reference face-to-face energy curves from Zhang et al. [13] in Section 2.2.
  • Buckingham repulsion coefficient A, outer-outer = 600 kcal/mol
    Fitted to reference edge-to-edge energy curves from Zhang et al. [13].
  • Buckingham repulsion coefficient A, inner-outer = 400 kcal/mol
    Fitted to reference face-to-edge energy curves from Zhang et al. [13].
  • Buckingham length parameter rho, inner-inner = 0.78 Å
    Fitted parameter from Section 2.2, Table 1.
  • Buckingham length parameter rho, outer-outer = 0.70 Å
    Fitted parameter from Section 2.2, Table 1.
  • Buckingham length parameter rho, inner-outer = 1.2 Å
    Fitted parameter from Section 2.2, Table 1.
  • Hertzian stiffness k_n, inner-inner = 0.2 kcal/(mol·Å^3/2)
    Fitted to reference face-to-face energy curves from Zhang et al. [13].
  • Hertzian stiffness k_n, outer-outer = 0.08 kcal/(mol·Å^3/2)
    Fitted to reference edge-to-edge energy curves from Zhang et al. [13].
  • Hertzian stiffness k_n, inner-outer = 1.3 kcal/(mol·Å^3/2)
    Fitted to reference face-to-edge energy curves from Zhang et al. [13].
  • Buckingham attraction coefficient C = 0 (assigned)
    Set to zero because the authors state there are no significant attractive interactions; this modeling choice removes the C/r^6 term and is not sensitivity-tested.
assumptions (5)
  • domain assumption The reference energy curves from Zhang et al. [13] are accurate and transferable targets for sodium montmorillonite interactions.
    All conservative parameters are fitted to these curves in Section 2.2; if the reference is wrong or not transferable, the baseline compressibility and subsequent results lose their foundation.
  • ad hoc to paper Omitting the attractive C/r^6 term in the Buckingham potential does not materially change the mechanical response.
    Section 2.2 sets C=0 because attractions are described as comparatively small, but no sensitivity test is performed. The model is repulsive-only in the Buckingham term, which may bias fabric and compressibility.
  • domain assumption A single random initial configuration of 1000 platelets with one set of initial velocities is representative of clay assemblies.
    Sections 2.3 and 3 use one initial configuration; the authors note directional stochasticity but do not run replicate simulations, so statistical uncertainty is absent.
  • domain assumption The viscoelastic Hertzian contact with Coulomb friction law (Eqs. 2-3) adequately represents hydration-layer friction before true steric contact.
    Invoked in Section 2.1 to justify the dual-component force field; no direct validation of the friction law against lubrication experiments or atomistic friction simulations is given.
  • domain assumption LAMMPS granular package with deform/pressure and rigid body integration correctly implements the intended thermodynamics for anisotropic platelets.
    The authors rely on LAMMPS built-in contact detection and rigid-body integration without reporting convergence tests for timestep, neighbor list, or rigid-body solver settings.
invented entities (1)
  • Inner and outer sphere interaction classes in the CG-MMT platelet
    purpose: Assign different potential parameters to platelet edge and interior regions to capture edge effects during fitting
    These are model constructs introduced in Section 2.2; they are fitted to synthetic face/edge energy curves and have no independent experimental or atomistic validation outside the paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Unveiling the Role of Friction in Coarse-Grained Clay: A Hybrid Framework Integrating Long-Range Interactions and Granular Contact Mechanics." pith.science (2026). https://pith.science/paper/HTOUQBJR

@misc{pith2026260804359,
  author       = {Pith},
  title        = {Pith review of: Unveiling the Role of Friction in Coarse-Grained Clay: A Hybrid Framework Integrating Long-Range Interactions and Granular Contact Mechanics},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HTOUQBJR}},
  note         = {Machine review of arXiv:2608.04359}
}
read the original abstract

Given the predominant role of inter-particle physicochemical forces in governing clay behavior, researchers have increasingly utilized coarse-grained molecular dynamics (CGMD) simulations. However, inter-particle friction has been historically overlooked due to methodological limitations, and the extent to which this omission influences simulation accuracy remains an unresolved question. This study proposes a novel hybrid CGMD framework explicitly coupling long-range Buckingham potential with Hertzian granular contact mechanics. A baseline model was validated via isotropic compression, where the resulting compressibility and derived compression index (Cc) aligned with macroscopic geotechnical observations. Parametric analyses revealed that viscoelastic damping of particle contacts governs structural evolution. Elevated damping suppresses densification, trapping platelets in disorganized, high-void-ratio configurations. Furthermore, evaluating the interplay with thermal fluctuations underscores the necessity of precise temperature control to prevent such unphysical kinetic trapping. Finally, uniaxial compression tests demonstrate the critical importance of inter-particle friction. Explicit friction locks sliding interfaces and sustains significantly higher loads compared to frictionless models; the latter rely solely on geometric interlocking and ultimately exhibit unphysical fluid-like yielding. By bridging atomistic potentials with contact mechanics, this framework highlights the fundamental role of the inter-particle friction and offers essential guidelines for future multi-scale simulations of clay assemblies.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

53 extracted references · 53 canonical work pages

  1. [1]

    Water, Air, & Soil Pollution, 2020

    Wang, F., et al., CTMAB-Modified Bentonite–Based PRB in Remediating Cr(VI) Contaminated Groundwater. Water, Air, & Soil Pollution, 2020. 231(1): p. 20

  2. [2]

    Zhou, and W.-M

    Cui, L.-Y ., C. Zhou, and W.-M. Ye, Modified hydraulic conductivity equations of bentonite -based materials under saturated and unsaturated conditions. Canadian Geotechnical Journal, 2025. 62: p. 1-17

  3. [3]

    Journal of Physical Chemistry C, 2021

    Brochard, L., Swelling of Montmorillonite from Molecular Simulations: Hydration Diagram and Confined Water Properties. Journal of Physical Chemistry C, 2021. 125(28): p. 15527-15543

  4. [4]

    Materials Chemistry and Physics, 2022

    Al-Zaoari, K., et al., Early stage of swelling process of dehydrated montmorillonite through molecular dynamics simulation. Materials Chemistry and Physics, 2022. 283: p. 126015

  5. [5]

    Clay Minerals, 2022

    Yuan, R., et al., Molecular dynamics modelling of Na-montmorillonite subjected to uniaxial compression and unidirectional shearing. Clay Minerals, 2022. 57(3-4): p. 241-252

  6. [6]

    Chemical Physics Letters, 2022

    Han, Z., et al., A molecular dynamics study on the structural and mechanical properties of pyrophyllite and M -Montmorillonites (M = Na, K, Ca, and Ba). Chemical Physics Letters, 2022. 803

  7. [7]

    Applied Clay Science, 2022

    Wei, P., et al., Effect of water content and structural anisotropy on tensile mechanical properties of montmorillonite using molecular dynamics. Applied Clay Science, 2022. 228

  8. [8]

    Coordination Chemistry Reviews, 2025

    Zhou, A., et al., Molecular modeling of clay minerals: A thirty-year journey and future perspectives. Coordination Chemistry Reviews, 2025. 526

Show all 53 references
  1. [9]

    The Journal of Chemical Physics, 2013

    Noid, W.G., Perspective: Coarse-grained models for biomolecular systems. The Journal of Chemical Physics, 2013. 139(9): p. 090901

  2. [10]

    Whittle, and R.J.M

    Ebrahimi, D., A.J. Whittle, and R.J.M. Pellenq, Mesoscale properties of clay aggregates from potential of mean force representation of interactions between nanoplatelets. Journal of Chemical Physics, 2014. 140(15)

  3. [11]

    Whittle, and R.- M

    Ebrahimi, D., A. Whittle, and R.- M. Pellenq, Effect of Polydispersity of Clay Platelets on the Aggregation And Mechanical Properties of Clay at the Mesoscale. Clays and Clay Minerals, 2016. 64(4): p. 425-437

  4. [12]

    Pellenq, and A.J

    Ebrahimi, D., R.J.-M. Pellenq, and A.J. Whittle, Mesoscale simulation of clay aggregate formation and mechanical properties. Granular Matter, 2016. 18(3)

  5. [13]

    Langmuir, 2022

    Zhang, Y ., et al., A Coarse-Grained Interaction Model for Sodium Dominant Montmorillonite. Langmuir, 2022. 38(43): p. 13226-13237

  6. [14]

    Applied Clay Science, 2023

    Zhang, Y ., et al., Mechanical properties and pore network connectivity of sodium montmorillonite as predicted by a coarse -grained molecular model. Applied Clay Science, 2023. 243

  7. [15]

    Journal of Molecular Liquids, 2020

    Sun, H.-m., et al., A coarse-grained water model for mesoscale simulation of clay-water interaction. Journal of Molecular Liquids, 2020. 318. 22

  8. [16]

    Langmuir, 2022

    Ghazanfari, S., et al., A Coarse-Grained Model for the Mechanical Behavior of Na-Montmorillonite Clay. Langmuir, 2022. 38(16): p. 4859-4869

  9. [17]

    Shen, and I.C

    Zheng, X., X. Shen, and I.C. Bourg, Coarse-grained simulation of colloidal self- assembly, cation exchange, and rheology in Na/Ca smectite clay gels. J Colloid Interface Sci, 2025. 693: p. 137573

  10. [18]

    Tasche, and M.R

    Potter, T.D., J. Tasche, and M.R. Wilson, Assessing the transferability of common top- down and bottom -up coarse -grained molecular models for molecular mixtures. Physical Chemistry Chemical Physics, 2019. 21(4): p. 1912-1927

  11. [19]

    Journal of Chemical Theory and Computation, 2022

    Jin, J., et al., Bottom-up Coarse-Graining: Principles and Perspectives. Journal of Chemical Theory and Computation, 2022. 18(10): p. 5759-5791

  12. [20]

    International Journal of Geomechanics, 2024

    Bandera, S., et al., Effects of the Absence of Friction in Coarse -Grained Molecular Dynamics Simulations of Clay. International Journal of Geomechanics, 2024. 24(10)

  13. [21]

    Cement and Concrete Composites, 2024

    Tang, Z.- Q., et al., A novel mesoscale modelling method for steel fibre - reinforced concrete with the combined finite-discrete element method. Cement and Concrete Composites, 2024. 149: p. 105479

  14. [22]

    Yin, and C

    Liu, D., Z.-Y . Yin, and C. O’Sullivan, DEM Exploration of Stress Transmission and Small Strain Behavior of Rubber Sand Mixtures. Journal of Geotechnical and Geoenvironmental Engineering, 2025. 151(5): p. 04025031

  15. [23]

    Computers and Geotechnics, 2024

    Shen, Z., et al., DEM simulation of microscopic structure and macroscopic mechanical behavior of clay in oedometer and triaxial compression tests. Computers and Geotechnics, 2024. 173: p. 106544

  16. [24]

    Yin, and P.-Y

    Liu, B., Z.- Y . Yin, and P.-Y . Hicher, The compressibility of kaolin clay: a micromechanical perspective based on DEM. Canadian Geotechnical Journal,

  17. [25]

    Wang, Y . and N. Lu, A Closed -Form Equation for Soil Sorptive Potential. Journal of Geotechnical and Geoenvironmental Engineering, 2025. 151(9)

  18. [26]

    Yin, and Y .-Y

    Xu, W.-Q., Z.- Y . Yin, and Y .-Y . Zheng, Investigating silica interface rate - dependent friction behavior under dry and lubricated conditions with molecular dynamics. Acta Geotechnica, 2023. 18: p. 3543–3554

  19. [27]

    Computer Physics Communications, 2022

    Thompson, A.P., et al., LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Computer Physics Communications, 2022. 271: p. 108171

  20. [28]

    Physical Review E, 2001

    Silbert, L.E., et al., Granular flow down an inclined plane: Bagnold scaling and rheology. Physical Review E, 2001. 64(5): p. 051302

  21. [29]

    Zheng, X. and I.C. Bourg, Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse -Grained Model with Explicit Counterions. The Journal of Physical Chemistry C, 2026. 130(10): p. 3990-4004

  22. [30]

    Zheng, and I.C

    Shen, X., X. Zheng, and I.C. Bourg, A coarse -grained model of clay colloidal aggregation and consolidation with explicit representation of the electrical 23 double layer. J Colloid Interface Sci, 2025. 683(Pt 1): p. 1188-1196

  23. [31]

    Polymer, 2020

    Fan, J., et al., Molecular dynamics predictions of thermomechanical properties of an epoxy thermosetting polymer. Polymer, 2020. 196

  24. [32]

    Yin, and Y .Y

    Xu, W.Q., Z.Y . Yin, and Y .Y . Zheng, FRP–soil interfacial mechanical properties with molecular dynamics simulations: Insights into friction and creep behavior. International Journal for Numerical and Analytical Methods in Geomechanics, 2023(47): p. 2951-2967

  25. [33]

    Colloids and Surfaces A: Physicochemical and Engineering Aspects, 2023

    Wei, P., et al., Nanoscale friction at the quartz -quartz/kaolinite interface. Colloids and Surfaces A: Physicochemical and Engineering Aspects, 2023. 676

  26. [34]

    Computers and Geotechnics, 2024

    Wei, P., et al., Nanoscale Stick-Slip Behavior and Hydration of Hydrated Illite Clay. Computers and Geotechnics, 2024. 166

  27. [35]

    Moore, and D.A

    Morrow, C.A., D.E. Moore, and D.A. Lockner, Frictional strength of wet and dry montmorillonite. Journal of Geophysical Research: Solid Earth, 2017. 122(5): p. 3392-3409

  28. [36]

    Choong, and R

    Miari, M., K.K. Choong, and R. Jankowski, Seismic pounding between adjacent buildings: Identification of parameters, soil interaction issues and mitigation measures. Soil Dynamics and Earthquake Engineering, 2019. 121: p. 135-150

  29. [37]

    Tripathy, and T

    Baille, W., S. Tripathy, and T. Schanz, Swelling pressures and one-dimensional compressibility behaviour of bentonite at large pressures. Applied Clay Science,

  30. [38]

    Delage, and Y .J

    Marcial, D., P. Delage, and Y .J. Cui, On the high stress compression of bentonites. Canadian Geotechnical Journal, 2002. 39(4): p. 812-820

  31. [39]

    Tiwari, B. and B. Ajmera, Consolidation and swelling behavior of major clay minerals and their mixtures. Applied Clay Science, 2011. 54(3-4): p. 264-273

  32. [40]

    Sustainability, 2020

    Cheng, G., et al., Experimental Investigation of Consolidation Properties of Nano-Bentonite Mixed Clayey Soil. Sustainability, 2020. 12(2)

  33. [41]

    ACS Applied Materials & Interfaces, 2024

    Chahal, R., et al., Deep -Learning Interatomic Potential Connects Molecular Structural Ordering to the Macroscale Properties of Polyacrylonitrile. ACS Applied Materials & Interfaces, 2024. 16(28): p. 36878-36891

  34. [42]

    Macromolecules, 2026

    Huang, W., et al., Isomorphic α-Crystal and Amorphous Multiscale Structural Evolution Enables Heat -Resistant and High- Strength Semiaromatic Copolyamide Fibers. Macromolecules, 2026. 59(3): p. 1581-1590

  35. [43]

    Molecular Physics, 2019

    Zhu, H., et al., Mesoscale simulation of aggregation of imogolite nanotubes from potential of mean force interactions. Molecular Physics, 2019. 117(22): p. 3445-3455

  36. [44]

    Vescovi, and B

    Berzi, D., D. Vescovi, and B. Nadler, Shaking into order: Q -tensor/kinetic theory of vibrated non-spherical grains in a confined geometry. Journal of Fluid Mechanics, 2025. 1024: p. A32

  37. [45]

    Basconi, J.E. and M.R. Shirts, Effects of Temperature Control Algorithms on Transport Properties and Kinetics in Molecular Dynamics Simulations. Journal of Chemical Theory and Computation, 2013. 9(7): p. 2887-2899

  38. [46]

    MacAinsh, and H.-X

    Hicks, A., M. MacAinsh, and H.-X. Zhou, Removing Thermostat Distortions of 24 Protein Dynamics in Constant-Temperature Molecular Dynamics Simulations. Journal of Chemical Theory and Computation, 2021. 17(9): p. 5920-5932

  39. [47]

    Appl Clay Sci, 2021

    Zhang, L.-L., et al., Nanoscale mechanical behavior of kaolinite under uniaxial strain conditions. Appl Clay Sci, 2021. 201

  40. [48]

    Applied Clay Science, 2023

    Xiao, C., et al., Mechanical properties of defective kaolinite in tension and compression: A molecular dynamics study. Applied Clay Science, 2023. 246: p. 107164

  41. [49]

    Katti, and D.R

    Thapa, K.B., K.S. Katti, and D.R. Katti, Compression of Na–Montmorillonite Swelling Clay Interlayer Is Influenced by Fluid Polarity: A Steered Molecular Dynamics Study. Langmuir, 2020. 36(40): p. 11742-11753

  42. [50]

    Katti, and D.R

    Faisal, H.M.N., K.S. Katti, and D.R. Katti, Molecular mechanics of the swelling clay tactoid under compression, tension and shear. Applied Clay Science, 2021. 200

  43. [51]

    Marine and Petroleum Geology, 2007

    Mondol, N.H., et al., Experimental mechanical compaction of clay mineral aggregates—Changes in physical properties of mudstones during burial. Marine and Petroleum Geology, 2007. 24(5): p. 289-311

  44. [52]

    Journal of Building Engineering, 2022

    Anglade, E., et al., Physical and mechanical properties of clay –sand mixes to assess the performance of earth construction materials. Journal of Building Engineering, 2022. 51: p. 104229

  45. [53]

    Li Xue, P., et al., Compressive Modulus of Soil -Bentonite Mixtures for Cutoff Walls from CPTU Data, in Geo-Chicago 2016. 2016. p. 528-536. 25 Tables Table 1. Summary of Hybrid Potential Parameters for CG-MMT Interactions. Pair Buckingham Hertz Diameter (Å) A (kcal/mol) ρ (Å) ...

Pith tools

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