Size-extensive wave functions for QMC A linear-scaling GVB approach

Size: px
Start display at page:

Download "Size-extensive wave functions for QMC A linear-scaling GVB approach"

Transcription

1 Size-extensive wave functions for QMC A linear-scaling GVB approach Claudia Filippi, University of Twente, The Netherlands Francesco Fracchia, University of Pisa, Italy Claudio Amovilli, University of Pisa, Italy ES12 Workshop, Wake Forest University, June 5-8, 2012

2 Wave functions and QMC: A healthy obsession Advantages of QMC Favorable scaling: N 4 versus for instance N 7 in CCSD(T) Freedom on the functional form of the wave functions... Inclusion of r ij terms, non-orthogonal orbitals etc. However No optimal recipe for the construction of the trial function Multi-determinant wave function is often solution Computational cost remains high even for small systems

3 In search of different functional forms Active area of research Geminals, pfaffians, backflow... The trial wavefunction should meet the following requirements Accurate Compact functional form Good scaling with respect to the size of the system Size-consistency and size-extensivity

4 Size-extensivity and size-consistency Size-extensivity (Bartlett, Ann. Rev. Phys. Chem. 1981) Computed energy scales linearly with number of particles lim (E n+1 E n ) = C n Size-consistency: Additivity of separate chemical units E A+B = E A + E B

5 Common choice in QMC Jastrow-Slater wave functions Ψ(r 1,..., r N ) = J (r 1,..., r N ) i c i D i (r 1,..., r N ) J Jastrow correlation factor Positive function of inter-particle distances Explicit dependence on electron-electron distances r ij ci D i Determinants of single-particle orbitals Few and not millions of determinants Determines the nodal surface

6 Determinants: A complete active space expansion CAS wave functions Good for isomerization, excitations,...!"#$%&' ()*"&' +++' In principle, size-consistent and size-extensive but scales as N! Small active space Not all electrons correlated Small portion of PES described Large active space Truncation Not size-extensive/consistent... dangerous

7 CAS wave functions and atomization energies G2 set with full-valence CAS (Petruzielo et al. JCP, 2012) Nearly chemical accuracy on atomization energies MAD of 1.2 kcal/mol in DMC Very careful truncation (small thereshold + inclusion of doubles) However, C 2 H determinants!

8 CAS wave functions and excitations Cyanines dyes (Send, Valsson, Filippi, JCTC 2011) Starting determinantal wave function: Already compromises CAS (4,6) (6,10) (8,14) (10,9) (12,11) Det QMC setup: Truncation Thr Det

9 A localized description of a many-body wave function Desire From molecular orbital to local description of molecule Bonding (over 2 or more centers), antibonding, lone-pairs...! "#$%&'#(%#)!!!! *$+,! &'#(%#)!!!!! In the context of a multi-determinant wave function!

10 GVB wave functions as inspiration Inspiration from multi-determinantal expression of GVB-PP Ψ GVB PP = Φ N i=1 Limited variational flexibility Our starting point q i Φ a i a i b i b i + Ψ LGVB1 = c 0 Φ N i<j 1 2 N i=1 q i q j Φ a i a i,a j a j b i b i,b j b j +... c i Φ a i a i b i b i

11 Our starting wave function: Emerging of localization in MCSCF LGVB1 wave function Ψ LGVB1 = c 0 Φ N i=1 c i Φ a i a i b i b i

12 Our starting wave function: Emerging of localization LGVB1 wave function Ψ LGVB1 = c 0 Φ 0 + Correlated all valence electrons 1 2 N i=1 c i Φ b i b i a i a i Good description of single-bond dissociation in a molecule Compact and scales linearly with system size Localization emerges naturally in MCSCF optimization

13 Improving on our starting LGVB1 wave function Tests: Use these localized orbitals in CISDTQ 1) Correlation mainly among close pairs (one atom in common)... Blocks of CAS(4,4) wave functions

14 Improving on our starting LGVB1 wave function Tests: Use these localized orbitals in CISDTQ 2) Hierarchy of excitations within close pair, completing CAS(4,4) E1 E2

15 Hierarchy of excitations: Correlating adjacent electron pairs E3 E4 E5 E6

16 Hierarchy of excitations: Correlating adjacent electron pairs E7 E8 E9 E10

17 Modular truncation Only excitations that couple adjacent electron pairs Approach is size-extensive Number of determinants scales linearly with size system Additional hierarchy within excitations of each pair

18 Building the LGVB wave function Correlating all 18 adjacent electron pairs

19 Building the LGVB wave function CAS(4,4) Correlating all 18 adjacent electron pairs

20 Building the LGVB wave function CAS(4,4) Correlating all 18 adjacent electron pairs

21 Building the LGVB wave function CAS(4,4) Correlating all 18 adjacent electron pairs

22 Building the LGVB wave function CAS(4,4) Correlating all 18 adjacent electron pairs

23 LGVBn wave functions: Correlating adjecent electron pairs with n carbon atoms Ψ CSFs N(CSF) LGVB1 E0 + E1 3n + 2 LGVB2 LGVB1 + E2 15n + 2 LGVB3 LGVB2 + E3 18n + 3 LGVB4 LGVB3 + E4 30n + 3 LGVB5 LGVB4 + E5 36n + 3 LGVB6 LGVB5 + E6 48n + 3 LGVB7 LGVB6 + E7 60n + 3 LGVB8 LGVB7 + E8 72n + 3 LGVB9 LGVB8 + E9 84n + 3 LGVB10 LGVB9 + E10 96n + 3

24 Size of LGVBn wave functions Ψ N(CSF) LGVB1 11 LGVB2 47 LGVB4 93 LGVB CAS(20,20) >10 8

25 Back to quantum Monte Carlo: Jastrow-LGVBn Ψ(r 1,..., r N ) = J (r 1,..., r N ) LGVBn Generation of localized orbitals in MCSCF with LGVB1 Optimization of orbitals, CI, Jastrow in J-LGVB1 Further increase the wave function as J-LGVBn (n> 1)

26 Fragmentation tests Homolytic bond breaking Breaking of one electron pair Difficult to balance description of correlation in various fragments We must respect size-consistency and size-extensivity Set of fragmentation reactions: N 2 H 4 2NH 2 HNO 2 OH + NO CH 3 NH 2 CH 3 + NH 2 CH 3 OH CH 3 + OH

27 N 2 H 4 : VMC total energy VMC energy (Hartree) J-LGVB1 J-LGVB J-LGVBn order

28 N 2 H 4 : VMC total energy det VMC energy (Hartree) J-LGVB1 J-LGVB J-LGVBn order

29 N 2 H 4 : VMC total energy VMC energy (Hartree) J-LGVB1 1 det CAS(2,2) J-LGVB J-LGVBn order

30 N 2 H 4 : VMC total energy VMC energy (Hartree) J-LGVB1 1 det CAS(2,2) CAS(6,6) J-LGVB J-LGVBn order

31 N 2 H 4 : VMC total energy VMC energy (Hartree) J-LGVB1 1 det CAS(2,2) CAS(6,6) CAS(14,14) J-LGVB J-LGVBn order

32 N 2 H 4 : VMC total energy VMC energy (Hartree) J-LGVB det CAS(2,2) CAS(6,6) CAS(14,14) J-LGVBn order J-LGVB10

33 DMC fragmentation energy Exp. DMC binding energy (kcal/mol) J-LGVB1 J-LGVB J-LGVBn order

34 DMC fragmentation energy + DMC binding energy (kcal/mol) Exp. J-LGVB1 1 det J-LGVB J-LGVBn order

35 DMC fragmentation energy + DMC binding energy (kcal/mol) Exp. J-LGVB1 CAS(2,2) 1 det J-LGVB J-LGVBn order

36 DMC fragmentation energy Exp. DMC binding energy (kcal/mol) J-LGVB1 CAS(2,2) CAS(6,6) 1 det J-LGVB J-LGVBn order

37 DMC fragmentation energy + DMC binding energy (kcal/mol) Exp. J-LGVB1 CAS(14,14) CAS(2,2) CAS(6,6) 1 det J-LGVB J-LGVBn order

38 Fragmentation energy: DMC and CC at comparison Binding energy (kcal/mol) J-LGVB1 Exp. J-LGVB J-LGVBn order

39 Fragmentation energy: DMC and CC at comparison Binding energy (kcal/mol) J-LGVB1 Exp. J-LGVB CCSD(T)/cc-pVTZ J-LGVBn order

40 Fragmentation energy: DMC and CC at comparison Binding energy (kcal/mol) J-LGVB1 Exp. CCSD(T)/cc-pVQZ J-LGVB CCSD(T)/cc-pVTZ J-LGVBn order

41 Fragmentation energy: DMC and CC at comparison Binding energy (kcal/mol) J-LGVB1 Exp. CCSD(T)/cc-pV5Z CCSD(T)/cc-pVQZ J-LGVB CCSD(T)/cc-pVTZ J-LGVBn order

42 Fragmentations: Summary Method MAD MAX (kcal/mol) DMC 1 det DMC J-LGVB DMC J-LGVB CCSD(T)/cc-pVTZ CCSD(T)/cc-pVQZ CCSD(T)/cc-pV5Z B3LYP PBE M06-2X M08-HX B2PLYP DMC/pVTZ; DFT/aug-cc-pVQZ

43 Barrier heigths in reactions: In progress (1) Example: CH 4 + NH 2 CH 3 + NH 3 New orbital type: Bonding/antibonding of 2 elec on 3 centers Transition state Method BH for BH rev E DMC a 1 det J-LGVB J-LGVB G CCSD(T) a Reference value ? a DMC/aug-pVTZ; CC/aug-cc-pVTZ

44 Barrier heigths in reactions: In progress (2) Example: NH 2 + H 2 O NH 3 + OH Again: Bonding/antibonding of 2 elec on 3 centers Transition state Method BH for BH rev E DMC a 1 det J-LGVB J-LGVB G CCSD(T) a Reference value ??? a DMC/aug-pVTZ; CC/aug-cc-pVTZ

45 Challenges: The classical case of benzene 1) Multiple resonances J-(M)LGBV wave functions? antibonding bonding Problematic for MCSCF but not for QMC

46 Challenges: The classical case of benzene 2) Hybrid localized-delocalized description Localized description would be preferable

47 Conclusions Development of Jastrow-Slater Linear-GVB wave functions based on localized orbital, allow modular truncation correlate all valence electrons, compact, and size-extensive number of determinants scales linearly with system size Excellent performance on fragmentation of homolytic bonds J-LGVB2 satisfies chemical accuracy Performs as well as CCSD(T) with large cc-pv5z Barrier heights in progress Open challenge: Multiple resonances

48 Collaborators Francesco Fracchia, University of Pisa, Italy Claudio Amovilli, University of Pisa, Italy Fracchia, Filippi, Amovilli, JCTC (2012)

Recent advances in quantum Monte Carlo for quantum chemistry: optimization of wave functions and calculation of observables

Recent advances in quantum Monte Carlo for quantum chemistry: optimization of wave functions and calculation of observables Recent advances in quantum Monte Carlo for quantum chemistry: optimization of wave functions and calculation of observables Julien Toulouse 1, Cyrus J. Umrigar 2, Roland Assaraf 1 1 Laboratoire de Chimie

More information

Quantum Monte Carlo wave functions and their optimization for quantum chemistry

Quantum Monte Carlo wave functions and their optimization for quantum chemistry Quantum Monte Carlo wave functions and their optimization for quantum chemistry Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France CEA Saclay, SPhN Orme des Merisiers April 2015 Outline

More information

QMC dissociation energy of the water dimer: Time step errors and backflow calculations

QMC dissociation energy of the water dimer: Time step errors and backflow calculations QMC dissociation energy of the water dimer: Time step errors and backflow calculations Idoia G. de Gurtubay and Richard J. Needs TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay.

More information

Optimization of quantum Monte Carlo (QMC) wave functions by energy minimization

Optimization of quantum Monte Carlo (QMC) wave functions by energy minimization Optimization of quantum Monte Carlo (QMC) wave functions by energy minimization Julien Toulouse, Cyrus Umrigar, Roland Assaraf Cornell Theory Center, Cornell University, Ithaca, New York, USA. Laboratoire

More information

Ab initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier

Ab initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier Ab initio calculations for potential energy surfaces D. Talbi GRAAL- Montpellier A theoretical study of a reaction is a two step process I-Electronic calculations : techniques of quantum chemistry potential

More information

Quantum Monte Carlo backflow calculations of benzene dimers

Quantum Monte Carlo backflow calculations of benzene dimers Quantum Monte Carlo backflow calculations of benzene dimers Kathleen Schwarz*, Cyrus Umrigar**, Richard Hennig*** *Cornell University Department of Chemistry, **Cornell University Department of Physics,

More information

Optimization of quantum Monte Carlo wave functions by energy minimization

Optimization of quantum Monte Carlo wave functions by energy minimization Optimization of quantum Monte Carlo wave functions by energy minimization Julien Toulouse, Roland Assaraf, Cyrus J. Umrigar Laboratoire de Chimie Théorique, Université Pierre et Marie Curie and CNRS, Paris,

More information

Quantum Monte Carlo for excited state calculations

Quantum Monte Carlo for excited state calculations Quantum Monte Carlo for excited state calculations Claudia Filippi MESA+ Institute for Nanotechnology, Universiteit Twente, The Netherlands Winter School in Theoretical Chemistry, Helsinki, Finland, Dec

More information

Fixed-Node quantum Monte Carlo for Chemistry

Fixed-Node quantum Monte Carlo for Chemistry Fixed-Node quantum Monte Carlo for Chemistry Michel Caffarel Lab. Physique et Chimie Quantiques, CNRS-IRSAMC, Université de Toulouse e-mail : caffarel@irsamc.ups-tlse.fr. p.1/29 The N-body problem of Chemistry

More information

QMC dissociation energies of three-electron hemibonded radical cation dimers... and water clusters

QMC dissociation energies of three-electron hemibonded radical cation dimers... and water clusters QMC dissociation energies of three-electron hemibonded radical cation dimers... and water clusters Idoia G. de Gurtubay TCM group. Cavendish Laboratory University of Cambridge Idoia G. de Gurtubay. Quantum

More information

arxiv: v2 [physics.chem-ph] 11 Mar 2011

arxiv: v2 [physics.chem-ph] 11 Mar 2011 Quantum Monte Carlo with Jastrow-Valence-Bond wave functions Benoît Braïda 1, Julien Toulouse 1, Michel Caffarel 2, and C. J. Umrigar 3 1 Laboratoire de Chimie Théorique, Université Pierre et Marie Curie

More information

Wave Function Optimization and VMC

Wave Function Optimization and VMC Wave Function Optimization and VMC Jeremy McMinis The University of Illinois July 26, 2012 Outline Motivation History of Wave Function Optimization Optimization in QMCPACK Multideterminant Example Motivation:

More information

Ab-initio molecular dynamics for High pressure Hydrogen

Ab-initio molecular dynamics for High pressure Hydrogen Ab-initio molecular dynamics for High pressure Hydrogen Claudio Attaccalite Institut d'electronique, Microélectronique et Nanotechnologie (IEMN), Lille Outline A brief introduction to Quantum Monte Carlo

More information

Orbital-dependent backflow transformations in quantum Monte Carlo

Orbital-dependent backflow transformations in quantum Monte Carlo transformations in quantum Monte Carlo P. Seth, P. López Ríos, and R. J. Needs TCM group, Cavendish Laboratory, University of Cambridge 5 December 2012 VMC and DMC Optimization Wave functions Variational

More information

Variational Monte Carlo Optimization and Excited States

Variational Monte Carlo Optimization and Excited States Variational Monte Carlo Optimization and Excited States Eric Neuscamman August 9, 2018 motivation charge transfer core spectroscopy double excitations the menu aperitif: number counting Jastrows main course:

More information

Introduction to Computational Chemistry

Introduction to Computational Chemistry Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry Chemicum 4th floor vesa.hanninen@helsinki.fi September 10, 2013 Lecture 3. Electron correlation methods September

More information

Quantum Monte Carlo methods

Quantum Monte Carlo methods Quantum Monte Carlo methods Lubos Mitas North Carolina State University Urbana, August 2006 Lubos_Mitas@ncsu.edu H= 1 2 i i 2 i, I Z I r ii i j 1 r ij E ion ion H r 1, r 2,... =E r 1, r 2,... - ground

More information

Multiconfiguration wave functions for quantum Monte Carlo calculations of first-row diatomic molecules

Multiconfiguration wave functions for quantum Monte Carlo calculations of first-row diatomic molecules Multiconfiguration wave functions for quantum Monte Carlo calculations of first-row diatomic molecules Claudia Filippi Laboratory of Atomic and Solid State Physics and Theory Center, Cornell University,

More information

A guide to the perplexed

A guide to the perplexed Wave function embedding methods and excited states: A guide to the perplexed Csaba Daday and Claudia Filippi MESA+ Institute for Nanotechnology, University of Twente, the Netherlands Carolin König and

More information

Auxiliary-field quantum Monte Carlo calculations of excited states and strongly correlated systems

Auxiliary-field quantum Monte Carlo calculations of excited states and strongly correlated systems Auxiliary-field quantum Monte Carlo calculations of excited states and strongly correlated systems Formally simple -- a framework for going beyond DFT? Random walks in non-orthogonal Slater determinant

More information

Time-dependent linear-response variational Monte Carlo.

Time-dependent linear-response variational Monte Carlo. Time-dependent linear-response variational Monte Carlo. Bastien Mussard bastien.mussard@colorado.edu https://mussard.github.io/ Julien Toulouse julien.toulouse@upmc.fr Sorbonne University, Paris (web)

More information

Anion-π and π-π cooperative interactions

Anion-π and π-π cooperative interactions Chapter 5 Anion-π and π-π cooperative interactions 5.1 Introduction The design of selective receptors of anionic species is a very active area of research within supramolecular chemistry due to the potential

More information

Electron Correlation - Methods beyond Hartree-Fock

Electron Correlation - Methods beyond Hartree-Fock Electron Correlation - Methods beyond Hartree-Fock how to approach chemical accuracy Alexander A. Auer Max-Planck-Institute for Chemical Energy Conversion, Mülheim September 4, 2014 MMER Summerschool 2014

More information

Resonating Valence Bond wave function with molecular orbitals: application to diatomic molecules

Resonating Valence Bond wave function with molecular orbitals: application to diatomic molecules Resonating Valence Bond wave function with molecular orbitals: application to diatomic molecules M. Marchi 1,2, S. Azadi 2, M. Casula 3, S. Sorella 1,2 1 DEMOCRITOS, National Simulation Center, 34014,

More information

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)

More information

4 Post-Hartree Fock Methods: MPn and Configuration Interaction

4 Post-Hartree Fock Methods: MPn and Configuration Interaction 4 Post-Hartree Fock Methods: MPn and Configuration Interaction In the limit of a complete basis, the Hartree-Fock (HF) energy in the complete basis set limit (ECBS HF ) yields an upper boundary to the

More information

Valence electronic structure of isopropyl iodide investigated by electron momentum spectroscopy. --- Influence of intramolecular interactions

Valence electronic structure of isopropyl iodide investigated by electron momentum spectroscopy. --- Influence of intramolecular interactions Valence electronic structure of isopropyl iodide investigated by electron momentum spectroscopy --- Influence of intramolecular interactions Minfu Zhao, Xu Shan, Shanshan Niu, Yaguo Tang, Zhaohui Liu,

More information

Beyond the Hartree-Fock Approximation: Configuration Interaction

Beyond the Hartree-Fock Approximation: Configuration Interaction Beyond the Hartree-Fock Approximation: Configuration Interaction The Hartree-Fock (HF) method uses a single determinant (single electronic configuration) description of the electronic wavefunction. For

More information

Performance of Hartree-Fock and Correlated Methods

Performance of Hartree-Fock and Correlated Methods Chemistry 460 Fall 2017 Dr. Jean M. Standard December 4, 2017 Performance of Hartree-Fock and Correlated Methods Hartree-Fock Methods Hartree-Fock methods generally yield optimized geomtries and molecular

More information

Electron Correlation Methods

Electron Correlation Methods Electron Correlation Methods HF method: electron-electron interaction is replaced by an average interaction E HF c = E 0 E HF E 0 exact ground state energy E HF HF energy for a given basis set HF E c

More information

Solution of the Electronic Schrödinger Equation. Using Basis Sets to Solve the Electronic Schrödinger Equation with Electron Correlation

Solution of the Electronic Schrödinger Equation. Using Basis Sets to Solve the Electronic Schrödinger Equation with Electron Correlation Solution of the Electronic Schrödinger Equation Using Basis Sets to Solve the Electronic Schrödinger Equation with Electron Correlation Errors in HF Predictions: Binding Energies D e (kcal/mol) HF Expt

More information

Density Functional Theory

Density Functional Theory Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function

More information

Development and application of statistical and quantum mechanical methods for modelling molecular ensembles

Development and application of statistical and quantum mechanical methods for modelling molecular ensembles Development and application of statistical and quantum mechanical methods for modelling molecular ensembles By Ellen T. Swann A thesis submitted for the degree of Doctor of Philosophy of the Australian

More information

Lec20 Fri 3mar17

Lec20 Fri 3mar17 564-17 Lec20 Fri 3mar17 [PDF]GAUSSIAN 09W TUTORIAL www.molcalx.com.cn/wp-content/uploads/2015/01/gaussian09w_tutorial.pdf by A Tomberg - Cited by 8 - Related articles GAUSSIAN 09W TUTORIAL. AN INTRODUCTION

More information

Properties of reactive oxygen species by quantum Monte Carlo

Properties of reactive oxygen species by quantum Monte Carlo Properties of reactive oxygen species by quantum Monte Carlo The MIT Faculty has made this article openly available. Please share how this access benefits you. Your story matters. Citation As Published

More information

Accurate barrier heights using diffusion Monte Carlo

Accurate barrier heights using diffusion Monte Carlo Accurate barrier heights using diffusion Monte Carlo Kittithat Krongchon, Brian Busemeyer, and Lucas K. Wagner Department of Physics, University of Illinois at Urbana-Champaign (Dated: March 9, 217) Fixed

More information

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah 1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations

More information

Electronic structure quantum Monte Carlo methods and variable spins: beyond fixedphase/node

Electronic structure quantum Monte Carlo methods and variable spins: beyond fixedphase/node Electronic structure quantum Monte Carlo methods and variable spins: beyond fixedphase/node approximations Cody Melton, M. Chandler Bennett, L. Mitas, with A. Ambrosetti, F. Pederiva North Carolina State

More information

Building a wavefunction within the Complete-Active. Cluster with Singles and Doubles formalism: straightforward description of quasidegeneracy

Building a wavefunction within the Complete-Active. Cluster with Singles and Doubles formalism: straightforward description of quasidegeneracy Building a wavefunction within the Complete-Active Active-Space Coupled-Cluster Cluster with Singles and Doubles formalism: straightforward description of quasidegeneracy Dmitry I. Lyakh (Karazin Kharkiv

More information

Many-body electronic structure theory: Is Quantum Monte-Carlo affordable for molecular crystals?

Many-body electronic structure theory: Is Quantum Monte-Carlo affordable for molecular crystals? Many-body electronic structure theory: Is Quantum Monte-Carlo affordable for molecular crystals? Gerit Brandenburg 23 rd of November 2017 CHEMISTRY DEPARTMENT, UNIVERSITY OF TURIN,

More information

arxiv:physics/ v2 [physics.chem-ph] 9 Apr 2005

arxiv:physics/ v2 [physics.chem-ph] 9 Apr 2005 Electronic Quantum Monte Carlo Calculations of Atomic Forces, Vibrations, and Anharmonicities arxiv:physics/0411209v2 [physics.chem-ph] 9 Apr 2005 Myung Won Lee a), Massimo Mella b), and Andrew M. Rappe

More information

OVERVIEW OF QUANTUM CHEMISTRY METHODS

OVERVIEW OF QUANTUM CHEMISTRY METHODS OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density

More information

QUANTUM CHEMISTRY FOR TRANSITION METALS

QUANTUM CHEMISTRY FOR TRANSITION METALS QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline

More information

Convergence of many-body wavefunction expansions using a plane wave basis: From the homogeneous electron gas to the solid state

Convergence of many-body wavefunction expansions using a plane wave basis: From the homogeneous electron gas to the solid state Convergence of many-body wavefunction expansions using a plane wave basis: From the homogeneous electron gas to the solid state TCM Electronic Structure Discussion Group James Shepherd (CUC3, Alavi Group)

More information

Ab-initio simulation of liquid water by quantum Monte Carlo

Ab-initio simulation of liquid water by quantum Monte Carlo Ab-initio simulation of liquid water by quantum Monte Carlo Sandro Sorella G. Mazzola & Y. Luo SISSA, IOM DEMOCRITOS, Trieste A. Zen, L. Guidoni U. of L Aquila, L Aquila 28 July 2014, Mike Towler Institute,

More information

Module 6 1. Density functional theory

Module 6 1. Density functional theory Module 6 1. Density functional theory Updated May 12, 2016 B A DDFT C K A bird s-eye view of density-functional theory Authors: Klaus Capelle G http://arxiv.org/abs/cond-mat/0211443 R https://trac.cc.jyu.fi/projects/toolbox/wiki/dft

More information

PAPER:2, PHYSICAL CHEMISTRY-I QUANTUM CHEMISTRY. Module No. 34. Hückel Molecular orbital Theory Application PART IV

PAPER:2, PHYSICAL CHEMISTRY-I QUANTUM CHEMISTRY. Module No. 34. Hückel Molecular orbital Theory Application PART IV Subject PHYSICAL Paper No and Title TOPIC Sub-Topic (if any), PHYSICAL -II QUANTUM Hückel Molecular orbital Theory Module No. 34 TABLE OF CONTENTS 1. Learning outcomes. Hückel Molecular Orbital (HMO) Theory

More information

3: Many electrons. Orbital symmetries. l =2 1. m l

3: Many electrons. Orbital symmetries. l =2 1. m l 3: Many electrons Orbital symmetries Atomic orbitals are labelled according to the principal quantum number, n, and the orbital angular momentum quantum number, l. Electrons in a diatomic molecule experience

More information

Uncertainty in Molecular Photoionization!

Uncertainty in Molecular Photoionization! Uncertainty in Molecular Photoionization! Robert R. Lucchese! Department of Chemistry! Texas A&M University Collaborators:! At Texas A&M: R. Carey, J. Lopez, J. Jose! At ISMO, Orsay, France: D. Dowek and

More information

T. Helgaker, Department of Chemistry, University of Oslo, Norway. T. Ruden, University of Oslo, Norway. W. Klopper, University of Karlsruhe, Germany

T. Helgaker, Department of Chemistry, University of Oslo, Norway. T. Ruden, University of Oslo, Norway. W. Klopper, University of Karlsruhe, Germany 1 The a priori calculation of molecular properties to chemical accuarcy T. Helgaker, Department of Chemistry, University of Oslo, Norway T. Ruden, University of Oslo, Norway W. Klopper, University of Karlsruhe,

More information

Noncollinear spins in QMC: spiral Spin Density Waves in the HEG

Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Noncollinear spins in QMC: spiral Spin Density Waves in the HEG Zoltán Radnai and Richard J. Needs Workshop at The Towler Institute July 2006 Overview What are noncollinear spin systems and why are they

More information

Hückel Molecular orbital Theory Application PART III

Hückel Molecular orbital Theory Application PART III Subject PHYSICAL Paper No and Title TOPIC Sub-Topic (if any) 2, PHYSICAL -I QUANTUM Hückel Molecular orbital Theory Module No. 33 PAPER: 2, PHYSICAL -I TABLE OF CONTENTS 1. Learning outcomes 2. Hückel

More information

AN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University

AN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,

More information

Practical Issues on the Use of the CASPT2/CASSCF Method in Modeling Photochemistry: the Selection and Protection of an Active Space

Practical Issues on the Use of the CASPT2/CASSCF Method in Modeling Photochemistry: the Selection and Protection of an Active Space Practical Issues on the Use of the CASPT2/CASSCF Method in Modeling Photochemistry: the Selection and Protection of an Active Space Roland Lindh Dept. of Chemistry Ångström The Theoretical Chemistry Programme

More information

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah 1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations

More information

Chapter 2 Quantum chemistry using auxiliary field Monte Carlo

Chapter 2 Quantum chemistry using auxiliary field Monte Carlo Chapter 2 Quantum chemistry using auxiliary field Monte Carlo 1. The Hubbard-Stratonovich Transformation 2. Neuhauser s shifted contour 3. Calculation of forces and PESs 4. Multireference AFMC 5. Examples

More information

2~:J~ -ryej- r- 2 Jr. A - f3. sr(djk nv~tor rn~ +~ rvjs (::-CJ) ::;-1-.'--~ -. rhd. ('-.Ji.L.~ )- r'-d)c, -r/~ JJr - 2~d ~2-Jr fn'6.

2~:J~ -ryej- r- 2 Jr. A - f3. sr(djk nv~tor rn~ +~ rvjs (::-CJ) ::;-1-.'--~ -. rhd. ('-.Ji.L.~ )- r'-d)c, -r/~ JJr - 2~d ~2-Jr fn'6. .~, ~ I, sr(djk nv~tor rn~ +~ rvjs (::-CJ) ::;-1-.'--~ -. rhd. ('-.Ji.L.~ )- r'-d)c, -r/~ JJr - 2~d ~2-Jr fn'6.)1e'" 21t-ol Je C'...-------- lj-vi, J? Jr Jr \Ji 2~:J~ -ryej- r- 2 Jr A - f3 c _,~,= ~,.,w._..._.

More information

Towards gas-phase accuracy for condensed phase problems

Towards gas-phase accuracy for condensed phase problems Towards gas-phase accuracy for condensed phase problems Fred Manby Centre for Computational Chemistry, School of Chemistry University of Bristol STC 2006: Quantum Chemistry Methods and Applications Erkner,

More information

Quantum Monte Carlo Benchmarks Density Functionals: Si Defects

Quantum Monte Carlo Benchmarks Density Functionals: Si Defects Quantum Monte Carlo Benchmarks Density Functionals: Si Defects K P Driver, W D Parker, R G Hennig, J W Wilkins (OSU) C J Umrigar (Cornell), R Martin, E Batista, B Uberuaga (LANL), J Heyd, G Scuseria (Rice)

More information

Coupled-Cluster Perturbative Triples for Bond Breaking

Coupled-Cluster Perturbative Triples for Bond Breaking Coupled-Cluster Perturbative Triples for Bond Breaking Andrew G. Taube and Rodney J. Bartlett Quantum Theory Project University of Florida INT CC Meeting Seattle July 8, 2008 Why does chemistry need triples?

More information

Renner-Teller Effect in Tetra-Atomic Molecules

Renner-Teller Effect in Tetra-Atomic Molecules Groupe de Chimie Théorique du MSME Renner-Teller Effect in Tetra-Atomic Molecules Laurent Jutier, G. Dhont, H. Khalil and C. Léonard jutier@univ-mlv.fr (non linear) Outline General Presentation Structure

More information

Computational Chemistry I

Computational Chemistry I Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of

More information

Module:32,Huckel Molecular Orbital theory- Application Part-II PAPER: 2, PHYSICAL CHEMISTRY-I QUANTUM CHEMISTRY

Module:32,Huckel Molecular Orbital theory- Application Part-II PAPER: 2, PHYSICAL CHEMISTRY-I QUANTUM CHEMISTRY Subject PHYSICAL Paper No and Title TOPIC Sub-Topic (if any) 2, PHYSICAL -I QUANTUM Hückel Molecular orbital Theory Application PART II Module No. 32 TABLE OF CONTENTS 1. Learning outcomes 2. Hückel Molecular

More information

Basis sets for electron correlation

Basis sets for electron correlation Basis sets for electron correlation Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry

More information

Molecular Orbital Theory. Molecular Orbital Theory: Electrons are located in the molecule, not held in discrete regions between two bonded atoms

Molecular Orbital Theory. Molecular Orbital Theory: Electrons are located in the molecule, not held in discrete regions between two bonded atoms Molecular Orbital Theory Valence Bond Theory: Electrons are located in discrete pairs between specific atoms Molecular Orbital Theory: Electrons are located in the molecule, not held in discrete regions

More information

High-level Quantum Chemistry Methods and Benchmark Datasets for Molecules

High-level Quantum Chemistry Methods and Benchmark Datasets for Molecules High-level Quantum Chemistry Methods and Benchmark Datasets for Molecules Markus Schneider Fritz Haber Institute of the MPS, Berlin, Germany École Polytechnique Fédérale de Lausanne, Switzerland دانشگاه

More information

Exercise 1: Structure and dipole moment of a small molecule

Exercise 1: Structure and dipole moment of a small molecule Introduction to computational chemistry Exercise 1: Structure and dipole moment of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the dipole moment of a small

More information

The wavefunction that describes a bonding pair of electrons:

The wavefunction that describes a bonding pair of electrons: 4.2. Molecular Properties from VB Theory a) Bonding and Bond distances The wavefunction that describes a bonding pair of electrons: Ψ b = a(h 1 ) + b(h 2 ) where h 1 and h 2 are HAOs on adjacent atoms

More information

Introduction to Path Integral Monte Carlo. Part I.

Introduction to Path Integral Monte Carlo. Part I. Introduction to Path Integral Monte Carlo. Part I. Alexey Filinov, Jens Böning, Michael Bonitz Institut für Theoretische Physik und Astrophysik, Christian-Albrechts-Universität zu Kiel, D-24098 Kiel, Germany

More information

Introduction to multiconfigurational quantum chemistry. Emmanuel Fromager

Introduction to multiconfigurational quantum chemistry. Emmanuel Fromager Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Notations

More information

The QMC Petascale Project

The QMC Petascale Project The QMC Petascale Project Richard G. Hennig What will a petascale computer look like? What are the limitations of current QMC algorithms for petascale computers? How can Quantum Monte Carlo algorithms

More information

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah 1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations

More information

Supporting Information

Supporting Information Supporting Information Conflict in the Mechanism and Kinetics of the Barrierless Reaction between SH and NO 2 Radicals Ramanpreet Kaur and Vikas * Quantum Chemistry Group, Department of Chemistry & Centre

More information

Fractional-charge and Fractional-spin errors in Range-Separated Density-Functional Theory with Random Phase Approximation.

Fractional-charge and Fractional-spin errors in Range-Separated Density-Functional Theory with Random Phase Approximation. Fractional-charge and Fractional-spin errors in Range-Separated Density-Functional Theory with Random Phase Approximation. Bastien Mussard, Julien Toulouse Laboratoire de Chimie Théorique Institut du Calcul

More information

MOLECULAR ORBITAL THEORY Chapter 10.8, Morrison and Boyd

MOLECULAR ORBITAL THEORY Chapter 10.8, Morrison and Boyd MOLECULAR ORBITAL THEORY Chapter 10.8, Morrison and Boyd more understanding: why oxygen is paramagnetic, why H2 + exists; explanation of excited electronic states (e.g., visible spectra) eliminates need

More information

Computational Methods. Chem 561

Computational Methods. Chem 561 Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational

More information

Methods for Treating Electron Correlation CHEM 430

Methods for Treating Electron Correlation CHEM 430 Methods for Treating Electron Correlation CHEM 430 Electron Correlation Energy in the Hartree-Fock approximation, each electron sees the average density of all of the other electrons two electrons cannot

More information

Introduction to Electronic Structure Theory

Introduction to Electronic Structure Theory CSC/PRACE Spring School in Computational Chemistry 2017 Introduction to Electronic Structure Theory Mikael Johansson http://www.iki.fi/~mpjohans Objective: To get familiarised with the, subjectively chosen,

More information

Atomic Structure and Bonding. Chapter 1 Organic Chemistry, 8 th Edition John McMurry

Atomic Structure and Bonding. Chapter 1 Organic Chemistry, 8 th Edition John McMurry Atomic Structure and Bonding Chapter 1 Organic Chemistry, 8 th Edition John McMurry 1 Common Elements Groups First row Second row In most organic molecules carbon is combined with relatively few elements

More information

Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid

Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid Quantum Monte Carlo study of the ground state of the two-dimensional Fermi fluid N. D. Drummond and R. J. Needs TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge

More information

Chapter 3. Orbitals and Bonding

Chapter 3. Orbitals and Bonding Chapter 3. Orbitals and Bonding What to master Assigning Electrons to Atomic Orbitals Constructing Bonding and Antibonding Molecular Orbitals with Simple MO Theory Understanding Sigma and Pi Bonds Identifying

More information

Department of Physics and NCSA University of Illinois, Urbana-Champaign, IL, 61801, USA

Department of Physics and NCSA University of Illinois, Urbana-Champaign, IL, 61801, USA UNDERSTANDING ELECTRONIC WAVE FUNCTIONS D. M. Ceperley Department of Physics and NCSA University of Illinois, Urbana-Champaign, IL, 61801, USA INTRODUCTION In this article I discuss some aspects of what

More information

Quantum Monte Carlo calculations of two neutrons in finite volume

Quantum Monte Carlo calculations of two neutrons in finite volume Quantum Monte Carlo calculations of two neutrons in finite volume Philipp Klos with J. E. Lynn, I. Tews, S. Gandolfi, A. Gezerlis, H.-W. Hammer, M. Hoferichter, and A. Schwenk Nuclear Physics from Lattice

More information

Correlation in correlated materials (mostly transition metal oxides) Lucas K. Wagner University of Illinois at Urbana-Champaign

Correlation in correlated materials (mostly transition metal oxides) Lucas K. Wagner University of Illinois at Urbana-Champaign Correlation in correlated materials (mostly transition metal oxides) Lucas K. Wagner University of Illinois at Urbana-Champaign Understanding of correlated materials is mostly phenomenological FN- DMC

More information

MO theory is better for spectroscopy (Exited State Properties; Ionization)

MO theory is better for spectroscopy (Exited State Properties; Ionization) CHEM 2060 Lecture 25: MO Theory L25-1 Molecular Orbital Theory (MO theory) VB theory treats bonds as electron pairs. o There is a real emphasis on this point (over-emphasis actually). VB theory is very

More information

CHEM 2010 Symmetry, Electronic Structure and Bonding Winter Numbering of Chapters and Assigned Problems

CHEM 2010 Symmetry, Electronic Structure and Bonding Winter Numbering of Chapters and Assigned Problems CHEM 2010 Symmetry, Electronic Structure and Bonding Winter 2011 Numbering of Chapters and Assigned Problems The following table shows the correspondence between the chapter numbers in the full book (Physical

More information

Shapes of Molecules. Lewis structures are useful but don t allow prediction of the shape of a molecule.

Shapes of Molecules. Lewis structures are useful but don t allow prediction of the shape of a molecule. Shapes of Molecules Lewis structures are useful but don t allow prediction of the shape of a molecule. H O H H O H Can use a simple theory based on electron repulsion to predict structure (for non-transition

More information

PART 3 Chemical Bonds, Valence Bond Method, and Molecular Shapes. Reference: Chapter 9 10 in textbook

PART 3 Chemical Bonds, Valence Bond Method, and Molecular Shapes. Reference: Chapter 9 10 in textbook PART 3 Chemical Bonds, Valence Bond Method, and Molecular Shapes Reference: Chapter 9 10 in textbook 1 Valence Electrons Valence ae Electron Define: the outer shell electrons Important for determination

More information

Chem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.

Chem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components. Chem 44 Review for Exam Hydrogenic atoms: The Coulomb energy between two point charges Ze and e: V r Ze r Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative

More information

Many-body wavefunctions for normal liquid 3 He

Many-body wavefunctions for normal liquid 3 He Many-body wavefunctions for normal liquid 3 He Markus Holzmann, 1 Bernard Bernu, 1 and D. M. Ceperley 2 1 LPTMC, UMR 7600 of CNRS, Université Pierre et Marie Curie, Paris, France 2 Department of Physics

More information

More quantum chemistry with BOINC

More quantum chemistry with BOINC More quantum chemistry with BOINC Martin Korth Theory of Condensed Matter group, Cavendish Laboratory University of Cambridge The 6th BOINC Workshop 01/09/2010 From Chemistry to Quantum Chemistry Chemistry

More information

Molecular Magnetism. Magnetic Resonance Parameters. Trygve Helgaker

Molecular Magnetism. Magnetic Resonance Parameters. Trygve Helgaker Molecular Magnetism Magnetic Resonance Parameters Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Laboratoire de Chimie Théorique,

More information

We can model covalent bonding in molecules in essentially two ways:

We can model covalent bonding in molecules in essentially two ways: CHEM 2060 Lecture 22: VB Theory L22-1 PART FIVE: The Covalent Bond We can model covalent bonding in molecules in essentially two ways: 1. Localized Bonds (retains electron pair concept of Lewis Structures)

More information

Molecular Magnetic Properties

Molecular Magnetic Properties Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry and

More information

Molecular Geometry and Chemical Bonding Theory

Molecular Geometry and Chemical Bonding Theory Molecular Geometry and Chemical Bonding Theory The Valence -Shell Electron -Pair Repulsion (VSEPR) Model predicts the shapes of the molecules and ions by assuming that the valence shell electron pairs

More information

Random Estimates in QMC

Random Estimates in QMC Page 1 Random Estimates in QMC J.R. Trail and R.J. Needs Theory of Condensed Matter Group, Cavendish Laboratory, University of Cambridge, UK February 2007 Page 2 Numerical Integration Unkown function,

More information

New applications of Diffusion Quantum Monte Carlo

New applications of Diffusion Quantum Monte Carlo New applications of Diffusion Quantum Monte Carlo Paul R. C. Kent (ORNL) Graphite: P. Ganesh, J. Kim, C. Park, M. Yoon, F. A. Reboredo (ORNL) Platinum: W. Parker, A. Benali, N. Romero (ANL), J. Greeley

More information

The challenges. Challenge 1: Give me insight and good numbers, on another unusual bonding

The challenges. Challenge 1: Give me insight and good numbers, on another unusual bonding The challenges hallenge 1: Give me insight and good numbers, on another unusual bonding hallenge 2: can VB be a simpler, easier, and faster method than MO-I to provides accurate numbers? Stop me at any

More information

Dispersion Interactions from the Exchange-Hole Dipole Moment

Dispersion Interactions from the Exchange-Hole Dipole Moment Dispersion Interactions from the Exchange-Hole Dipole Moment Erin R. Johnson and Alberto Otero-de-la-Roza Chemistry and Chemical Biology, University of California, Merced E. R. Johnson (UC Merced) Dispersion

More information

Atomic Wave Function Forms. Received 3 June 1996; revised 24 December 1997; accepted 2 January 1997

Atomic Wave Function Forms. Received 3 June 1996; revised 24 December 1997; accepted 2 January 1997 Atomic Wave Function Forms S. A. ALEXANDER, 1 R. L. COLDWELL 2 1 Department of Physics, University of Texas at Arlington, Arlington, Texas 76019 2 Department of Physics, University of Florida, Gainesville,

More information