1
|
Anderson TA, Umrigar CJ. Nonlocal pseudopotentials and time-step errors in diffusion Monte Carlo. J Chem Phys 2021; 154:214110. [PMID: 34240977 DOI: 10.1063/5.0052838] [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/15/2022] Open
Abstract
We present a version of the T-moves approach for treating nonlocal pseudopotentials in diffusion Monte Carlo, which has much smaller time-step errors than the existing T-moves approaches, while at the same time preserving desirable features such as the upper-bound property for the energy. In addition, we modify the reweighting factor of the projector used in diffusion Monte Carlo to reduce the time-step error. The latter is applicable not only to pseudopotential calculations but also to all-electron calculations.
Collapse
Affiliation(s)
- Tyler A Anderson
- Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, New York 14853, USA
| | - C J Umrigar
- Laboratory of Atomic and Solid State Physics, Cornell University, Ithaca, New York 14853, USA
| |
Collapse
|
2
|
Varandas AJC. Extrapolation in quantum chemistry: Insights on energetics and reaction dynamics. JOURNAL OF THEORETICAL & COMPUTATIONAL CHEMISTRY 2020. [DOI: 10.1142/s0219633620300013] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Since there is no exact solution for problems in physics and chemistry, extrapolation methods may assume a key role in quantitative quantum chemistry. Two topics where it bears considerable impact are addressed, both at the heart of computational quantum chemistry: electronic structure and reaction dynamics. In the first, the problem of extrapolating the energy obtained by solving the electronic Schrödinger equation to the limit of the complete one-electron basis set is addressed. With the uniform-singlet-and-triplet-extrapolation (USTE) scheme at the focal point, the emphasis is on recent updates covering from the energy itself to other molecular properties. The second topic refers to extrapolation of quantum mechanical reactive scattering probabilities from zero total angular momentum to any of the values that it may assume when running quasiclassical trajectories, QCT/QM-[Formula: see text]J. With the extrapolation guided in both cases by physically motivated asymptotic theories, realism is seeked by avoiding unsecure jumps into the unknown. Although, mostly review oriented, a few issues are addressed for the first time here and there. Prospects for future work conclude the overview.
Collapse
Affiliation(s)
- A. J. C. Varandas
- School of Physics and Physical Engineering, Qufu Normal University, Qufu 273165, P. R. China
- Department of Physics, Universidade Federal do Espírito Santo, Vitória 29075-910, Brazil
- Department of Chemistry and Chemistry Centre, University of Coimbra, Coimbra 3004-535, Portugal
| |
Collapse
|
3
|
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
|
4
|
Zen A, Coccia E, Gozem S, Olivucci M, Guidoni L. Quantum Monte Carlo Treatment of the Charge Transfer and Diradical Electronic Character in a Retinal Chromophore Minimal Model. J Chem Theory Comput 2016; 11:992-1005. [PMID: 25821414 PMCID: PMC4357234 DOI: 10.1021/ct501122z] [Citation(s) in RCA: 36] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2014] [Indexed: 01/22/2023]
Abstract
![]()
The
penta-2,4-dieniminium cation (PSB3) displays similar ground
state and first excited state potential energy features as those of
the retinal protonated Schiff base (RPSB) chromophore in rhodopsin.
Recently, PSB3 has been used to benchmark several electronic structure
methods, including highly correlated multireference wave function
approaches, highlighting the necessity to accurately describe the
electronic correlation in order to obtain reliable properties even
along the ground state (thermal) isomerization paths. In this work,
we apply two quantum Monte Carlo approaches, the variational Monte
Carlo and the lattice regularized diffusion Monte Carlo, to study
the energetics and electronic properties of PSB3 along representative
minimum energy paths and scans related to its thermal cis–trans isomerization. Quantum Monte Carlo
is used in combination with the Jastrow antisymmetrized geminal power
ansatz, which guarantees an accurate and balanced description of the
static electronic correlation thanks to the multiconfigurational nature
of the antisymmetrized geminal power term, and of the dynamical correlation,
due to the presence of the Jastrow factor explicitly depending on
electron–electron distances. Along the two ground state isomerization
minimum energy paths of PSB3, CASSCF calculations yield wave functions
having either charge transfer or diradical character in proximity
of the two transition state configurations. Here, we observe that
at the quantum Monte Carlo level of theory, only the transition state
with charge transfer character can be located. The conical intersection,
which becomes highly sloped, is observed only if the path connecting
the two original CASSCF transition states is extended beyond the diradical
one, namely by increasing the bond-length-alternation (BLA). These
findings are in good agreement with the results obtained by MRCISD+Q
calculations, and they demonstrate the importance of having an accurate
description of the static and dynamical correlation when studying
isomerization and transition states of conjugated systems.
Collapse
|
5
|
Sorella S, Devaux N, Dagrada M, Mazzola G, Casula M. Geminal embedding scheme for optimal atomic basis set construction in correlated calculations. J Chem Phys 2015; 143:244112. [DOI: 10.1063/1.4938089] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- S. Sorella
- International School for Advanced Studies (SISSA), Via Beirut 2-4, 34014 Trieste, Italy and INFM Democritos National Simulation Center, Trieste, Italy
| | - N. Devaux
- Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, Case 115, 4 Place Jussieu, 75252 Paris Cedex 05, France
| | - M. Dagrada
- Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, Case 115, 4 Place Jussieu, 75252 Paris Cedex 05, France
| | - G. Mazzola
- Theoretische Physik, ETH Zurich, 8093 Zurich, Switzerland
| | - M. Casula
- CNRS and Institut de Minéralogie, de Physique des Matériaux et de Cosmochimie, Université Pierre et Marie Curie, Case 115, 4 Place Jussieu, 75252 Paris Cedex 05, France
| |
Collapse
|
6
|
Wu Y, Wagner LK, Aluru NR. The interaction between hexagonal boron nitride and water from first principles. J Chem Phys 2015; 142:234702. [DOI: 10.1063/1.4922491] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Yanbin Wu
- Department of Mechanical Science and Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - Lucas K. Wagner
- Department of Physics, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801-3080, USA
| | - Narayana R. Aluru
- Department of Mechanical Science and Engineering, Beckman Institute for Advanced Science and Technology, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| |
Collapse
|
7
|
Zen A, Luo Y, Mazzola G, Guidoni L, Sorella S. Ab initio molecular dynamics simulation of liquid water by quantum Monte Carlo. J Chem Phys 2015; 142:144111. [PMID: 25877566 DOI: 10.1063/1.4917171] [Citation(s) in RCA: 55] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article, we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous density functional theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab initio simulations of complex chemical systems.
Collapse
Affiliation(s)
- Andrea Zen
- Dipartimento di Fisica, “La Sapienza” - Università di Roma, piazzale Aldo Moro 5, 00185 Rome, Italy
- London Centre for Nanotechnology, University College London, London WC1E 6BT, United Kingdom
| | - Ye Luo
- SISSA–International School for Advanced Studies, Via Bonomea 26, 34136 Trieste, Italy
- Democritos Simulation Center CNR–IOM Istituto Officina dei Materiali, 34151 Trieste, Italy
| | - Guglielmo Mazzola
- SISSA–International School for Advanced Studies, Via Bonomea 26, 34136 Trieste, Italy
- Democritos Simulation Center CNR–IOM Istituto Officina dei Materiali, 34151 Trieste, Italy
| | - Leonardo Guidoni
- Dipartimento di Fisica, “La Sapienza” - Università di Roma, piazzale Aldo Moro 5, 00185 Rome, Italy
- Dipartimento di Scienze Fisiche e Chimiche, Università degli Studi dell’ Aquila, via Vetoio, 67100 L’ Aquila, Italy
| | - Sandro Sorella
- SISSA–International School for Advanced Studies, Via Bonomea 26, 34136 Trieste, Italy
- Democritos Simulation Center CNR–IOM Istituto Officina dei Materiali, 34151 Trieste, Italy
| |
Collapse
|
8
|
Zen A, Coccia E, Luo Y, Sorella S, Guidoni L. Static and Dynamical Correlation in Diradical Molecules by Quantum Monte Carlo Using the Jastrow Antisymmetrized Geminal Power Ansatz. J Chem Theory Comput 2014; 10:1048-61. [DOI: 10.1021/ct401008s] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Affiliation(s)
- Andrea Zen
- Dipartimento
di Fisica, Sapienza-Università di Roma, Piazzale Aldo Moro
5, 00185 Rome, Italy
| | - Emanuele Coccia
- Dipartimento
di Scienze Fisiche e Chimiche, Università degli Studi dell’Aquila, Via Vetoio 2, 67010 L’Aquila, Italy
| | - Ye Luo
- Scuola
Internazionale Superiore di Studi Avanzati (SISSA) and Democritos
National Simulation Center, Istituto Officina dei Materiali del CNR, via Bonomea 265, 34136 Trieste, Italy
| | - Sandro Sorella
- Scuola
Internazionale Superiore di Studi Avanzati (SISSA) and Democritos
National Simulation Center, Istituto Officina dei Materiali del CNR, via Bonomea 265, 34136 Trieste, Italy
| | - Leonardo Guidoni
- Dipartimento
di Scienze Fisiche e Chimiche, Università degli Studi dell’Aquila, Via Vetoio 2, 67010 L’Aquila, Italy
| |
Collapse
|
9
|
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
|
10
|
Coccia E, Guidoni L. Quantum monte carlo study of the retinal minimal model C5H6NH2+. J Comput Chem 2012; 33:2332-9. [DOI: 10.1002/jcc.23071] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2012] [Revised: 06/22/2012] [Accepted: 06/26/2012] [Indexed: 12/22/2022]
|
11
|
Morales MA, McMinis J, Clark BK, Kim J, Scuseria GE. Multideterminant Wave Functions in Quantum Monte Carlo. J Chem Theory Comput 2012; 8:2181-8. [DOI: 10.1021/ct3003404] [Citation(s) in RCA: 92] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Miguel A. Morales
- Lawrence Livermore National
Laboratory, Livermore, California 94550, United States
| | - Jeremy McMinis
- Department of Physics, University
of Illinois at Urbana−Champaign, Urbana, Illinois 61801, United
States
| | - Bryan K. Clark
- Princeton Center For
Theoretical
Science and Department of Physics, Joseph Henry Laboratories, Princeton
University, Princeton, New Jersey 08544, United States
| | - Jeongnim Kim
- National Center for Supercomputing
Applications, University of Illinois at Urbana−Champaign, Urbana,
Illinois 61801, United States and Materials Science and Technology
Division and Computational Chemistry and Materials Division, Oak Ridge
National Laboratory, Oak Ridge, Tennessee 37830, United States
| | - Gustavo E. Scuseria
- Department of Chemistry and Department of Physics & Astronomy, Rice University, Houston, Texas 77005-1892, United States
| |
Collapse
|
12
|
Coccia E, Chernomor O, Barborini M, Sorella S, Guidoni L. Molecular Electrical Properties from Quantum Monte Carlo Calculations: Application to Ethyne. J Chem Theory Comput 2012; 8:1952-62. [DOI: 10.1021/ct300171q] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Emanuele Coccia
- Dipartimento di Scienze Fisiche
e Chimiche, Universitá degli Studi dell'Aquila, via Vetoio (Coppito), 67100, L'Aquila, Italy
| | - Olga Chernomor
- Università Degli Studi de L’Aquila, Dipartimento di Matematica
Pura ed Applicata, via Vetoio (Coppito 1), 67100 L’Aquila,
Italy
| | - Matteo Barborini
- Università Degli Studi de L’Aquila, Dipartimento di Matematica
Pura ed Applicata, via Vetoio (Coppito 1), 67100 L’Aquila,
Italy
| | - Sandro Sorella
- Scuola Internazionale
Superiore
di Studi Avanzati (SISSA) and Democritos National Simulation Center, Istituto Officina dei Materiali del CNR, via Bonomea
265, 34136 Trieste, Italy
| | - Leonardo Guidoni
- Dipartimento di Scienze Fisiche
e Chimiche, Universitá degli Studi dell'Aquila, via Vetoio (Coppito), 67100, L'Aquila, Italy
| |
Collapse
|
13
|
Petruzielo FR, Toulouse J, Umrigar CJ. Approaching chemical accuracy with quantum Monte Carlo. J Chem Phys 2012; 136:124116. [DOI: 10.1063/1.3697846] [Citation(s) in RCA: 63] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|