1
|
Zhang P, Feng M, Xu X. Double-Layer Distribution of Hydronium and Hydroxide Ions in the Air-Water Interface. ACS PHYSICAL CHEMISTRY AU 2024; 4:336-346. [PMID: 39069983 PMCID: PMC11274287 DOI: 10.1021/acsphyschemau.3c00076] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/28/2023] [Revised: 04/04/2024] [Accepted: 04/05/2024] [Indexed: 07/30/2024]
Abstract
The acid-base nature of the aqueous interface has long been controversial. Most macroscopic experiments suggest that the air-water interface is basic based on the detection of negative charges at the interface that indicates the enrichment of hydroxides (OH-), whereas microscopic studies mostly support the acidic air-water interface with the observation of hydronium (H3O+) accumulation in the top layer of the interface. It is crucial to clarify the interfacial preference of OH- and H3O+ ions for rationalizing the debate. In this work, we perform deep potential molecular dynamics simulations to investigate the preferential distribution of OH- and H3O+ ions at the aqueous interfaces. The neural network potential energy surface is trained based on density functional theory calculations with the SCAN functional, which can accurately describe the diffusion of these two ions both in the interface and in the bulk water. In contrast to the previously reported single ion enrichment, we show that both OH- and H3O+ surprisingly prefer to accumulate in interfaces but at different interfacial depths, rendering a double-layer ionic distribution within ∼1 nm near the Gibbs dividing surface. The H3O+ preferentially resides in the topmost layer of the interface, but the OH-, which is enriched in the deeper interfacial layer, has a higher equilibrium concentration due to the more negative free energy of interfacial stabilization [-0.90 (OH-) vs -0.56 (H3O+) kcal/mol]. The present finding of the ionic double-layer distribution may qualitatively offer a self-consistent explanation for the long-term controversy about the acid-base nature of the air-water interface.
Collapse
Affiliation(s)
- Pengchao Zhang
- Center
for Combustion Energy, Department of Energy and Power Engineering,
and Key Laboratory for Thermal Science and Power Engineering of Ministry
of Education, Tsinghua University, Beijing 100084, China
| | - Muye Feng
- School
of Mechanical and Power Engineering, Nanjing
Tech University, Nanjing 211816, China
| | - Xuefei Xu
- Center
for Combustion Energy, Department of Energy and Power Engineering,
and Key Laboratory for Thermal Science and Power Engineering of Ministry
of Education, Tsinghua University, Beijing 100084, China
| |
Collapse
|
2
|
Zhai Y, Rashmi R, Palos E, Paesani F. Many-body interactions and deep neural network potentials for water. J Chem Phys 2024; 160:144501. [PMID: 38587225 DOI: 10.1063/5.0203682] [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/2024] [Accepted: 03/23/2024] [Indexed: 04/09/2024] Open
Abstract
We present a detailed assessment of deep neural network potentials developed within the Deep Potential Molecular Dynamics (DeePMD) framework and trained on the MB-pol data-driven many-body potential energy function. Specific focus is directed at the ability of DeePMD-based potentials to correctly reproduce the accuracy of MB-pol across various water systems. Analyses of bulk and interfacial properties as well as many-body interactions characteristic of water elucidate inherent limitations in the transferability and predictive accuracy of DeePMD-based potentials. These limitations can be traced back to an incomplete implementation of the "nearsightedness of electronic matter" principle, which may be common throughout machine learning potentials that do not include a proper representation of self-consistently determined long-range electric fields. These findings provide further support for the "short-blanket dilemma" faced by DeePMD-based potentials, highlighting the challenges in achieving a balance between computational efficiency and a rigorous, physics-based representation of the properties of water. Finally, we believe that our study contributes to the ongoing discourse on the development and application of machine learning models in simulating water systems, offering insights that could guide future improvements in the field.
Collapse
Affiliation(s)
- Yaoguang Zhai
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
- Department of Computer Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
| | - Richa Rashmi
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Etienne Palos
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
- Halicioğlu Data Science Institute, University of California San Diego, La Jolla, California 92093, USA
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
3
|
Shen H, Shen X, Wu Z. Simulating the isotropic Raman spectra of O-H stretching mode in liquid H 2O based on a machine learning potential: the influence of vibrational couplings. Phys Chem Chem Phys 2023; 25:28180-28188. [PMID: 37819214 DOI: 10.1039/d3cp03035k] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/13/2023]
Abstract
In this study, we trained a deep potential (DP) for H2O, an accurate machine learning (ML) potential. We performed molecular dynamics (MD) simulations of liquid water using the DP model (or DeePMD simulations). Our results showed that the DP model exhibits DFT-level accuracy, and the DeePMD simulation is a promising approach for modeling the structural properties of liquid water. Based on the DeePMD simulation trajectories, we calculated the isotropic Raman spectra of the O-H stretching mode using the surface-specific velocity-velocity correlation function (ssVVCF), showing that the DeePMD/ssVVCF approach can correctly capture the bimodal characteristics of the experimental Raman spectra, with one peak located near 3400 cm-1 and the other near 3250 cm-1. The success of the DeePMD/ssVVCF approach should be credited to (1) the DFT-level accuracy of the DP model for H2O, (2) the ssVVCF formulation considering the coupling between vibrational modes, and (3) non-Condon effects. Furthermore, the DeePMD simulations revealed that the anharmonic interactions between the coupled water molecules in the first and second hydration shells should play an essential role in the strong mixing of the H-O-H bending mode and the O-H stretching mode, leading to the delocalization of the O-H stretching band. In particular, increasing the strength of hydrogen bonds would enhance the bend-stretch coupling, leading to the red-shifting of the O-H vibrational spectra and the increase in the intensity of the shoulder around 3250 cm-1.
Collapse
Affiliation(s)
- Hujun Shen
- Guizhou Provincial Key Laboratory of Computational Nano-Material Science, Guizhou Education University, Guiyang 550018, China.
| | - Xu Shen
- National Center of Technology Innovation for Intelligent Design and Numerical Control, Huazhong University of Science and Technology, Wuhan, 430074, China
| | - Zhenhua Wu
- Department of Big Data and Artificial Intelligence, Guizhou Vocational Technology College of Electronics & Information, Kaili, 556000, China
| |
Collapse
|
4
|
Zeng J, Zhang D, Lu D, Mo P, Li Z, Chen Y, Rynik M, Huang L, Li Z, Shi S, Wang Y, Ye H, Tuo P, Yang J, Ding Y, Li Y, Tisi D, Zeng Q, Bao H, Xia Y, Huang J, Muraoka K, Wang Y, Chang J, Yuan F, Bore SL, Cai C, Lin Y, Wang B, Xu J, Zhu JX, Luo C, Zhang Y, Goodall REA, Liang W, Singh AK, Yao S, Zhang J, Wentzcovitch R, Han J, Liu J, Jia W, York DM, E W, Car R, Zhang L, Wang H. DeePMD-kit v2: A software package for deep potential models. J Chem Phys 2023; 159:054801. [PMID: 37526163 PMCID: PMC10445636 DOI: 10.1063/5.0155600] [Citation(s) in RCA: 41] [Impact Index Per Article: 41.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2023] [Accepted: 07/03/2023] [Indexed: 08/02/2023] Open
Abstract
DeePMD-kit is a powerful open-source software package that facilitates molecular dynamics simulations using machine learning potentials known as Deep Potential (DP) models. This package, which was released in 2017, has been widely used in the fields of physics, chemistry, biology, and material science for studying atomistic systems. The current version of DeePMD-kit offers numerous advanced features, such as DeepPot-SE, attention-based and hybrid descriptors, the ability to fit tensile properties, type embedding, model deviation, DP-range correction, DP long range, graphics processing unit support for customized operators, model compression, non-von Neumann molecular dynamics, and improved usability, including documentation, compiled binary packages, graphical user interfaces, and application programming interfaces. This article presents an overview of the current major version of the DeePMD-kit package, highlighting its features and technical details. Additionally, this article presents a comprehensive procedure for conducting molecular dynamics as a representative application, benchmarks the accuracy and efficiency of different models, and discusses ongoing developments.
Collapse
Affiliation(s)
- Jinzhe Zeng
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | | | - Denghui Lu
- HEDPS, CAPT, College of Engineering, Peking University, Beijing 100871, People’s Republic of China
| | - Pinghui Mo
- College of Electrical and Information Engineering, Hunan University, Changsha, People’s Republic of China
| | - Zeyu Li
- Yuanpei College, Peking University, Beijing 100871, People’s Republic of China
| | - Yixiao Chen
- Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08540, USA
| | - Marián Rynik
- Department of Experimental Physics, Comenius University, Mlynská Dolina F2, 842 48 Bratislava, Slovakia
| | - Li’ang Huang
- Center for Quantum Information, Institute for Interdisciplinary Information Sciences, Tsinghua University, Beijing 100084, People’s Republic of China
| | | | - Shaochen Shi
- ByteDance Research, Zhonghang Plaza, No. 43, North 3rd Ring West Road, Haidian District, Beijing, People’s Republic of China
| | | | - Haotian Ye
- Yuanpei College, Peking University, Beijing 100871, People’s Republic of China
| | - Ping Tuo
- AI for Science Institute, Beijing 100080, People’s Republic of China
| | - Jiabin Yang
- Baidu, Inc., Beijing, People’s Republic of China
| | | | - Yifan Li
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | | | - Qiyu Zeng
- Department of Physics, National University of Defense Technology, Changsha, Hunan 410073, People’s Republic of China
| | | | - Yu Xia
- ByteDance Research, Zhonghang Plaza, No. 43, North 3rd Ring West Road, Haidian District, Beijing, People’s Republic of China
| | | | - Koki Muraoka
- Department of Chemical System Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan
| | - Yibo Wang
- DP Technology, Beijing 100080, People’s Republic of China
| | | | - Fengbo Yuan
- DP Technology, Beijing 100080, People’s Republic of China
| | - Sigbjørn Løland Bore
- Hylleraas Centre for Quantum Molecular Sciences and Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, 0315 Oslo, Norway
| | | | - Yinnian Lin
- Wangxuan Institute of Computer Technology, Peking University, Beijing 100871, People’s Republic of China
| | - Bo Wang
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Key Laboratory of Green Chemistry and Chemical Process, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, People’s Republic of China
| | - Jiayan Xu
- School of Chemistry and Chemical Engineering, Queen’s University Belfast, Belfast BT9 5AG, United Kingdom
| | - Jia-Xin Zhu
- State Key Laboratory of Physical Chemistry of Solid Surfaces, iChEM, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, People’s Republic of China
| | - Chenxing Luo
- Department of Applied Physics and Applied Mathematics, Columbia University, New York, New York 10027, USA
| | - Yuzhi Zhang
- DP Technology, Beijing 100080, People’s Republic of China
| | | | - Wenshuo Liang
- DP Technology, Beijing 100080, People’s Republic of China
| | - Anurag Kumar Singh
- Department of Data Science, Indian Institute of Technology, Palakkad, Kerala, India
| | - Sikai Yao
- DP Technology, Beijing 100080, People’s Republic of China
| | - Jingchao Zhang
- NVIDIA AI Technology Center (NVAITC), Santa Clara, California 95051, USA
| | | | - Jiequn Han
- Center for Computational Mathematics, Flatiron Institute, New York, New York 10010, USA
| | - Jie Liu
- College of Electrical and Information Engineering, Hunan University, Changsha, People’s Republic of China
| | | | - Darrin M. York
- Laboratory for Biomolecular Simulation Research, Institute for Quantitative Biomedicine and Department of Chemistry and Chemical Biology, Rutgers University, Piscataway, New Jersey 08854, USA
| | | | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | | | - Han Wang
- Author to whom correspondence should be addressed:
| |
Collapse
|
5
|
Ko HY, Calegari Andrade MF, Sparrow ZM, Zhang JA, DiStasio RA. High-Throughput Condensed-Phase Hybrid Density Functional Theory for Large-Scale Finite-Gap Systems: The SeA Approach. J Chem Theory Comput 2023. [PMID: 37385014 DOI: 10.1021/acs.jctc.2c00827] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/01/2023]
Abstract
High-throughput electronic structure calculations (often performed using density functional theory (DFT)) play a central role in screening existing and novel materials, sampling potential energy surfaces, and generating data for machine learning applications. By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal DFT and furnish a more accurate description of the underlying electronic structure, albeit at a computational cost that often prohibits such high-throughput applications. To address this challenge, we have constructed a robust, accurate, and computationally efficient framework for high-throughput condensed-phase hybrid DFT and implemented this approach in the PWSCF module of Quantum ESPRESSO (QE). The resulting SeA approach (SeA = SCDM + exx + ACE) combines and seamlessly integrates: (i) the selected columns of the density matrix method (SCDM, a robust noniterative orbital localization scheme that sidesteps system-dependent optimization protocols), (ii) a recently extended version of exx (a black-box linear-scaling EXX algorithm that exploits sparsity between localized orbitals in real space when evaluating the action of the standard/full-rank V^xx operator), and (iii) adaptively compressed exchange (ACE, a low-rank V^xx approximation). In doing so, SeA harnesses three levels of computational savings: pair selection and domain truncation from SCDM + exx (which only considers spatially overlapping orbitals on orbital-pair-specific and system-size-independent domains) and low-rank V^xx approximation from ACE (which reduces the number of calls to SCDM + exx during the self-consistent field (SCF) procedure). Across a diverse set of 200 nonequilibrium (H2O)64 configurations (with densities spanning 0.4-1.7 g/cm3), SeA provides a 1-2 order-of-magnitude speedup in the overall time-to-solution, i.e., ≈8-26× compared to the convolution-based PWSCF(ACE) implementation in QE and ≈78-247× compared to the conventional PWSCF(Full) approach, and yields energies, ionic forces, and other properties with high fidelity. As a proof-of-principle high-throughput application, we trained a deep neural network (DNN) potential for ambient liquid water at the hybrid DFT level using SeA via an actively learned data set with ≈8,700 (H2O)64 configurations. Using an out-of-sample set of (H2O)512 configurations (at nonambient conditions), we confirmed the accuracy of this SeA-trained potential and showcased the capabilities of SeA by computing the ground-truth ionic forces in this challenging system containing >1,500 atoms.
Collapse
Affiliation(s)
- Hsin-Yu Ko
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Marcos F Calegari Andrade
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
- Quantum Simulations Group, Materials Science Division, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Zachary M Sparrow
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Ju-An Zhang
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
6
|
Zhu X, Xu J, Zhang Y, Li B, Tian Y, Wu Y, Liu Z, Ma C, Tan S, Wang B. Revealing Intramolecular Isotope Effects with Chemical-Bond Precision. J Am Chem Soc 2023. [PMID: 37338304 DOI: 10.1021/jacs.3c02728] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/21/2023]
Abstract
Isotope substitution of a molecule not only changes its vibrational frequencies but also changes its vibrational distributions in real-space. Quantitatively measuring the isotope effects inside a polyatomic molecule requires both energy and spatial resolutions at the single-bond level, which has been a long-lasting challenge in macroscopic techniques. By achieving ångström resolution in tip-enhanced Raman spectroscopy (TERS), we record the corresponding local vibrational modes of pentacene and its fully deuterated form, enabling us to identify and measure the isotope effect of each vibrational mode. The measured frequency ratio νH/νD varies from 1.02 to 1.33 in different vibrational modes, indicating different isotopic contributions of H/D atoms, which can be distinguished from TERS maps in real-space and well described by the potential energy distribution simulations. Our study demonstrates that TERS can serve as a non-destructive and highly sensitive methodology for isotope detection and recognition with chemical-bond precision.
Collapse
Affiliation(s)
- Xiang Zhu
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information & Quantum Physics, New Cornerstone Science Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Jiayu Xu
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information & Quantum Physics, New Cornerstone Science Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yao Zhang
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information & Quantum Physics, New Cornerstone Science Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230088, China
| | - Bin Li
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information & Quantum Physics, New Cornerstone Science Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230088, China
| | - Yishu Tian
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information & Quantum Physics, New Cornerstone Science Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Yingying Wu
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information & Quantum Physics, New Cornerstone Science Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Zhiwei Liu
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information & Quantum Physics, New Cornerstone Science Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
| | - Chuanxu Ma
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information & Quantum Physics, New Cornerstone Science Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230088, China
| | - Shijing Tan
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information & Quantum Physics, New Cornerstone Science Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230088, China
| | - Bing Wang
- Hefei National Research Center for Physical Sciences at the Microscale and Synergetic Innovation Center of Quantum Information & Quantum Physics, New Cornerstone Science Laboratory, University of Science and Technology of China, Hefei, Anhui 230026, China
- Hefei National Laboratory, University of Science and Technology of China, Hefei, Anhui 230088, China
| |
Collapse
|
7
|
Zhu X, Huang C, Li N, Ma X, Li Z, Fan J. Distinct roles of graphene and graphene oxide nanosheets in regulating phospholipid flip-flop. J Colloid Interface Sci 2023; 637:112-122. [PMID: 36689797 DOI: 10.1016/j.jcis.2023.01.080] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2022] [Revised: 12/30/2022] [Accepted: 01/15/2023] [Indexed: 01/19/2023]
Abstract
Two-dimensional (2D) nanomaterials, such as graphene nanosheets (GNs) and graphene oxide nanosheets (GOs), could adhere onto or insert into a biological membrane, leading to a change in membrane properties and biological activities. Consequently, GN and GO become potential candidates for mediating interleaflet phospholipid transfer. In this work, molecular dynamics (MD) simulations were employed to investigate the effects of GN and GO on lipid flip-flop behavior and the underlying molecular mechanisms. Of great interest is that GN and GO work in opposite directions. The inserted GN can induce the formation of an ordered nanodomain, which dramatically elevates the free energy barrier of flipping phospholipids from one leaflet to the other, thus leading to a decreased lipid flip-flop rate. In contrast, the embedded GO can catalyze the transport of phospholipids between membrane leaflets by facilitating the formation of water pores. These results suggest that GN may work as an inhibitor of the interleaflet lipid translocation, while GO may play the role of scramblases. These findings are expected to expand promising biomedical applications of 2D nanomaterials.
Collapse
Affiliation(s)
- Xiaohong Zhu
- Department of Materials Science and Engineering, City University of Hong Kong, Hong Kong, China
| | - Changxiong Huang
- Department of Materials Science and Engineering, City University of Hong Kong, Hong Kong, China
| | - Na Li
- Department of Materials Science and Engineering, City University of Hong Kong, Hong Kong, China
| | - Xinyao Ma
- Department of Materials Science and Engineering, City University of Hong Kong, Hong Kong, China
| | - Zhen Li
- School of Materials Science and Engineering, China University of Petroleum (East China), Qingdao 266580, China.
| | - Jun Fan
- Department of Materials Science and Engineering, City University of Hong Kong, Hong Kong, China; Center for Advanced Nuclear Safety and Sustainable Development, City University of Hong Kong, Hong Kong, China.
| |
Collapse
|
8
|
Plé T, Mauger N, Adjoua O, Inizan TJ, Lagardère L, Huppert S, Piquemal JP. Routine Molecular Dynamics Simulations Including Nuclear Quantum Effects: From Force Fields to Machine Learning Potentials. J Chem Theory Comput 2023; 19:1432-1445. [PMID: 36856658 DOI: 10.1021/acs.jctc.2c01233] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/02/2023]
Abstract
We report the implementation of a multi-CPU and multi-GPU massively parallel platform dedicated to the explicit inclusion of nuclear quantum effects (NQEs) in the Tinker-HP molecular dynamics (MD) package. The platform, denoted Quantum-HP, exploits two simulation strategies: the Ring-Polymer Molecular Dynamics (RPMD) that provides exact structural properties at the cost of a MD simulation in an extended space of multiple replicas and the adaptive Quantum Thermal Bath (adQTB) that imposes the quantum distribution of energy on a classical system via a generalized Langevin thermostat and provides computationally affordable and accurate (though approximate) NQEs. We discuss some implementation details, efficient numerical schemes, and parallelization strategies and quickly review the GPU acceleration of our code. Our implementation allows an efficient inclusion of NQEs in MD simulations for very large systems, as demonstrated by scaling tests on water boxes with more than 200,000 atoms (simulated using the AMOEBA polarizable force field). We test the compatibility of the approach with Tinker-HP's recently introduced Deep-HP machine learning potentials module by computing water properties using the DeePMD potential with adQTB thermostatting. Finally, we show that the platform is also compatible with the alchemical free energy estimation capabilities of Tinker-HP and fast enough to perform simulations. Therefore, we study how NQEs affect the hydration free energy of small molecules solvated with the recently developed Q-AMOEBA water force field. Overall, the Quantum-HP platform allows users to perform routine quantum MD simulations of large condensed-phase systems and will help to shed new light on the quantum nature of important interactions in biological matter.
Collapse
Affiliation(s)
- Thomas Plé
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Nastasia Mauger
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Olivier Adjoua
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | | | - Louis Lagardère
- Sorbonne Université, LCT, UMR 7616 CNRS, F-75005 Paris, France
| | - Simon Huppert
- Institut des Nanosciences de Paris (INSP), CNRS UMR 7588, and Sorbonne Université, 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
| |
Collapse
|
9
|
Chen BW. Equilibrium and kinetic isotope effects in heterogeneous catalysis: A density functional theory perspective. CATAL COMMUN 2023. [DOI: 10.1016/j.catcom.2023.106654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/29/2023] Open
|
10
|
Sours T, Kulkarni AR. Predicting Structural Properties of Pure Silica Zeolites Using Deep Neural Network Potentials. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2023; 127:1455-1463. [PMID: 36733763 PMCID: PMC9885523 DOI: 10.1021/acs.jpcc.2c08429] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/07/2022] [Revised: 12/20/2022] [Indexed: 06/18/2023]
Abstract
Machine learning potentials (MLPs) capable of accurately describing complex ab initio potential energy surfaces (PESs) have revolutionized the field of multiscale atomistic modeling. In this work, using an extensive density functional theory (DFT) data set (denoted as Si-ZEO22) consisting of 219 unique zeolite topologies (350,000 unique DFT calculations) found in the International Zeolite Association (IZA) database, we have trained a DeePMD-kit MLP to model the dynamics of silica frameworks. The performance of our model is evaluated by calculating various properties that probe the accuracy of the energy and force predictions. This MLP demonstrates impressive agreement with DFT for predicting zeolite structural properties, energy-volume trends, and phonon density of states. Furthermore, our model achieves reasonable predictions for stress-strain relationships without including DFT stress data during training. These results highlight the ability of MLPs to capture the flexibility of zeolite frameworks and motivate further MLP development for nanoporous materials with near-ab initio accuracy.
Collapse
|
11
|
Wisesa P, Andolina CM, Saidi WA. Development and Validation of Versatile Deep Atomistic Potentials for Metal Oxides. J Phys Chem Lett 2023; 14:468-475. [PMID: 36623167 DOI: 10.1021/acs.jpclett.2c03445] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
Machine learning interatomic potentials powered by neural networks have been shown to readily model a gradient of compositions in metallic systems. However, their application to date on ionic systems tends to focus on specific compositions and oxidation states owing to their more heterogeneous chemical nature. Herein we show that a deep neural network potential (DNP) can model various properties of metal oxides with different oxidation states without additional charge information. We created and validated DNPs for AgxOy, CuxOy MgxOy, PtxOy, and ZnxOy, whereby each system was trained without any limitations on oxidation states. We illustrate how the database can be augmented to enhance the DNP transferability for a new polymorph, surface energies, and thermal expansion. In addition, we show that these potentials can correctly interpolate significant pressure and temperature ranges, exhibit stability over long molecular dynamics simulation time scales, and replicate nonharmonic thermal expansion, consistent with experimental results.
Collapse
Affiliation(s)
- Pandu Wisesa
- Department of Mechanical Engineering and Materials Science, University of Pittsburgh, Pittsburgh, Pennsylvania 15216, United States
| | - Christopher M Andolina
- Department of Mechanical Engineering and Materials Science, University of Pittsburgh, Pittsburgh, Pennsylvania 15216, United States
| | - Wissam A Saidi
- Department of Mechanical Engineering and Materials Science, University of Pittsburgh, Pittsburgh, Pennsylvania 15216, United States
| |
Collapse
|
12
|
Liu J, Liu R, Cao Y, Chen M. Solvation structures of calcium and magnesium ions in water with the presence of hydroxide: a study by deep potential molecular dynamics. Phys Chem Chem Phys 2023; 25:983-993. [PMID: 36519362 DOI: 10.1039/d2cp04105g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
The solvation structures of calcium (Ca2+) and magnesium (Mg2+) ions with the presence of hydroxide (OH-) ion in water are essential for understanding their roles in biological and chemical processes but have not been fully explored. Ab initio molecular dynamics (AIMD) is an important tool to address this issue, but two challenges exist. First, an accurate description of OH- from AIMD needs an appropriate exchange-correlation functional. Second, a long trajectory is needed to reach an equilibrium state for the Ca2+-OH- and Mg2+-OH- ion pairs in aqueous solutions. Herein, we adopt a deep potential molecular dynamics (DPMD) method to simulate 1 ns trajectories for the Ca2+-OH- and Mg2+-OH- ion pairs in water; the DPMD method provides efficient machine-learning-based models that have the accuracy of the SCAN exchange-correlation functional within the framework of density functional theory. The solvation structures of the cations and the OH- in terms of three different species have been systematically investigated. On the one hand, we find that OH- have more significant effects on the solvation structure of Ca2+ than that of Mg2+. We observe that the OH- substantially affects the orientation angles of water molecules surrounding the cation. Through the time correlation functions, we conclude that the water molecules in the first solvation shell of Ca2+ change their preferred orientation faster than those of Mg2+. On the other hand, with the presence of the cation in the first solvation shell of OH-, we find that the hydrogen bonds of OH- are severely altered, and the adjacent water molecules of OH- are squeezed. The two cations have substantially different effects on the solvation structure of OH-. Our work provides new insight into the solvation structures of Ca2+ and Mg2+ in water with the presence of OH-.
Collapse
Affiliation(s)
- Jianchuan Liu
- HEDPS, CAPT, College of Engineering and School of Physics, Peking University, Beijing, 100871, China.
| | - Renxi Liu
- HEDPS, CAPT, College of Engineering and School of Physics, Peking University, Beijing, 100871, China. .,Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, 100871, China
| | - Yu Cao
- HEDPS, CAPT, College of Engineering and School of Physics, Peking University, Beijing, 100871, China.
| | - Mohan Chen
- HEDPS, CAPT, College of Engineering and School of Physics, Peking University, Beijing, 100871, China. .,Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, 100871, China
| |
Collapse
|
13
|
Li W, Ou Q, Chen Y, Cao Y, Liu R, Zhang C, Zheng D, Cai C, Wu X, Wang H, Chen M, Zhang L. DeePKS + ABACUS as a Bridge between Expensive Quantum Mechanical Models and Machine Learning Potentials. J Phys Chem A 2022; 126:9154-9164. [DOI: 10.1021/acs.jpca.2c05000] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Affiliation(s)
- Wenfei Li
- AI for Science Institute, Beijing100080, P. R. China
| | - Qi Ou
- AI for Science Institute, Beijing100080, P. R. China
| | - Yixiao Chen
- Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey08544, United States
| | - Yu Cao
- HEDPS, CAPT, College of Engineering and School of Physics, Peking University, Beijing100871, P. R. China
| | - Renxi Liu
- HEDPS, CAPT, College of Engineering and School of Physics, Peking University, Beijing100871, P. R. China
| | - Chunyi Zhang
- Department of Physics, Temple University, Philadelphia, Pennsylvania19122, United States
| | - Daye Zheng
- AI for Science Institute, Beijing100080, P. R. China
| | - Chun Cai
- AI for Science Institute, Beijing100080, P. R. China
- DP Technology, Beijing100080, P. R. China
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania19122, United States
| | - Han Wang
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Huayuan Road 6, Beijing100088, P. R. China
| | - Mohan Chen
- HEDPS, CAPT, College of Engineering and School of Physics, Peking University, Beijing100871, P. R. China
| | - Linfeng Zhang
- AI for Science Institute, Beijing100080, P. R. China
- DP Technology, Beijing100080, P. R. China
| |
Collapse
|
14
|
Lu D, Jiang W, Chen Y, Zhang L, Jia W, Wang H, Chen M. DP Compress: A Model Compression Scheme for Generating Efficient Deep Potential Models. J Chem Theory Comput 2022; 18:5559-5567. [PMID: 35926122 DOI: 10.1021/acs.jctc.2c00102] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Machine-learning-based interatomic potential energy surface (PES) models are revolutionizing the field of molecular modeling. However, although much faster than electronic structure schemes, these models suffer from costly computations via deep neural networks to predict the energy and atomic forces, resulting in lower running efficiency as compared to the typical empirical force fields. Herein, we report a model compression scheme for boosting the performance of the Deep Potential (DP) model, a deep learning-based PES model. This scheme, we call DP Compress, is an efficient postprocessing step after the training of DP models (DP Train). DP Compress combines several DP-specific compression techniques, which typically speed up DP-based molecular dynamics simulations by an order of magnitude faster and consume an order of magnitude less memory. We demonstrate that DP Compress is sufficiently accurate by testing a variety of physical properties of Cu, H2O, and Al-Cu-Mg systems. DP Compress applies to both CPU and GPU machines and is publicly available online.
Collapse
Affiliation(s)
- Denghui Lu
- HEDPS, CAPT, College of Engineering, Peking University, Beijing 100871, P. R. China
| | - Wanrun Jiang
- Songshan Lake Materials Laboratory, Dongguan, Guangdong 523808, P. R. China.,Institute of Physics, Chinese Academy of Sciences, Beijing 100190, P. R. China
| | - Yixiao Chen
- Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, United States
| | - Linfeng Zhang
- Beijing Institute of Big Data Research, Beijing 100871, P. R. China
| | - Weile Jia
- Institute of Computing Technology, Chinese Academy of Sciences, Beijing 100190, P. R. China.,University of Chinese Academy of Sciences, Beijing 100049, P. R. China
| | - Han Wang
- HEDPS, CAPT, College of Engineering, Peking University, Beijing 100871, P. R. China.,Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Fenghao East Road 2, Beijing 100094, P. R. China
| | - Mohan Chen
- HEDPS, CAPT, College of Engineering, Peking University, Beijing 100871, P. R. China
| |
Collapse
|
15
|
Vinson J. Advances in the OCEAN-3 spectroscopy package. Phys Chem Chem Phys 2022; 24:12787-12803. [PMID: 35608324 PMCID: PMC9844114 DOI: 10.1039/d2cp01030e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
The OCEAN code for calculating valence- and core-level spectra using the Bethe-Salpeter equation is briefly reviewed. OCEAN is capable of calculating optical absorption, near-edge X-ray absorption or non-resonant scattering, and resonant inelastic X-ray scattering, requiring only the structure of the material as input. Improved default behavior and reduced input requirements are detailed as well as new capabilities, such as incorporation of final-state-dependent broadening, finite-temperature dependence, and flexibility in the density-functional theory exchange-correlation potentials. OCEAN is built on top of a plane-wave, pseudopotential, density-functional theory foundation, and calculations are shown for systems ranging in size up to 7 nm3.
Collapse
Affiliation(s)
- John Vinson
- Material Measurement Laboratory, National Institute of Standards and Technology, Gaithersburg, MD 20899
| |
Collapse
|
16
|
Thomsen B, Shiga M. Structures of liquid and aqueous water isotopologues at ambient temperature from ab initio path integral simulations. Phys Chem Chem Phys 2022; 24:10851-10859. [PMID: 35504275 DOI: 10.1039/d2cp00499b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The heavy hydrogen isotopes D and T are found in trace amounts in water, and when their concentration increases they can play an intricate role in modulating the physical properties of the liquid. We present an analysis of the microscopic structures of ambient light water (H2O(l)), heavy water (D2O(l)), T2O(l), HDO(aq) and HTO(aq) studied by ab initio path integral molecular dynamics (PIMD). Unlike previous ab initio PIMD investigations of H2O(l) and D2O(l) [Chen et al., Phys. Rev. Lett., 2003, 91, 215503] [Machida et al., J. Chem. Phys., 2017, 148, 102324] we find that D2O(l) is more structured than H2O(l), as is predicted by the experiment. The agreement between the experiment and our simulation for H2O(l) and D2O(l) allows us to accurately predict the intra- and intermolecular structures of T2O(l) HDO(aq) and HTO(aq). T2O(l) is found to have a similar intermolecular structure to that of D2O(l), while the intramolecular structure is more compact, giving rise to a smaller dipole moment than those of H2O(l) and D2O(l). For the mixed isotope species, HDO(aq) and HTO(aq), we find smaller dipole moments and fewer hydrogen bonds when compared with the pure species H2O and D2O. We can attribute this effect to the relative compactness of the mixed isotope species, which results in a lower dipole moment than that of the pure species.
Collapse
Affiliation(s)
- Bo Thomsen
- CCSE, Japan Atomic Energy Agency, 178-4-4, Wakashiba, Kashiwa, Chiba, 277-0871, Japan.
| | - Motoyuki Shiga
- CCSE, Japan Atomic Energy Agency, 178-4-4, Wakashiba, Kashiwa, Chiba, 277-0871, Japan.
| |
Collapse
|
17
|
Batzner S, Musaelian A, Sun L, Geiger M, Mailoa JP, Kornbluth M, Molinari N, Smidt TE, Kozinsky B. E(3)-equivariant graph neural networks for data-efficient and accurate interatomic potentials. Nat Commun 2022; 13:2453. [PMID: 35508450 PMCID: PMC9068614 DOI: 10.1038/s41467-022-29939-5] [Citation(s) in RCA: 237] [Impact Index Per Article: 118.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2021] [Accepted: 04/07/2022] [Indexed: 11/08/2022] Open
Abstract
This work presents Neural Equivariant Interatomic Potentials (NequIP), an E(3)-equivariant neural network approach for learning interatomic potentials from ab-initio calculations for molecular dynamics simulations. While most contemporary symmetry-aware models use invariant convolutions and only act on scalars, NequIP employs E(3)-equivariant convolutions for interactions of geometric tensors, resulting in a more information-rich and faithful representation of atomic environments. The method achieves state-of-the-art accuracy on a challenging and diverse set of molecules and materials while exhibiting remarkable data efficiency. NequIP outperforms existing models with up to three orders of magnitude fewer training data, challenging the widely held belief that deep neural networks require massive training sets. The high data efficiency of the method allows for the construction of accurate potentials using high-order quantum chemical level of theory as reference and enables high-fidelity molecular dynamics simulations over long time scales.
Collapse
Affiliation(s)
- Simon Batzner
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA.
| | - Albert Musaelian
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
| | - Lixin Sun
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
| | - Mario Geiger
- École Polytechnique Fédérale de Lausanne, 1015, Lausanne, Switzerland
- Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
| | - Jonathan P Mailoa
- Robert Bosch Research and Technology Center, Cambridge, MA, 02139, USA
| | | | - Nicola Molinari
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA
| | - Tess E Smidt
- Computational Research Division and Center for Advanced Mathematics for Energy Research Applications, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA
- Department of Electrical Engineering and Computer Science and Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, MA, 02139, USA
| | - Boris Kozinsky
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, 02138, USA.
- Robert Bosch Research and Technology Center, Cambridge, MA, 02139, USA.
| |
Collapse
|
18
|
Yang Y, Peltier CR, Zeng R, Schimmenti R, Li Q, Huang X, Yan Z, Potsi G, Selhorst R, Lu X, Xu W, Tader M, Soudackov AV, Zhang H, Krumov M, Murray E, Xu P, Hitt J, Xu L, Ko HY, Ernst BG, Bundschu C, Luo A, Markovich D, Hu M, He C, Wang H, Fang J, DiStasio RA, Kourkoutis LF, Singer A, Noonan KJT, Xiao L, Zhuang L, Pivovar BS, Zelenay P, Herrero E, Feliu JM, Suntivich J, Giannelis EP, Hammes-Schiffer S, Arias T, Mavrikakis M, Mallouk TE, Brock JD, Muller DA, DiSalvo FJ, Coates GW, Abruña HD. Electrocatalysis in Alkaline Media and Alkaline Membrane-Based Energy Technologies. Chem Rev 2022; 122:6117-6321. [PMID: 35133808 DOI: 10.1021/acs.chemrev.1c00331] [Citation(s) in RCA: 89] [Impact Index Per Article: 44.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Hydrogen energy-based electrochemical energy conversion technologies offer the promise of enabling a transition of the global energy landscape from fossil fuels to renewable energy. Here, we present a comprehensive review of the fundamentals of electrocatalysis in alkaline media and applications in alkaline-based energy technologies, particularly alkaline fuel cells and water electrolyzers. Anion exchange (alkaline) membrane fuel cells (AEMFCs) enable the use of nonprecious electrocatalysts for the sluggish oxygen reduction reaction (ORR), relative to proton exchange membrane fuel cells (PEMFCs), which require Pt-based electrocatalysts. However, the hydrogen oxidation reaction (HOR) kinetics is significantly slower in alkaline media than in acidic media. Understanding these phenomena requires applying theoretical and experimental methods to unravel molecular-level thermodynamics and kinetics of hydrogen and oxygen electrocatalysis and, particularly, the proton-coupled electron transfer (PCET) process that takes place in a proton-deficient alkaline media. Extensive electrochemical and spectroscopic studies, on single-crystal Pt and metal oxides, have contributed to the development of activity descriptors, as well as the identification of the nature of active sites, and the rate-determining steps of the HOR and ORR. Among these, the structure and reactivity of interfacial water serve as key potential and pH-dependent kinetic factors that are helping elucidate the origins of the HOR and ORR activity differences in acids and bases. Additionally, deliberately modulating and controlling catalyst-support interactions have provided valuable insights for enhancing catalyst accessibility and durability during operation. The design and synthesis of highly conductive and durable alkaline membranes/ionomers have enabled AEMFCs to reach initial performance metrics equal to or higher than those of PEMFCs. We emphasize the importance of using membrane electrode assemblies (MEAs) to integrate the often separately pursued/optimized electrocatalyst/support and membranes/ionomer components. Operando/in situ methods, at multiscales, and ab initio simulations provide a mechanistic understanding of electron, ion, and mass transport at catalyst/ionomer/membrane interfaces and the necessary guidance to achieve fuel cell operation in air over thousands of hours. We hope that this Review will serve as a roadmap for advancing the scientific understanding of the fundamental factors governing electrochemical energy conversion in alkaline media with the ultimate goal of achieving ultralow Pt or precious-metal-free high-performance and durable alkaline fuel cells and related technologies.
Collapse
Affiliation(s)
- Yao Yang
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Cheyenne R Peltier
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Rui Zeng
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Roberto Schimmenti
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Qihao Li
- College of Chemistry and Molecular Sciences, Wuhan University, Wuhan 430072, China
| | - Xin Huang
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, United States
| | - Zhifei Yan
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Georgia Potsi
- Department of Materials Science and Engineering, Cornell University, Ithaca, New York 14853, United States
| | - Ryan Selhorst
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Xinyao Lu
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Weixuan Xu
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Mariel Tader
- Department of Physics, Cornell University, Ithaca, New York 14853, United States
| | - Alexander V Soudackov
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Hanguang Zhang
- Materials Physics and Applications Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Mihail Krumov
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Ellen Murray
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Pengtao Xu
- Department of Materials Science and Engineering, Cornell University, Ithaca, New York 14853, United States
| | - Jeremy Hitt
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Linxi Xu
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Hsin-Yu Ko
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Brian G Ernst
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Colin Bundschu
- Department of Physics, Cornell University, Ithaca, New York 14853, United States
| | - Aileen Luo
- Department of Materials Science and Engineering, Cornell University, Ithaca, New York 14853, United States
| | - Danielle Markovich
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, United States
| | - Meixue Hu
- College of Chemistry and Molecular Sciences, Wuhan University, Wuhan 430072, China
| | - Cheng He
- Chemical and Materials Science Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Hongsen Wang
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Jiye Fang
- Department of Chemistry, State University of New York at Binghamton, Binghamton, New York 13902, United States
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Lena F Kourkoutis
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, United States.,Kavli Institute at Cornell for Nanoscale Science, Cornell University, Ithaca, New York 14853, United States
| | - Andrej Singer
- Department of Materials Science and Engineering, Cornell University, Ithaca, New York 14853, United States
| | - Kevin J T Noonan
- Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| | - Li Xiao
- College of Chemistry and Molecular Sciences, Wuhan University, Wuhan 430072, China
| | - Lin Zhuang
- College of Chemistry and Molecular Sciences, Wuhan University, Wuhan 430072, China
| | - Bryan S Pivovar
- Chemical and Materials Science Center, National Renewable Energy Laboratory, Golden, Colorado 80401, United States
| | - Piotr Zelenay
- Materials Physics and Applications Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Enrique Herrero
- Instituto de Electroquímica, Universidad de Alicante, Alicante E-03080, Spain
| | - Juan M Feliu
- Instituto de Electroquímica, Universidad de Alicante, Alicante E-03080, Spain
| | - Jin Suntivich
- Department of Materials Science and Engineering, Cornell University, Ithaca, New York 14853, United States.,Kavli Institute at Cornell for Nanoscale Science, Cornell University, Ithaca, New York 14853, United States
| | - Emmanuel P Giannelis
- Department of Materials Science and Engineering, Cornell University, Ithaca, New York 14853, United States
| | | | - Tomás Arias
- Department of Physics, Cornell University, Ithaca, New York 14853, United States
| | - Manos Mavrikakis
- Department of Chemical and Biological Engineering, University of Wisconsin-Madison, Madison, Wisconsin 53706, United States
| | - Thomas E Mallouk
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Joel D Brock
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, United States
| | - David A Muller
- School of Applied and Engineering Physics, Cornell University, Ithaca, New York 14853, United States.,Kavli Institute at Cornell for Nanoscale Science, Cornell University, Ithaca, New York 14853, United States
| | - Francis J DiSalvo
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Geoffrey W Coates
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Héctor D Abruña
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.,Center for Alkaline Based Energy Solutions (CABES), Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
19
|
|
20
|
Gallo P, Bachler J, Bove LE, Böhmer R, Camisasca G, Coronas LE, Corti HR, de Almeida Ribeiro I, de Koning M, Franzese G, Fuentes-Landete V, Gainaru C, Loerting T, de Oca JMM, Poole PH, Rovere M, Sciortino F, Tonauer CM, Appignanesi GA. Advances in the study of supercooled water. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2021; 44:143. [PMID: 34825973 DOI: 10.1140/epje/s10189-021-00139-1] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Accepted: 10/17/2021] [Indexed: 06/13/2023]
Abstract
In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.
Collapse
Affiliation(s)
- Paola Gallo
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy.
| | - Johannes Bachler
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | - Livia E Bove
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale A. Moro 5, 00185, Roma, Italy
- Sorbonne Université, CNRS UMR 7590, IMPMC, 75005, Paris, France
| | - Roland Böhmer
- Fakultät Physik, Technische Universität Dortmund, 44221, Dortmund, Germany
| | - Gaia Camisasca
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy
| | - Luis E Coronas
- Secció de Física Estadística i Interdisciplinària-Departament de Física de la Matèria Condensada, Universitat de Barcelona, & Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, C. Martí i Franquès 1, 08028, Barcelona, Spain
| | - Horacio R Corti
- Departamento de Física de la Materia Condensada, Centro Atómico Constituyentes, Comisión Nacional de Energía Atómica, B1650LWP, Buenos Aires, Argentina
| | - Ingrid de Almeida Ribeiro
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859, Campinas, São Paulo, Brazil
| | - Maurice de Koning
- Instituto de Física "Gleb Wataghin", Universidade Estadual de Campinas, UNICAMP, 13083-859, Campinas, São Paulo, Brazil
- Center for Computing in Engineering & Sciences, Universidade Estadual de Campinas, UNICAMP, 13083-861, Campinas, São Paulo, Brazil
| | - Giancarlo Franzese
- Secció de Física Estadística i Interdisciplinària-Departament de Física de la Matèria Condensada, Universitat de Barcelona, & Institut de Nanociència i Nanotecnologia (IN2UB), Universitat de Barcelona, C. Martí i Franquès 1, 08028, Barcelona, Spain
| | - Violeta Fuentes-Landete
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | - Catalin Gainaru
- Fakultät Physik, Technische Universität Dortmund, 44221, Dortmund, Germany
| | - Thomas Loerting
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | | | - Peter H Poole
- Department of Physics, St. Francis Xavier University, Antigonish, NS, B2G 2W5, Canada
| | - Mauro Rovere
- Dipartimento di Matematica e Fisica, Università degli Studi Roma Tre, Via della Vasca Navale 84, 00146, Roma, Italy
| | - Francesco Sciortino
- Dipartimento di Fisica, Sapienza Università di Roma, Piazzale A. Moro 5, 00185, Roma, Italy
| | - Christina M Tonauer
- Institute of Physical Chemistry, University of Innsbruck, Innrain 52c, A-6020, Innsbruck, Austria
| | - Gustavo A Appignanesi
- INQUISUR, Departamento de Química, Universidad Nacional del Sur (UNS)-CONICET, Avenida Alem 1253, 8000, Bahía Blanca, Argentina
| |
Collapse
|
21
|
Ko HY, Santra B, DiStasio RA. Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory-Based Ab Initio Molecular Dynamics II: Extensions to the Isobaric-Isoenthalpic and Isobaric-Isothermal Ensembles. J Chem Theory Comput 2021; 17:7789-7813. [PMID: 34775753 DOI: 10.1021/acs.jctc.0c01194] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In the previous paper of this series [Ko, H.-Y. et al. J. Chem. Theory Comput. 2020, 16, 3757-3785], we presented a theoretical and algorithmic framework based on a localized representation of the occupied space that exploits the inherent sparsity in the real-space evaluation of the exact exchange (EXX) interaction in finite-gap systems. This was accompanied by a detailed description of exx, a massively parallel hybrid message-passing interface MPI/OpenMP implementation of this approach in Quantum ESPRESSO (QE) that enables linear scaling hybrid density functional theory (DFT)-based ab initio molecular dynamics (AIMD) in the microcanonical/canonical (NVE/NVT) ensembles of condensed-phase systems containing 500-1000 atoms (in fixed orthorhombic cells) with a wall time cost comparable to semi-local DFT. In this work, we extend the current capabilities of exx to enable hybrid DFT-based AIMD simulations of large-scale condensed-phase systems with general and fluctuating cells in the isobaric-isoenthalpic/isobaric-isothermal (NpH/NpT) ensembles. The theoretical extensions to this approach include an analytical derivation of the EXX contribution to the stress tensor for systems in general simulation cells with a computational complexity that scales linearly with system size. The corresponding algorithmic extensions to exx include optimized routines that (i) handle both static and fluctuating simulation cells with non-orthogonal lattice symmetries, (ii) solve Poisson's equation in general/non-orthogonal cells via an automated selection of the auxiliary grid directions in the Natan-Kronik representation of the discrete Laplacian operator, and (iii) evaluate the EXX contribution to the stress tensor. Using this approach, we perform a case study on a variety of condensed-phase systems (including liquid water, a benzene molecular crystal polymorph, and semi-conducting crystalline silicon) and demonstrate that the EXX contributions to the energy and stress tensor simultaneously converge with an appropriate choice of exx parameters. This is followed by a critical assessment of the computational performance of the extended exx module across several different high-performance computing architectures via case studies on (i) the computational complexity due to lattice symmetry during NpT simulations of three different ice polymorphs (i.e., ice Ih, II, and III) and (ii) the strong/weak parallel scaling during large-scale NpT simulations of liquid water. We demonstrate that the robust and highly scalable implementation of this approach in the extended exx module is capable of evaluating the EXX contribution to the stress tensor with negligible cost (<1%) as well as all other EXX-related quantities needed during NpT simulations of liquid water (with a very tight 150 Ry planewave cutoff) in ≈5.2 s ((H2O)128) and ≈6.8 s ((H2O)256) per AIMD step. As such, the extended exx module presented in this work brings us another step closer to routinely performing hybrid DFT-based AIMD simulations of sufficient duration for large-scale condensed-phase systems across a wide range of thermodynamic conditions.
Collapse
Affiliation(s)
- Hsin-Yu Ko
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Biswajit Santra
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
22
|
Zhang C, Tang F, Chen M, Xu J, Zhang L, Qiu DY, Perdew JP, Klein ML, Wu X. Modeling Liquid Water by Climbing up Jacob's Ladder in Density Functional Theory Facilitated by Using Deep Neural Network Potentials. J Phys Chem B 2021; 125:11444-11456. [PMID: 34533960 DOI: 10.1021/acs.jpcb.1c03884] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
Within the framework of Kohn-Sham density functional theory (DFT), the ability to provide good predictions of water properties by employing a strongly constrained and appropriately normed (SCAN) functional has been extensively demonstrated in recent years. Here, we further advance the modeling of water by building a more accurate model on the fourth rung of Jacob's ladder with the hybrid functional, SCAN0. In particular, we carry out both classical and Feynman path-integral molecular dynamics calculations of water with the SCAN0 functional and the isobaric-isothermal ensemble. To generate the equilibrated structure of water, a deep neural network potential is trained from the atomic potential energy surface based on ab initio data obtained from SCAN0 DFT calculations. For the electronic properties of water, a separate deep neural network potential is trained by using the Deep Wannier method based on the maximally localized Wannier functions of the equilibrated trajectory at the SCAN0 level. The structural, dynamic, and electric properties of water were analyzed. The hydrogen-bond structures, density, infrared spectra, diffusion coefficients, and dielectric constants of water, in the electronic ground state, are computed by using a large simulation box and long simulation time. For the properties involving electronic excitations, we apply the GW approximation within many-body perturbation theory to calculate the quasiparticle density of states and bandgap of water. Compared to the SCAN functional, mixing exact exchange mitigates the self-interaction error in the meta-generalized-gradient approximation and further softens liquid water toward the experimental direction. For most of the water properties, the SCAN0 functional shows a systematic improvement over the SCAN functional. However, some important discrepancies remain. The H-bond network predicted by the SCAN0 functional is still slightly overstructured compared to the experimental results.
Collapse
Affiliation(s)
- Chunyi Zhang
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Fujie Tang
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Mohan Chen
- HEDPS, Center for Applied Physics and Technology, College of Engineering, Peking University, Beijing 100871, China
| | - Jianhang Xu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, United States
| | - Diana Y Qiu
- Department of Mechanical Engineering and Materials Science, Yale University, New Haven, Connecticut 06520, United States
| | - John P Perdew
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States.,Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Michael L Klein
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States.,Institute for Computational Molecular Science, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| |
Collapse
|
23
|
Keith JA, Vassilev-Galindo V, Cheng B, Chmiela S, Gastegger M, Müller KR, Tkatchenko A. Combining Machine Learning and Computational Chemistry for Predictive Insights Into Chemical Systems. Chem Rev 2021; 121:9816-9872. [PMID: 34232033 PMCID: PMC8391798 DOI: 10.1021/acs.chemrev.1c00107] [Citation(s) in RCA: 190] [Impact Index Per Article: 63.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Indexed: 12/23/2022]
Abstract
Machine learning models are poised to make a transformative impact on chemical sciences by dramatically accelerating computational algorithms and amplifying insights available from computational chemistry methods. However, achieving this requires a confluence and coaction of expertise in computer science and physical sciences. This Review is written for new and experienced researchers working at the intersection of both fields. We first provide concise tutorials of computational chemistry and machine learning methods, showing how insights involving both can be achieved. We follow with a critical review of noteworthy applications that demonstrate how computational chemistry and machine learning can be used together to provide insightful (and useful) predictions in molecular and materials modeling, retrosyntheses, catalysis, and drug design.
Collapse
Affiliation(s)
- John A. Keith
- Department
of Chemical and Petroleum Engineering Swanson School of Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Valentin Vassilev-Galindo
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Bingqing Cheng
- Accelerate
Programme for Scientific Discovery, Department
of Computer Science and Technology, 15 J. J. Thomson Avenue, Cambridge CB3 0FD, United Kingdom
| | - Stefan Chmiela
- Department
of Software Engineering and Theoretical Computer Science, Technische Universität Berlin, 10587, Berlin, Germany
| | - Michael Gastegger
- Department
of Software Engineering and Theoretical Computer Science, Technische Universität Berlin, 10587, Berlin, Germany
| | - Klaus-Robert Müller
- Machine
Learning Group, Technische Universität
Berlin, 10587, Berlin, Germany
- Department
of Artificial Intelligence, Korea University, Anam-dong, Seongbuk-gu, Seoul, 02841, Korea
- Max-Planck-Institut für Informatik, 66123 Saarbrücken, Germany
- Google Research, Brain Team, 10117 Berlin, Germany
| | - Alexandre Tkatchenko
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
24
|
Zhang L, Wang H, Car R, E W. Phase Diagram of a Deep Potential Water Model. PHYSICAL REVIEW LETTERS 2021; 126:236001. [PMID: 34170175 DOI: 10.1103/physrevlett.126.236001] [Citation(s) in RCA: 91] [Impact Index Per Article: 30.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Accepted: 04/28/2021] [Indexed: 06/13/2023]
Abstract
Using the Deep Potential methodology, we construct a model that reproduces accurately the potential energy surface of the SCAN approximation of density functional theory for water, from low temperature and pressure to about 2400 K and 50 GPa, excluding the vapor stability region. The computational efficiency of the model makes it possible to predict its phase diagram using molecular dynamics. Satisfactory overall agreement with experimental results is obtained. The fluid phases, molecular and ionic, and all the stable ice polymorphs, ordered and disordered, are predicted correctly, with the exception of ice III and XV that are stable in experiments, but metastable in the model. The evolution of the atomic dynamics upon heating, as ice VII transforms first into ice VII^{''} and then into an ionic fluid, reveals that molecular dissociation and breaking of the ice rules coexist with strong covalent fluctuations, explaining why only partial ionization was inferred in experiments.
Collapse
Affiliation(s)
- Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA
| | - Han Wang
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Fenghao East Road 2, Beijing 100094, People's Republic of China
| | - Roberto Car
- Department of Chemistry, Department of Physics, Program in Applied and Computational Mathematics, Princeton Institute for the Science and Technology of Materials, Princeton University, Princeton, New Jersey 08544, USA
| | - Weinan E
- Department of Mathematics and Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA and Beijing Institute of Big Data Research, Beijing 100871, People's Republic of China
| |
Collapse
|
25
|
Liu X, Zhang L, Liu J. Machine learning phase space quantum dynamics approaches. J Chem Phys 2021; 154:184104. [PMID: 34241027 DOI: 10.1063/5.0046689] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023] Open
Abstract
Derived from phase space expressions of the quantum Liouville theorem, equilibrium continuity dynamics is a category of trajectory-based phase space dynamics methods, which satisfies the two critical fundamental criteria: conservation of the quantum Boltzmann distribution for the thermal equilibrium system and being exact for any thermal correlation functions (even of nonlinear operators) in the classical and harmonic limits. The effective force and effective mass matrix are important elements in the equations of motion of equilibrium continuity dynamics, where only the zeroth term of an exact series expansion of the phase space propagator is involved. We introduce a machine learning approach for fitting these elements in quantum phase space, leading to a much more efficient integration of the equations of motion. Proof-of-concept applications to realistic molecules demonstrate that machine learning phase space dynamics approaches are possible as well as competent in producing reasonably accurate results with a modest computation effort.
Collapse
Affiliation(s)
- Xinzijian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| | - Linfeng Zhang
- Beijing Institute of Big Data Research, Beijing 100871, China
| | - Jian Liu
- Beijing National Laboratory for Molecular Sciences, Institute of Theoretical and Computational Chemistry, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
26
|
Yang M, Bonati L, Polino D, Parrinello M. Using metadynamics to build neural network potentials for reactive events: the case of urea decomposition in water. Catal Today 2021. [DOI: 10.1016/j.cattod.2021.03.018] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
|
27
|
Liang W, Lu G, Yu J. Machine-Learning-Driven Simulations on Microstructure and Thermophysical Properties of MgCl 2-KCl Eutectic. ACS APPLIED MATERIALS & INTERFACES 2021; 13:4034-4042. [PMID: 33430593 DOI: 10.1021/acsami.0c20665] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Theoretical studies on the MgCl2-KCl eutectic heavily rely on ab initio calculations based on density functional theory (DFT). However, neither large-scale nor long-time calculations are feasible in the framework of the ab initio method, which makes it challenging to accurately predict some properties. To address this issue, a scheme based on ab initio calculation, deep neural networks, and machine learning is introduced. By training on high-quality data sets generated by ab initio calculations, a deep potential (DP) is constructed to describe the interaction between atoms. This work shows that the DP enables higher efficiency and similar accuracy relative to DFT. By performing molecular dynamics simulations with DP, the microstructure and thermophysical properties of the MgCl2-KCl eutectic (32:68 mol %) are investigated. The structural evolution with temperature is analyzed through partial radial distribution functions, coordination numbers, angular distribution functions, and structural factors. Meanwhile, the estimated thermophysical properties are discussed, including density, thermal expansion coefficient, shear viscosity, self-diffusion coefficient, and specific heat capacity. It reveals that the Mg2+ ions in this system have a distorted tetrahedral geometry rather than an octahedral one (with vacancies). The microstructure of the MgCl2-KCl eutectic shows the feature of medium-range order, and this feature will be enhanced at a higher temperature. All predicted thermophysical properties are in good agreement with the experimental results. The hydrodynamic radius determined from the shear viscosity and self-diffusion coefficient shows that the Mg2+ ions have a strong local structure and diffuse as if with an intact coordination shell. Overall, this work provides a thorough understanding of the microstructure and enriches the data of the thermophysical properties of the MgCl2-KCl eutectic.
Collapse
Affiliation(s)
- Wenshuo Liang
- School of Resources and Environmental Engineering, East China University of Science and Technology, Shanghai 200237, China
- National Engineering Research Center for Integrated Utilization of Salt Lake Resource, East China University of Science and Technology, Shanghai 200237, China
| | - Guimin Lu
- School of Resources and Environmental Engineering, East China University of Science and Technology, Shanghai 200237, China
- National Engineering Research Center for Integrated Utilization of Salt Lake Resource, East China University of Science and Technology, Shanghai 200237, China
| | - Jianguo Yu
- School of Resources and Environmental Engineering, East China University of Science and Technology, Shanghai 200237, China
- National Engineering Research Center for Integrated Utilization of Salt Lake Resource, East China University of Science and Technology, Shanghai 200237, China
| |
Collapse
|
28
|
Yue S, Muniz MC, Calegari Andrade MF, Zhang L, Car R, Panagiotopoulos AZ. When do short-range atomistic machine-learning models fall short? J Chem Phys 2021; 154:034111. [DOI: 10.1063/5.0031215] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Affiliation(s)
- Shuwen Yue
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | - Maria Carolina Muniz
- Department of Chemical and Biological Engineering, Princeton University, Princeton, New Jersey 08544, USA
| | | | - Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University, Princeton, New Jersey 08544, USA
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, USA
| | | |
Collapse
|
29
|
Vinson J. Faster exact exchange in periodic systems using single-precision arithmetic. J Chem Phys 2020; 153:204106. [PMID: 33261476 DOI: 10.1063/5.0030493] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
Density-functional theory simplifies many-electron calculations by approximating the exchange and correlation interactions with a one-electron operator that is a functional of the density. Hybrid functionals incorporate some amount of exact exchange, improving agreement with measured electronic and structural properties. However, calculations with hybrid functionals require substantial computational resources, limiting their use. By calculating the exchange interaction of periodic systems with single-precision arithmetic, the computation time is cut nearly in half with a negligible loss in accuracy. This improvement makes exact exchange calculations quicker and more feasible, especially for high-throughput calculations. Example hybrid density-functional theory calculations of band energies, forces, and x-ray absorption spectra show that this single-precision implementation maintains accuracy with significantly reduced runtime and memory requirements.
Collapse
Affiliation(s)
- John Vinson
- Material Measurement Laboratory, National Institute of Standards and Technology, 100 Bureau Drive, Gaithersburg, Maryland 20899, USA
| |
Collapse
|
30
|
Liang W, Lu G, Yu J. Molecular Dynamics Simulations of Molten Magnesium Chloride Using Machine‐Learning‐Based Deep Potential. ADVANCED THEORY AND SIMULATIONS 2020. [DOI: 10.1002/adts.202000180] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Affiliation(s)
- Wenshuo Liang
- School of Resources and Environmental Engineering East China University of Science and Technology Shanghai 200237 China
- National Engineering Research Center for Integrated Utilization of Salt Lake Resource East China University of Science and Technology Shanghai 200237 China
| | - Guimin Lu
- School of Resources and Environmental Engineering East China University of Science and Technology Shanghai 200237 China
- National Engineering Research Center for Integrated Utilization of Salt Lake Resource East China University of Science and Technology Shanghai 200237 China
| | - Jianguo Yu
- School of Resources and Environmental Engineering East China University of Science and Technology Shanghai 200237 China
- National Engineering Research Center for Integrated Utilization of Salt Lake Resource East China University of Science and Technology Shanghai 200237 China
| |
Collapse
|
31
|
Martelli F, Leoni F, Sciortino F, Russo J. Connection between liquid and non-crystalline solid phases in water. J Chem Phys 2020; 153:104503. [DOI: 10.1063/5.0018923] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Affiliation(s)
- Fausto Martelli
- IBM Research Europe, Hartree Centre, Daresbury WA4 4AD, United Kingdom
- School of Mathematics, University of Bristol, Bristol BS8 1UG, United Kingdom
| | - Fabio Leoni
- School of Mathematics, University of Bristol, Bristol BS8 1UG, United Kingdom
- Department of Physics, Sapienza University of Rome, P.le Aldo Moro 5, 00185 Rome, Italy
| | - Francesco Sciortino
- Department of Physics, Sapienza University of Rome, P.le Aldo Moro 5, 00185 Rome, Italy
| | - John Russo
- School of Mathematics, University of Bristol, Bristol BS8 1UG, United Kingdom
- Department of Physics, Sapienza University of Rome, P.le Aldo Moro 5, 00185 Rome, Italy
| |
Collapse
|
32
|
Ko HY, Jia J, Santra B, Wu X, Car R, DiStasio RA. Enabling Large-Scale Condensed-Phase Hybrid Density Functional Theory Based Ab Initio Molecular Dynamics. 1. Theory, Algorithm, and Performance. J Chem Theory Comput 2020; 16:3757-3785. [PMID: 32045232 DOI: 10.1021/acs.jctc.9b01167] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
By including a fraction of exact exchange (EXX), hybrid functionals reduce the self-interaction error in semilocal density functional theory (DFT) and thereby furnish a more accurate and reliable description of the underlying electronic structure in systems throughout biology, chemistry, physics, and materials science. However, the high computational cost associated with the evaluation of all required EXX quantities has limited the applicability of hybrid DFT in the treatment of large molecules and complex condensed-phase materials. To overcome this limitation, we describe a linear-scaling approach that utilizes a local representation of the occupied orbitals (e.g., maximally localized Wannier functions (MLWFs)) to exploit the sparsity in the real-space evaluation of the quantum mechanical exchange interaction in finite-gap systems. In this work, we present a detailed description of the theoretical and algorithmic advances required to perform MLWF-based ab initio molecular dynamics (AIMD) simulations of large-scale condensed-phase systems of interest at the hybrid DFT level. We focus our theoretical discussion on the integration of this approach into the framework of Car-Parrinello AIMD, and highlight the central role played by the MLWF-product potential (i.e., the solution of Poisson's equation for each corresponding MLWF-product density) in the evaluation of the EXX energy and wave function forces. We then provide a comprehensive description of the exx algorithm implemented in the open-source Quantum ESPRESSO program, which employs a hybrid MPI/OpenMP parallelization scheme to efficiently utilize the high-performance computing (HPC) resources available on current- and next-generation supercomputer architectures. This is followed by a critical assessment of the accuracy and parallel performance (e.g., strong and weak scaling) of this approach when AIMD simulations of liquid water are performed in the canonical (NVT) ensemble. With access to HPC resources, we demonstrate that exx enables hybrid DFT-based AIMD simulations of condensed-phase systems containing 500-1000 atoms (e.g., (H2O)256) with a wall time cost that is comparable to that of semilocal DFT. In doing so, exx takes us one step closer to routinely performing AIMD simulations of complex and large-scale condensed-phase systems for sufficiently long time scales at the hybrid DFT level of theory.
Collapse
Affiliation(s)
- Hsin-Yu Ko
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States.,Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States
| | - Junteng Jia
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| | - Biswajit Santra
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States.,Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Xifan Wu
- Department of Physics, Temple University, Philadelphia, Pennsylvania 19122, United States
| | - Roberto Car
- Department of Chemistry, Princeton University, Princeton, New Jersey 08544, United States.,Department of Physics, Princeton University, Princeton, New Jersey 08544, United States
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, United States
| |
Collapse
|
33
|
Calegari Andrade MF, Ko HY, Zhang L, Car R, Selloni A. Free energy of proton transfer at the water-TiO 2 interface from ab initio deep potential molecular dynamics. Chem Sci 2020; 11:2335-2341. [PMID: 34084393 PMCID: PMC8157430 DOI: 10.1039/c9sc05116c] [Citation(s) in RCA: 85] [Impact Index Per Article: 21.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
TiO2 is a widely used photocatalyst in science and technology and its interface with water is important in fields ranging from geochemistry to biomedicine. Yet, it is still unclear whether water adsorbs in molecular or dissociated form on TiO2 even for the case of well-defined crystalline surfaces. To address this issue, we simulated the TiO2-water interface using molecular dynamics with an ab initio-based deep neural network potential. Our simulations show a dynamical equilibrium of molecular and dissociative adsorption of water on TiO2. Water dissociates through a solvent-assisted concerted proton transfer to form a pair of short-lived hydroxyl groups on the TiO2 surface. Molecular adsorption of water is ΔF = 8.0 ± 0.9 kJ mol-1 lower in free energy than the dissociative adsorption, giving rise to a 5.6 ± 0.5% equilibrium water dissociation fraction at room temperature. Due to the relevance of surface hydroxyl groups to the surface chemistry of TiO2, our model might be key to understanding phenomena ranging from surface functionalization to photocatalytic mechanisms.
Collapse
Affiliation(s)
| | - Hsin-Yu Ko
- Department of Chemistry, Princeton University Princeton NJ 08544 USA
| | - Linfeng Zhang
- Program in Applied and Computational Mathematics, Princeton University Princeton NJ 08544 USA
| | - Roberto Car
- Department of Chemistry, Princeton University Princeton NJ 08544 USA
| | - Annabella Selloni
- Department of Chemistry, Princeton University Princeton NJ 08544 USA
| |
Collapse
|