Applying Quantum Monte Carlo to the Electronic Structure
Applying Quantum Monte Carlo to the Electronic Structure Problem: Potential Energy Curves of N2 and CO with Diffusion Monte Carlo Andrew D Powell, Richard Dawes Department of Chemistry, Missouri University of Science and Technology, Rolla, MO, USA. 1 Outline Motivation Scaling and Accuracy Introduction Background of QMC methodology Results
Comparison with spectroscopically accurate PECs Cost analysis Conclusion 2 Motivation Scaling Quantum Monte Carlo (QMC) methods are efficiently parallelized with a large number of cores n3 with the number of electrons as n7 or worse with traditional quantum methods With its cost prefactor, how practical is it for routine application? Accuracy Capable of capturing large fractions of dynamic correlation energy How does it compare with MRCI or CCSD(T) for small molecules? Can QMC produce all data for global PESs?
3 Introduction Quantum Monte Carlo (QMC) is an alternative method to solve the Schrdinger equation. CASINO1 QMCPACK2 For traditional quantum chemistry methods, there are many codes available. MOLPRO Gaussian
GAMESS CFOUR R.J. Needs, M.D. Towler, N.D. Drummond and P. Lpez Ros, J. Phys.: Condensed Matter 22, 023201 (2010). 2 J. Kim, K.P. Esler, J. McMinis, M.A. Morales, B.K. Clark, L. Shulenburger, and D.M. Ceperley, J. Phys.: Conf. Ser. 402, 012008 (2012). 1 4 Introduction Test systems: diatomic PECs for N2 and CO CASINO code Background of QMC Variational Monte Carlo (VMC) and Diffusion Monte Carlo (DMC)
5 Background of VMC VMC is designed to sample a trial wave function and to calculate the expectation value of the Hamiltonian using Monte Carlo numerical integration. VMC is primarily used to optimize parameters for dynamic electron correlation. 6 Background of DMC The DMC method attempts to simultaneously create and sample the unknown exact ground state wave function.
Constrained from above the exact solution by the fixed-node approximation Fixed-node approximation: the nodal surface of the wave function is constrained to be that of the trial wave function. (1) (2) (3) Due to the statistical nature of the approach, there is uncertainty associated with a calculated energy. E E 7 Preparation of trial wave function A trial wave function is used as an initial reference for the method. It can be prepared by methods such as DFT, HF, CASSCF, etc. We use MCSCF trial wave functions from GAMESS.3
Molecular orbitals and configuration coefficients are prepared by scripts for use in the CASINO program. aug-cc-pwCV5Z* basis set High angular momentum functions (l f) were removed 3 M.W.Schmidt, K.K.Baldridge, J.A.Boatz, S.T.Elbert, M.S.Gordon, J.H.Jensen, S.Koseki N.Matsunaga, K.A.Nguyen, S.J.Su, T.L.Windus, M.Dupuis, J.A.Montgomery 8 Trial Wave Function Multi-determinant expansion is used to describe static electron correlation Includes a Jastrow factor to capture dynamic correlation
Makes the trial wave function depend explicitly on the interparticle separations Maintains nodal structure of the trial wave function. Enforces the cusp conditions. Tightness of VMC optimization Procedure 1: 20,000 CPU-hrs. Procedure 2: 2,000 CPU-hrs. 9 Single-reference breakdown Divergence of CCSD(T) occurs for rNN > 2 .4 4 X. Li and J. Paldus, J. Chem. Phys. 129, 054104 (2008) 10
N(4S) + N(4S) N2 (X1g+) ~80 cm-1 uncertainties after 105 samples (1500-2000 cpu-hrs for DMC) 11 A DW-MRCI benchmark PEC for CO 1-state MRCI calculation produces discontinuity 9-states are degenerate asymptotically 12 A DW-MRCI benchmark PEC for CO DW-MRCI/CBS with 9 states SO- and SR-Coupling
13 Accuracy of benchmark PEC for CO CO J=0 Vibrational Levels v Calculated Experiment Error 0 1 2 3 4
0.18 0.60 1.17 1.88 14 C(3P) + O(3P) CO(X1+) ~80 cm-1 uncertainties after 105 samples (1500-2000 cpu-hrs for DMC) 15 Conclusions Cost of QMC Tightening of the VMC optimization gave no benefit in the subsequent DMC. Not competitive with traditional chemistry (N2 and CO take ~5 min per point with MRCI/CBS. Unlikely to be used to generate all points of a PES. n3 scaling ensures crossover point.
Accuracy Comparable to the best ab initio or empirical PECs. Future Arbiter in difficult cases? Testing reaction barriers for F + H2O, O + O2, and O2 + propargyl. Using QMCPACK (more efficient for MCSCF trial wave functions, capable of open-shell systems, allows spdfg functions) Testing FCIQMC in same applications 16 17 Dynamically-weighted state-averaged multireference electronic structure theory E SA MCSCF
n w i EiMCSCF i 0 W n W wi i 0 HCO Dynamic weighting1 2 wi cosh (Ei 0 / )
M.P. Deskevich, D.J. Nesbitt and H-J. Werner, JCP 120 7281 (2004) 6 R. Dawes, A.W. JasperS.A. Reid et al. J. Phys. Chem. Lett. 1 641 (2010) 5 Importance of Dynamic Weighting: BeOBe BeO + Be 1 state 23 states equal weights In single state calc (at left), avoided crossing with higher lying state (blue oval, at right) causes discontinuity. In fixed weights calc (at right) crossing with higher lying state (red
oval) causes discontinuity cascade affecting ground state. Six lowest A1 states (co-linear asymmetric stretch C2V) Importance of Dynamic Weighting: BeOBe BeO + Be 23 states equal weights 23 states dynamic weights With dynamic weights (at right), crossing still occurs (blue oval), but due to negligible weight, no disruption results.
Selects between two and only another is OK Flat Rayleigh fading channel -overlap -windows -equal Pfa: fair comparison !! Code length 64 M=64, R varies M=window size R=hop size Pfa set by simulations for windowing & overlap cases Windowing causes...
A "New World" One Italian sailor, Christopher Columbus, believed that he could reach Asia by sailing west across the Atlantic. The king and queen of Spain gave Columbus the money to make his journey. In October 1492, Columbus landed on...
Arial Calibri Harlow Solid Italic Haettenschweiler Default Design Telling Time Slide 2 Get the answer in your head (or on paper) before clicking for it: You know how we say at 2:10- "It's 10 after 2:00," or at 4:20 "It's...
By capturing bats in the summer and putting small transmitters on their backs, we have been able to track them to caves that we did not know existed. And, with the emergence of WNS, we have also discovered new caves...
The Story of an Hour Kate Chopin Main characters Mrs. Mallard Wife to Brently Mallard Main character in the story Josephine Friend to Mrs. Mallard Person who tells of the "disaster" Richards Friend of Brently, who shares in the retelling...
The Universal Declaration of Human Rights of 1948 is a landmark in the history of the United Nation's which upholds the dignity of every human beings in this earth. Religious freedom is an inherited right of every human being and...
Ready to download the document? Go ahead and hit continue!