Molecular Magnetic Properties
|
|
- Alberta Wells
- 5 years ago
- Views:
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 Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Raman Centre for Atomic, Molecular and Optical
More informationMolecular Magnetism. Magnetic Resonance Parameters. Trygve Helgaker
Molecular Magnetism Magnetic Resonance Parameters Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Laboratoire de Chimie Théorique,
More informationMolecules in Magnetic Fields
Molecules in Magnetic Fields Trygve Helgaker Hylleraas Centre, Department of Chemistry, University of Oslo, Norway and Centre for Advanced Study at the Norwegian Academy of Science and Letters, Oslo, Norway
More informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Hylleraas Centre, Department of Chemistry, University of Oslo, Norway and Centre for Advanced Study at the Norwegian Academy of Science and Letters, Oslo,
More informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry and
More informationMolecular Magnetic Properties. The 11th Sostrup Summer School. Quantum Chemistry and Molecular Properties July 4 16, 2010
1 Molecular Magnetic Properties The 11th Sostrup Summer School Quantum Chemistry and Molecular Properties July 4 16, 2010 Trygve Helgaker Centre for Theoretical and Computational Chemistry, Department
More informationMolecular Magnetic Properties ESQC 07. Overview
1 Molecular Magnetic Properties ESQC 07 Trygve Helgaker Department of Chemistry, University of Oslo, Norway Overview the electronic Hamiltonian in an electromagnetic field external and nuclear magnetic
More informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway European Summer School in Quantum Chemistry
More informationMolecular 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 informationMolecular Magnetism. Molecules in an External Magnetic Field. Trygve Helgaker
Molecular Magnetism Molecules in an External Magnetic Field Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Laboratoire de
More informationMolecules in strong magnetic fields
Molecules in strong magnetic fields Trygve Helgaker, Kai Lange, Alessandro Soncini, and Erik Tellgren Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More informationMolecular electronic structure in strong magnetic fields
Trygve Helgaker 1, Mark Hoffmann 1,2 Kai Molecules Lange in strong 1, Alessandro magnetic fields Soncini 1,3, andcctcc Erik Jackson Tellgren 211 1 (CTCC, 1 / 23 Uni Molecular electronic structure in strong
More informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
More informationMolecules 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 informationChemical 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 informationDerivatives 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 informationMolecular 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 informationDalton Quantum Chemistry Program
1 Quotation from home page: Dalton Quantum Chemistry Program Dalton QCP represents a powerful quantum chemistry program for the calculation of molecular properties with SCF, MP2, MCSCF or CC wave functions.
More informationMolecules in strong magnetic fields
Molecules in strong magnetic fields Trygve Helgaker, Kai Lange, Alessandro Soncini, and Erik Tellgren Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway
More informationBasis sets for electron correlation
Basis sets for electron correlation Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry
More informationT. Helgaker, Department of Chemistry, University of Oslo, Norway. T. Ruden, University of Oslo, Norway. W. Klopper, University of Karlsruhe, Germany
1 The a priori calculation of molecular properties to chemical accuarcy T. Helgaker, Department of Chemistry, University of Oslo, Norway T. Ruden, University of Oslo, Norway W. Klopper, University of Karlsruhe,
More informationDiamagnetism and Paramagnetism in Atoms and Molecules
Diamagnetism and Paramagnetism in Atoms and Molecules Trygve Helgaker Alex Borgoo, Maria Dimitrova, Jürgen Gauss, Florian Hampe, Christof Holzer, Wim Klopper, Trond Saue, Peter Schwerdtfeger, Stella Stopkowicz,
More informationIntroduction to Computational Chemistry
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 informationThe calculation of the universal density functional by Lieb maximization
The calculation of the universal density functional by Lieb maximization Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry,
More informationElectric properties of molecules
Electric properties of molecules For a molecule in a uniform electric fielde the Hamiltonian has the form: Ĥ(E) = Ĥ + E ˆµ x where we assume that the field is directed along the x axis and ˆµ x is the
More informationNMR and IR spectra & vibrational analysis
Lab 5: NMR and IR spectra & vibrational analysis A brief theoretical background 1 Some of the available chemical quantum methods for calculating NMR chemical shifts are based on the Hartree-Fock self-consistent
More informationQuantum chemical modelling of molecular properties - parameters of EPR spectra
Quantum chemical modelling of molecular properties - parameters of EPR spectra EPR ( electric paramagnetic resonance) spectra can be obtained only for open-shell systems, since they rely on transitions
More informationConvergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems
1 Convergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems T. Helgaker Centre for Theoretical and Computational Chemistry, Department of
More informationThe Rigorous Calculation of Molecular Properties to Chemical Accuracy. T. Helgaker, Department of Chemistry, University of Oslo, Norway
1 The Rigorous Calculation of Molecular Properties to Chemical Accuracy T. Helgaker, Department of Chemistry, University of Oslo, Norway A. C. Hennum and T. Ruden, University of Oslo, Norway S. Coriani,
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationAb-initio 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 informationHighly accurate quantum-chemical calculations
1 Highly accurate quantum-chemical calculations T. Helgaker Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, Norway A. C. Hennum and T. Ruden, University
More informationSpin 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 informationMD 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 informationImporting ab-initio theory into DFT: Some applications of the Lieb variation principle
Importing ab-initio theory into DFT: Some applications of the Lieb variation principle Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department
More informationElectron Correlation - Methods beyond Hartree-Fock
Electron Correlation - Methods beyond Hartree-Fock how to approach chemical accuracy Alexander A. Auer Max-Planck-Institute for Chemical Energy Conversion, Mülheim September 4, 2014 MMER Summerschool 2014
More informationTime-independent molecular properties
Time-independent molecular properties Trygve Helgaker Hylleraas Centre, Department of Chemistry, University of Oslo, Norway and Centre for Advanced Study at the Norwegian Academy of Science and Letters,
More informationDensity functional theory 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 informationcharges 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 informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
More informationQuantum 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 informationChem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.
Chem 44 Review for Exam Hydrogenic atoms: The Coulomb energy between two point charges Ze and e: V r Ze r Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative
More informationLecture 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 informationParameters, 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 informationComputational Chemistry I
Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of
More informationBeyond the Hartree-Fock Approximation: Configuration Interaction
Beyond the Hartree-Fock Approximation: Configuration Interaction The Hartree-Fock (HF) method uses a single determinant (single electronic configuration) description of the electronic wavefunction. For
More informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More informationThe adiabatic connection
The adiabatic connection Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Dipartimento di Scienze
More informationHyperfine 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 informationQUANTUM 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 information2m 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 informationThe Accurate Calculation of Molecular Energies and Properties: A Tour of High-Accuracy Quantum-Chemical Methods
1 The Accurate Calculation of Molecular Energies and Properties: A Tour of High-Accuracy Quantum-Chemical Methods T. Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry,
More informationPerformance of Hartree-Fock and Correlated Methods
Chemistry 460 Fall 2017 Dr. Jean M. Standard December 4, 2017 Performance of Hartree-Fock and Correlated Methods Hartree-Fock Methods Hartree-Fock methods generally yield optimized geomtries and molecular
More informationDensity-functional theory in quantum chemistry. Trygve Helgaker. From Quarks to the Nuclear Many-Body Problem
1 Density-functional theory in quantum chemistry Trygve Helgaker Centre for Theoretical and Computational Chemistry, University of Oslo, Norway From Quarks to the Nuclear Many-Body Problem A conference
More informationChemistry 4560/5560 Molecular Modeling Fall 2014
Final Exam Name:. User s guide: 1. Read questions carefully and make sure you understand them before answering (if not, ask). 2. Answer only the question that is asked, not a different question. 3. Unless
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationAn 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 information6 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 informationElectric 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 informationExercise 1: Structure and dipole moment of a small molecule
Introduction to computational chemistry Exercise 1: Structure and dipole moment of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the dipole moment of a small
More informationP. W. Atkins and R. S. Friedman. Molecular Quantum Mechanics THIRD EDITION
P. W. Atkins and R. S. Friedman Molecular Quantum Mechanics THIRD EDITION Oxford New York Tokyo OXFORD UNIVERSITY PRESS 1997 Introduction and orientation 1 Black-body radiation 1 Heat capacities 2 The
More informationThe 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 informationThe 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 information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationQUANTUM 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 informationThis 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 informationWhat 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.
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 informationPreliminary 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 informationAnswers 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 information14. 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 informationLecture 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 informationMOLECULAR 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 informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationI. 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 informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle
More informationH 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 informationQuantum 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 informationIV. 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 informationCE 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 informationMethods for Treating Electron Correlation CHEM 430
Methods for Treating Electron Correlation CHEM 430 Electron Correlation Energy in the Hartree-Fock approximation, each electron sees the average density of all of the other electrons two electrons cannot
More informationAN INTRODUCTION TO QUANTUM CHEMISTRY. Mark S. Gordon Iowa State University
AN INTRODUCTION TO QUANTUM CHEMISTRY Mark S. Gordon Iowa State University 1 OUTLINE Theoretical Background in Quantum Chemistry Overview of GAMESS Program Applications 2 QUANTUM CHEMISTRY In principle,
More informationv(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 informationQuantum mechanics can be used to calculate any property of a molecule. The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is,
Chapter : Molecules Quantum mechanics can be used to calculate any property of a molecule The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is, E = Ψ H Ψ Ψ Ψ 1) At first this seems like
More informationI. 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 informationJoint 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 informationStudents 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 informationCHEM-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.
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 informationModule 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 informationChapter 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 informationLecture 5: More about one- Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory.
Lecture 5: More about one- determinant wave functions Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory. Items from Lecture 4 Could the Koopmans theorem
More informationH 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 informationPhysics 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 informationIntroduction to Electronic Structure Theory
CSC/PRACE Spring School in Computational Chemistry 2017 Introduction to Electronic Structure Theory Mikael Johansson http://www.iki.fi/~mpjohans Objective: To get familiarised with the, subjectively chosen,
More information6.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 informationQUANTUM 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 informationBonding 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 informationSpin 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