Time-independent molecular properties
|
|
- Lenard Rogers
- 6 years ago
- Views:
Transcription
1 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, Oslo, Norway European Summer School in Quantum Chemistry (ESQC) 2017 Torre Normanna, Sicily, Italy September 10 23, 2017 Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
2 Sections 1 Introduction 2 Energy Functions 3 Derivatives for Variational Wave Functions 4 Derivatives for Nonvariational Wave Functions 5 Examples of derivatives 6 Derivatives in second quantization 7 Geometrical Properties Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
3 Introduction Section 1 Introduction Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
4 Introduction Table of Contents 1 Introduction 2 Energy Functions Variational and nonvariational wave functions 3 Derivatives for Variational Wave Functions Molecular gradients Molecular Hessians 4 Derivatives for Nonvariational Wave Functions Nonvariational Wave Functions Lagrangian method 5 Examples of derivatives Hartree Fock Molecular Gradient FCI Lagrangian and Derivatives 6 Derivatives in second quantization Second-quantization Hamiltonian Molecular gradients 7 Geometrical Properties Bond Distances Harmonic and Anharmonic Constants Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
5 Introduction Time-Independent Molecular Properties When a molecular system is perturbed, its total energy changes E(µ) = E (0) + E (1) µ E(2) µ 2 + The expansion coefficients are characteristic of the molecule and its quantum state we refer to these coefficients as molecular properties When the perturbation is static, the properties may be calculated by differentiation E (1) = de dµ µ=0 E (2) = d2 E dµ 2 µ=0 such properties are said to be time independent We do not here consider time-dependent molecular properties if periodic, these can be calculated in similar way using the quasi-energy Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
6 Introduction Examples of Derivatives Responses to geometrical perturbations forces and force constants spectroscopic constants Responses to external electromagnetic fields permanent and induced moments polarizabilities and magnetizabilities optical activity Responses to external magnetic fields and nuclear magnetic moments NMR shielding and indirect spin spin coupling constants EPR hyperfine coupling constants and g values Responses to nuclear quadrupole moments nuclear field gradients, quadrupole coupling constants Responses to molecular rotation spin rotation constants and molecular g values Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
7 Introduction Numerical vs. analytical differentiation Numerical differentiation (finite differences and polynomial fitting) often simple to implement (at least for real singlet perturbations) difficulties related to numerical accuracy and computational efficiency Analytical differentiation (derivatives calculated from analytical expressions) considerable programming effort required greater speed, precision, and convenience Implementations of analytical techniques first-order properties (dipole moments and gradients) second-order properties (polarizabilities and Hessians, NMR parameters) Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
8 Energy Functions Section 2 Energy Functions Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
9 Energy Functions Table of Contents 1 Introduction 2 Energy Functions Variational and nonvariational wave functions 3 Derivatives for Variational Wave Functions Molecular gradients Molecular Hessians 4 Derivatives for Nonvariational Wave Functions Nonvariational Wave Functions Lagrangian method 5 Examples of derivatives Hartree Fock Molecular Gradient FCI Lagrangian and Derivatives 6 Derivatives in second quantization Second-quantization Hamiltonian Molecular gradients 7 Geometrical Properties Bond Distances Harmonic and Anharmonic Constants Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
10 Energy Functions Electronic Energy Function The electronic energy function contains the Hamiltonian and the wave function: It depends on two distinct sets of parameters: E(x, λ) = λ H(x) λ x: external (perturbation) parameters (geometry, external field) λ: electronic (wave-function) parameters (MOs, cluster amplitudes) The Hamiltonian (here in second quantization) H(x) = h pq(x)e pq + 1 g 2 pqrs(x)e pqrs + h nuc(x) pq pqrs depends explicitly on the external parameters: h pq(x) = φ p(x) h(x) φ q(x) The wave function λ depends implicitly on the external parameters λ(x). Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
11 Energy Functions Electronic Energy and its Derivatives The electronic energy E(x) is obtained by optimizing the energy function E(x, λ) with respect to λ for each value of x: E(x) = E(x, λ ) note: the optimization is not necessarily a variational minimization Our task is to calculate derivatives of E(x) with respect to x: de(x) E(x, λ ) E(x, λ) = + λ dx x λ }{{} λ=λ x λ=λ }{{} explicit dependence implicit dependence the implicit as well as explicit dependence must be accounted for The quantity λ/ x is the wave-function response it tells us how the electronic structure changes when the system is perturbed To proceed, we need to make a distinction between variationally determined wave functions nonvariationally determined wave functions Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
12 Energy Functions Variational and nonvariational wave functions Variational and Nonvariational Wave Functions Variational wave functions the optimized energy fulfils the stationary (variational) condition: E(x, λ) = 0 (for all x) λ the Hartree Fock energy in an unconstrained exponential parameterization HF = exp( κ) 0, κ = κ the energy of the full CI (FCI) wave function FCI = i c i i as an expectation value: E FCI (x, c) = 0 c Nonvariational wave functions wave functions whose energy does not fulfil the stationary condition: E(x, λ) 0 λ the Hartree Fock and Kohn Sham energies in a constrained LCAO parameterization (orthonormality) HF = 1 N! det φ 1, φ 2,... φ N, φ p(r; x) = µ Cµpχµ(r; x), φp φ q = δpq the truncated CI energy with respect to orbital rotations: E CI (x, c, κ) = 0, c E CI (x, c, κ) 0 κ Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
13 Derivatives for Variational Wave Functions Section 3 Derivatives for Variational Wave Functions Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
14 Derivatives for Variational Wave Functions Table of Contents 1 Introduction 2 Energy Functions Variational and nonvariational wave functions 3 Derivatives for Variational Wave Functions Molecular gradients Molecular Hessians 4 Derivatives for Nonvariational Wave Functions Nonvariational Wave Functions Lagrangian method 5 Examples of derivatives Hartree Fock Molecular Gradient FCI Lagrangian and Derivatives 6 Derivatives in second quantization Second-quantization Hamiltonian Molecular gradients 7 Geometrical Properties Bond Distances Harmonic and Anharmonic Constants Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
15 Derivatives for Variational Wave Functions Molecular gradients Molecular Gradients for Variational Wave Functions Applying the chain rule, we obtain for the total derivative of the energy: de(x) dx = E(x, λ) x + E(x, λ) λ λ x the first term accounts for the explicit dependence on x the last term accounts for the implicit dependence on x We now invoke the stationary condition: E(x, λ) λ The molecular gradient then simplifies to = 0 (zero electronic gradient for all x) de(x) dx = E(x, λ) x examples: HF/KS and MCSCF molecular gradients (exponential parameterization) For variational wave functions, we do not need the response of the wave function λ/ x to calculate the molecular gradient de/dx. Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
16 Derivatives for Variational Wave Functions Molecular gradients Hellmann Feynman Theorem Assume that the (stationary) energy is an expectation value: E(x, λ) = λ H(x) λ The gradient is then given by the expectation-value expression: de(x) dx = E(x, λ) x = λ H x λ Relationship to first-order perturbation theory: E (1) = 0 H (1) 0 the Hellmann Feynman theorem The Hellmann Feynman theorem was originally stated for geometrical distortions: de Z K r ik = λ dr K r 3 λ + Z I Z K R IK i ik R 3 I K IK Classical interpretation: integration over the force operator Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
17 Derivatives for Variational Wave Functions Molecular Hessians Molecular Hessians for Variational Wave Functions Differentiating the molecular gradient, we obtain the molecular Hessian: d 2 E(x) dx 2 = d dx E(x, λ) x ( = x + λ x ) E(x, λ) λ x = 2 E(x, λ) x E(x, λ) λ x λ x we need the first-order response λ/ x to calculate the Hessian but we do not need the second-order response 2 λ/ x 2 for stationary energies To determine the response, we differentiate the stationary condition: E(x, λ) λ = 0 (all x) = d dx These are the first-order response equations: E(x, λ) λ = 2 E(x, λ) x λ 2 E λ λ }{{ 2 x = 2 E x λ }}{{} electronic right-hand Hessian side = E(x, λ) λ λ 2 x = 0 Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
18 Derivatives for Variational Wave Functions Molecular Hessians Response Equations The molecular Hessian for stationary energies: The response equations: d 2 E dx 2 = 2 E x E λ x λ x electronic Hessian 2 E λ 2 λ x = 2 E λ x perturbed electronic gradient the electronic Hessian is a Hermitian matrix, independent of the perturbation its dimensions are usually large and it cannot be constructed explicitly the response equations are typically solved by iterative techniques key step: multiplication of the Hessian with a trial vector Analogy with Hooke s law: force constant kx = F force the wave function relaxes by an amount proportional to the perturbation Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
19 Derivatives for Variational Wave Functions Molecular Hessians 2n + 1 Rule For molecular gradients and Hessians, we have the expressions de dx = E x d 2 E dx 2 = 2 E x E λ x λ x In general, we have the 2n + 1 rule: zero-order response needed first-order response needed For variational wave functions, the derivatives of the wave function to order n determine the derivatives of the energy to order 2n + 1. Examples: wave-function responses needed to fourth order: energy E (0) E (1) E (2) E (3) E (4) wave function λ (0) λ (0) λ (0), λ (1) λ (0), λ (1) λ (0), λ (1), λ (2) Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
20 Derivatives for Nonvariational Wave Functions Section 4 Derivatives for Nonvariational Wave Functions Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
21 Derivatives for Nonvariational Wave Functions Table of Contents 1 Introduction 2 Energy Functions Variational and nonvariational wave functions 3 Derivatives for Variational Wave Functions Molecular gradients Molecular Hessians 4 Derivatives for Nonvariational Wave Functions Nonvariational Wave Functions Lagrangian method 5 Examples of derivatives Hartree Fock Molecular Gradient FCI Lagrangian and Derivatives 6 Derivatives in second quantization Second-quantization Hamiltonian Molecular gradients 7 Geometrical Properties Bond Distances Harmonic and Anharmonic Constants Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
22 Derivatives for Nonvariational Wave Functions Nonvariational Wave Functions Nonvariational Wave Functions The 2n + 1 rule simplifies property evaluation for variational wave functions What about the nonvariational wave functions? any energy may be made stationary by Lagrange s method of undetermined multipliers the 2n + 1 rule is therefore of general interest Example: the CI energy the CI energy function is given by: E CI (x, c, κ) CI parameters c orbital-rotation parameters κ it is nonstationary with respect to the orbital-rotation parameters: We shall now consider its molecular gradient: E CI (x, c, κ) = 0 stationary c E CI (x, c, κ) 0 nonstationary κ 1 by straightforward differentiation of the CI energy 2 by differentiation of the CI Lagrangian In coupled-cluster theory, all parameters are nonvariationally determined Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
23 Derivatives for Nonvariational Wave Functions Nonvariational Wave Functions CI Molecular Gradients the Straightforward Way Straightforward differentiation of E CI (x, c, κ) gives the expression de CI dx = E CI x = E CI x + E CI c + E CI κ c x + E CI κ κ x κ x κ contribution does not vanish it appears that we need the first-order response of the orbitals The HF orbitals used in CI theory fulfil the following conditions at all geometries: E HF κ = 0 HF stationary conditions we obtain the orbital responses by differentiating this equation with respect to x: 2 E HF κ 2 κ x = 2 E HF x κ 1st-order response equations one such set of equations must be solved for each perturbation Calculated in this manner, the CI gradient becomes expensive Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
24 Derivatives for Nonvariational Wave Functions Lagrangian method Lagrange s Method of Undetermined Multipliers To calculate the CI energy, we minimize E CI with respect to c and κ: min E CI(x, c, κ) c,κ subject to the constraints Use Lagrange s method of undetermined multipliers: E HF (x, κ) κ construct the CI Lagrangian by adding these constraints with multipliers to the energy: ( ) EHF (x, κ) L CI (x, c, κ, κ) = E CI (x, c, κ) + κ 0 κ adjust the Lagrange multipliers κ such that the Lagrangian becomes stationary: = 0 L CI c L CI κ L CI κ = 0 = E CI c = 0 = E CI κ = 0 = E HF κ = 0 CI conditions + κ 2 E HF κ 2 = 0 linear set of equations for κ = 0 HF conditions note the duality between κ and κ Note that E CI = L CI when the Lagrangian is stationary we now have a stationary CI energy expression L CI Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
25 Derivatives for Nonvariational Wave Functions Lagrangian method CI Molecular Gradients the Easy Way The CI Lagrangian is given by L CI = E CI + κ E HF κ stationary with respect to all variables Since the Lagrangian is stationary, we may invoke the 2n + 1 rule: de CI dx = dl CI dx = L CI x zero-order response equations = E CI x κ 2 E HF κ 2 This result should be contrasted with the original expression first-order response equations de CI dx = E CI x 2 E HF κ 2 + κ 2 E HF κ x = E CI κ + E CI κ κ x κ x = 2 E HF κ x We have greatly reduced the number of response equations to be solved Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
26 Derivatives for Nonvariational Wave Functions Lagrangian method Lagrange s Method Summarized (1) Establish the energy function E(x, λ) and identify conditions on the variables (2) Set up the Lagrangian energy function: L(x, λ, λ) }{{} Lagrangian e(x, λ) = 0 = E(x, λ) }{{} + λ (e(x, λ) 0) }{{} energy function constraints (3) Satisfy the stationary conditions for the variables and their multipliers: L = e(x, λ) = 0 λ L λ = E λ + λ e λ = 0 condition for λ determines λ condition for λ determines λ note the duality between λ and λ! (4) Calculate derivatives from the stationary Lagrangian The Lagrangian approach is generally applicable: it gives the Hylleraas functional when applied to a perturbation expression it may be generalized to time-dependent properties Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
27 Derivatives for Nonvariational Wave Functions Lagrangian method 2n + 1 and 2n + 2 Rules For variational wave functions, we have the 2n + 1 rule: λ (n) determines the energy to order 2n + 1. The Lagrangian technique extends this rule to nonvariational wave functions For the new variables the multipliers the stronger 2n + 2 rule applies: λ (n) determines the energy to order 2n + 2. Responses required to order 10: E (n) λ (k) λ (k) Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
28 Examples of derivatives Section 5 Examples of derivatives Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
29 Examples of derivatives Table of Contents 1 Introduction 2 Energy Functions Variational and nonvariational wave functions 3 Derivatives for Variational Wave Functions Molecular gradients Molecular Hessians 4 Derivatives for Nonvariational Wave Functions Nonvariational Wave Functions Lagrangian method 5 Examples of derivatives Hartree Fock Molecular Gradient FCI Lagrangian and Derivatives 6 Derivatives in second quantization Second-quantization Hamiltonian Molecular gradients 7 Geometrical Properties Bond Distances Harmonic and Anharmonic Constants Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
30 Examples of derivatives Hartree Fock Molecular Gradient Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
31 Examples of derivatives Hartree Fock Molecular Gradient Hartree Fock Energy The MOs are expanded in atom-fixed AOs φ p(r; x) = Cµpχµ(r; x) µ The HF energy may be written in the general form E HF = pq Dpqhpq pqrs dpqrsgpqrs + K>L where h pq(x) = ( φ p(r, x) ) Z K K rk φ q(r, x) dr g pqrs(x) = φ p(r 1,x)φ q(r 1,x)φ r (r 2,x)φ s (r 2,x) r 12 dr 1 dr 2 the integrals depend explicitly on the geometry In closed-shell restricted HF (RHF) theory, the energy is given by E RHF = 2 h ii + ( ) 2giijj g ijji + i ij K>L summations over doubly occupied orbitals Z K Z L R KL Z K Z L R KL Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
32 Examples of derivatives Hartree Fock Molecular Gradient Hartree Fock Equations The HF energy is minimized subject to orthonormality constraints We therefore introduce the HF Lagrangian: L HF = E HF ε ( ) ij Sij δ ij ij S ij = φ i φ j = δij = D ij h ij + 1 ij 2 d ijkl g ijkl ε ( ) ij Sij δ ij + ijkl ij K>L The stationary conditions on the Lagrangian become: Z K Z L R KL L HF ε ij = S ij δ ij = 0 L HF C µi = E HF C µi ε S kl kl = 0 kl C µi the multiplier conditions are the orthonormality constraints the MO stationary conditions are the Roothaan Hall equations E HF C µi = ε S kl kl = F AO C = S AO Cε kl C µi Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
33 Examples of derivatives Hartree Fock Molecular Gradient Hartree Fock Molecular Gradient From the 2n + 1 rule, we obtain the RHF molecular gradient: de HF = dl HF = L HF = E HF dx dx x x S ij ε ij x ij In terms of MO integrals and density-matrix elements, we obtain the expression de HF dx = ij D ij h ij x We then transform to the AO basis: de HF dx = µν Dµν AO hµν AO x µνρσ ijkl d ijkl g ijkl x ij dµνρσ AO gµνρσ AO x µν ε ij S ij x + Fnuc ε AO Sµν AO µν + F nuc x density matrices transformed to AO basis derivative integrals added directly to gradient elements Important points: the gradient does not involve MO differentiation because of the 2n + 1 rule the time-consuming step is integral differentiation Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
34 Examples of derivatives FCI Lagrangian and Derivatives FCI Energy Consider a normalized CI wave function: c = n=0 cn n, ct c = 1, m n = δ mn The basis functions n are the normalized CI eigenstates of the unperturbed problem: m H n = δ mne n, E 0 E 1 E 2 We assume that the ground-state energy function depends on two external parameters: c H(x, y) c = cm m H(x, y) n cn, mn n We construct a variational CI Lagrangian: L(x, y, c, µ) = ( c m m H(x, y) n c n µ mn n The stationary conditions are given by c 2 n = 1 ) cn 2 1 L(x, y, c, µ) = 0 c n = 2 n H(x, y) c 2µc n = 0 = H(x, y)c = E 0 (x, y)c L(x, y, c, µ) = 0 µ = cn 2 1 = 0 n = ct c = 1 the first condition is the CI eigenvalue problem with ground-state energy E 0(x, y) = µ the second condition is the CI normalization condition Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
35 Examples of derivatives FCI Lagrangian and Derivatives FCI Molecular Gradient and Hessian Using the CI Lagrangian, we calculate CI energy derivative in the usual way: de dx = L x, d 2 E dxdy = 2 L x y + n 2 L x c n c n x, 2 L c n c n m c n x = 2 L x c m By inverting the electronic Hessian, we obtain the more compact expression: d 2 E dxdy = 2 L x y 2 [ L 2 ] 1 L 2 L x c mn m c m c n c n y We next evaluate the various partial derivatives at x = y = 0 where c = 0 : L x = 0 H x 0, 2 L x y = 0 2 H x y 0, 2 L c m c n = 2 m H E 0 n = 2(E n E 0 )δ mn 2 L = 2 n H x c x 0 n Inserted above, we recover Rayleigh Schrödinger perturbation theory to second order: de dx = 0 H d 2 E x 0, dxdy = 0 2 H x y H x n n H y 0 E n n E 0 Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
36 Derivatives in second quantization Section 6 Derivatives in second quantization Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
37 Derivatives in second quantization Table of Contents 1 Introduction 2 Energy Functions Variational and nonvariational wave functions 3 Derivatives for Variational Wave Functions Molecular gradients Molecular Hessians 4 Derivatives for Nonvariational Wave Functions Nonvariational Wave Functions Lagrangian method 5 Examples of derivatives Hartree Fock Molecular Gradient FCI Lagrangian and Derivatives 6 Derivatives in second quantization Second-quantization Hamiltonian Molecular gradients 7 Geometrical Properties Bond Distances Harmonic and Anharmonic Constants Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
38 Derivatives in second quantization Second-quantization Hamiltonian Second-Quantization Hamiltonian In second quantization, the Hamiltonian operator is given by: H = pq h pqa pa q + pqrs g pqrsa pa r a sa q + h nuc h pq = φ p (r) h(r) φq(r) g pqrs = φ p (r 1 )φ r (r 2 ) r 1 12 φq(r 1)φ s(r 2 ) Its construction assumes an orthonormal basis of MOs φ p: [a p, a q] + = 0, [a p, a q] + = 0, [a p, a q] + = δ pq The MOs are expanded in AOs, which often depend explicitly on the perturbation such basis sets are said to be perturbation-dependent: φ p(r) = µ C pµ χ µ(r, x) we must make sure that the MOs remain orthonormal for all x this introduces complications as we take derivatives with respect to x Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
39 Derivatives in second quantization Second-quantization Hamiltonian MOs and Hamiltonian at Distorted geometries 1. Orthonormal MOs at the reference geometry: 2. Geometrical distortion x = x 0 + x: φ(x 0 ) = C (0) χ(x 0 ) S(x 0 ) = φ(x 0 ) φ (x 0 ) = I φ(x) = C (0) χ(x) S(x) = φ(x) φ (x) I note: this basis is nonorthogonal and not useful for setting up the Hamiltonian. 3. Orthonormalize the basis set (e.g., by Löwdin orthonormalization): ψ(x) = S 1/2 (x)φ(x) S(x) = S 1/2 (x)s(x)s 1/2 (x) = I 4. From these orthonormalized MOs (OMOs) ψ p, construct Hamiltonian in the usual manner H = h pq pqa pa q + pqrs gpqrsa pa r a sa q + h nuc h pq = ψ p (r) h(r) ψ q(r) g pqrs = ψ p (r 1 )ψ r (r 2 ) r 1 12 ψq(r 1)ψ s(r 2 ) Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
40 Derivatives in second quantization Second-quantization Hamiltonian Hamiltonian at all Geometries The Hamiltonian is now well defined at all geometries: H(x) = h pq(x)e pq(x) + 1 g 2 pqrs(x)e pqrs(x) + h nuc(x) pq pqrs The OMO integrals are given by h pq(x) = h mn(x)[s 1/2 ] mp(x)[s 1/2 ] nq(x) mn in terms of the usual MO integrals h mn(x) = µν C mµc (0) nν (0) hµν AO (x), Smn(x) = C mµc (0) nν (0) Sµν AO (x) µν and similarly for the two-electron integrals. What about the geometry dependence of the excitation operators? this may be neglected when calculating derivatives since, for all geometries, ] [a p(x), a q (x) = S pq(x) = δ pq + Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
41 Derivatives in second quantization Molecular gradients HF Molecular Gradients in Second Quantization The molecular gradient now follows from the Hellmann Feynman theorem: E (1) = 0 H (1) 0 = D h(1) pq pq + 1 d 2 pqrs g pqrs (1) + h nuc (1) pq pqrs We need the derivatives of the OMO integrals: h pq (1) = [ hmn(s 1/2 ) mp(s 1/2 ] (1) (1) ) nq = h pq 1 2 mn m S (1) pmh (0) mq 1 2 The gradient may therefore be written in the form E (1) = D pqh (1) pq + 1 d 2 pqrsg pqrs (1) F pqs pq (1) + h nuc, (1) pq pqrs pq m h pms (0) mq (1) where the generalized Fock matrix is given by: F pq = n D pnh qn + d pnrsg qnrs nrs For RHF theory, this result is equivalent to that derived in first quantization Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
42 Geometrical Properties Section 7 Geometrical Properties Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
43 Geometrical Properties Table of Contents 1 Introduction 2 Energy Functions Variational and nonvariational wave functions 3 Derivatives for Variational Wave Functions Molecular gradients Molecular Hessians 4 Derivatives for Nonvariational Wave Functions Nonvariational Wave Functions Lagrangian method 5 Examples of derivatives Hartree Fock Molecular Gradient FCI Lagrangian and Derivatives 6 Derivatives in second quantization Second-quantization Hamiltonian Molecular gradients 7 Geometrical Properties Bond Distances Harmonic and Anharmonic Constants Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
44 Geometrical Properties Geometrical Derivatives In the Born Oppenheimer approximation, the nuclei move on the electronic potential-energy surface E(x), which is a function of the nuclear geometry: E(x) = E 0 + E (1) x E(2) x 2 + expansion around the reference geometry The derivatives of this surface are therefore important: The geometrical derivatives are E (1) = de dx E (2) = d2 E dx 2 molecular gradient molecular Hessian used for locating and characterizing critical points related to spectroscopic constants, vibrational frequencies, and intensities Usually, only a few terms are needed in the expansions in some cases low-order expansions are inadequate or useless Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
45 Geometrical Properties Uses of Geometrical Derivatives To explore molecular potential-energy surfaces (3N 6 dimensions) localization and characterization of stationary points localization of avoided crossings and conical intersections calculation of reaction paths and reaction-path Hamiltonians application to direct dynamics To calculate spectroscopic constants molecular structure quadratic force constants and harmonic frequencies cubic and quartic force constants; fundamental frequencies partition functions dipole gradients and vibrational infrared intensities polarizability gradients and Raman intensities Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
46 Geometrical Properties Bond Distances Bond Distances Mean and mean abs. errors for 28 distances at the all-el. cc-pvx Z level (pm): 1-1 pvdz MP4 CCSD(T) MP2 CCSD MP3 DZ TZ QZ CCSD CCSD(T) pvtz pvqz CISD HF Bonds shorten with increasing basis: HF: DZ TZ 0.8 pm; TZ QZ 0.1 pm corr.: DZ TZ 1.6 pm; TZ QZ pm Bonds lengthen with improvements in the N-electron model: singles < doubles < triples < There is considerable scope for error cancellation: CISD/DZ, MP3/DZ Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
47 Geometrical Properties Bond Distances Bond Distances R e of BH, CO, N 2, HF, and F 2 (pm) CCSD(T) cc-pcvdz 2.1 CCSD(T) cc-pcvtz 0.2 CCSD(T) cc-pcvqz 0.1 CCSD(T) cc-pcv5z CCSD cc-pcvdz 1.2 CCSD cc-pcvtz 0.6 CCSD cc-pcvqz 0.9 CCSD cc-pcv5z MP2 cc-pcvdz 1.5 MP2 cc-pcvtz 0.9 MP2 cc-pcvqz 0.8 MP2 cc-pcv5z SCF cc-pcvdz 2.5 SCF cc-pcvtz 3.4 SCF cc-pcvqz 3.5 SCF cc-pcv5z Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
48 Geometrical Properties Bond Distances Contributions to Equilibrium Bond Distances (pm) RHF SD T Q 5 rel. adia. theory exp. err. HF N F CO We have agreement with experiment to within 0.01 pm except for F 2 Hartree Fock theory underestimates bond distances by up to 8.6 pm (for F 2 ) All correlation contributions are positive approximately linear convergence, slowest for F 2 triples contribute up to 2.0, quadruples up to 0.4, and quintuples 0.03 pm sextuples are needed for convergence to within 0.01 pm Relativistic corrections are small except for F 2 (0.05 pm) of the same magnitude and direction as the quintuples Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
49 Geometrical Properties Harmonic and Anharmonic Constants Harmonic Constants ω e of BH, CO, N 2, HF, and F 2 (cm 1 ) CCSD(T) cc-pcvdz 42 CCSD(T) cc-pcvtz 14 CCSD(T) cc-pcvqz 9 CCSD(T) cc-pcv5z CCSD cc-pcvdz 34 CCSD cc-pcvtz 64 CCSD cc-pcvqz 71 CCSD cc-pcv5z MP2 cc-pcvdz 68 MP2 cc-pcvtz 81 MP2 cc-pcvqz 73 MP2 cc-pcv5z SCF cc-pcvdz 269 SCF cc-pcvtz 288 SCF cc-pcvqz 287 SCF cc-pcv5z Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
50 Geometrical Properties Harmonic and Anharmonic Constants Contributions to Harmonic Frequencies ω e (cm 1 ) RHF SD T Q 5 rel. adia. theory exp. err. HF N F CO We have agreement with experiment to within 1 cm 1 except for F 2 Hartree Fock theory overestimates harmonic frequencies by up to 38% (in F 2 ). All correlation contributions are large and negative triples contribute up to 95 cm 1, quadruples 20 cm 1, and quintuples 4 cm 1 sextuples are sometimes needed for convergence to within 1 cm 1 The relativistic corrections are of the order of 1 cm 1 of the same magnitude and direction as the quadruples or quintuples Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
51 Geometrical Properties Harmonic and Anharmonic Constants Higher-Order Connected Contributions to ω e in N 2 (cm 1 ) There are substantial higher-order corrections: HF CCSD FC CCSD T FC CCSD T CCSDT CCSDTQ CCSDTQ5 connected triples relaxation contributes 9.7 cm 1 (total triples 70.5 cm 1 ) connected quadruples contribute 18.8 cm 1 connected quintuples contribute 3.9 cm 1 Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
52 Geometrical Properties Harmonic and Anharmonic Constants Anharmonic Constants ω e x e of BH, CO, N 2, HF, and F 2 (cm 1 ) CCSD(T) cc-pcvdz 1 CCSD(T) cc-pcvtz 1 CCSD(T) cc-pcvqz 0 CCSD(T) cc-pcv5z CCSD cc-pcvdz 2 CCSD cc-pcvtz 2 CCSD cc-pcvqz 2 CCSD cc-pcv5z MP2 cc-pcvdz 3 MP2 cc-pcvtz 3 MP2 cc-pcvqz 3 MP2 cc-pcv5z SCF cc-pcvdz 4 SCF cc-pcvtz 5 SCF cc-pcvqz 4 SCF cc-pcv5z Trygve Helgaker (University of Oslo) Time-independent molecular properties ESQC / 51
Derivatives and Properties
Derivatives and Properties Trygve Helgaker Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Norway Summer School: Modern Wavefunction Methods in Electronic
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 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 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 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 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 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 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 informationGradient Theory. unfamiliar with this particular formalism, the essential elements of second quantization are described in the Appendix (Section 9).
Gradient Theory Trygve Helgaker University of Oslo, Norway 1 Introduction 1157 2 Molecular Properties Calculated as erivatives 1157 3 Numerical Versus Analytical Evaluation of Molecular Properties 1159
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 informationMolecules in Magnetic Fields
Molecules in Magnetic Fields 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, Norway
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
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 informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Raman Centre for Atomic, Molecular and Optical
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 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 informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry and
More informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
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 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 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 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 informationMolecules in strong magnetic fields
Molecules in strong magnetic fields Trygve Helgaker, Kai Lange, Alessandro Soncini, and Erik Tellgren Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway
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 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 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 informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
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 informationMolecular Magnetic Properties ESQC 07. Overview
1 Molecular Magnetic Properties ESQC 07 Trygve Helgaker Department of Chemistry, University of Oslo, Norway Overview the electronic Hamiltonian in an electromagnetic field external and nuclear magnetic
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 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 informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway European Summer School in Quantum Chemistry
More informationPure and zero-point vibrational corrections to molecular properties
Pure and zero-point vibrational corrections to molecular properties Kenneth Ruud UiT The Arctic University of Norway I UNIVERSITETET TROMSØ June 30 2015 Outline Why consider vibrational effects? General
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 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 informationIntroduction to the calculation of molecular properties by response theory
Introduction to the calculation of molecular properties by response theory Julien Toulouse Laboratoire de Chimie Théorique Université Pierre et Marie Curie et CNRS, 75005 Paris, France julien.toulouse@upmc.fr
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 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 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 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 informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway European Summer School in Quantum Chemistry
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 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 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 informationDalton Quantum Chemistry Program
1 Quotation from home page: Dalton Quantum Chemistry Program Dalton QCP represents a powerful quantum chemistry program for the calculation of molecular properties with SCF, MP2, MCSCF or CC wave functions.
More informationNMR and IR spectra & vibrational analysis
Lab 5: NMR and IR spectra & vibrational analysis A brief theoretical background 1 Some of the available chemical quantum methods for calculating NMR chemical shifts are based on the Hartree-Fock self-consistent
More informationMolecular electronic structure in strong magnetic fields
Trygve Helgaker 1, Mark Hoffmann 1,2 Kai Molecules Lange in strong 1, Alessandro magnetic fields Soncini 1,3, andcctcc Erik Jackson Tellgren 211 1 (CTCC, 1 / 23 Uni Molecular electronic structure in strong
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 adiabatic connection
The adiabatic connection Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Dipartimento di Scienze
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 informationDensity-functional theory in quantum chemistry. Trygve Helgaker. From Quarks to the Nuclear Many-Body Problem
1 Density-functional theory in quantum chemistry Trygve Helgaker Centre for Theoretical and Computational Chemistry, University of Oslo, Norway From Quarks to the Nuclear Many-Body Problem A conference
More informationMolecules in strong magnetic fields
Molecules in strong magnetic fields Trygve Helgaker, Kai Lange, Alessandro Soncini, and Erik Tellgren Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway
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 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 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 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 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 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 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 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 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 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 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 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 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 informationP. W. Atkins and R. S. Friedman. Molecular Quantum Mechanics THIRD EDITION
P. W. Atkins and R. S. Friedman Molecular Quantum Mechanics THIRD EDITION Oxford New York Tokyo OXFORD UNIVERSITY PRESS 1997 Introduction and orientation 1 Black-body radiation 1 Heat capacities 2 The
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 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 informationGeometry optimization
Geometry optimization Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway European Summer School in Quantum Chemistry (ESQC) 211 Torre
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 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 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 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 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 informationAccurate multireference configuration interaction calculations on the lowest 1 and 3 electronic states of C 2,CN, BN, and BO
Accurate multireference configuration interaction calculations on the lowest 1 and 3 electronic states of C 2,CN, BN, and BO Kirk A. Peterson a) Department of Chemistry, Washington State University and
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 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 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 informationOn the Uniqueness of Molecular Orbitals and limitations of the MO-model.
On the Uniqueness of Molecular Orbitals and limitations of the MO-model. The purpose of these notes is to make clear that molecular orbitals are a particular way to represent many-electron wave functions.
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 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 informationQuantum chemical modelling of molecular properties - parameters of EPR spectra
Quantum chemical modelling of molecular properties - parameters of EPR spectra EPR ( electric paramagnetic resonance) spectra can be obtained only for open-shell systems, since they rely on transitions
More informationSection 9 Variational Method. Page 492
Section 9 Variational Method Page 492 Page 493 Lecture 27: The Variational Method Date Given: 2008/12/03 Date Revised: 2008/12/03 Derivation Section 9.1 Variational Method: Derivation Page 494 Motivation
More informationHighly Accurate Ab Initio Computation of Thermochemical Data
Chapter 1 Highly Accurate Ab Initio Computation of Thermochemical Data Trygve Helgaker Department of Chemistry, University of Oslo, P. O. Box 1033 Blindern, N-0315 Oslo, Norway Wim Klopper and Asger Halkier
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 informationMolecular Magnetism. Molecules in an External Magnetic Field. Trygve Helgaker
Molecular Magnetism Molecules in an External Magnetic Field Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Laboratoire de
More informationDiamagnetism and Paramagnetism in Atoms and Molecules
Diamagnetism and Paramagnetism in Atoms and Molecules Trygve Helgaker Alex Borgoo, Maria Dimitrova, Jürgen Gauss, Florian Hampe, Christof Holzer, Wim Klopper, Trond Saue, Peter Schwerdtfeger, Stella Stopkowicz,
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 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 informationDirect optimization of the atomic-orbital density matrix using the conjugate-gradient method with a multilevel preconditioner
JOURNAL OF CHEMICAL PHYSICS VOLUME 115, NUMBER 21 1 DECEMBER 2001 Direct optimization of the atomic-orbital density matrix using the conjugate-gradient method with a multilevel preconditioner Helena Larsen,
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 informationLecture 9. Hartree Fock Method and Koopman s Theorem
Lecture 9 Hartree Fock Method and Koopman s Theorem Ψ(N) is approximated as a single slater determinant Φ of N orthogonal One electron spin-orbitals. One electron orbital φ i = φ i (r) χ i (σ) χ i (σ)
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 informationChem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.
Chem 44 Review for Exam Hydrogenic atoms: The Coulomb energy between two point charges Ze and e: V r Ze r Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative
More informationQuantum mechanics can be used to calculate any property of a molecule. The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is,
Chapter : Molecules Quantum mechanics can be used to calculate any property of a molecule The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is, E = Ψ H Ψ Ψ Ψ 1) At first this seems like
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 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 information