1
|
Wittmann L, Gordiy I, Friede M, Helmich-Paris B, Grimme S, Hansen A, Bursch M. Extension of the D3 and D4 London dispersion corrections to the full actinides series. Phys Chem Chem Phys 2024; 26:21379-21394. [PMID: 39092890 DOI: 10.1039/d4cp01514b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/04/2024]
Abstract
Efficient dispersion corrections are an indispensable component of modern density functional theory, semi-empirical quantum mechanical, and even force field methods. In this work, we extend the well established D3 and D4 London dispersion corrections to the full actinides series, francium, and radium. To keep consistency with the existing versions, the original parameterization strategy of the D4 model was only slightly modified. This includes improved reference Hirshfeld atomic partial charges at the ωB97M-V/ma-def-TZVP level to fit the required electronegativity equilibration charge (EEQ) model. In this context, we developed a new actinide data set called AcQM, which covers the most common molecular actinide compound space. Furthermore, the efficient calculation of dynamic polarizabilities that are needed to construct CAB6 dispersion coefficients was implemented into the ORCA program package. The extended models are assessed for the computation of dissociation curves of actinide atoms and ions, geometry optimizations of crystal structure cutouts, gas-phase structures of small uranium compounds, and an example extracted from a small actinide complex protein assembly. We found that the novel parameterizations perform on par with the computationally more demanding density-dependent VV10 dispersion correction. With the presented extension, the excellent cost-accuracy ratio of the D3 and D4 models can now be utilized in various fields of computational actinide chemistry and, e.g., in efficient composite DFT methods such as r2SCAN-3c. They are implemented in our freely available standalone codes (dftd4, s-dftd3) and the D4 version will be also available in the upcoming ORCA 6.0 program package.
Collapse
Affiliation(s)
- Lukas Wittmann
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Igor Gordiy
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Marvin Friede
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstr. 4, 53115 Bonn, Germany.
| | - Markus Bursch
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany.
- FACCTs GmbH, 50677, Köln, Germany
| |
Collapse
|
2
|
Pham HQ, Ouyang R, Lv D. Scalable Quantum Monte Carlo with Direct-Product Trial Wave Functions. J Chem Theory Comput 2024; 20:3524-3534. [PMID: 38700513 DOI: 10.1021/acs.jctc.3c00769] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/15/2024]
Abstract
The computational demand posed by applying multi-Slater determinant trials in phaseless auxiliary-field quantum Monte Carlo methods (MSD-AFQMC) is particularly significant for molecules exhibiting strong correlations. Here, we propose using direct-product wave functions as trials for MSD-AFQMC, aiming to reduce computational overhead by leveraging the compactness of multi-Slater determinant trials in direct-product form (DP-MSD). This efficiency arises when the active space can be divided into noncoupling subspaces, a condition we term "decomposable active space". By employing localized-active space self-consistent field wave functions as an example of such trials, we demonstrate our proposed approach across a range of molecular systems, each exhibiting varying degrees of complexity in their electronic structures. Our findings indicate that the compact DP-MSD trials can reduce computational costs substantially, by up to 36 times for the C2H6N4 molecule where the two double bonds between nitrogen N=N are clearly separated by a C-C single bond, while maintaining accuracy when active spaces are decomposable. In the case of larger systems such as the benzene dimer, characterized by weak coupling between the two monomers, we observed a decrease in computational cost compared to using a complete active space trial, yet we retained the same level of accuracy. However, for systems where these active subspaces strongly couple, a scenario we refer to as "strong subspace coupling", the method's accuracy decreases compared to that achieved with a complete active space approach. We anticipate that our method will be beneficial for systems with noncoupling to weakly coupling subspaces that require local multireference treatments.
Collapse
Affiliation(s)
- Hung Q Pham
- ByteDance Research, San Jose, California 95110, United States
| | | | - Dingshun Lv
- ByteDance Research, Zhonghang Plaza, No. 43, North third Ring West Road, 100098 Beijing, China
| |
Collapse
|
3
|
Weber JL, Vuong H, Friesner RA, Reichman DR. Expanding the Design Space of Constraints in Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2023; 19:7567-7576. [PMID: 37889331 DOI: 10.1021/acs.jctc.3c00654] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/28/2023]
Abstract
We formulate and characterize a new constraint for auxiliary-field quantum Monte Carlo (AFQMC) applicable for general fermionic systems, which allows for the accumulation of phase in the random walk but disallows walkers with a magnitude of phase greater than π with respect to the trial wave function. For short imaginary times, before walkers accumulate sizable phase values, this approach is equivalent to exact free projection, allowing one to observe the accumulation of bias associated with the constraint and thus estimate its magnitude a priori. We demonstrate the stability of this constraint over arbitrary imaginary times and system sizes, highlighting the removal of noise due to the fermionic sign problem. Benchmark total energies for a variety of weakly and strongly correlated molecular systems reveal a distinct bias with respect to standard phaseless AFQMC, with a comparative increase in accuracy given sufficient quality of the trial wave function for the set of studied cases. We then take this constraint, termed linecut AFQMC (lc-AFQMC), and systematically release it (lcR-AFQMC), providing a route to obtain a smooth bridge between constrained AFQMC and the exact free projection results.
Collapse
Affiliation(s)
- John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Hung Vuong
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
4
|
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
|
5
|
Neugebauer H, Vuong HT, Weber JL, Friesner RA, Shee J, Hansen A. Toward Benchmark-Quality Ab Initio Predictions for 3d Transition Metal Electrocatalysts: A Comparison of CCSD(T) and ph-AFQMC. J Chem Theory Comput 2023; 19:6208-6225. [PMID: 37655473 DOI: 10.1021/acs.jctc.3c00617] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/02/2023]
Abstract
Generating accurate ab initio ionization energies for transition metal complexes is an important step toward the accurate computational description of their electrocatalytic reactions. Benchmark-quality data is required for testing existing theoretical methods and developing new ones but is complicated to obtain for many transition metal compounds due to the potential presence of both strong dynamical and static electron correlation. In this regime, it is questionable whether the so-called gold standard, coupled cluster with singles, doubles, and perturbative triples (CCSD(T)), provides the desired level of accuracy─roughly 1-3 kcal/mol. In this work, we compiled a test set of 28 3d metal-containing molecules relevant to homogeneous electrocatalysis (termed 3dTMV) and computed their vertical ionization energies (ionization potentials) with CCSD(T) and phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) in the def2-SVP basis set. A substantial effort has been made to converge away the phaseless bias in the ph-AFQMC reference values. We assess a wide variety of multireference diagnostics and find that spin-symmetry breaking of the CCSD wave function and the PBE0 density functional correlate well with our analysis of multiconfigurational wave functions. We propose quantitative criteria based on symmetry breaking to delineate correlation regimes inside of which appropriately performed CCSD(T) can produce mean absolute deviations from the ph-AFQMC reference values of roughly 2 kcal/mol or less and outside of which CCSD(T) is expected to fail. We also present a preliminary assessment of density functional theory (DFT) functionals on the 3dTMV set.
Collapse
Affiliation(s)
- Hagen Neugebauer
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, D-53115 Bonn, Germany
| | - Hung T Vuong
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, Rice University, Houston, Texas 77005, United States
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Clausius Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, D-53115 Bonn, Germany
| |
Collapse
|
6
|
Shee J, Weber JL, Reichman DR, Friesner RA, Zhang S. On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond. J Chem Phys 2023; 158:140901. [PMID: 37061483 PMCID: PMC10089686 DOI: 10.1063/5.0134009] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/17/2023] Open
Abstract
Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within "chemical accuracy" (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob's ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 - N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.
Collapse
Affiliation(s)
- James Shee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, USA
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, USA
| |
Collapse
|
7
|
Lee J, Pham HQ, Reichman DR. Twenty Years of Auxiliary-Field Quantum Monte Carlo in Quantum Chemistry: An Overview and Assessment on Main Group Chemistry and Bond-Breaking. J Chem Theory Comput 2022; 18:7024-7042. [PMID: 36255074 DOI: 10.1021/acs.jctc.2c00802] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
In this work, we present an overview of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) approach from a computational quantum chemistry perspective and present a numerical assessment of its performance on main group chemistry and bond-breaking problems with a total of 1004 relative energies. While our benchmark study is somewhat limited, we make recommendations for the use of ph-AFQMC for general main-group chemistry applications. For systems where single determinant wave functions are qualitatively accurate, we expect the accuracy of ph-AFQMC in conjunction with a single-determinant trial wave function to be between that of coupled-cluster with singles and doubles (CCSD) and CCSD with perturbative triples (CCSD(T)). For these applications, ph-AFQMC should be a method of choice when canonical CCSD(T) is too expensive to run. For systems where multireference (MR) wave functions are needed for qualitative accuracy, ph-AFQMC is far more accurate than MR perturbation theory methods and competitive with MR configuration interaction (MRCI) methods. Due to the computational efficiency of ph-AFQMC compared to MRCI, we recommended ph-AFQMC as a method of choice for handling dynamic correlation in MR problems. We conclude with a discussion of important directions for future development of the ph-AFQMC approach.
Collapse
Affiliation(s)
- Joonho Lee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Hung Q Pham
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
8
|
Liebermann N, Ghanem K, Alavi A. Importance-sampling FCIQMC: solving weak sign-problem systems. J Chem Phys 2022; 157:124111. [DOI: 10.1063/5.0107317] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We investigate the exact FCIQMC algorithm (without the initiator approximation) applied to weak sign-problem fermionic systems, namely systems in which the energy gap to the corresponding sign-free or ``stoquastized" state is small. We show that the minimum number of walkers required to exactly overcome the sign problem can be significantly reduced via an importance-sampling similarity transformation, even though the similarity-transformed Hamiltonian has the same stoquastic gap as the untransformed one. Furthermore, we show that in the off-half-filling Hubbard model at $U/t=8$, the real-space (site) representation has a much weaker sign problem compared to the momentum space representation. By applying importance sampling using a Gutzwiller-like guiding wavefunction, we are able to reduce the minimum number of walkers substantially in the case of $2 \times \ell$ Hubbard ladders, enabling us to get exact energies for sizeable ladders. With these results, we calculate the fundamental charge gap $\Delta E^{\mathrm{fund}}=E(N+1)+E(N-1)-2 E(N)$ for the ladder systems compared to strictly one-dimensional Hubbard chains, and show that the ladder systems have a reduced fundamental gap compared to the 1D chains.
Collapse
Affiliation(s)
- Niklas Liebermann
- Electronic Structure Theory, Max-Planck-Institute for Solid State Research, Germany
| | | | - Ali Alavi
- Max-Planck-Institute for Solid State Research, Germany
| |
Collapse
|
9
|
Iyer GR, Rubenstein BM. Finite-Size Error Cancellation in Diffusion Monte Carlo Calculations of Surface Chemistry. J Phys Chem A 2022; 126:4636-4646. [PMID: 35820033 DOI: 10.1021/acs.jpca.2c01957] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The accurate prediction of reaction mechanisms in heterogeneous (surface) catalysis is one of the central challenges in computational chemistry. Quantum Monte Carlo methods─Diffusion Monte Carlo (DMC) in particular─are being recognized as higher-accuracy, albeit more computationally expensive, alternatives to Density Functional Theory (DFT) for energy predictions of catalytic systems. A major computational bottleneck in the broader adoption of DMC for catalysis is the need to perform finite-size extrapolations by simulating increasingly large periodic cells (supercells) to eliminate many-body finite-size effects and obtain energies in the thermodynamic limit. Here, we show that it is possible to significantly reduce this computational cost by leveraging the cancellation of many-body finite-size errors that accompanies the evaluation of energy differences when calculating quantities like adsorption (binding) energies and mapping potential energy surfaces. We analyze the cancellation and convergence of many-body finite-size errors in two well-known adsorbate/slab systems, H2O/LiH(001) and CO/Pt(111). Based on this analysis, we identify strategies for obtaining binding energies in the thermodynamic limit that optimally utilize error cancellation to balance accuracy and computational efficiency. Using one such strategy, we then predict the correct order of adsorption site preference on CO/Pt(111), a challenging problem for a wide range of density functionals. Our accurate and inexpensive DMC calculations are found to unambiguously recover the top > bridge > hollow site order, in agreement with experimental observations. We proceed to use this DMC method to map the potential energy surface of CO hopping between Pt(111) adsorption sites. This reveals the existence of an L-shaped top-bridge-hollow diffusion trajectory characterized by energy barriers that provide an additional kinetic justification for experimental observations of CO/Pt(111) adsorption. Overall, this work demonstrates that it is routinely possible to achieve order-of-magnitude speedups and memory savings in DMC calculations by taking advantage of error cancellation in the calculation of energy differences that are ubiquitous in heterogeneous catalysis and surface chemistry more broadly.
Collapse
Affiliation(s)
- Gopal R Iyer
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Brenda M Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
10
|
Rudshteyn B, Weber JL, Coskun D, Devlaminck PA, Zhang S, Reichman DR, Shee J, Friesner RA. Calculation of Metallocene Ionization Potentials via Auxiliary Field Quantum Monte Carlo: Toward Benchmark Quantum Chemistry for Transition Metals. J Chem Theory Comput 2022; 18:2845-2862. [PMID: 35377642 PMCID: PMC9123894 DOI: 10.1021/acs.jctc.1c01071] [Citation(s) in RCA: 15] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
The accurate ab initio prediction of ionization energies is essential to understanding the electrochemistry of transition metal complexes in both materials science and biological applications. However, such predictions have been complicated by the scarcity of gas phase experimental data, the relatively large size of the relevant molecules, and the presence of strong electron correlation effects. In this work, we apply all-electron phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) utilizing multideterminant trial wave functions to six metallocene complexes to compare the computed adiabatic and vertical ionization energies with experimental results. We find that ph-AFQMC yields mean absolute errors (MAEs) of 1.69 ± 1.02 kcal/mol for the adiabatic energies and 2.85 ± 1.13 kcal/mol for the vertical energies. We also carry out density functional theory (DFT) calculations using a variety of functionals, which yields MAEs of 3.62-6.98 kcal/mol and 3.31-9.88 kcal/mol, as well as one variant of localized coupled cluster calculations (DLPNO-CCSD(T0) with moderate PNO cutoffs), which has MAEs of 4.96 and 6.08 kcal/mol, respectively. We also test the reliability of DLPNO-CCSD(T0) and DFT on acetylacetonate (acac) complexes for adiabatic energies measured in the same manner experimentally, and we find higher MAEs, ranging from 4.56 to 10.99 kcal/mol (with a different ordering) for DFT and 6.97 kcal/mol for DLPNO-CCSD(T0). Finally, by utilizing experimental solvation energies, we show that accurate reduction potentials in solution for the metallocene series can be obtained from the AFQMC gas phase results.
Collapse
Affiliation(s)
- Benjamin Rudshteyn
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Dilek Coskun
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Pierre A Devlaminck
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 Fifth Avenue, New York, New York 10010, United States
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
11
|
Weber JL, Vuong H, Devlaminck PA, Shee J, Lee J, Reichman DR, Friesner RA. A Localized-Orbital Energy Evaluation for Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2022; 18:3447-3459. [PMID: 35507769 DOI: 10.1021/acs.jctc.2c00111] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) has recently emerged as a promising method for the production of benchmark-level simulations of medium- to large-sized molecules because of its accuracy and favorable polynomial scaling with system size. Unfortunately, the memory footprints of standard energy evaluation algorithms are nontrivial, which can significantly impact timings on graphical processing units (GPUs) where memory is limited. Previous attempts to reduce scaling by taking advantage of the low-rank structure of the Coulombic integrals have been successful but exhibit high prefactors, making their utility limited to very large systems. Here we present a complementary cubic-scaling route to reduce memory and computational scaling based on the low rank of the Coulombic interactions between localized orbitals, focusing on the application to ph-AFQMC. We show that the error due to this approximation, which we term localized-orbital AFQMC (LO-AFQMC), is systematic and controllable via a single variable and that the method is computationally favorable even for small systems. We present results demonstrating robust retention of accuracy versus both experiment and full ph-AFQMC for a variety of test cases chosen for their potential difficulty for localized-orbital-based methods, including the singlet-triplet gaps of the polyacenes benzene through pentacene, the heats of formation for a set of Platonic hydrocarbon cages, and the total energy of ferrocene, Fe(Cp)2. Finally, we reproduce our previous result for the gas-phase ionization energy of Ni(Cp)2, agreeing with full ph-AFQMC to within statistical error while using less than 1/15th of the computer time.
Collapse
Affiliation(s)
- John L Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Hung Vuong
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Pierre A Devlaminck
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Joonho Lee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - David R Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
12
|
Shen T, Liu Y, Yu Y, Rubenstein BM. Finite temperature auxiliary field quantum Monte Carlo in the canonical ensemble. J Chem Phys 2020; 153:204108. [DOI: 10.1063/5.0026606] [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)
- Tong Shen
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, USA
| | - Yuan Liu
- Center for Ultracold Atoms, Research Laboratory of Electronics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Yang Yu
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, USA
| | | |
Collapse
|
13
|
Weber JL, Churchill EM, Jockusch S, Arthur EJ, Pun AB, Zhang S, Friesner RA, Campos LM, Reichman DR, Shee J. In silico prediction of annihilators for triplet-triplet annihilation upconversion via auxiliary-field quantum Monte Carlo. Chem Sci 2020; 12:1068-1079. [PMID: 34163873 PMCID: PMC8179011 DOI: 10.1039/d0sc03381b] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2020] [Accepted: 11/09/2020] [Indexed: 12/31/2022] Open
Abstract
The energy of the lowest-lying triplet state (T1) relative to the ground and first-excited singlet states (S0, S1) plays a critical role in optical multiexcitonic processes of organic chromophores. Focusing on triplet-triplet annihilation (TTA) upconversion, the S0 to T1 energy gap, known as the triplet energy, is difficult to measure experimentally for most molecules of interest. Ab initio predictions can provide a useful alternative, however low-scaling electronic structure methods such as the Kohn-Sham and time-dependent variants of Density Functional Theory (DFT) rely heavily on the fraction of exact exchange chosen for a given functional, and tend to be unreliable when strong electronic correlation is present. Here, we use auxiliary-field quantum Monte Carlo (AFQMC), a scalable electronic structure method capable of accurately describing even strongly correlated molecules, to predict the triplet energies for a series of candidate annihilators for TTA upconversion, including 9,10 substituted anthracenes and substituted benzothiadiazole (BTD) and benzoselenodiazole (BSeD) compounds. We compare our results to predictions from a number of commonly used DFT functionals, as well as DLPNO-CCSD(T0), a localized approximation to coupled cluster with singles, doubles, and perturbative triples. Together with S1 estimates from absorption/emission spectra, which are well-reproduced by TD-DFT calculations employing the range-corrected hybrid functional CAM-B3LYP, we provide predictions regarding the thermodynamic feasibility of upconversion by requiring (a) the measured T1 of the sensitizer exceeds that of the calculated T1 of the candidate annihilator, and (b) twice the T1 of the annihilator exceeds its S1 energetic value. We demonstrate a successful example of in silico discovery of a novel annihilator, phenyl-substituted BTD, and present experimental validation via low temperature phosphorescence and the presence of upconverted blue light emission when coupled to a platinum octaethylporphyrin (PtOEP) sensitizer. The BTD framework thus represents a new class of annihilators for TTA upconversion. Its chemical functionalization, guided by the computational tools utilized herein, provides a promising route towards high energy (violet to near-UV) emission.
Collapse
Affiliation(s)
- John L Weber
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - Emily M Churchill
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - Steffen Jockusch
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - Evan J Arthur
- Schrodinger Inc 120 West 45th Street New York NY 1003 USA
| | - Andrew B Pun
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute 162 5th Avenue New York NY 10010 USA
- Department of Physics, College of William and Mary Williamsburg VA 23187 USA
| | - Richard A Friesner
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - Luis M Campos
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - David R Reichman
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| | - James Shee
- Department of Chemistry, Columbia University 3000 Broadway New York NY 10027 USA
| |
Collapse
|
14
|
Malone FD, Zhang S, Morales MA. Accelerating Auxiliary-Field Quantum Monte Carlo Simulations of Solids with Graphical Processing Units. J Chem Theory Comput 2020; 16:4286-4297. [DOI: 10.1021/acs.jctc.0c00262] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Shuai Zhang
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| | - Miguel A. Morales
- Quantum Simulations Group, Lawrence Livermore National Laboratory, Livermore, California 94550, United States
| |
Collapse
|
15
|
Rudshteyn B, Coskun D, Weber JL, Arthur EJ, Zhang S, Reichman DR, Friesner RA, Shee J. Predicting Ligand-Dissociation Energies of 3d Coordination Complexes with Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2020; 16:3041-3054. [DOI: 10.1021/acs.jctc.0c00070] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Affiliation(s)
- Benjamin Rudshteyn
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Dilek Coskun
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - John L. Weber
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Evan J. Arthur
- Schrodinger Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, 162 5th Avenue, New York, New York 10010, United States
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187, United States
| | - David R. Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - James Shee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
16
|
Lee J, Malone FD, Morales MA. An auxiliary-Field quantum Monte Carlo perspective on the ground state of the dense uniform electron gas: An investigation with Hartree-Fock trial wavefunctions. J Chem Phys 2019. [DOI: 10.1063/1.5109572] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Joonho Lee
- College of Chemistry, University of California, Berkeley, California 94720, USA
| | - Fionn D. Malone
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| | - Miguel A. Morales
- Quantum Simulations Group, Lawrence Livermore National Laboratory, 7000 East Avenue, Livermore, California 94551, USA
| |
Collapse
|
17
|
Motta M, Shee J, Zhang S, Chan GKL. Efficient Ab Initio Auxiliary-Field Quantum Monte Carlo Calculations in Gaussian Bases via Low-Rank Tensor Decomposition. J Chem Theory Comput 2019; 15:3510-3521. [DOI: 10.1021/acs.jctc.8b00996] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Mario Motta
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - James Shee
- Department of Chemistry, Columbia University, New York, New York 10027, United States
| | - Shiwei Zhang
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| |
Collapse
|
18
|
On Achieving High Accuracy in Quantum Chemical Calculations of 3d Transition Metal-Containing Systems: A Comparison of Auxiliary-Field Quantum Monte Carlo with Coupled Cluster, Density Functional Theory, and Experiment for Diatomic Molecules. J Chem Theory Comput 2019; 15:2346-2358. [DOI: 10.1021/acs.jctc.9b00083] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022]
|
19
|
Lee J, Small DW, Head-Gordon M. Open-shell coupled-cluster valence-bond theory augmented with an independent amplitude approximation for three-pair correlations: Application to a model oxygen-evolving complex and single molecular magnet. J Chem Phys 2018; 149:244121. [DOI: 10.1063/1.5052667] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - David W. Small
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
20
|
Hao H, Shee J, Upadhyay S, Ataca C, Jordan KD, Rubenstein BM. Accurate Predictions of Electron Binding Energies of Dipole-Bound Anions via Quantum Monte Carlo Methods. J Phys Chem Lett 2018; 9:6185-6190. [PMID: 30299101 DOI: 10.1021/acs.jpclett.8b02733] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Neutral molecules with sufficiently large dipole moments can bind electrons in diffuse nonvalence orbitals with most of their charge density far from the nuclei, forming so-called dipole-bound anions. Because long-range correlation effects play an important role in the binding of an excess electron and overall binding energies are often only on the order of 10s-100s of wave numbers, predictively modeling dipole-bound anions remains a challenge. Here, we demonstrate that quantum Monte Carlo methods can accurately characterize molecular dipole-bound anions with near-threshold dipole moments. We also show that correlated sampling Auxiliary Field Quantum Monte Carlo is particularly well-suited for resolving the fine energy differences between the neutral and anionic species. These results shed light on the fundamental limitations of quantum Monte Carlo methods and pave the way toward using them for the study of weakly bound species that are too large to model using traditional electron structure methods.
Collapse
Affiliation(s)
- Hongxia Hao
- Department of Chemistry , Brown University , Providence , Rhode Island 02912 , United States
| | - James Shee
- Department of Chemistry , Columbia University , New York , New York 10027 , United States
| | - Shiv Upadhyay
- Department of Chemistry , University of Pittsburgh , Pittsburgh , Pennsylvania 15260 , United States
| | - Can Ataca
- Department of Physics , University of Maryland-Baltimore County , Baltimore , Maryland 21250 , United States
| | - Kenneth D Jordan
- Department of Chemistry , University of Pittsburgh , Pittsburgh , Pennsylvania 15260 , United States
| | - Brenda M Rubenstein
- Department of Chemistry , Brown University , Providence , Rhode Island 02912 , United States
| |
Collapse
|
21
|
Tiihonen J, Kylänpää I, Rantala TT. Computation of Dynamic Polarizabilities and van der Waals Coefficients from Path-Integral Monte Carlo. J Chem Theory Comput 2018; 14:5750-5763. [DOI: 10.1021/acs.jctc.8b00859] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Juha Tiihonen
- Laboratory of Physics, Tampere University of Technology, P. O. Box 692, FI-33101 Tampere, Finland
| | - Ilkka Kylänpää
- Laboratory of Physics, Tampere University of Technology, P. O. Box 692, FI-33101 Tampere, Finland
| | - Tapio T. Rantala
- Laboratory of Physics, Tampere University of Technology, P. O. Box 692, FI-33101 Tampere, Finland
| |
Collapse
|
22
|
Shee J, Arthur EJ, Zhang S, Reichman DR, Friesner RA. Phaseless Auxiliary-Field Quantum Monte Carlo on Graphical Processing Units. J Chem Theory Comput 2018; 14:4109-4121. [DOI: 10.1021/acs.jctc.8b00342] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- James Shee
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Evan J. Arthur
- Schrödinger
Inc., 120 West 45th Street, New York, New York 10036, United States
| | - Shiwei Zhang
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - David R. Reichman
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| | - Richard A. Friesner
- Department of Chemistry, Columbia University, 3000 Broadway, New York, New York 10027, United States
| |
Collapse
|
23
|
Ab initio computations of molecular systems by the auxiliary‐field quantum Monte Carlo method. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1364] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
24
|
Motta M, Zhang S. Communication: Calculation of interatomic forces and optimization of molecular geometry with auxiliary-field quantum Monte Carlo. J Chem Phys 2018; 148:181101. [DOI: 10.1063/1.5029508] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Mario Motta
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Shiwei Zhang
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187-8795, USA
| |
Collapse
|
25
|
Motta M, Zhang S. Computation of Ground-State Properties in Molecular Systems: Back-Propagation with Auxiliary-Field Quantum Monte Carlo. J Chem Theory Comput 2017; 13:5367-5378. [DOI: 10.1021/acs.jctc.7b00730] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Mario Motta
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| | - Shiwei Zhang
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187-8795, United States
| |
Collapse
|
26
|
Takeshita TY, de Jong WA, Neuhauser D, Baer R, Rabani E. Stochastic Formulation of the Resolution of Identity: Application to Second Order Møller–Plesset Perturbation Theory. J Chem Theory Comput 2017; 13:4605-4610. [DOI: 10.1021/acs.jctc.7b00343] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Tyler Y. Takeshita
- Department
of Chemistry, University of California Berkeley, Berkeley, California 94720, United States
- Materials
Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Wibe A. de Jong
- Computational
Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Daniel Neuhauser
- Department
of Chemistry and Biochemistry, University of California, Los Angeles, California 90095, United States
| | - Roi Baer
- Fritz
Harber Center for Molecular Dynamics, Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Eran Rabani
- Department
of Chemistry, University of California Berkeley, Berkeley, California 94720, United States
- Materials
Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
- The Sackler
Center for Computational Molecular Science, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|