Wave function methods for the electronic Schrödinger equation
|
|
- Virgil Harrison
- 6 years ago
- Views:
Transcription
1 Wave function methods for the electronic Schrödinger equation Zürich 2008
2 DFG Reseach Center Matheon: Mathematics in Key Technologies A7: Numerical Discretization Methods in Quantum Chemistry DFG Priority Program: Modern and Universal First-Principle-Methods for Many-Electron Systems in Chemistry and Physics DFG Priority Program: Extraction of Essential Information from Complex Systems
3 Basic model - electronic Schrödinger equation Electronic Schrödinger equation N nonrelativistic electrons + Born Oppenheimer approximation HΨ = EΨ The Hamilton operator H = 1 2 i i N i K ν=1 Z ν x i a ν N i,j 1 x i x j acts on anti-symmetric wave functions Ψ H 1 (R 3 {± 1 2 })N, Ψ(x 1, s 1,..., x N, s N ) R, x i = (x i, s i ) R 3 {± 1 2 }.
4 Goals To go beyond the accuracy of Density Functional Models one tries to approximate the eigenfuntion Ψ, of HΨ = EΨ. high precision: Correlation energy. E corr := E HF E 0 dynamic correlation: (closed shell R) HF is a good approximation, but the resolution of the electron- electron cusp hempers good convergence static correlation: a single Slater deteminant cannot provide a sufficiently good reference approximation, due e.g. open shells, requires multi configurational models excited states degenerate ground states
5 CI Configuration Interaction Method Approximation space for (spin) orbitals (x j, s j ) ϕ(x j, s j ) X h := span {ϕ i : i = 1,..., N }, ϕ i, ϕ j = δ i,j E.g. Canonical orbitals ϕ i, i = 1,..., N are eigenfunctions of the discretized Fock operator F := F h = N k=1 F k : V FCI V FCI Fϕ i λ i ϕ i, φ h = 0 φ h X h Full CI (for benchmark computations N = 18) is a Galerkin method w.r.t. the subspace V FCI = N X h = span{ψ SL = Ψ[ν 1,..ν N ] = 1 N! det(ϕ ν i (x j, s j )) N i,j=1, } i=1
6 CI Configuration Interaction Method N (discrete) eigenfunction ϕ i, Fϕ i λ i ϕ i, φ h = 0 φ h X h The first N eigenfunctions ϕ i are called occupied orbitals the others are called unoccupied orbitals (traditionally) ϕ 1,..., ϕ N, ϕ N+1,..., ϕ N Galerkin ansatz: Ψ = c 0 Ψ 0 + ν J c νψ ν H = ( Ψ ν, HΨ ν ), Hc = Ec, dim V h = N N Theorem (Brillouins theorem) Let Ψ 0 = Ψ HF and Ψ a i = Ψ S then Ψ a i, HΨ 0 = 0
7 Configuration Interaction Method Assumption: For any ϕ H 1 there exist ϕ h X h such that ϕ h ϕ H 1 0, if h 0 (roughly: lim h 0 X h = H 1 (R 3 {± 1 2 })) Theorem Let E 0 be a single eigenvalue and HΨ = E 0 Ψ and E 0,h, Ψ h V h V FCI be the Galerkin solution, then for h < h 0 hold Ψ Ψ h V c inf φ h V h Ψ φ h V E 0,h E 0 C inf φ h V h Ψ φ h 2 V. Since dim V h = O(N N ), (curse of dimension), the full CI method is infeasible for large N or N!!!!
8 Multi-Configuration Self Consistent Field Method - MCSCF Method Optimization of the orbital basis functions ϕ i, N M << N Y h := span {ϕ i : i = 1,..., M} X h := span {ϕ i : i = 1,..., N }, N V MCSCF = Y h = span{ψ SL = Ψ[ν 1,..ν N ] : ϕ ν Y h } i=1 For Ψ = c 0 Ψ 0 + ν J c νψ ν and E = Ψ, HΨ = J MCSCF (c, Φ) E 0 min{j MCSCF (c, Φ ) : c = 1, Φ S M } unknowns c, Φ := (ϕ i ) i=1,...,m S V M on Stiefel manifold S Exisitence results: Friesecke, Lewin(03)
9 Conclusions MCSCF can be solved with 2nd order optimization methods, e.g. trust region Newton methods Instead of FCI one can use a Complete Active Space MCSCF is for multi-configurational problems, ( where RHF is rather bad) static correlation Due to the e-e cusp, M ɛ 1/2 is large! FCI and MCSCF methods are scaling exponentially with N O( M N ) = O(e N )!!! restriction to single-double excitations etc. are not size consistent (see below)
10 Rayleigh-Schrödinger perturbation theory H = H 0 + U, Solution for H 0 is known H 0 ψ 0 = E 0 ψ 0. Adiabatic perturbation H(λ) = H 0 + λu. Assumption: E 0 is a simple eigenvalue of H 0 and λ E(λ), λ ψ(λ) are analytic in {λ C : λ 1 + ɛ} for some ɛ > 0. E(λ) = E k λ k, ψ(λ) = ψ k λ k k=0 k=0
11 Rayleigh-Schrödinger perturbation theory Intermediate normalization ψ = ψ(1) ψ, ψ 0 = 1 ψ l, ψ 0 = 0, l > 0. The projection P 0 onto the orthogonal complement of ψ 0 is defined by P 0 u = I ψ 0, u ψ 0. Inserting the ansatz into the Schrödinger equation k (H 0 E 0 )ψ k λ k = λuψ k λ k + ( E l ψ k l )λ k. k=0 k=0 k=1 l=1 for all λ 1 + ɛ, + sorting w.r.t. to powers of λ gives k (H 0 E 0 )ψ k = Uψ k 1 + E l ψ k l, k = 1,.... l=1
12 Rayleigh-Schrödinger perturbation theory Testing against ψ 0 gives 0 = ψ 0, (H 0 E 0 )ψ k = ψ 0, ( Uψ k 1 + ψ 0, Uψ k 1 + E k. k E l ψ k l ) = l=1 E k = ψ 0, Uψ k 1 k ψ k = P 0 (H 0 E 0 ) 1 P 0 ( Uψ k 1 + E l ψ k l ). For the computation of the energy contribution E k we need all ψ l up to l = k 1. l=1
13 Møller-Plesset Perturbation theory In the sequel we restrict ourselves finite dimensional subspace to V h = V FCI! N H = F + U = F i + ( 1 r i r j N (J i K i )) i=1 i<j i=1 where F is the Fock operator and U is called the fluctuation potential. Proposition Let Ψ = Ψ SL = Ψ[ν 1,..., ν N ] be a Slater determinant of canonical orbitals, then N N FΨ = εψ, ε = λ νi, E 0 = λ i. i=1 i=1
14 Let Ψ 0 be the HF wave function, the first order energy contribution is E 1 = Ψ 0, UΨ 0 = 1 N ij ij 2 i,j=1 Hartree-Fock energy E MP1 = E 0 + E 1 = E HF. Ψ 1 = P 0 (F E 0 ) 1 P 0 UΨ 0 and the second order contribution to the energy MP2: E 2 = Ψ 0, UΨ 1 = Ψ 0, UP 0 (F E 0 ) 1 P 0 UΨ 0, E MP2 = E 0 + E 1 + E 2 = E HF + E 2. Higher order contributions can be computed, e.g. MP3, MP4 etc. too. The convergence of the expansion is not guaranteed.
15 Second quantization Second quantization: annihilation operators: a j Ψ[j, 1,..., N] = Ψ[1,..., N] and = 0 if j not apparent in Ψ[...]. sign-normalization: j appears in the first place in Ψ[j, 1,..., N]. The adjoint of a b is a creation operator a b a b Ψ[1,..., N] = Ψ[b, 1,..., N] = ( 1)N Ψ[1,..., N, b] Lemma a k a l = a l a k, a k a l = a l a k, a k a l + a l a k = δ k.l
16 Excitation operators Single excitation operator e.g. X k 1 ( 1) p Ψ k 1 = Ψ[k, 2,..., N] = X k 1 Ψ 0 = X k j Ψ[1,...,..., N] = a k a 1Ψ 0 higher excitation operator X µ := X b 1,...,b k l 1,...,l k = k i=1 X b i l i, 1 l i < l i+1 N, N < b i < b i+1. A CI solution Ψ = c 0 Ψ 0 + µ J c µψ µ can be written by Ψ = c 0 + c µ X µ Ψ 0 = (I + T )Ψ 0 T = T 1 + T 2 + T µ J
17 Coupled Cluster Method - Exponential-ansatz Theorem (S. 06) Let Ψ 0 be a reference Slater determinant, e.g. Ψ 0 = Ψ HF and Ψ V FCI, or V, satisfying Ψ, Ψ 0 = 1 intermediate normalization. Then there exists an excitation operator (T 1 - single-, T 2 - double-,... excitation operators) N T = T i = t µ X µ i=1 µ J such that Ψ = e T Ψ 0
18 Baker-Campell-Hausdorff expansion We recall the Baker-Campell-Hausdorff formula e T Ae T = A + [A, T ] + 1 2! [[A, T ], T ] + 1 [[[A, T ], T ], T ] +... = 3! A + k=1 1 k! [A, T ] k. For Ψ V h the above series terminates, exercise** e T He T = H+[H, T ]+ 1 2! [[H, T ], T ]+ 1 3! [[[H, T ], T ], T ]+ 1 4! [H, T ] 4 e.g. for a single particle operator e.g. F there holds e T Fe T = F + [F, T ] + [[F, T ], T ]
19 Coupled Cluster energy Ψ V or V FCI = V FCI = span{ψ ν : ν J } Φ, (H E 0 )Ψ = 0 Φ V, V FCI due to Slater Condon rules and normalization Ψ, Ψ 0 = 1 E = Ψ 0, HΨ = E Ψ 0, H(I + T T )Ψ 0 = Ψ 0, H(I + T 1 +T 2 + T T )Ψ 0 Proposition E = Ψ 0, HΨ = Ψ 0, H(I + T 1 +T T 1 2 )Ψ 0
20 Projected Coupled Cluster Method amplitude equations Ψ µ, e T He T Ψ 0 = Ψ µ, e T HΨ = E Ψ µ, e T Ψ = 0 µ J The Projected Coupled Cluster Method consists in the ansatz l T = T k = t µ X µ µ J h k=1, 0 µ J h J, i.e. Ψ µ V h V FCI e.g. CCSD T = T 1 + T 2 = T (t) satisfying 0 = Ψ µ, e T He T Ψ 0 =: f µ (t), t = (t ν ) ν Jh, µ, ν J h These are L = J h << dimv FCI nonlinear equations for L unknown excitation amplitudes t µ.
21 Iteration method to solve CC amplitude equations Quasi-Newton method to solve CC amplitude equations We decompose the (discretized) Hamiltonian H = F + U, F - Fock operator, U - fluctuation potential. There holds [F, X µ ] = [F, X a 1,...,a k l 1,...,l k ] = ( k (λ aj λ lj ))X µ =: ɛ µ X µ. j=1 and [[F, X µ ], X µ ] = 0 together with (Bach-Lieb-Solojev) ɛ µ λ N+1 λ N > 0
22 Iteration method to solve CC amplitude equations The amplitude function t f(t) = (f µ (t)) µ Jh must be zero f µ (t) = Ψ µ, e T He T Ψ 0 = Ψ µ, [F, X µ ]Ψ 0 + Ψ µ, [U, T ]Ψ 0 = 0. The nonlinear amplitude equation f(t) = 0 is solved by Algorithm (quasi Newton-scheme) 1 Choose t 0, e.g. t 0 = 0. 2 Compute t n+1 = t n A 1 f(t n ), where A = diag (ɛ µ ) µ J > 0. The Coupled Cluster Method is size consistent!: H AB = H A +H B, e (T A+T B ) (H A +H B )e Ta+T B = e T A H A e T A +e T B H B e T B
23 Analysis of the Coupled Cluster Method We consider the projected CC as an approximation of the full CI solution! If h 0, then M and max ɛ µ! We need estimates uniformly w.r.t. h, N Definition Let M := dim V FCI dimensional parameter space V = R M equipped with the norm t 2 V := ɛ µ t µ Ψ µ 2 L 2 ((R 3 {± 1 2 })N ) = ɛ µ t µ 2 µ J µ J
24 Analysis of the Coupled Cluster Method Lemma (S.06) There holds t V T Ψ 0 H 1 ((R 3 {± 1 2 })N ) T Ψ 0 V. Lemma (S.06) For t l 2 (J ), the operator T := ν J t νx ν maps N T Ψ L2 t l2 Ψ L 2 Ψ V FCI L 2 (R 3 {± 1 2 }) i=1
25 Analysis of the Coupled Cluster Method Lemma (S.06) For t V, the operator T := ν J t νx ν maps T Ψ H 1 t V Ψ H 1 Ψ V FCI Corollary (S06) The function f : V V is differentiable at t V with the Frechet derivative f [t] : V V given by (f [t]) ν,µ = Ψ ν, e T [H, X µ ]e T Ψ 0 = ɛ ν δ ν,µ + Ψ ν, e T [U, X µ ]e T Ψ 0 All Frechet derivatives t f (k) [t] : V V, are Lipschitz continuous. In particular f (5) 0.
26 Convergence of the Coupled Cluster Method A function f : is called strictly monotone at t if f(t) f(t ), (t t γ t t 2 V for some γ > 0 and all t t V < δ. Theorem (S.06 (a priori estimate)) If f(t) = 0 and f is strictly monotone at t, then t h, resp. Ψ h = e T hψ 0 satisfy t t h V inf t v h V. v R J h P µ J v Ψ Ψ h H 1 inf Ψ e h µx µ Ψ 0 H 1. v R L
27 Convergence of the CC Method - Duality estimates Energy functional J(t) := E(t) := Ψ 0, H(1 + T T 1 2 )Ψ 0 where t solves the amplitude equations (f(t)) ν = Ψ ν, e T He T Ψ 0 = 0, ν J. Let us further consider the Lagrange functional L(t, a) := J(t) f(t), a, t V, a V. and its stationary points L t [t, a](r, b) := J [t]r f [t]r, a = 0, for all a V. and L a [t, a](r, b) = f(t), b = 0 for all b V.
28 Duality estimates for the Coupled Cluster Method The dual solution a satisfies f [t] a = (J [t]) V. Its Galerkin approximation is given by a h = (a µ ) µ Jh V h, f [t h ] a h, v h = (J [t h ]), v h, v h V h. The discrete primal solution t h = (t ν ) ν Jh V h solves f[t h ], b h = 0 for all b h J h. We define the corresponding residual r, r V (f µ ), µ J h (r(t h )) µ = 0, µ J h together Reinhold withschneider, the dual MATHEON residual TU Berlin
29 Duality estimates for the Coupled Cluster Method The derivatives J [t] and f [t] can be explicitly computed Ψ 0, HX µ e T Ψ 0 µ J 1 single (J [t]) = Ψ 0, UΨ µ, µ J 2 double 0, otherwise where T = t ν X ν, and (f [t]) µ,ν = ɛ µ δ µ,ν + Ψ µ, e T [U, X ν ]e T Ψ 0, µ, ν J. Lemma (dual weighted residual, Rannacher) Let x := (t, a) V V, x h := (t h, a h ) V h V h and e h = x x h. L (x) = 0, L (x h )y h = 0 y h V h V h Then L(x) L(x h ) = L [x h ](x y h ) + R 3, y h V h V h
30 Convergence of the Coupled Cluster Energies Theorem (S. 06 a priori estimate) The error in the energy J(t) J(t h ) can be estimated by E E h t t h V a a h V + ( t t h V ) 2 inf t u h V inf a b h V + u h V h b h V +( inf u h V h t u h V ) 2. E E h t t h V a a h V inf t u h V inf a b h V u h V h b h V All constants involved above are uniform w.r.t. N.
31 Conclusions Improvement by adding the e-e-cusp singularity explicitely, r 1,2, f 1,2 methods (Kutzelnigg-Klopper... ) CCSD and CCSD(T) are standard CCSDT; CCSDTQ etc. only for extremely accurate computations CC is the most powerful tool for computing dynamical correlation not good for multi-configurational problems, ( where RHF is rather bad) How to do Multi Reference Coupled Cluster???
Direct Minimization in Density Functional Theory
Direct Minimization in Density Functional Theory FOCM Hongkong 2008 Partners joint paper with: J. Blauert, T. Rohwedder (TU Berlin), A. Neelov (U Basel) joint EU NEST project BigDFT together with Dr. Thierry
More informationI ) Direct minimization in Hartree Fock and DFT - II) Converence o
I ) Direct minimization in Hartree Fock and DFT - II) Converence of Coupled Cluster Computational methods in many electron systems : MP2 Workshop Göteborg 2009 Ab initio computation Reliable computation
More informationTensorization of Electronic Schrödinger Equation and Second Quantization
Tensorization of Electronic Schrödinger Equation and Second Quantization R. Schneider (TUB Matheon) Mariapfarr, 2014 Introduction Electronic Schrödinger Equation Quantum mechanics Goal: Calculation of
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 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 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 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 informationMathematical representations of quantum states
Mathematical representations of quantum states Eric CANCES Ecole des Ponts and INRIA, Paris, France IPAM, Los Angeles, Septembre 13-16, 2016 Outline of the tutorial 2 1. Quantum states of many-particle
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 informationReinhold Schneider, MATHEON TU Berlin. Direct Minimization for effective single particle models (DFT)
Direct Minimization for effective single particle models (DFT) Partners joint work with: J. Blauert, T. Rohwedder (TU Berlin), A. Neelov (U Basel) joint EU NEST project BigDFT together with Dr. Thierry
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 information0 belonging to the unperturbed Hamiltonian H 0 are known
Time Independent Perturbation Theory D Perturbation theory is used in two qualitatively different contexts in quantum chemistry. It allows one to estimate (because perturbation theory is usually employed
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 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 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 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 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 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 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 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 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 Hartree-Fock approximation
Contents The Born-Oppenheimer approximation Literature Quantum mechanics 2 - Lecture 7 November 21, 2012 Contents The Born-Oppenheimer approximation Literature 1 The Born-Oppenheimer approximation 2 3
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 informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Chemistry 3502/4502 Final Exam Part I May 14, 2005 1. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle (e) The
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle
More information1 Time reversal. 1.1 Without spin. Time-dependent Schrödinger equation: 2m + V (r) ψ (r, t) (7) Local time-reversal transformation, T :
1 Time reversal 1.1 Without spin Time-dependent Schrödinger equation: i t ψ (r, t = Local time-reversal transformation, T : Transformed Schrödinger equation d (f T (t dt ] 2m + V (r ψ (r, t (1 t 1 < t
More informationLecture 5. Hartree-Fock Theory. WS2010/11: Introduction to Nuclear and Particle Physics
Lecture 5 Hartree-Fock Theory WS2010/11: Introduction to Nuclear and Particle Physics Particle-number representation: General formalism The simplest starting point for a many-body state is a system 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 informationAttempts at relativistic QM
Attempts at relativistic QM based on S-1 A proper description of particle physics should incorporate both quantum mechanics and special relativity. However historically combining quantum mechanics and
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 informationLecture 3: Quantum Satis*
Lecture 3: Quantum Satis* Last remarks about many-electron quantum mechanics. Everything re-quantized! * As much as needed, enough. Electron correlation Pauli principle Fermi correlation Correlation energy
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 information2~: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 informationMean-Field Theory: HF and BCS
Mean-Field Theory: HF and BCS Erik Koch Institute for Advanced Simulation, Jülich N (x) = p N! ' (x ) ' 2 (x ) ' N (x ) ' (x 2 ) ' 2 (x 2 ) ' N (x 2 )........ ' (x N ) ' 2 (x N ) ' N (x N ) Slater determinant
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 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 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 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 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 informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
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 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 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 informationThe Hartree-Fock Method
1/41 18. December, 2015 Problem Statement First Quantization Second Quantization References Overview 1. Problem Statement 2. First Quantization 3. Second Quantization 4. 5. References 2/41 Problem Statement
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 informationTime reversible Born Oppenheimer molecular dynamics
Time reversible Born Oppenheimer molecular dynamics Jianfeng Lu Mathematics Department Department of Physics Duke University jianfeng@math.duke.edu KI-Net Conference, CSCAMM, University of Maryland, May
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 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 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 informationAb-initio studies of the adiabatic connection in density-functional theory
Ab-initio studies of the adiabatic connection in density-functional theory Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry,
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 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 informationDoes Møller Plesset perturbation theory converge? A look at two electron systems
Does Møller Plesset perturbation theory converge? A look at two electron systems Mark S. Herman George A. Hagedorn Department of Mathematics and Center for Statistical Mechanics, Mathematical Physics,
More informationHartree, Hartree-Fock and post-hf methods
Hartree, Hartree-Fock and post-hf methods MSE697 fall 2015 Nicolas Onofrio School of Materials Engineering DLR 428 Purdue University nonofrio@purdue.edu 1 The curse of dimensionality Let s consider a multi
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 informationwhere P a is a projector to the eigenspace of A corresponding to a. 4. Time evolution of states is governed by the Schrödinger equation
1 Content of the course Quantum Field Theory by M. Srednicki, Part 1. Combining QM and relativity We are going to keep all axioms of QM: 1. states are vectors (or rather rays) in Hilbert space.. observables
More information(1.1) In particular, ψ( q 1, m 1 ; ; q N, m N ) 2 is the probability to find the first particle
Chapter 1 Identical particles 1.1 Distinguishable particles The Hilbert space of N has to be a subspace H = N n=1h n. Observables Ân of the n-th particle are self-adjoint operators of the form 1 1 1 1
More informationNon-relativistic Quantum Electrodynamics
Rigorous Aspects of Relaxation to the Ground State Institut für Analysis, Dynamik und Modellierung October 25, 2010 Overview 1 Definition of the model Second quantization Non-relativistic QED 2 Existence
More informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
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 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 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 informationThe very basic: QF methods and many body perturbation theory
The very basic: QF methods and many body perturbation theory Georg Kresse, Felix Hummel Faculty of Physics, Universität Wien Funded by the Austrian FWF SFB ViCoM FOR 1346 Outline Intro 2nd PT Diag. Introduction
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 informationIntroduction to Electronic Structure Theory
Introduction to Electronic Structure Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology June 2002 Last Revised: June 2003 1 Introduction The purpose of these
More informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
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 informationDRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS
DRIFT OF SPECTRALLY STABLE SHIFTED STATES ON STAR GRAPHS ADILBEK KAIRZHAN, DMITRY E. PELINOVSKY, AND ROY H. GOODMAN Abstract. When the coefficients of the cubic terms match the coefficients in the boundary
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 informationApplications of Matrix Functions Part III: Quantum Chemistry
Applications of Matrix Functions Part III: Quantum Chemistry Emory University Department of Mathematics and Computer Science Atlanta, GA 30322, USA Prologue The main purpose of this lecture is to present
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 informationYingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2
Purdue University CHM 67300 Computational Quantum Chemistry REVIEW Yingwei Wang October 10, 2013 Review: Prof Slipchenko s class, Fall 2013 Contents 1 Hartree energy 2 2 Many-body system 2 3 Born-Oppenheimer
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationPart 5: Penalty and augmented Lagrangian methods for equality constrained optimization. Nick Gould (RAL)
Part 5: Penalty and augmented Lagrangian methods for equality constrained optimization Nick Gould (RAL) x IR n f(x) subject to c(x) = Part C course on continuoue optimization CONSTRAINED MINIMIZATION x
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 informationSimulation Methods II
Simulation Methods II Maria Fyta Institute for Computational Physics Universität Stuttgart Summer Term 2018 SM II - contents First principles methods Hartree-Fock and beyond Density-funtional-theory Ab
More informationConvex Optimization and Modeling
Convex Optimization and Modeling Duality Theory and Optimality Conditions 5th lecture, 12.05.2010 Jun.-Prof. Matthias Hein Program of today/next lecture Lagrangian and duality: the Lagrangian the dual
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 information2 Electronic structure theory
Electronic structure theory. Generalities.. Born-Oppenheimer approximation revisited In Sec..3 (lecture 3) the Born-Oppenheimer approximation was introduced (see also, for instance, [Tannor.]). We are
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 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 informationAlgorithms and Computational Aspects of DFT Calculations
Algorithms and Computational Aspects of DFT Calculations Part I Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational
More informationWe will use C1 symmetry and assume an RHF reference throughout. Index classes to be used include,
1 Common Notation We will use C1 symmetry and assume an RHF reference throughout. Index classes to be used include, µ, ν, λ, σ: Primary basis functions, i.e., atomic orbitals (AOs) {φ µ ( r)}. The size
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 informationGround-State Energies of Coulomb Systems and Reduced One-Particle Density Matrices
Ground-State Energies of Coulomb Systems and Reduced One-Particle Density Matrices Heinz Siedentop Chennai, August 16, 2010 I. Introduction Atomic Schrödinger operator H N,Z := N n=1 self-adjointly realized
More informationCoupled cluster theory: Analytic derivatives, molecular properties, and response theory
Coupled cluster theory: Analytic derivatives, molecular properties, and response theory Wim Klopper and David P. Tew Lehrstuhl für Theoretische Chemie Institut für Physikalische Chemie Universität Karlsruhe
More informationQMI PRELIM Problem 1. All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work.
QMI PRELIM 013 All problems have the same point value. If a problem is divided in parts, each part has equal value. Show all your work. Problem 1 L = r p, p = i h ( ) (a) Show that L z = i h y x ; (cyclic
More informationComputational Chemistry. Ab initio methods seek to solve the Schrödinger equation.
Theory Computational Chemistry Ab initio methods seek to solve the Schrödinger equation. Molecular orbital theory expresses the solution as a linear combination of atomic orbitals. Density functional theory
More information2.1 Green Functions in Quantum Mechanics
Chapter 2 Green Functions and Observables 2.1 Green Functions in Quantum Mechanics We will be interested in studying the properties of the ground state of a quantum mechanical many particle system. We
More informationExercises : Questions
Exercises 18.05.2017: Questions Problem 1 where Calculate the following commutators: a) [ Ĥ, ˆp ], b) [ Ĥ, ˆr ], Ĥ = 1 2m ˆp2 + V ˆr), 1) ˆp 2 = ˆp 2 x + ˆp 2 y + ˆp 2 z and V ˆr) = V ˆx, ŷ, ẑ) is an arbitrary
More informationComputational Chemistry
Computational Chemistry Physical Chemistry Course Autumn 2015 Lecturers: Dos. Vesa Hänninen and Dr Garold Murdachaew vesa.hanninen@helsinki.fi Room B407 http://www.helsinki.fi/kemia/fysikaalinen/opetus/
More informationDiagonalization of bosonic quadratic Hamiltonians by Bogoliubov transformations 1
Diagonalization of bosonic quadratic Hamiltonians by Bogoliubov transformations 1 P. T. Nam 1 M. Napiórkowski 1 J. P. Solovej 2 1 Institute of Science and Technology Austria 2 Department of Mathematics,
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 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 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 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 informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
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 information5. Duality. Lagrangian
5. Duality Convex Optimization Boyd & Vandenberghe Lagrange dual problem weak and strong duality geometric interpretation optimality conditions perturbation and sensitivity analysis examples generalized
More information