arxiv: 2605.11604 · v1 · submitted 2026-05-12 · ⚛️ physics.chem-ph · cond-mat.mtrl-sci· physics.comp-ph

Recognition: 2 theorem links

· Lean Theorem

Gel-Chemistry-Dependent Heavy-Metal Ion Transport and Immobilization in Cementitious Nanopores: A Molecular Dynamics Study

Kai Gong, Qiyao He, Weiqiang Chen

Pith reviewed 2026-05-13 01:27 UTC · model grok-4.3

classification ⚛️ physics.chem-ph cond-mat.mtrl-sciphysics.comp-ph

keywords molecular dynamicscementitious gelsheavy-metal ionsion immobilizationnanoporesbinding strengthC-S-Hgel chemistry

0 comments

The pith

A relative total binding strength descriptor correlates with how strongly different cement gels immobilize heavy-metal ions in nanopores.

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

Classical molecular dynamics simulations track the diffusion of Pb2+, Ba2+, and Cs+ ions inside idealized nanopores formed by three cement gel types: C-S-H, aluminum-substituted C-(N)-A-S-H, and N-A-S-H. Mobility drops in all confined spaces compared with bulk water, yet the degree of slowing depends on the gel's surface chemistry, with aluminum-rich compositions showing the largest effect through accumulation near specific oxygen sites. The work introduces a single scalar called relative total binding strength that aggregates ion coordination, residence times, and density profiles, and reports that this scalar rises in step with measured immobilization across gel compositions, ion types, and pore diameters. The findings separate the binding roles of hydroxyl oxygens and calcium bridges in C-S-H from ion-exchange processes involving sodium in the aluminum-bearing gel.

Core claim

Ion mobility is substantially reduced in all gel nanopores relative to bulk solutions, but the extent and mechanism of suppression vary strongly with gel chemistry. C-(N)-A-S-H with higher Al/Si ratios exhibits the strongest retention, driven by ion accumulation around Al-linked oxygen species via an ion-exchange-like mechanism with charge-balancing Na+. C-S-H immobilizes ions primarily through surface hydroxyl oxygens and Ca-mediated linkages, whereas N-A-S-H exhibits more distributed binding environments. Pb2+ and Ba2+ exhibit broadly similar immobilization mechanisms, whereas Cs+ shows more distinct, gel-dependent interactions with silicate and aluminosilicate oxygen sites. A relative (rT

What carries the argument

The relative total binding strength (rTBS) descriptor, constructed from pore-averaged diffusivity, residence-time distributions, ion-density profiles, radial distribution functions, and coordination numbers to quantify how gel surface chemistry suppresses ion mobility.

If this is right

Higher Al/Si ratios strengthen ion retention in C-(N)-A-S-H through Na+-mediated exchange at Al-linked oxygens.
Pb2+ and Ba2+ share similar surface-binding pathways that differ from the more variable silicate-site interactions of Cs+.
The rTBS descriptor maintains its positive correlation with immobilization across the examined range of pore sizes.
N-A-S-H produces the most spatially distributed binding environments among the three gels.

Where Pith is reading between the lines

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

Formulations that increase aluminum content in cement gels could be screened for improved heavy-metal encapsulation performance using the rTBS metric.
The same descriptor might be adapted to compare ion retention in other nanoporous silicate materials beyond cement.
If the correlation survives in larger-scale or hydrated models, it could reduce the need for full trajectory simulations when ranking new gel chemistries.

Load-bearing premise

Classical molecular dynamics force fields and idealized nanopore models capture the real ion-gel surface interactions and charge-exchange processes without significant quantum or many-body effects that would change the mobility trends.

What would settle it

Direct experimental measurements of Pb2+, Ba2+, or Cs+ diffusivity inside real cement nanopores that reverse the ordering of retention strengths predicted for C-S-H versus C-(N)-A-S-H versus N-A-S-H gels.

read the original abstract

Cementitious materials are widely used for hazardous-waste encapsulation, yet the molecular mechanisms governing heavy-metal ion retention across different gel chemistries remain insufficiently resolved. Here, classical molecular dynamics simulations were employed to investigate the adsorption-controlled mobility of representative heavy-metal ions (Pb2+, Ba2+, and Cs+) within nanopores of C-S-H, C-(N)-A-S-H, and N-A-S-H gels. By combining pore-averaged diffusivity, spatially resolved diffusivity and residence-time analysis, ion-density profiles, two-dimensional adsorption maps, radial distribution functions, coordination analysis, and interfacial binding-strength descriptors, this study establishes a comparative atomistic framework linking gel surface chemistry to ion mobility suppression under nanoconfinement. Ion mobility is substantially reduced in all gel nanopores relative to bulk solutions, but the extent and mechanism of suppression vary strongly with gel chemistry. C-(N)-A-S-H with higher Al/Si ratios exhibits the strongest retention, driven by ion accumulation around Al-linked oxygen species via an ion-exchange-like mechanism with charge-balancing Na+. C-S-H immobilizes ions primarily through surface hydroxyl oxygens and Ca-mediated linkages, whereas N-A-S-H exhibits more distributed binding environments. Pb2+ and Ba2+ exhibit broadly similar immobilization mechanisms, whereas Cs+ shows more distinct, gel-dependent interactions with silicate and aluminosilicate oxygen sites. A relative total binding strength (rTBS) descriptor is introduced, showing a strong positive correlation with the extent of ion immobilization across gel types, ion species, and pore sizes examined. These results clarify gel-specific and ion-specific mechanisms controlling heavy-metal retention in idealized cementitious nanopores.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This MD study compares Pb2+, Ba2+, and Cs+ mobility across three cement gel nanopores and ties the differences to a new rTBS binding metric, but the results depend on how well the classical potentials rank the actual surface interactions.

read the letter

The main thing to know is that higher Al content in the gels strengthens ion retention, especially for Pb2+ and Ba2+, through accumulation at Al-linked oxygens and Na+ exchange, while the rTBS descriptor tracks that effect across the systems they simulated. Cs+ behaves differently and more variably with the silicate sites. The paper lays this out with consistent use of diffusivity profiles, residence times, RDFs, coordination numbers, and density maps in idealized C-S-H, C-(N)-A-S-H, and N-A-S-H pores. That side-by-side view is the clearest addition over earlier single-gel MD work on ion transport in C-S-H. The analyses are standard but applied evenly enough to make the gel-chemistry and ion-specific patterns readable. The rTBS correlation with immobilization holds in their data and gives a compact way to summarize the binding-mobility link. The soft spots sit in the methods. Classical force fields for these cations with aluminosilicate oxygens can miss polarization and charge transfer that DFT and ab initio MD often reorder, particularly for Pb2+ near Al sites. If the potentials get the relative binding strengths wrong, both the gel ranking and the rTBS-immobilization relation move. The perfect cylindrical pores also leave out defects and long-range electrostatics that real cement surfaces show matter for retention. rTBS is built from the same trajectories as the diffusivity numbers, so the correlation is partly by construction, though still useful to see quantified. This paper is aimed at materials chemists and engineers working on hazardous-waste encapsulation or nuclear-waste barriers in cement. Anyone already running MD on nanoporous gels or needing atomistic reasons for why one formulation retains ions better than another will find the comparative data worth reading. It shows straightforward engagement with the existing C-S-H simulation literature and produces concrete, falsifiable trends, so it deserves a serious referee. I would send it to peer review with requests for force-field validation against quantum benchmarks and some sensitivity checks on the pore model.

Referee Report

2 major / 2 minor

Summary. The manuscript employs classical molecular dynamics simulations to examine the adsorption-controlled mobility and immobilization of Pb2+, Ba2+, and Cs+ ions inside idealized nanopores of C-S-H, C-(N)-A-S-H, and N-A-S-H gels. It reports pore-averaged and spatially resolved diffusivities, residence times, ion-density profiles, RDFs, coordination numbers, and two-dimensional adsorption maps, then introduces a relative total binding strength (rTBS) descriptor that exhibits a strong positive correlation with the extent of ion immobilization across gel chemistries, ion species, and pore sizes.

Significance. If the rTBS-immobilization correlation proves robust, the work supplies a useful atomistic comparative framework for understanding how gel surface chemistry (particularly Al-linked oxygens and ion-exchange with Na+) modulates heavy-metal retention under nanoconfinement. This is relevant to hazardous-waste encapsulation in cementitious materials. The multi-metric analysis across three gel types and three ions is a clear strength, though the classical force-field foundation limits transferability.

major comments (2)

[Results section on rTBS descriptor] The rTBS descriptor (introduced in the results) is constructed directly from residence times, RDF peak heights, and coordination numbers extracted from the identical MD trajectories used to compute diffusivity suppression and immobilization metrics. This creates a moderate circularity risk: the reported strong correlation may partly reflect internal consistency of the simulation data rather than an independent, predictive link between binding strength and mobility.
[Methods] The methods section provides no explicit force-field parameters for the heavy-metal ions (especially Pb2+ and Ba2+ interactions with aluminosilicate oxygen sites), no validation against DFT or ab initio MD binding energies, and no error analysis or convergence checks on the reported diffusivities and residence times. Classical potentials are known to mis-rank ion-oxygen affinities relative to quantum results for Pb2+ with Al-linked oxygens; if such reordering occurs, both the gel-chemistry ranking and the rTBS correlation would shift.

minor comments (2)

[Abstract and Results] The abstract refers to 'two-dimensional adsorption maps' and 'interfacial binding-strength descriptors' without indicating which figures present them or whether quantitative error bars are shown; the main text should ensure these visualizations are clearly labeled and reproducible.
[Results] The notation for rTBS should be formalized as an explicit equation (including the weighting of residence time, RDF, and coordination contributions) rather than described only qualitatively, to facilitate independent verification.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed comments on our manuscript. These have prompted us to clarify the construction and interpretation of the rTBS descriptor and to strengthen the methods section with additional parameter details, limitations discussion, and convergence information. We address each major comment below.

read point-by-point responses

Referee: [Results section on rTBS descriptor] The rTBS descriptor (introduced in the results) is constructed directly from residence times, RDF peak heights, and coordination numbers extracted from the identical MD trajectories used to compute diffusivity suppression and immobilization metrics. This creates a moderate circularity risk: the reported strong correlation may partly reflect internal consistency of the simulation data rather than an independent, predictive link between binding strength and mobility.

Authors: We acknowledge the referee's valid concern about potential circularity. The rTBS is indeed a composite metric assembled from binding indicators (residence times, RDF peak intensities, and coordination numbers) drawn from the same trajectories that yield the diffusivity-based immobilization measures. While diffusivity suppression is a dynamic observable and the binding metrics are more structural, they are not fully independent. The observed correlation across gel types, ions, and pore sizes nevertheless provides a useful internal consistency check and a compact way to compare retention trends. In the revised manuscript we have added an explicit paragraph in the results section describing how rTBS is assembled, its relationship to the mobility metrics, and the fact that it functions as a descriptive rather than strictly independent predictor. We have also noted that future validation against independent datasets would be valuable. revision: partial
Referee: [Methods] The methods section provides no explicit force-field parameters for the heavy-metal ions (especially Pb2+ and Ba2+ interactions with aluminosilicate oxygen sites), no validation against DFT or ab initio MD binding energies, and no error analysis or convergence checks on the reported diffusivities and residence times. Classical potentials are known to mis-rank ion-oxygen affinities relative to quantum results for Pb2+ with Al-linked oxygens; if such reordering occurs, both the gel-chemistry ranking and the rTBS correlation would shift.

Authors: We agree that the original methods section was insufficiently detailed on these points. In the revised version we have inserted the full set of Lennard-Jones parameters and partial charges employed for Pb2+ and Ba2+ (taken from established literature sets compatible with the ClayFF-based cement force field), together with the precise mixing rules used for cross-interactions with aluminosilicate oxygen sites. We have added a dedicated limitations paragraph acknowledging that classical force fields can mis-rank Pb2+ affinities relative to ab initio results and that any such reordering would affect the reported gel rankings and rTBS correlation; we cite relevant quantum studies on this issue. For convergence and error analysis we have expanded the methods with simulation equilibration and production lengths, block-averaging procedures for diffusivities, and standard-error estimates on residence times and coordination numbers; error bars have been added to the corresponding figures and tables. revision: yes

Circularity Check

0 steps flagged

No circularity: rTBS-immobilization correlation is an observed empirical relationship from distinct MD-derived quantities

full rationale

The paper performs classical MD simulations and computes multiple independent observables from the trajectories: pore-averaged and spatially resolved diffusivities (for immobilization extent), residence times, RDFs, coordination numbers, ion-density profiles, and binding-strength descriptors. The rTBS is introduced as a composite descriptor from the binding metrics and then correlated with the immobilization metrics. These are physically distinct quantities (binding strength vs. mobility suppression) whose correlation is reported as a finding, not derived by definition or by fitting one to predict the other. No load-bearing self-citations, uniqueness theorems, or ansatzes are invoked; the work is self-contained simulation analysis without tautological reduction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claims rest on standard classical MD assumptions and the newly introduced rTBS metric derived directly from the trajectories.

axioms (1)

domain assumption Classical force fields and rigid water models sufficiently describe heavy-metal ion interactions with silicate and aluminosilicate surfaces under nanoconfinement.
Invoked throughout the MD setup and analysis without additional validation mentioned.

invented entities (1)

relative total binding strength (rTBS) no independent evidence
purpose: Quantitative descriptor to correlate ion immobilization across conditions
New metric introduced from simulation outputs; no independent experimental validation provided.

pith-pipeline@v0.9.0 · 5615 in / 1296 out tokens · 34298 ms · 2026-05-13T01:27:01.668155+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear
A relative total binding strength (rTBS) descriptor is introduced, showing a strong positive correlation with the extent of ion immobilization across gel types, ion species, and pore sizes examined.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages

[1]

immobile layer

> N –A–S–H > C –S–H ( Figure S3 c–f). Increasing the Al/Si ratio produces only a marginal additional reduction in diffusivity. 0 0.5 1 1.5 2 Pb in PbCl2 Ba in BaCl2 Cs in CsCl Diffusion coefficient (x10−9 m2/s) Solution species Bulk C–(N)–A–S–H (0) C–(N)–A–S–H (0.059) C–(N)–A–S–H (0.125) C–S–H N–A–S–H (a) Cation 0% 25% 50% 75% 100% Pb in PbCl2 Ba in BaCl2...

work page
[2]

Pedone, G

A. Pedone, G. Malavasi, M.C. Menziani, A.N. Cormack, U. Segre, A new self -consistent empirical interatomic potential model for oxides, silicates, and silica -based glasses, J. Phys. Chem. B, 110 (2006) 11780–11795

work page 2006
[3]

White, K

C.E. White, K. Page, N.J. Henson, J.L. Provis, In situ synchrotron X-ray pair distribution function analysis of the early stages of gel formation in metakaolin -based geopolymers, Appl. Clay Sci., 73 (2013) 17–25

work page 2013
[4]

Farrow, P

C. Farrow, P. Juhas, J. Liu, D. Bryndin, E. Božin, J. Bloch, T. Proffen, S. Billinge, PDFfit2 and PDFgui: computer programs for studying nanostructure in crystals, J. Phys. Condens. Matter, 19 (2007) 335219

work page 2007
[5]

Cygan, J.- J

R.T. Cygan, J.- J. Liang, A.G. Kalinichev, Molecular models of hydroxide, oxyhydroxide, and clay phases and the development of a general force field, J. Phys. Chem. B, 108 (2004) 1255–1266

work page 2004
[6]

Koneshan, J.C

S. Koneshan, J.C. Rasaiah, R. Lynden-Bell, S. Lee, Solvent structure, dynamics, and ion mobility in aqueous solutions at 25 °C, J. Phys. Chem. B, 102 (1998) 4193–4204

work page 1998
[7]

Smith, L.X

D.E. Smith, L.X. Dang, Computer simulations of NaCl association in polarizable water, J. Chem. Phys., 100 (1994) 3757–3766

work page 1994
[8]

Aqvist, Ion -water interaction potentials derived from free energy perturbation simulations, J

J. Aqvist, Ion -water interaction potentials derived from free energy perturbation simulations, J. Phys. Chem., 94 (1990) 8021–8024

work page 1990
[9]

Smith, L.X

D.E. Smith, L.X. Dang, Computer simulations of cesium –water clusters: Do ion –water clusters form gas‐phase clathrates?, J. Chem. Phys., 101 (1994) 7873–7881

work page 1994
[10]

Honorio, H

T. Honorio, H. Carasek, O. Cascudo, Water self -diffusion in CSH: Effect of confinement and temperature studied by molecular dynamics, Cem. Concr. Res., 155 (2022) 106775

work page 2022
[11]

W. Chen, K. Gong, Insights into ionic diffusion in C –S–H gel pore from molecular dynamics simulations: spatial distributions, energy barriers, and structural descriptor, J. Phys. Chem. B, 129 (2025) 10550–10567

work page 2025
[12]

MATLAB version: 9.14.0 (R2023a), The MathWorks Inc., Natick, Massachusetts, United States, 2023

The MathWorks Inc. MATLAB version: 9.14.0 (R2023a), The MathWorks Inc., Natick, Massachusetts, United States, 2023

work page 2023
[13]

B. Rode, C. Schwenk, T. Hofer, B. Randolf, Coordination and ligand exchange dynamics of solvated metal ions, Coord. Chem. Rev., 249 (2005) 2993–3006

work page 2005
[14]

Bhattacharjee, T.S

A. Bhattacharjee, T.S. Hofer, A.B. Pribil, B.R. Randolf, L.H.V. Lim, A.F. Lichtenberger, B.M. Rode, Revisiting the hydration of Pb (II): a QMCF MD approach, J. Phys. Chem. B, 113 (2009) 13007–13013

work page 2009
[15]

Mehta, P

P.K. Mehta, P. Monteiro, Concrete: microstructure, properties, and materials, 4th ed., McGraw- Hill Education, New York, 2014

work page 2014
[16]

Jiang, Y.-J

Z.-L. Jiang, Y.-J. Pan, J.-F. Lu, Y.-C. Wang, Pore structure characterization of cement paste by different experimental methods and its influence on permeability evaluation, Cem. Concr. Res., 159 (2022) 106892

work page 2022