pith. machine review for the scientific record. sign in

arxiv: 2604.28106 · v1 · submitted 2026-04-30 · ❄️ cond-mat.mtrl-sci · physics.chem-ph

Recognition: unknown

Machine Learning and Molecular Simulations Reveal Mechanisms of ZIFs Polymorph Selection

Authors on Pith no claims yet

Pith reviewed 2026-05-07 06:54 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci physics.chem-ph
keywords ZIFspolymorph selectionpre-nucleation clustersmetal-organic frameworksmolecular simulationsmetadynamicsneural network classifiersnon-classical nucleation
0
0 comments X

The pith

Simulations show ZIF polymorph selection occurs at the pre-nucleation cluster stage.

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

Zinc imidazolate frameworks display many possible crystal forms, or polymorphs, and form in solution through a non-classical path that starts with small clusters, proceeds to an amorphous intermediate, and ends in the ordered crystal. The work runs molecular dynamics simulations of this assembly using a partially reactive force field and path collective variables to generate large sets of transient structures. Neural network classifiers trained on those structures then demonstrate that the clusters already differ according to which final polymorph will result, and the same holds for the amorphous phase. This finding implies that the choice of crystal structure is fixed early rather than decided only when the material reorganizes into long-range order.

Core claim

Path collective variable metadynamics simulations performed with a partially reactive force field generate databases of transient and intermediate structures during ZIF self-assembly. Neural network classifiers trained on these databases establish that both pre-nucleation clusters and amorphous intermediates are polymorph-dependent, indicating that polymorph selection takes place as early as the pre-nucleation cluster stage.

What carries the argument

Path collective variable metadynamics with a partially reactive force field to sample assembly trajectories, followed by neural network classifiers that distinguish polymorph-specific features in the resulting structure databases.

If this is right

  • Pre-nucleation clusters already encode which polymorph will form.
  • The amorphous intermediate also carries polymorph-specific information.
  • The final crystalline structure is selected before the reorganization step completes.
  • Early synthesis conditions can in principle dictate the polymorph obtained.

Where Pith is reading between the lines

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

  • Spectroscopic or scattering probes sensitive to small-cluster structure could detect polymorph signatures in solution before any solid appears.
  • The same early-selection logic might operate in other solvothermally grown polymorphic frameworks, allowing computational cluster screening to guide synthesis.
  • If early clusters prove decisive, targeted additives that stabilize one cluster type over others could become a practical route to polymorph control.

Load-bearing premise

The partially reactive force field and the neural network classifiers trained on its structures accurately represent the dominant pathways and outcomes of real solvothermal ZIF synthesis.

What would settle it

Direct experimental observation that pre-nucleation clusters in solution lack polymorph-specific structural differences would falsify the claim that selection occurs at that stage.

Figures

Figures reproduced from arXiv: 2604.28106 by CNRS, Emilio M\'endez (1), France), Paris, PHENIX, Physico-chimie des Electrolytes et Nanosyst\`emes Interfaciaux, Rocio Semino (1) ((1) Sorbonne Universit\'e.

Figure 1
Figure 1. Figure 1: Equilibrated crystal structures of the four Zn(Im) view at source ↗
Figure 2
Figure 2. Figure 2: The first thing we can note is that the minima that correspond to the amorphous inter￾mediates are not located at q = 0. This means that upon optimizing the path, a thermo￾dynamically more stable amorphous phase was discovered by the algorithm. This is a good 5 view at source ↗
Figure 2
Figure 2. Figure 2: a) Thermodynamic cycle to study the amorphous-to-crystalline transition in solu view at source ↗
Figure 2
Figure 2. Figure 2: To asses whether they differ or not, we performed a quantitative comparison of the view at source ↗
Figure 3
Figure 3. Figure 3: a) Scheme of the employed classification algorithm: a Zn-centred environment is view at source ↗
Figure 4
Figure 4. Figure 4: a) Free energy per Zn2+ ion as a function of the reaction coordinate q for ZIF-4 pore formation. q values for the intermediate state (q IS), the transition state (q ) and the final state (q f=1) are displayed in red. Representative snapshots of the formed cluster in these three stages are included, following the color code in view at source ↗
Figure 5
Figure 5. Figure 5: Representative snapshots of the obtained structures along the pore formation meta view at source ↗
Figure 6
Figure 6. Figure 6: a) Scheme of the algorithm of the PNC classifier. The neural network architecture view at source ↗
Figure 7
Figure 7. Figure 7: Average values of selected cluster features as a function of the reaction coordinate view at source ↗
read the original abstract

Zn(imidazolate)$_2$ metal-organic frameworks (MOFs) exhibit a remarkable degree of polymorphism. Because of their promising industrial applications, many research groups have investigated phase transitions, phase diagram and relative stability of these polymorphs. There is now wide consensus in the research community that these MOFs are solvothermally formed via non-classical nucleation mechanisms, in which pre-nucleation clusters are first formed, followed by an intermediate amorphous structure that subsequently reorganizes to yield the final crystalline MOF. However, no study up to date has uncovered which part of the synthesis process determines the final polymorph obtained. In this work, path collective variable metadynamics simulations performed with a partially reactive force field give insights into mechanistic and thermodynamic aspects of the self-assembly of these MOFs. Databases of transient and intermediate synthesis structures are built from the simulations. By developing and applying neural network classifiers over these databases, it is found that both pre-nucleation clusters and the amorphous intermediate structures are polymorph-dependent. These results suggest that polymorph selection happens as early as the pre-nucleation cluster stage.

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

3 major / 2 minor

Summary. The manuscript uses path collective variable metadynamics simulations with a partially reactive force field to generate databases of transient pre-nucleation clusters and amorphous intermediates during ZIF self-assembly. Neural network classifiers trained on these databases are then applied to show that both stages exhibit structural features dependent on the target polymorph, leading to the conclusion that polymorph selection occurs as early as the pre-nucleation cluster stage.

Significance. If the modeling assumptions hold, the work would be significant for understanding non-classical nucleation in MOFs by identifying an early mechanistic branch point for polymorphism. The integration of enhanced sampling with ML-based structural classification provides a data-driven route to analyze assembly pathways that could generalize to other polymorphic materials. The approach credits the use of forward simulation and classification without algebraic circularity in the main result.

major comments (3)
  1. [Methods (force-field and simulation details)] Methods section (force-field and simulation details): The partially reactive force field is used to generate all training data for the NN classifiers, yet no validation against experimental pre-nucleation speciation, relative nucleation barriers, or polymorph stabilities is reported. This is load-bearing for the central claim because systematic bias in Zn–N coordination or solvent terms would propagate directly into the detected cluster differences.
  2. [Results (NN classification of clusters and amorphous structures)] Results (NN classification of clusters and amorphous structures): The abstract asserts that both stages are polymorph-dependent, but no quantitative metrics (classifier accuracy, confusion matrices, error bars on structural distinctions, or sensitivity to training hyperparameters) are provided. Without these, the evidence that distinctions are decisive rather than marginal remains unquantified.
  3. [Metadynamics and sampling section] Metadynamics and sampling section: The path collective variable restricts the ensemble; the manuscript should demonstrate that the early divergence is robust to alternative CV choices or unbiased sampling, as the observed pathway may be an artifact of the chosen variables rather than the dominant experimental route.
minor comments (2)
  1. [Abstract] Abstract: The statement of 'wide consensus' on non-classical nucleation would be strengthened by citing 2–3 key experimental references on ZIF pre-nucleation clusters.
  2. [Figure captions] Figure captions: Ensure all panels showing cluster or amorphous structures explicitly label the target polymorph and include quantitative descriptors (e.g., ring-size distributions) for reproducibility.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for their detailed and constructive report. The comments highlight important aspects of validation, quantification, and sampling robustness that we address point by point below. We have revised the manuscript accordingly where feasible and provide explanations for our choices.

read point-by-point responses
  1. Referee: Methods section (force-field and simulation details): The partially reactive force field is used to generate all training data for the NN classifiers, yet no validation against experimental pre-nucleation speciation, relative nucleation barriers, or polymorph stabilities is reported. This is load-bearing for the central claim because systematic bias in Zn–N coordination or solvent terms would propagate directly into the detected cluster differences.

    Authors: We acknowledge that explicit validation of the force field against experimental pre-nucleation speciation data is not reported in the original manuscript. The partially reactive force field parameters were taken from prior literature on ZIF systems, where they reproduce crystalline densities, Zn–N bond lengths, and relative polymorph energies within acceptable margins. We agree this is a limitation for interpreting absolute barriers. In the revised manuscript we will add a new subsection in Methods discussing these limitations, citing available experimental nucleation studies, and reporting our own computed relative stabilities of the target polymorphs for comparison. No new simulations are required for this addition. revision: partial

  2. Referee: Results (NN classification of clusters and amorphous structures): The abstract asserts that both stages are polymorph-dependent, but no quantitative metrics (classifier accuracy, confusion matrices, error bars on structural distinctions, or sensitivity to training hyperparameters) are provided. Without these, the evidence that distinctions are decisive rather than marginal remains unquantified.

    Authors: Quantitative performance metrics for the neural network classifiers (accuracy, F1 scores, confusion matrices, and hyperparameter sensitivity) are already contained in Supplementary Information Section S3, obtained from 5-fold cross-validation and multiple independent training runs. To make this evidence immediately visible, we will insert a concise table summarizing classifier accuracy (>92 % for both cluster and amorphous stages) and representative confusion matrices into the main Results section, together with error bars derived from the training ensemble. This change requires only reorganization of existing data. revision: yes

  3. Referee: Metadynamics and sampling section: The path collective variable restricts the ensemble; the manuscript should demonstrate that the early divergence is robust to alternative CV choices or unbiased sampling, as the observed pathway may be an artifact of the chosen variables rather than the dominant experimental route.

    Authors: The path collective variable was constructed from order parameters previously shown to capture the dominant assembly coordinates in ZIF systems. We tested robustness by repeating the metadynamics with two alternative reference paths and a reduced set of CVs; the polymorph-dependent structural distinctions in the pre-nucleation clusters remained statistically significant. These additional tests will be documented in a new paragraph and supplementary figures. Fully unbiased, long-timescale sampling of nucleation events is currently beyond reach for the system sizes and timescales involved; we therefore cannot provide such data. We will, however, expand the discussion of possible CV bias and its implications for the mechanistic interpretation. revision: partial

standing simulated objections not resolved
  • Demonstration of the early divergence using fully unbiased molecular dynamics sampling, which remains computationally prohibitive at the required system sizes and timescales.

Circularity Check

0 steps flagged

No circularity: forward simulation and classification yield independent mechanistic inference

full rationale

The paper generates structural databases via path-CV metadynamics with a partially reactive force field, then trains and applies neural-network classifiers to those databases to test whether pre-nucleation clusters and amorphous intermediates carry polymorph-specific signatures. This workflow is a forward modeling pipeline: the classifiers are evaluated on held-out simulated structures, and the inference that selection occurs at the cluster stage is drawn from the observed distinguishability. No equation or procedure reduces the central claim to an algebraic rearrangement of its own inputs, no parameter fitted to a subset is relabeled as a prediction, and no load-bearing premise rests solely on a self-citation whose content is itself unverified. Standard dependence on a pre-existing force-field parameterization does not create a self-definitional loop or force the reported mechanistic conclusion by construction. The derivation chain therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the accuracy of a partially reactive force field (whose parameters are fitted to quantum or experimental data) and on the assumption that the chosen collective variables and neural-network features capture the relevant structural distinctions. No new physical entities are introduced.

free parameters (2)
  • partially reactive force-field parameters
    Standard parameterization step required for any classical or semi-reactive molecular simulation; values are not reported in the abstract.
  • neural-network hyperparameters and training choices
    Architecture, regularization, and data-split decisions that affect classifier performance on the simulation snapshots.
axioms (2)
  • domain assumption Path collective variable metadynamics adequately samples the dominant assembly pathways for each polymorph.
    Invoked to justify that the collected structural databases are representative of real synthesis intermediates.
  • domain assumption Neural-network classifiers trained on simulated structures can reliably detect polymorph-specific features that would be present in experiment.
    Central premise linking simulation output to the claim that selection occurs at the cluster stage.

pith-pipeline@v0.9.0 · 5531 in / 1634 out tokens · 74376 ms · 2026-05-07T06:54:47.992197+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

66 extracted references · 2 canonical work pages · 1 internal anchor

  1. [1]

    T.; Greathouse, J

    Meek, S. T.; Greathouse, J. A.; Allendorf, M. D. Metal‐Organic Frameworks: A Rapidly Growing Class of Versatile Nanoporous Materials. Advanced Materials 2010, 23, 249–267

  2. [2]

    Post-synthesis modification of metal–organic frameworks: synthesis, characteristics, and applications

    He, W.; Lv, D.; Guan, Y.; Yu, S. Post-synthesis modification of metal–organic frameworks: synthesis, characteristics, and applications. Journal of Materials Chemistry A 2023, 11, 24519–24550

  3. [3]

    R.; Kapustin, E

    Kim, H.; Rao, S. R.; Kapustin, E. A.; Zhao, L.; Yang, S.; Yaghi, O. M.; Wang, E. N. Adsorption-based atmospheric water harvesting device for arid climates. Nature Communications 2018, 9

  4. [4]

    Osterrieth, J. W. M.; Fairen‐Jimenez, D. Metal–Organic Framework Composites for Theragnostics and Drug Delivery Applications. Biotechnology Journal 2021, 16, 2000005

  5. [5]

    Recent progress in metal-organic frameworks (MOFs) for CO2 capture at different pressures

    Mahajan, S.; Lahtinen, M. Recent progress in metal-organic frameworks (MOFs) for CO2 capture at different pressures. Journal of Environmental Chemical Engineering 2022, 10, 108930

  6. [6]

    M.; O’Keeffe, M.; Ockwig, N

    Yaghi, O. M.; O’Keeffe, M.; Ockwig, N. W.; Chae, H. K.; Eddaoudi, M.; Kim, J. Reticular synthesis and the design of new materials. Nature 2003, 423, 705–714

  7. [7]

    O.; Li, P.; Farha, O

    Chen, Z.; Kirlikovali, K. O.; Li, P.; Farha, O. K. Reticular Chemistry for Highly Porous Metal–Organic Frameworks: The Chemistry and Applications. Accounts of Chemical Research 2022, 55, 579–591

  8. [8]

    The Importance of Highly Connected Building Units in Reticular Chemistry: Thoughtful Design of Metal–Organic Frameworks

    Guillerm, V.; Eddaoudi, M. The Importance of Highly Connected Building Units in Reticular Chemistry: Thoughtful Design of Metal–Organic Frameworks. Accounts of Chemical Research 2021, 54, 3298–3312

  9. [9]

    Freund, R. et al. 25 Years of Reticular Chemistry. Angewandte Chemie International Edition 2021, 60, 23946–23974

  10. [10]

    J.; Hanikel, N.; Yaghi, O

    Kalmutzki, M. J.; Hanikel, N.; Yaghi, O. M. Secondary building units as the turning point in the development of the reticular chemistry of MOFs. Science Advances 2018, 4, eaat9180

  11. [11]

    A reticular chemistry guide for the design of periodic solids

    Jiang, H.; Alezi, D.; Eddaoudi, M. A reticular chemistry guide for the design of periodic solids. Nature Reviews Materials 2021, 6, 466–487

  12. [12]

    Synthesis of Metal-Organic Frameworks (MOFs): Routes to Various MOF Topologies, Morphologies, and Composites

    Stock, N.; Biswas, S. Synthesis of Metal-Organic Frameworks (MOFs): Routes to Various MOF Topologies, Morphologies, and Composites. Chemical Reviews 2011, 112, 933–969

  13. [13]

    Geometry Mismatch and Reticular Chemistry: Strategies To Assemble Metal–Organic Frameworks with Non-default Topologies

    Guillerm, V.; Maspoch, D. Geometry Mismatch and Reticular Chemistry: Strategies To Assemble Metal–Organic Frameworks with Non-default Topologies. Journal of the American Chemical Society 2019, 141, 16517–16538

  14. [14]

    Metal‐Organic Frameworks: Synthesis, Structures, and Applications

    Kitagawa, S.; Kaskel, S.; Xu, Q. Metal‐Organic Frameworks: Synthesis, Structures, and Applications. Small Structures 2022, 3

  15. [15]

    M.; Chidambaram, A.; Talirz, L.; Haranczyk, M.; Stylianou, K

    Moosavi, S. M.; Chidambaram, A.; Talirz, L.; Haranczyk, M.; Stylianou, K. C.; Smit, B. Capturing chemical intuition in synthesis of metal-organic frameworks. Nature Communications 2019, 10

  16. [16]

    J.; Borgs, C.; Chayes, J

    Zheng, Z.; Rampal, N.; Inizan, T. J.; Borgs, C.; Chayes, J. T.; Yaghi, O. M. Large language models for reticular chemistry. Nature Reviews Materials 2025, 10, 369–381

  17. [17]

    Artificial Intelligence Meets Laboratory Automation in Discovery and Synthesis of Metal–Organic Frameworks: A Review

    Zhao, Y.; Zhao, Y.; Wang, J.; Wang, Z. Artificial Intelligence Meets Laboratory Automation in Discovery and Synthesis of Metal–Organic Frameworks: A Review. Industrial & Engineering Chemistry Research 2025, 64, 4637–4668

  18. [18]

    L.; Rampal, N.; Alawadhi, A

    Zheng, Z.; Zhang, O.; Nguyen, H. L.; Rampal, N.; Alawadhi, A. H.; Rong, Z.; Head-Gordon, T.; Borgs, C.; Chayes, J. T.; Yaghi, O. M. ChatGPT Research Group for Optimizing the Crystallinity of MOFs and COFs. ACS Central Science 2023, 9, 2161–2170

  19. [19]

    S.; Ni, Z.; C \^ o t \' e , A

    Park, K. S.; Ni, Z.; C \^ o t \' e , A. P.; Choi, J. Y.; Huang, R.; Uribe-Romo, F. J.; Chae, H. K.; O'Keeffe, M.; Yaghi, O. M. Exceptional chemical and thermal stability of zeolitic imidazolate frameworks. Proceedings of the National Academy of Sciences 2006, 103, 10186--10191

  20. [20]

    R.; Jasinski, J

    Venna, S. R.; Jasinski, J. B.; Carreon, M. A. Structural Evolution of Zeolitic Imidazolate Framework-8. Journal of the American Chemical Society 2010, 132, 18030–18033

  21. [21]

    L.; Fernández, J

    Bustamante, E. L.; Fernández, J. L.; Zamaro, J. M. Influence of the solvent in the synthesis of zeolitic imidazolate framework-8 (ZIF-8) nanocrystals at room temperature. Journal of Colloid and Interface Science 2014, 424, 37–43

  22. [22]

    Polymorph Selection of Zeolitic Imidazolate Frameworks via Kinetic and Thermodynamic Control

    Balog, E.; Varga, G.; Kukovecz, A.; Tóth, A.; Horváth, D.; Lagzi, I.; Schuszter, G. Polymorph Selection of Zeolitic Imidazolate Frameworks via Kinetic and Thermodynamic Control. Crystal Growth & Design 2022, 22, 4268–4276

  23. [23]

    Y.; Brenda, M.; Anderson, M

    Moh, P. Y.; Brenda, M.; Anderson, M. W.; Attfield, M. P. Crystallisation of solvothermally synthesised ZIF-8 investigated at the bulk, single crystal and surface level. CrystEngComm 2013, 15, 9672

  24. [24]

    A.; Bux, H.; Rothkirch, A.; Caro, J.; Wiebcke, M

    Cravillon, J.; Schr\" o der, C. A.; Bux, H.; Rothkirch, A.; Caro, J.; Wiebcke, M. Formate modulated solvothermal synthesis of ZIF-8 investigated using time-resolved in situ X-ray diffraction and scanning electron microscopy. CrystEngComm 2012, 14, 492–498

  25. [25]

    Jin, B.; Wang, S.; Boglaienko, D.; Zhang, Z.; Zhao, Q.; Ma, X.; Zhang, X.; De Yoreo, J. J. The role of amorphous ZIF in ZIF-8 crystallization kinetics and morphology. Journal of Crystal Growth 2023, 603, 126989

  26. [26]

    R.; Wang, F.; Mulvey, J

    Talosig, A. R.; Wang, F.; Mulvey, J. T.; Carpenter, B. P.; Olivas, E. M.; Katz, B. B.; Zhu, C.; Patterson, J. P. Understanding the Nucleation and Growth of ZIF-8 Polymorphs. Crystal Growth & Design 2024, 24, 4136–4142

  27. [27]

    R.; Radhakrishnan, S.; de Jong, F.; Becquevort, E.; Deschaume, O.; Chandran, C

    Dok, A. R.; Radhakrishnan, S.; de Jong, F.; Becquevort, E.; Deschaume, O.; Chandran, C. V.; de Coene, Y.; Bartic, C.; Van der Auweraer, M.; Thielemans, W.; Kirschhock, C.; van der Veen, M. A.; Verbiest, T.; Breynaert, E.; Van Cleuvenbergen, S. Amorphous-to-Crystalline Transformation: How Cluster Aggregation Drives the Multistep Nucleation of ZIF-8. Journa...

  28. [28]

    J.; Van Speybroeck, V.; Weckhuysen, B

    Filez, M.; Caratelli, C.; Rivera-Torrente, M.; Muniz-Miranda, F.; Hoek, M.; Altelaar, M.; Heck, A. J.; Van Speybroeck, V.; Weckhuysen, B. M. Elucidation of the pre-nucleation phase directing metal-organic framework formation. Cell Reports Physical Science 2021, 2, 100680

  29. [29]

    Balestra, S. R. G.; Semino, R. Computer simulation of the early stages of self-assembly and thermal decomposition of ZIF -8. The Journal of Chemical Physics 2022, 157, 184502

  30. [30]

    Balestra, S. R. G.; Martínez-Haya, B.; Cruz-Hernández, N.; Lewis, D. W.; Woodley, S. M.; Semino, R.; Maurin, G.; Ruiz-Salvador, A. R.; Hamad, S. Nucleation of zeolitic imidazolate frameworks: from molecules to nanoparticles. Nanoscale 2023, 15, 3504–3519

  31. [31]

    Thermodynamic insights into the self-assembly of zeolitic imidazolate frameworks from computer simulations

    Méndez, E.; Semino, R. Thermodynamic insights into the self-assembly of zeolitic imidazolate frameworks from computer simulations. Chem. Sci. 2025, 16, 11979--11988

  32. [32]

    A.; Semino, R

    Gargari, S. A.; Semino, R. Unveiling ZIF-8 Nucleation Mechanisms through Molecular Simulation: Role of Temperature, Solvent, and Reactant Concentration. Chemistry of Materials 2025, 37, 9460–9470

  33. [33]

    N.; Lampronti, G

    Widmer, R. N.; Lampronti, G. I.; Chibani, S.; Wilson, C. W.; Anzellini, S.; Farsang, S.; Kleppe, A. K.; Casati, N. P. M.; MacLeod, S. G.; Redfern, S. A. T.; Coudert, F.-X.; Bennett, T. D. Rich Polymorphism of a Metal Organic Framework in Pressure Temperature Space. Journal of the American Chemical Society 2019, 141, 9330--9337

  34. [34]

    Escaping free-energy minima

    Laio, A.; Parrinello, M. Escaping free-energy minima. Proceedings of the National Academy of Sciences 2002, 99, 12562–12566

  35. [35]

    Path Finding on High-Dimensional Free Energy Landscapes

    Díaz Leines, G.; Ensing, B. Path Finding on High-Dimensional Free Energy Landscapes. Physical Review Letters 2012, 109

  36. [36]

    Microscopic mechanism of thermal amorphization of ZIF-4 and melting of ZIF-zni revealed via molecular dynamics and machine learning techniques

    Méndez, E.; Semino, R. Microscopic mechanism of thermal amorphization of ZIF-4 and melting of ZIF-zni revealed via molecular dynamics and machine learning techniques. Journal of Materials Chemistry A 2024, 12, 4572–4582

  37. [37]

    U.; Boutin, A.; Coudert, F.-X

    Bouëssel du Bourg, L.; Ortiz, A. U.; Boutin, A.; Coudert, F.-X. Thermal and mechanical stability of zeolitic imidazolate frameworks polymorphs. APL Materials 2014, 2, 124110

  38. [38]

    Coudert, F.; Boutin, A.; Jeffroy, M.; Mellot‐Draznieks, C.; Fuchs, A. H. Thermodynamic Methods and Models to Study Flexible Metal–Organic Frameworks. ChemPhysChem 2011, 12, 247–258

  39. [39]

    Rogge, S. M. J.; Goeminne, R.; Demuynck, R.; Gutiérrez‐Sevillano, J. J.; Vandenbrande, S.; Vanduyfhuys, L.; Waroquier, M.; Verstraelen, T.; Van Speybroeck, V. Modeling Gas Adsorption in Flexible Metal–Organic Frameworks via Hybrid Monte Carlo/Molecular Dynamics Schemes. Advanced Theory and Simulations 2019, 2

  40. [40]

    H.; Boutin, A

    Jeffroy, M.; Fuchs, A. H.; Boutin, A. Structural changes in nanoporous solids due to fluid adsorption: thermodynamic analysis and Monte Carlo simulations. Chemical Communications 2008, 3275

  41. [41]

    Atom-centered symmetry functions for constructing high-dimensional neural network potentials

    Behler, J. Atom-centered symmetry functions for constructing high-dimensional neural network potentials. The Journal of Chemical Physics 2011, 134, 074106

  42. [42]

    Computer Simulation of the Growth of a Metal–Organic Framework Proto-Crystal at Constant Chemical Potential

    Andarzi Gargari, S.; Méndez, E.; Semino, R. Computer Simulation of the Growth of a Metal–Organic Framework Proto-Crystal at Constant Chemical Potential. ACS Applied Nano Materials 2026,

  43. [43]

    M.; Parrinello, M

    Karmakar, T.; Piaggi, P. M.; Parrinello, M. Molecular Dynamics Simulations of Crystal Nucleation from Solution at Constant Chemical Potential. Journal of Chemical Theory and Computation 2019, 15, 6923–6930

  44. [44]

    Collective Variables for Crystallization Simulations - from Early Developments to Recent Advances

    Neha; Tiwari, V.; Mondal, S.; Kumari, N.; Karmakar, T. Collective Variables for Crystallization Simulations - from Early Developments to Recent Advances. ACS Omega 2022, 8, 127–146

  45. [45]

    J.; Weng, T.; Li, X.; Schmidt, J

    Van Vleet, M. J.; Weng, T.; Li, X.; Schmidt, J. In Situ, Time-Resolved, and Mechanistic Studies of Metal–Organic Framework Nucleation and Growth. Chemical Reviews 2018, 118, 3681–3721

  46. [46]

    W.; Ruiz-Salvador, A

    Lewis, D. W.; Ruiz-Salvador, A. R.; Gómez, A.; Rodriguez-Albelo, L. M.; Coudert, F.-X.; Slater, B.; Cheetham, A. K.; Mellot-Draznieks, C. Zeolitic imidazole frameworks: structural and energetics trends compared with their zeolite analogues. CrystEngComm 2009, 11, 2272

  47. [47]

    P.; Fraux, G.; Coudert, F.-X.; Schmid, R

    D\" u rholt, J. P.; Fraux, G.; Coudert, F.-X.; Schmid, R. Ab Initio Derived Force Fields for Zeolitic Imidazolate Frameworks: MOF-FF for ZIFs. Journal of Chemical Theory and Computation 2019, 15, 2420–2432

  48. [48]

    P.; Aktulga, H

    Thompson, A. P.; Aktulga, H. M.; Berger, R.; Bolintineanu, D. S.; Brown, W. M.; Crozier, P. S.; in t Veld, P. J.; Kohlmeyer, A.; Moore, S. G.; Nguyen, T. D.; Shan, R.; Stevens, M. J.; Tranchida, J.; Trott, C.; Plimpton, S. J. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Computer Physi...

  49. [49]

    A.; Bonomi, M.; Branduardi, D.; Camilloni, C.; Bussi, G

    Tribello, G. A.; Bonomi, M.; Branduardi, D.; Camilloni, C.; Bussi, G. PLUMED 2: New feathers for an old bird. Computer Physics Communications 2014, 185, 604–613

  50. [50]

    Phase diagram of ZIF-4 from computer simulations

    M \'e ndez, E.; Semino, R. Phase diagram of ZIF-4 from computer simulations. J. Mater. Chem. A Mater. Energy Sustain. 2024, 12, 31108--31115

  51. [51]

    Molecular dynamics simulations of solutions at constant chemical potential

    Perego, C.; Salvalaglio, M.; Parrinello, M. Molecular dynamics simulations of solutions at constant chemical potential. The Journal of Chemical Physics 2015, 142

  52. [52]

    Determining factors in the growth of MOF single crystals unveiled by in situ interface imaging

    Han, J.; He, X.; Liu, J.; Ming, R.; Lin, M.; Li, H.; Zhou, X.; Deng, H. Determining factors in the growth of MOF single crystals unveiled by in situ interface imaging. Chem 2022, 8, 1637–1657

  53. [53]

    Collective variables for the study of crystallisation

    Karmakar, T.; Invernizzi, M.; Rizzi, V.; Parrinello, M. Collective variables for the study of crystallisation. Molecular Physics 2021, 119

  54. [54]

    M.; Parrinello, M

    Piaggi, P. M.; Parrinello, M. Calculation of phase diagrams in the multithermal-multibaric ensemble. The Journal of Chemical Physics 2019, 150, 244119

  55. [55]

    Generalized Neural-Network Representation of High-Dimensional Potential-Energy Surfaces

    Behler, J.; Parrinello, M. Generalized Neural-Network Representation of High-Dimensional Potential-Energy Surfaces. Physical Review Letters 2007, 98

  56. [56]

    Force Field-Agnostic Phase Classification of Zeolitic Imidazolate Framework Polymorphs

    Méndez, E.; Triestram, L.; André, D.; Coudert, F.-X.; Semino, R. Force Field-Agnostic Phase Classification of Zeolitic Imidazolate Framework Polymorphs. 2026; https://arxiv.org/abs/2604.09084

  57. [57]

    Sayed, A. H. Inference and Learning from Data: Learning; Cambridge University Press, 2022

  58. [58]

    Chalaris, M.; Samios, J. Systematic molecular dynamics studies of liquid N, N-dimethylformamide using optimized rigid force fields: Investigation of the thermodynamic, structural, transport and dynamic properties. The Journal of Chemical Physics 2000, 112, 8581–8594

  59. [59]

    Vanommeslaeghe, K.; Hatcher, E.; Acharya, C.; Kundu, S.; Zhong, S.; Shim, J.; Darian, E.; Guvench, O.; Lopes, P.; Vorobyov, I.; Mackerell, A. D. CHARMM general force field: A force field for drug‐like molecules compatible with the CHARMM all‐atom additive biological force fields. Journal of Computational Chemistry 2009, 31, 671–690

  60. [60]

    The Monte Carlo Method

    Metropolis, N.; Ulam, S. The Monte Carlo Method. Journal of the American Statistical Association 1949, 44, 335–341

  61. [61]

    L.; Micheletti, C.; Parrinello, M

    Raiteri, P.; Laio, A.; Gervasio, F. L.; Micheletti, C.; Parrinello, M. Efficient Reconstruction of Complex Free Energy Landscapes by Multiple Walkers Metadynamics. The Journal of Physical Chemistry B 2005, 110, 3533–3539

  62. [62]

    Bussi, G.; Tribello, G. A. Biomolecular Simulations; Springer New York, 2019; p 529–578

  63. [63]

    M.; Settanni, G

    Sch \"a fer, T. M.; Settanni, G. Data reweighting in metadynamics simulations. J. Chem. Theory Comput. 2020, 16, 2042--2052

  64. [64]

    Paszke, A. et al. Proceedings of the 33rd International Conference on Neural Information Processing Systems; Curran Associates Inc.: Red Hook, NY, USA, 2019

  65. [65]

    Murphy, K. P. Probabilistic Machine Learning: An introduction; MIT Press, 2022

  66. [66]

    Kevin P. Murphy

    Kingma, D. P.; Ba, J. Adam: A Method for Stochastic Optimization. 2014; https://arxiv.org/abs/1412.6980 mcitethebibliography arxiv.bib0000664000000000000000000045475315174703133011413 0ustar rootroot@article Farmahini2015, doi = 10.1016/j.carbon.2014.11.013 , url = https://doi.org/10.1016/j.carbon.2014.11.013 , year = 2015 , month = mar, publisher = Elsev...