-
Accurate nuclear quantum statistics on machine-learned classical effective potentials
Authors:
Iryna Zaporozhets,
Félix Musil,
Venkat Kapil,
Cecilia Clementi
Abstract:
The contribution of nuclear quantum effects (NQEs) to the properties of various hydrogen-bound systems, including biomolecules, is increasingly recognized. Despite the development of many acceleration techniques, the computational overhead of incorporating NQEs in complex systems is sizable, particularly at low temperatures. In this work, we leverage deep learning and multiscale coarse-graining te…
▽ More
The contribution of nuclear quantum effects (NQEs) to the properties of various hydrogen-bound systems, including biomolecules, is increasingly recognized. Despite the development of many acceleration techniques, the computational overhead of incorporating NQEs in complex systems is sizable, particularly at low temperatures. In this work, we leverage deep learning and multiscale coarse-graining techniques to mitigate the computational burden of path integral molecular dynamics (PIMD). Specifically, we employ a machine-learned potential to accurately represent corrections to classical potentials, thereby significantly reducing the computational cost of simulating NQEs. We validate our approach using four distinct systems: Morse potential, Zundel cation, single water molecule, and bulk water. Our framework allows us to accurately compute position-dependent static properties, as demonstrated by the excellent agreement obtained between the machine-learned potential and computationally intensive PIMD calculations, even in the presence of strong NQEs. This approach opens the way to the development of transferable machine-learned potentials capable of accurately reproducing NQEs in a wide range of molecular systems.
△ Less
Submitted 3 July, 2024;
originally announced July 2024.
-
The Wetting of H$_2$O by CO$_2$
Authors:
Samuel G. H. Brookes,
Venkat Kapil,
Christoph Schran,
Angelos Michaelides
Abstract:
Biphasic interfaces are complex but fascinating regimes that display a number of properties distinct from those of the bulk. The CO$_2$-H$_2$O interface, in particular, has been the subject of a number of studies on account of its importance for the carbon life cycle as well as carbon capture and sequestration schemes. Despite this attention, there remain a number of open questions on the nature o…
▽ More
Biphasic interfaces are complex but fascinating regimes that display a number of properties distinct from those of the bulk. The CO$_2$-H$_2$O interface, in particular, has been the subject of a number of studies on account of its importance for the carbon life cycle as well as carbon capture and sequestration schemes. Despite this attention, there remain a number of open questions on the nature of the CO$_2$-H$_2$O interface, particularly concerning the interfacial tension and phase behavior of CO$_2$ at the interface. In this paper, we seek to address these ambiguities using ab initio-quality simulations. Harnessing the benefits of machine-learned potentials and enhanced statistical sampling methods, we present an ab initio-level description of the CO$_2$-H$_2$O interface. Interfacial tensions are predicted from 1-500 bar and found to be in close agreement with experiment at the pressures for which experimental data is available. Structural analyses indicate the build-up of an adsorbed, saturated CO$_2$ film forming at low pressure (20 bar) with properties similar to those of the bulk liquid, but preferential perpendicular alignment with respect to the interface. CO$_2$ monolayer build-up coincides with a reduced structuring of water molecules close to the interface. This study highlights the predictive nature of machine-learned potentials for complex macroscopic properties of biphasic interfaces, and the mechanistic insight obtained into carbon dioxide aggregation at the water interface is of high relevance for geoscience, climate research, and materials science.
△ Less
Submitted 21 June, 2024;
originally announced June 2024.
-
Data-efficient fine-tuning of foundational models for first-principles quality sublimation enthalpies
Authors:
Harveen Kaur,
Flaviano Della Pia,
Ilyes Batatia,
Xavier R. Advincula,
Benjamin X. Shi,
Jinggang Lan,
Gábor Csányi,
Angelos Michaelides,
Venkat Kapil
Abstract:
Calculating sublimation enthalpies of molecular crystal polymorphs is relevant to a wide range of technological applications. However, predicting these quantities at first-principles accuracy -- even with the aid of machine learning potentials -- is a challenge that requires sub-kJ/mol accuracy in the potential energy surface and finite-temperature sampling. We present an accurate and data-efficie…
▽ More
Calculating sublimation enthalpies of molecular crystal polymorphs is relevant to a wide range of technological applications. However, predicting these quantities at first-principles accuracy -- even with the aid of machine learning potentials -- is a challenge that requires sub-kJ/mol accuracy in the potential energy surface and finite-temperature sampling. We present an accurate and data-efficient protocol based on fine-tuning of the foundational MACE-MP-0 model and showcase its capabilities on sublimation enthalpies and physical properties of ice polymorphs. Our approach requires only a few tens of training structures to achieve sub-kJ/mol accuracy in the sublimation enthalpies and sub 1 % error in densities for polymorphs at finite temperature and pressure. Exploiting this data efficiency, we explore simulations of hexagonal ice at the random phase approximation level of theory at experimental temperatures and pressures, calculating its physical properties, like pair correlation function and density, with good agreement with experiments. Our approach provides a way forward for predicting the stability of molecular crystals at finite thermodynamic conditions with the accuracy of correlated electronic structure theory.
△ Less
Submitted 30 May, 2024;
originally announced May 2024.
-
i-PI 3.0: a flexible and efficient framework for advanced atomistic simulations
Authors:
Yair Litman,
Venkat Kapil,
Yotam M. Y. Feldman,
Davide Tisi,
Tomislav Begušić,
Karen Fidanyan,
Guillaume Fraux,
Jacob Higer,
Matthias Kellner,
Tao E. Li,
Eszter S. Pós,
Elia Stocco,
George Trenins,
Barak Hirshberg,
Mariana Rossi,
Michele Ceriotti
Abstract:
Atomic-scale simulations have progressed tremendously over the past decade, largely due to the availability of machine-learning interatomic potentials. These potentials combine the accuracy of electronic structure calculations with the ability to reach extensive length and time scales. The i-PI package facilitates integrating the latest developments in this field with advanced modeling techniques,…
▽ More
Atomic-scale simulations have progressed tremendously over the past decade, largely due to the availability of machine-learning interatomic potentials. These potentials combine the accuracy of electronic structure calculations with the ability to reach extensive length and time scales. The i-PI package facilitates integrating the latest developments in this field with advanced modeling techniques, thanks to a modular software architecture based on inter-process communication through a socket interface. The choice of Python for implementation facilitates rapid prototyping but can add computational overhead. In this new release, we carefully benchmarked and optimized i-PI for several common simulation scenarios, making such overhead negligible when i-PI is used to model systems up to tens of thousands of atoms using widely adopted machine learning interatomic potentials, such as Behler-Parinello, DeePMD and MACE neural networks. We also present the implementation of several new features, including an efficient algorithm to model bosonic and fermionic exchange, a framework for uncertainty quantification to be used in conjunction with machine-learning potentials, a communication infrastructure that allows deeper integration with electronic-driven simulations, and an approach to simulate coupled photon-nuclear dynamics in optical or plasmonic cavities.
△ Less
Submitted 10 July, 2024; v1 submitted 24 May, 2024;
originally announced May 2024.
-
A foundation model for atomistic materials chemistry
Authors:
Ilyes Batatia,
Philipp Benner,
Yuan Chiang,
Alin M. Elena,
Dávid P. Kovács,
Janosh Riebesell,
Xavier R. Advincula,
Mark Asta,
Matthew Avaylon,
William J. Baldwin,
Fabian Berger,
Noam Bernstein,
Arghya Bhowmik,
Samuel M. Blau,
Vlad Cărare,
James P. Darby,
Sandip De,
Flaviano Della Pia,
Volker L. Deringer,
Rokas Elijošius,
Zakariya El-Machachi,
Fabio Falcioni,
Edvin Fako,
Andrea C. Ferrari,
Annalena Genreith-Schriever
, et al. (51 additional authors not shown)
Abstract:
Machine-learned force fields have transformed the atomistic modelling of materials by enabling simulations of ab initio quality on unprecedented time and length scales. However, they are currently limited by: (i) the significant computational and human effort that must go into development and validation of potentials for each particular system of interest; and (ii) a general lack of transferabilit…
▽ More
Machine-learned force fields have transformed the atomistic modelling of materials by enabling simulations of ab initio quality on unprecedented time and length scales. However, they are currently limited by: (i) the significant computational and human effort that must go into development and validation of potentials for each particular system of interest; and (ii) a general lack of transferability from one chemical system to the next. Here, using the state-of-the-art MACE architecture we introduce a single general-purpose ML model, trained on a public database of 150k inorganic crystals, that is capable of running stable molecular dynamics on molecules and materials. We demonstrate the power of the MACE-MP-0 model - and its qualitative and at times quantitative accuracy - on a diverse set problems in the physical sciences, including the properties of solids, liquids, gases, chemical reactions, interfaces and even the dynamics of a small protein. The model can be applied out of the box and as a starting or "foundation model" for any atomistic system of interest and is thus a step towards democratising the revolution of ML force fields by lowering the barriers to entry.
△ Less
Submitted 1 March, 2024; v1 submitted 29 December, 2023;
originally announced January 2024.
-
MACE-OFF23: Transferable Machine Learning Force Fields for Organic Molecules
Authors:
Dávid Péter Kovács,
J. Harry Moore,
Nicholas J. Browning,
Ilyes Batatia,
Joshua T. Horton,
Venkat Kapil,
William C. Witt,
Ioan-Bogdan Magdău,
Daniel J. Cole,
Gábor Csányi
Abstract:
Classical empirical force fields have dominated biomolecular simulation for over 50 years. Although widely used in drug discovery, crystal structure prediction, and biomolecular dynamics, they generally lack the accuracy and transferability required for predictive modelling. In this paper, we introduce MACE-OFF23, a transferable force field for organic molecules created using state-of-the-art mach…
▽ More
Classical empirical force fields have dominated biomolecular simulation for over 50 years. Although widely used in drug discovery, crystal structure prediction, and biomolecular dynamics, they generally lack the accuracy and transferability required for predictive modelling. In this paper, we introduce MACE-OFF23, a transferable force field for organic molecules created using state-of-the-art machine learning technology and first-principles reference data computed with a high level of quantum mechanical theory. MACE-OFF23 demonstrates the remarkable capabilities of local, short-range models by accurately predicting a wide variety of gas and condensed phase properties of molecular systems. It produces accurate, easy-to-converge dihedral torsion scans of unseen molecules, as well as reliable descriptions of molecular crystals and liquids, including quantum nuclear effects. We further demonstrate the capabilities of MACE-OFF23 by determining free energy surfaces in explicit solvent, as well as the folding dynamics of peptides. Finally, we simulate a fully solvated small protein, observing accurate secondary structure and vibrational spectrum. These developments enable first-principles simulations of molecular systems for the broader chemistry community at high accuracy and low computational cost.
△ Less
Submitted 29 December, 2023; v1 submitted 23 December, 2023;
originally announced December 2023.
-
Quantum dynamics using path integral coarse-graining
Authors:
Félix Musil,
Iryna Zaporozhets,
Frank Noé,
Cecilia Clementi,
Venkat Kapil
Abstract:
Vibrational spectra of condensed and gas-phase systems containing light nuclei are influenced by their quantum-mechanical behaviour. The quantum dynamics of light nuclei can be approximated by the imaginary time path integral (PI) formulation, but still at a large computational cost that increases sharply with decreasing temperature. By leveraging advances in machine-learned coarse-graining, we de…
▽ More
Vibrational spectra of condensed and gas-phase systems containing light nuclei are influenced by their quantum-mechanical behaviour. The quantum dynamics of light nuclei can be approximated by the imaginary time path integral (PI) formulation, but still at a large computational cost that increases sharply with decreasing temperature. By leveraging advances in machine-learned coarse-graining, we develop a PI method with the reduced computational cost of a classical simulation. We also propose a simple temperature elevation scheme to significantly attenuate the artefacts of standard PI approaches and also eliminate the unfavourable temperature scaling of the computational cost.We illustrate the approach, by calculating vibrational spectra using standard models of water molecules and bulk water, demonstrating significant computational savings and dramatically improved accuracy compared to more expensive reference approaches. We believe that our simple, efficient and accurate method could enable routine calculations of vibrational spectra including nuclear quantum effects for a wide range of molecular systems.
△ Less
Submitted 23 September, 2022; v1 submitted 12 August, 2022;
originally announced August 2022.
-
General embedded cluster protocol for accurate modeling of oxygen vacancies in metal-oxides
Authors:
Benjamin Xu Shi,
Venkat Kapil,
Andrea Zen,
Ji Chen,
Ali Alavi,
Angelos Michaelides
Abstract:
The O vacancy (Ov) formation energy, $E_\textrm{Ov}$, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g. the generalized gradient approximation) are insufficient for modeling th…
▽ More
The O vacancy (Ov) formation energy, $E_\textrm{Ov}$, is an important property of a metal-oxide, governing its performance in applications such as fuel cells or heterogeneous catalysis. These defects are routinely studied with density functional theory (DFT). However, it is well-recognized that standard DFT formulations (e.g. the generalized gradient approximation) are insufficient for modeling the Ov, requiring higher levels of theory. The embedded cluster method offers a promising approach to compute $E_\textrm{Ov}$ accurately, giving access to all electronic structure methods. Central to this approach is the construction of quantum(-mechanically treated) clusters placed within suitable embedding environments. Unfortunately, current approaches to constructing the quantum clusters either require large system sizes, preventing application of high-level methods, or require significant manual input, preventing investigations of multiple systems simultaneously. In this work, we present a systematic and general quantum cluster design protocol that can determine small converged quantum clusters for studying the Ov in metal-oxides with accurate methods such as local coupled cluster with single, double and perturbative triple excitations [CCSD(T)]. We apply this protocol to study the Ov in the bulk and surface planes of rutile TiO2 and rocksalt MgO, producing the first accurate and well-converged determinations of $E_\textrm{Ov}$ with this method. These reference values are used to benchmark exchange-correlation functionals in DFT and we find that all studied functionals underestimate $E_\textrm{Ov}$, with the average error decreasing along the rungs of Jacob's ladder. This protocol is automatable for high-throughput calculations and can be generalized to study other point defects or adsorbates.
△ Less
Submitted 31 March, 2022; v1 submitted 9 February, 2022;
originally announced February 2022.
-
Efficient Quantum Vibrational Spectroscopy of Water with High-Order Path Integrals: from Bulk to Interfaces
Authors:
Sam Shepherd,
Jinggang Lan,
David M. Wilkins,
Venkat Kapil
Abstract:
Vibrational spectroscopy is key for probing the interplay between the structure and dynamics of aqueous systems. In order 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 quantum vibrati…
▽ More
Vibrational spectroscopy is key for probing the interplay between the structure and dynamics of aqueous systems. In order 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 quantum vibrational spectra of aqueous systems can be dramatically reduced compared to 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 of the infrared and Raman spectra of bulk water, but also of 2D correlation and 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.
△ Less
Submitted 6 August, 2021;
originally announced August 2021.
-
The importance of nuclear quantum effects for NMR crystallography
Authors:
Edgar A. Engel,
Venkat Kapil,
Michele Ceriotti
Abstract:
The resolving power of solid-state nuclear magnetic resonance (NMR) crystallography depends heavily on the accuracy of computational predictions of NMR chemical shieldings of candidate structures, which are usually taken to be local minima in the potential energy. To test the limits of this approximation, we systematically study the importance of finite-temperature and quantum nuclear fluctuations…
▽ More
The resolving power of solid-state nuclear magnetic resonance (NMR) crystallography depends heavily on the accuracy of computational predictions of NMR chemical shieldings of candidate structures, which are usually taken to be local minima in the potential energy. To test the limits of this approximation, we systematically study the importance of finite-temperature and quantum nuclear fluctuations for $^1$H, $^{13}$C, and $^{15}$N shieldings in polymorphs of three paradigmatic molecular crystals -- benzene, glycine, and succinic acid. The effect of quantum fluctuations is comparable to the typical errors of shielding predictions for static nuclei with respect to experiments, and their inclusion to improve the agreement with measurements, translating to more reliable assignment of the NMR spectra to the correct candidate structure. The use of integrated machine-learning models, trained on first-principles energies and shieldings, renders rigorous sampling of nuclear fluctuations affordable, setting a new standard for the calculations underlying NMR structure determinations.
△ Less
Submitted 9 January, 2022; v1 submitted 27 June, 2021;
originally announced June 2021.
-
Uncertainty estimation for molecular dynamics and sampling
Authors:
Giulio Imbalzano,
Yongbin Zhuang,
Venkat Kapil,
Kevin Rossi,
Edgar A. Engel,
Federico Grasselli,
Michele Ceriotti
Abstract:
Machine learning models have emerged as a very effective strategy to sidestep time-consuming electronic-structure calculations, enabling accurate simulations of greater size, time scale and complexity. Given the interpolative nature of these models, the reliability of predictions depends on the position in phase space, and it is crucial to obtain an estimate of the error that derives from the fini…
▽ More
Machine learning models have emerged as a very effective strategy to sidestep time-consuming electronic-structure calculations, enabling accurate simulations of greater size, time scale and complexity. Given the interpolative nature of these models, the reliability of predictions depends on the position in phase space, and it is crucial to obtain an estimate of the error that derives from the finite number of reference structures included during the training of the model. When using a machine-learning potential to sample a finite-temperature ensemble, the uncertainty on individual configurations translates into an error on thermodynamic averages, and provides an indication for the loss of accuracy when the simulation enters a previously unexplored region. Here we discuss how uncertainty quantification can be used, together with a baseline energy model, or a more robust although less accurate interatomic potential, to obtain more resilient simulations and to support active-learning strategies. Furthermore, we introduce an on-the-fly reweighing scheme that makes it possible to estimate the uncertainty in the thermodynamic averages extracted from long trajectories. We present examples covering different types of structural and thermodynamic properties, and systems as diverse as water and liquid gallium.
△ Less
Submitted 14 January, 2021; v1 submitted 9 November, 2020;
originally announced November 2020.
-
Inexpensive modelling of quantum dynamics using path integral generalized Langevin equation thermostats
Authors:
Venkat Kapil,
David M. Wilkins,
Jinggang Lan,
Michele Ceriotti
Abstract:
The properties of molecules and materials containing light nuclei are affected by their quantum mechanical nature. Modelling these quantum nuclear effects accurately requires computationally demanding path integral techniques. Considerable success has been achieved in reducing the cost of such simulations by using generalized Langevin dynamics to induce frequency-dependent fluctuations. Path integ…
▽ More
The properties of molecules and materials containing light nuclei are affected by their quantum mechanical nature. Modelling these quantum nuclear effects accurately requires computationally demanding path integral techniques. Considerable success has been achieved in reducing the cost of such simulations by using generalized Langevin dynamics to induce frequency-dependent fluctuations. Path integral generalized Langevin equation methods, however, have this far been limited to the study of static, thermodynamic properties due to the large perturbation to the system's dynamics induced by the aggressive thermostatting. Here we introduce a post-processing scheme, based on analytical estimates of the dynamical perturbation induced by the generalized Langevin dynamics, that makes it possible to recover meaningful time correlation properties from a thermostatted trajectory. We show that this approach yields spectroscopic observables for model and realistic systems which have an accuracy comparable to much more demanding approximate quantum dynamics techniques based on full path integral simulations.
△ Less
Submitted 6 December, 2019;
originally announced December 2019.
-
Modeling the structural and thermal properties of loaded metal-organic frameworks. An interplay of quantum and anharmonic fluctuations
Authors:
Venkat Kapil,
Jelle Wieme,
Steven Vandenbrande,
Aran Lamaire,
Veronique Van Speybroeck,
Michele Ceriotti
Abstract:
Metal-organic frameworks show both fundamental interest and great promise for applications in adsorption-based technologies, such as the separation and storage of gases. The flexibility and complexity of the molecular scaffold poses a considerable challenge to atomistic modeling, especially when also considering the presence of guest molecules. We investigate the role played by quantum and anharmo…
▽ More
Metal-organic frameworks show both fundamental interest and great promise for applications in adsorption-based technologies, such as the separation and storage of gases. The flexibility and complexity of the molecular scaffold poses a considerable challenge to atomistic modeling, especially when also considering the presence of guest molecules. We investigate the role played by quantum and anharmonic fluctuations in the archetypical case of MOF-5, comparing the material at various levels of methane loading. Accurate path integral simulations of such effects are made affordable by the introduction of an accelerated simulation scheme and the use of an optimized force field based on first-principles reference calculations. We find that the level of statistical treatment that is required for predictive modeling depends significantly on the property of interest. The thermal properties of the lattice are generally well described by a quantum harmonic treatment, with the adsorbate behaving in a classical but strongly anharmonic manner. The heat capacity of the loaded framework - which plays an important role in the characterization of the framework and in determining its stability to thermal fluctuations during adsorption/desorption cycles - requires, however, a full quantum and anharmonic treatment, either by path integral methods or by a simple but approximate scheme. We also present molecular-level insight into the nanoscopic interactions contributing to the material's properties and suggest design principles to optimize them
△ Less
Submitted 11 January, 2019;
originally announced January 2019.
-
i-PI 2.0: A Universal Force Engine for Advanced Molecular Simulations
Authors:
Venkat Kapil,
Mariana Rossi,
Ondrej Marsalek,
Riccardo Petraglia,
Yair Litman,
Thomas Spura,
Bingqing Cheng,
Alice Cuzzocrea,
Robert H. Meißner,
David M. Wilkins,
Przemyslaw Juda,
Sébastien P. Bienvenue,
Wei Fang,
Jan Kessler,
Igor Poltavsky,
Steven Vandenbrande,
Jelle Wieme,
Clemence Corminboeuf,
Thomas D. Kühne,
David E. Manolopoulos,
Thomas E. Markland,
Jeremy O. Richardson,
Alexandre Tkatchenko,
Gareth A. Tribello,
Veronique Van Speybroeck
, et al. (1 additional authors not shown)
Abstract:
Progress in the atomic-scale modelling of matter over the past decade has been tremendous. This progress has been brought about by improvements in methods for evaluating interatomic forces that work by either solving the electronic structure problem explicitly, or by computing accurate approximations of the solution and by the development of techniques that use the Born-Oppenheimer (BO) forces to…
▽ More
Progress in the atomic-scale modelling of matter over the past decade has been tremendous. This progress has been brought about by improvements in methods for evaluating interatomic forces that work by either solving the electronic structure problem explicitly, or by computing accurate approximations of the solution and by the development of techniques that use the Born-Oppenheimer (BO) forces to move the atoms on the BO potential energy surface. As a consequence of these developments it is now possible to identify stable or metastable states, to sample configurations consistent with the appropriate thermodynamic ensemble, and to estimate the kinetics of reactions and phase transitions. All too often, however, progress is slowed down by the bottleneck associated with implementing new optimization algorithms and/or sampling techniques into the many existing electronic-structure and empirical-potential codes. To address this problem, we are thus releasing a new version of the i-PI software. This piece of software is an easily extensible framework for implementing advanced atomistic simulation techniques using interatomic potentials and forces calculated by an external driver code. While the original version of the code was developed with a focus on path integral molecular dynamics techniques, this second release of i-PI not only includes several new advanced path integral methods, but also offers other classes of algorithms. In other words, i-PI is moving towards becoming a universal force engine that is both modular and tightly coupled to the driver codes that evaluate the potential energy surface and its derivatives.
△ Less
Submitted 17 September, 2018; v1 submitted 11 August, 2018;
originally announced August 2018.
-
The Anisotropy of the Proton Momentum Distribution in Water
Authors:
Venkat Kapil,
Alice Cuzzocrea,
Michele Ceriotti
Abstract:
One of the many peculiar properties of water is the pronounced deviation of the proton momentum distribution from Maxwell-Boltzmann behaviour. This deviation from the classical limit is a manifestation of the quantum mechanical nature of protons. Its extent, which can be probed directly by Deep Inelastic Neutron Scattering (DINS) experiments, gives important insight on the potential of mean force…
▽ More
One of the many peculiar properties of water is the pronounced deviation of the proton momentum distribution from Maxwell-Boltzmann behaviour. This deviation from the classical limit is a manifestation of the quantum mechanical nature of protons. Its extent, which can be probed directly by Deep Inelastic Neutron Scattering (DINS) experiments, gives important insight on the potential of mean force felt by H atoms. The determination of the full distribution of particle momenta, however, is a real tour de force for both experiments and theory, which has led to unresolved discrepancies between the two. In this Letter we present comprehensive, fully-converged momentum distributions for water at several thermodynamic state points, focusing on the components that cannot be described in terms of a scalar contribution to the quantum kinetic energy, and providing a benchmark that can serve as a reference for future simulations and experiments. In doing so, we also introduce a number of technical developments that simplify and accelerate greatly the calculation of momentum distributions by means of atomistic simulations.
△ Less
Submitted 3 May, 2018;
originally announced May 2018.
-
Fine Tuning Classical and Quantum Molecular Dynamics using a Generalized Langevin Equation
Authors:
Mariana Rossi,
Venkat Kapil,
Michele Ceriotti
Abstract:
Generalized Langevin Equation (GLE) thermostats have been used very effectively as a tool to manipulate and optimize the sampling of thermodynamic ensembles and the associated static properties. Here we show that a similar, exquisite level of control can be achieved for the dynamical properties computed from thermostatted trajectories. By developing quantitative measures of the disturbance induced…
▽ More
Generalized Langevin Equation (GLE) thermostats have been used very effectively as a tool to manipulate and optimize the sampling of thermodynamic ensembles and the associated static properties. Here we show that a similar, exquisite level of control can be achieved for the dynamical properties computed from thermostatted trajectories. By developing quantitative measures of the disturbance induced by the GLE to the Hamiltonian dynamics of a harmonic oscillator, we show that these analytical results accurately predict the behavior of strongly anharmonic systems. We also show that it is possible to correct, to a significant extent, the effects of the GLE term onto the corresponding microcanonical dynamics, which puts on more solid grounds the use of non-equilibrium Langevin dynamics to approximate quantum nuclear effects and could help improve the prediction of dynamical quantities from techniques that use a Langevin term to stabilize dynamics. Finally we address the use of thermostats in the context of approximate path-integral-based models of quantum nuclear dynamics. We demonstrate that a custom-tailored GLE can alleviate some of the artifacts associated with these techniques, improving the quality of results for the modelling of vibrational dynamics of molecules, liquids and solids.
△ Less
Submitted 17 April, 2017;
originally announced April 2017.
-
High Order Path Integrals Made Easy
Authors:
Venkat Kapil,
Jörg Behler,
Michele Ceriotti
Abstract:
The precise description of quantum nuclear fluctuations in atomistic modelling is possible by employing path integral techniques, which involve a considerable computational overhead due to the need of simulating multiple replicas of the system. Many approaches have been suggested to reduce the required number of replicas. Among these, high-order factorizations of the Boltzmann operator are particu…
▽ More
The precise description of quantum nuclear fluctuations in atomistic modelling is possible by employing path integral techniques, which involve a considerable computational overhead due to the need of simulating multiple replicas of the system. Many approaches have been suggested to reduce the required number of replicas. Among these, high-order factorizations of the Boltzmann operator are particularly attractive for high-precision and low-temperature scenarios. Unfortunately, to date several technical challenges have prevented a widespread use of these approaches to study nuclear quantum effects in condensed-phase systems. Here we introduce an inexpensive molecular dynamics scheme that overcomes these limitations, thus making it possible to exploit the improved convergence of high-order path integrals without having to sacrifice the stability, convenience and flexibility of conventional second-order techniques. The capabilities of the method are demonstrated by simulations of liquid water and ice, as described by a neural-network potential fitted to dispersion-corrected hybrid density functional theory calculations.
△ Less
Submitted 30 December, 2016; v1 submitted 2 June, 2016;
originally announced June 2016.
-
Accurate molecular dynamics and nuclear quantum effects at low cost by multiple steps in real and imaginary time: using density functional theory to accelerate wavefunction methods
Authors:
Venkat Kapil,
Joost VandeVondele,
Michele Ceriotti
Abstract:
The development and implementation of increasingly accurate methods for electronic structure calculations mean that, for many atomistic simulation problems, treating light nuclei as classical particles is now one of the most serious approximations. Even though recent developments have significantly reduced the overhead for modelling the quantum nature of the nuclei, the cost is still prohibitive w…
▽ More
The development and implementation of increasingly accurate methods for electronic structure calculations mean that, for many atomistic simulation problems, treating light nuclei as classical particles is now one of the most serious approximations. Even though recent developments have significantly reduced the overhead for modelling the quantum nature of the nuclei, the cost is still prohibitive when combined with advanced electronic structure methods. Here we present how multiple time step integrators can be combined with ring-polymer contraction techniques (effectively, multiple time stepping in imaginary time) to reduce virtually to zero the overhead of modelling nuclear quantum effects, while describing inter-atomic forces at high levels of electronic structure theory. This is demonstrated for a combination of MP2 and semi-local DFT applied to the Zundel cation. The approach can be seamlessly combined with other methods to reduce the computational cost of path integral calculations, such as high-order factorizations of the Boltzmann operator, or generalized Langevin equation thermostats.
△ Less
Submitted 1 December, 2015;
originally announced December 2015.
-
Sampling Free Energy Surfaces as Slices by Combining Umbrella Sampling and Metadynamics
Authors:
Shalini Awasthi,
Venkat Kapil,
Nisanth N. Nair
Abstract:
Metadynamics (MTD) is a very powerful technique to sample high-dimensional free energy landscapes, and due to its self-guiding property, the method has been successful in studying complex reactions and conformational changes. MTD sampling is based on filling the free energy basins by biasing potentials and thus for cases with flat, broad and unbound free energy wells, the computational time to sam…
▽ More
Metadynamics (MTD) is a very powerful technique to sample high-dimensional free energy landscapes, and due to its self-guiding property, the method has been successful in studying complex reactions and conformational changes. MTD sampling is based on filling the free energy basins by biasing potentials and thus for cases with flat, broad and unbound free energy wells, the computational time to sample them becomes very large. To alleviate this problem, we combine the standard Umbrella Sampling (US) technique with MTD to sample orthogonal collective variables (CVs) in a simultaneous way. Within this scheme, we construct the equilibrium distribution of CVs from biased distributions obtained from independent MTD simulations with umbrella potentials. Reweighting is carried out by a procedure that combines US reweighting and Tiwary-Parrinello MTD reweighting within the Weighted Histogram Analysis Method (WHAM). The approach is ideal for a controlled sampling of a CV in a MTD simulation, making it computationally efficient in sampling flat, broad and unbound free energy surfaces. This technique also allows for a distributed sampling of a high-dimensional free energy surface, further increasing the computational efficiency in sampling. We demonstrate the application of this technique in sampling high-dimensional surface for various chemical reactions using ab initio and QM/MM hybrid molecular dynamics simulations. Further, in order to carry out MTD bias reweighting for computing forward reaction barriers in ab initio or QM/MM simulations, we propose a computationally affordable approach that does not require recrossing trajectories.
△ Less
Submitted 4 February, 2016; v1 submitted 24 July, 2015;
originally announced July 2015.