1
|
Zheng J, Frisch MJ. Multiple-time scale integration method based on an interpolated potential energy surface for ab initio path integral molecular dynamics. J Chem Phys 2024; 160:144111. [PMID: 38597307 DOI: 10.1063/5.0196634] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2024] [Accepted: 03/25/2024] [Indexed: 04/11/2024] Open
Abstract
A new multiple-time scale integration method is presented that propagates ab initio path integral molecular dynamics (PIMD). This method uses a large time step to generate an approximate geometrical configuration whose energy and gradient are evaluated at the level of an ab initio method, and then, a more precise integration scheme, e.g., the Bulirsch-Stoer method or velocity Verlet integration with a smaller time step, is used to integrate from the previous step using the computationally efficient interpolated potential energy surface constructed from two consecutive points. This method makes the integration of PIMD more efficient and accurate compared with the velocity Verlet integration. A Nosé-Hoover chain thermostat combined with this new multiple-time scale method has good energy conservation even with a large time step, which is usually challenging in velocity Verlet integration for PIMD due to the very small chain mass when a large number of beads are used. The new method is used to calculate infrared spectra and free energy profiles to demonstrate its accuracy and capabilities.
Collapse
Affiliation(s)
- Jingjing Zheng
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| | - Michael J Frisch
- Gaussian, Inc., 340 Quinnipiac St. Bldg. 40, Wallingford, Connecticut 06492, USA
| |
Collapse
|
2
|
London N, Limbu DK, Momeni MR, Shakib FA. DL_POLY Quantum 2.0: A modular general-purpose software for advanced path integral simulations. J Chem Phys 2024; 160:132501. [PMID: 38557837 DOI: 10.1063/5.0197822] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2024] [Accepted: 03/12/2024] [Indexed: 04/04/2024] Open
Abstract
DL_POLY Quantum 2.0, a vastly expanded software based on DL_POLY Classic 1.10, is a highly parallelized computational suite written in FORTRAN77 with a modular structure for incorporating nuclear quantum effects into large-scale/long-time molecular dynamics simulations. This is achieved by presenting users with a wide selection of state-of-the-art dynamics methods that utilize the isomorphism between a classical ring polymer and Feynman's path integral formalism of quantum mechanics. The flexible and user-friendly input/output handling system allows the control of methodology, integration schemes, and thermostatting. DL_POLY Quantum is equipped with a module specifically assigned for calculating correlation functions and printing out the values for sought-after quantities, such as dipole moments and center-of-mass velocities, with packaged tools for calculating infrared absorption spectra and diffusion coefficients.
Collapse
Affiliation(s)
- Nathan London
- Division of Energy, Matter and Systems, School of Science and Engineering, University of Missouri-Kansas City, Kansas City, Missouri 64110, USA
| | - Dil K Limbu
- Division of Energy, Matter and Systems, School of Science and Engineering, University of Missouri-Kansas City, Kansas City, Missouri 64110, USA
- Department of Chemistry and Environmental Science, New Jersey Institute of Technology, Newark, New Jersey 07102, USA
| | - Mohammad R Momeni
- Division of Energy, Matter and Systems, School of Science and Engineering, University of Missouri-Kansas City, Kansas City, Missouri 64110, USA
| | - Farnaz A Shakib
- Department of Chemistry and Environmental Science, New Jersey Institute of Technology, Newark, New Jersey 07102, USA
| |
Collapse
|
3
|
Lan J, Chergui M, Pasquarello A. Dynamics of the charge transfer to solvent process in aqueous iodide. Nat Commun 2024; 15:2544. [PMID: 38514610 DOI: 10.1038/s41467-024-46772-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/21/2023] [Accepted: 03/05/2024] [Indexed: 03/23/2024] Open
Abstract
Charge-transfer-to-solvent states in aqueous halides are ideal systems for studying the electron-transfer dynamics to the solvent involving a complex interplay between electronic excitation and solvent polarization. Despite extensive experimental investigations, a full picture of the charge-transfer-to-solvent dynamics has remained elusive. Here, we visualise the intricate interplay between the dynamics of the electron and the solvent polarization occurring in this process. Through the combined use of ab initio molecular dynamics and machine learning methods, we investigate the structure, dynamics and free energy as the excited electron evolves through the charge-transfer-to-solvent process, which we characterize as a sequence of states denoted charge-transfer-to-solvent, contact-pair, solvent-separated, and hydrated electron states, depending on the distance between the iodine and the excited electron. Our assignment of the charge-transfer-to-solvent states is supported by the good agreement between calculated and measured vertical binding energies. Our results reveal the charge transfer process in terms of the underlying atomic processes and mechanisms.
Collapse
Affiliation(s)
- Jinggang Lan
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland.
- Department of Chemistry, New York University, New York, NY, 10003, USA.
- Simons Center for Computational Physical Chemistry at New York University, New York, NY, 10003, USA.
| | - Majed Chergui
- Lausanne Centre for Ultrafast Science (LACUS), ISIC, Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
- Elettra - Sincrotrone Trieste, Area Science Park I - 34149, Trieste, Italy
| | - Alfredo Pasquarello
- Chaire de Simulation à l'Echelle Atomique (CSEA), Ecole Polytechnique Fédérale de Lausanne (EPFL), CH-1015, Lausanne, Switzerland
| |
Collapse
|
4
|
Kurnikov IV, Pereyaslavets L, Kamath G, Sakipov SN, Voronina E, Butin O, Illarionov A, Leontyev I, Nawrocki G, Darkhovskiy M, Olevanov M, Ivahnenko I, Chen Y, Lock CB, Levitt M, Kornberg RD, Fain B. Neural Network Corrections to Intermolecular Interaction Terms of a Molecular Force Field Capture Nuclear Quantum Effects in Calculations of Liquid Thermodynamic Properties. J Chem Theory Comput 2024; 20:1347-1357. [PMID: 38240485 PMCID: PMC11042917 DOI: 10.1021/acs.jctc.3c00921] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/19/2024]
Abstract
We incorporate nuclear quantum effects (NQE) in condensed matter simulations by introducing short-range neural network (NN) corrections to the ab initio fitted molecular force field ARROW. Force field NN corrections are fitted to average interaction energies and forces of molecular dimers, which are simulated using the Path Integral Molecular Dynamics (PIMD) technique with restrained centroid positions. The NN-corrected force field allows reproduction of the NQE for computed liquid water and methane properties such as density, radial distribution function (RDF), heat of evaporation (HVAP), and solvation free energy. Accounting for NQE through molecular force field corrections circumvents the need for explicit computationally expensive PIMD simulations in accurate calculations of the properties of chemical and biological systems. The accuracy and locality of pairwise NN NQE corrections indicate that this approach could be applicable to complex heterogeneous systems, such as proteins.
Collapse
Affiliation(s)
- Igor V Kurnikov
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - Leonid Pereyaslavets
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - Ganesh Kamath
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - Serzhan N Sakipov
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - Ekaterina Voronina
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - Oleg Butin
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - Alexey Illarionov
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - Igor Leontyev
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - Grzegorz Nawrocki
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - Mikhail Darkhovskiy
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - Michael Olevanov
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - Ilya Ivahnenko
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - YuChun Chen
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| | - Christopher B Lock
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
- Department of Neurology and Neurological Sciences, Stanford University School of Medicine, Palo Alto, California 94304, United States
| | - Michael Levitt
- Department of Structural Biology, Stanford University School of Medicine, Stanford, California 94305, United States
| | - Roger D Kornberg
- Department of Structural Biology, Stanford University School of Medicine, Stanford, California 94305, United States
| | - Boris Fain
- InterX Inc., (a Subsidiary of NeoTX Therapeutics Ltd.), 805 Allston Way, Berkeley, California 94710, United States
| |
Collapse
|
5
|
Kapil V, Kovács DP, Csányi G, Michaelides A. First-principles spectroscopy of aqueous interfaces using machine-learned electronic and quantum nuclear effects. Faraday Discuss 2024; 249:50-68. [PMID: 37799072 PMCID: PMC10845015 DOI: 10.1039/d3fd00113j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Accepted: 07/18/2023] [Indexed: 10/07/2023]
Abstract
Vibrational spectroscopy is a powerful approach to visualising interfacial phenomena. However, extracting structural and dynamical information from vibrational spectra is a challenge that requires first-principles simulations, including non-Condon and quantum nuclear effects. We address this challenge by developing a machine-learning enhanced first-principles framework to speed up predictive modelling of infrared, Raman, and sum-frequency generation spectra. Our approach uses machine learning potentials that encode quantum nuclear effects to generate quantum trajectories using simple molecular dynamics efficiently. In addition, we reformulate bulk and interfacial selection rules to express them unambiguously in terms of the derivatives of polarisation and polarisabilities of the whole system and predict these derivatives efficiently using fully-differentiable machine learning models of dielectric response tensors. We demonstrate our framework's performance by predicting the IR, Raman, and sum-frequency generation spectra of liquid water, ice and the water-air interface by achieving near quantitative agreement with experiments at nearly the same computational efficiency as pure classical methods. Finally, to aid the experimental discovery of new phases of nanoconfined water, we predict the temperature-dependent vibrational spectra of monolayer water across the solid-hexatic-liquid phases transition.
Collapse
Affiliation(s)
- Venkat Kapil
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK.
| | | | - Gábor Csányi
- Engineering Laboratory, University of Cambridge, Cambridge, CB2 1PZ, UK
| | - Angelos Michaelides
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, UK.
| |
Collapse
|
6
|
Brezina K, Beck H, Marsalek O. Reducing the Cost of Neural Network Potential Generation for Reactive Molecular Systems. J Chem Theory Comput 2023; 19:6589-6604. [PMID: 37747971 PMCID: PMC10569056 DOI: 10.1021/acs.jctc.3c00391] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/06/2023] [Indexed: 09/27/2023]
Abstract
Although machine learning potentials have recently had a substantial impact on molecular simulations, the construction of a robust training set can still become a limiting factor, especially due to the requirement of a reference ab initio simulation that covers all the relevant geometries of the system. Recognizing that this can be prohibitive for certain systems, we develop the method of transition tube sampling that mitigates the computational cost of training set and model generation. In this approach, we generate classical or quantum thermal geometries around a transition path describing a conformational change or a chemical reaction using only a sparse set of local normal mode expansions along this path and select from these geometries by an active learning protocol. This yields a training set with geometries that characterize the whole transition without the need for a costly reference trajectory. The performance of the method is evaluated on different molecular systems with the complexity of the potential energy landscape increasing from a single minimum to a double proton-transfer reaction with high barriers. Our results show that the method leads to training sets that give rise to models applicable in classical and path integral simulations alike that are on par with those based directly on ab initio calculations while providing the computational speedup we have come to expect from machine learning potentials.
Collapse
Affiliation(s)
- Krystof Brezina
- Charles University, Faculty of Mathematics
and Physics, Ke Karlovu
3, 121 16, Prague
2, Czech Republic
| | - Hubert Beck
- Charles University, Faculty of Mathematics
and Physics, Ke Karlovu
3, 121 16, Prague
2, Czech Republic
| | - Ondrej Marsalek
- Charles University, Faculty of Mathematics
and Physics, Ke Karlovu
3, 121 16, Prague
2, Czech Republic
| |
Collapse
|
7
|
Khuu T, Schleif T, Mohamed A, Mitra S, Johnson MA, Valdiviezo J, Heindel JP, Head-Gordon T. Intra-cluster Charge Migration upon Hydration of Protonated Formic Acid Revealed by Anharmonic Analysis of Cold Ion Vibrational Spectra. J Phys Chem A 2023; 127:7501-7509. [PMID: 37669457 DOI: 10.1021/acs.jpca.3c03971] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/07/2023]
Abstract
The rates of many chemical reactions are accelerated when carried out in micron-sized droplets, but the molecular origin of the rate acceleration remains unclear. One example is the condensation reaction of 1,2-diaminobenzene with formic acid to yield benzimidazole. The observed rate enhancements have been rationalized by invoking enhanced acidity at the surface of methanol solvent droplets with low water content to enable protonation of formic acid to generate a cationic species (protonated formic acid or PFA) formed by attachment of a proton to the neutral acid. Because PFA is the key feature in this reaction mechanism, vibrational spectra of cryogenically cooled, microhydrated PFA·(H2O)n=1-6 were acquired to determine how the extent of charge localization depends on the degree of hydration. Analysis of these highly anharmonic spectra with path integral ab initio molecular dynamics simulations reveals the gradual displacement of the excess proton onto the water network in the microhydration regime at low temperatures with n = 3 as the tipping point for intra-cluster proton transfer.
Collapse
Affiliation(s)
- Thien Khuu
- Sterling Chemistry Laboratory, Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Tim Schleif
- Sterling Chemistry Laboratory, Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Ahmed Mohamed
- Sterling Chemistry Laboratory, Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Sayoni Mitra
- Sterling Chemistry Laboratory, Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Mark A Johnson
- Sterling Chemistry Laboratory, Department of Chemistry, Yale University, New Haven, Connecticut 06511, United States
| | - Jesús Valdiviezo
- Pitzer Theory Center, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Joseph P Heindel
- Pitzer Theory Center, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Teresa Head-Gordon
- Pitzer Theory Center, Department of Chemistry, University of California, Berkeley, California 94720, United States
| |
Collapse
|
8
|
Atsango AO, Morawietz T, Marsalek O, Markland TE. Developing machine-learned potentials to simultaneously capture the dynamics of excess protons and hydroxide ions in classical and path integral simulations. J Chem Phys 2023; 159:074101. [PMID: 37581418 DOI: 10.1063/5.0162066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2023] [Accepted: 07/31/2023] [Indexed: 08/16/2023] Open
Abstract
The transport of excess protons and hydroxide ions in water underlies numerous important chemical and biological processes. Accurately simulating the associated transport mechanisms ideally requires utilizing ab initio molecular dynamics simulations to model the bond breaking and formation involved in proton transfer and path-integral simulations to model the nuclear quantum effects relevant to light hydrogen atoms. These requirements result in a prohibitive computational cost, especially at the time and length scales needed to converge proton transport properties. Here, we present machine-learned potentials (MLPs) that can model both excess protons and hydroxide ions at the generalized gradient approximation and hybrid density functional theory levels of accuracy and use them to perform multiple nanoseconds of both classical and path-integral proton defect simulations at a fraction of the cost of the corresponding ab initio simulations. We show that the MLPs are able to reproduce ab initio trends and converge properties such as the diffusion coefficients of both excess protons and hydroxide ions. We use our multi-nanosecond simulations, which allow us to monitor large numbers of proton transfer events, to analyze the role of hypercoordination in the transport mechanism of the hydroxide ion and provide further evidence for the asymmetry in diffusion between excess protons and hydroxide ions.
Collapse
Affiliation(s)
- Austin O Atsango
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Tobias Morawietz
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Ondrej Marsalek
- Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic
| | - Thomas E Markland
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
9
|
Schiltz C, Rappoport D, Mandelshtam VA. Implementation of the self-consistent phonons method with ab initio potentials (AI-SCP). J Chem Phys 2023; 158:2890485. [PMID: 37184023 DOI: 10.1063/5.0146682] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Accepted: 04/24/2023] [Indexed: 05/16/2023] Open
Abstract
The self-consistent phonon (SCP) method allows one to include anharmonic effects when treating a many-body quantum system at thermal equilibrium. The system is then described by an effective temperature-dependent harmonic Hamiltonian, which can be used to estimate its various dynamic and static properties. In this paper, we combine SCP with ab initio (AI) potential energy evaluation in which case the numerical bottleneck of AI-SCP is the evaluation of Gaussian averages of the AI potential energy and its derivatives. These averages are computed efficiently by the quasi-Monte Carlo method utilizing low-discrepancy sequences leading to a fast convergence with respect to the number, S, of the AI energy evaluations. Moreover, a further substantial (an-order-of-magnitude) improvement in efficiency is achieved once a numerically cheap approximation of the AI potential is available. This is based on using a perturbation theory-like (the two-grid) approach in which it is the average of the difference between the AI and the approximate potential that is computed. The corresponding codes and scripts are provided.
Collapse
Affiliation(s)
- Colin Schiltz
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | - Dmitrij Rappoport
- Department of Chemistry, University of California, Irvine, California 92697, USA
| | | |
Collapse
|
10
|
Chen Z, Yang Y. Incorporating Nuclear Quantum Effects in Molecular Dynamics with a Constrained Minimized Energy Surface. J Phys Chem Lett 2023; 14:279-286. [PMID: 36595586 DOI: 10.1021/acs.jpclett.2c02905] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The accurate incorporation of nuclear quantum effects in large-scale molecular dynamics (MD) simulations remains a significant challenge. Recently, we combined constrained nuclear-electronic orbital (CNEO) theory with classical MD and obtained a new approach (CNEO-MD) that can accurately and efficiently incorporate nuclear quantum effects into classical simulations. In this Letter, we provide the theoretical foundation for CNEO-MD by developing an alternative formulation of the equations of motion for MD. In this new formulation, the expectation values of quantum nuclear positions evolve classically on an effective energy surface that is obtained from a constrained energy minimization procedure when solving for the quantum nuclear wave function, thus enabling the incorporation of nuclear quantum effects in classical MD simulations. For comparison with other existing approaches, we examined a series of model systems and found that this new MD approach is significantly more accurate than the conventional way of performing classical MD and generally outperforms centroid MD and ring-polymer MD in describing vibrations in these model systems.
Collapse
Affiliation(s)
- Zehua Chen
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin53706, United States
| | - Yang Yang
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin53706, United States
| |
Collapse
|
11
|
Lan J, Rybkin VV, Pasquarello A. Temperature Dependent Properties of the Aqueous Electron. Angew Chem Int Ed Engl 2022; 61:e202209398. [PMID: 35849110 PMCID: PMC9541610 DOI: 10.1002/anie.202209398] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2022] [Indexed: 11/07/2022]
Abstract
The temperature‐dependent properties of the aqueous electron have been extensively studied using mixed quantum‐classical simulations in a wide range of thermodynamic conditions based on one‐electron pseudopotentials. While the cavity model appears to explain most of the physical properties of the aqueous electron, only a non‐cavity model has so far been successful in accounting for the temperature dependence of the absorption spectrum. Here, we present an accurate and efficient description of the aqueous electron under various thermodynamic conditions by combining hybrid functional‐based molecular dynamics, machine learning techniques, and multiple time‐step methods. Our advanced simulations accurately describe the temperature dependence of the absorption maximum in the presence of cavity formation. Specifically, our work reveals that the red shift of the absorption maximum results from an increasing gyration radius with temperature, rather than from global density variations as previously suggested.
Collapse
Affiliation(s)
- Jinggang Lan
- Chaire de Simulation àl'Echelle Atomique (CSEA)Ecole Polytechnique Fédérale de Lausanne (EPFL)CH-1015LausanneSwitzerland
| | | | - Alfredo Pasquarello
- Chaire de Simulation àl'Echelle Atomique (CSEA)Ecole Polytechnique Fédérale de Lausanne (EPFL)CH-1015LausanneSwitzerland
| |
Collapse
|
12
|
Lan J, Rybkin VV, Pasquarello A. Temperature Dependent Properties of the Aqueous Electron. Angew Chem Int Ed Engl 2022. [DOI: 10.1002/ange.202209398] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Jinggang Lan
- EPFL: Ecole Polytechnique Federale de Lausanne Chaire de Simulation à l’Echelle Atomique 1015 Lausanne SWITZERLAND
| | | | - Alfredo Pasquarello
- EPFL: Ecole Polytechnique Federale de Lausanne Chaire de Simulation à l’Echelle Atomique SWITZERLAND
| |
Collapse
|
13
|
Impact of nuclear quantum effects on the $$^{13}$$C nuclear shielding of linear carboxylates in water. Theor Chem Acc 2022. [DOI: 10.1007/s00214-022-02894-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
14
|
Juraskova V, Celerse F, Laplaza R, Corminboeuf C. Assessing the persistence of chalcogen bonds in solution with neural network potentials. J Chem Phys 2022; 156:154112. [DOI: 10.1063/5.0085153] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Non-covalent bonding patterns are commonly harvested as a design principle in the field of catalysis, supramolecular chemistry and functional materials to name a few. Yet, their computational description generally neglects finite temperature and environment effects, which promote competing interactions and alter their static gas-phase properties. Recently, neural network potentials (NNPs) trained on Density Functional Theory (DFT) data have become increasingly popular to simulate molecular phenomena in condensed phase with an accuracy comparable to ab initio methods. To date, most applications have centered on solid-state materials or fairly simple molecules made of a limited number of elements. Herein, we focus on the persistence and strength of chalcogen bonds involving a benzotelluradiazole in condensed phase. While the tellurium-containing heteroaromatic molecules are known to exhibit pronounced interactions with anions and lone pairs of different atoms, the relevance of competing intermolecular interactions, notably with the solvent, is complicated to monitor experimentally but also challenging to model at an accurate electronic structure level. Here, we train direct and baselined NNPs to reproduce hybrid DFT energies and forces in order to identify what are the most prevalent non-covalent interactions occurring in a solute-Cl$^-$-THF mixture. The simulations in explicit solvent highlight competition with chalcogen bonds formed with the solvent and the short-range directionality of the interaction with direct consequences for the molecular properties in the solution. The comparison with other potentials (e.g., AMOEBA, direct NNP and continuum solvent model) also demonstrates that baselined NNPs offer a reliable picture of the non-covalent interaction interplay occurring in solution.
Collapse
|
15
|
Li C, Voth GA. Using Machine Learning to Greatly Accelerate Path Integral Ab Initio Molecular Dynamics. J Chem Theory Comput 2022; 18:599-604. [PMID: 34982562 PMCID: PMC8864787 DOI: 10.1021/acs.jctc.1c01085] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Ab initio molecular dynamics (AIMD) has become one of the most popular and robust approaches for modeling complicated chemical, liquid, and material systems. However, the formidable computational cost often limits its widespread application in simulations of the largest-scale systems. The situation becomes even more severe in cases where the hydrogen nuclei may be better described as quantized particles using a path integral representation. Here, we present a computational approach that combines machine learning with recent advances in path integral contraction schemes, and we achieve a 2 orders of magnitude acceleration over direct path integral AIMD simulation while at the same time maintaining its accuracy.
Collapse
|
16
|
Xue Y, Wang JN, Hu W, Zheng J, Li Y, Pan X, Mo Y, Shao Y, Wang L, Mei Y. Affordable Ab Initio Path Integral for Thermodynamic Properties via Molecular Dynamics Simulations Using Semiempirical Reference Potential. J Phys Chem A 2021; 125:10677-10685. [PMID: 34894680 PMCID: PMC9108008 DOI: 10.1021/acs.jpca.1c07727] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
Path integral molecular dynamics (PIMD) is becoming a routinely applied method for incorporating the nuclear quantum effect in computer simulations. However, direct PIMD simulations at an ab initio level of theory are formidably expensive. Using the protonated 1,8-bis(dimethylamino)naphthalene molecule as an example, we show in this work that the computational expense for the intramolecular proton transfer between the two nitrogen atoms can be remarkably reduced by implementing the idea of reference-potential methods. The simulation time can be easily extended to a scale of nanoseconds while maintaining the accuracy on an ab initio level of theory for thermodynamic properties. In addition, postprocessing can be carried out in parallel on massive computer nodes. A 545-fold reduction in the total CPU time can be achieved in this way as compared to a direct PIMD simulation at the same ab initio level of theory.
Collapse
Affiliation(s)
- Yuanfei Xue
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Jia-Ning Wang
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China
| | - Wenxin Hu
- The Computer Center, School of Data Science & Engineering, East China Normal University, Shanghai 200062, China
| | - Jun Zheng
- The Computer Center, School of Data Science & Engineering, East China Normal University, Shanghai 200062, China
| | - Yongle Li
- Department of Physics, International Center of Quantum and Molecular Structure, and Shanghai Key Laboratory of High Temperature Superconductors, Shanghai University, Shanghai 200444, China
| | - Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Yan Mo
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China,NYU–ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China,Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, Norman, Oklahoma 73019, United States
| | - Lu Wang
- Department of Chemistry and Chemical Biology, Institute for Quantitative Biomedicine, Rutgers University, Piscataway, New Jersey 08854, United States
| | - Ye Mei
- State Key Laboratory of Precision Spectroscopy, School of Physics and Electronic Science, East China Normal University, Shanghai 200062, China,NYU–ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China,Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| |
Collapse
|
17
|
Gaigeot MP. Some opinions on MD-based vibrational spectroscopy of gas phase molecules and their assembly: An overview of what has been achieved and where to go. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2021; 260:119864. [PMID: 34052762 DOI: 10.1016/j.saa.2021.119864] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/14/2021] [Revised: 04/13/2021] [Accepted: 04/18/2021] [Indexed: 06/12/2023]
Abstract
We hereby review molecular dynamics simulations for anharmonic gas phase spectroscopy and provide some of our opinions of where the field is heading. With these new directions, the theoretical IR/Raman spectroscopy of large (bio)-molecular systems will be more easily achievable over longer time-scale MD trajectories for an increase in accuracy of the MD-IR and MD-Raman calculated spectra. With the new directions presented here, the high throughput 'decoding' of experimental IR/Raman spectra into 3D-structures should thus be possible, hence advancing e.g. the field of MS-IR for structural characterization by spectroscopy. We also review the assignment of vibrational spectra in terms of anharmonic molecular modes from the MD trajectories, and especially introduce our recent developments based on Graph Theory algorithms. Graph Theory algorithmic is also introduced in this review for the identification of the molecular 3D-structures sampled over MD trajectories.
Collapse
Affiliation(s)
- Marie-Pierre Gaigeot
- Université Paris-Saclay, Univ Evry, CNRS, LAMBE UMR8587, 91025 Evry-Courcouronnes, France.
| |
Collapse
|
18
|
Shepherd S, Lan J, Wilkins DM, Kapil V. Efficient Quantum Vibrational Spectroscopy of Water with High-Order Path Integrals: From Bulk to Interfaces. J Phys Chem Lett 2021; 12:9108-9114. [PMID: 34523941 DOI: 10.1021/acs.jpclett.1c02574] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Vibrational spectroscopy is key in probing the interplay between the structure and dynamics of aqueous systems. To map different regions of experimental spectra to the microscopic structure of a system, it is important to combine them with first-principles atomistic simulations that incorporate the quantum nature of nuclei. Here we show that the large cost of calculating the quantum vibrational spectra of aqueous systems can be dramatically reduced compared with standard path integral methods by using approximate quantum dynamics based on high-order path integrals. Together with state-of-the-art machine-learned electronic properties, our approach gives an excellent description not only of the infrared and Raman spectra of bulk water but also of the 2D correlation and the more challenging sum-frequency generation spectra of the water-air interface. This paves the way for understanding complex interfaces such as water encapsulated between or in contact with hydrophobic and hydrophilic materials through robust and inexpensive surface-sensitive and multidimensional spectra with first-principles accuracy.
Collapse
Affiliation(s)
- Sam Shepherd
- Atomistic Simulation Centre, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, Northern Ireland, United Kingdom
| | - Jinggang Lan
- Department of Chemistry, University of Zürich, Zürich 8057, Switzerland
| | - David M Wilkins
- Atomistic Simulation Centre, School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 1NN, Northern Ireland, United Kingdom
| | - Venkat Kapil
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW,United Kingdom
| |
Collapse
|
19
|
Mauger N, Plé T, Lagardère L, Bonella S, Mangaud É, Piquemal JP, Huppert S. Nuclear Quantum Effects in Liquid Water at Near Classical Computational Cost Using the Adaptive Quantum Thermal Bath. J Phys Chem Lett 2021; 12:8285-8291. [PMID: 34427440 DOI: 10.1021/acs.jpclett.1c01722] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We demonstrate the accuracy and efficiency of a recently introduced approach to account for nuclear quantum effects (NQEs) in molecular simulations: the adaptive quantum thermal bath (adQTB). In this method, zero-point energy is introduced through a generalized Langevin thermostat designed to precisely enforce the quantum fluctuation-dissipation theorem. We propose a refined adQTB algorithm with improved accuracy and report adQTB simulations of liquid water. Through extensive comparison with reference path integral calculations, we demonstrate that it provides excellent accuracy for a broad range of structural and thermodynamic observables as well as infrared vibrational spectra. The adQTB has a computational cost comparable to that of classical molecular dynamics, enabling simulations of up to millions of degrees of freedom.
Collapse
Affiliation(s)
- Nastasia Mauger
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Thomas Plé
- CNRS, Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588, 4 Place Jussieu, F-75005 Paris, France
| | - Louis Lagardère
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Sara Bonella
- CECAM Centre Européen de Calcul Atomique et Moléculaire, École Polytechnique Fédérale de Lausanne, Batochimie, Avenue Forel 2, 1015 Lausanne, Switzerland
| | - Étienne Mangaud
- CNRS, Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588, 4 Place Jussieu, F-75005 Paris, France
| | - Jean-Philip Piquemal
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
- Institut Universitaire de France, 75005 Paris, France
- Department of Biomedical Engineering, The University of Texas at Austin, Austin, Texas 78712, United States
| | - Simon Huppert
- CNRS, Sorbonne Université, Institut des NanoSciences de Paris, UMR 7588, 4 Place Jussieu, F-75005 Paris, France
| |
Collapse
|
20
|
Yao Y, Kanai Y. Nuclear Quantum Effect and Its Temperature Dependence in Liquid Water from Random Phase Approximation via Artificial Neural Network. J Phys Chem Lett 2021; 12:6354-6362. [PMID: 34231366 DOI: 10.1021/acs.jpclett.1c01566] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
We report structural and dynamical properties of liquid water described by the random phase approximation (RPA) correlation together with the exact exchange energy (EXX) within density functional theory. By utilizing thermostated ring polymer molecular dynamics, we examine the nuclear quantum effects and their temperature dependence. We circumvent the computational limitation of performing direct first-principles molecular dynamics simulation at this high level of electronic structure theory by adapting an artificial neural network model. We show that the EXX+RPA level of theory accurately describes liquid water in terms of both dynamical and structural properties.
Collapse
Affiliation(s)
- Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Durham, North Carolina 27599, United States
- Department of Mechanical Engineering and Materials Science, Duke University, Durham, North Carolina 27708, United States
| | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Durham, North Carolina 27599, United States
| |
Collapse
|
21
|
Imbalzano G, Zhuang Y, Kapil V, Rossi K, Engel EA, Grasselli F, Ceriotti M. Uncertainty estimation for molecular dynamics and sampling. J Chem Phys 2021; 154:074102. [DOI: 10.1063/5.0036522] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Giulio Imbalzano
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Yongbin Zhuang
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Collaborative Innovation Center of Chemistry for Energy Materials, Xiamen University, Xiamen 361005, China
| | - Venkat Kapil
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Kevin Rossi
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
- Laboratory of Nanochemistry for Energy, ISIC, École Polytechnique Fédérale de Lausanne, 1950 Sion, Switzerland
| | - Edgar A. Engel
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Federico Grasselli
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
22
|
Pure non-local machine-learned density functional theory for electron correlation. Nat Commun 2021; 12:344. [PMID: 33436595 PMCID: PMC7804195 DOI: 10.1038/s41467-020-20471-y] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Accepted: 12/02/2020] [Indexed: 01/02/2023] Open
Abstract
Density-functional theory (DFT) is a rigorous and (in principle) exact framework for the description of the ground state properties of atoms, molecules and solids based on their electron density. While computationally efficient density-functional approximations (DFAs) have become essential tools in computational chemistry, their (semi-)local treatment of electron correlation has a number of well-known pathologies, e.g. related to electron self-interaction. Here, we present a type of machine-learning (ML) based DFA (termed Kernel Density Functional Approximation, KDFA) that is pure, non-local and transferable, and can be efficiently trained with fully quantitative reference methods. The functionals retain the mean-field computational cost of common DFAs and are shown to be applicable to non-covalent, ionic and covalent interactions, as well as across different system sizes. We demonstrate their remarkable possibilities by computing the free energy surface for the protonated water dimer at hitherto unfeasible gold-standard coupled cluster quality on a single commodity workstation. Semilocal density functionals, while computationally efficient, do not account for non-local correlation. Here, the authors propose a machine-learning approach to DFT that leads to non-local and transferable functionals applicable to non-covalent, ionic and covalent interactions across system of different sizes.
Collapse
|
23
|
Kamat K, Guo R, Reutzel-Edens SM, Price SL, Peters B. Diabat method for polymorph free energies: Extension to molecular crystals. J Chem Phys 2020; 153:244105. [PMID: 33380078 DOI: 10.1063/5.0024727] [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
Lattice-switch Monte Carlo and the related diabat methods have emerged as efficient and accurate ways to compute free energy differences between polymorphs. In this work, we introduce a one-to-one mapping from the reference positions and displacements in one molecular crystal to the positions and displacements in another. Two features of the mapping facilitate lattice-switch Monte Carlo and related diabat methods for computing polymorph free energy differences. First, the mapping is unitary so that its Jacobian does not complicate the free energy calculations. Second, the mapping is easily implemented for molecular crystals of arbitrary complexity. We demonstrate the mapping by computing free energy differences between polymorphs of benzene and carbamazepine. Free energy calculations for thermodynamic cycles, each involving three independently computed polymorph free energy differences, all return to the starting free energy with a high degree of precision. The calculations thus provide a force field independent validation of the method and allow us to estimate the precision of the individual free energy differences.
Collapse
Affiliation(s)
- Kartik Kamat
- Department of Chemical Engineering, University of California-Santa Barbara, Santa Barbara, California 93106, USA
| | - Rui Guo
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Susan M Reutzel-Edens
- Small Molecule Design and Development, Eli Lilly and Company, Indianapolis, Indiana 46285, USA
| | - Sarah L Price
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Baron Peters
- Department of Chemical and Biomolecular Engineering, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
24
|
Brieuc F, Schran C, Uhl F, Forbert H, Marx D. Converged quantum simulations of reactive solutes in superfluid helium: The Bochum perspective. J Chem Phys 2020; 152:210901. [DOI: 10.1063/5.0008309] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Affiliation(s)
- Fabien Brieuc
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Christoph Schran
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Felix Uhl
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Harald Forbert
- Center for Solvation Science ZEMOS, Ruhr-Universität Bochum, 44780 Bochum, Germany
| | - Dominik Marx
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, 44780 Bochum, Germany
| |
Collapse
|
25
|
Poltavsky I, Kapil V, Ceriotti M, Kim KS, Tkatchenko A. Accurate Description of Nuclear Quantum Effects with High-Order Perturbed Path Integrals (HOPPI). J Chem Theory Comput 2020; 16:1128-1135. [PMID: 31913625 DOI: 10.1021/acs.jctc.9b00881] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Imaginary time path-integral (PI) simulations that account for nuclear quantum effects (NQE) beyond the harmonic approximation are increasingly employed together with modern electronic-structure calculations. Existing PI methods are applicable to molecules, liquids, and solids; however, the computational cost of such simulations increases dramatically with decreasing temperature. To address this challenge, here, we propose to combine high-order PI factorization with perturbation theory (PT). Already for conventional second-order PI simulations, the PT ansatz increases the accuracy 2-fold compared to fourth-order schemes with the same settings. In turn, applying PT to high-order path integrals (HOPI) further improves the efficiency of simulations for molecular and condensed matter systems especially at low temperatures. We present results for bulk liquid water, the aspirin molecule, and the CH5+ molecule. Perturbed HOPI simulations remain both efficient and accurate down to 20 K and provide a convenient method to estimate the convergence of quantum-mechanical observables.
Collapse
Affiliation(s)
- Igor Poltavsky
- Physics and Materials Science Research Unit , University of Luxembourg , L-1511 Luxembourg City , Luxembourg
| | - Venkat Kapil
- Laboratory of Computational Science and Modelling, Institute of Materials , Ecole Polytechnique Fédérale de Lausanne , Lausanne , Switzerland
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials , Ecole Polytechnique Fédérale de Lausanne , Lausanne , Switzerland
| | - Kwang S Kim
- Department of Chemistry, School of Natural Science , Ulsan National Institute of Science and Technology , Ulsan 44919 , Korea
| | - Alexandre Tkatchenko
- Physics and Materials Science Research Unit , University of Luxembourg , L-1511 Luxembourg City , Luxembourg
| |
Collapse
|
26
|
Liu C, Piquemal JP, Ren P. Implementation of Geometry-Dependent Charge Flux into the Polarizable AMOEBA+ Potential. J Phys Chem Lett 2020; 11:419-426. [PMID: 31865706 PMCID: PMC7384396 DOI: 10.1021/acs.jpclett.9b03489] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Molecular dynamics (MD) simulations employing classical force fields (FFs) have been widely used to model molecular systems. The important ingredient of the current FFs, atomic charge, remains fixed during MD simulations despite the atomic environment or local geometry changes. This approximation hinders the transferability of the potential being used in multiple phases. Here we implement a geometry-dependent charge flux (GDCF) model into the multipole-based AMOEBA+ polarizable potential. The CF in the current work explicitly depends on the local geometry (bond and angle) of the molecule. To our knowledge, this is the first study that derives energy and force expressions due to GDCF in a multipole-based polarizable FF framework. Due to the inclusion of GDCF, the AMOEBA+ water model is noticeably improved in terms of describing the monomer properties, cluster binding/interaction energy, and a variety of liquid properties, including the infrared spectra that previous flexible water models were not able to capture.
Collapse
Affiliation(s)
- Chengwen Liu
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
| | - Jean-Philip Piquemal
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
- Laboratoire de Chimie Théorique , Sorbonne Université, UMR7616 CNRS , 75252 Paris , France
- Institut Universitaire de France , 75005 , Paris , France
| | - Pengyu Ren
- Department of Biomedical Engineering , The University of Texas at Austin , Austin , Texas 78712 , United States
| |
Collapse
|
27
|
Wieme J, Vandenbrande S, Lamaire A, Kapil V, Vanduyfhuys L, Van Speybroeck V. Thermal Engineering of Metal-Organic Frameworks for Adsorption Applications: A Molecular Simulation Perspective. ACS APPLIED MATERIALS & INTERFACES 2019; 11:38697-38707. [PMID: 31556593 PMCID: PMC6818952 DOI: 10.1021/acsami.9b12533] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/17/2019] [Accepted: 09/26/2019] [Indexed: 05/29/2023]
Abstract
Thermal engineering of metal-organic frameworks for adsorption-based applications is very topical in view of their industrial potential, in particular, since heat management and thermal stability have been identified as important obstacles. Hence, a fundamental understanding of the structural and chemical features underpinning their intrinsic thermal properties is highly sought-after. Herein, we investigate the nanoscale behavior of a diverse set of frameworks using molecular simulation techniques and critically compare properties such as thermal conductivity, heat capacity, and thermal expansion with other classes of materials. Furthermore, we propose a hypothetical thermodynamic cycle to estimate the temperature rise associated with adsorption for the most important greenhouse and energy-related gases (CO2 and CH4). This macroscopic response on the heat of adsorption connects the intrinsic thermal properties with the adsorption properties and allows us to evaluate their importance.
Collapse
Affiliation(s)
- Jelle Wieme
- Center for Molecular
Modeling, Ghent University, Tech Lane Ghent Science Park Campus
A, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Steven Vandenbrande
- Center for Molecular
Modeling, Ghent University, Tech Lane Ghent Science Park Campus
A, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Aran Lamaire
- Center for Molecular
Modeling, Ghent University, Tech Lane Ghent Science Park Campus
A, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Venkat Kapil
- Laboratory
of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Louis Vanduyfhuys
- Center for Molecular
Modeling, Ghent University, Tech Lane Ghent Science Park Campus
A, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Veronique Van Speybroeck
- Center for Molecular
Modeling, Ghent University, Tech Lane Ghent Science Park Campus
A, Technologiepark 46, 9052 Zwijnaarde, Belgium
| |
Collapse
|
28
|
Ko HY, Zhang L, Santra B, Wang H, E W, DiStasio Jr RA, Car R. Isotope effects in liquid water via deep potential molecular dynamics. Mol Phys 2019. [DOI: 10.1080/00268976.2019.1652366] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Affiliation(s)
- Hsin-Yu Ko
- Department of Chemistry, Princeton University, Princeton, NJ, USA
| | - Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA
| | - Biswajit Santra
- Department of Physics, Temple University, Philadelphia, PA, USA
| | - Han Wang
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing, People's Republic of China
| | - Weinan E
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA
- Department of Mathematics, Princeton University, Princeton, NJ, USA
| | | | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, NJ, USA
- Program in Applied and Computational Mathematics, Princeton University, Princeton, NJ, USA
- Department of Physics and Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, NJ, USA
| |
Collapse
|
29
|
Clark T, Heske J, Kühne TD. Opposing Electronic and Nuclear Quantum Effects on Hydrogen Bonds in H 2 O and D 2 O. Chemphyschem 2019; 20:2461-2465. [PMID: 31449714 PMCID: PMC6790677 DOI: 10.1002/cphc.201900839] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2019] [Indexed: 11/09/2022]
Abstract
The effect of extending the O-H bond length(s) in water on the hydrogen-bonding strength has been investigated using static ab initio molecular orbital calculations. The "polar flattening" effect that causes a slight σ-hole to form on hydrogen atoms is strengthened when the bond is stretched, so that the σ-hole becomes more positive and hydrogen bonding stronger. In opposition to this electronic effect, path-integral ab initio molecular-dynamics simulations show that the nuclear quantum effect weakens the hydrogen bond in the water dimer. Thus, static electronic effects strengthen the hydrogen bond in H2 O relative to D2 O, whereas nuclear quantum effects weaken it. These quantum fluctuations are stronger for the water dimer than in bulk water.
Collapse
Affiliation(s)
- Timothy Clark
- Computer-Chemie-Centrum, Department Chemie und PharmazieFriedrich-Alexander-Universität Erlangen-NürnbergNägelsbachstr. 2591052ErlangenGermany
| | - Julian Heske
- Chair of Theoretical Chemistry, Dynamics of Condensed Matter, Department of ChemistryUniversity of PaderbornWarburger Str. 10033098PaderbornGermany
| | - Thomas D. Kühne
- Chair of Theoretical Chemistry, Dynamics of Condensed Matter, Department of ChemistryUniversity of PaderbornWarburger Str. 10033098PaderbornGermany
| |
Collapse
|
30
|
Kapil V, Wieme J, Vandenbrande S, Lamaire A, Van Speybroeck V, Ceriotti M. Modeling the Structural and Thermal Properties of Loaded Metal–Organic Frameworks. An Interplay of Quantum and Anharmonic Fluctuations. J Chem Theory Comput 2019; 15:3237-3249. [DOI: 10.1021/acs.jctc.8b01297] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Venkat Kapil
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Jelle Wieme
- Center for Molecular Modeling, Ghent University, Tech Lane Ghent Science Park Campus A, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Steven Vandenbrande
- Center for Molecular Modeling, Ghent University, Tech Lane Ghent Science Park Campus A, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Aran Lamaire
- Center for Molecular Modeling, Ghent University, Tech Lane Ghent Science Park Campus A, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Veronique Van Speybroeck
- Center for Molecular Modeling, Ghent University, Tech Lane Ghent Science Park Campus A, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| |
Collapse
|
31
|
Affiliation(s)
- Wei Fang
- School of Physics and Collaborative Innovation Centre of Quantum Matter, Peking University, Beijing, People's Republic of China
- Thomas Young Centre, London Centre for Nanotechnology, and Department of Physics and Astronomy, University College London, London, UK
- Laboratory of Physical Chemistry, ETH Zurich, Zurich, Switzerland
| | - Ji Chen
- Department of Electronic Structure Theory, Max Plank Institute for Solid State Research, Stuttgart, Germany
| | - Yexin Feng
- School of Physics and Electronics, Hunan University, Changsha, People's Republic of China
| | - Xin-Zheng Li
- School of Physics and Collaborative Innovation Centre of Quantum Matter, Peking University, Beijing, People's Republic of China
- State Key Laboratory for Artificial Microstructure and Mesoscopic Physics, Peking University, Beijing, People's Republic of China
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology, and Department of Physics and Astronomy, University College London, London, UK
| |
Collapse
|
32
|
Lamaire A, Wieme J, Rogge SMJ, Waroquier M, Van Speybroeck V. On the importance of anharmonicities and nuclear quantum effects in modelling the structural properties and thermal expansion of MOF-5. J Chem Phys 2019; 150:094503. [DOI: 10.1063/1.5085649] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Affiliation(s)
- Aran Lamaire
- Center for Molecular Modeling (CMM), Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Jelle Wieme
- Center for Molecular Modeling (CMM), Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Sven M. J. Rogge
- Center for Molecular Modeling (CMM), Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Michel Waroquier
- Center for Molecular Modeling (CMM), Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| | - Veronique Van Speybroeck
- Center for Molecular Modeling (CMM), Ghent University, Technologiepark 46, 9052 Zwijnaarde, Belgium
| |
Collapse
|
33
|
Elton DC, Fritz M, Fernández-Serra M. Using a monomer potential energy surface to perform approximate path integral molecular dynamics simulation of ab initio water at near-zero added cost. Phys Chem Chem Phys 2018; 21:409-417. [PMID: 30534683 DOI: 10.1039/c8cp06077k] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
It is now established that nuclear quantum motion plays an important role in determining water's hydrogen bonding, structure, and dynamics. Such effects are important to include in density functional theory (DFT) based molecular dynamics simulation of water. The standard way of treating nuclear quantum effects, path integral molecular dynamics (PIMD), multiplies the number of energy/force calculations by the number of beads required. In this work we introduce a method whereby PIMD can be incorporated into a DFT simulation with little extra cost and little loss in accuracy. The method is based on the many body expansion of the energy and has the benefit of including a monomer level correction to the DFT energy. Our method calculates intramolecular forces using the highly accurate monomer potential energy surface developed by Partridge-Schwenke, which is cheap to evaluate. Intermolecular forces and energies are calculated with DFT only once per timestep using the centroid positions. We show how our method may be used in conjunction with a multiple time step algorithm for an additional speedup and how it relates to ring polymer contraction and other schemes that have been introduced recently to speed up PIMD simulations. We show that our method, which we call "monomer PIMD", correctly captures changes in the structure of water found in a full PIMD simulation but at much lower computational cost.
Collapse
Affiliation(s)
- Daniel C Elton
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794-3800, USA.
| | | | | |
Collapse
|
34
|
Ruiz Pestana L, Marsalek O, Markland TE, Head-Gordon T. The Quest for Accurate Liquid Water Properties from First Principles. J Phys Chem Lett 2018; 9:5009-5016. [PMID: 30118601 DOI: 10.1021/acs.jpclett.8b02400] [Citation(s) in RCA: 53] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Developing accurate ab initio molecular dynamics (AIMD) models that capture both electronic reorganization and nuclear quantum effects associated with hydrogen bonding is key to quantitative understanding of bulk water and its anomalies as well as its role as a universal solvent. For condensed phase simulations, AIMD has typically relied on the generalized gradient approximation (GGA) of density functional theory (DFT) as the underlying model chemistry for the potential energy surface, with nuclear quantum effects (NQEs) sometimes modeled by performing classical molecular dynamics simulations at elevated temperatures. Here we show that the properties of liquid water obtained from the meta-GGA B97M-rV functional, when evaluated using accelerated path integral molecular dynamics simulations, display accuracy comparable to a computationally expensive dispersion-corrected hybrid functional, revPBE0-D3. We show that the meta-GGA DFT functional reproduces bulk water properties including radial distribution functions, self-diffusion coefficients, and infrared spectra with comparable accuracy of a much more expensive hybrid functional. This work demonstrates that the underlying quality of a good DFT functional requires evaluation with quantum nuclei and that high-temperature simulations are a poor proxy for properly treating NQEs.
Collapse
Affiliation(s)
| | - Ondrej Marsalek
- Charles University , Faculty of Mathematics and Physics , Ke Karlovu 3 , 121 16 Prague 2 , Czech Republic
| | - Thomas E Markland
- Department of Chemistry , Stanford University , Stanford , California 94305 , United States
| | | |
Collapse
|
35
|
Drużbicki K, Krzystyniak M, Hollas D, Kapil V, Slavíček P, Romanelli G, Fernandez-Alonso F. Hydrogen dynamics in solid formic acid: insights from simulations with quantum colored-noise thermostats. ACTA ACUST UNITED AC 2018. [DOI: 10.1088/1742-6596/1055/1/012003] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
36
|
Litman Y, Donadio D, Ceriotti M, Rossi M. Decisive role of nuclear quantum effects on surface mediated water dissociation at finite temperature. J Chem Phys 2018; 148:102320. [PMID: 29544260 DOI: 10.1063/1.5002537] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
Water molecules adsorbed on inorganic substrates play an important role in several technological applications. In the presence of light atoms in adsorbates, nuclear quantum effects (NQEs) influence the structural stability and the dynamical properties of these systems. In this work, we explore the impact of NQEs on the dissociation of water wires on stepped Pt(221) surfaces. By performing ab initio molecular dynamics simulations with van der Waals corrected density functional theory, we note that several competing minima for both intact and dissociated structures are accessible at finite temperatures, making it important to assess whether harmonic estimates of the quantum free energy are sufficient to determine the relative stability of the different states. We thus perform ab initio path integral molecular dynamics (PIMD) in order to calculate these contributions taking into account the conformational entropy and anharmonicities at finite temperatures. We propose that when adsorption is weak and NQEs on the substrate are negligible, PIMD simulations can be performed through a simple partition of the system, resulting in considerable computational savings. We then calculate the full contribution of NQEs to the free energies, including also anharmonic terms. We find that they result in an increase of up to 20% of the quantum contribution to the dissociation free energy compared with the harmonic estimates. We also find that the dissociation process has a negligible contribution from tunneling but is dominated by zero point energies, which can enhance the rate of dissociation by three orders of magnitude. Finally we highlight how both temperature and NQEs indirectly impact dipoles and the redistribution of electron density, causing work function changes of up to 0.4 eV with respect to static estimates. This quantitative determination of the change in the work function provides a possible approach to determine experimentally the most stable configurations of water oligomers on the stepped surfaces.
Collapse
Affiliation(s)
- Yair Litman
- Fritz Haber Institute of the Max Planck Society, Faradayweg 4-6, 14195 Berlin, Germany
| | - Davide Donadio
- Department of Chemistry, University of California Davis, One Shields Ave., Davis, California 95616, USA
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Mariana Rossi
- Fritz Haber Institute of the Max Planck Society, Faradayweg 4-6, 14195 Berlin, Germany
| |
Collapse
|
37
|
Kapil V, Cuzzocrea A, Ceriotti M. Anisotropy of the Proton Momentum Distribution in Water. J Phys Chem B 2018; 122:6048-6054. [DOI: 10.1021/acs.jpcb.8b03896] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Affiliation(s)
- Venkat Kapil
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| | - Alice Cuzzocrea
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
| |
Collapse
|
38
|
|
39
|
Buxton SJ, Habershon S. Accelerated path-integral simulations using ring-polymer interpolation. J Chem Phys 2017; 147:224107. [DOI: 10.1063/1.5006465] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Affiliation(s)
- Samuel J. Buxton
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
40
|
Unravelling the influence of quantum proton delocalization on electronic charge transfer through the hydrogen bond. Chem Phys Lett 2017. [DOI: 10.1016/j.cplett.2017.04.034] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
41
|
Mouhat F, Sorella S, Vuilleumier R, Saitta AM, Casula M. Fully Quantum Description of the Zundel Ion: Combining Variational Quantum Monte Carlo with Path Integral Langevin Dynamics. J Chem Theory Comput 2017; 13:2400-2417. [DOI: 10.1021/acs.jctc.7b00017] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Félix Mouhat
- IMPMC, Sorbonne Universités, Université Pierre et Marie Curie Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, 4 place Jussieu, 75252 Paris, France
| | - Sandro Sorella
- International
School for Advanced Studies (SISSA), Via Bonomea 26, 34136 Trieste, Italy
- INFM Democritos
National Simulation Center, 34151 Trieste, Italy
| | - Rodolphe Vuilleumier
- PASTEUR,
Département de chimie, École normale supérieure,
UPMC Univ Paris 06, CNRS, PSL Research University, 75005 Paris, France
- Sorbonne Universités, UPMC Univ Paris 06, École normale supérieure, CNRS, Processus d’activation
sélective par transfert d’énergie uni-électronique
ou radiatif (PASTEUR), 75005 Paris, France
| | - Antonino Marco Saitta
- IMPMC, Sorbonne Universités, Université Pierre et Marie Curie Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, 4 place Jussieu, 75252 Paris, France
| | - Michele Casula
- IMPMC, Sorbonne Universités, Université Pierre et Marie Curie Paris 06, CNRS UMR 7590, IRD UMR 206, MNHN, 4 place Jussieu, 75252 Paris, France
| |
Collapse
|
42
|
Marsalek O, Markland TE. Quantum Dynamics and Spectroscopy of Ab Initio Liquid Water: The Interplay of Nuclear and Electronic Quantum Effects. J Phys Chem Lett 2017; 8:1545-1551. [PMID: 28296422 DOI: 10.1021/acs.jpclett.7b00391] [Citation(s) in RCA: 137] [Impact Index Per Article: 19.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Understanding the reactivity and spectroscopy of aqueous solutions at the atomistic level is crucial for the elucidation and design of chemical processes. However, the simulation of these systems requires addressing the formidable challenges of treating the quantum nature of both the electrons and nuclei. Exploiting our recently developed methods that provide acceleration by up to 2 orders of magnitude, we combine path integral simulations with on-the-fly evaluation of the electronic structure at the hybrid density functional theory level to capture the interplay between nuclear quantum effects and the electronic surface. Here we show that this combination provides accurate structure and dynamics, including the full infrared and Raman spectra of liquid water. This allows us to demonstrate and explain the failings of lower-level density functionals for dynamics and vibrational spectroscopy when the nuclei are treated quantum mechanically. These insights thus provide a foundation for the reliable investigation of spectroscopy and reactivity in aqueous environments.
Collapse
Affiliation(s)
- Ondrej Marsalek
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| | - Thomas E Markland
- Department of Chemistry, Stanford University , Stanford, California 94305, United States
| |
Collapse
|
43
|
Rybkin VV, VandeVondele J. Nuclear Quantum Effects on Aqueous Electron Attachment and Redox Properties. J Phys Chem Lett 2017; 8:1424-1428. [PMID: 28296416 DOI: 10.1021/acs.jpclett.7b00386] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/06/2023]
Abstract
Nuclear quantum effects (NQEs) on the reduction and oxidation properties of small aqueous species (CO2, HO2, and O2) are quantified and rationalized by first-principles molecular dynamics and thermodynamic integration. Vertical electron attachment, or electron affinity, and detachment energies (VEA and VDE) are strongly affected by NQEs, decreasing in absolute value by 0.3 eV going from a classical to a quantum description of the nuclei. The effect is attributed to NQEs that lessen the solvent response upon oxidation/reduction. The reduction of solvent reorganization energy is expected to be general for small solutes in water. In the thermodynamic integral that yields the free energy of oxidation/reduction, these large changes enter with opposite sign, and only a small net effect (0.1 eV) remains. This is not obvious for CO2, where the integrand is strongly influenced by NQEs due to the onset of interaction of the reduced orbital with the conduction band of the liquid during thermodynamic integration. We conclude that NQEs might not have to be included in the computation of redox potentials, unless high accuracy is needed, but are important for VEA and VDE calculations.
Collapse
Affiliation(s)
- Vladimir V Rybkin
- Nanoscale Simulations, Department of Materials, ETH Zürich , Wolfgang-Pauli-Str. 27, CH-8093 Zürich, Switzerland
| | - Joost VandeVondele
- Nanoscale Simulations, Department of Materials, ETH Zürich , Wolfgang-Pauli-Str. 27, CH-8093 Zürich, Switzerland
| |
Collapse
|
44
|
Hirshberg B, Sagiv L, Gerber RB. Approximate Quantum Dynamics using Ab Initio Classical Separable Potentials: Spectroscopic Applications. J Chem Theory Comput 2017; 13:982-991. [DOI: 10.1021/acs.jctc.6b01129] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Barak Hirshberg
- Institute
of Chemistry and the Fritz Haber Center for Molecular Dynamics, The Hebrew University, Jerusalem 9190401, Israel
| | - Lior Sagiv
- Institute
of Chemistry and the Fritz Haber Center for Molecular Dynamics, The Hebrew University, Jerusalem 9190401, Israel
| | - R. Benny Gerber
- Institute
of Chemistry and the Fritz Haber Center for Molecular Dynamics, The Hebrew University, Jerusalem 9190401, Israel
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| |
Collapse
|
45
|
Affiliation(s)
- Venkat Kapil
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| | - Jörg Behler
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, Bochum, Germany
| | - Michele Ceriotti
- Laboratory of Computational Science and Modelling, Institute of Materials, Ecole Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
| |
Collapse
|
46
|
Suleimanov YV, Aoiz FJ, Guo H. Chemical Reaction Rate Coefficients from Ring Polymer Molecular Dynamics: Theory and Practical Applications. J Phys Chem A 2016; 120:8488-8502. [DOI: 10.1021/acs.jpca.6b07140] [Citation(s) in RCA: 92] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Yury V. Suleimanov
- Computation-based Science
and Technology Research Center, Cyprus Institute, 20 Kavafi Street, Nicosia 2121, Cyprus
- Department
of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - F. Javier Aoiz
- Departamento de Química
Física I, Facultad de CC. Químicas, Universidad Complutense de Madrid, 28040 Madrid, Spain
| | - Hua Guo
- Department of Chemistry and
Chemical Biology, University of New Mexico, Albuquerque, New Mexico 87131, United States
| |
Collapse
|
47
|
Rossi M, Gasparotto P, Ceriotti M. Anharmonic and Quantum Fluctuations in Molecular Crystals: A First-Principles Study of the Stability of Paracetamol. PHYSICAL REVIEW LETTERS 2016; 117:115702. [PMID: 27661700 DOI: 10.1103/physrevlett.117.115702] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/28/2016] [Indexed: 06/06/2023]
Abstract
Molecular crystals often exist in multiple competing polymorphs, showing significantly different physicochemical properties. Computational crystal structure prediction is key to interpret and guide the search for the most stable or useful form, a real challenge due to the combinatorial search space, and the complex interplay of subtle effects that work together to determine the relative stability of different structures. Here we take a comprehensive approach based on different flavors of thermodynamic integration in order to estimate all contributions to the free energies of these systems with density-functional theory, including the oft-neglected anharmonic contributions and nuclear quantum effects. We take the two main stable forms of paracetamol as a paradigmatic example. We find that anharmonic contributions, different descriptions of van der Waals interactions, and nuclear quantum effects all matter to quantitatively determine the stability of different phases. Our analysis highlights the many challenges inherent in the development of a quantitative and predictive framework to model molecular crystals. However, it also indicates which of the components of the free energy can benefit from a cancellation of errors that can redeem the predictive power of approximate models, and suggests simple steps that could be taken to improve the reliability of ab initio crystal structure prediction.
Collapse
Affiliation(s)
- Mariana Rossi
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Piero Gasparotto
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Michele Ceriotti
- Laboratory of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| |
Collapse
|
48
|
Tikhonov DS, Otlyotov AA, Rybkin VV. The effect of molecular dynamics sampling on the calculated observable gas-phase structures. Phys Chem Chem Phys 2016; 18:18237-45. [PMID: 27331660 DOI: 10.1039/c6cp02973f] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
In this study, we compare the performance of various ab initio molecular dynamics (MD) sampling methods for the calculation of the observable vibrationally-averaged gas-phase structures of benzene, naphthalene and anthracene molecules. Nose-Hoover (NH), canonical and quantum generalized-Langevin-equation (GLE) thermostats as well as the a posteriori quantum correction to the classical trajectories have been tested and compared to the accurate path-integral molecular dynamics (PIMD), static anharmonic vibrational calculations as well as to the experimental gas electron diffraction data. Classical sampling methods neglecting quantum effects (NH and canonical GLE thermostats) dramatically underestimate vibrational amplitudes for the bonded atom pairs, both C-H and C-C, the resulting radial distribution functions exhibit nonphysically narrow peaks. This deficiency is almost completely removed by taking the quantum effects on the nuclei into account. The quantum GLE thermostat and a posteriori correction to the canonical GLE and NH thermostatted trajectories capture most vibrational quantum effects and closely reproduce computationally expensive PIMD and experimental radial distribution functions. These methods are both computationally feasible and accurate and are therefore recommended for calculations of the observable gas-phase structures. A good performance of the quantum GLE thermostat for the gas-phase calculations is encouraging since its parameters have been originally fitted for the condensed-phase calculations. Very accurate molecular structures can be predicted by combining the equilibrium geometry obtained at a high level of electronic structure theory with vibrational amplitudes and corrections calculated using MD driven by a lower level of electronic structure theory.
Collapse
Affiliation(s)
- Denis S Tikhonov
- Universität Bielefeld, Lehrstuhl für Anorganische Chemie und Strukturchemie, Universitätsstrasse 25, 33615, Bielefeld, Germany.
| | | | | |
Collapse
|
49
|
Ceriotti M, Fang W, Kusalik PG, McKenzie RH, Michaelides A, Morales MA, Markland TE. Nuclear Quantum Effects in Water and Aqueous Systems: Experiment, Theory, and Current Challenges. Chem Rev 2016; 116:7529-50. [DOI: 10.1021/acs.chemrev.5b00674] [Citation(s) in RCA: 339] [Impact Index Per Article: 42.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Michele Ceriotti
- Laboratory
of Computational Science and Modeling, Institute of Materials, École Polytechnique Fédérale de Lausanne, 1015 Lausanne, Switzerland
| | - Wei Fang
- Thomas
Young Centre, London Centre for Nanotechnology and Department of Physics
and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Peter G. Kusalik
- Department
of Chemistry, University of Calgary, 2500 University Drive NW, Calgary, Alberta T2N 1N4, Canada
| | - Ross H. McKenzie
- School
of Mathematics and Physics, University of Queensland, Brisbane, 4072 Queensland Australia
| | - Angelos Michaelides
- Thomas
Young Centre, London Centre for Nanotechnology and Department of Physics
and Astronomy, University College London, London WC1E 6BT, United Kingdom
| | - Miguel A. Morales
- Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Thomas E. Markland
- Department
of Chemistry, Stanford University, 333 Campus Drive, Stanford, California 94305, United States
| |
Collapse
|
50
|
John C, Spura T, Habershon S, Kühne TD. Quantum ring-polymer contraction method: Including nuclear quantum effects at no additional computational cost in comparison to ab initio molecular dynamics. Phys Rev E 2016; 93:043305. [PMID: 27176426 DOI: 10.1103/physreve.93.043305] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2016] [Indexed: 06/05/2023]
Abstract
We present a simple and accurate computational method which facilitates ab initio path-integral molecular dynamics simulations, where the quantum-mechanical nature of the nuclei is explicitly taken into account, at essentially no additional computational cost in comparison to the corresponding calculation using classical nuclei. The predictive power of the proposed quantum ring-polymer contraction method is demonstrated by computing various static and dynamic properties of liquid water at ambient conditions using density functional theory. This development will enable routine inclusion of nuclear quantum effects in ab initio molecular dynamics simulations of condensed-phase systems.
Collapse
Affiliation(s)
- Christopher John
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| | - Thomas Spura
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Thomas D Kühne
- Dynamics of Condensed Matter, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany and Paderborn Center for Parallel Computing and Institute for Lightweight Design, Department of Chemistry, University of Paderborn, Warburger Strasse 100, D-33098 Paderborn, Germany
| |
Collapse
|