Molecular Magnetic Properties

Size: px
Start display at page:

Download "Molecular Magnetic Properties"

Transcription

1 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 (ESQC) 2013 Torre Normanna, Sicily, Italy September 8 21, 2013 Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 1 / 67

2 Molecular Magnetism Part 1: The electronic Hamiltonian Hamiltonian mechanics and quantization electromagnetic fields scalar and vector potentials electron spin Part 2: Molecules in an external magnetic field Hamiltonian in an external magnetic field gauge transformations and London orbitals magnetizabilities diamagnetism and paramagnetism induced currents molecules in strong fields Part 3: NMR parameters Zeeman and hyperfine operators nuclear shielding constants indirect nuclear spin spin coupling constants Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 2 / 67

3 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 such 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 (Hamilton s equations of motion) Example: a single particle of mass m in a conservative force field F (q) the Hamiltonian function is constructed from the corresponding scalar potential: H(q, p) = p2 2m + V (q), V (q) F (q) = q Note: Hamilton s equations 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 of motion) = q q Newton s equations are second-order differential equations Hamilton s equations are first-order differential equations the Hamiltonian function is not unique! Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 3 / 67

4 Quantization of a particle in 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 (the 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 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 (CTCC, University of Oslo) Magnetic Properties ESQC12 4 / 67

5 The 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 and B satisfy Maxwell s equations ( ): Note: E = ρ/ε 0 B ε 0 µ 0 E/ t = µ 0 J B = 0 E + B/ t = 0 Coulomb s law Ampère s law with Maxwell s correction Faraday s law of induction when the charge and current densities ρ(r, t) and J(r, t) are known, Maxwell s equations can be solved for E(r, t) and B(r, t) on the other hand, since the charges (particles) are driven by the Lorentz force, ρ(r, t) and J(r, t) are functions of E(r, t) and B(r, t) We shall consider the motion of particles in a given (fixed) electromagnetic field. Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 5 / 67

6 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). They are obtained by solving the first, inhomogeneous pair of Maxwell s equations, which contains ρ and J. Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 6 / 67

7 Gauge transformations The scalar and vector potentials φ and A are not unique. Consider the following transformation of the potentials: φ = φ f } t A f = f (q, t) gauge function of position and time = A + f This gauge transformation of the potentials does not affect the physical fields: E = φ A t = φ + f t A t f t B = (A + f ) = B + f = B = E We are free to choose f (q, t) so as to make φ and A satisfy additional conditions. In the Coulomb gauge, the gauge function is chosen to make vector potential divergenceless: A = 0 Coulomb gauge Note: the Hamiltonian changes in the following manner upon a gauge transformation: H = H z f t However, the equations of motion are unaffected! Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 7 / 67

8 The 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 (CTCC, University of Oslo) Magnetic Properties ESQC12 8 / 67

9 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 only coupled in the presence of an external magnetic field Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 9 / 67

10 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 but we note that, in the nonrelativistic limit, all magnetic fields disappear... We interpret σ by associating an intrinsic angular momentum (spin) with the electron: s = σ/2 Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 10 / 67

11 Molecular 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 (CTCC, University of Oslo) Magnetic Properties ESQC12 11 / 67

12 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 (CTCC, University of Oslo) Magnetic Properties ESQC12 12 / 67

13 Review We have studied the nonrelativistic electronic Hamiltonian: H = H 0 + H (1) + H (2) = H 0 + A (r) p + B (r) s A (r)2 where the vector potential in the Coulomb gauge is not unique but satisfies the relations A (r) = B (r), A (r) = 0 We shall consider uniform external field and the nuclear magnetic moments: A (r) = 1 2 B r O, A K (r) = α 2 M K r K rk 3, 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 We shall first consider external magnetic fields and next nuclear magnetic fields Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 13 / 67

14 Molecular Magnetic Properties Part 1: The electronic Hamiltonian Hamiltonian mechanics and quantization electromagnetic fields scalar and vector potentials electron spin Part 2: Molecules in an external magnetic field Hamiltonian in an external magnetic field gauge transformations and London orbitals magnetizabilities diamagnetism and paramagnetism induced currents molecules in strong fields Part 3: NMR parameters Zeeman and hyperfine operators nuclear shielding constants indirect nuclear spin spin coupling constants Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 14 / 67

15 Review We have studied the nonrelativistic electronic Hamiltonian: H = H 0 + H (1) + H (2) = H 0 + A (r) p + B (r) s A (r)2 where the vector potential in the Coulomb gauge is not unique but satisfies the relations A (r) = B (r), A (r) = 0 We shall consider uniform external field and the nuclear magnetic moments: A (r) = 1 2 B r O, A K (r) = α 2 M K r K rk 3, 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 We shall first consider external magnetic fields and next nuclear magnetic fields Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 15 / 67

16 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 (static) fields 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: 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: L O = r O p 1 2 A2 (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] a change of the origin is a gauge transformation Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 16 / 67

17 Gauge transformation of the 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 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 new wave function must undergo a compensating unitary transformation: Ψ = exp ( if ) Ψ All observable properties such as the electron density are then unaffected: ρ = (Ψ ) Ψ = [Ψ exp( if )] [exp( if )Ψ] = Ψ Ψ = ρ Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 17 / 67

18 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 (CTCC, University of Oslo) Magnetic Properties ESQC12 18 / 67

19 London atomic orbitals The exact wave function transforms in the following manner: Ψ exact G = exp [ i 1 2 B (G O) r] Ψ exact O this behaviour cannot easily be modeled 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 centered 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 (CTCC, University of Oslo) Magnetic Properties ESQC12 19 / 67

20 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 (CTCC, University of Oslo) Magnetic Properties ESQC12 20 / 67

21 Molecular magnetic moments and magnetizabilities Expand the molecular electronic energy in the external magnetic field B: E (B) = E 0 B T M perm 1 2 BT ξb + We have here introduced the permanent magnetic moment and the magnetizability: M perm = de db permanent magnetic moment B=0 describes the first-order change in the energy but vanishes for closed-shell systems ξ = d2 E db 2 molecular a) magnetizabilityb) x B= describes the second-order energy and the first-order induced magnetic 10 moment 0.06 The permanent magnetic moment vanishes for closed-shell systems: 0.04 c.c. ˆΩ 0.02 imaginary c.c. 0 The magnetizability describes the curvature at zero magnetic field: c) d) a) b) x left: diamagnetic dependence on the field (ξ < 0); right: paramagnetic dependence on the field (ξ > 0) c) Trygve Helgaker (CTCC, University of Oslo) d) Magnetic Properties 0.02 ESQC12 21 /

22 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 (CTCC, University of Oslo) Magnetic Properties ESQC12 22 / 67

23 Normal distributions of errors for magnetizabilities Normal distributions of magnetizability errors for 27 molecules in the aug-cc-pcvqz basis relative to CCSD(T)/aug-cc-pCV[TQ]Z values (Lutnæs et al., JCP 131, (2009)) Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 23 / 67

24 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 (CTCC, University of Oslo) Magnetic Properties ESQC12 24 / 67

25 Molecular diamagnetism and paramagnetism The Hamiltonian has paramagnetic and diamagnetic parts: H = H BLz + Bsz B2 (x 2 + y 2 ) para- and diamagnetic parts Second-order perturbation theory gives the magnetizability expression: ξ zz = d2 E db 2 = x 2 + y L z n n L z 0 a) dia- andb) paramagnetic parts 2 x 10 E n n E Most closed-shell molecules are diamagnetic with 0.08 ξ zz < 0 (below left) their energy increases in an applied magnetic field the induced currents oppose the field Some closed-shell systems are paramagnetic with ξ zz > 0 (below right) their energy decreases in a magnetic field relaxation of the wave function lowers the energy, dominating the diamagnetic term RHF calculations of the field dependence of the energy for two closed-shell systems: a) c) b) 14 0 x 10 3 d) Trygve Helgaker (CTCC, University of c) Oslo) Magnetic Properties d) ESQC12 25 / 67

26 Closed-shell paramagnetism explained The Hamiltonian has a linear and quadratic dependence on the magnetic field: H = H BLz + Bsz B2 (x 2 + y 2 ) linear and quadratic B terms Consider the effect of B on the ground state and on the first excited state: the quadratic term raises the energy of both states, creating diabatic energy curves depending on the states, diabatic curve crossings may occur at a given field strength the linear term couples the two states, creating adiabatic states and avoiding crossings A sufficiently strong coupling gives a paramagnetic ground-state with a double minimum Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 26 / 67

27 C 20 in a perpendicular magnetic field All systems become diamagnetic in sufficiently strong fields: The transition occurs at a characteristic stabilizing critical field strength B c B c 0.01 for C 20 (ring conformation) above B c is inversely proportional to the area of the molecule normal to the field we estimate that B c should be observable for C 72 H 72 We may in principle separate such molecules by applying a field gradient Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 27 / 67

28 Induced electron rotation The magnetic field induces a rotation of the electrons about the field direction: the derivative of the energy with the field is half the kinetic angular-momentum operator E (B) = r π 0, π = p + A Paramagnetic closed-shell molecules (here BH): E(B x ) Energy Energy!25.11!25.11!25.12!25.12!25.13!25.13 E(B x )!25.14!25.14 L x (B x ) L (B ) x x!25.15!25.15!25.16! Angular momentum 1 Angular momentum !0.5! L x / / r r!! C C x nuc nuc Nuclear shielding Nuclear shielding in Boron Boron Hydrog Hydrogen ! there is no rotation at the field-free energy maximum: B = 0 the onset of paramagnetic rotation Orbital (against energies the field) reduces the energy HOMO!LUMO for B > 0gap the strongest paramagnetic rotation occurs at the energy inflexion 0.5 point the rotation comes to a halt 0.1 at the stabilizing field strength: B = B c the onset of diamagnetic rotation in the opposite direction increases 0.48 the energy for B > B c Diamagnetic closed-shell molecules: induced diamagnetic rotation!0.1 always increases LUMOthe energy HOMO!(B x )! Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties 0.4 ESQC12 28 / 67! gap (B x ) 0.44 "(B x ) Singlet excitation en

29 Molecules in strong magnetic fields The non-relativistic electronic Hamiltonian (a.u.) in a magnetic field B along the z axis: H = H BLz + Bsz B2 (x 2 + y 2 ) linear and quadratic B terms one atomic unit of B corresponds to T = G Coulomb regime: B 1 a.u. earth-like conditions: Coulomb interactions dominate magnetic interactions are treated perturbatively earth magnetism 10 10, NMR 10 4 ; pulsed laboratory field 10 3 a.u. familiar chemistry Intermediate regime: B 1 a.u. white dwarves: up to about 1 a.u. the Coulomb and magnetic interactions are equally important complicated behaviour resulting from an interplay of linear and quadratic terms exotic chemistry Landau regime: B 1 a.u. neutron stars: a.u. magnetic interactions dominate (Landau levels) Coulomb interactions are treated perturbatively relativity becomes important for B α a.u. alien chemistry Trygve Helgaker (CTCC, University of Oslo) Magnetic Properties ESQC12 29 / 67

30 The H2 molecule in a strong magnetic field I Singlet H2 1σg2 aligns itself with the field and becomes more strongly bound I Triplet H2 1σg 1σu i becomes bound in a perpendicular field orientation B I Re 74 pm 66 pm singlet De 459 kj/mol 594 kj/mol Erot 0 83 kj/mol Re 136 pm triplet De 0 12 kj/mol Erot 0 12 kj/mol Polar plots of E (R, Θ) in the singlet (left) and triplet (right) states I red for high energy; blue for low energy Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 30 / 67

31 Perpendicular paramagnetic bonding The magnetic field modifies the MO energy level diagram, creating a zero-bond-order bond perpendicular orientation: antibonding MO is stabilized, while bonding MO is (less) destabilized parallel orientation: both MOs are unaffected relative to the AOs EêkJmol -1 EêkJmol Σ H 2 1Σ u He HF s A HF 1s B FCI Rêpm Σ g FCI Rêpm In the helium limit, the bonding and antibonding MOs transform into 1s and 2p AOs: lim 1σg = 1s, R 0 { lim R 0 1σ u = 2p 0, 2p 1, all orientations parallel orientation perpendicular orientation (stabilized) Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 31 / 67

32 Molecular Magnetic Properties Part 1: the electronic Hamiltonian Hamiltonian mechanics and quantization electromagnetic fields scalar and vector potentials electron spin Part 2: molecules in an external magnetic field Hamiltonian in an external magnetic field gauge transformations and London orbitals magnetizabilities diamagnetism and paramagnetism induced currents molecules bonding in strong fields Part 3: NMR parameters Zeeman and hyperfine operators nuclear shielding constants indirect nuclear spin spin coupling constants Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 32 / 67

33 Review: magnetic perturbations In atomic units, the molecular Hamiltonian is given by H = H 0 + A i (r i ) p i + B i (r i ) s i i i i }{{}}{{} orbital paramagnetic spin paramagnetic There are two kinds of magnetic perturbation operators: φ i (r i ) + 1 A 2 i 2 (r i ) 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 (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 33 / 67

34 Zeeman interactions The scalar and vector potentials of the uniform (static) fields E and B are given by: E = φ A } { = const φ(r) = E r t B = A = const A(r) = 1 2 B r Electrostatic interaction with a uniform electric field: φ(r i ) = E r i = E d e, d e = i i i Orbital Zeeman interaction with a uniform magnetic field: A p i = 1 B r 2 i p i = 1 2 B L, L = r i p i i i i r i electric dipole operator Spin Zeeman interaction with a uniform magnetic field: B s i = B S, S = s i spin ang. mom. op. i i Total Zeeman interaction with a uniform magnetic field: orbital ang. mom. op. H z = B d m, d m = 1 (L + 2S) magnetic dipole operator 2 Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 34 / 67

35 Orbital hyperfine interaction The nuclear magnetic moment M K at R K sets up the magnetic vector potential A K (r): A K (r) = α 2 M K r K rk 3, r K = r R K Compare this vector potential with that of a uniform magnetic field: A = 1 2 B r O, r O = r O 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. It gives rise to the paramagnetic hyperfine interaction: A K p = M K h PSO K, h PSO K = α 2 r K p rk 3 = α 2 L K rk 3 Compare with the orbital Zeeman interaction: paramagnetic SO (PSO) A p = 1 2 B (r O p) = 1 2 B L O The operator h PSO K is often called the paramagnetic spin orbit (PSO) operator Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 35 / 67

36 Spin hyperfine interactions The nuclear magnetic moment of nucleus K has the vector potential A K (r) = α 2 M K r K rk 3 Taking the curl of this vector potential, we obtain the corresponding magnetic field: B K (r) = A K (r) = 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 (r) thus gives rise to two spin hyperfine interactions: B K (r) 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) these are the Fermi-contact and spin dipole interactions they are triplet operators, which change the spin of the wave function Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 36 / 67

37 Perturbation theory The nonrelativistic electronic Hamiltonian: H = H 0 + H (1) + H (2) = H 0 + A (r) p + B (r) s A (r)2 Second-order Rayleigh Schrödinger perturbation theory: 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 Vector potentials of the uniform external field and the nuclear magnetic moments: A (r) = 1 2 B r O, A K (r) = α 2 M K r K r 3 K Orbital and spin Zeeman interactions with the external magnetic field: H (1) Z = 1 2 B L O + B s Orbital and spin hyperfine interactions with the nuclear magnetic moments: H (1) hf = α 2 M K L K } rk 3 {{ } PSO + 8πα2 3 δ (r K ) M K s }{{} FC + α 2 3(s r K )(r K M K ) (M K s)rk 2 } rk 5 {{ } SD Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 37 / 67

38 Zeeman and hyperfine interactions Zeeman SO SS, SO, OO SO PSO FC+SD hyperfine FC+SD PSO Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 38 / 67

39 Taylor expansion of the 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 (already dealt with) 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 (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 39 / 67

40 First-order molecular properties The first-order properties are expectation values of H (1) Permanent magnetic moment M = 0 (1) H Z 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 (1) H hf 0 = 8πα δ (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 (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 40 / 67

41 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 (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 41 / 67

42 Simulated 200 MHz NMR spectra of vinyllithium 12 C 2 H 3 6 Li experiment RHF MCSCF B3LYP Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 42 / 67

43 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 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 (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 43 / 67

44 Zeeman and hyperfine interactions Zeeman SO SS, SO, OO SO PSO FC+SD hyperfine FC+SD PSO Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 44 / 67

45 Ramsey s expression for the 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 L 2 O n n r 3 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 (B) A (M 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 (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 45 / 67

46 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 (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 46 / 67

47 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 (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 47 / 67

48 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 (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 48 / 67

49 Kohn Sham shielding constants (ppm) HF LDA BLYP B3LYP KT2 CCSD(T) exp. HF F ± 6 H 2O O ± 6 NH 3 N CH 4 C F 2 F N 2 N ± 0.2 CO C ± 0.9 (eq) O ± 6 (eq) Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 49 / 67

50 NMR shieldings: normal distributions of errors LDA 0.06 BLYP 0.06 PBE KT B3LYP B B PBE CAMB3LYP OB3LYP OB OB OPBE OCAMB3LYP CCSD Normal distributions of the errors in NMR shielding constants calculated in the aug-cc-pcvqz basis relative to the CCSD(T)/aug-cc-pCV[TQ]Z benchmark data set for 26 molecules (Teale et al., JCP) Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 50 / 67

51 Mean absolute NMR shielding errors relative to empirical equilibrium values RHF LDA MAE / ppm CC GGA Hybrid OEP-Hybrid H C N O F All 10 0 Mean absolute errors (in ppm) for NMR shielding constants relative to empirical equilibrium values for H (white), C (grey), N (blue), O (red), and F (yellow). The total mean absolute errors over all nuclear types are shown by the purple bars. The DFT methodologies are arranged in the categories LDA, GGA, hybrid and OEP-hybrid. (Teal et al., JCP) Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 51 / 67

52 MAE (Exp.) MAE (Emp. Eq.) NMR shieldings: mean absolute errors relative to experiment MAE / ppm MAE (Exp.) MAE (Emp. Eq.) 0 Mean absolute errors (in ppm) for NMR shielding constants relative to experimental (blue) and empirical equilibrium values (red). The inclusion of vibrational corrections in the empirical equilibrium values leads to a degradation of the quality of the RHF and DFT results but to a notable improvement for the CCSD and CCSD(T) methods. (Teale et al., JCP) Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 52 / 67

53 Direct and indirect nuclear spin spin couplings The last term in the expansion of the molecular electronic energy in B and M K E (B, M) = E BT E (20) B K BT E (11) K M K KL MT K E(02) KL M L + describes the coupling of the nuclear magnetic moments in the presence of electrons There are two distinct contributions to the coupling: the direct and indirect contributions E (02) KL = D KL + K KL The direct coupling occurs by a classical dipole mechanism: D KL = α 2 R 5 ( KL R 2 KL I 3 3R KL R T ) KL a.u. it is anisotropic and vanishes in isotropic media such as gases and liquids The indirect coupling arises from hyperfine interactions with the surrounding electrons: it is exceedingly small: K KL a.u. 1 Hz it does not vanish in isotropic media it gives the fine structure of high-resolution NMR spectra Experimentalists usually work in terms of the (nonreduced) spin spin couplings J KL = h γ K 2π γ L 2π K KL isotope dependent Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 53 / 67

54 Zeeman and hyperfine interactions Zeeman SO SS, SO, OO SO PSO FC+SD hyperfine FC+SD PSO Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 54 / 67

55 Ramsey s expression for indirect nuclear spin spin coupling tensors The indirect nuclear spin spin coupling tensors of a closed-shell system are given by: K KL = d2 E el = 0 2 H dm K dm L M K M H M n n H K M 0 L L E n n E 0 Carrying out the differentiation of the Hamiltonian, we obtain Ramsey s expression: K KL = α 0 4 rk T r LI 3 r K rl T rk 3 r L 3 0 2α 0 r 3 4 K L K n n r 3 L LT L 0 E n n E 0 }{{}}{{} diamagnetic spin orbit (DSO) paramagnetic spin orbit (PSO) 0 8π 3 δ(r K )s + 3r K rk T r K 2 I 3 s n n 8π 3 δ(r L)s T + 3r LrL T r L 2 I 3 s T 0 r 5 K r 5 L 2α 4 E n n E 0 }{{} Fermi contact (FC) and spin dipole (SD) the isotropic FC/FC term often dominates short-range coupling constants the FC/SD and SD/FC terms often dominate the anisotropic part of K KL the orbital contributions (especially DSO) are usually but not invariably small for large internuclear separations, the DSO and PSO contributions cancel Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 55 / 67

56 Relative importance of the contributions to spin spin coupling constants The isotropic indirect spin spin coupling constants can be uniquely decomposed as: J KL = J DSO KL + J PSO KL + J FC KL + JSD KL The spin spin coupling constants are often dominated by the FC term Since the FC term is relatively easy to calculate, it is tempting to ignore the other terms. However, none of the contributions can be a priori neglected (N 2 and CO)! 200 PSO FC PSO FC FC FC FC FC SD FC FC PSO FC SD FC -100 SD H2 HF H2O OH NH3 NH CH4 C2H4 HCN CH CC NC N2 CO C2H2 CC Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 56 / 67

57 Calculation of indirect nuclear spin spin coupling constants The calculation of spin spin coupling constants is a challenging task Spin spin coupling constants depend on many coupling mechanisms: 3 singlet response equations and 7 triplet equations for each nucleus for shieldings, only 3 equations are required, for molecules of all sizes Spin spin coupling constants require a proper description of static correlation the Hartree Fock model fails abysmally MCSCF theory treats static correlation properly but is expensive Spin spin couplings are sensitive to the basis set the FC contribution requires an accurate electron density at the nuclei steep s functions must be included in the basis Spin spin couplings are sensitive to the molecular geometry equilibrium structures must be chosen carefully large vibrational corrections (often 5% 10%) For heavy elements, a relativistic treatment may be necessary. However, there is no need for London orbitals since no external magnetic field is involved. Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 57 / 67

58 Restricted Hartree Fock theory and triplet instabilities The correct description of triplet excitations is important for spin spin coupling constants In restricted Hartree Fock (RHF) theory, triplet excitations are often poorly described upon H 2 dissociation, RHF does not describe the singlet ground state correctly but the lowest triplet state dissociates correctly, leading to triplet instabilities 1 D2p g FCI 1 P1s2p 3 P1s2p 1 1 g 1Σ u 2 1 u 1Σg1Σu R ionic covion 3 u 1Σg1Σu covalent 2 1 g 1Σ g 2 1 S1s 2 g FCI Near such instabilities, the RHF description of spin interactions becomes unphysical C 2H 4/Hz 1 J CC 1 J CH 2 J CH 2 J HH 3 J cis 3 J trans exp RHF CAS B3LYP Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 58 / 67

59 Reduced spin spin coupling constants by wave-function theory RHF CAS RAS SOPPA CCSD CC3 exp vib HF 1 K HF CO 1 K CO N 2 1 K NN H 2O 1 K OH K HH NH 3 1 K NH K HH C 2H 4 1 K CC K CH K CH K HH K cis K tns abs at R e % SOPPA: second-order polarization-propagator approximation Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 59 / 67

60 Reduced spin spin coupling constants by density-functional theory LDA BLYP B3LYP PBE B97-3 RAS exp vib HF 1 K HF CO 1 K CO N 2 1 K NN H 2O 1 K OH K HH NH 3 1 K NH K HH C 2H 4 1 K CC K CH K CH K HH K cis K tns abs at R e % Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 60 / 67

61 Comparison of density-functional and wave-function theory Normal distributions of errors for indirect nuclear spin spin coupling constants for the same molecules as on the previous slides HF LDA BLYP B3LYP CAS RAS SOPPA CCSD Some observations: LDA underestimates only slightly, but has a large standard deviation BLYP reduces the LDA errors by a factor of two B3LYP errors are similar to those of CASSCF The CCSD method is slightly better than the SOPPA method Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 61 / 67

62 The Karplus curve Vicinal (three-bond) spin spin coupling constants depend critically on the dihedral angle: 3 J HH in ethane as a function of the dihedral angle: empirical DFT Good agreement with the (empirically constructed) Karplus curve Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 62 / 67

63 Valinomycin C 54 H 90 N 8 O 18 DFT can be applied to large molecular systems such as valinomycin (168 atoms) there are a total of 7587 spin spin couplings to the carbon atoms in valinomycin below, we have plotted the magnitude of the reduced LDA/6-31G coupling constants on a logarithmic scale, as a function of the internuclear distance: the coupling constants decay in characteristic fashion, which we shall examine most of the indirect couplings beyond 500 pm are small and cannot be detected Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 63 / 67

64 Valinomycin C 54 H 90 N 8 O 18 One-bond spin spin couplings to CH, CO, CN, CC greater than 0.01 Hz Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 64 / 67

65 Valinomycin C 54 H 90 N 8 O 18 Two-bond spin spin couplings to CH, CO, CN, CC greater than 0.01 Hz Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 65 / 67

66 Valinomycin C 54 H 90 N 8 O 18 Three-bond spin spin couplings to CH, CO, CN, CC greater than 0.01 Hz Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 66 / 67

67 Valinomycin C 54 H 90 N 8 O 18 Four-bond spin spin couplings to CH, CO, CN, CC greater than 0.01 Hz Trygve Helgaker (CTCC, University of Oslo) Bonding in magnetic fields ESQC12 67 / 67

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

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

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

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

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

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

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

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

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

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

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

Derivatives and Properties

Derivatives and Properties Derivatives and Properties Trygve Helgaker Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Norway Summer School: Modern Wavefunction Methods in Electronic

More information

Molecular bonding in strong magnetic fields

Molecular bonding in strong magnetic fields Trygve Helgaker 1, Mark Hoffmann 1,2 Kai Molecules Lange in strong 1, Alessandro magnetic fields Soncini 1,3 AMAP,, andzurich, ErikJune Tellgren 1 4 212 1 (CTCC, 1 / 35 Uni Molecular bonding in strong

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

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

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

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

Diamagnetism and Paramagnetism in Atoms and Molecules

Diamagnetism and Paramagnetism in Atoms and Molecules Diamagnetism and Paramagnetism in Atoms and Molecules Trygve Helgaker Alex Borgoo, Maria Dimitrova, Jürgen Gauss, Florian Hampe, Christof Holzer, Wim Klopper, Trond Saue, Peter Schwerdtfeger, Stella Stopkowicz,

More 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

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

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

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

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

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

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

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

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

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

Spin Interactions. Giuseppe Pileio 24/10/2006

Spin Interactions. Giuseppe Pileio 24/10/2006 Spin Interactions Giuseppe Pileio 24/10/2006 Magnetic moment µ = " I ˆ µ = " h I(I +1) " = g# h Spin interactions overview Zeeman Interaction Zeeman interaction Interaction with the static magnetic field

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

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

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

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

Density functional theory in magnetic fields

Density functional theory in magnetic fields Density functional theory in magnetic fields U. E. Ekström, T. Helgaker, S. Kvaal, E. Sagvolden, E. Tellgren, A. M. Teale Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry,

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

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

Quantum Chemistry in Magnetic Fields

Quantum Chemistry in Magnetic Fields Quantum Chemistry in Magnetic Fields Trygve Helgaker Hylleraas Centre of Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Norway 11th Triennial Congress of the World Association

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

Lecture 11: calculation of magnetic parameters, part II

Lecture 11: calculation of magnetic parameters, part II TO DO IS TO BE SOCRATES TO BE IS TO DO SARTRE OO BE DO BE DO SINATRA Lecture 11: calculation of magnetic parameters, part II classification of magnetic perturbations, nuclear quadrupole interaction, J-coupling,

More information

Parameters, Calculation of Nuclear Magnetic Resonance

Parameters, Calculation of Nuclear Magnetic Resonance Parameters, Calculation of Nuclear Magnetic Resonance Cynthia J. Jameson University of Illinois at Chicago, Chicago, IL, USA 1 Introduction 1 1.1 Absolute Shielding Tensor and Nuclear Magnetic Resonance

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

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

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

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

Hyperfine interaction

Hyperfine interaction Hyperfine interaction The notion hyperfine interaction (hfi) comes from atomic physics, where it is used for the interaction of the electronic magnetic moment with the nuclear magnetic moment. In magnetic

More information

QUANTUM CHEMISTRY FOR TRANSITION METALS

QUANTUM CHEMISTRY FOR TRANSITION METALS QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline

More 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

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

Performance of Hartree-Fock and Correlated Methods

Performance of Hartree-Fock and Correlated Methods Chemistry 460 Fall 2017 Dr. Jean M. Standard December 4, 2017 Performance of Hartree-Fock and Correlated Methods Hartree-Fock Methods Hartree-Fock methods generally yield optimized geomtries and molecular

More information

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

Chemistry 4560/5560 Molecular Modeling Fall 2014

Chemistry 4560/5560 Molecular Modeling Fall 2014 Final Exam Name:. User s guide: 1. Read questions carefully and make sure you understand them before answering (if not, ask). 2. Answer only the question that is asked, not a different question. 3. Unless

More 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

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

6 NMR Interactions: Zeeman and CSA

6 NMR Interactions: Zeeman and CSA 6 NMR Interactions: Zeeman and CSA 6.1 Zeeman Interaction Up to this point, we have mentioned a number of NMR interactions - Zeeman, quadrupolar, dipolar - but we have not looked at the nature of these

More information

Electric and magnetic multipoles

Electric and magnetic multipoles Electric and magnetic multipoles Trond Saue Trond Saue (LCPQ, Toulouse) Electric and magnetic multipoles Virginia Tech 2017 1 / 22 Multipole expansions In multipolar gauge the expectation value of the

More information

Exercise 1: Structure and dipole moment of a small molecule

Exercise 1: Structure and dipole moment of a small molecule Introduction to computational chemistry Exercise 1: Structure and dipole moment of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the dipole moment of a small

More 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

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

The hydrogen molecule in magnetic fields: The ground states of the Σ manifold of the parallel configuration

The hydrogen molecule in magnetic fields: The ground states of the Σ manifold of the parallel configuration The hydrogen molecule in magnetic fields: The ground states of the Σ manifold of the parallel configuration T.Detmer, P. Schmelcher, F. K. Diakonos and L. S. Cederbaum Theoretische Chemie, Physikalisch

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

QUANTUM MECHANICS AND MOLECULAR STRUCTURE

QUANTUM MECHANICS AND MOLECULAR STRUCTURE 6 QUANTUM MECHANICS AND MOLECULAR STRUCTURE 6.1 Quantum Picture of the Chemical Bond 6.2 Exact Molecular Orbital for the Simplest Molecule: H + 2 6.3 Molecular Orbital Theory and the Linear Combination

More information

This is a simple input file for the calculation of NMR chemical shieldings for a given molecule using the B3LYP functional and def2-tzvpp basis set:

This is a simple input file for the calculation of NMR chemical shieldings for a given molecule using the B3LYP functional and def2-tzvpp basis set: Computing NMR parameters using ORCA This practical comes with a short lecture on the basics of the computation of NMR parameters using standard electronic structure theory methods. By now you should have

More information

What Do Molecules Look Like?

What Do Molecules Look Like? What Do Molecules Look Like? The Lewis Dot Structure approach provides some insight into molecular structure in terms of bonding, but what about 3D geometry? Recall that we have two types of electron pairs:

More information

ψ s a ˆn a s b ˆn b ψ Hint: Because the state is spherically symmetric the answer can depend only on the angle between the two directions.

ψ s a ˆn a s b ˆn b ψ Hint: Because the state is spherically symmetric the answer can depend only on the angle between the two directions. 1. Quantum Mechanics (Fall 2004) Two spin-half particles are in a state with total spin zero. Let ˆn a and ˆn b be unit vectors in two arbitrary directions. Calculate the expectation value of the product

More information

Preliminary Quantum Questions

Preliminary Quantum Questions Preliminary Quantum Questions Thomas Ouldridge October 01 1. Certain quantities that appear in the theory of hydrogen have wider application in atomic physics: the Bohr radius a 0, the Rydberg constant

More information

Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014

Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde

More information

14. Structure of Nuclei

14. Structure of Nuclei 14. Structure of Nuclei Particle and Nuclear Physics Dr. Tina Potter Dr. Tina Potter 14. Structure of Nuclei 1 In this section... Magic Numbers The Nuclear Shell Model Excited States Dr. Tina Potter 14.

More information

Lecture 9 Electronic Spectroscopy

Lecture 9 Electronic Spectroscopy Lecture 9 Electronic Spectroscopy Molecular Orbital Theory: A Review - LCAO approximaton & AO overlap - Variation Principle & Secular Determinant - Homonuclear Diatomic MOs - Energy Levels, Bond Order

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

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

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

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

H B. θ = 90 o. Lecture notes Part 4: Spin-Spin Coupling. θ θ

H B. θ = 90 o. Lecture notes Part 4: Spin-Spin Coupling. θ θ Lecture notes Part 4: Spin-Spin Coupling F. olger Försterling October 4, 2011 So far, spins were regarded spins isolated from each other. owever, the magnetic moment of nuclear spins also have effect on

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

IV. Electronic Spectroscopy, Angular Momentum, and Magnetic Resonance

IV. Electronic Spectroscopy, Angular Momentum, and Magnetic Resonance IV. Electronic Spectroscopy, Angular Momentum, and Magnetic Resonance The foundation of electronic spectroscopy is the exact solution of the time-independent Schrodinger equation for the hydrogen atom.

More information

CE 530 Molecular Simulation

CE 530 Molecular Simulation 1 CE 530 Molecular Simulation Lecture 14 Molecular Models David A. Kofke Department of Chemical Engineering SUNY Buffalo kofke@eng.buffalo.edu 2 Review Monte Carlo ensemble averaging, no dynamics easy

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

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

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

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

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

Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012

Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data January 2012 2327-3 Joint ICTP-IAEA Workshop on Fusion Plasma Modelling using Atomic and Molecular Data 23-27 January 2012 Qunatum Methods for Plasma-Facing Materials Alain ALLOUCHE Univ.de Provence, Lab.de la Phys.

More information

Students are required to pass a minimum of 15 AU of PAP courses including the following courses:

Students are required to pass a minimum of 15 AU of PAP courses including the following courses: School of Physical and Mathematical Sciences Division of Physics and Applied Physics Minor in Physics Curriculum - Minor in Physics Requirements for the Minor: Students are required to pass a minimum of

More information

CHEM-UA 127: Advanced General Chemistry I

CHEM-UA 127: Advanced General Chemistry I CHEM-UA 7: Advanced General Chemistry I I. LINEAR COMBINATION OF ATOMIC ORBITALS Linear combination of atomic orbitals (LCAO) is a simple method of quantum chemistry that yields a qualitative picture of

More information

(relativistic effects kinetic energy & spin-orbit coupling) 3. Hyperfine structure: ) (spin-spin coupling of e & p + magnetic moments) 4.

(relativistic effects kinetic energy & spin-orbit coupling) 3. Hyperfine structure: ) (spin-spin coupling of e & p + magnetic moments) 4. 4 Time-ind. Perturbation Theory II We said we solved the Hydrogen atom exactly, but we lied. There are a number of physical effects our solution of the Hamiltonian H = p /m e /r left out. We already said

More information

Module 6 1. Density functional theory

Module 6 1. Density functional theory Module 6 1. Density functional theory Updated May 12, 2016 B A DDFT C K A bird s-eye view of density-functional theory Authors: Klaus Capelle G http://arxiv.org/abs/cond-mat/0211443 R https://trac.cc.jyu.fi/projects/toolbox/wiki/dft

More information

Chapter 7. Nuclear Magnetic Resonance Spectroscopy

Chapter 7. Nuclear Magnetic Resonance Spectroscopy Chapter 7 Nuclear Magnetic Resonance Spectroscopy I. Introduction 1924, W. Pauli proposed that certain atomic nuclei have spin and magnetic moment and exposure to magnetic field would lead to energy level

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

H 2 in the minimal basis

H 2 in the minimal basis H 2 in the minimal basis Alston J. Misquitta Centre for Condensed Matter and Materials Physics Queen Mary, University of London January 27, 2016 Overview H 2 : The 1-electron basis. The two-electron basis

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

Introduction to Electronic Structure Theory

Introduction to Electronic Structure Theory CSC/PRACE Spring School in Computational Chemistry 2017 Introduction to Electronic Structure Theory Mikael Johansson http://www.iki.fi/~mpjohans Objective: To get familiarised with the, subjectively chosen,

More information

6.1 Nondegenerate Perturbation Theory

6.1 Nondegenerate Perturbation Theory 6.1 Nondegenerate Perturbation Theory Analytic solutions to the Schrödinger equation have not been found for many interesting systems. Fortunately, it is often possible to find expressions which are analytic

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

Bonding forces and energies Primary interatomic bonds Secondary bonding Molecules

Bonding forces and energies Primary interatomic bonds Secondary bonding Molecules Chapter 2. Atomic structure and interatomic bonding 2.1. Atomic structure 2.1.1.Fundamental concepts 2.1.2. Electrons in atoms 2.1.3. The periodic table 2.2. Atomic bonding in solids 2.2.1. Bonding forces

More information

Spin densities and related quantities in paramagnetic defects

Spin densities and related quantities in paramagnetic defects Spin densities and related quantities in paramagnetic defects Roberto Orlando Dipartimento di Scienze e Tecnologie Avanzate Università del Piemonte Orientale Via G. Bellini 25/G, Alessandria roberto.orlando@unipmn.it

More information