Derivatives and Properties

Size: px
Start display at page:

Download "Derivatives and Properties"

Transcription

1 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 Structure Theory Wissenschaftspark Gelsenkirchen, Gelsenkirchen, Germany September 30 October 5, 2018 Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

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 8 Electronic Hamiltonian 9 London Orbitals and Gauge-Origin Transformations 10 Zeeman and Hyperfine Interactions 11 First-Order Magnetic Properties 12 Molecular Magnetizabilities 13 High-Resolution NMR Spectra Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

3 Introduction Section 1 Introduction Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

4 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

5 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

6 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

7 Energy Functions Section 2 Energy Functions Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

8 Energy Functions Energy and 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

9 Energy Functions Energy and 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

10 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 the optimized energy does not fulfil the stationary (variational) 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

11 Derivatives for Variational Wave Functions Section 3 Derivatives for Variational Wave Functions Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

12 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

13 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

14 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

15 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

16 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

17 Derivatives for Nonvariational Wave Functions Section 4 Derivatives for Nonvariational Wave Functions Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

18 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: E CI (x, c, κ) = 0 stationary c E CI (x, c, κ) 0 nonstationary κ We shall now consider its molecular gradient: 1 by differentiation of the CI energy 2 by differentiation of the CI Lagrangian In coupled-cluster theory, all parameters are nonvariationally determined Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

19 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

20 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

21 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

22 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

23 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

24 Examples of derivatives Section 5 Examples of derivatives Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

25 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 Z K Z L R KL where the one- and two-electron integrals are given by ( h pq(x) = φ p(r, x) ) Z K K rk φ q(r, x) dr φp(r1, x)φ q(r 1, x)φ r (r 2, x)φ s(r 2, x) g pqrs(x) = dr 1 dr 2 r 12 note: all 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 Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

26 Examples of derivatives Hartree Fock Molecular Gradient Hartree Fock Equations The HF energy is minimized subject to orthonormality constraints S ij = φ i φ j = δij We therefore introduce the HF Lagrangian: L HF = E HF ε ( ) ij Sij δ ij 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 Note: 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

27 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 ijkl d ijkl g ijkl x We then express the gradient in terms of AO integrals: ij ε ij S ij x + Fnuc de HF dx = µν Dµν AO hµν AO x µνρσ dµνρσ AO gµνρσ AO x µν ε 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

28 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

29 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

30 Derivatives in second quantization Section 6 Derivatives in second quantization Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

31 Derivatives in second quantization Second-quantization Hamiltonian Second-Quantization Hamiltonian In second quantization, the Hamiltonian operator is given by: H = pq h pqa p aq + g pqrsa p a r asaq + hnuc pqrs 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, [ap, 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

32 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 p aq + pqrs gpqrsa p a r asaq + hnuc 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

33 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

34 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 pq h (1) 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

35 Geometrical Properties Section 7 Geometrical Properties Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

36 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

37 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

38 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

39 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

40 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

41 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

42 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

43 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

44 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 (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

45 Geometrical Properties Harmonic and anharmonic constants 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 8 Electronic Hamiltonian 9 London Orbitals and Gauge-Origin Transformations 10 Zeeman and Hyperfine Interactions 11 First-Order Magnetic Properties 12 Molecular Magnetizabilities 13 High-Resolution NMR Spectra Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

46 Electronic Hamiltonian Section 8 Electronic Hamiltonian Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

47 Electronic Hamiltonian Hamiltonian Mechanics In classical Hamiltonian mechanics, a system of particles is described in terms their positions q i and conjugate momenta p i For each system, there exists a scalar Hamiltonian function H(q i, p i ) such that the classical equations of motion are given by: q i = H p i, ṗ i = H q i note: the Hamiltonian H is not unique! (Hamilton s equations) Example: a single particle of mass m in a conservative force field F (q) the Hamiltonian is constructed from the corresponding scalar potential: H(q, p) = p2 2m + V (q), V (q) F (q) = q Hamilton s equations of motion are equivalent to Newton s equations: q = H(q,p) = p } p m ṗ = H(q,p) V (q) = m q = F (q) (Newton s equations) = q q Hamilton s equations are first-order differential equations Newton s are second-order Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

48 Electronic Hamiltonian Quantization of a Particle in a Conservative Force Field The Hamiltonian formulation is more general than the Newtonian formulation: it is invariant to coordinate transformations it provides a uniform description of matter and field it constitutes the springboard to quantum mechanics The Hamiltonian function (total energy) of a particle in a conservative force field: H(q, p) = p2 2m + V (q) Standard rule for quantization (in Cartesian coordinates): carry out the operator substitutions p i, H i t multiply the resulting expression by the wave function Ψ(q) from the right: i Ψ(q) ] = [ 2 t 2m 2 + V (q) Ψ(q) This approach is sufficient for a treatment of electrons in an electrostatic field it is insufficient for nonconservative systems it is therefore inappropriate for systems in a general electromagnetic field Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

49 Electronic Hamiltonian Particle in a Lorentz Force Field Lorentz Force and Maxwell s Equations In the presence of an electric field E and a magnetic field (magnetic induction) B, a classical particle of charge z experiences the Lorentz force: F = z (E + v B) since this force depends on the velocity v of the particle, it is not conservative The electric and magnetic fields E(r, t) and B(r, t) satisfy Maxwell s equations ( ): E = ρ/ε 0 Coulomb s law B ε 0 µ 0 E/ t = µ 0 J Ampère s law with Maxwell s correction B = 0 E + B/ t = 0 Faraday s law of induction where ρ(r, t) and J(r, t) are the charge and current densities, respectively Note: when ρ and J are known, Maxwell s equations can be solved for E and B but the particles are driven by the Lorentz force, so ρ and J are functions of E and B We here consider the motion of particles in a given (fixed) electromagnetic field Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

50 Electronic Hamiltonian Particle in a Lorentz Force Field Scalar and Vector Potentials The second, homogeneous pair of Maxwell s equations involves only E and B: 1 Eq. (1) is satisfied by introducing the vector potential A: B = 0 (1) E + B t = 0 (2) B = 0 = B = A vector potential (3) 2 inserting Eq. (3) in Eq. (2) and introducing a scalar potential φ, we obtain ( E + A ) = 0 = E + A t t = φ scalar potential The second pair of Maxwell s equations is thus automatically satisfied by writing E = φ A t B = A The potentials (φ, A) contain four rather than six components as in (E, B). φ and A are obtained by solving the inhomogeneous pair of Maxwell s equations, containing ρ and J Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

51 Electronic Hamiltonian Particle in a Lorentz Force Field Gauge Transformations Consider the following gauge transformation of the potentials: φ = φ f } t A with f = f (q, t) gauge function of position and time = A + f Such a transformation of the potentials does not affect the physical fields: E = φ A t = φ + f t A t f t B = A = (A + f ) = B + f = B Conclusion: the scalar and vector potentials φ and A are not unique we are free to choose f (q, t) to make the potentials satisfy additional conditions typically, we require the vector potential to be divergenceless: A = 0 = (A + f ) = 0 = 2 f = A = E Coulomb gauge We shall always assume that the vector potential satisfies the Coulomb gauge: A = B, A = 0 Coulomb gauge note: A is still not uniquely determined, the following transformation being allowed: A = A + f, 2 f = 0 Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

52 Electronic Hamiltonian Particle in a Lorentz Force Field Hamiltonian in an Electromagnetic Field We must construct a Hamiltonian function such that Hamilton s equations are equivalent to Newton s equation with the Lorentz force: q i = H p i & ṗ i = H q i ma = z (E + v B) To this end, we introduce scalar and vector potentials φ and A such that E = φ A t, B = A In terms of these potentials, the classical Hamiltonian function becomes H = π2 + zφ, π = p za kinetic momentum 2m Quantization is then accomplished in the usual manner, by the substitutions p i, H i t The time-dependent Schrödinger equation for a particle in an electromagnetic field: i Ψ t = 1 ( i za) ( i za) Ψ + zφ Ψ 2m Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

53 Electronic Hamiltonian Electron Spin Electron Spin The nonrelativistic Hamiltonian for an electron in an electromagnetic field is then given by: H = π2 2m eφ, π = i + ea However, this description ignores a fundamental property of the electron: spin. Spin was introduced by Pauli in 1927, to fit experimental observations: (σ π)2 π2 H = eφ = 2m 2m + e 2m B σ eφ where σ contains three operators, represented by the two-by-two Pauli spin matrices ( 0 1 σ x = 1 0 ) ( 0 i, σ y = i 0 ) ( 1 0, σ z = 0 1 ) The Schrödinger equation now becomes a two-component equation: ( ) π 2 e eφ + 2m 2m Bz e (Ψα ) ( ) (Bx iby ) 2m Ψα e 2m (Bx + iby ) π 2 = E e eφ 2m 2m Bz Ψ β Ψ β Note: the two components are coupled only in the presence of an external magnetic field Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

54 Electronic Hamiltonian Electron Spin Spin and Relativity The introduction of spin by Pauli in 1927 may appear somewhat ad hoc By contrast, spin arises naturally from Dirac s relativistic treatment in 1928 is spin a relativistic effect? However, reduction of Dirac s equation to nonrelativistic form yields the Hamiltonian H = (σ π)2 2m π2 eφ = 2m + e 2m B σ eφ π2 2m eφ in this sense, spin is not a relativistic property of the electron on the other hand, in the nonrelativistic limit, all magnetic fields disappear... We interpret σ by associating an intrinsic angular momentum (spin) with the electron: s = σ/2 Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

55 Electronic Hamiltonian Molecular Electronic Hamiltonian Molecular Electronic Hamiltonian The nonrelativistic Hamiltonian for an electron in an electromagnetic field is therefore H = π2 2m + e B s eφ, π = p + ea, p = i m expanding π 2 and assuming the Coulomb gauge A = 0, we obtain π 2 Ψ = (p + ea) (p + ea) Ψ = p 2 Ψ + ep AΨ + ea pψ + e 2 A 2 Ψ = p 2 Ψ + e(p A)Ψ + 2eA pψ + e 2 A 2 Ψ = ( p 2 + 2eA p + e 2 A 2) Ψ in molecules, the dominant electromagnetic contribution is from the nuclear charges: φ = 1 Z K e 4πɛ 0 K + φ r ext K Summing over all electrons and adding pairwise Coulomb interactions, we obtain H = 1 2m p2 i e2 Z K + e2 r 1 ij 4πɛ i 0 r Ki ik 4πɛ 0 i>j + e A i p i + e B i s i e φ i m m i i i + e2 A 2 i 2m i zero-order Hamiltonian first-order Hamiltonian second-order Hamiltonian Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

56 Electronic Hamiltonian Molecular Electronic Hamiltonian Magnetic Perturbations In atomic units, the molecular Hamiltonian is given by H = H 0 + A(r i ) p i + B(r i ) s i i i i }{{}}{{} orbital paramagnetic spin paramagnetic There are two kinds of magnetic perturbation operators: φ(r i ) + 1 A 2 (r i ) 2 i }{{} diamagnetic the paramagnetic operator is linear and may lower or raise the energy the diamagnetic operator is quadratic and always raises the energy There are two kinds of paramagnetic operators: the orbital paramagnetic operator couples the field to the electron s orbital motion the spin paramagnetic operator couples the field to the electron s spin In the study of magnetic properties, we are interested in two types of perturbations: uniform external magnetic field B, with vector potential A ext(r) = 1 2 B r leads to Zeeman interactions nuclear magnetic moments M K, with vector potential A nuc(r) = α 2 K M K r K r 3 K leads to hyperfine interactions where α 1/137 is the fine-structure constant Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

57 London Orbitals and Gauge-Origin Transformations Section 9 London Orbitals and Gauge-Origin Transformations Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

58 London Orbitals and Gauge-Origin Transformations Hamiltonian in a Uniform Magnetic Field The nonrelativistic electronic Hamiltonian (implied summation over electrons): H = H 0 + A (r) p + B (r) s A (r)2 The vector potential of the uniform field B is given by: B = A = const = A O (r) = 1 2 B (r O) = 1 2 B r O note: the gauge origin O is arbitrary! The orbital paramagnetic interaction becomes: A O (r) p = 1 2 B (r O) p = 1 2 B (r O) p = 1 2 B L O where we have introduced the angular momentum relative to the gauge origin: The diamagnetic interaction becomes: L O = r O p 1 2 A2 O (r) = 1 8 (B r O) (B r O ) = 1 8 [ B 2 r 2 O (B r O) 2] The electronic Hamiltonian in a uniform magnetic field depends on the gauge origin: H = H B L O + B s [ B 2 r 2 O (B r O) 2] as we shall see, a change of the origin is a gauge transformation Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

59 London Orbitals and Gauge-Origin Transformations Gauge-Origin Transformations Gauge Transformation of Schrödinger Equation What is the effect of a gauge transformation on the wave function? Consider a general gauge transformation for the electron (atomic units): A = A + f, φ = φ f t It can be shown that this represents a unitary transformation of H i / t: ( H i ) ( = exp ( if ) H i ) exp (if ) t t In order that the Schrödinger equation is still satisfied ( H i ) Ψ t ( H i t ) Ψ, the wave function undergoes a compensating unitary transformation: Ψ = exp ( if ) Ψ All observable properties such as the electron density are then unaffected: ρ = (Ψ ) Ψ = [Ψ exp( if )] [exp( if )Ψ] = Ψ Ψ = ρ Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

60 London Orbitals and Gauge-Origin Transformations Gauge-Origin Transformations Gauge-Origin Transformations Different choices of gauge origin in the external vector potential are related by gauge transformations: A O (r) = 1 B (r O) 2 A G (r) = A O (r) A O (G) = A O (r) + f, f (r) = A O (G) r The exact wave function transforms accordingly and gives gauge-invariant results: Ψ exact G = exp [ if (r)] Ψ exact O = exp [ia O (G) r] Ψ exact O (rapid) oscillations Illustration: H 2 on the z axis in a magnetic field B = 0.2 a.u. in the y direction wave function with gauge origin at O = (0, 0, 0) (left) and G = (100, 0, 0) (right) Wave function, ψ Gauge transformed wave function, ψ" Re(ψ) Im(ψ) ψ Re(ψ") Im(ψ") ψ" Space coordinate, x (along the bond) Space coordinate, x (along the bond) Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

61 London Orbitals and Gauge-Origin Transformations London Orbitals London Orbitals The exact wave function transforms in the following manner: Ψ exact G (r) = exp [ i 1 2 B (G O) r] Ψ exact O (r) this behaviour cannot easily be modelled by standard atomic orbitals Let us build this behaviour directly into the atomic orbitals: ω lm (r K, B, G) = exp [ i 1 2 B (G K) r] χ lm (r K ) χ lm (r K ) is a normal atomic orbital centred at K and quantum numbers lm ω lm (r K, B, G) is a field-dependent orbital at K with field B and gauge origin G Each AO now responds in a physically sound manner to an applied magnetic field indeed, all AOs are now correct to first order in B, for any gauge origin G the calculations become rigorously gauge-origin independent uniform (good) quality follows, independent of molecule size These are the London orbitals after Fritz London (1937) Questions: also known as GIAOs (gauge-origin independent AOs or gauge-origin including AOs) are London orbitals needed in atoms? why not attach the phase factor to the total wave function instead? Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

62 London Orbitals and Gauge-Origin Transformations London Orbitals Dissociation With and Without London Orbitals Let us consider the FCI dissociation of H 2 in a magnetic field full lines: with London atomic orbitals dashed lines: without London atomic orbitals B B B 0.0 Without London orbitals, the FCI method is not size extensive in magnetic fields Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

63 Zeeman and Hyperfine Interactions Section 10 Zeeman and Hyperfine Interactions Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

64 Zeeman and Hyperfine Interactions Hamiltonian in Magnetic Field In atomic units, the molecular Hamiltonian is given by H = H 0 + A(r) p + B(r) s + 1 }{{}}{{} 2 A2 (r) }{{} orbital paramagnetic spin paramagnetic diamagnetic There are two kinds of magnetic perturbation operators: paramagnetic (may lower or raise energy) and diamagnetic (always raises energy) There are two kinds of paramagnetic operators: the orbital paramagnetic and spin paramagnetic First- and second-order Rayleigh Schrödinger perturbation theory gives: E (1) = 0 A p + B s 0 E (2) = 1 0 A n 0 A p + B s n n A p + B s 0 E n E 0 In the study of magnetic properties, we are interested in two types of perturbations: externally applied uniform magnetic fields B fields generated internally by nuclear magnetic moments M K Both fields are weak well described by perturbation (response) theory Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

65 Zeeman and Hyperfine Interactions Paramagnetic Operators Orbital Paramagnetic Interactions: A p Vector potentials corresponding to uniform fields and nuclear magnetic moments: A O = 1 2 B r O, A K = α 2 M K r K rk 3, α 1/137 the external field is typically about 10 4 a.u. (NMR experiments) the nuclear vector potential is exceedingly small (about 10 8 a.u.) since: α 2 = c a.u., M K = γ K I K 10 4 a.u. We obtain the following orbital paramagnetic operators: A O p = 1 2 B r O p = A K p = α 2 M K r K p r 3 K 1 2 B r O p = = α 2 M K r K p r 3 K 1 2 B L O orbital Zeeman = α 2 M K LK r 3 K orbital hyperfine interactions depend on angular momenta L O and L K relative to O and R K, respectively orbital hyperfine interaction expressed in terms of the paramagnetic spin orbit operator: A K p = M K h PSO K, hpso K = α 2 L K rk 3 These are imaginary singlet operators they have zero expectation values of closed-shell states they generate complex wave functions Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

66 Zeeman and Hyperfine Interactions Paramagnetic Operators Spin Paramagnetic Interactions: B s The spin interaction with the external uniform field B is trivial: B s spin Zeeman interaction should be compared with the orbital Zeeman interaction 1 2 B L O (different prefactor!) Taking the curl of A K, we obtain the nuclear magnetic field: B K = A K = 8πα2 3 δ(r K )M K + α 2 3r K (r K M K ) rk 2 M K rk 5 the first term is a contact interaction and contributes only at the nucleus the second term is a classical dipole field and contributes at a distance This magnetic field B K thus gives rise to two spin hyperfine interactions: B K s = M K (h FC K + hsd K ), h FC K h SD = 8πα2 3 δ(r K ) s Fermi contact (FC) K = α2 3r K rk T r K 2 I 3 r K 5 s spin dipole (SD) the FC operator contributes when the electron passes through the nucleus the SD operator is a classical dipolar interaction, decaying as r 3 K These are real triplet operators, which change the spin of the wave function they have zero expectation values of closed-shell states they couple closed-shell states to triplet states Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

67 Zeeman and Hyperfine Interactions Hamiltonian with Zeeman and Hyperfine Operators Perturbation Theory with Zeeman and Hyperfine Operators Hamiltonian with a uniform external field and with nuclear magnetic moments: H = H 0 + H (1) + H (2) = H 0 + H z (1) + H (1) hf A2 Zeeman interactions with the external magnetic field B: H (1) z = 1 2 B L O + B s 10 4 hyperfine interactions with the nuclear magnetic moments M K : H (1) hf = K M K h PSO K + K ( M K h FC K ) + hsd K 10 8 h PSO K = α 2 L K rk 3, h FC K = 8πα2 3 δ(r K ) s, h SD K = 3r α2 K rk T r K 2 I 3 rk 5 s Second-order Rayleigh Schrödinger perturbation theory: E (1) = 0 H z (1) + H (1) hf 0 E (2) = 1 0 A 2 0 H (1) z + H (1) 0 hf n n H (1) z + H (1) hf 0 2 E n n E 0 Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

68 Zeeman and Hyperfine Interactions Hamiltonian with Zeeman and Hyperfine Operators Zeeman and Hyperfine Interactions Zeeman SO SS, SO, OO SO PSO FC+SD hyperfine FC+SD PSO Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

69 Zeeman and Hyperfine Interactions Diamagnetic Operators Diamagnetic Operators: 1 2 A2 From A O and A K, we obtain three diamagnetic operators: A = A O + A K = A 2 = A O A O + 2A O A K + A K A K Their explicit forms and typical magnitudes (atomic units) are given by A O A O = 1 4 (B r O) (B r O ) 10 8 A O A K = α2 2 These are all real singlet operators (B r O ) (M K r K ) r 3 K A K A L = α 4 (M K r K ) (M L r L ) r 3 K r 3 L their expectation values contribute to second-order magnetic properties they are all exceedingly small but nonetheless all observable Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

70 First-Order Magnetic Properties Section 11 First-Order Magnetic Properties Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

71 First-Order Magnetic Properties Taylor Expansion of Energy Expand the energy in the presence of an external magnetic field B and nuclear magnetic moments M K around zero field and zero moments: E (B, M) = E 0 + perm. magnetic moments {}}{ B T E (10) BT E (20) B }{{} K magnetizability hyperfine coupling {}}{ K B T E (11) K M K } {{ } shieldings + 1 M T K E(01) K + 1 M T K 2 E(02) KL M L + KL }{{} spin spin couplings First-order terms vanish for closed-shell systems because of symmetry they shall be considered only briefly here Second-order terms are important for many molecular properties magnetizabilities nuclear shieldings constants of NMR nuclear spin spin coupling constants of NMR electronic g tensors of EPR (not dealt with here) Higher-order terms are negligible since the perturbations are tiny: 1) the magnetic induction B is weak ( 10 4 a.u.) 2) the nuclear magnetic moments M K couple weakly (µ 0 µ N 10 8 a.u.) Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

72 First-Order Magnetic Properties First-Order Magnetic Properties First-Order Molecular Properties The first-order properties are expectation values of H (1) Permanent magnetic moment M = 0 H z (1) 0 = L O + s 0 permanent magnetic moment dominates the magnetism of molecules the molecule reorients itself and enters the field such molecules are therefore paramagnetic Hyperfine coupling constants A K = 0 H (1) hf 0 = 8πα2 3 0 δ (rk ) s 0 MK + measure spin density at the nucleus important in electron paramagnetic resonance (EPR) recall: there are three hyperfine mechanisms: FC, SD and PSO Note: there are no first-order Zeeman or hyperfine couplings for closed-shell molecules all expectation values vanish for imaginary operators and triplet operators: c.c. ˆΩ imaginary c.c. c.c. ˆΩ triplet c.c. 0 Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

73 First-Order Magnetic Properties First-Order Magnetic Properties Atoms in a Magnetic Field Lowest states of the fluorine atom (left) and neon atom (right) in a magnetic field ( H = H 0 + Bs z BLz B2 x 2 + y 2) the ground-state neon atom is closed shell, with zero initial slope all other states are open shells, with a nonzero initial slope one atomic unit magnetic field strength is B 0 = T CCSD(T) calculations in uncontracted aug-cc-pcvqz basis (atomic units) Stopkowicz, Gauss, Lange, Tellgren, and Helgaker, JCP 143, (2015) Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

74 Molecular Magnetizabilities Section 12 Molecular Magnetizabilities Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

75 Molecular Magnetizabilities Molecular Magnetizabilities Expand the molecular electronic energy in the external magnetic field: E (B) = E 0 B T M 1 2 BT ξb + The magnetizability describes the second-order energy: 0.1 ξ = d2 E db 2 = 0 2 H B H B n n H 0.08 B E n n E 0 = 1 0 ro T 4 r OI 3 r O ro T L O n n L T O E }{{} n n E 0 } {{ } 0.1 diamagnetic term paramagnetic term The magnetizability describes the curvature at zero magnetic field: a) a) c) b) x b) d) x left: diamagnetic dependence on the field (ξ < 0); right: paramagnetic dependence on the field (ξ > 0) c) d) Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

76 Molecular Magnetizabilities Zeeman and Hyperfine Interactions Zeeman SO SS, SO, OO SO PSO FC+SD hyperfine FC+SD PSO Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

77 Molecular Magnetizabilities Basis-Set Convergence of Hartree Fock Magnetizabilities London orbitals are correct to first-order in the external magnetic field For this reason, basis-set convergence is usually improved RHF magnetizabilities of benzene: basis set χ xx χ yy χ zz London STO-3G G cc-pvdz aug-cc-pvdz origin CM STO-3G G cc-pvdz aug-cc-pvdz origin H STO-3G G cc-pvdz aug-cc-pvdz Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

78 Molecular Magnetizabilities Mean Absolute Errors for Magnetizabilities Mean relative errors (MREs, %) in magnetizabilities of 27 molecules relative to the CCSD(T)/aug-cc-pCV[TQ]Z values. The DFT results are grouped by functional type. The heights of the bars correspond to the largest MRE in each category. (Lutnæs et al., JCP 131, (2009)) Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

79 Molecular Magnetizabilities C 20 in a Perpendicular Magnetic Field All systems become diamagnetic in sufficiently strong fields: Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

80 High-Resolution NMR Spectra Section 13 High-Resolution NMR Spectra Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

81 High-Resolution NMR Spectra NMR Spin Hamiltonian High-Resolution NMR Spin Hamiltonian Consider a molecule in an external magnetic field B along the z axis and with nuclear spins I K related to the nuclear magnetic moments M K as: M K = γ K I K 10 4 a.u. where γ K is the magnetogyric ratio of the nucleus. Assuming free molecular rotation, the nuclear magnetic energy levels can be reproduced by the following high-resolution NMR spin Hamiltonian: where we have introduced H NMR = γ K (1 σ K )BI K z + γ K γ L 2 K KL I K I L K K>L }{{}}{{} nuclear Zeeman interaction nuclear spin spin interaction the nuclear shielding constants σ K the (reduced) indirect nuclear spin spin coupling constants K KL This is an effective nuclear spin Hamiltonian: it reproduces NMR spectra without considering the electrons explicitly the spin parameters σ K and K KL are adjusted to fit the observed spectra we shall consider their evaluation from molecular electronic-structure theory Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

82 High-Resolution NMR Spectra NMR Spin Hamiltonian Simulated 200 MHz NMR spectra of Vinyllithium 12 C 2 H 3 6 Li Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

83 NMR Shielding Constants Section 14 NMR Shielding Constants Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

84 NMR Shielding Constants Nuclear Shielding Constants Nuclear Shielding Constants Expansion of closed-shell energy in an external field B and nuclear magnetic moments M K : Here E (11) K E (B, M) = E BT E (20) B K B T E (11) K M K KL M T K E(02) KL M L + describes the coupling between the applied field and the nuclear moments: in the absence of electrons (i.e., in vacuum), this coupling is identical to I 3 : Hz nuc = B M K the purely nuclear Zeeman interaction K in the presence of electrons (i.e., in a molecule), the coupling is modified slightly: E (11) K = I 3 + σ K the nuclear shielding tensor Shielding constants arise from a hyperfine interaction between the electrons and the nuclei and Zeeman interaction between the electrons and the field they are of the order of α and are measured in ppm The nuclear Zeeman interaction does not enter the electronic problem compare with the nuclear nuclear Coulomb repulsion Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

85 NMR Shielding Constants Nuclear Shielding Constants Zeeman and Hyperfine Interactions Zeeman SO SS, SO, OO SO PSO FC+SD hyperfine FC+SD PSO Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

86 NMR Shielding Constants Nuclear Shielding Constants Ramsey s Expression for Nuclear Shielding Tensors Ramsey s expression for nuclear shielding tensors of a closed-shell system: σ K = d2 E el = 0 2 H dbdm K B M H B n n H M 0 K K E n n E 0 = α2 ro T 0 r K I 3 r O rk T 2 rk 3 0 α 0 LO n n r 3 2 K LT K 0 E n n E 0 }{{}}{{} diamagnetic term paramagnetic term The (usually) dominant diamagnetic term arises from differentiation of the operator: A O A K = 1 2 α2 r 3 K (B r O) (M K r K ) As for the magnetizability, there is no spin contribution for singlet states: S 0 0 singlet state For 1 S systems (closed-shell atoms), the paramagnetic term vanishes completely and the shielding is given by (assuming gauge origin at the nucleus): σ Lamb = α2 S r 1 K 1 S Lamb formula Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

87 NMR Shielding Constants Nuclear Shielding Constants Benchmark Calculations of BH Shieldings (ppm) σ( 11 B) σ( 11 B) σ( 1 H) σ( 1 H) HF MP CCSD CCSD(T) CCSDT CCSDTQ CISD CISDT CISDTQ FCI TZP+ basis, R BH = pm, all electrons correlated J. Gauss and K. Ruud, Int. J. Quantum Chem. S29 (1995) 437 M. Kállay and J. Gauss, J. Chem. Phys. 120 (2004) 6841 Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

88 NMR Shielding Constants Nuclear Shielding Constants Coupled-Cluster Convergence of Shielding Constants in CO (ppm) CCSD CCSD(T) CCSDT CCSDTQ CCSDTQ5 FCI σ( 13 C) σ( 13 C) σ( 17 O) σ( 17 O) All calculations in the cc-pvdz basis and with a frozen core. Kállay and Gauss, J. Chem. Phys. 120 (2004) Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

89 NMR Shielding Constants Nuclear Shielding Constants Calculated and Experimental Shielding Constants (ppm) HF CAS MP2 CCSD CCSD(T) exp. HF F ± 6 (300K) H ± 0.2 (300K) H 2O O ± 6 (300K) H ± 0.02 NH 3 N H ± 1.0 CH 4 C H F 2 F N 2 N ± 0.2 (300K) CO C ± 0.9 (eq) O ± 6 (eq) For references and details, see Chem. Rev. 99 (1999) 293. for exp. CO and H 2 O values, see Wasylishen and Bryce, JCP 117 (2002) Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

90 NMR Shielding Constants Nuclear Shielding Constants MAE (Exp.) MAE (Emp. Eq.) NMR: Mean Absolute Errors Relative to Experiment Mean absolute errors relative to experimental (blue) and empirical equilibrium values (red) MAE / ppm MAE (Exp.) MAE (Emp. Eq.) 0 Kohn Sham calculations give shielding constants of uneven quality errors increase when vibrational corrections are applied Teale et al. JCP 138, (2013) Trygve Helgaker (Hylleraas, University of Oslo) Derivatives and Properties MWM Summer School / 99

Time-independent molecular properties

Time-independent molecular properties Time-independent molecular properties Trygve Helgaker Hylleraas Centre, Department of Chemistry, University of Oslo, Norway and Centre for Advanced Study at the Norwegian Academy of Science and Letters,

More information

Molecules in Magnetic Fields

Molecules 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 information

Molecular Magnetic Properties

Molecular 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 information

Molecular Magnetic Properties

Molecular 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 information

Molecular Magnetic Properties ESQC 07. Overview

Molecular 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 information

Molecular Magnetism. Magnetic Resonance Parameters. Trygve Helgaker

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

More information

Molecular Magnetic Properties

Molecular 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 information

Molecular Magnetic Properties

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

More information

Molecular Magnetic Properties. The 11th Sostrup Summer School. Quantum Chemistry and Molecular Properties July 4 16, 2010

Molecular 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 information

Highly accurate quantum-chemical calculations

Highly 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 information

Molecular Magnetic Properties

Molecular 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 information

Basis sets for electron correlation

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

More information

Molecular Magnetic Properties. Workshop on Theoretical Chemistry. Mariapfarr, Salzburg, Austria February 14 17, Overview

Molecular Magnetic Properties. Workshop on Theoretical Chemistry. Mariapfarr, Salzburg, Austria February 14 17, Overview 1 Molecular Magnetic Properties Workshop on Theoretical Chemistry Mariapfarr, Salzburg, Austria February 14 17, 2006 Trygve Helgaker Department of Chemistry, University of Oslo, Norway Overview The electronic

More information

The Accurate Calculation of Molecular Energies and Properties: A Tour of High-Accuracy Quantum-Chemical Methods

The 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 information

DFT calculations of NMR indirect spin spin coupling constants

DFT 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 information

Convergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems

Convergence 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 information

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

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

More information

Oslo node. Highly accurate calculations benchmarking and extrapolations

Oslo 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 information

The Rigorous Calculation of Molecular Properties to Chemical Accuracy. T. Helgaker, Department of Chemistry, University of Oslo, Norway

The 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 information

Gradient Theory. unfamiliar with this particular formalism, the essential elements of second quantization are described in the Appendix (Section 9).

Gradient 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 information

Molecules in strong magnetic fields

Molecules 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 information

Introduction to Electronic Structure Theory

Introduction 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 information

Quantum Chemistry Methods

Quantum 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 information

P. W. Atkins and R. S. Friedman. Molecular Quantum Mechanics THIRD EDITION

P. 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 information

Molecular Magnetism. Molecules in an External Magnetic Field. Trygve Helgaker

Molecular 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 information

Lecture 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- 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 information

Molecular electronic structure in strong magnetic fields

Molecular 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 information

Introduction to Computational Chemistry

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

More information

Computational Methods. Chem 561

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

More information

Electric properties of molecules

Electric 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 information

Molecules in strong magnetic fields

Molecules 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 information

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

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

More information

Dalton Quantum Chemistry Program

Dalton 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 information

OVERVIEW OF QUANTUM CHEMISTRY METHODS

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

More information

Chemistry 3502/4502. Final Exam Part I. May 14, 2005

Chemistry 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 information

2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements

2. 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 information

The calculation of the universal density functional by Lieb maximization

The 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 information

NMR and IR spectra & vibrational analysis

NMR 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 information

Chemistry 334 Part 2: Computational Quantum Chemistry

Chemistry 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 information

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

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

More information

Coupled cluster theory: Analytic derivatives, molecular properties, and response theory

Coupled 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 information

Quantum chemical modelling of molecular properties - parameters of EPR spectra

Quantum 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 information

Wave function methods for the electronic Schrödinger equation

Wave 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 information

charges q r p = q 2mc 2mc L (1.4) ptles µ e = g e

charges q r p = q 2mc 2mc L (1.4) ptles µ e = g e APAS 5110. Atomic and Molecular Processes. Fall 2013. 1. Magnetic Moment Classically, the magnetic moment µ of a system of charges q at positions r moving with velocities v is µ = 1 qr v. (1.1) 2c charges

More information

Importing 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 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 information

2 Electronic structure theory

2 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 information

MO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )

MO 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 information

I. CSFs Are Used to Express the Full N-Electron Wavefunction

I. CSFs Are Used to Express the Full N-Electron Wavefunction Chapter 11 One Must be Able to Evaluate the Matrix Elements Among Properly Symmetry Adapted N- Electron Configuration Functions for Any Operator, the Electronic Hamiltonian in Particular. The Slater-Condon

More information

v(r i r j ) = h(r i )+ 1 N

v(r i r j ) = h(r i )+ 1 N Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j

More information

Attempts at relativistic QM

Attempts 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

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

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

More information

The adiabatic connection

The 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 information

Introduction to multiconfigurational quantum chemistry. Emmanuel Fromager

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

More information

The Particle-Field Hamiltonian

The Particle-Field Hamiltonian The Particle-Field Hamiltonian For a fundamental understanding of the interaction of a particle with the electromagnetic field we need to know the total energy of the system consisting of particle and

More information

Electron Correlation - Methods beyond Hartree-Fock

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

More information

Pure and zero-point vibrational corrections to molecular properties

Pure 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 information

Section 9 Variational Method. Page 492

Section 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 information

2m 2 Ze2. , where δ. ) 2 l,n is the quantum defect (of order one but larger

2m 2 Ze2. , where δ. ) 2 l,n is the quantum defect (of order one but larger PHYS 402, Atomic and Molecular Physics Spring 2017, final exam, solutions 1. Hydrogenic atom energies: Consider a hydrogenic atom or ion with nuclear charge Z and the usual quantum states φ nlm. (a) (2

More information

MOLECULAR SPECTROSCOPY

MOLECULAR SPECTROSCOPY MOLECULAR SPECTROSCOPY First Edition Jeanne L. McHale University of Idaho PRENTICE HALL, Upper Saddle River, New Jersey 07458 CONTENTS PREFACE xiii 1 INTRODUCTION AND REVIEW 1 1.1 Historical Perspective

More information

Methods for Treating Electron Correlation CHEM 430

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

More information

Ab-initio studies of the adiabatic connection in density-functional theory

Ab-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 information

MD simulation: output

MD simulation: output Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square

More information

Quantum 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. 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 information

Chemistry 3502/4502. Final Exam Part I. May 14, 2005

Chemistry 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 information

Density Functional Theory. Martin Lüders Daresbury Laboratory

Density 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 information

QUANTUM MECHANICS. Franz Schwabl. Translated by Ronald Kates. ff Springer

QUANTUM MECHANICS. Franz Schwabl. Translated by Ronald Kates. ff Springer Franz Schwabl QUANTUM MECHANICS Translated by Ronald Kates Second Revised Edition With 122Figures, 16Tables, Numerous Worked Examples, and 126 Problems ff Springer Contents 1. Historical and Experimental

More information

Quantum mechanics can be used to calculate any property of a molecule. The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is,

Quantum 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 information

We also deduced that transformations between Slater determinants are always of the form

We 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 information

Electron Correlation

Electron 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 information

Quantum Mechanical Simulations

Quantum 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 information

An Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine

An 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 information

Session 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education

Session 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 information

Physics 221A Fall 1996 Notes 13 Spins in Magnetic Fields

Physics 221A Fall 1996 Notes 13 Spins in Magnetic Fields Physics 221A Fall 1996 Notes 13 Spins in Magnetic Fields A nice illustration of rotation operator methods which is also important physically is the problem of spins in magnetic fields. The earliest experiments

More information

Density functional calculation of nuclear magnetic resonance chemical shifts

Density functional calculation of nuclear magnetic resonance chemical shifts Density functional calculation of nuclear magnetic resonance chemical shifts Christoph van Wüllen Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany Received 22 August

More information

Chemical bonding in strong magnetic fields

Chemical bonding in strong magnetic fields Trygve Helgaker 1, Mark Hoffmann 1,2 Chemical Kai Lange bonding in 1, strong Alessandro magnetic fields Soncini 1,3, CMS212, and Erik JuneTellgren 24 27 212 1 (CTCC, 1 / 32 Uni Chemical bonding in strong

More information

Phonons and lattice dynamics

Phonons and lattice dynamics Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed

More information

Basis Set for Molecular Orbital Theory

Basis 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 information

Molecular Simulation I

Molecular 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 information

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

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

More information

Introduction to the calculation of molecular properties by response theory

Introduction 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 information

LECTURES ON QUANTUM MECHANICS

LECTURES ON QUANTUM MECHANICS LECTURES ON QUANTUM MECHANICS GORDON BAYM Unitsersity of Illinois A II I' Advanced Bock Progrant A Member of the Perseus Books Group CONTENTS Preface v Chapter 1 Photon Polarization 1 Transformation of

More information

Fundamentals of Spectroscopy for Optical Remote Sensing. Course Outline 2009

Fundamentals of Spectroscopy for Optical Remote Sensing. Course Outline 2009 Fundamentals of Spectroscopy for Optical Remote Sensing Course Outline 2009 Part I. Fundamentals of Quantum Mechanics Chapter 1. Concepts of Quantum and Experimental Facts 1.1. Blackbody Radiation and

More information

I. INTRODUCTION JOURNAL OF CHEMICAL PHYSICS VOLUME 118, NUMBER 3 15 JANUARY ; Electronic mail:

I. INTRODUCTION JOURNAL OF CHEMICAL PHYSICS VOLUME 118, NUMBER 3 15 JANUARY ; Electronic mail: JOURNAL OF CHEMICAL PHYSICS VOLUME 118, NUMBER 3 15 JANUARY 2003 Quasirelativistic theory for the magnetic shielding constant. I. Formulation of Douglas Kroll Hess transformation for the magnetic field

More information

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

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

More information

Lecture 4: Hartree-Fock Theory

Lecture 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 information

Molecules in strong magnetic fields

Molecules in strong magnetic fields Trygve Helgaker 1, Mark Hoffmann 1,2 Kai Molecules Lange in strong 1, Alessandro magnetic fields Soncini 1,3 Nottingham,, and Erik 3th October Tellgren 213 1 (CTCC, 1 / 27 Uni Molecules in strong magnetic

More information

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

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

More information

Quantum Mechanics: Fundamentals

Quantum Mechanics: Fundamentals Kurt Gottfried Tung-Mow Yan Quantum Mechanics: Fundamentals Second Edition With 75 Figures Springer Preface vii Fundamental Concepts 1 1.1 Complementarity and Uncertainty 1 (a) Complementarity 2 (b) The

More information

An Introduction to Hyperfine Structure and Its G-factor

An Introduction to Hyperfine Structure and Its G-factor An Introduction to Hyperfine Structure and Its G-factor Xiqiao Wang East Tennessee State University April 25, 2012 1 1. Introduction In a book chapter entitled Model Calculations of Radiation Induced Damage

More information

Coupled-Cluster Perturbative Triples for Bond Breaking

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

More information

Computational Chemistry I

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

More information

Non-relativistic scattering

Non-relativistic scattering Non-relativistic scattering Contents Scattering theory 2. Scattering amplitudes......................... 3.2 The Born approximation........................ 5 2 Virtual Particles 5 3 The Yukawa Potential

More information

( r) = 1 Z. e Zr/a 0. + n +1δ n', n+1 ). dt ' e i ( ε n ε i )t'/! a n ( t) = n ψ t = 1 i! e iε n t/! n' x n = Physics 624, Quantum II -- Exam 1

( r) = 1 Z. e Zr/a 0. + n +1δ n', n+1 ). dt ' e i ( ε n ε i )t'/! a n ( t) = n ψ t = 1 i! e iε n t/! n' x n = Physics 624, Quantum II -- Exam 1 Physics 624, Quantum II -- Exam 1 Please show all your work on the separate sheets provided (and be sure to include your name) You are graded on your work on those pages, with partial credit where it is

More information

Introduction to Computational Chemistry

Introduction to Computational Chemistry Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry room B430, Chemicum 4th floor vesa.hanninen@helsinki.fi September 3, 2013 Introduction and theoretical backround September

More information

Electron Correlation Methods

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

More information

Density-functional theory in quantum chemistry. Trygve Helgaker. From Quarks to the Nuclear Many-Body Problem

Density-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 information

Lecture 10. The Dirac equation. WS2010/11: Introduction to Nuclear and Particle Physics

Lecture 10. The Dirac equation. WS2010/11: Introduction to Nuclear and Particle Physics Lecture 10 The Dirac equation WS2010/11: Introduction to Nuclear and Particle Physics The Dirac equation The Dirac equation is a relativistic quantum mechanical wave equation formulated by British physicist

More information

Lecture 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 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 information

Beyond the Hartree-Fock Approximation: Configuration Interaction

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

More information

Approximation Methods in QM

Approximation Methods in QM Chapter 3 Approximation Methods in QM Contents 3.1 Time independent PT (nondegenerate)............... 5 3. Degenerate perturbation theory (PT)................. 59 3.3 Time dependent PT and Fermi s golden

More information