Coupled-Cluster Perturbative Triples for Bond Breaking
|
|
- Marcia Warren
- 6 years ago
- Views:
Transcription
1 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
2 Why does chemistry need triples? Method/ Molecule Hartree-Fock H2O N2 HF (99.7%) (99.7%) (99.9%) Total energies in hartree CCSD (99.995%) (99.987%) (99.998%) FCI hartree = 27.2 ev Scale (for these molecules) of total energies: 10 2 hartree = 10 0 kev Equilibrium structures. H2O in cc-pvdz, N2 in cc-pvdz, HF in DZ
3 Atomization Energies 16 Small Molecules, cc-pcvqz Bak, et. al., JCP 112, 9229 (2000) H2O 2H + O Scale of AE ~ kj/mol ~ 10 0 ev ev ;<=>5&:??@A??@A5>: 1 ev = 96.5 kj/mol Errors in AE ~ kj/mol ~ 10 1 mev mev * *!!""!#"!$"!%"!&" " &" %" $" #"!"" '(()(*+,*-.)/+01.+),*',2(34*5678/)9:
4 19 Rxns of those 16 small molecules, Ibid. Reactions <= BBCD BBCD6A; * Scale of ΔErxn ~ kj/mol ~ 10 2 mev ev Errors in ΔErxn ~ kj/mol ~ 10 0 mev mev *!!""!#"!$"!%"!&" " &" %" $" #"!"" '(()(*+,*-./01+),*',12/34+.5*6789:)3;
5 19 Rxns of those 16 small molecules, Ibid. Reactions <= BBCD BBCD6A; * Scale of ΔErxn ~ kj/mol ~ 10 2 mev ev Errors in ΔErxn ~ kj/mol ~ 10 0 mev mev Total Energy ~ 1 kev want errors in ΔErxn ~ 1 mev *!!""!#"!$"!%"!&" " &" %" $" #"!"" '(()(*+,*-./01+),*',12/34+.5*6789:)3;
6 What can we (reasonably) expect from high-level quantum chemistry Identification of and calculation of properties at stationary points of moderately-sized molecules Parameterization / reference for less accurate methods for larger molecules Driver for Molecular Dynamics / (Classical) Monte Carlo Simulations for small molecules
7 What can we (reasonably) expect from high-level quantum chemistry Identification of and calculation of properties at stationary points of moderately-sized molecules Parameterization / reference for less accurate methods for larger molecules Driver for Molecular Dynamics / (Classical) Monte Carlo Simulations for small molecules For the latter two, need data away from stationary points For MC: Energies For MD: Forces
8 Computational Scaling CCSD T=T1+T2 CCSD(T) T=T1+T2+T3 [2] CCSDT T=T1+T2+T3 O(N 6 ) ( o 3 v 3 + o 2 v 4) N it O(N 7 ) ( o 3 v 3 + o 2 v 4) N it + o 3 v 4 + o 4 v 3 O(N 8 ) ( o 4 v 4 + o 3 v 5) N it
9 The Gold Standard E CCSD = E ref + 0 F ov T W ( T ) 2 T E CCSD(T) = E CCSD + 0 ( ) T 1 + T 2 (F ov + W ) R 3 (W T 2 ) C 0 R 3 = ( ) 2 1 3! ijk abc abc ijk abc ijk ɛ abc ijk [2] = [2] [2] [2]
10 CCSD(T)
11 CCSD(T)
12 CCSD(T)
13 There are two different ways that CCSD(T) fails for bond breaking
14 Bond Length (n*re) There are two different ways that CCSD(T) fails for bond breaking Collapse of the spin-restricted (RHF) wavefunction RHF CCSD RHF CCSD(T) DMRG Full CI 376.5!E (kcal/mol)
15 There are two different ways that CCSD(T) fails for bond breaking Collapse of the spin-restricted (RHF) wavefunction Spin symmetry breaking RHF CCSD(T) UHF CCSD(T) Absolute Error from Full CI (kcal/mol) Bond Length (n*re)
16 F2 at 2Re Denominatorless T3 amplitudes δe t abc ijk RHF λ ijk abc [2] ɛ abc ijk δe t abc ijk UHF λ ijk abc ɛ abc ijk The T3 amplitudes themselves are bad - not denominators ( E CCSD = E ref + 0 F ov T W T ) 2 T O(T 2 ) ( ) E CCSD(T) = E CCSD + 0 T 1 + T 2 (F ov + W ) R 3 (W T 2 ) C 0 O(T2 2 ) If the CCSD wavefunction is bad then so is the CCSD(T) energy
17 Ways to fix this problem
18 Ways to fix this problem Symmetry-broken Reference Doesn t fully fix the problem
19 Ways to fix this problem Symmetry-broken Reference Multireference Techniques Not Black-box, Extensivity, Computational Cost
20 Ways to fix this problem Symmetry-broken Reference Multireference Techniques Higher excitations Computational Cost
21 Ways to fix this problem Symmetry-broken Reference Multireference Techniques Higher excitations Use CI Extensivity, Accuracy
22 Ways to fix this problem Symmetry-broken Reference Multireference Techniques Higher excitations Use CI Improved Functionals Our Choice
23 Properties of CCSD(T) that are important to keep Extensivity Invariance to OO and VV rotations Availability of analytical gradients O[(o 3 v 3 + o 2 v 4 )N it + o 3 v 4 + o 4 v 3 ] Applicability to arbitrary reference functions
24 Löwdin Partitioning H P P c P + H P Q c Q H QP c P + H QQ c Q = Ec P = Ec Q c Q = ˆQ ( E H QQ ) 1 ˆQ HQP c P H P P + H P Q ˆQ ( E HQQ ) 1 ˆQ HQP c P = Ec P H is a similarity transformation of the Hamiltonian same eigenvalue spectrum
25 Perturbative Expansion E = 0 L P HP P R P L P HP Q ˆQ ( E HQQ ) 1 ˆQ HQP R P 0 Want corrections to CCSD Expand perturbatively around CCSD solution H [0] P P = ˆP H CCSD ˆP E [0] = E CCSD H[0] P Q L [0] P = (1 + Λ CCSD) R [0] P = 1 = H [0] QP = 0 CCSD(T) Choice ) L [0] P (1 + T CCSD
26 Partitioning Choices H [0] Λ(T) = ˆP H CCSD ˆP + ˆQ H [0] (2) T = ˆP H CCSD ˆP + ˆQ ( E CCSD + H (0)) ˆQ ( E CCSD + H (0)) ˆQ H [0] CR(2,3) = ˆP H CCSD ˆP + ˆQ ( ECCSD + H CCSD ) ˆQ H [1] Λ(T) = H (1) CCSD H [1] (2) T = H CCSD H [1] CR(2,3) = H CCSD H [0] Λ(T) H [0] (2) T H [0] CR(2,3) H (0) = F oo + F vv H (1) = F ov + F vo + W
27 (T)-Like Functionals ( ) 2 E [3] 1 CCSD(T) = 3! ijk E [3] ΛCCSD(T) = ( 1 3! E [2] CCSD(2) T = E [2] CR(2,3) = ( 1 3! ( 1 3! abc ) 2 ijk abc ) 2 ijk abc ) 2 ijk abc 0 T (F ov + W ) abc ijk 0 Λ (F ov + W ) abc ijk 0 Λ H SD abc ijk 0 Λ H SD abc ijk 1 ɛ abc ijk 1 ɛ abc ijk 1 ɛ abc ijk 1 ɛ abc ijk abc ijk (W T 2 ) C 0 abc ijk (W T 2 ) C 0 abc ijk H SD 0 abc ijk H SD 0 Formally should have quadruples CCSD(T): Many CCSD(2)T: Hirata ΛCCSD(T): Bartlett CR(2,3) : Piecuch
28 Non-iterative O(N 7 ) Triples ΛCCSD(T) is the minimal extension to CCSD(T), yet... ΛCCSD(T) CCSD(T) CCSD(2)T CR(2,3) Improved RHF bond breaking Orbital Invariance Non-HF Reference Analytical Gradients * CR-CCSD(T) is not extensive
29 ΛCCSD(T) and CCSD(T) E ΛCCSD(T) = E CCSD + 0 (Λ 1 + Λ 2 ) (F ov + W ) R 3 (W T 2 ) C 0 [2] [2] [2] E CCSD(T) = E CCSD + 0 ( ) T 1 + T 2 (F ov + W ) R 3 (W T 2 ) C 0 [2] [2] [2]
30 Stationary Formulation of ΛCCSD(T) E ΛCCSD(T) = 0 (1 + Λ + Π) H (1 + Λ) ( HΣ )C Λ 1 W T [2] Λ 2 (F ov + W ) T [2] Λ [2] 3 (F oo + F vv ) T [2] Λ [2] 3 (W T 2) C 0 Allows for (Generalized) Hellman-Feynman Theorem to be used for gradients
31 Stationary Formulation of ΛCCSD(T) E ΛCCSD(T) = 0 (1 + Λ + Π) H (1 + Λ) ( HΣ )C 0 Combination of the Stat. Functional for CCSD(T) and EOM-CCSD + 0 Λ 1 W T [2] Λ 2 (F ov + W ) T [2] Λ [2] 3 (F oo + F vv ) T [2] Λ [2] 3 (W T 2) C 0 Allows for (Generalized) Hellman-Feynman Theorem to be used for gradients
32 Reaction Energies CCSD(T)!CCSD(T) Error in Reaction Enthalpy (kcal/mol) ~40 well-characterized reactions from Lynch, Zhao, Truhlar data set
33 Lithium Fluoride RHF CCSD RHF CCSD(T) RHF "CCSD(T) DMRG Full CI!E (kcal/mol) double-ζ Bond Length (n*re) FCI: Bauschlicher, Langoff, JCP 89, 4246 (1988). DMRG: Legeza, Roder, Hess, Mol. Phys. 101, 2019 (2003).
34 Carbon Monoxide RHF CCSD RHF CCSD(T) RHF "CCSD(T) cc-pvtz 284.7!E (kcal/mol) De ~ 8 kcal/mol too high Bond Length (n*re)
35 Carbon Monoxide RHF CCSD RHF CCSD(T) RHF "CCSD(T) cc-pvtz 284.7!E (kcal/mol) De ~ 8 kcal/mol too high Bond Length (n*re)
36 H2O RHF 0 cc-pvdz Error from Full CI (kcal/mol) RHF CCSD(T) RHF!CCSD(T) RHF CR-CCSD(T) RHF CCSD(2)T RHF CR-CCSD(T)L Bond Length (n*re) FCI: Olsen, et. al., JCP 104, 8007 (1996). CR-CCSD(T)/CCSD(2)T: Hirata, et. al., JCP 121, (2004 ). CR-CCSD(T) L: Włoch, et. al., MP 104, 2149 (2006)
37 Average Non-parallelity Error CCSD CCSD(T)!CCSD(T) Weighted Non-Parallelity Error (millihartree) HF: Dutta, Sherrill, JCP 118, 1610 (2003). H2O: Olsen, et al., JCP 104, 8007 (1996). C2: Sherrill, Piecuch, JCP 122, (2005). N2: Chan, et al., JCP 121, 6110 (2004).
38 HF Force Curve RHF CCSD(T) RHF!CCSD(T) 6-31G** Force (E h /a 0 ) Bond Length (n*r e )
39 FCI: Olsen, et. al., JCP 104, 8007 (1996). CR-CCSD(T)/CCSD(2)T: Hirata, et. al., JCP 121, (2004 ). H2O UHF 7.53 cc-pvdz UHF CCSD(T) UHF!CCSD(T) UHF CCSD(2)T Error from Full CI (kcal/mol) Bond Length (n*re)
40 [2] Denominatorless T3 amplitudes δe t abc ijk UHF λ ijk abc T3 amplitudes not source of the UHF problem ɛ abc ijk Error from RHF CCSDT (kcal/mol) F2 at 2Re RHF CCSD(T) RHF!CCSD(T) UHF CCSD(T) UHF!CCSD(T) Bond Length (n*r e )
41 RHF: The ΛCCSD(T) Story So Far Improvements for intermediate bond lengths Asymptotics not quantitative UHF: No Improvements for intermediate bond lengths Asymptotics still ok ΛCCSD(T) fixes spin-restricted PT failure Not spin-symmetry breaking problem AGT, RJB. JCP 128, & (2008)
42 Can spin-restricted refs be good longer? Yes... Brueckner Orbitals Best single-determinant approximation to the exact wavefunction 0 B min 0 0 Ψ 2 equivalent to: 0 B max 0 0 Ψ * For normalized wfns and determinants Brueckner orbitals should maintain symmetry better than Hartree-Fock orbitals
43 CC Condition What about the case using a CC wfn? Intermediate normalization: Ω = e κ 1 0 Ψ Ψ 1 κ 1 = ia κ a i ( {a i} {i a} ) Define f(κ 1 ) = Ω e κ 1 Ψ Ψ Ψ 1/2 0 B max κ 1 f(κ 1 )
44 CC Condition, continued f(κ 1 ) κ a i Expanding κ1 =0 = Ω ( i a a i ) e κ 1 Ψ Ψ Ψ 1/2 = Ωa i Ψ Ψ Ψ 1/2 = 0 κ1 =0 Ψ = e T Ω = Ω + ia t a i Ω a i +
45 CC Condition, continued f(κ 1 ) κ a i Expanding κ1 =0 = Ω ( i a a i ) e κ 1 Ψ Ψ Ψ 1/2 = Ωa i Ψ Ψ Ψ 1/2 = 0 κ1 =0 Ψ = e T Ω = Ω + ia t a i Ω a i + t a i = 0
46 Brueckner Instability Condition 2 f(κ 1 ) κ a i κb j κ1 =0 = Ω ( i a a i ) ( j b b j ) e κ 1 Ψ Ψ Ψ 1/2 Ω ab ij = Ψ δ ij δ ab Ω Ψ Ψ Ψ 1/2 = Ψ Ψ 1/2 ( t ab ) ij δ ij δ ab κ1 =0 Want this to be negative definite (maximum) Diagonalize A = I T 2 Any negative eigenvalues are instabilities
47 Onset of Instability for BCCD Molecule N2 F2 C2 HF RHF UHF 1.1Re Re 0.5Re 1.5Re RB UB 1.9Re >5Re 2.1Re >5Re
48 Computational Procedure Brueckner Instability Analysis Brueckner DIIS (Padé Approximants for matrix exponentials and logarithms) Regularization of CC equations Initial guess using prior converged brueckner result
49 An interlude on solving the CC equations [n] t ab ij = 1 ɛ ab ij F ( [n 1] t, [n 1] t 2, [n 1] t 3, [n 1] t 4) Note the possibility of a (near-)singularity due to the denominator [n] t ab ij = ( ɛ ab ij + ɛ ab ij ) 2 + δ 2 F ( [n 1] t, [n 1] t 2, [n 1] t 3, [n 1] t 4) ( ɛ ab ij δ 2 ) 2 + δ 2 [n 1] t ab ij Based on an approximate Tikhonov regularization. Will converge to the same result as unregularized methods
50 t ab ij = At convergence of CC equations... [ (ɛ ab ij ɛ ab ij ( ( ) ɛ ab F t, t 2, t 3, t 4) δ 2 ij ( ɛ ab ij δ 2 ) 2 + δ 2 tab ) ] 2 + δ 2 t ab ij = ɛ ab ( ij F t, t 2, t 3, t 4) + δ 2 t ab ij ( ) ɛ ab 2 ij t ab ij = ɛ ab ( ij F t, t 2, t 3, t 4) t ab ij = 1 ɛ ab ij F ( t, t 2, t 3, t 4) ij
51 Brueckner Functional E BCC = 0 (1 + Λ) [ (He κ 1 ) e T ] C C 0 T = n=2 T n κ 1 = κ 1 1 Λ = n=1 Λ n BΛCCSD(T) same as ΛCCSD(T) with modified Λ Singles are needed for Λ, not for T equations Reformulation of Handy, Pople, Head-Gordon, Raghavachari, Trucks. CPL 164, 185 (1989)
52 HF RHF / UHF / RB G** Error from Full CI (kcal/mol) RHF CCSD(T) RHF!CCSD(T) RB CCD(T) RB!CCSD(T) UHF CCSD(T) UHF!CCSD(T) Bond Length (n*re) FCI: Dutta, Sherrill, JCP 118, 1610 (2003).
53 H2O UHF / UB Error from Full CI (kcal/mol) cc-pvdz UHF CCSD(T) UHF!CCSD(T) UB CCD(T) UB!CCSD(T) Bond Length (n*re) FCI:Olsen, et. al., JCP 104, 8007 (1996).
54 F2 RHF / UHF / RB 5 cc-pvdz Error from RHF CCSDT (kcal/mol) RHF CCSD(T) RHF!CCSD(T) RB CCD(T) RB!CCSD(T) UHF CCSD(T) UHF!CCSD(T) Bond Length (n*r e ) CCSDT: Kowalski, Piecuch, CPL 344, 165 (2001).
55 N2 UHF / UB 15 cc-pvdz 10 Error from FCI (kcal/mol) UHF CCSD(T) UHF!CCSD(T) UB CCD(T) UB!CCSD(T) Bond Length (n*r e ) FCI: Chan, et al., JCP 121, 6110 (2004).
56 Conclusions about BΛCCSD(T) Turned Spin-Symmetry Breaking (T) Problem into Spin-Restricted (T) Problem Allows BΛCCSD(T) PES to be Substantially improved over HF CCSD(T) or BCCD(T) Brueckner orbitals improve perturbative approximations; no evidence this is true for iterative CC methods
57 Acknowledgments Symposium Organizers and Staff Prof. Rod Bartlett and the Bartlett Group All of you Funding: National Defense Science and Engineering Fellowship AFOSR, ARL
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 informationElectron 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 informationCOUPLED-CLUSTER CALCULATIONS OF GROUND AND EXCITED STATES OF NUCLEI
COUPLED-CLUSTER CALCULATIONS OF GROUND AND EXCITED STATES OF NUCLEI Marta Włoch, a Jeffrey R. Gour, a and Piotr Piecuch a,b a Department of Chemistry,Michigan State University, East Lansing, MI 48824 b
More informationFull configuration interaction potential energy curves for breaking bonds to hydrogen: An assessment of single-reference correlation methods
JOURNAL OF CHEMICAL PHYSICS VOLUME 118, NUMBER 4 22 JANUARY 2003 Full configuration interaction potential energy curves for breaking bonds to hydrogen: An assessment of single-reference correlation methods
More informationNWChem: Coupled Cluster Method (Tensor Contraction Engine)
NWChem: Coupled Cluster Method (ensor Contraction Engine) What we want to solve H Ψ = E Ψ Many Particle Systems Molecular/Atomic Physics, Quantum Chemistry (electronic Schrödinger equations) Solid State
More informationIntroduction 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 informationMethods 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 informationOther methods to consider electron correlation: Coupled-Cluster and Perturbation Theory
Other methods to consider electron correlation: Coupled-Cluster and Perturbation Theory Péter G. Szalay Eötvös Loránd University Institute of Chemistry H-1518 Budapest, P.O.Box 32, Hungary szalay@chem.elte.hu
More informationAN 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 informationChapter 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 informationMolecular properties in quantum chemistry
Molecular properties in quantum chemistry Monika Musiał Department of Theoretical Chemistry Outline 1 Main computational schemes for correlated calculations 2 Developmentoftheabinitiomethodsforthe calculation
More informationThe Accurate Calculation of Molecular Energies and Properties: A Tour of High-Accuracy Quantum-Chemical Methods
1 The Accurate Calculation of Molecular Energies and Properties: A Tour of High-Accuracy Quantum-Chemical Methods T. Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry,
More informationLecture 5: More about one- Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory.
Lecture 5: More about one- determinant wave functions Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory. Items from Lecture 4 Could the Koopmans theorem
More informationHighly accurate quantum-chemical calculations
1 Highly accurate quantum-chemical calculations T. Helgaker Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, Norway A. C. Hennum and T. Ruden, University
More informationNWChem: Coupled Cluster Method (Tensor Contraction Engine)
NWChem: Coupled Cluster Method (Tensor Contraction Engine) Why CC is important?! Correlation effects are important!! CC is size-extensive theory: can be used to describe dissociation processes.! Higher-order
More informationCoupled-cluster and perturbation methods for macromolecules
Coupled-cluster and perturbation methods for macromolecules So Hirata Quantum Theory Project and MacroCenter Departments of Chemistry & Physics, University of Florida Contents Accurate electronic structure
More informationConvergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems
1 Convergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems T. Helgaker Centre for Theoretical and Computational Chemistry, Department of
More informationComparison of Nuclear Configuration Interaction Calculations and Coupled Cluster Calculations
Comparison of Nuclear Configuration Interaction Calculations and Coupled Cluster Calculations Mihai Horoi Department of Physics, Central Michigan University, Mount Pleasant, Michigan 48859, USA Support
More informationTime-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 informationOVERVIEW 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 informationJack 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 information4 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 informationPredictive Computing for Solids and Liquids
Predictive Computing for Solids and Liquids So Hirata Department of Chemistry May 214 Blue Waters Symposium 1 Schrödinger equation for a water molecule 1-particle, 3-dimensional partial differential equation
More informationFaddeev Random Phase Approximation (FRPA) Application to Molecules
Faddeev Random Phase Approximation (FRPA) Application to Molecules Matthias Degroote Center for Molecular Modeling (CMM) Ghent University INT 2011 Spring Program Fermions from Cold Atoms to Neutron Stars:
More informationBenchmark studies of variational, unitary and extended coupled cluster methods
Benchmark studies of variational, unitary and extended coupled cluster methods Bridgette Cooper and Peter J. Knowles Citation: The Journal of Chemical Physics 133, 23412 (21); doi: 1.163/1.352564 View
More informationInternally Contracted Multireference Coupled Cluster Theory. Computer Aided Implementation of Advanced Electronic Structure Methods.
Internally Contracted Multireference Coupled Cluster Theory. Computer Aided Implementation of Advanced Electronic Structure Methods. Marcel Nooijen Chemistry Department University of Waterloo nooijen@uwaterloo.ca
More informationElectron 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 informationChemistry 4560/5560 Molecular Modeling Fall 2014
Final Exam Name:. User s guide: 1. Read questions carefully and make sure you understand them before answering (if not, ask). 2. Answer only the question that is asked, not a different question. 3. Unless
More informationBeyond 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 informationAb 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 informationTowards 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 information1 Rayleigh-Schrödinger Perturbation Theory
1 Rayleigh-Schrödinger Perturbation Theory All perturbative techniques depend upon a few simple assumptions. The first of these is that we have a mathematical expression for a physical quantity for which
More informationWave function methods for the electronic Schrödinger equation
Wave function methods for the electronic Schrödinger equation Zürich 2008 DFG Reseach Center Matheon: Mathematics in Key Technologies A7: Numerical Discretization Methods in Quantum Chemistry DFG Priority
More informationT. 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 informationMulti-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005
Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson
More informationApproximate variational coupled cluster theory
Approximate variational coupled cluster theory James B. Robinson and Peter J. Knowles Citation: The Journal of Chemical Physics 135, 044113 (2011); doi: 10.1063/1.3615060 View online: http://dx.doi.org/10.1063/1.3615060
More informationExercise 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 informationTheoretical and Computational Studies of Interstellar C2nH and SiC2m+1H. Ryan Fortenberry
Theoretical and Computational Studies of Interstellar C2nH and SiC2m+1H Ryan Fortenberry 1 Introduction Astrobiology Pillars of Creation Titan Interstellar Spectra DIBs 1.0 SiC3H Comparison Spectrum 0.8
More informationIntroduction 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 informationCOUPLED-CLUSTER BASED METHODS FOR EXCITED STATE ENERGIES AND GRADIENTS STEVEN RAY GWALTNEY
COUPLED-CLUSTER BASED METHODS FOR EXCITED STATE ENERGIES AND GRADIENTS By STEVEN RAY GWALTNEY A DISSERTATION PRESENTED TO THE GRADUATE SCHOOL OF THE UNIVERSITY OF FLORIDA IN PARTIAL FULFILLMENT OF THE
More informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More informationAb Initio Electromagnetic Transitions with the IMSRG
Ab Initio Electromagnetic Transitions with the IMSRG Nathan Parzuchowski Michigan State University March, 017 1 / 1 Outline IMSRG IMSRG rotates the Hamiltonian into a coordinate system where simple methods
More information$ ARO $ University of Florida: Quantum Theory Project
Approaching Multi-reference (and Open-Shell Problems) in Electronic Structure Theory from the Equation-of-Motion Coupled-Cluster Perspective INT Workshop on Advances in Many-Body Theory: From Nuclei to
More informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More informationQuantum Chemistry Methods
1 Quantum Chemistry Methods T. Helgaker, Department of Chemistry, University of Oslo, Norway The electronic Schrödinger equation Hartree Fock theory self-consistent field theory basis functions and basis
More informationOptimization 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 informationBasis 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 informationRecent 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 informationThe Rigorous Calculation of Molecular Properties to Chemical Accuracy. T. Helgaker, Department of Chemistry, University of Oslo, Norway
1 The Rigorous Calculation of Molecular Properties to Chemical Accuracy T. Helgaker, Department of Chemistry, University of Oslo, Norway A. C. Hennum and T. Ruden, University of Oslo, Norway S. Coriani,
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationIntroduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule. Vesa Hänninen
Introduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the electronic energy
More informationPutting a spin on time dependent electronic structure theory! Joshua Goings! General Exam! Thursday, May 14, 10:00am CHB 339
Putting a spin on time dependent electronic structure theory Joshua Goings General Exam Thursday, May 14, 10:00am CHB 339 Thank you to General Exam Committee Xiaosong Li (chair) Jim Pfaendtner (GSR) Matthew
More informationCoupled Cluster Theory and Tensor Decomposition Techniques
Coupled Cluster Theory and Tensor Decomposition Techniques Alexander A. Auer Atomistic Modelling Group Interface Chemistry and Surface Engineering Hierarchy of quantum chemical methods Analysis of orbitals
More informationSpin-restricted Brueckner orbitals for coupled-cluster wavefunctions
Spin-restricted Brueckner orbitals for coupled-cluster wavefunctions T. Daniel Crawford Center for Computational Quantum Chemistry University of Georgia thens Georgia 30602-2556 Timothy J. Lee NS mes Research
More informationTHE JOURNAL OF CHEMICAL PHYSICS 128,
THE JOURNAL OF CHEMICAL PHYSICS 128, 044108 2008 A comparative assessment of the perturbative and renormalized coupled cluster theories with a noniterative treatment of triple excitations for thermochemical
More informationAdvanced 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 informationChemistry 334 Part 2: Computational Quantum Chemistry
Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical
More informationTime-independent molecular properties
Time-independent molecular properties Trygve Helgaker Hylleraas Centre, Department of Chemistry, University of Oslo, Norway and Centre for Advanced Study at the Norwegian Academy of Science and Letters,
More informationCoupled-Cluster Theory. Nuclear Structure
Coupled-Cluster Theory! for Nuclear Structure!!!! Sven Binder INSTITUT FÜR KERNPHYSIK! 1 Nuclear Interactions from Chiral EFT NN 3N 4N NLO LO N 2 LO +... N 3 LO +... +... +... 2 Nuclear Interactions from
More informationPractical Advice for Quantum Chemistry Computations. C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology
Practical Advice for Quantum Chemistry Computations C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Choice of Basis Set STO-3G is too small 6-31G* or 6-31G** 6 probably
More informationTheoretical determination of the heat of formation of methylene
Theoretical determination of the heat of formation of methylene Nikos L. Doltsinis and Peter J. Knowles School of Chemistry, University of Birmingham, Edgbaston, Birmingham B5 2TT, United Kingdom The heat
More informationCoupled-cluster theory for medium-mass nuclei
Coupled-cluster theory for medium-mass nuclei Thomas Papenbrock and G. Hagen (ORNL) D. J. Dean (ORNL) M. Hjorth-Jensen (Oslo) A. Nogga (Juelich) A. Schwenk (TRIUMF) P. Piecuch (MSU) M. Wloch (MSU) Seattle,
More informationAssessment of the accuracy of coupled cluster. perturbation theory for open-shell systems. I. Triples expansions
Assessment of the accuracy of coupled cluster perturbation theory for open-shell systems. I. Triples expansions Janus J. Eriksen,, Devin A. Matthews, Poul Jørgensen, and Jürgen Gauss qleap Center for Theoretical
More informationSupporting Information
Supporting Information Computational Evidence of Inversion of 1 L a and 1 L b -Derived Excited States in Naphthalene Excimer Formation from ab Initio Multireference Theory with Large Active Space: DMRG-CASPT2
More informationRenormalization of Tensor- Network States Tao Xiang
Renormalization of Tensor- Network States Tao Xiang Institute of Physics/Institute of Theoretical Physics Chinese Academy of Sciences txiang@iphy.ac.cn Physical Background: characteristic energy scales
More informationQuantum 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 informationUNEDF- PACK FOREST. RODNEY J. BARTLETT Quantum Theory Project, Departments of Chemistry and Physics University of Florida, Gainesville, Florida USA
UNEDF- PACK FOREST Ab-initio dft: The seamless connection with wave-function theory June 25, 2008 RODNEY J. BARTLETT Quantum Theory Project, Departments of Chemistry and Physics University of Florida,
More informationChemistry Publications
Chemistry Publications Chemistry 2007 Accurate Ab Initio Potential Energy Curve of F2. I. Nonrelativistic Full Valence Configuration Interaction Energies Using the Correlation Energy Extrapolation by Intrinsic
More information( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r)
Born Oppenheimer Approximation: Ĥ el ( R)Ψ el ( r;r) = E el ( R)Ψ el ( r;r) For a molecule with N electrons and M nuclei: Ĥ el What is E el (R)? s* potential surface Reaction Barrier Unstable intermediate
More informationMolecular Simulation I
Molecular Simulation I Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences
More informationPerformance 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 informationQUANTUM CHEMISTRY PROJECT 3: ATOMIC AND MOLECULAR STRUCTURE
Chemistry 460 Fall 2017 Dr. Jean M. Standard November 1, 2017 QUANTUM CHEMISTRY PROJECT 3: ATOMIC AND MOLECULAR STRUCTURE OUTLINE In this project, you will carry out quantum mechanical calculations of
More informationQUANTUM 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 informationFine Structure Calculations of Atomic Data for Ar XVI
Journal of Modern Physics, 2015, 6, 1609-1630 Published Online September 2015 in SciRes. http://www.scirp.org/journal/jmp http://dx.doi.org/10.4236/jmp.2015.611163 Fine Structure Calculations of Atomic
More informationSelf-consistent intermediate Hamiltonians: A coupled cluster type formulation of the singles and doubles configuration interaction matrix dressing
Self-consistent intermediate Hamiltonians: A coupled cluster type formulation of the singles and doubles configuration interaction matrix dressing I. Nebot-Gil and J. Sánchez-Marín Departament de Química
More informationIntroduction 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 informationImporting ab-initio theory into DFT: Some applications of the Lieb variation principle
Importing ab-initio theory into DFT: Some applications of the Lieb variation principle Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department
More informationJack 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 informationElectronic structure theory: Fundamentals to frontiers. VI. Analysis and more.
Electronic structure theory: Fundamentals to frontiers. VI. Analysis and more. MARTIN HEAD-GORDON Department of Chemistry, University of California, Berkeley, and, Chemical Sciences Division, Lawrence
More informationOrbital dependent correlation potentials in ab initio density functional theory
Orbital dependent correlation potentials in ab initio density functional theory noniterative - one step - calculations Ireneusz Grabowski Institute of Physics Nicolaus Copernicus University Toruń, Poland
More informationComputational 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 informationIntroduction to Hartree-Fock Molecular Orbital Theory
Introduction to Hartree-Fock Molecular Orbital Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Origins of Mathematical Modeling in Chemistry Plato (ca. 428-347
More informationElectric properties of molecules
Electric properties of molecules For a molecule in a uniform electric fielde the Hamiltonian has the form: Ĥ(E) = Ĥ + E ˆµ x where we assume that the field is directed along the x axis and ˆµ x is the
More informationLocal Approaches to the Simulation of Electron Correlation in complex systems
Local Approaches to the Simulation of Electron Correlation in complex systems Martin Schütz Institut für Physikalische und Theoretische Chemie, Universität Regensburg Universitätsstraße 31, D-93040 Regensburg
More informationComputational 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 informationAb initio asymptotic-expansion coefficients for pair energies in Møller-Plesset perturbation theory for atoms
Ab initio asymptotic-expansion coefficients for pair energies in Møller-Plesset perturbation theory for atoms K. JANKOWSKI a, R. SŁUPSKI a, and J. R. FLORES b a Nicholas Copernicus University 87-100 Toruń,
More informationManuel Díaz-Tinoco and J. V. Ortiz Department of Chemistry and Biochemistry Auburn University Auburn AL Abstract
JCP Comment on Are polynuclear superhalogens without halogen atoms probable? A high level ab initio case study on triple bridged binuclear anions with cyanide ligands [J. Chem. Phys. 140, 094301 (2014)]
More informationRenner-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 informationSize-extensive wave functions for QMC A linear-scaling GVB approach
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
More informationQMC 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 informationComputational Material Science Part II. Ito Chao ( ) Institute of Chemistry Academia Sinica
Computational Material Science Part II Ito Chao ( ) Institute of Chemistry Academia Sinica Ab Initio Implementations of Hartree-Fock Molecular Orbital Theory Fundamental assumption of HF theory: each electron
More informationToward Coupled-Cluster Implementations in Nuclear Structure
Toward Coupled-Cluster Implementations in Nuclear Structure D.J. Dean and M. Hjorth-Jensen 2 Physics Division, Oak Ridge National Laboratory P.O. Box 2008, Oak Ridge, TN 3783-6373 USA 2 Department of Physics,
More informationThe Nuclear Many Body Problem Lecture 4. Coupled Cluster theory and its application to medium sized nuclei.
The Nuclear Many Body Problem Lecture 4 Coupled Cluster theory and its application to medium sized nuclei. Extending the Ab Initio program beyond the lightest nuclei. Need a theory which scales softly
More informationElectronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
More informationuse the backs of pages as needed
CHEMISTRY 4021/8021 Q1) Propose a simple, united-atom molecular mechanics force-field needed to generate a potential energy surface for an isolated molecule of acetone (Me 2 CO). I.e., provide an energy
More informationAccurate description of potential energy surfaces by ab initio methods : a review and application to ozone
Accurate description of potential energy surfaces by ab initio methods : a review and application to ozone Péter G. Szalay Laboratory of Theoretical Chemistry Institute of Chemistry Eötvös Loránd University,
More informationLecture 4: Hartree-Fock Theory
Lecture 4: Hartree-Fock Theory One determinant to rule them all, One determinant to find them, One determinant to bring them all and in the darkness bind them Second quantization rehearsal The formalism
More informationDensity functional theory predictions of anharmonicity and spectroscopic constants for diatomic molecules
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 6 8 AUGUST 2001 Density functional theory predictions of anharmonicity and spectroscopic constants for diatomic molecules Mutasem Omar Sinnokrot and C. David
More informationAuxiliary-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 informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
More information