REVIEW 2 major objections 6 minor 30 references
A Dataset of Equilibrium State Configurations of Adsorption in Zeolites
T0 review · 2 major / 6 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read AdsZeo records equilibrium methane configurations, not just averaged loadings, for 4,775 zeolite framework variants.
desk verdict Coordinate-resolved GCMC adsorption dataset that is genuinely new and internally validated; accept after fixing a text/figure mismatch and being clear about single-seed convergence limits. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing mechanism is a coordinate-preserving grand canonical Monte Carlo protocol: each framework realisation is equilibrated for 5,000 cycles and run for 200,000 production cycles, with one full configuration frame saved every 1,000 cycles. Each saved frame stores Cartesian coordinates of united-atom methane pseudo-atoms and mobile sodium cations together with per-frame counts and energy components, and all frames are linked by run and structure identifiers to scalar isotherm records in a processed database. This pairing of coordinates with scalar averages is what lets the dataset be checked for consistency and what makes it usable for configuration-level machine learning.
What would settle it
Rerun the state point with the largest reported block deviation, 7.935 methane per unit cell, from new random seeds and compare the per-frame distribution of methane counts and energies; if the new runs do not fall within the reported block statistics, the equilibrium-equivalence claim for that framework realisation fails.
Extended reading notes
Core claim
The central claim is that AdsZeo provides coordinate-resolved adsorption data across variations in framework topology, aluminium content and distribution, sodium cation arrangement, pressure, and methane loading. The processed release contains 12,415,000 saved production frames and 1,245,376,215 saved particle-coordinate records. The paper validates these records by showing that loadings reconstructed from saved frames agree with the simulation software's scalar loadings to a mean absolute deviation of 0.1794 methane per unit cell, that block statistics show no systematic drift, that pair-distance checks show no unphysical overlaps, and that all isotherms increase monotonically with pressure. External comparisons with experimental methane isotherms for three sodium-exchanged zeolites reproduce the overall uptake trends, supporting the claim that the stored configurations are realistic equilibrium samples.
Load-bearing premise
The force field and the single-run GCMC protocol produce true equilibrium ensembles for every one of the 4,775 framework realisations, including those with aluminium distributions and Al-O-Al environments the parameters were not specifically fitted to.
Editorial extensions
If this is right
- Machine-learning models can be trained to generate full adsorption configurations in charged zeolite pores, not just to predict scalar loadings.
- Spatial statistics and density estimation of methane and sodium cation distributions become possible across 191 topologies and a 0.1 to 100 bar pressure range.
- The saved frames can be used as equilibrium samples in benchmarks for configuration generation, provided data splits respect framework topology and run identity.
- Zero-loading frames are retained as valid equilibrium samples, preserving the low-pressure onset of adsorption for modelling.
- Because each topology has 25 aluminium-substitution realisations, the dataset supports studying how aluminium distribution and cation arrangement affect local adsorption geometry.
Reading between the lines
- The saved mobile sodium coordinates make the database a potential resource for studying how cation positions respond to methane loading, a coupling the paper notes but does not itself analyse.
- The force field is validated experimentally on only a few framework compositions, so a natural extension is targeted validation on additional aluminosilicate frameworks, especially those containing Al-O-Al environments.
- Users modelling low-loading regimes should weigh the reported block-deviation statistics, with a largest absolute deviation of 7.935 methane per unit cell, when deciding which runs are converged enough for their task.
- Coordinate-resolved datasets of this kind could be used to train conditional generative samplers that combine pressure, topology, and cation arrangement as inputs, potentially reducing the cost of future GCMC screening.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper presents AdsZeo, a coordinate-resolved dataset of grand-canonical Monte Carlo (GCMC) methane adsorption configurations in aluminium-substituted, sodium-exchanged zeolites. The authors generate 4,775 framework realisations from 191 IZA topologies, simulate each at 13 pressures between 0.1 and 100 bar at 298 K with RASPA, and store scalar adsorption records (isotherms, frame-level energy/counts) together with saved production-frame Cartesian coordinates of methane pseudo-atoms and Na+ cations in a DuckDB database. The paper describes the structure-generation workflow, the assembled force field, the GCMC protocol, the database schema, and four internal validation checks (block statistics, coordinate-derived loading consistency, pair-distance geometry, monotonicity of isotherms) plus a comparison with experimental methane isotherms for Na-ZSM-5, NaY, and NaX. The central claim is that the released coordinates represent equilibrated GCMC samples across a broad range of framework chemistry and pressure, suitable for training machine-learning models for configuration generation.
Significance. If the dataset meets its claims, it fills a genuine gap in the high-throughput adsorption literature, which has mostly archived scalar properties; the coordinate-level records would enable generative and density-estimation models for charged nanoporous hosts. The paper's internal validation is appropriately matched to the dataset claim: the coordinate-derived loadings reconstruct the RASPA scalar loadings with a mean absolute deviation of 0.18 CH4 per unit cell, no short-range coordinate artefacts are reported, all isotherms are monotonic, and the dataset and generation code are publicly available. The main risk to the equilibrium claim is the per-run convergence tail documented in Section 4.1, which is disclosed but not resolved at the level of individual runs. On balance, the resource is valuable and the weaknesses are largely addressable.
major comments (2)
- [§2.3, §4.1] The paper labels all 62,075 runs as equilibrium samples, but each run used a single random seed and the block statistics in Section 4.1 show a non-negligible tail of runs with drift: the 95th percentile absolute relative block deviation is 5.26% and the largest absolute block deviation is 7.935 CH4 per unit cell. Because the central claim is that the saved coordinates are equilibrium configurations, the release should include a per-run convergence diagnostic (e.g., block-deviation values or a pass/fail flag) and the text should state the fraction of runs satisfying a defined convergence criterion; otherwise the 'equilibrium' label is not verifiable for the outlier runs.
- [§2.2, §4.5] The force field assembled from Refs. 24-26 is applied to all 191 topologies, including Al-O-Al environments, but the external validation in Section 4.5 covers only Na-ZSM-5, NaY, and NaX. The paper should either provide additional transferability evidence for representative topologies outside these three systems or explicitly state that the 'realistic methane adsorption behaviour' claim at the end of Section 4.5 is established only for the tested systems and is a transferability assumption elsewhere.
minor comments (6)
- [§4.1, Figure 4] The body text gives example mean methane counts of 32.3, 50.5, and 60.0 CH4 for the 1, 10, and 100 bar traces, while the Figure 4 labels show 36.0, 86.8, and 120.9 CH4; these numbers must be reconciled so that the primary production-stability figure is self-consistent.
- [§2.1, §4.5] Section 2.1 states that Löwenstein's rule was not explicitly enforced, while Section 4.5 says Al substitutions were generated 'subject to the non-Loewenstein constraint'; please clarify whether this means the same thing, since the phrase 'non-Loewenstein constraint' could be read as enforcing absence of Al-O-Al rather than allowing it.
- [§4.4] There is a typo in 'near-zere methan loading'; it should read 'near-zero methane loading'.
- [§4.1] The sentence 'do not show a appreciable, systematic first-to-last-block drift' contains a grammar error ('a appreciable'); it should read 'no appreciable systematic drift'.
- [Figure 2, Table 3] Figure 2 lists 'Structures / isotherms 4,775' while Table 3 reports 1,489,800 rows in the isotherms table; please clarify that the former counts distinct isotherms while the latter counts scalar records, to avoid an apparent inconsistency.
- [§2.3] The 13 pressures are only described as 'logarithmically spaced between 0.1 and 100 bar'; for reproducibility, please list the exact pressure values or point to where they are stored in the database schema.
Circularity Check
No circularity: the dataset claim rests on an inherited force field and GCMC runs with clearly disclosed validation limits; no prediction reduces to a fit or to the dataset's own outputs.
full rationale
AdsZeo makes no parameter fit and derives no quantitative result from its own outputs. Its central claim is that the released coordinates are equilibrated GCMC samples under the stated force field. The force field is assembled from Refs. 24-26, which are prior published parameter sets by overlapping groups; the paper does not refit them, and those references carry their own external validation, so this is independent support rather than circularity. The block-stability checks (Sec. 4.1) and coordinate-vs-scalar consistency check (Sec. 4.2) are internal consistency validations of parsing and stationarity, not predictions derived from the data they validate; the paper explicitly frames Sec. 4.2 as consistency 'rather than identity between two independently stored quantities.' The experimental comparison (Sec. 4.5) is presented as an external consistency check, and the paper discloses that quantitative agreement is not expected because simulated structures are generated from ideal frameworks. No uniqueness theorem, ansatz, or fitted-parameter prediction is invoked. The paper itself discloses the main validity limits: single-seed runs (Sec. 2.3), the largest block deviation of 7.935 CH4/uc and the 4 CH4/uc cutoff for relative statistics (Sec. 4.1), and transferability assessed with only three experimental comparisons (Sec. 4.5). These affect evidence strength, not circularity. One non-circular consistency issue should be corrected: Sec. 4.1 gives example means 32.3/50.5/60.0 CH4 while Fig. 4 labels show 36.0/86.8/120.9; the figure and text need reconciliation, but this does not indicate that any claimed result reduces to its inputs.
Assumptions & free parameters
free parameters (5)
- Al substitution count per topology =
interpolated from 1 Al to about 25 percent of tetrahedral sites
- Number of framework realisations per topology =
25
- Supercell lattice-vector threshold =
all lattice vectors longer than 24 angstroms
- Production length and saving interval =
200,000 production cycles; one saved frame per 1,000 cycles
- Pressure grid =
13 log-spaced pressures from 0.1 to 100 bar
assumptions (5)
- domain assumption Framework atoms are rigid during adsorption.
- domain assumption The force field parameters from references 24-26 transfer to all 191 topologies, including Al-O-Al environments.
- domain assumption 200,000 production cycles without repeated seeds are sufficient for equilibration in every run.
- domain assumption Loewenstein's rule is not enforced, so Al-O-Al environments are allowed.
- domain assumption United-atom methane and point-charge Na+ are adequate representations for the intended machine-learning training purpose.
Cite this review
Pith. "Pith review of A Dataset of Equilibrium State Configurations of Adsorption in Zeolites." pith.science (2026). https://pith.science/paper/FNA3PYAE
@misc{pith2026260808848,
author = {Pith},
title = {Pith review of: A Dataset of Equilibrium State Configurations of Adsorption in Zeolites},
year = {2026},
howpublished = {\url{https://pith.science/paper/FNA3PYAE}},
note = {Machine review of arXiv:2608.08848}
}
abstract
Zeolites are crystalline nanoporous materials widely used in adsorption, separation, and catalytic processes. Molecular simulations are commonly used to predict adsorption properties, but most high-throughput adsorption datasets report only ensemble-averaged quantities such as loadings or isotherms, rather than the molecular configurations from which these averages are obtained. Here, we present AdsZeo, a coordinate-resolved dataset of equilibrium methane adsorption configurations in aluminium-substituted, sodium-containing zeolite frameworks. The processed release contains 4,775 framework realisations derived from 191 zeolite topologies. Each framework realisation was simulated at 13 methane pressures between 0.1 and 100 bar at 298 K using grand canonical Monte Carlo simulations, giving 62,075 production simulations in total. In addition to scalar adsorption records, the dataset stores production-frame methane pseudo-atom coordinates, mobile Na$^+$ cation coordinates, framework atomic coordinates, per-frame loading and energy statistics, and simulation metadata in a processed DuckDB database. The release contains 12,415,000 saved production-frame records and 1,245,376,215 saved particle-coordinate records. AdsZeo provides coordinate-resolved adsorption data across variations in framework topology, aluminium content and distribution, sodium cation arrangement, pressure, and methane loading, enabling reuse for adsorption analysis, spatial statistics, density estimation, and machine-learning models for molecular configuration generation in charged zeolite pores.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
Eduardo P´ erez-Botella, Susana Valencia, and Fernando Rey. Zeolites in adsorption processes: State of the art and future prospects.Chemical reviews, 122(24):17647–17695, 2022
work page 2022
-
[2]
Metal–organic framework materials as catalysts.Chemical Society Reviews, 38(5):1450–1459, 2009
JeongYong Lee, Omar K Farha, John Roberts, Karl A Scheidt, SonBinh T Nguyen, and Joseph T Hupp. Metal–organic framework materials as catalysts.Chemical Society Reviews, 38(5):1450–1459, 2009
work page 2009
-
[3]
Hao Li, Libo Li, Rui-Biao Lin, Wei Zhou, Zhangjing Zhang, Shengchang Xiang, and Banglin Chen. Porous metal-organic frameworks for gas storage and separation: Status and challenges.Energy- Chem, 1(1):100006, 2019
work page 2019
-
[4]
Tina D¨ uren, Youn-Sang Bae, and Randall Q Snurr. Using molecular simulation to characterise metal–organic frameworks for adsorption applications.Chemical Society Reviews, 38(5):1237–1247, 2009
work page 2009
-
[5]
David Dubbeldam, Sof ´ ıa Calero, Donald E Ellis, and Randall Q Snurr. RASPA: molecular simulation software for adsorption and diffusion in flexible nanoporous materials.Molecular Simulation, 42(2): 81–101, 2016
work page 2016
-
[6]
On the inner workings of monte carlo codes.Molecular Simulation, 39(14-15):1253–1292, 2013
David Dubbeldam, Ariana Torres-Knoop, and Krista S Walton. On the inner workings of monte carlo codes.Molecular Simulation, 39(14-15):1253–1292, 2013. 13
work page 2013
-
[7]
Daniel W Siderius, Harold W Hatch, and Vincent K Shen. Flat-histogram monte carlo simulation of water adsorption in metal–organic frameworks.The Journal of Physical Chemistry B, 128(19): 4830–4845, 2024
work page 2024
-
[8]
Zhenzi Yu, Dylan M Anstine, Salah Eddine Boulfelfel, Chenkai Gu, Coray M Colina, and David S Sholl. Incorporating flexibility effects into metal–organic framework adsorption simulations using different models.ACS Applied Materials & Interfaces, 13(51):61305–61315, 2021
work page 2021
Show all 30 references
-
[9]
Language models are few-shot learners.Advances in neural information processing systems, 33:1877–1901, 2020
Tom Brown, Benjamin Mann, Nick Ryder, Melanie Subbiah, Jared D Kaplan, Prafulla Dhariwal, Arvind Neelakantan, Pranav Shyam, Girish Sastry, Amanda Askell, et al. Language models are few-shot learners.Advances in neural information processing systems, 33:1877–1901, 2020
1901
-
[10]
High- resolution image synthesis with latent diffusion models
Robin Rombach, Andreas Blattmann, Dominik Lorenz, Patrick Esser, and Bj¨ orn Ommer. High- resolution image synthesis with latent diffusion models. InProceedings of the IEEE/CVF conference on computer vision and pattern recognition, pages 10684–10695, 2022
2022
-
[11]
De novo design of protein structure and function with rfdiffusion.Nature, 620(7976):1089–1100, 2023
Joseph L Watson, David Juergens, Nathaniel R Bennett, Brian L Trippe, Jason Yim, Helen E Eisenach, Woody Ahern, Andrew J Borst, Robert J Ragotte, Lukas F Milles, et al. De novo design of protein structure and function with rfdiffusion.Nature, 620(7976):1089–1100, 2023
2023
-
[12]
Boltzmann generators: Sampling equilibrium states of many-body systems with deep learning.Science, 365(6457):eaaw1147, 2019
Frank No´ e, Simon Olsson, Jonas K¨ ohler, and Hao Wu. Boltzmann generators: Sampling equilibrium states of many-body systems with deep learning.Science, 365(6457):eaaw1147, 2019
2019
-
[13]
Torsional dif- fusion for molecular conformer generation.Advances in neural information processing systems, 35: 24240–24253, 2022
Bowen Jing, Gabriele Corso, Jeffrey Chang, Regina Barzilay, and Tommi Jaakkola. Torsional dif- fusion for molecular conformer generation.Advances in neural information processing systems, 35: 24240–24253, 2022
2022
-
[14]
Equivariant flow matching.Advances in Neural Information Processing Systems, 36:59886–59910, 2023
Leon Klein, Andreas Kr¨ amer, and Frank No´ e. Equivariant flow matching.Advances in Neural Information Processing Systems, 36:59886–59910, 2023
2023
-
[15]
Transferable boltzmann generators.Advances in Neural Information Processing Systems, 37:45281–45314, 2024
Leon Klein and Frank No´ e. Transferable boltzmann generators.Advances in Neural Information Processing Systems, 37:45281–45314, 2024
2024
-
[16]
MOFX-DB: An online database of computational adsorption data for nanoporous materials.Journal of Chemical & Engineering Data, 68(2):483–498, 2023
N Scott Bobbitt, Kaihang Shi, Benjamin J Bucior, Haoyuan Chen, Nathaniel Tracy-Amoroso, Zhao Li, Yangzesheng Sun, Julia H Merlin, J Ilja Siepmann, Daniel W Siderius, et al. MOFX-DB: An online database of computational adsorption data for nanoporous materials.Journal of Chemica...
2023
-
[17]
Data-driven design of metal–organic frameworks for wet flue gas CO2 capture
Peter G Boyd, Arunraj Chidambaram, Enrique Garc ´ ıa-D ´ ıez, Christopher P Ireland, Thomas D Daff, Richard Bounds, Andrzej G ladysiak, Pascal Schouwink, Seyed Mohamad Moosavi, M Mercedes Maroto-Valer, et al. Data-driven design of metal–organic frameworks for wet flue gas CO2 ...
2019
-
[18]
Database for CO2 separation performances of MOFs based on computational materials screening.ACS applied materials & interfaces, 10(20):17257–17268, 2018
Cigdem Altintas, Gokay Avci, Hilal Daglar, Ayda Nemati Vesali Azar, Sadiye Velioglu, Ilknur Eru- car, and Seda Keskin. Database for CO2 separation performances of MOFs based on computational materials screening.ACS applied materials & interfaces, 10(20):17257–17268, 2018
2018
-
[19]
Crafted: An exploratory database of simulated adsorption isotherms of metal-organic frameworks.Scientific Data, 10(1):230, 2023
Felipe Lopes Oliveira, Conor Cleeton, Rodrigo Neumann Barros Ferreira, Binquan Luan, Amir H Farmahini, Lev Sarkisov, and Mathias Steiner. Crafted: An exploratory database of simulated adsorption isotherms of metal-organic frameworks.Scientific Data, 10(1):230, 2023
2023
-
[20]
The open DAC 2023 dataset and challenges for sorbent discovery in direct air capture, 2024
Anuroop Sriram, Sihoon Choi, Xiaohan Yu, Logan M Brabson, Abhishek Das, Zachary Ulissi, Matt Uyttendaele, Andrew J Medford, and David S Sholl. The open DAC 2023 dataset and challenges for sorbent discovery in direct air capture, 2024
2023
-
[21]
The open DAC 2025 dataset for sorbent discovery in direct air capture.arXiv preprint arXiv:2508.03162, 2025
Anuroop Sriram, Logan M Brabson, Xiaohan Yu, Sihoon Choi, Kareem Abdelmaqsoud, Elias Moubarak, Pim de Haan, Sindy L¨ owe, Johann Brehmer, John R Kitchin, et al. The open DAC 2025 dataset for sorbent discovery in direct air capture.arXiv preprint arXiv:2508.03162, 2025
2025
-
[22]
Elsevier, 2007
Christian Baerlocher, Lynne B McCusker, and David H Olson.Atlas of zeolite framework types. Elsevier, 2007
2007
-
[23]
marko-petkovic/porran: Porran, March 2025
Marko Petkovic and Koen Wortelboer. marko-petkovic/porran: Porran, March 2025. URLhttps: //doi.org/10.5281/zenodo.15050436. 14
2025 doi
-
[24]
Adsorption of carbon dioxide in non-l¨ owenstein zeolites.Chemistry of Materials, 35(13):5222–5231, 2023
Pablo Romero-Marimon, Juan Jose Gutierrez-Sevillano, and Sofia Calero. Adsorption of carbon dioxide in non-l¨ owenstein zeolites.Chemistry of Materials, 35(13):5222–5231, 2023
2023
-
[25]
Understanding the role of sodium during adsorption: a force field for alkanes in sodium-exchanged faujasites.Journal of the American Chemical Society, 126(36):11377–11386, 2004
Sof ´ ıa Calero, David Dubbeldam, Rajamani Krishna, Berend Smit, Thijs JH Vlugt, Joeri FM De- nayer, Johan A Martens, and Theo LM Maesen. Understanding the role of sodium during adsorption: a force field for alkanes in sodium-exchanged faujasites.Journal of the American Chemic...
2004
-
[26]
Transferable force field for carbon dioxide adsorption in zeolites.The Journal of Physical Chemistry C, 113(20):8814–8820, 2009
Almudena Garcia-Sanchez, Conchi O Ania, Jos´ e B Parra, David Dubbeldam, Thijs JH Vlugt, Rajamani Krishna, and Sofia Calero. Transferable force field for carbon dioxide adsorption in zeolites.The Journal of Physical Chemistry C, 113(20):8814–8820, 2009
2009
-
[27]
AdsZeo: A dataset of equilibrium state configurations of adsorption in zeolite
Marko Petkovic, Rachna Ramesh, Vlado Menkovski, and Sofia Calero. AdsZeo: A dataset of equilibrium state configurations of adsorption in zeolite. Zenodo, July 2026. URLhttps: //doi.org/10.5281/zenodo.21445386
2026 doi
-
[28]
Calorimetric heats of adsorption and adsorption isotherms
JA Dunne, M Rao, S Sircar, RJ Gorte, and AL Myers. Calorimetric heats of adsorption and adsorption isotherms. 2. O2, N2, Ar, CO2, CH4, C2H6, and SF6 on NaX, H-ZSM-5, and Na-ZSM-5 zeolites.Langmuir, 12(24):5896–5904, 1996
1996
-
[29]
Effect of cations on methane adsorption by nay, mgy, cay, sry, and bay zeolites.The Journal of Physical Chemistry, 97(49):12894–12898, 1993
Orhan Talu, Si Yang Zhang, and David T Hayhurst. Effect of cations on methane adsorption by nay, mgy, cay, sry, and bay zeolites.The Journal of Physical Chemistry, 97(49):12894–12898, 1993
1993
-
[30]
Adsorption of methane on nax zeolite in the subcritical and supercritical regions.Bull
EV Chkhaidze, AA Fomkin, VV Serpinskii, and GV Tsitsishuili. Adsorption of methane on nax zeolite in the subcritical and supercritical regions.Bull. Acad. Sci. USSR, Div. Chem. Sci.(Engl. Transl.);(United States), 34(5 PT 1), 1985. 15
1985
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.