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 Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry and Molecular Properties July Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 1 / 44

2 Magnetic perturbations In atomic units, the molecular Hamiltonian is given by H = H + 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 2 / 44

3 Molecules in an external magnetic field Part 1: Molecules in an external magnetic field Hamiltonian in an external magnetic field gauge transformations and London orbitals diamagnetism and paramagnetism induced currents Part 2: Magnetic resonance parameters Zeeman and hyperfine operators magnetizabilities nuclear shielding constants indirect nuclear spin spin coupling constants Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 3 / 44

4 Hamiltonian in a uniform magnetic field The nonrelativistic electronic Hamiltonian (implied summation over electrons): H = H + 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 (B) = 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 4 / 44

5 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 5 / 44

6 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 =.2 a.u. in the y direction wave function with gauge origin at O = (,, ) (left) and G = (1,, ) (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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 6 / 44

7 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 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 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 7 / 44

8 London orbitals Let us consider the FCI dissociation of H 2 in a magnetic field full lines: London atomic orbitals dashed lines: AOs with gauge origin between atoms dotted lines: AOs with gauge origin on one of the atoms Without London orbitals, the FCI method is not size extensive in magnetic fields Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 8 / 44

9 Diamagnetism and paramagnetism The Hamiltonian has paramagnetic and diamagnetic parts: H = H BLz + Bsz B2 (x 2 + y 2 ) linear and quadratic B terms Most closed-shell molecules are diamagnetic their energy increases in an applied magnetic field.4 induced currents oppose the field according to Lenz s law Some closed-shell systems are paramagnetic their energy decreases in a magnetic field relaxation of the wave function lowers the energy RHF calculations of the field-dependence for two closed-shell systems: a) a) c) b) x b) d) x left: benzene: diamagnetic dependence on an out-of-plane field, ξ < c) d) right: BH: paramagnetic dependence on a perpendicular field, ξ >.2.2 Trygve Helgaker (CTCC, University of Oslo).4 Molecular Magnetic Properties 12th Sostrup Summer School (212) 9 / 44

10 Transition from para- to diamagnetism However, all closed-shell systems become diamagnetic in sufficiently strong fields: One atomic unit field strength is T highest fields created in laboratories is 1 3 a.u. The transition occurs at a characteristic stabilizing critical field strength B c B c.1 for C 2 (ring conformation) above B c is inversely proportional to the area of the molecule normal to the field B c should be observable for C 72 H 72 Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 1 / 44

11 Closed-shell paramagnetic molecules Ground and (singlet) excited states of BH along the z axis zz = 1sB 2 2σ2 BH 2p2 z, zx = 1sB 2 2σ2 BH 2pz 2px, zy = 1s2 B 2σ2 BH2pz 2py All expectation values increase quadratically in a perpendicular field in the y direction: n H BLy B2 (x 2 + z 2 ) n = En n x 2 + z 2 n B 2 = E n 1 2 χnb2 The zz ground state is coupled to the low-lying zx excited state by this field: zz H BLy B2 (x 2 + z 2 ) xz = 1 zz Ly xz B A paramagnetic ground-state with a double minimum is generated by strong coupling Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 11 / 44

12 Induced electron rotation The magnetic field induces a rotation of the electrons about the field direction: the amount of rotation is the expectation value of the kinetic angular-momentum operator Λ = 2E (B), Λ = r π, π = p + A Paramagnetic closed-shell molecules (here BH): Energy Angular momentum Nuclear shielding Energy 1 Angular momentum Nuclear shielding in! Boron! Boron Hydrog.4 Hydrogen!25.12! ! !25.13!25.14!25.14!25.15!25.15!.2!25.16!25.16!.5! there is no rotation at the field-free energy maximum: B = Orbital energies HOMO!LUMO gap the onset of paramagnetic rotation (against the field) reduces the energy for B > Singlet excitation en the strongest paramagnetic rotation occurs at the energy.5 inflexion point.1.4 the rotation comes to a halt at the stabilizing field.48 strength: B = B c the onset of diamagnetic rotation (with the field) increases the energy for B > B.35 c.46.3!.1 LUMO.44 diamagnetic rotation always increases HOMO the energy according to Lenz s law.25! E(B x ) E(B x ) Diamagnetic closed-shell molecules:!(b x ) Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties.4 12th Sostrup Summer School (212) 12 / 44 L x (B x ) L (B ) x x! gap (B x ) L x / / r r!! C C x nuc nuc "(B x )

13 Expansion of the molecular energy in a magnetic field Expand the molecular electronic energy in the external magnetic induction B: E (B) = E + B T E (1) BT E (2) B + The molecular magnetic moment at B is now given by M mol (B) def de (B) = db = E(1) E (2) B + = M perm + ξb +, where we have introduced the permanent magnetic moment and the magnetizability: M perm = E (1) = de db permanent magnetic moment B= describes the first-order change in the energy but vanishes for closed-shell systems ξ = E (2) = d2 E db 2 molecular magnetizability B= describes the second-order energy and the first-order induced magnetic moment First-order energies for imaginary and triplet perturbations vanish for closed-shell systems: c.c. ˆΩ imaginary c.c. c.c. ˆΩ triplet c.c. Such molecules therefore do not have a permanent magnetic moment. Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 13 / 44

14 The magnetizability The electronic Hamiltonian in a uniform magnetic field: H = H B L O + B s + 1 [ 8 B 2 ro 2 (B r O) 2] The molecular magnetizability of a closed-shell system: ξ = d2 E db 2 = 1 ) r O ro (r T 4 O T r O I L O n n L T O 2 E }{{} n n E diamagnetic term The isotropic part of the diamagnetic term is given by: ξ dia = 1 3 Tr ξ dia = 1 6 xo 2 + y O 2 + O z2 = 1 6 } {{ } paramagnetic term r 2 O Only the orbital Zeeman interaction contributions to the paramagnetic term: S singlet state system surface for 1 S systems (closed-shell atoms), the paramagnetic term vanishes altogether: 1 2 L O 1 S gauge origin at nucleus In most (but not all) systems the diamagnetic term dominates: 1 3 JT 2 RHF exp. diff. H 2O % NH % CH % CO % PH % Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 14 / 44

15 The H2 molecule I I I Earth magnetism is 1 1 a.u., whereas NMR uses 1 3 a.u. In stellar atmospheres, much stronger fields exist Polar plots of the singlet (left) and triplet (right) energy E (R, Θ) at B = 1 a.u I Bond distance Re (pm), orientation Θe ( ), diss. energy De, and rot. barrier E (kj/mol) B. 1. Re Trygve Helgaker (CTCC, University of Oslo) θe singlet De E 83 Re 136 Molecular Magnetic Properties triplet θe De 9 12 E 12 12th Sostrup Summer School (212) 15 / 44

16 Part 2: Magnetic resonance parameters Part 1: Molecules in an external magnetic field Hamiltonian in an external magnetic field gauge transformations and London orbitals diamagnetism and paramagnetism induced currents Part 2: Magnetic resonance parameters Zeeman and hyperfine operators magnetizabilities nuclear shielding constants indirect nuclear spin spin coupling constants Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 16 / 44

17 Magnetic perturbations In atomic units, the molecular Hamiltonian is given by H = H + 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 17 / 44

18 Perturbation theory The nonrelativistic electronic Hamiltonian: H = H + H (1) + H (2) = H + A (r) p + B (r) s A (r)2, B (r) = A (r) Second-order Rayleigh Schrödinger perturbation theory: E (1) = A p + B s E (2) = 1 2 A 2 n A p + B s n n A p + B s E n E Zeeman interactions with uniform external field: A (r) = 1 2 B r O H (1) Z = 1 2 B L O + B s Hyperfine interactions with nuclear magnetic moments (1 8 a.u.) A K (r) = α 2 M K r K r 3 K H (1) hf = α 2 M K L K r 3 K }{{} paramagnetic spin orbit (PSO) + 8πα2 3 δ (r K ) M K s + α 2 3(s r K )(r K M K ) (M K s)rk 2 r }{{} K 5 }{{} Fermi contact (FC) spin dipole (SD) Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 18 / 44

19 Zeeman and hyperfine interactions Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 19 / 44

20 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 + perm. magnetic moments {}}{ B T E (1) BT E (2) B }{{} K magnetizability hyperfine coupling {}}{ K B T E (11) K M K } {{ } shieldings + 1 M T K E(1) K + 1 M T K 2 E(2) 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 Seond-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 ( 1 4 a.u.) 2) the nuclear magnetic moments M K couple weakly (µ µ N 1 8 a.u.) Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 2 / 44

21 First-order molecular properties The first-order properties are expectation values of H (1) Permanent magnetic moment M = (1) H Z = 1 2 L O + s 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 = H (1) hf = 8πα2 3 δ (r K ) s M K + measure spin density at the nucleus important in electron paramagnetic resonance (EPR) recall: there are three hyperfine mechanisms: FC, SD and PSO H (1) hf = 8πα2 δ (r 3 K ) M K s + α 2 3(s r K )(r K M K ) (M K s)rk 2 rk 5 + α 2 M K L K rk 3 Note: there are no first-order Zeeman or hyperfine couplings for closed-shell molecules c.c. ˆΩ imaginary c.c. c.c. ˆΩ triplet c.c. Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 21 / 44

22 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 1 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 22 / 44

23 Simulated 2 MHz NMR spectra of vinyllithium 12 C 2 H 3 6 Li experiment RHF MCSCF B3LYP Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 23 / 44

24 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 (2) B K B T E (11) K M K KL M T K E(2) 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 24 / 44

25 Zeeman and hyperfine interactions Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 25 / 44

26 Ramsey s expression for the nuclear shielding tensors Ramsey s expression for nuclear shielding tensors of a closed-shell system: σ K = d2 E el = 2 H dbdm K B M 2 H B n n H M K K E n n E = α2 ro T r K I 3 r O rk T 2 rk 3 α L 2 O n n r 3 K LT K E n n E }{{}}{{} 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 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 26 / 44

27 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. 12 (24) 6841 Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 27 / 44

28 Calculated and experimental shielding constants (ppm) HF CAS MP2 CCSD CCSD(T) exp. HF F ± 6 (3K) H ±.2 (3K) H 2O O ± 6 (3K) H ±.2 NH 3 N H ± 1. CH 4 C H F 2 F N 2 N ±.2 (3K) CO C ±.9 (eq) O ± 6 (eq) mean mean abs For references and details, see Chem. Rev. 99 (1999) 293. for exp. CO and H 2 O values, see Wasylishen and Bryce, JCP 117 (22) 161 Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 28 / 44

29 Kohn Sham shielding constants (ppm) HF LDA BLYP B3LYP KT2 CCSD(T) exp. HF F ± 6 H ±.2 H 2O O ± 6 H ±.2 NH 3 N H ± 1. CH 4 C H F 2 F N 2 N ±.2 CO C ±.9 (eq) O ± 6 (eq) mean mean abs Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 29 / 44

30 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 (2) B K BT E (11) K M K KL MT K E(2) 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 (2) 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 1 12 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 1 16 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 3 / 44

31 Zeeman and hyperfine interactions Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 31 / 44

32 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 = 2 H dm K dm L M K M 2 H M n n H K M L L E n n E Carrying out the differentiation of the Hamiltonian, we obtain Ramsey s expression: K KL = α 4 rk T r LI 3 r K rl T rk 3 r L 3 2α r 3 4 K L K n n r 3 L LT L E n n E }{{}}{{} diamagnetic spin orbit (DSO) paramagnetic spin orbit (PSO) 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 r 5 K r 5 L 2α 4 E n n E }{{} 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 32 / 44

33 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)! 2 PSO 1 FC PSO FC FC FC FC FC SD FC FC PSO FC SD FC -1 SD H2 HF H2O OH NH3 NH CH4 C2H4 HCN CH CC NC N2 CO C2H2 CC Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 33 / 44

34 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 more generally, the lowest RHF triplet excitations are underestimated 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 34 / 44

35 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 propertly 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% 1%) 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 35 / 44

36 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 36 / 44

37 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 37 / 44

38 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 38 / 44

39 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) Molecular Magnetic Properties 12th Sostrup Summer School (212) 39 / 44

40 Valinomycin C 54 H 9 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 5 pm are small and cannot be detected Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 4 / 44

41 Valinomycin C 54 H 9 N 8 O 18 One-bond spin spin couplings to CH, CO, CN, CC greater than.1 Hz Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 41 / 44

42 Valinomycin C 54 H 9 N 8 O 18 Two-bond spin spin couplings to CH, CO, CN, CC greater than.1 Hz Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 42 / 44

43 Valinomycin C 54 H 9 N 8 O 18 Three-bond spin spin couplings to CH, CO, CN, CC greater than.1 Hz Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 43 / 44

44 Valinomycin C 54 H 9 N 8 O 18 Four-bond spin spin couplings to CH, CO, CN, CC greater than.1 Hz Trygve Helgaker (CTCC, University of Oslo) Molecular Magnetic Properties 12th Sostrup Summer School (212) 44 / 44

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

An introduction to Solid State NMR and its Interactions

An introduction to Solid State NMR and its Interactions An introduction to Solid State NMR and its Interactions From tensor to NMR spectra CECAM Tutorial September 9 Calculation of Solid-State NMR Parameters Using the GIPAW Method Thibault Charpentier - CEA

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

Density Functional Theory Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function

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

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

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

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

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

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

55416 Introduction to computational chemistry Exercise on calculation of excitation energies and second-order molecular properties

55416 Introduction to computational chemistry Exercise on calculation of excitation energies and second-order molecular properties 55416 Introduction to computational chemistry Exercise on calculation of excitation energies and second-order molecular properties Juha Vaara Laboratory of Physical Chemistry, Department of Chemistry,

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

Quantum chemistry and vibrational spectra

Quantum chemistry and vibrational spectra Chapter 3 Quantum chemistry and vibrational spectra This chapter presents the quantum chemical results for the systems studied in this work, FHF (Section 3.) and OHF (Section 3.3). These triatomic anions

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

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

Nuclear magnetic resonance spin spin coupling constants from coupled perturbed density functional theory

Nuclear magnetic resonance spin spin coupling constants from coupled perturbed density functional theory JOURNAL OF CHEMICAL PHYSICS VOLUME 113, NUMBER 9 1 SEPTEMBER 2000 Nuclear magnetic resonance spin spin coupling constants from coupled perturbed density functional theory Vladimír Sychrovský, Jürgen Gräfenstein,

More information

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby

Advanced Electronic Structure Theory Density functional theory. Dr Fred Manby Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)

More 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

Chemistry 431. NC State University. Lecture 17. Vibrational Spectroscopy

Chemistry 431. NC State University. Lecture 17. Vibrational Spectroscopy Chemistry 43 Lecture 7 Vibrational Spectroscopy NC State University The Dipole Moment Expansion The permanent dipole moment of a molecule oscillates about an equilibrium value as the molecule vibrates.

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

Chemistry 2. Lecture 1 Quantum Mechanics in Chemistry

Chemistry 2. Lecture 1 Quantum Mechanics in Chemistry Chemistry 2 Lecture 1 Quantum Mechanics in Chemistry Your lecturers 8am Assoc. Prof Timothy Schmidt Room 315 timothy.schmidt@sydney.edu.au 93512781 12pm Assoc. Prof. Adam J Bridgeman Room 222 adam.bridgeman@sydney.edu.au

More information

Lecture B6 Molecular Orbital Theory. Sometimes it's good to be alone.

Lecture B6 Molecular Orbital Theory. Sometimes it's good to be alone. Lecture B6 Molecular Orbital Theory Sometimes it's good to be alone. Covalent Bond Theories 1. VSEPR (valence shell electron pair repulsion model). A set of empirical rules for predicting a molecular geometry

More information

Chemistry 431. Lecture 23

Chemistry 431. Lecture 23 Chemistry 431 Lecture 23 Introduction The Larmor Frequency The Bloch Equations Measuring T 1 : Inversion Recovery Measuring T 2 : the Spin Echo NC State University NMR spectroscopy The Nuclear Magnetic

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

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

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

Lecture 4: methods and terminology, part II

Lecture 4: methods and terminology, part II So theory guys have got it made in rooms free of pollution. Instead of problems with the reflux, they have only solutions... In other words, experimentalists will likely die of cancer From working hard,

More 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

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

(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

Chapter 5. Molecular Orbitals

Chapter 5. Molecular Orbitals Chapter 5. Molecular Orbitals MO from s, p, d, orbitals: - Fig.5.1, 5.2, 5.3 Homonuclear diatomic molecules: - Fig. 5.7 - Para- vs. Diamagnetic Heteronuclear diatomic molecules: - Fig. 5.14 - ex. CO Hybrid

More information

Non perturbative properties of

Non perturbative properties of Non perturbative properties of molecules in strong magnetic fields Alessandro Soncini INPAC Institute for NanoscalePhysics and Chemistry, University of Leuven, elgium Lb ti Nti ld Ch M éti I t Laboratoire

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

wbt Λ = 0, 1, 2, 3, Eq. (7.63)

wbt Λ = 0, 1, 2, 3, Eq. (7.63) 7.2.2 Classification of Electronic States For all diatomic molecules the coupling approximation which best describes electronic states is analogous to the Russell- Saunders approximation in atoms The orbital

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

Chapter IV: Electronic Spectroscopy of diatomic molecules

Chapter IV: Electronic Spectroscopy of diatomic molecules Chapter IV: Electronic Spectroscopy of diatomic molecules IV.2.1 Molecular orbitals IV.2.1.1. Homonuclear diatomic molecules The molecular orbital (MO) approach to the electronic structure of diatomic

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

Polarised Nucleon Targets for Europe, 2nd meeting, Bochum 2005

Polarised Nucleon Targets for Europe, 2nd meeting, Bochum 2005 Polarised Nucleon Targets for Europe, nd meeting, Bochum Temperature dependence of nuclear spin-lattice relaxations in liquid ethanol with dissolved TEMPO radicals H. Štěpánková, J. Englich, J. Kohout,

More information

one ν im: transition state saddle point

one ν im: transition state saddle point Hypothetical Potential Energy Surface Ethane conformations Hartree-Fock theory, basis set stationary points all ν s >0: minimum eclipsed one ν im: transition state saddle point multiple ν im: hilltop 1

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

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

Unusual long-range spin-spin coupling in fluorinated polyenes: A mechanistic analysis

Unusual long-range spin-spin coupling in fluorinated polyenes: A mechanistic analysis THE JOURNAL OF CHEMICAL PHYSICS 127, 174704 2007 Unusual long-range spin-spin coupling in fluorinated polyenes: A mechanistic analysis Jürgen Gräfenstein a Department of Chemistry, Göteborg University,

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

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

The Physical Basis of the NMR Experiment

The Physical Basis of the NMR Experiment The Physical Basis of the NMR Experiment 1 Interaction of Materials with Magnetic Fields F F S N S N Paramagnetism Diamagnetism 2 Microscopic View: Single Spins an electron has mass and charge in addition

More information

HONOUR SCHOOL OF NATURAL SCIENCE. Final Examination GENERAL PHYSICAL CHEMISTRY I. Answer FIVE out of nine questions

HONOUR SCHOOL OF NATURAL SCIENCE. Final Examination GENERAL PHYSICAL CHEMISTRY I. Answer FIVE out of nine questions HONOUR SCHOOL OF NATURAL SCIENCE Final Examination GENERAL PHYSICAL CHEMISTRY I Monday, 12 th June 2000, 9.30 a.m. - 12.30 p.m. Answer FIVE out of nine questions The numbers in square brackets indicate

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

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

Electronic Spectra of Complexes

Electronic Spectra of Complexes Electronic Spectra of Complexes Interpret electronic spectra of coordination compounds Correlate with bonding Orbital filling and electronic transitions Electron-electron repulsion Application of MO theory

More information

The symmetry properties & relative energies of atomic orbitals determine how they react to form molecular orbitals. These molecular orbitals are then

The symmetry properties & relative energies of atomic orbitals determine how they react to form molecular orbitals. These molecular orbitals are then 1 The symmetry properties & relative energies of atomic orbitals determine how they react to form molecular orbitals. These molecular orbitals are then filled with the available electrons according to

More information

Origin of Chemical Shifts BCMB/CHEM 8190

Origin of Chemical Shifts BCMB/CHEM 8190 Origin of Chemical Shifts BCMB/CHEM 8190 Empirical Properties of Chemical Shift υ i (Hz) = γb 0 (1-σ i ) /2π The Larmor frequencies of nuclei depend on the electronic structure of the molecule and the

More information

Origin of the first Hund rule in He-like atoms and 2-electron quantum dots

Origin of the first Hund rule in He-like atoms and 2-electron quantum dots in He-like atoms and 2-electron quantum dots T Sako 1, A Ichimura 2, J Paldus 3 and GHF Diercksen 4 1 Nihon University, College of Science and Technology, Funabashi, JAPAN 2 Institute of Space and Astronautical

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

Origin of Chemical Shifts BCMB/CHEM 8190

Origin of Chemical Shifts BCMB/CHEM 8190 Origin of Chemical Shifts BCMB/CHEM 8190 Empirical Properties of Chemical Shift υ i (Hz) = γb 0 (1-σ i ) /2π σ i, shielding constant dependent on electronic structure, is ~ 10-6. Measurements are made

More information