1
|
Endo K, Yuhara D, Yasuoka K. Efficient Monte Carlo Sampling for Molecular Systems Using Continuous Normalizing Flow. J Chem Theory Comput 2022; 18:1395-1405. [PMID: 35175774 DOI: 10.1021/acs.jctc.1c01047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Monte Carlo molecular simulation is a powerful computational method for simulating molecular behavior. It generates samples of the possible states of molecular systems. To generate a sample efficiently, it is advantageous to avoid suggesting extremely high-energy states that would never become possible states. In this study, we propose a new sampling method for Monte Carlo molecular simulation, that is, a continuous normalizing molecular flow (CNMF) method, which can create various probabilistic distributions of molecular states from some initial distribution. The CNMF method generates samples by solving a first-order differential equation with two-body intermolecular interaction terms. We also develop specific probabilistic distributions using CNMF called inverse square flow, which yields distributions with zero probability density when molecule pairs are in close proximity, whereas probability densities are compressed uniformly from the initial distribution in all other cases. Using inverse square flow, we demonstrate that Monte Carlo molecular simulation is more efficient than the standard simulation. Although the increased computational costs of the CNMF method are non-negligible, this method is feasible for parallel computation and has the potential for expansion.
Collapse
Affiliation(s)
- Katsuhiro Endo
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama, 223-8522, Japan
| | - Daisuke Yuhara
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama, 223-8522, Japan.,Science & Innovation Center, Mitsubishi Chemical Corporation, Yokohama, 227-8502, Japan
| | - Kenji Yasuoka
- Department of Mechanical Engineering, Keio University, 3-14-1 Hiyoshi, Kohoku-ku, Yokohama, 223-8522, Japan
| |
Collapse
|
2
|
The role of structural symmetry on proton tautomerization: A DFTB/Meta-Dynamics computational study. Chem Phys 2021. [DOI: 10.1016/j.chemphys.2021.111222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
3
|
Jadrich RB, Ticknor C, Leiding JA. First principles reactive simulation for equation of state prediction. J Chem Phys 2021; 154:244307. [PMID: 34241343 DOI: 10.1063/5.0050676] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The high cost of density functional theory (DFT) has hitherto limited the ab initio prediction of the equation of state (EOS). In this article, we employ a combination of large scale computing, advanced simulation techniques, and smart data science strategies to provide an unprecedented ab initio performance analysis of the high explosive pentaerythritol tetranitrate (PETN). Comparison to both experiment and thermochemical predictions reveals important quantitative limitations of DFT for EOS prediction and thus the assessment of high explosives. In particular, we find that DFT predicts the energy of PETN detonation products to be systematically too high relative to the unreacted neat crystalline material, resulting in an underprediction of the detonation velocity, pressure, and temperature at the Chapman-Jouguet state. The energetic bias can be partially accounted for by high-level electronic structure calculations of the product molecules. We also demonstrate a modeling strategy for mapping chemical composition across a wide parameter space with limited numerical data, the results of which suggest additional molecular species to consider in thermochemical modeling.
Collapse
Affiliation(s)
- Ryan B Jadrich
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Christopher Ticknor
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - Jeffery A Leiding
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| |
Collapse
|
4
|
Maerzke KA, Yoon TJ, Jadrich RB, Leiding JA, Currier RP. First-Principles Simulations of CuCl in High-Temperature Water Vapor. J Phys Chem B 2021; 125:4794-4807. [PMID: 33938730 DOI: 10.1021/acs.jpcb.1c00083] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Experimental data suggest that the solubility of copper in high-temperature water vapor is controlled by the formation of hydrated clusters of the form CuCl(H2O)n, where the average number of water molecules in the cluster generally increases with increasing density [Migdisov, A. A.; et al. Geochim. Cosmochim. Acta 2014, 129, 33-53]. However, the precise nature of these clusters is difficult to probe experimentally. Moreover, there are some discrepancies between experimental estimates of average cluster size and prior simulation work [Mei, Y. Geofluids 2018, 2018, 4279124]. We have performed first-principles Monte Carlo (MC) and molecular dynamics (MD) simulations to explore these clusters in finer detail. We find that molecular dynamics is not the most appropriate technique for studying aggregation in vapor phases, even at relatively high temperatures. Specifically, our MD simulations exhibit substantial problems in adequately sampling the equilibrium cluster size distribution. In contrast, MC simulations with specialized cluster moves are able to accurately sample the phase space of hydrogen-bonding vapors. At all densities, we find a stable, slightly distorted linear H2O-Cu-Cl structure, which is in agreement with the earlier simulations, surrounded by a variable number of water molecules. The surrounding water molecules do not form a well-defined second solvation shell but rather a loose network of hydrogen-bonded water with molecular CuCl on the outside edge of the water cluster. We also find a broad distribution of hydration numbers, especially at higher densities. In contrast to previous simulation work but in agreement with experimental data, we find that the average hydration number substantially increases with increasing density. Moreover, the value of the hydration number depends on the choice of cluster definition.
Collapse
Affiliation(s)
- Katie A Maerzke
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Tae Jun Yoon
- Materials Physics and Applications Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Ryan B Jadrich
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States.,Center for Nonlinear Studies, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Jeffery A Leiding
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Robert P Currier
- Chemistry Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| |
Collapse
|
5
|
Gilabert JF, Lecina D, Estrada J, Guallar V. Monte Carlo Techniques for Drug Design: The Success Case of PELE. BIOMOLECULAR SIMULATIONS IN STRUCTURE-BASED DRUG DISCOVERY 2018. [DOI: 10.1002/9783527806836.ch5] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Joan F. Gilabert
- Barcelona Supercomputing Center (BSC); Life Science Department; Jordi Girona 29 08034 Barcelona Spain
| | - Daniel Lecina
- Barcelona Supercomputing Center (BSC); Life Science Department; Jordi Girona 29 08034 Barcelona Spain
| | - Jorge Estrada
- Barcelona Supercomputing Center (BSC); Life Science Department; Jordi Girona 29 08034 Barcelona Spain
| | - Victor Guallar
- Barcelona Supercomputing Center (BSC); Life Science Department; Jordi Girona 29 08034 Barcelona Spain
- ICREA; Passeig Lluís Companys 23 08010 Barcelona Spain
| |
Collapse
|
6
|
Mitchell I, Aradi B, Page AJ. Density functional tight binding-based free energy simulations in the DFTB+ program. J Comput Chem 2018; 39:2452-2458. [DOI: 10.1002/jcc.25583] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2018] [Revised: 08/14/2018] [Accepted: 08/15/2018] [Indexed: 01/11/2023]
Affiliation(s)
- Izaac Mitchell
- School of Environmental and Life Sciences, University of Newcastle
| | - Bálint Aradi
- Bremen Center for Computational Materials Science, University of Bremen
| | - Alister J. Page
- School of Environmental and Life Sciences, University of Newcastle
| |
Collapse
|
7
|
Aussems DUB, Bal KM, Morgan TW, van de Sanden MCM, Neyts EC. Atomistic simulations of graphite etching at realistic time scales. Chem Sci 2017; 8:7160-7168. [PMID: 29081947 PMCID: PMC5635421 DOI: 10.1039/c7sc02763j] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Accepted: 08/23/2017] [Indexed: 11/21/2022] Open
Abstract
Hydrogen-graphite interactions are relevant to a wide variety of applications, ranging from astrophysics to fusion devices and nano-electronics. In order to shed light on these interactions, atomistic simulation using Molecular Dynamics (MD) has been shown to be an invaluable tool. It suffers, however, from severe time-scale limitations. In this work we apply the recently developed Collective Variable-Driven Hyperdynamics (CVHD) method to hydrogen etching of graphite for varying inter-impact times up to a realistic value of 1 ms, which corresponds to a flux of ∼1020 m-2 s-1. The results show that the erosion yield, hydrogen surface coverage and species distribution are significantly affected by the time between impacts. This can be explained by the higher probability of C-C bond breaking due to the prolonged exposure to thermal stress and the subsequent transition from ion- to thermal-induced etching. This latter regime of thermal-induced etching - chemical erosion - is here accessed for the first time using atomistic simulations. In conclusion, this study demonstrates that accounting for long time-scales significantly affects ion bombardment simulations and should not be neglected in a wide range of conditions, in contrast to what is typically assumed.
Collapse
Affiliation(s)
- D U B Aussems
- DIFFER - Dutch Institute for Fundamental Energy Research , De Zaale 20 , 5612 AJ Eindhoven , The Netherlands .
| | - K M Bal
- University of Antwerp , Department of Chemistry , PLASMANT Research Group , Universiteitsplein 1 , 2610 Antwerp , Belgium
| | - T W Morgan
- DIFFER - Dutch Institute for Fundamental Energy Research , De Zaale 20 , 5612 AJ Eindhoven , The Netherlands .
| | - M C M van de Sanden
- DIFFER - Dutch Institute for Fundamental Energy Research , De Zaale 20 , 5612 AJ Eindhoven , The Netherlands .
- Eindhoven University of Technology , PO Box 513 , 5600 MB Eindhoven , The Netherlands
| | - E C Neyts
- University of Antwerp , Department of Chemistry , PLASMANT Research Group , Universiteitsplein 1 , 2610 Antwerp , Belgium
| |
Collapse
|
8
|
|
9
|
Mitchell I, Irle S, Page AJ. A global reaction route mapping-based kinetic Monte Carlo algorithm. J Chem Phys 2016; 145:024105. [DOI: 10.1063/1.4954660] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Affiliation(s)
- Izaac Mitchell
- Newcastle Institute for Energy and Resources, The University of Newcastle, Callaghan 2308, Australia
| | - Stephan Irle
- Institute of Transformative Bio-Molecules (WPI-ITbM) and Department of Chemistry, Graduate School of Science, Nagoya University, Nagoya 464-8602, Japan
| | - Alister J. Page
- Newcastle Institute for Energy and Resources, The University of Newcastle, Callaghan 2308, Australia
| |
Collapse
|
10
|
Atomic scale simulation of carbon nanotube nucleation from hydrocarbon precursors. Nat Commun 2015; 6:10306. [PMID: 26691537 PMCID: PMC4703880 DOI: 10.1038/ncomms10306] [Citation(s) in RCA: 54] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2015] [Accepted: 11/27/2015] [Indexed: 11/08/2022] Open
Abstract
Atomic scale simulations of the nucleation and growth of carbon nanotubes is essential for understanding their growth mechanism. In spite of over twenty years of simulation efforts in this area, limited progress has so far been made on addressing the role of the hydrocarbon growth precursor. Here we report on atomic scale simulations of cap nucleation of single-walled carbon nanotubes from hydrocarbon precursors. The presented mechanism emphasizes the important role of hydrogen in the nucleation process, and is discussed in relation to previously presented mechanisms. In particular, the role of hydrogen in the appearance of unstable carbon structures during in situ experimental observations as well as the initial stage of multi-walled carbon nanotube growth is discussed. The results are in good agreement with available experimental and quantum-mechanical results, and provide a basic understanding of the incubation and nucleation stages of hydrocarbon-based CNT growth at the atomic level. Atomic scale simulation of the nucleation and growth of carbon nanotubes is essential for understanding their growth mechanism. Here, the authors look at cap nucleation of nanotubes from hydrocarbon precursors, specifically probing the role of hydrogen in the early stages of growth.
Collapse
|
11
|
Neumann T, Danilov D, Wenzel W. Multiparticle moves in acceptance rate optimized monte carlo. J Comput Chem 2015; 36:2236-45. [DOI: 10.1002/jcc.24205] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2015] [Revised: 08/31/2015] [Accepted: 09/02/2015] [Indexed: 11/11/2022]
Affiliation(s)
- Tobias Neumann
- Institute of Nanotechnology, Karlsruhe Institute of Technology; PO Box 3640, D-76021 Karlsruhe Germany
| | - Denis Danilov
- Institute of Nanotechnology, Karlsruhe Institute of Technology; PO Box 3640, D-76021 Karlsruhe Germany
| | | |
Collapse
|
12
|
Bal KM, Neyts EC. Merging Metadynamics into Hyperdynamics: Accelerated Molecular Simulations Reaching Time Scales from Microseconds to Seconds. J Chem Theory Comput 2015; 11:4545-54. [DOI: 10.1021/acs.jctc.5b00597] [Citation(s) in RCA: 66] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Kristof M. Bal
- Department
of Chemistry, University of Antwerp, Universiteitsplein 1, 2610 Wilrijk, Antwerp, Belgium
| | - Erik C. Neyts
- Department
of Chemistry, University of Antwerp, Universiteitsplein 1, 2610 Wilrijk, Antwerp, Belgium
| |
Collapse
|
13
|
Bal KM, Neyts EC. On the time scale associated with Monte Carlo simulations. J Chem Phys 2014; 141:204104. [DOI: 10.1063/1.4902136] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Affiliation(s)
- Kristof M. Bal
- Department of Chemistry, University of Antwerp, Research Group PLASMANT, Universiteitsplein 1, 2610 Wilrijk, Antwerp, Belgium
| | - Erik C. Neyts
- Department of Chemistry, University of Antwerp, Research Group PLASMANT, Universiteitsplein 1, 2610 Wilrijk, Antwerp, Belgium
| |
Collapse
|
14
|
Engelmann Y, Bogaerts A, Neyts EC. Thermodynamics at the nanoscale: phase diagrams of nickel-carbon nanoclusters and equilibrium constants for phase transitions. NANOSCALE 2014; 6:11981-11987. [PMID: 25177915 DOI: 10.1039/c4nr02354d] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
Using reactive molecular dynamics simulations, the melting behavior of nickel-carbon nanoclusters is examined. The phase diagrams of icosahedral and Wulff polyhedron clusters are determined using both the Lindemann index and the potential energy. Formulae are derived for calculating the equilibrium constants and the solid and liquid fractions during a phase transition, allowing more rational determination of the melting temperature with respect to the arbitrary Lindemann value. These results give more insight into the properties of nickel-carbon nanoclusters in general and can specifically be very useful for a better understanding of the synthesis of carbon nanotubes using the catalytic chemical vapor deposition method.
Collapse
Affiliation(s)
- Yannick Engelmann
- Research Group PLASMANT, Department of Chemistry, University of Antwerp, Universiteitsplein 1, 2610, Antwerp, Belgium.
| | | | | |
Collapse
|
15
|
Combining molecular dynamics with Monte Carlo simulations: implementations and applications. HIGHLIGHTS IN THEORETICAL CHEMISTRY 2014. [DOI: 10.1007/978-3-642-41315-5_23] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
|
16
|
Neyts EC, van Duin ACT, Bogaerts A. Formation of single layer graphene on nickel under far-from-equilibrium high flux conditions. NANOSCALE 2013; 5:7250-7255. [PMID: 23695014 DOI: 10.1039/c3nr00153a] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
Abstract
We investigate the theoretical possibility of single layer graphene formation on a nickel surface at different substrate temperatures under far-from-equilibrium high precursor flux conditions, employing state-of-the-art hybrid reactive molecular dynamics/uniform acceptance force bias Monte Carlo simulations. It is predicted that under these conditions, the formation of a single layer graphene-like film may proceed through a combined deposition-segregation mechanism on a nickel substrate, rather than by pure surface segregation as is typically observed for metals with high carbon solubility. At 900 K and above, nearly continuous graphene layers are obtained. These simulations suggest that single layer graphene deposition is theoretically possible on Ni under high flux conditions.
Collapse
Affiliation(s)
- Erik C Neyts
- Department of Chemistry, University of Antwerp, Research group PLASMANT, Universiteitsplein 1, 2610 Antwerp, Belgium.
| | | | | |
Collapse
|
17
|
Neyts EC, Bogaerts A. Combining molecular dynamics with Monte Carlo simulations: implementations and applications. Theor Chem Acc 2012. [DOI: 10.1007/s00214-012-1320-x] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|