1
|
Zhou X, Huang Z, He X. Diffusion Monte Carlo method for barrier heights of multiple proton exchanges and complexation energies in small water, ammonia, and hydrogen fluoride clusters. J Chem Phys 2024; 160:054103. [PMID: 38310472 DOI: 10.1063/5.0182164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2023] [Accepted: 01/09/2024] [Indexed: 02/05/2024] Open
Abstract
Proton exchange reactions are of key importance in many processes in water. However, it is nontrivial to achieve reliable barrier heights for multiple proton exchanges and complexation energies in hydrogen-bonded systems theoretically. Performance of the fixed-node diffusion quantum Monte Carlo (FN-DMC) with the single-Slater-Jastrow trial wavefunction on total energies, barrier heights of multiple proton exchanges, and complexation energies of small water, ammonia, and hydrogen fluoride clusters is investigated in this study. Effects of basis sets and those of locality approximation (LA), T-move approximation (T-move), and determinant localization approximation (DLA) schemes in dealing with the nonlocal part of pseudopotentials on FN-DMC results are evaluated. According to our results, diffuse basis functions are important in achieving reliable barrier heights and complexation energies with FN-DMC, although the cardinal number of the basis set is more important than diffuse basis functions on total energies of these systems. Our results also show that the time step bias with DLA and LA is smaller than T-move; however, the time step bias of DMC energies with respect to time steps using the T-move is roughly linear up to 0.06 a.u., while this is not the case with LA and DLA. Barrier heights and complexation energies with FN-DMC using these three schemes are always within chemical accuracy. Taking into account the fact that T-move and DLA are typically more stable than LA, FN-DMC calculations with the T-move or DLA scheme and basis sets containing diffuse basis functions are suggested for barrier heights of multiple proton exchanges and complexation energies of hydrogen-bonded clusters.
Collapse
Affiliation(s)
- Xiaojun Zhou
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, People's Republic of China
- School of Physics & Information Science, Shaanxi University of Science and Technology, Xi'an 710021, People's Republic of China
| | - Zhiru Huang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu, People's Republic of China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, Shanghai Frontiers Science Center of Molecule Intelligent Syntheses, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, People's Republic of China
- Chongqing Key Laboratory of Precision Optics, Chongqing Institute of East China Normal University, Chongqing 401120, People's Republic of China
- New York University-East China Normal University Center for Computational Chemistry, New York University Shanghai, Shanghai 200062, People's Republic of China
| |
Collapse
|
2
|
Herman KM, Xantheas SS. An extensive assessment of the performance of pairwise and many-body interaction potentials in reproducing ab initio benchmark binding energies for water clusters n = 2-25. Phys Chem Chem Phys 2023; 25:7120-7143. [PMID: 36853239 DOI: 10.1039/d2cp03241d] [Citation(s) in RCA: 10] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
We assess the performance of 7 pairwise additive (TIP3P, TIP4P, TIP4P-ice, TIP5P, OPC, SPC, SPC/E) and 8 families of many-body potentials (q-AQUA, HIPPO, AMOEBA, EFP, TTM, WHBB, MB-pol, MB-UCB) in reproducing high-level ab initio benchmark values, CCSD(T) or MP2 at the complete basis set (CBS) limit for the binding energy and the many-body expansion (MBE) of water clusters n = 2-11, 16-17, 20, 25. By including a large range of cluster sizes having dissimilar hydrogen bonding networks, we obtain an understanding of how these potentials perform for different hydrogen bonding arrangements that are mostly outside of their parameterization range. While it is appropriate to compare the results of ab initio based many-body potentials directly to the electronic binding energies (De's), the pairwise additive ones are compared to the enthalpies at T = 298 K, ΔH(298 K), as the latter class of force fields are parametrized to reproduce enthalpies (implicitly accounting for zero-point energy corrections) rather than binding energies. We find that all pairwise additive potentials considered overestimate the reference ΔH values for the n = 2-25 clusters by >13%. For the water dimer (n = 2) in particular, the errors are in the range 83-119% for the pairwise additive potentials studied since these are based on an effective rather than the true 2-body interaction specifically designed as a means of partially accounting for the missing many-body terms. This stronger 2-body interaction is achieved by an enhanced monomer dipole moment that mimics its increase from the gas phase monomer to the condensed phase value. Indeed, for cluster sizes n ≥ 4 the percent deviations become slightly smaller (albeit all exceeding 13%). In contrast, we find that the many-body potentials perform more accurately in reproducing the electronic binding energies (De's) throughout the entire cluster range (n = 2-25), all reproducing the ab initio benchmark binding energies within ±7% of the respective CBS values. We further assess the ability of a subset of the many-body potentials (MB-UCB, q-AQUA, MB-pol, and TTM2.1-F) to also reproduce the magnitude of the ab initio many-body energy terms for water cluster sizes n = 7, 10, 16 and 17. The potentials show an overall good agreement with the available benchmark values. However, we identify characteristic differences upon comparing the many-body terms at both the ab initio-optimized geometries and the respective potential-optimized geometries to the reference ab initio values. Additionally, by applying this analysis to a wide range of cluster sizes, trends in the MBE of the potentials with increasing cluster size can be identified. Finally, in an attempt to draw a parallel between the pairwise additive and many-body potentials, we report the analysis of the individual molecular dipole moments for water clusters with 1 to ∼4 solvation shells with the TTM2.1-F potential. We find that the internally solvated water molecules have in general a larger molecular dipole moment ranging from 2.6-3.0 D. This justifies the use of an enhanced, with respect to the gas-phase value, molecular dipole moment for the pairwise additive potentials, which is intended to fold in the many body terms into an effective (enhanced) pairwise interaction through the choice of the charges. These results have important implications for the development of future generations of efficient, transferable, and highly accurate classical interaction potentials in both the pairwise additive and many-body categories.
Collapse
Affiliation(s)
- Kristina M Herman
- Department of Chemistry, University of Washington, Seattle, WA 98195, USA.
| | - Sotiris S Xantheas
- Department of Chemistry, University of Washington, Seattle, WA 98195, USA. .,Advanced Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, 902 Battelle Boulevard, P.O. Box 999, MS K1-83, WA, 99352, USA.
| |
Collapse
|
3
|
Hill AN, Meijer AJHM, Hill JG. Correlation Consistent Basis Sets and Core Polarization Potentials for Al-Ar with ccECP Pseudopotentials. J Phys Chem A 2022; 126:5853-5863. [PMID: 35976118 PMCID: PMC9442647 DOI: 10.1021/acs.jpca.2c04446] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
New correlation consistent basis sets for the second-row
atoms
(Al–Ar) to be used with the neon-core correlation consistent
effective core potentials (ccECPs) have been developed. The basis
sets, denoted cc-pV(n+d)Z-ccECP (n = D, T, Q), include the “tight”-d functions that are
known to be important for second-row elements. Sets augmented with
additional diffuse functions are also reported. Effective core polarization
potentials (CPPs) to account for the effect of core–valence
correlation have been adjusted for the same elements, and two different
forms of the CPP cutoff function have been analyzed. The accuracy
of both the basis sets and the CPPs is assessed through benchmark
calculations at the coupled-cluster level of theory for atomic and
molecular properties. Agreement with all-electron results is much
improved relative to the basis sets that originally accompanied the
ccECPs; moreover, the combination of cc-pV(n+d)Z-ccECP
and CPPs is found to be a computationally efficient and accurate alternative
to including core electrons in the correlation treatment.
Collapse
Affiliation(s)
- Adam N Hill
- Department of Chemistry, University of Sheffield, Sheffield S3 7HF, United Kingdom
| | - Anthony J H M Meijer
- Department of Chemistry, University of Sheffield, Sheffield S3 7HF, United Kingdom
| | - J Grant Hill
- Department of Chemistry, University of Sheffield, Sheffield S3 7HF, United Kingdom
| |
Collapse
|
4
|
Yang DC, Kim DY, Kim KS. Quantum Monte Carlo Study of the Water Dimer Binding Energy and Halogen-π Interactions. J Phys Chem A 2019; 123:7785-7791. [PMID: 31418568 DOI: 10.1021/acs.jpca.9b04072] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Halogen-π systems are involved with competition between halogen bonding and π-interaction. Using the diffusion quantum Monte Carlo (DMC) method, we compare the interaction of benzene with diatomic halogens (X2: Cl2/Br2) with the typical hydrogen bonding in the water dimer, taking into account explicit correlations of up to three bodies. The benzene-Cl2/Br2 binding energies (13.07 ± 0.42/16.62 ± 0.02 kJ/mol) attributed to both halogen bonding and dispersion are smaller than but comparable to the typical hydrogen bonding in the water dimer binding energy (20.88 ± 0.27 kJ/mol). All of the above values are in good agreement with those from the coupled-cluster with single, double, and noniterative triple excitations (CCSD(T)) results at the complete basis set limit (benzene-Cl2/Br2: 12.78/16.17 kJ/mol; water dimer: 21.0 kJ/mol).
Collapse
Affiliation(s)
- D ChangMo Yang
- Center for Superfunctional Materials, Department of Chemistry , Ulsan National Institute of Science and Technology , Ulsan 44919 , Republic of Korea
| | - Dong Yeon Kim
- Center for Superfunctional Materials, Department of Chemistry , Ulsan National Institute of Science and Technology , Ulsan 44919 , Republic of Korea
| | - Kwang S Kim
- Center for Superfunctional Materials, Department of Chemistry , Ulsan National Institute of Science and Technology , Ulsan 44919 , Republic of Korea
| |
Collapse
|
5
|
Dubecký M. Noncovalent Interactions by Fixed-Node Diffusion Monte Carlo: Convergence of Nodes and Energy Differences vs Gaussian Basis-Set Size. J Chem Theory Comput 2017; 13:3626-3635. [DOI: 10.1021/acs.jctc.7b00537] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Matúš Dubecký
- Department of Physics, Faculty of Science, University of Ostrava, 30. dubna 22, 701
03 Ostrava, Czech Republic
- ATRI, Faculty of Materials
Science and Technology, Slovak University of Technology, Paulínska
16, 917 24 Trnava, Slovakia
| |
Collapse
|
6
|
Trail JR, Needs RJ. Shape and energy consistent pseudopotentials for correlated electron systems. J Chem Phys 2017; 146:204107. [PMID: 28571391 PMCID: PMC5444922 DOI: 10.1063/1.4984046] [Citation(s) in RCA: 31] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2017] [Accepted: 05/11/2017] [Indexed: 11/14/2022] Open
Abstract
A method is developed for generating pseudopotentials for use in correlated-electron calculations. The paradigms of shape and energy consistency are combined and defined in terms of correlated-electron wave-functions. The resulting energy consistent correlated electron pseudopotentials (eCEPPs) are constructed for H, Li-F, Sc-Fe, and Cu. Their accuracy is quantified by comparing the relaxed molecular geometries and dissociation energies which they provide with all electron results, with all quantities evaluated using coupled cluster singles, doubles, and triples calculations. Errors inherent in the pseudopotentials are also compared with those arising from a number of approximations commonly used with pseudopotentials. The eCEPPs provide a significant improvement in optimised geometries and dissociation energies for small molecules, with errors for the latter being an order-of-magnitude smaller than for Hartree-Fock-based pseudopotentials available in the literature. Gaussian basis sets are optimised for use with these pseudopotentials.
Collapse
Affiliation(s)
- J R Trail
- Theory of Condensed Matter Group, Cavendish Laboratory, J J Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - R J Needs
- Theory of Condensed Matter Group, Cavendish Laboratory, J J Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
7
|
Zhou X, Wang F. Barrier heights of hydrogen-transfer reactions with diffusion quantum monte carlo method. J Comput Chem 2017; 38:798-806. [DOI: 10.1002/jcc.24750] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2016] [Revised: 01/11/2017] [Accepted: 01/13/2017] [Indexed: 01/17/2023]
Affiliation(s)
- Xiaojun Zhou
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University; Chengdu 610065 People's Republic of China
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University; Chengdu 610065 People's Republic of China
| |
Collapse
|
8
|
Huang C, Muñoz-García AB, Pavone M. Effective scheme for partitioning covalent bonds in density-functional embedding theory: From molecules to extended covalent systems. J Chem Phys 2016; 145:244103. [DOI: 10.1063/1.4972012] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
|
10
|
Affiliation(s)
- Matúš Dubecký
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Lubos Mitas
- Department
of Physics and CHiPS, North Carolina State University, Raleigh, North Carolina 27695, United States
| | - Petr Jurečka
- Regional
Centre of Advanced Technologies and Materials, Department of Physical
Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| |
Collapse
|
11
|
Gasperich K, Jordan KD. Diffusion Monte Carlo Study of the Parallel Displaced Form of the Benzene Dimer. ACS SYMPOSIUM SERIES 2016. [DOI: 10.1021/bk-2016-1234.ch007] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Affiliation(s)
- Kevin Gasperich
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Kenneth D. Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
12
|
Exploration of Brueckner orbital trial wave functions in diffusion Monte Carlo calculations. Chem Phys Lett 2016. [DOI: 10.1016/j.cplett.2015.11.053] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
13
|
Deible MJ, Tuguldur O, Jordan KD. Theoretical Study of the Binding Energy of a Methane Molecule in a (H2O)20 Dodecahedral Cage. J Phys Chem B 2014; 118:8257-63. [DOI: 10.1021/jp501592h] [Citation(s) in RCA: 33] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Affiliation(s)
- Michael J. Deible
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Odbadrakh Tuguldur
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Kenneth D. Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
14
|
|
15
|
Zen A, Luo Y, Sorella S, Guidoni L. Molecular Properties by Quantum Monte Carlo: An Investigation on the Role of the Wave Function Ansatz and the Basis Set in the Water Molecule. J Chem Theory Comput 2013; 9:4332-4350. [PMID: 24526929 PMCID: PMC3920371 DOI: 10.1021/ct400382m] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Quantum Monte Carlo methods are accurate and promising many body techniques for electronic structure calculations which, in the last years, are encountering a growing interest thanks to their favorable scaling with the system size and their efficient parallelization, particularly suited for the modern high performance computing facilities. The ansatz of the wave function and its variational flexibility are crucial points for both the accurate description of molecular properties and the capabilities of the method to tackle large systems. In this paper, we extensively analyze, using different variational ansatzes, several properties of the water molecule, namely, the total energy, the dipole and quadrupole momenta, the ionization and atomization energies, the equilibrium configuration, and the harmonic and fundamental frequencies of vibration. The investigation mainly focuses on variational Monte Carlo calculations, although several lattice regularized diffusion Monte Carlo calculations are also reported. Through a systematic study, we provide a useful guide to the choice of the wave function, the pseudopotential, and the basis set for QMC calculations. We also introduce a new method for the computation of forces with finite variance on open systems and a new strategy for the definition of the atomic orbitals involved in the Jastrow-Antisymmetrised Geminal power wave function, in order to drastically reduce the number of variational parameters. This scheme significantly improves the efficiency of QMC energy minimization in case of large basis sets.
Collapse
Affiliation(s)
- Andrea Zen
- Dipartimento di Fisica, La Sapienza—Università di Roma, Piazzale Aldo Moro 2, 00185 Rome, Italy
| | - Ye Luo
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), via Bonomea 265, 34136 Trieste, Italy
| | - Sandro Sorella
- Scuola Internazionale Superiore di Studi Avanzati (SISSA), via Bonomea 265, 34136 Trieste, Italy
| | - Leonardo Guidoni
- Dipartimento di Scienze Fisiche e Chimiche, Università degli studi de L’Aquila, Via Vetoio, 67100 Coppito, L’Aquila, Italy
| |
Collapse
|
16
|
Dubecký M, Jurečka P, Derian R, Hobza P, Otyepka M, Mitas L. Quantum Monte Carlo Methods Describe Noncovalent Interactions with Subchemical Accuracy. J Chem Theory Comput 2013; 9:4287-92. [DOI: 10.1021/ct4006739] [Citation(s) in RCA: 76] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Matúš Dubecký
- Regional Centre of Advanced
Technologies and Materials, Department of
Physical Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Petr Jurečka
- Regional Centre of Advanced
Technologies and Materials, Department of
Physical Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - René Derian
- Institute
of Physics, Slovak Academy of Sciences, Dúbravská cesta 9, 845
11 Bratislava, Slovakia
| | - Pavel Hobza
- Regional Centre of Advanced
Technologies and Materials, Department of
Physical Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
- Institute
of Organic Chemistry and Biochemistry, Academy of Sciences of the Czech Republic, Flemingovo
nám. 2, 166 10 Prague 6, Czech Republic
| | - Michal Otyepka
- Regional Centre of Advanced
Technologies and Materials, Department of
Physical Chemistry, Faculty of Science, Palacký University Olomouc, tř.
17 listopadu 12, 771 46 Olomouc, Czech Republic
| | - Lubos Mitas
- Department
of Physics and CHiPS, North Carolina State University, Raleigh, North Carolina 27695, United States
| |
Collapse
|
17
|
Pozun ZD, Su X, Jordan KD. Establishing the Ground State of the Disjoint Diradical Tetramethyleneethane with Quantum Monte Carlo. J Am Chem Soc 2013; 135:13862-9. [DOI: 10.1021/ja406002n] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Zachary D. Pozun
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Xiaoge Su
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| | - Kenneth D. Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, United States
| |
Collapse
|
18
|
Wang FF, Deible MJ, Jordan KD. Benchmark Study of the Interaction Energy for an (H2O)16 Cluster: Quantum Monte Carlo and Complete Basis Set Limit MP2 Results. J Phys Chem A 2013; 117:7606-11. [DOI: 10.1021/jp404541c] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Fang-Fang Wang
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania
15260, United States
| | - Michael J. Deible
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania
15260, United States
| | - Kenneth D. Jordan
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania
15260, United States
| |
Collapse
|