Coupled Cluster Theory and Tensor Decomposition Techniques
|
|
- Hilary Dawson
- 6 years ago
- Views:
Transcription
1 Coupled Cluster Theory and Tensor Decomposition Techniques Alexander A. Auer Atomistic Modelling Group Interface Chemistry and Surface Engineering
2 Hierarchy of quantum chemical methods
3 Analysis of orbitals and electron density
4 predicting 13 C nmr chemical shifts deviation from experiment [ppm] + vibr. corr. + large basis CCSD(T) CCSD MP2 DFT HF-SCF benchmark for 16 compounds : A.A.A., J.Gauss, J.F.Stanton J. Chem. Phys., 118, (2003)
5 problem from synthesis / spectroscopy NMR chemical shifts [ppm] : C 1 C 2 C 3 C S.G.D yachkova et al. Russ. Chem. Bull., Int. Ed., 4, 751 (2001)
6 the right testcase NMR chemical shifts [ppm] : C 1 C 2 C 3 C 4 C 1 C 2 C 3 C K.Banert in Science of Synthesis, Houben-Weyl, Methods of Molecular Transformations, Bd. 24, S. 1059, Thieme Verlag Stuttgart (2006)
7 benchmarks geometry: MP2 / pvtz - NMR: B3-LYP / TZVP C 1 C 2 C 2 C 3 C 4 C 3 C [ppm]
8 benchmarks geometry: MP2 / pvtz - NMR: B3-LYP / TZVP C 1 C 2 C 2 C 3 C 4 C 3 C geometry: CCSD(T) [ppm] / pvtz - NMR: MP2 / tz2p C 1 C 2 C 2 C 3 C 4 C 4 C
9 benchmarks geometry: MP2 / pvtz - NMR: B3-LYP / TZVP C 1 C 2 C 2 C 3 C 4 C 3 C geometry: CCSD(T) [ppm] / pvtz - NMR: MP2 / tz2p C 1 C 2 C 2 C 3 C 4 C 4 C geometry: CCSD(T) / pvtz - NMR: CCSD(T) / tz2p C 1 C 2 C 2 C 3 C 3 C
10 ZPV and solvent effects solvent effects: chem. shift [ppm] ε = 1 ε = 2 ε = 20 ε = 78 C C experimental observation: change in chemical shifts for EtOCCH in different solvents at most 1 ppm benchmark for ZVP effects for 13 C shifts: in the order of 1-2 ppm
11 CCSD(T) prediction for azide compound theory geometry: MP2 / pvtz experiment NMR: CCSD(T) / qz2p C 1 C 2 C 2 C 3 C 4 C 3 C [ppm]
12 E. Prochnow, A. A. Auer and K. Banert, J. Phys. Chem. A, 111, 9945 (2007) CCSD(T) prediction for azide compound theory geometry: MP2 / pvtz experiment NMR: CCSD(T) / qz2p C 1 C 2 C 2 C 3 C 4 C 3 C [ppm]
13 ab-initio methods : configuration interaction exact wavefunction is a many electron function!
14 ab-initio methods : configuration interaction exact wavefunction is a many electron function! idea : expand exact wavefunction in complete basis of many electron functions
15 ab-initio methods : configuration interaction exact wavefunction is a many electron function! idea : expand exact wavefunction in complete basis of many electron functions many electron function : HF-SCF slater determinant complete basis of many electron functions : all possible determinants from occupied and virtuals Configuration Interaction (CI) : linear expansion, one parameter per determinant, optimized variationally
16 ab-initio methods : configuration interaction wavefunction expression as linear CI expansion : Ψ CI = Ψ 0 + ci a Ψ a i + 1 ia 4 ijab c ab ij Ψ ab ij +...
17 ab-initio methods : configuration interaction wavefunction expression as linear CI expansion : Ψ CI = Ψ 0 + ci a Ψ a i + 1 ia 4 ijab c ab ij Ψ ab ij +... Ψ CI = Ψ 0 + Ĉ 1 Ψ 0 + Ĉ 2 Ψ diagonalization of Hamiltonian matrix yields energies for ground and excited states
18 ab-initio methods : configuration interaction use all determinants that can be constructed expansion in complete basis (number of determinants grows exponentially!!)
19 ab-initio methods : configuration interaction use all determinants that can be constructed expansion in complete basis (number of determinants grows exponentially!!) exact solution of electronic problem in given AO basis (Full Configuration Interaction FCI)
20 ab-initio methods : configuration interaction use all determinants that can be constructed expansion in complete basis (number of determinants grows exponentially!!) exact solution of electronic problem in given AO basis (Full Configuration Interaction FCI) in practice only truncated expansion feasible : CIS : only singly substituted determinants CISD : singly and double substituted determinants CISDT : single, double and triple substitutions...
21 configuration interaction and size constancy But : extensivity of energy not given!! (CI is not size-extensive) unphysical description leads to wrong asymptotic of energy for extended systems the larger the system, the smaller the fraction of correlation energy from a truncated CI
22 configuration interaction and size constancy But : extensivity of energy not given!! (CI is not size-extensive) unphysical description leads to wrong asymptotic of energy for extended systems the larger the system, the smaller the fraction of correlation energy from a truncated CI recommendation : if possible don t use CI!
23 coupled cluster theory ansatz : instead of linear expansion, use exponential expansion! Ψ CC = e ˆT Ψ 0
24 coupled cluster theory ansatz : instead of linear expansion, use exponential expansion! Ψ CC = e ˆT Ψ 0 cluster operator ˆT analogous to CI-Operator ˆT = ˆT 1 + ˆT 2 + ˆT
25 coupled cluster theory ansatz : instead of linear expansion, use exponential expansion! Ψ CC = e ˆT Ψ 0 cluster operator ˆT analogous to CI-Operator coefficients are called amplitudes: ˆT = ˆT 1 + ˆT 2 + ˆT e ˆT = 1 + ˆT + 1 2! ˆT ! ˆT
26 coupled cluster theory product terms resulting from exponential expansion Ψ CC = (1 + ˆT ! ˆT ! ˆT ˆT ! ˆT ! ˆT ˆT 1 ˆT ! ˆT 1 2 ˆT ! ˆT 1 ˆT ) Ψ 0
27 coupled cluster theory Energy equations : Amplitude equations : Ψ 0 e ˆT Ĥe ˆT Ψ 0 = E Ψ ab.. ij.. e ˆT Ĥe ˆT Ψ 0 = 0 energy as expectation value of similarity transformed Hamiltonian projection trick to derive amplitude equations iterative solution of nonlinear set of equations! state selective (ground state) quite expensive!
28 coupled cluster theory CI based method with inclusion of higher excitations strictly size consistent faster convergence towards FCI highly accurate! approximate methods : CCSD : Coupled Cluster Singles and Doubles CCSD(T):Singles, Doubles + perturbative approximation to Triples CCSDT : Singles, Doubles and Triples...
29 CCD T 2 equations T ab ij v ab ij fi m tmj ab f b teb ij t fa mnvef mn e tij ea tab tef ij t ab mnvij mn mnvef mn + t ea nj vie nb 1 2 tab mit ef nj vef mn tef ij vef ab tmit eb nj fa vef mn
30 typical problem: 50 occ, 400 virt T ab ij amplitudes (3 GB) W(abef) (200 GB) (less if stored in AO-Basis) T abc ijk (1.5 TB) one a<bc block (250 MB) memory requirements 10 GB typical nr. of iterations : typical calculation times : days, weeks (.. months) molecular properties: factor of >2 more expensive
31 parallel CCSD(T) implementation geometry optimization of Ferrocene (staggered and eclipsed), ae-ccsd(t)/cc-pwcvtz, 96 e, 672 Bf., 9 days per step on 14 nodes spin-rotation tensor of HC 7 N, ae-ccsd(t)/cc-pvqz, 470 Bf., 4 days on 7 nodes adamantyl cation (C 10 H + 15 ), ae-ccsd(t)/cc-pvtz, 510 Bf., 5 days per step on 9 nodes M. Harding, T. Metzroth, A. A. Auer and J. Gauss J. Chem. Theory Comput. 4, 64 (2008)
32 hierarchy of post-hf ab-initio methods MP2 / CC2 : non-iterative / iterative N 5 CCSD : iterative NoccN 2 vrt 4 (formal N 6 ) CC3, CCSD(T) : iterative / non-iterative NoccN 3 vrt 4 (formal N 7 ) CCSDT : iterative NoccN 3 vrt 5 (formal N 8 ) CCSDTQ : iterative NoccN 4 vrt 6 (formal N 10 )... T.U. Helgaker, W. Klopper and D.P. Tew Mol. Phys., 106, 2107 (2008)
33 task for future methods
34 tensor decomposition in quantum chemistry related techniques have been used for DFT since 70ies: resolution of the identity / density fitting Laplace-MP2 Cholesky decomposition various applications of SVD techniques
35 tensor decomposition SVD for n-dimensions is not well-defined but: any n-dimensional tensor can be approximated by decomposition find convenient format for decomposition construct trivial decomposition and then apply rank reduction rank reduction : algorithm for finding a low-rank approximation with control of error - used as back box here (one threshold!)
36 tensor decomposition format chosen for application in Coupled Cluster Theory: canonical format X ab ij k l =1 (x a ) l (x b ) l (x i ) l (x j ) l advantages: - most compact form - leads to most general contraction scheme - best reduction of scaling in contractions task: devise CC algorithm in decomposed format
37 notes on decomposition algorithm convert original tensor to canonical format by trivial decomposition find large entry in original tensor to use as new rank optimize rank to improve A-residual iterative procedure to adjust all ranks to minimize à A. if difference is larger ε, pick up more ranks U. Benedikt, A. A. Auer, M. Espig, W. Hackbusch, J. Chem. Phys., accepted
38 tensor decomposition for chemists decomposition of arbitrary tensor: a 1 a 2. a r br T r a 1 ( b 1 b 2... ) 1 + a 2. = A ( b 1 b 2... ) = 2 x 11 x x 21 x
39 tensor decomposition for chemists decomposition of arbitrary tensor: a 1 a 2. ( b 1 b 2... ) 1 + a 1 a 2. ( b 1 b 2... ) =
40 tensor decomposition for chemists decomposition of arbitrary tensor: ( ) = ( ) ( )
41 tensor decomposition for chemists rank reduction from 3 to 2 : 1 2 ( ) + 3 = ( )
42 notes on decomposition algorithm 1: Choose initial à 0 V R d and parameter ε R >0. Define iteration count k := 0, compute the gradient G 0 := J(à 0 ) and D 0 := G 0. 2: while the gradient G k > ε do 3: Compute the smallest root α k [0, 1] of the polynomial p(α) := { } α k := min α R 0 : p(α) := J(à k + αd k ), D k = 0. 4: Update the representation system of Ã, i.e. à k+1 := à k + α k D k. 5: Compute the gradient for the updated system, i.e. G k+1 := J(à k+1 ). 6: Compute β k := Gk+1 G k,g k+1, γ G k 2 k := max{0, β k }. 7: Update the new search direction, i.e. D k+1 := G k+1 + γ k D k. 8: k k : end while J(à k + αd k ), D k, i.e.
43 improved decomposition algorithm
44 coupled cluster in a nutshell calculation of HF yields MO coefficients C µ,p, energies ε i and integrals like µν σρ transform integrals to MO basis ab ij = µνσρ Cµ a Cν b Ci σ C ρ j µν σρ iteratively solve CC amplitude equations by performing contractions and updating residual obtain energy from amplitudes
45 two electron integrals decomposition of two electron integrals: µν σρ r l=1 quantity of 4 #AOs r benefit if rank r is smaller than #AO 3 (χ µ ) l (χ ν ) l (χ σ ) l (χ ρ ) l
46 integral transformation ab ij = Cµ a Cν b Ci σ C ρ j µν σρ µνσρ = = µνσρ r l=1 r l=1 Cµ a Cν b Ci σ C ρ j ( µ r l=1 C a µ (χ µ ) l ) ( ν (χ µ ) l (χ ν ) l (χ σ ) l (χ ρ ) l (v a ) l (v b ) l (v i ) l (v j ) l ) ( Cν b (χ ν ) l C σ σ i (χ σ ) l )( ρ C ρ j (χρ ) l )
47 integral transformation ab ij = Cµ a Cν b Ci σ C ρ j µν σρ µνσρ = r l=1 (v a ) l (v b ) l (v i ) l (v j ) l note : rank stays the same upon transformation rank is the same irrespective of basis used! computational effort for transformation : 4 #MOs #AOs r benefit if rank is smaller than #AO 3
48 perturbation theory E [2] = 1 4 ij ab 2 abij ε i + ε j ε a ε b use integrals in decomposed form factorizable denominator approximation
49 perturbation theory E [2] = 1 4 ij ab 2 abij ε i + ε j ε a ε b use integrals in decomposed form factorizable denominator approximation 1 ε i + ε j ε a ε b = N n=1 N ω n exp( α n (ε i + ε j ε a ε b )) n=1 ω n exp( α n ε i ) exp( α n ε j ) exp(α n ε a ) exp(α n ε b )
50 post HF-methods in spirit of coupled cluster - rewrite MP2 energy in terms of amplitudes : E [2] = 1 4 tij ab ij ab with tij ab = abij ij ab ε i + ε j ε a ε b if integrals and denominator are decomposed, amplitudes will always automatically be obtained in decomposed form t k l =1 (t a ) l (t b ) l (t i ) l (t j ) l
51 post HF-methods E [2] = 1 4 abij tij ab ij ab MP2 energy can be calculated as inner products ( k ) ( ) E [2] 1 4 r (t a ) l (t b ) l (t i ) l (t j ) l (v a ) l (v b ) l (v i ) l (v j ) l abij l=1 = 1 4 k r l =1 l=1 ( a=1 l =1 ) (t a ) l (v a ) l )( (t b ) l (v b ) l )( (t i ) l (v i ) l )( (t j ) l (v j ) l b=1 i=1 j=1
52 coupled cluster theory amplitude equations Ψ ab.. ij.. e ˆT Ĥe ˆT Ψ 0 = 0 example of T 2 equations for CCD: T ab ij v ab ij teb ij t fa fi m tmj ab f b mnvef mn e tij ea tef ij t ab tab mnvef mn mnvij mn + t ea 1 2 tab mit ef nj vie nb nj vef mn express all (!) quantities in decomposed form tef ij vef ab tmit eb nj fa vef mn
53 coupled cluster theory r ab ij tij ef vef ab = all tensors in decomposed representation t ef ij v ab = k r ef = l =1 l=1e=1 f =1 k r l =1 l=1 ( e=1 tij ef vef ab e=1 f =1 [ ] (t e ) l (t f ) l (t i ) l (t j ) l ][(v a ) l (v b ) l (v e ) l (v f ) l ) ( ) (t e ) l (v e ) l (t f ) l (v f ) l (v a ) l (v b ) l (t i ) l (t j ) l f =1 complexity is reduced from occ 2 virt 4 to k r N benefit if product of ranks is smaller than N 5
54 benchmarks and results decomposition of AO integrals H 2 O - different basis size (STO-3G, 6-31G, cc-pvdz, aug-cc-pvdz, cc-pvtz) initial rank eps = 10-2 eps = 10-4 eps = 10-6 slice = 250 slice = rank number of basis functions (N) scaling of rank with basis functions: n 2.5
55 benchmarks and results decomposition of AO integrals initial ranks final ranks (10-2 ) final ranks (10-4 ) final ranks (10-6 ) LiH-chain (6-31G) rank (LiH) n scaling of rank with system size: N 2
56 benchmarks and results decomposition of MO integrals scaling of v abcd -ranks in C n H 2n 6-31G scaling of v ijkl -ranks in C n H 2n 6-31G initial (n 2.7 ) eps=10-2 (n 1.3 ) eps=10-4 (n 2.1 ) eps=10-6 (n x ) initial (n 2.8 ) eps=10-2 (n 1.8 ) eps=10-4 (n 2.0 ) eps=10-6 (n 2.3 ) rank rank n n
57 benchmarks and results decomposition MP2 amplitudes 2000 initial eps=10-2 eps=10-4 eps=10-6 t(abij) in H 2 O using different basis size 1500 rank number of basis functions (N) scaling of rank with basis functions: n 1.2
58 benchmarks and results decomposition MP2 amplitudes t(abij) in (LiH) n 6-31G 2000 initial eps=10-2 eps=10-4 eps= rank (LiH) n scaling of rank with system size: N 1.2
59 benchmarks and results error in MP2 energy -1.0e e-01 E MP2 for different eps in C n H 2n 6-31G CH 4 C 2 H 6 C 3 H 8 C 4 H e-02 E MP2 [Hartree] -1.0e e e e e eps mh accuracy at thresholds around
60 CC algorithm - first tests test-lccd (not all terms yet...) only use decomposed integrals and amplitudes rank reduction after every contraction apply same threshold for all quantities yet no DIIS error in correlation energy for H 2 O using 6-31G correlation energy ε = 10-2 ε = error [Hartree] energy [Hartree] iteration
61 how compact is decomposed representation? rank does not change upon transformation no guarantee to find lowest possible rank does algorithms know about sparsity? Test : dissociate O 2 and investigate ranks MP2 energy in O 2 using 6 31G ε=10 4 interaction energy energy error rank E or E [µhartree] rank O O distance [Å]
62 improvement for CCSD(T) expected? Terms like : T abc ijk = ˆP can be cast into the form: ( ) R L (t m ) r (v m ) l r=1 m=1 l=1 Tim ab Vjk cm m + e Tij ae Vek bc (t a ) r (t b ) r (v c ) l (t i ) r (v j ) l (v k ) l current estimate : R N 1.5, L N 2.0 complexity : R L orb (with current estimate N 4.5 )
63 summary - tensor decomposition integral storage and transformation scales as r N coupled cluster iterations scale as r ampl. r int. N actual ranks scale as N 1.5 (ampl.) and N 2.0 (MO int.) decomposition scales with initial rank and quickly becomes main bottleneck! next steps: CCSD, CCSD(T) (N 5?) and CCSDT more efficient decomposition algorithm (trivial!) parallelization molecular properties
64 acknowledgments Udo Benedikt (MPIE) DFG : SPP 1145 AU 206/2-1 Eric Prochnow (TUC) AU 206/1-1 Prof. J. Gauss (Uni Mainz) BMBF M. Harding, T. Metzroth Fonds der Chemischen Industrie Prof. K. Banert (TUC) Prof. W. Hackbusch (MPI-MIS) M. Espig TU Chemnitz MPIE CES (RUB) M. Nooijen (Uni Waterloo) S. Hirata (QTP) TCE Group : D. Bernholt, R. Pitzer, P. Saddayappan, G. Baumgartner...
65 thanks for your attention!!
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 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 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 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 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 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 informationCoupled-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 informationHigh Accuracy Local Correlation Methods: Computer Aided Implementation
High Accuracy Local Correlation Methods: Computer Aided Implementation Marcel Nooijen Alexander Auer Princeton University University of Waterloo USA Canada So Hirata, PNNL Supported by: NSF ITR (Information
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 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 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 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 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 informationCyclops Tensor Framework: reducing communication and eliminating load imbalance in massively parallel contractions
Cyclops Tensor Framework: reducing communication and eliminating load imbalance in massively parallel contractions Edgar Solomonik 1, Devin Matthews 3, Jeff Hammond 4, James Demmel 1,2 1 Department of
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 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 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 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 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 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 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 informationCoupled-Cluster Theory. Jürgen Gauss
Coupled-Cluster Theory Jürgen Gauss Electron Correlation, Size Consistency, and Exponential Ansatz Electron Correlation: the Helium Atom nucleus electron 2 correlation effect wave function Ψ exact (r 1,r
More informationIntroduction to Computational Quantum Chemistry: Theory
Introduction to Computational Quantum Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC 3108 Course Lectures 2007 Introduction Hartree Fock Theory Configuration Interaction Lectures 1 Introduction
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 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 informationBeyond Hartree-Fock: MP2 and Coupled-Cluster Methods for Large Systems
John von Neumann Institute for Computing Beyond Hartree-Fock: MP2 and Coupled-Cluster Methods for Large Systems Christof Hättig published in Computational Nanoscience: Do It Yourself, J. Grotendorst, S.
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 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 informationBuilding 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 informationIntroduction to Computational Chemistry: Theory
Introduction to Computational Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC andrew.gilbert@anu.edu.au 3023 Course Lectures Introduction Hartree Fock Theory Basis Sets Lecture 1 1 Introduction
More informationGeminal product wave function ansatz for the description of correlated systems
Geminal product wave function ansatz for the description of correlated systems D. Van Neck, S. De Baerdemacker, P. Bultinck Ghent University P.A. Limacher, P. Johnson, P.W. Ayers McMaster University April
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 informationTDDFT in Chemistry and Biochemistry III
TDDFT in Chemistry and Biochemistry III Dmitrij Rappoport Department of Chemistry and Chemical Biology Harvard University TDDFT Winter School Benasque, January 2010 Dmitrij Rappoport (Harvard U.) TDDFT
More informationMethod development - new algorithms and approximations for post HF Methods
Method development - new algorithms and approximations for post HF Methods The great advantage of high level ab-initio methods like from Coupled Cluster theory is that they can be improved systematically
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 informationMODERN ELECTRONIC STRUCTURE THEORY: Electron Correlation
5.61 Physical Chemistry Lecture #30 1 MODERN ELECTRONIC STRUCTURE THEORY: Electron Correlation In the previous lecture, we covered all the ingredients necessary to choose a good atomic orbital basis set.
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 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 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 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 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 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 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 informationSolution 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 informationHigh-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 informationPost Hartree-Fock: MP2 and RPA in CP2K
Post Hartree-Fock: MP2 and RPA in CP2K A tutorial Jan Wilhelm jan.wilhelm@chem.uzh.ch 4 September 2015 Further reading MP2 and RPA by Mauro Del Ben, Jürg Hutter, Joost VandeVondele Del Ben, M; Hutter,
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 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 informationLazy matrices for contraction-based algorithms
Lazy matrices for contraction-based algorithms Michael F. Herbst michael.herbst@iwr.uni-heidelberg.de https://michael-herbst.com Interdisziplinäres Zentrum für wissenschaftliches Rechnen Ruprecht-Karls-Universität
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 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 informationA Whirlwind Introduction to Coupled Cluster Response Theory. 1 Time-Independent Coupled Cluster Theory
A Whirlwind Introduction to Coupled Cluster Response Theory T. Daniel Crawford, Virginia Tech, Blacksburg, Virginia, U.S.A. 1 Time-Independent Coupled Cluster Theory If the Hamiltonian is independent of
More informationHeterogeneous programming for hybrid CPU-GPU systems: Lessons learned from computational chemistry
Heterogeneous programming for hybrid CPU-GPU systems: Lessons learned from computational chemistry and Eugene DePrince Argonne National Laboratory (LCF and CNM) (Eugene moved to Georgia Tech last week)
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 informationThe calculation of the universal density functional by Lieb maximization
The calculation of the universal density functional by Lieb maximization Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry,
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 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 informationWe also deduced that transformations between Slater determinants are always of the form
.3 Hartree-Fock The Hartree-Fock method assumes that the true N-body ground state wave function can be approximated by a single Slater determinant and minimizes the energy among all possible choices of
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 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 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 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 informationUsing BLIS for tensor computations in Q-Chem
Using BLIS for tensor computations in Q-Chem Evgeny Epifanovsky Q-Chem BLIS Retreat, September 19 20, 2016 Q-Chem is an integrated software suite for modeling the properties of molecular systems from first
More informationMolecular Magnetic Properties. The 11th Sostrup Summer School. Quantum Chemistry and Molecular Properties July 4 16, 2010
1 Molecular Magnetic Properties The 11th Sostrup Summer School Quantum Chemistry and Molecular Properties July 4 16, 2010 Trygve Helgaker Centre for Theoretical and Computational Chemistry, Department
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 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 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 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 informationCyclops Tensor Framework
Cyclops Tensor Framework Edgar Solomonik Department of EECS, Computer Science Division, UC Berkeley March 17, 2014 1 / 29 Edgar Solomonik Cyclops Tensor Framework 1/ 29 Definition of a tensor A rank r
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 informationSoftware implementation of correlated quantum chemistry methods. Exploiting advanced programming tools and new computer architectures
Software implementation of correlated quantum chemistry methods. Exploiting advanced programming tools and new computer architectures Evgeny Epifanovsky Q-Chem Septermber 29, 2015 Acknowledgments Many
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 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 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 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 informationDensity Matrix Renormalization Group Algorithm - Optimization in TT-format
Density Matrix Renormalization Group Algorithm - Optimization in TT-format R Schneider (TU Berlin MPI Leipzig 2010 Density Matrix Renormalisation Group (DMRG DMRG algorithm is promininent algorithm to
More informationConvergence 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 informationLow Rank Bilinear Algorithms for Symmetric Tensor Contractions
Low Rank Bilinear Algorithms for Symmetric Tensor Contractions Edgar Solomonik ETH Zurich SIAM PP, Paris, France April 14, 2016 1 / 14 Fast Algorithms for Symmetric Tensor Contractions Tensor contractions
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 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 informationMolecular Magnetic Properties
Molecular Magnetic 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, Oslo,
More informationMolecular 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 informationThis is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry.
1 Computational Chemistry (Quantum Chemistry) Primer This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry. TABLE OF CONTENTS Methods...1 Basis
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 informationExplicitly correlated wave functions
Explicitly correlated wave functions Wim Klopper and David P. Tew Lehrstuhl für Theoretische Chemie Institut für Physikalische Chemie Universität Karlsruhe (TH) C 4 Tutorial, Zürich, 2 4 October 2006 Outline
More informationLecture 4: methods and terminology, part II
So theory guys have got it made in rooms free of pollution. Instead of problems with the reflux, they have only solutions... In other words, experimentalists will likely die of cancer From working hard,
More informationIntroduction to DFTB. Marcus Elstner. July 28, 2006
Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD
More informationChem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 28, December 08, 2014
Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi Lecture 28, December 08, 2014 Solved Homework Water, H 2 O, involves 2 hydrogen atoms and an oxygen
More informationBasis Set for Molecular Orbital Theory
Basis Set for Molecular Orbital Theory! Different Types of Basis Functions! Different Types of Atom Center Basis Functions! Classifications of Gaussian Basis Sets! Pseudopotentials! Molecular Properties
More informationLec20 Wed 1mar17 update 3mar 10am
564-17 Lec20 Wed 1mar17 update 3mar 10am Figure 15.2 Shows that increasing the diversity of the basis set lowers The HF-SCF energy considerably, but comes nowhere near the exact experimental energy, regardless
More informationDensity-Fitting Approximations to the Electron Repulsion Integrals
Density-Fitting Approximations to the Electron Repulsion Integrals C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Created on 19 July 2010 Here we will follow some
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 informationCoupled Cluster Green function: model involving singles and. doubles excitations
Coupled Cluster Green function: model involving singles and doubles excitations Kiran Bhaskaran-Nair a,b, Karol Kowalski c, William A. Shelton a,b a Cain Department of Chemical Engineering, Louisiana State
More information4 October - 8 October 2010: Report on the article entitled R12 Methods in Explicitly correlated Molecular electronic structure Theory
Ohio University Department of Mathematics 1 Journal of Magagula Vusi Mpendulo: Fall 2010-2011 Magagula Vusi Mpendulo 4 October - 8 October 2010: Report on the article entitled R12 Methods in Explicitly
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
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 informationSCF calculation on HeH +
SCF calculation on HeH + Markus Meuwly Department of Chemistry, University of Basel, Basel, Switzerland Abstract This document describes the main steps involved in carrying out a SCF calculation on the
More informationA practical view on linear algebra tools
A practical view on linear algebra tools Evgeny Epifanovsky University of Southern California University of California, Berkeley Q-Chem September 25, 2014 What is Q-Chem? Established in 1993, first release
More informationTHE DEVELOPMENT OF THE CBS-Wes METHOD
Wesleyan University THE DEVELOPMENT OF THE CBS-Wes METHOD By Duminda Sampath Ranasinghe Faculty Advisor: Prof. George A. Petersson A Dissertation submitted to the Faculty of Wesleyan University in partial
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 information