1
|
France-Lanord A, Menon S, Lam J. Harvesting nucleating structures in nanoparticle crystallization: The example of gold, silver, and iron. J Chem Phys 2024; 161:044108. [PMID: 39046348 DOI: 10.1063/5.0200850] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2024] [Accepted: 07/10/2024] [Indexed: 07/25/2024] Open
Abstract
The thermodynamics and kinetics of nanoparticle crystallization, as opposed to bulk phases, may be influenced by surface and size effects. We investigate the importance of such factors in the crystallization process of gold, silver, and iron nanodroplets using numerical simulations in the form of molecular dynamics combined with path sampling. This modeling strategy is targeted at obtaining representative ensembles of structures located at the transition state of the crystallization process. A structural analysis of the transition state ensembles reveals that both the average size and location of the critical nucleation cluster are influenced by surface and nanoscale size effects. Furthermore, we also show that transition state structures in smaller nanodroplets exhibit a more ordered liquid phase, and differentiating between a well-ordered critical cluster and its surrounding disordered liquid phase becomes less evident. All in all, these findings demonstrate that crystallization mechanisms in nanoparticles go beyond the assumptions of classical nucleation theory.
Collapse
Affiliation(s)
- Arthur France-Lanord
- Sorbonne Université, Muséum National d'Histoire Naturelle, UMR CNRS 7590, Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, IMPMC, F-75005 Paris, France
| | - Sarath Menon
- Max-Planck-Institut für Eisenforschung GmbH, D-40237 Düsseldorf, Germany
| | - Julien Lam
- Univ. Lille, CNRS, INRA, ENSCL, UMR 8207, UMET, Unité Matériaux et Transformations, F 59000 Lille, France
- Centre d'élaboration des Matériaux et d'Etudes Structurales, CNRS (UPR 8011), 29 rue Jeanne Marvig, 31055 Toulouse Cedex 4, France
| |
Collapse
|
2
|
Al-Awad AS, Batet L, Rives R, Sedano L. Stochastic computer experiments of the thermodynamic irreversibility of bulk nanobubbles in supersaturated and weak gas-liquid solutions. J Chem Phys 2024; 161:024503. [PMID: 38984961 DOI: 10.1063/5.0204665] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2024] [Accepted: 06/24/2024] [Indexed: 07/11/2024] Open
Abstract
Spontaneous gas-bubble nucleation in weak gas-liquid solutions has been a challenging topic in theory, experimentation, and computer simulations. In analogy with recent advances in crystallization and droplet formation studies, the diffusive-shielding stabilization and thermodynamic irreversibility of bulk nanobubble (bNB) mechanisms are revisited and deployed to characterize nucleation processes in a stochastic framework of computer experiments using the large-scale atomic/molecular massively parallel simulator code. Theoretical bases, assumptions, and limitations underlying the irreversibility hypothesis of bNBs, and their computational counterparts, are extensively described and illustrated. In essence, it is established that the irreversibility hypothesis can be numerically investigated by converging the system volume (due to the finiteness of interatomic forces) and the initial dissolved-gas concentration in the solution (due to the single-bNB limitation). Helium nucleation in liquid Pb17Li alloy is selected as a representative case study, where it exhibits typical characteristics of noble-gas/liquid-metal systems. The proposed framework lays down the bases on which the stability of gas-bNBs in weak and supersaturated gas-liquid solutions can be inferred and explained from a novel perspective. In essence, it stochastically marches toward a unique irreversible state along out-of-equilibrium nucleation/growth trajectories. Moreover, it does not attempt to characterize the interface or any interface-related properties, neither theoretically nor computationally. It was concluded that bNBs of a few tens of He-atoms are irreversible when dissolved-He concentrations in the weak gas-liquid solution are at least ∼50 and ∼105 mol m-3 at 600 and 1000 K (and ∼80 MPa), respectively, whereas classical molecular dynamics -estimated solubilities are at least two orders of magnitude smaller.
Collapse
Affiliation(s)
- Abdulrahman S Al-Awad
- Department of Physics, Universitat Politècnica de Catalunya-BarcelonaTech (UPC), Barcelona 08028, Spain
| | - Lluis Batet
- Department of Physics, Universitat Politècnica de Catalunya-BarcelonaTech (UPC), Barcelona 08028, Spain
| | - Ronny Rives
- Department of Physics, Universitat Politècnica de Catalunya-BarcelonaTech (UPC), Barcelona 08028, Spain
| | - Luis Sedano
- Department of Physics, Universitat Politècnica de Catalunya-BarcelonaTech (UPC), Barcelona 08028, Spain
- Instituto de Ciencia de Materiales de Barcelona (ICMAB/CSIC), Bellaterra 08193, Spain
| |
Collapse
|
3
|
Domingues TS, Hussain S, Haji-Akbari A. Divergence among Local Structure, Dynamics, and Nucleation Outcome in Heterogeneous Nucleation of Close-Packed Crystals. J Phys Chem Lett 2024; 15:1279-1287. [PMID: 38284350 DOI: 10.1021/acs.jpclett.3c03561] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2024]
Abstract
Heterogeneous crystal nucleation is the dominant mechanism of crystallization in most systems, yet its underlying physics remains an enigma. While emergent interfacial crystalline order precedes heterogeneous nucleation, its importance in the nucleation mechanism is unclear. Here, we use path sampling simulations of two model systems to demonstrate that crystalline order in its traditional sense is not predictive of the outcome of the heterogeneous nucleation of close-packed crystals. Consequently, structure-based collective variables (CVs) that reliably describe homogeneous nucleation can be poor descriptors of heterogeneous nucleation. This divergence between structure and nucleation outcome is accompanied by an intriguing dynamical anomaly, wherein low-coordinated crystalline particles outpace their liquid-like counterparts. We use committor analysis, high-throughput screening, and machine learning to devise CV optimization strategies and present suitable structural heuristics within the metastable fluid for CV prescreening. Employing such optimized CVs is pivotal for properly characterizing the mechanism of heterogeneous nucleation in metallic and colloidal systems.
Collapse
Affiliation(s)
- Tiago S Domingues
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06511, United States
| | - Sarwar Hussain
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06511, United States
| | - Amir Haji-Akbari
- Department of Chemical and Environmental Engineering, Yale University, New Haven, Connecticut 06511, United States
| |
Collapse
|
4
|
Naleem N, Abreu CRA, Warmuz K, Tong M, Kirmizialtin S, Tuckerman ME. An exploration of machine learning models for the determination of reaction coordinates associated with conformational transitions. J Chem Phys 2023; 159:034102. [PMID: 37458344 DOI: 10.1063/5.0147597] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Accepted: 06/23/2023] [Indexed: 07/20/2023] Open
Abstract
Determining collective variables (CVs) for conformational transitions is crucial to understanding their dynamics and targeting them in enhanced sampling simulations. Often, CVs are proposed based on intuition or prior knowledge of a system. However, the problem of systematically determining a proper reaction coordinate (RC) for a specific process in terms of a set of putative CVs can be achieved using committor analysis (CA). Identifying essential degrees of freedom that govern such transitions using CA remains elusive because of the high dimensionality of the conformational space. Various schemes exist to leverage the power of machine learning (ML) to extract an RC from CA. Here, we extend these studies and compare the ability of 17 different ML schemes to identify accurate RCs associated with conformational transitions. We tested these methods on an alanine dipeptide in vacuum and on a sarcosine dipeptoid in an implicit solvent. Our comparison revealed that the light gradient boosting machine method outperforms other methods. In order to extract key features from the models, we employed Shapley Additive exPlanations analysis and compared its interpretation with the "feature importance" approach. For the alanine dipeptide, our methodology identifies ϕ and θ dihedrals as essential degrees of freedom in the C7ax to C7eq transition. For the sarcosine dipeptoid system, the dihedrals ψ and ω are the most important for the cisαD to transαD transition. We further argue that analysis of the full dynamical pathway, and not just endpoint states, is essential for identifying key degrees of freedom governing transitions.
Collapse
Affiliation(s)
- Nawavi Naleem
- Chemistry Program, Science Division, New York University, Abu Dhabi, UAE
| | - Charlles R A Abreu
- Chemical Engineering Department, Escola de Química, Universidade Federal do Rio de Janeiro, 21941-909 Rio de Janeiro, RJ, Brazil
| | - Krzysztof Warmuz
- Computer Science Program, Science Division, New York University, Abu Dhabi, UAE
| | - Muchen Tong
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
| | - Serdal Kirmizialtin
- Chemistry Program, Science Division, New York University, Abu Dhabi, UAE
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
- Center for Smart Engineering Materials, New York University, Abu Dhabi, UAE
| | - Mark E Tuckerman
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Rd. North, Shanghai 200062, China
- Simons Center for Computational Physical Chemistry at New York University, New York, New York 10003, USA
| |
Collapse
|
5
|
Rogal J, Díaz Leines G. Controlling crystallization: what liquid structure and dynamics reveal about crystal nucleation mechanisms. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2023; 381:20220249. [PMID: 37211029 DOI: 10.1098/rsta.2022.0249] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/03/2022] [Accepted: 12/06/2022] [Indexed: 05/23/2023]
Abstract
Over recent years, molecular simulations have provided invaluable insights into the microscopic processes governing the initial stages of crystal nucleation and growth. A key aspect that has been observed in many different systems is the formation of precursors in the supercooled liquid that precedes the emergence of crystalline nuclei. The structural and dynamical properties of these precursors determine to a large extent the nucleation probability as well as the formation of specific polymorphs. This novel microscopic view on nucleation mechanisms has further implications for our understanding of the nucleating ability and polymorph selectivity of nucleating agents, as these appear to be strongly linked to their ability in modifying structural and dynamical characteristics of the supercooled liquid, namely liquid heterogeneity. In this perspective, we highlight recent progress in exploring the connection between liquid heterogeneity and crystallization, including the effects of templates, and the potential impact for controlling crystallization processes. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.
Collapse
Affiliation(s)
- Jutta Rogal
- Department of Chemistry, New York University, New York, NY 10003, USA
- Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany
| | - Grisell Díaz Leines
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, UK
| |
Collapse
|
6
|
Hall SW, Díaz Leines G, Sarupria S, Rogal J. Practical guide to replica exchange transition interface sampling and forward flux sampling. J Chem Phys 2022; 156:200901. [DOI: 10.1063/5.0080053] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Path sampling approaches have become invaluable tools to explore the mechanisms and dynamics of the so-called rare events that are characterized by transitions between metastable states separated by sizable free energy barriers. Their practical application, in particular to ever more complex molecular systems, is, however, not entirely trivial. Focusing on replica exchange transition interface sampling (RETIS) and forward flux sampling (FFS), we discuss a range of analysis tools that can be used to assess the quality and convergence of such simulations, which is crucial to obtain reliable results. The basic ideas of a step-wise evaluation are exemplified for the study of nucleation in several systems with different complexities, providing a general guide for the critical assessment of RETIS and FFS simulations.
Collapse
Affiliation(s)
- Steven W. Hall
- Department of Chemical Engineering and Materials Science, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Grisell Díaz Leines
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| | - Sapna Sarupria
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
- Department of Chemical and Biomolecular Engineering, Clemson University, Clemson, South Carolina 29634, USA
| | - Jutta Rogal
- Department of Chemistry, New York University, New York, New York 10003, USA
- Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
7
|
Díaz Leines G, Rogal J. Template-Induced Precursor Formation in Heterogeneous Nucleation: Controlling Polymorph Selection and Nucleation Efficiency. PHYSICAL REVIEW LETTERS 2022; 128:166001. [PMID: 35522521 DOI: 10.1103/physrevlett.128.166001] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/26/2021] [Accepted: 03/29/2022] [Indexed: 06/14/2023]
Abstract
We present an atomistic study of heterogeneous nucleation in Ni employing transition path sampling, which reveals a template precursor-mediated mechanism of crystallization. Most notably, we find that the ability of tiny templates to modify the structural features of the liquid and promote the formation of precursor regions with enhanced bond-orientational order is key to determining their nucleation efficiency and the polymorphs that crystallize. Our results reveal an intrinsic link between structural liquid heterogeneity and the nucleating ability of templates, which significantly advances our understanding toward the control of nucleation efficiency and polymorph selection.
Collapse
Affiliation(s)
- Grisell Díaz Leines
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridgeshire CB2 1EW, United Kingdom
| | - Jutta Rogal
- Department of Chemistry, New York University, New York, New York 10003, USA and Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
8
|
Kikutsuji T, Mori Y, Okazaki KI, Mori T, Kim K, Matubayasi N. Explaining reaction coordinates of alanine dipeptide isomerization obtained from deep neural networks using Explainable Artificial Intelligence (XAI). J Chem Phys 2022; 156:154108. [DOI: 10.1063/5.0087310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A method for obtaining appropriate reaction coordinates is required to identify transition states distinguishing product and reactant in complex molecular systems. Recently, abundant research has been devoted to obtaining reaction coordinates using artificial neural networks from deep learning literature, where many collective variables are typically utilized in the input layer. However, it is difficult to explain the details of which collective variables contribute to the predicted reaction coordinates owing to the complexity of the nonlinear functions in deep neural networks. To overcome this limitation, we used Explainable Artificial Intelligence (XAI) methods of the Local Interpretable Model-agnostic Explanation (LIME) and the game theory-based framework known as Shapley Additive exPlanations (SHAP). We demonstrated that XAI enables us to obtain the degree of contribution of each collective variable to reaction coordinates that is determined by nonlinear regressions with deep learning for the committor of the alanine dipeptide isomerization in vacuum. In particular, both LIME and SHAP provide important features to the predicted reaction coordinates, which are characterized by appropriate dihedral angles consistent with those previously reported from the committor test analysis. The present study offers an AI-aided framework to explain the appropriate reaction coordinates, which acquires considerable significance when the number of degrees of freedom increases.
Collapse
Affiliation(s)
| | | | - Kei-ichi Okazaki
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Japan
| | - Toshifumi Mori
- Kyushu University Institute for Materials Chemistry and Engineering, Japan
| | - Kang Kim
- Graduate School of Engineering Science, Osaka University - Toyonaka Campus, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Japan
| |
Collapse
|
9
|
Bolhuis PG, Swenson DWH. Transition Path Sampling as Markov Chain Monte Carlo of Trajectories: Recent Algorithms, Software, Applications, and Future Outlook. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202000237] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Peter G. Bolhuis
- Amsterdam Center for Multiscale Modeling van 't Hoff Institute for Molecular Sciences University of Amsterdam PO Box 94157 1090 GD Amsterdam The Netherlands
| | - David W. H. Swenson
- Centre Blaise Pascal Ecole Normale Superieure 46, allée d'Italie 69364 Lyon Cedex 07 France
| |
Collapse
|
10
|
Díaz Leines G, Michaelides A, Rogal J. Interplay of structural and dynamical heterogeneity in the nucleation mechanism in Ni. Faraday Discuss 2021; 235:406-415. [DOI: 10.1039/d1fd00099c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Gaining fundamental understanding of crystal nucleation processes in metal alloys is crucial for the development and design of high-performance materials with targeted properties. Yet, crystallizationis a complex non-equilibrium process and,...
Collapse
|
11
|
Menon S, Díaz Leines G, Drautz R, Rogal J. Role of pre-ordered liquid in the selection mechanism of crystal polymorphs during nucleation. J Chem Phys 2020; 153:104508. [DOI: 10.1063/5.0017575] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sarath Menon
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Grisell Díaz Leines
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Ralf Drautz
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44801 Bochum, Germany
| | - Jutta Rogal
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44801 Bochum, Germany
| |
Collapse
|
12
|
Mori Y, Okazaki KI, Mori T, Kim K, Matubayasi N. Learning reaction coordinates via cross-entropy minimization: Application to alanine dipeptide. J Chem Phys 2020; 153:054115. [DOI: 10.1063/5.0009066] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Affiliation(s)
- Yusuke Mori
- Division of Chemical Engineering, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| | | | - Toshifumi Mori
- Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan
- The Graduate University for Advanced Studies, Okazaki, Aichi 444-8585, Japan
| | - Kang Kim
- Division of Chemical Engineering, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
- Institute for Molecular Science, Okazaki, Aichi 444-8585, Japan
| | - Nobuyuki Matubayasi
- Division of Chemical Engineering, Department of Materials Engineering Science, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
13
|
Liang Y, Díaz Leines G, Drautz R, Rogal J. Identification of a multi-dimensional reaction coordinate for crystal nucleation in Ni3Al. J Chem Phys 2020; 152:224504. [DOI: 10.1063/5.0010074] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
- Yanyan Liang
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Grisell Díaz Leines
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Ralf Drautz
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Jutta Rogal
- Interdisciplinary Centre for Advanced Materials Simulation, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|