1
|
Mahajan A, Kurian JS, Lee J, Reichman DR, Sharma S. Response properties in phaseless auxiliary field quantum Monte Carlo. J Chem Phys 2023; 159:184101. [PMID: 37937933 DOI: 10.1063/5.0171996] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Accepted: 10/11/2023] [Indexed: 11/09/2023] Open
Abstract
We present a method for calculating first-order response properties in phaseless auxiliary field quantum Monte Carlo by applying automatic differentiation (AD). Biases and statistical efficiency of the resulting estimators are discussed. Our approach demonstrates that AD enables the calculation of reduced density matrices with the same computational cost scaling per sample as energy calculations, accompanied by a cost prefactor of less than four in our numerical calculations. We investigate the role of self-consistency and trial orbital choice in property calculations. We find that orbitals obtained using density functional theory perform well for the dipole moments of selected molecules compared to those optimized self-consistently.
Collapse
Affiliation(s)
- Ankit Mahajan
- Department of Chemistry, Columbia University, New York, New York 10027, USA
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, USA
| | - Jo S Kurian
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, USA
| | - Joonho Lee
- Department of Chemistry, Columbia University, New York, New York 10027, USA
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - David R Reichman
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Sandeep Sharma
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, USA
| |
Collapse
|
2
|
Tiihonen J, Häkkinen H. Towards structural optimization of gold nanoclusters with quantum Monte Carlo. J Chem Phys 2023; 159:174301. [PMID: 37909449 DOI: 10.1063/5.0174383] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Accepted: 10/09/2023] [Indexed: 11/03/2023] Open
Abstract
We study the prospects of using quantum Monte Carlo techniques (QMC) to optimize the electronic wavefunctions and atomic geometries of gold compounds. Complex gold nanoclusters are widely studied for diverse biochemical applications, but the dynamic correlation and relativistic effects in gold set the bar high for reliable, predictive simulation methods. Here we study selected ground state properties of few-atom gold clusters by using density functional theory (DFT) and various implementations of the variational Monte Carlo (VMC) and diffusion Monte Carlo. We show that the QMC methods mitigate the exchange-correlation (XC) approximation made in the DFT approach: the average QMC results are more accurate and significantly more consistent than corresponding DFT results based on different XC functionals. Furthermore, we use demonstrate structural optimization of selected thiolated gold clusters with between 1 and 3 gold atoms using VMC forces. The optimization workflow is demonstrably consistent, robust, and its computational cost scales with nb, where b < 3 and n is the system size. We discuss the implications of these results while laying out steps for further developments.
Collapse
Affiliation(s)
- Juha Tiihonen
- Department of Physics, Nanoscience Center, University of Jyväskylä, Jyväskylä, Finland
| | - Hannu Häkkinen
- Department of Physics, Nanoscience Center, University of Jyväskylä, Jyväskylä, Finland
- Department of Chemistry, Nanoscience Center, University of Jyväskylä, Jyväskylä, Finland
| |
Collapse
|
3
|
Abstract
Diffusion Monte Carlo (DMC) is one of the most accurate techniques available for calculating the electronic properties of molecules and materials, yet it often remains a challenge to economically compute forces using this technique. As a result, ab initio molecular dynamics simulations and geometry optimizations that employ Diffusion Monte Carlo forces are often out of reach. One potential approach for accelerating the computation of "DMC forces" is to machine learn these forces from DMC energy calculations. In this work, we employ Behler-Parrinello Neural Networks to learn DMC forces from DMC energy calculations for geometry optimization and molecular dynamics simulations of small molecules. We illustrate the unique challenges that stem from learning forces without explicit force data and from noisy energy data by making rigorous comparisons of potential energy surface, dynamics, and optimization predictions among ab initio density functional theory (DFT) simulations and machine-learning models trained on DFT energies with forces, DFT energies without forces, and DMC energies without forces. We show for three small molecules─C2, H2O, and CH3Cl─that machine-learned DMC dynamics can reproduce average bond lengths and angles within a few percent of known experimental results at one hundredth of the typical cost. Our work describes a much-needed means of performing dynamics simulations on high-accuracy, DMC PESs and for generating DMC-quality molecular geometries given current algorithmic constraints.
Collapse
Affiliation(s)
- Cancan Huang
- Department of Chemistry, Brown University, Providence, Rhode Island02912, United States
| | - Brenda M Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island02912, United States
| |
Collapse
|
4
|
Jiang T, Fang W, Alavi A, Chen J. General Analytical Nuclear Forces and Molecular Potential Energy Surface from Full Configuration Interaction Quantum Monte Carlo. J Chem Theory Comput 2022; 18:7233-7242. [PMID: 36326847 DOI: 10.1021/acs.jctc.2c00440] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
The full configuration interaction quantum Monte Carlo (FCIQMC) is a state-of-the-art stochastic electronic structure method, providing a methodology to compute FCI-level state energies of molecular systems within a quantum chemical basis. However, especially to probe dynamics at the FCIQMC level, it is necessary to devise more efficient schemes to produce nuclear forces and potential energy surfaces (PES) from FCIQMC. In this work, we derive the general formula for nuclear forces from FCIQMC, and clarify different contributions of the total force. This method to obtain FCIQMC forces eliminates previous restrictions and can be used with frozen core approximation and free selection of orbitals, making it promising for more efficient nuclear forces calculations. After some numerical checks of this procedure on the binding curve of N2 molecule, we use the FCIQMC energy and force to obtain the full-dimensional ground state PES of the water molecule via Gaussian processes regression. The new water FCIQMC PES can be used as the basis for H2O ground state nuclear dynamics, structure optimization, and rotation-vibrational spectrum calculation.
Collapse
Affiliation(s)
- Tonghuan Jiang
- School of Physics, Peking University, Beijing100871, P. R. China
| | - Wei Fang
- State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian116023, P. R. China.,Department of Chemistry, Fudan University, Shanghai200438, P. R. China
| | - Ali Alavi
- Max Planck Institute for Solid State Research, Heisenbergstrasse 1, 70569Stuttgart, Germany.,University of Cambridge, Lensfield Road, CambridgeCB2 1EW, United Kingdom
| | - Ji Chen
- School of Physics, Peking University, Beijing100871, P. R. China.,Collaborative Innovation Center of Quantum Matter, Beijing100871, P. R. China.,Interdisciplinary Institute of Light-Element Quantum Materials and Research Center for Light-Element Advanced Materials, Peking University, Beijing100871, P. R. China.,Frontiers Science Center for Nano-Optoelectronics, Peking University, Beijing100871, P. R. China
| |
Collapse
|
5
|
Tiihonen J, Kent PRC, Krogel JT. Surrogate Hessian accelerated structural optimization for stochastic electronic structure theories. J Chem Phys 2022; 156:054104. [DOI: 10.1063/5.0079046] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Juha Tiihonen
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Paul R. C. Kent
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Tennessee 37831, USA
| | - Jaron T. Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| |
Collapse
|
6
|
van Rhijn J, Filippi C, De Palo S, Moroni S. Energy Derivatives in Real-Space Diffusion Monte Carlo. J Chem Theory Comput 2021; 18:118-123. [PMID: 34930005 PMCID: PMC8757439 DOI: 10.1021/acs.jctc.1c00496] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/04/2022]
Abstract
![]()
We present unbiased,
finite-variance estimators of energy derivatives
for real-space diffusion Monte Carlo calculations within the fixed-node
approximation. The derivative dλE is fully consistent with the dependence E(λ)
of the energy computed with the same time step. We address the issue
of the divergent variance of derivatives related to variations of
the nodes of the wave function both by using a regularization for
wave function parameter gradients recently proposed in variational
Monte Carlo and by introducing a regularization based on a coordinate
transformation. The essence of the divergent variance problem is distilled
into a particle-in-a-box toy model, where we demonstrate the algorithm.
Collapse
Affiliation(s)
- Jesse van Rhijn
- MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Claudia Filippi
- MESA+ Institute for Nanotechnology, University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands
| | - Stefania De Palo
- CNR-IOM DEMOCRITOS, Istituto Officina dei Materiali, and SISSA Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, I-34136 Trieste, Italy
| | - Saverio Moroni
- CNR-IOM DEMOCRITOS, Istituto Officina dei Materiali, and SISSA Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, I-34136 Trieste, Italy
| |
Collapse
|
7
|
Tiihonen J, Clay RC, Krogel JT. Toward quantum Monte Carlo forces on heavier ions: Scaling properties. J Chem Phys 2021; 154:204111. [PMID: 34241166 DOI: 10.1063/5.0052266] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Quantum Monte Carlo (QMC) forces have been studied extensively in recent decades because of their importance with spectroscopic observables and geometry optimization. Here, we benchmark the accuracy and computational cost of QMC forces. The zero-variance zero-bias (ZVZB) force estimator is used in standard variational and diffusion Monte Carlo simulations with mean-field based trial wavefunctions and atomic pseudopotentials. Statistical force uncertainties are obtained with a recently developed regression technique for heavy tailed QMC data [P. Lopez Rios and G. J. Conduit, Phys. Rev. E 99, 063312 (2019)]. By considering selected atoms and dimers with elements ranging from H to Zn (1 ≤ Zeff ≤ 20), we assess the accuracy and the computational cost of ZVZB forces as the effective pseudopotential valence charge, Zeff, increases. We find that the costs of QMC energies and forces approximately follow simple power laws in Zeff. The force uncertainty grows more rapidly, leading to a best case cost scaling relationship of approximately Zeff 6.5(3) for diffusion Monte Carlo. We find that the accessible system size at fixed computational cost scales as Zeff -2, insensitive to model assumptions or the use of the "space warp" variance-reduction technique. Our results predict the practical cost of obtaining forces for a range of materials, such as transition metal oxides where QMC forces have yet to be applied, and underscore the importance of further developing force variance-reduction techniques, particularly for atoms with high Zeff.
Collapse
Affiliation(s)
- Juha Tiihonen
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - Raymond C Clay
- Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| |
Collapse
|
8
|
Powell AD, Kroes GJ, Doblhoff-Dier K. Quantum Monte Carlo calculations on dissociative chemisorption of H2 + Al(110): Minimum barrier heights and their comparison to DFT values. J Chem Phys 2020; 153:224701. [DOI: 10.1063/5.0022919] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Andrew D. Powell
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| | - Geert-Jan Kroes
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| | - Katharina Doblhoff-Dier
- Leiden Institute of Chemistry, Gorlaeus Laboratories, Leiden University, P.O. Box 9502, 2300 RA Leiden, Netherlands
| |
Collapse
|
9
|
Needs RJ, Towler MD, Drummond ND, López Ríos P, Trail JR. Variational and diffusion quantum Monte Carlo calculations with the CASINO code. J Chem Phys 2020; 152:154106. [DOI: 10.1063/1.5144288] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- R. J. Needs
- TCM Group, Cavendish Laboratory, University of Cambridge, 19 J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - M. D. Towler
- University College London, London WC1E 6BT, United Kingdom
| | - N. D. Drummond
- Department of Physics, Lancaster University, Lancaster LA1 4YB, United Kingdom
| | - P. López Ríos
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - J. R. Trail
- TCM Group, Cavendish Laboratory, University of Cambridge, 19 J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
10
|
Ríos PL, Conduit GJ. Tail-regression estimator for heavy-tailed distributions of known tail indices and its application to continuum quantum Monte Carlo data. Phys Rev E 2019; 99:063312. [PMID: 31330629 DOI: 10.1103/physreve.99.063312] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2019] [Indexed: 11/07/2022]
Abstract
Standard statistical analysis is unable to provide reliable confidence intervals on expectation values of probability distributions that do not satisfy the conditions of the central limit theorem. We present a regression-based estimator of an arbitrary moment of a probability distribution with power-law heavy tails that exploits knowledge of the exponents of its asymptotic decay to bypass this issue entirely. Our method is applied to synthetic data and to energy and atomic force data from variational and diffusion quantum Monte Carlo calculations, whose distributions have known asymptotic forms [J. R. Trail, Phys. Rev. E 77, 016703 (2008)PLEEE81539-375510.1103/PhysRevE.77.016703; A. Badinski et al., J. Phys.: Condens. Matter 22, 074202 (2010)JCOMEL0953-898410.1088/0953-8984/22/7/074202]. We obtain convergent, accurate confidence intervals on the variance of the local energy of an electron gas and on the Hellmann-Feynman force on an atom in the all-electron carbon dimer. In each of these cases the uncertainty on our estimator is 45% and 60 times smaller, respectively, than the nominal (ill-defined) standard error.
Collapse
Affiliation(s)
- Pablo López Ríos
- Max-Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany.,Theory of Condensed Matter Group, Cavendish Laboratory, 19 J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - Gareth J Conduit
- Theory of Condensed Matter Group, Cavendish Laboratory, 19 J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
11
|
Liu YYF, Andrews B, Conduit GJ. Direct evaluation of the force constant matrix in quantum Monte Carlo. J Chem Phys 2019; 150:034104. [DOI: 10.1063/1.5070138] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Y. Y. F. Liu
- Theory of Condensed Matter Group, Cavendish Laboratory, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - B. Andrews
- Theory of Condensed Matter Group, Cavendish Laboratory, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| | - G. J. Conduit
- Theory of Condensed Matter Group, Cavendish Laboratory, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom
| |
Collapse
|
12
|
Crawford TD, Kumar A, Bazanté AP, Di Remigio R. Reduced‐scaling coupled cluster response theory: Challenges and opportunities. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1406] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Affiliation(s)
- T. Daniel Crawford
- Department of Chemistry Virginia Tech, Blacksburg Virginia
- The Molecular Sciences Software Institute Blacksburg Virginia
| | - Ashutosh Kumar
- Department of Chemistry Virginia Tech, Blacksburg Virginia
| | | | - Roberto Di Remigio
- Department of Chemistry Virginia Tech, Blacksburg Virginia
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry University of Tromsø ‐ The Arctic University of Norway Tromsø Norway
| |
Collapse
|
13
|
Cleland DM, Per MC. Performance of quantum Monte Carlo for calculating molecular bond lengths. J Chem Phys 2016; 144:124108. [DOI: 10.1063/1.4944826] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Deidre M. Cleland
- CSIRO Virtual Nanoscience Laboratory, 343 Royal Parade, Parkville, Victoria 3052, Australia
| | - Manolo C. Per
- CSIRO Virtual Nanoscience Laboratory, 343 Royal Parade, Parkville, Victoria 3052, Australia
| |
Collapse
|
14
|
Assaraf R, Caffarel M, Kollias AC. Chaotic versus nonchaotic stochastic dynamics in Monte Carlo simulations: a route for accurate energy differences in N-body systems. PHYSICAL REVIEW LETTERS 2011; 106:150601. [PMID: 21568537 DOI: 10.1103/physrevlett.106.150601] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/27/2010] [Indexed: 05/30/2023]
Abstract
We present a method to efficiently evaluate small energy differences of two close N-body systems by employing stochastic processes having a stability versus chaos property. By using the same random noise, energy differences are computed from close trajectories without reweighting procedures. The approach is presented for quantum systems but can be applied to classical N-body systems as well. It is exemplified with diffusion Monte Carlo simulations for long chains of hydrogen atoms and molecules for which it is shown that the long-standing problem of computing energy derivatives is solved.
Collapse
Affiliation(s)
- Roland Assaraf
- Laboratoire de Chimie Théorique, CNRS-UMR 7616 et Université Pierre et Marie Curie, 75252 Paris Cedex, France
| | | | | |
Collapse
|
15
|
Trail JR, Maezono R. Optimum and efficient sampling for variational quantum Monte Carlo. J Chem Phys 2010; 133:174120. [DOI: 10.1063/1.3488651] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|