Ab initio asymptotic-expansion coefficients for pair energies in Møller-Plesset perturbation theory for atoms
|
|
- Derek Bradley
- 5 years ago
- Views:
Transcription
1 Ab initio asymptotic-expansion coefficients for pair energies in Møller-Plesset perturbation theory for atoms K. JANKOWSKI a, R. SŁUPSKI a, and J. R. FLORES b a Nicholas Copernicus University Toruń, Poland b Universidad de Vigo Vigo, Spain
2 Outline OUTLINE 1 Slow-convergence problem of CI-type approaches. 2 Studies on the asymptotic expansion of second-order energies (He-like systems). 3 Second-order Møller-Plesset PT (MP2) for closed-shell atoms. 4 Asymptotic expansion (AE) of symmetry-adapted MP2 pair energies. 5 Disagreement of theoretical and computed AE coefficients. 6 Derivation of formulae for AE coefficients in the PW/m expansion. 7 Some examples. 8 Conclusions
3 Slow-convergence problem of CI-type approaches Slow convergence of many-electron wave functions in CI-type bases causes major problems in the description of correlation effects.
4 Slow-convergence problem of CI-type approaches Slow convergence of many-electron wave functions in CI-type bases causes major problems in the description of correlation effects. To overcome these problems reliable extrapolation methods are required in calculations of accurate energies. In the last decade much effort has been spent to develop methods of approximating the complete basis set (CBS) limit.
5 Slow-convergence problem of CI-type approaches Slow convergence of many-electron wave functions in CI-type bases causes major problems in the description of correlation effects. To overcome these problems reliable extrapolation methods are required in calculations of accurate energies. In the last decade much effort has been spent to develop methods of approximating the complete basis set (CBS) limit. The functions used in the energy extrapolation procedures are to some extend of empirical nature. However, there are several cases when guidance provided by formal mathematical results has been used.
6 Slow-convergence problem of CI-type approaches Slow convergence of many-electron wave functions in CI-type bases causes major problems in the description of correlation effects. To overcome these problems reliable extrapolation methods are required in calculations of accurate energies. In the last decade much effort has been spent to develop methods of approximating the complete basis set (CBS) limit. The functions used in the energy extrapolation procedures are to some extend of empirical nature. However, there are several cases when guidance provided by formal mathematical results has been used. Unfortunately, there is a real scarcity of such results. Until quite recently, formal information concerning the energy convergence was available only for He-like atomic systems.
7 Studies on the asymptotic expansion of second-order energies (He-like systems) C. Schwartz (1962), pioneered the studies on the rate of convergence correlation energies for atoms. For the He ground state within the 1/Z expansion Perturbation Theory (1/Z PT) the increments to the second-order energy E (2) l determined by orbitals corresponding to the angular momentum quantum number l decrease according to the asymptotic formula E (2) l = a 4 (l + 1/2) 4 + a 6 (l + 1/2) 6 + O [ (l + 1/2) 8],
8 Studies on the asymptotic expansion of second-order energies (He-like systems) C. Schwartz (1962), pioneered the studies on the rate of convergence correlation energies for atoms. For the He ground state within the 1/Z expansion Perturbation Theory (1/Z PT) the increments to the second-order energy E (2) l determined by orbitals corresponding to the angular momentum quantum number l decrease according to the asymptotic formula E (2) l = a 4 (l + 1/2) 4 + a 6 (l + 1/2) 6 + O [ (l + 1/2) 8], The basis-set truncation error, due to omission of functions of l > l max : δe (2) l max = l=l max +1 E (2) l l=l max +1/2 E (2) l dl = a 4 3 (l max + 1) 3 + a 6 5 (l max + 1) 5.
9 Studies on the asymptotic expansion of second-order energies (He-like systems) Kutzelnigg and Morgan (KM) 1992 generalized this result to arbitrary states of He-like systems in 1/Z PT.
10 Studies on the asymptotic expansion of second-order energies (He-like systems) Kutzelnigg and Morgan (KM) 1992 generalized this result to arbitrary states of He-like systems in 1/Z PT. They classified the singlet or triplet states two-electron systems with angular momentum L as natural- or unnatural-parity states: natural parity equal to ( 1) L, e.g., 1s2s 3 S, 2p3p 1 D, unnatural parity equal to ( 1) L+1, e.g., (2p) 2 3 P, 3p3d 1 D,
11 Studies on the asymptotic expansion of second-order energies (He-like systems) Kutzelnigg and Morgan (KM) 1992 generalized this result to arbitrary states of He-like systems in 1/Z PT. They classified the singlet or triplet states two-electron systems with angular momentum L as natural- or unnatural-parity states: natural parity equal to ( 1) L, e.g., 1s2s 3 S, 2p3p 1 D, unnatural parity equal to ( 1) L+1, e.g., (2p) 2 3 P, 3p3d 1 D, There are four types of states of He-like systems: 1 natural-parity singlets, 2 natural-parity triplets, 3 unnatural-parity triplets, 4 unnatural-parity singlets.
12 Studies on the asymptotic expansion of second-order energies (He-like systems) KATO S CORRELATION CUSP CONDITION ( ˆΨ r 12 )r 12 =0 = Ψ(r 12 = 0), where ˆΨ is Ψ averaged over a small sphere about R 12 = 0.
13 Studies on the asymptotic expansion of second-order energies (He-like systems) KATO S CORRELATION CUSP CONDITION ( ˆΨ r 12 )r 12 =0 = Ψ(r 12 = 0), where ˆΨ is Ψ averaged over a small sphere about R 12 = 0. First-order wave functions of He-like systems in 1/Z PT at r 12 = 0 Ψ (1) = α r 12Ψ (0) + O[r 2 12 ] = α (r 12) l P l (cos ϑ 12)Ψ (0) (1, 2) + O[r 2 12 ] (1) l=0 (r 12) l = 1 2l + 3 r l+2 < 1 r l+1 2l 1 > r l <, ϑ 12 angle between r 1 and r 2. r l 1 >
14 Studies on the asymptotic expansion of second-order energies (He-like systems) KATO S CORRELATION CUSP CONDITION ( ˆΨ r 12 )r 12 =0 = Ψ(r 12 = 0), where ˆΨ is Ψ averaged over a small sphere about R 12 = 0. First-order wave functions of He-like systems in 1/Z PT at r 12 = 0 Ψ (1) = α r 12Ψ (0) + O[r 2 12 ] = α (r 12) l P l (cos ϑ 12)Ψ (0) (1, 2) + O[r 2 12 ] (1) l=0 (r 12) l = 1 2l + 3 r l+2 < 1 r l+1 2l 1 > α = 1/2 natural-parity singlets (Kato,1957), r l <, ϑ 12 angle between r 1 and r 2. r l 1 > α = 1/4 natural- and unnatural-parity triplets (Pack, Byers-Brown, 1966), α = 1/6 unnatural-parity singlets (Kutzelnigg-Morgan, 1992).
15 Studies on the asymptotic expansion of second-order energies (He-like systems ASYMPTOTIC EXPANSIONS OF THE SECOND-ORDER ENERGY for He-like systems in 1/Z PT (K-M) 1 natural-parity singlets, E (2) l 2 natural-parity triplets, = a NS 4 (l + 1/2) 4 + a NS 6 (l + 1/2) 6 + O [ (l + 1/2) 8], E (2) l 3 unnatural-parity triplets, E (2) l 4 unnatural-parity singlets. E (2) l = a NT 6 (l + 1/2) 6 + O [ (l + 1/2) 8], = a UT 6 (l + 1/2) 6 + O [ (l + 1/2) 8], = a US 8 (l + 1/2) 8 + O [ (l + 1/2) 10], K-M derived formulae for a NS 4, ans 6, ant 6, a UT 6, a US 8 in terms of radial integrals. No comparison with results for computed energies made so far!
16 Studies on the asymptotic expansion of second-order energies (He-like systems Leading AE coefficients of the pair energies for He-like systems (KM notation) 1 a NS 4 = 3 8 NR(5) + F 0. [ 2 a NT 6 = N 3 a UT 6 = NF 1R (5) + [ 4 a US 8 = N N = (1 + δ ni n j δ li l j ) 1 F p [l 1, l 2] 2 min(l 1,l 2 ) n= F 0R (7) (F 1 G 1 ) R (5) F 1R (7) (F 2 G 2 ) R (5) + [ [n] n(n + 1) Radial integrals: R (5) + = ] p { l1 l 1 n l 2 l 2 L 0 ] ] }( l1 l 1 n ρ +(r, r)r 5 dr, R (7) = )( l2 l 2 n ρ (r, r)r 7 dr, where ρ ±(r, r ) = 1 Rn1 l 2 1 (r)r n2 l 2 (r ) ± R n1 l 1 (r )R n2 l 2 (r) 2 ρ (r, r ) = 2 ρ (r, r ), r 2 R ni l i (r) radial parts of hydrogenlike orbitals. )
17 Second-order Møller-Plesset PT (MP2) for closed-shell atoms Second-order Møller-Plesset perturbation theory (MP2) in Symmetry Adapted Pair Functions (SAPF) formulation H = N h(i) + r 1 ij, H 0 H HF = N f (i) i=1 i<j i=1 h(i) bare nuclear one-electron Hamiltonian N Fock operator: f (i) = h(i) + [J n(i) K n(i)], n=1 J n(i) and K n(i) Coulomb and exchange potentials
18 Second-order Møller-Plesset PT (MP2) for closed-shell atoms Second-order Møller-Plesset perturbation theory (MP2) in Symmetry Adapted Pair Functions (SAPF) formulation H = N h(i) + r 1 ij, H 0 H HF = N f (i) i=1 i<j i=1 h(i) bare nuclear one-electron Hamiltonian N Fock operator: f (i) = h(i) + [J n(i) K n(i)], n=1 J n(i) and K n(i) Coulomb and exchange potentials Perturbation operator: V = H H HF = i<j r 1 ij N n=1 N i=1 [J n(i) K n(i)], Ψ (0) = Φ HF consists of canonical HF spinorbitals f (1)ϕ k (1) = ɛ k ϕ k (1), (k = 1,..., N).
19 Second-order Møller-Plesset PT (MP2) for closed-shell atoms Equation for Ψ (1) (N-electron level) : (H HF E 0 )Ψ (1) + (V E (1) )Ψ HF = 0, Brillouin theorem Ψ (1) consists entirely of double excitations on Φ HF. Ψ (1) (1, 2,..., N) = 2 1/2 Â[u ij (i, j) ϕ k (k)], i<j k i,j  and ˆB N- and 2-electron antisymmetrizers, respectively. First-order spinorbital pair functions u ij - - strongly orthogonal to spinorbitals of Ψ HF : Ω(1, 2)u ij (1, 2) = u ij (1, 2), Ω(1, 2) strong-orthogonality projector Ω(1, 2)[J n(i) K n(i)]ψ HF = 0 Ω(1, 2)V Ψ HF = Ω(1, 2)r 1 12 Ψ HF
20 Second-order Møller-Plesset PT (MP2) for closed-shell atoms Equation for Ψ (1) (N-electron level) : (H HF E 0 )Ψ (1) + (V E (1) )Ψ HF = 0, Brillouin theorem Ψ (1) consists entirely of double excitations on Φ HF. Ψ (1) (1, 2,..., N) = 2 1/2 Â[u ij (i, j) ϕ k (k)], i<j k i,j  and ˆB N- and 2-electron antisymmetrizers, respectively. First-order spinorbital pair functions u ij - - strongly orthogonal to spinorbitals of Ψ HF : Ω(1, 2)u ij (1, 2) = u ij (1, 2), Ω(1, 2) strong-orthogonality projector Ω(1, 2)[J n(i) K n(i)]ψ HF = 0 Ω(1, 2)V Ψ HF = Ω(1, 2)r 1 12 Ψ HF Equations for Ψ (1) (2-electron level) : [f (1) + f (2) ɛ i ɛ j ]u ij (1, 2) + Ω(1, 2)r 1 12 ˆB[ϕ i (1)ϕ j (2)] = 0,
21 Second-order Møller-Plesset PT (MP2) for closed-shell atoms First-order symmetry-adapted pair functions (SAPF): u(t 1, 2) = ˆB [R T (r 1, r 2)Z l i l j (LM L δ 1, δ 2) S(SM S 1, 2)], where T (T, M L, M S ) with T (n i l i, n j l j, L, S), ( ) Z j 1 j 2 j (LM L δ 1, δ 2) 1 j 2 +M [L] 1/2 j1 j 2 L Y m 1 (δ m 1 m 2 M j 1 1)Y m 2 (δ j 2 2). m 1,m 2 ( 1)
22 Second-order Møller-Plesset PT (MP2) for closed-shell atoms First-order symmetry-adapted pair functions (SAPF): u(t 1, 2) = ˆB [R T (r 1, r 2)Z l i l j (LM L δ 1, δ 2) S(SM S 1, 2)], where T (T, M L, M S ) with T (n i l i, n j l j, L, S), ( ) Z j 1 j 2 j (LM L δ 1, δ 2) 1 j 2 +M [L] 1/2 j1 j 2 L Y m 1 (δ m 1 m 2 M j 1 1)Y m 2 (δ j 2 2). Equations for u(t 1, 2): m 1,m 2 ( 1) [f (1) + f (2) ɛ ni l i ɛ nj l j ]u(t 1, 2) + Ω(1, 2)r 1 12 u0(t 1, 2) = 0. u 0(T 1, 2) = ˆB [R ni l i (r 1)R nj l j (r 2)Z l i l j (LM L δ 1, δ 2) S(SM S 1, 2)] R nj l j (r) HF orbitals.
23 Second-order Møller-Plesset PT (MP2) for closed-shell atoms First-order symmetry-adapted pair functions (SAPF): u(t 1, 2) = ˆB [R T (r 1, r 2)Z l i l j (LM L δ 1, δ 2) S(SM S 1, 2)], where T (T, M L, M S ) with T (n i l i, n j l j, L, S), ( ) Z j 1 j 2 j (LM L δ 1, δ 2) 1 j 2 +M [L] 1/2 j1 j 2 L Y m 1 (δ m 1 m 2 M j 1 1)Y m 2 (δ j 2 2). Equations for u(t 1, 2): m 1,m 2 ( 1) [f (1) + f (2) ɛ ni l i ɛ nj l j ]u(t 1, 2) + Ω(1, 2)r 1 12 u0(t 1, 2) = 0. u 0(T 1, 2) = ˆB [R ni l i (r 1)R nj l j (r 2)Z l i l j (LM L δ 1, δ 2) S(SM S 1, 2)] R nj l j (r) HF orbitals. Main partial wave (PW/m) expansion of u(t 1, 2) u(t 1, 2) = u j1j2 (T 1, 2) : j 1 j 2 u j1j2(t 1, 2) = ˆB [R T j 1 j 2 (r 1, r 2)Z j 1 j 2 (LM L δ 1, δ 2) S(SM S 1, 2)].
24 Second-order Møller-Plesset PT (MP2) for closed-shell atoms MP2 energy of an N-electron atom (SAPF formulation) E (2) = T E (2) (T ); Second-order pair energy for the SAP specified by T : E (2) (T ) = u(t ) Ωr 1 12 u 0(T ). PW/m expansion of the second-order SAP energy: E (2) (T ) = j 1j 2 E (2) j 1j 2 (T ), E (2) j 1j 2 (T ) Second-order PW/m energy increment: E (2) j 1j 2 (T ) = u j1j2 (T ) Ωr 1 12 u 0(T )
25 Second-order Møller-Plesset PT (MP2) for closed-shell atoms Equations for the first-order SAPFs He-like systems 1/Z PT [h(1) + h(2) E ni l i E nj l j ]Ψ (1) (T 1, 2) + [r 1 12 E (1) ]Ψ (0) (T 1, 2) = 0, Electron pairs in atoms MP2 PT [f (1) + f (2) ɛ ni l i ɛ nj l j ]u(t 1, 2) + Ω(1, 2)r 1 12 u0(t 1, 2) = 0. Zero-order SAPFs: Ψ 0 (T 1, 2) = ˆB [ R H n i l i (r 1) R H n j l j (r 2)Z l i l j (LM L δ 1, δ 2) S(SM S 1, 2)]. u 0(T 1, 2) = ˆB [R ni l i (r 1)R nj l j (r 2)Z l i l j (LM L δ 1, δ 2) S(SM S 1, 2)]. R H n j l j (r) and R nj l j (r) hydrogenlike and HF radial orbitals, respectively.
26 Second-order Møller-Plesset PT (MP2) for closed-shell atoms Relevant differences between 1/Z-PT and MP-PT equations: 1 h(i) f (i) (no singular terms introduced by J n and K n terms) 2 r 1 12 Ω(1, 2)r 1 12 (Ω(1, 2) modifies, for any atom, only terms with l 3 in the PW expansion of r 1 12 ).
27 Second-order Møller-Plesset PT (MP2) for closed-shell atoms Relevant differences between 1/Z-PT and MP-PT equations: 1 h(i) f (i) (no singular terms introduced by J n and K n terms) 2 r 1 12 Ω(1, 2)r 1 12 (Ω(1, 2) modifies, for any atom, only terms with l 3 in the PW expansion of r 1 12 ). In the vicinity of r 12 = 0 (i.e., for l 3) Ω(1, 2)r 1 12 r 1 12, In vicinity of r 12 = 0, u(t 1, 2) should have an analogous form like Ψ (1), i.e., u(t 1, 2) α r 12 u 0 (T 1, 2) + O[r12 2 ] = ũ l (T 1, 2) + O[r12 2 ] l=0 where ũ l (T 1, 2) α (r 12 ) l P l (cos ϑ 12 ) u 0 (T 1, 2). (2) K-M formulae for the AE coefficients should be applicable after replacing the hydrogenic radial orbitals by the HF ones. Possibility of comparison of theoretical and computed values of AE coefficients
28 Disagreement of theoretical and computed AE coefficients Flores-Słupski-Jankowski-Malinowski (JCP 121,12334(2004)) For Zn symmetry adapted pairs (SAP) PW/m up to l=45! 28 natural-parity singlets 18 natural-parity triplets 7 unnatural-parity triplets 3 unnatural-parity singlets Very accurate AE coefficients a 1 and a 2 have been computed by a specially designed method from accurate E (2) j 1 j 2 (T ) PW/m increments (JCP 124,104107(2006)).
29 Disagreement of theoretical and computed AE coefficients Flores-Słupski-Jankowski-Malinowski (JCP 121,12334(2004)) For Zn symmetry adapted pairs (SAP) PW/m up to l=45! 28 natural-parity singlets 18 natural-parity triplets 7 unnatural-parity triplets 3 unnatural-parity singlets Very accurate AE coefficients a 1 and a 2 have been computed by a specially designed method from accurate E (2) j 1 j 2 (T ) PW/m increments (JCP 124,104107(2006)). We obtained AE coefficients which agree with the formulae: 1 a NS 4 = 3 8 NR(5) + F0. 2 a NT 6 = N [ 5 64 F0R(7) ( a UT 6 = ( 7 3 ) NF1R(5) + [ 4 a US 8 = N ( 9 5 ) F1R(7) + 77 ] ) (F1 G1) R(5) ( 9 2 ] ) (F2 G2) R(5) +
30 Disagreement of theoretical and computed AE coefficients Flores-Słupski-Jankowski-Malinowski (JCP 121,12334(2004)) For Zn symmetry adapted pairs (SAP) PW/m up to l=45! 28 natural-parity singlets 18 natural-parity triplets 7 unnatural-parity triplets 3 unnatural-parity singlets Very accurate AE coefficients a 1 and a 2 have been computed by a specially designed method from accurate E (2) j 1 j 2 (T ) PW/m increments (JCP 124,104107(2006)). We obtained AE coefficients which agree with the formulae: 1 a NS 4 = 3 8 NR(5) + F0. 2 a NT 6 = N [ 5 64 F0R(7) ( a UT 6 = ( 7 3 ) NF1R(5) + [ 4 a US 8 = N ( 9 5 ) F1R(7) + 77 ] ) (F1 G1) R(5) ( 9 2 ] ) (F2 G2) R(5) + WHAT IS THE REASON OF THIS DISAGREEMENT?
31 Derivation of formulae for AE coefficients in the PW/m expansion The MP2 SAP energies in the large-l regime E (2) (T ) = u(t ) Ωr 1 12 u 0(T ) = αr 12 u 0 (T ) }{{} α l=0 (r 12) l P l (cos ϑ 12 ) u 0 (T ) Ωr 1 12 u 0(T ) + O [ (l + 1/2) 8]
32 Derivation of formulae for AE coefficients in the PW/m expansion The MP2 SAP energies in the large-l regime E (2) (T ) = u(t ) Ωr 1 12 u 0(T ) = αr 12 u 0 (T ) }{{} α l=0 (r 12) l P l (cos ϑ 12 ) u 0 (T ) E (2) can be represented by the PW/a expansion: E (2) (T ) = E (2) (T ), l l=0 Ωr 1 12 u 0(T ) + O [ (l + 1/2) 8] [ E (2) l (T ) = α (r 12 ) l P l (cos ϑ 12 ) u 0 (T 1, 2) Ωr }{{} 1 12 u 0(T ) + O (l + 1/2) 8], ũ l (T 1,2) For l 3 we have: E (2) (T ) = ũ l l (T ) r 1 12 u 0(T ),
33 Derivation of formulae for AE coefficients in the PW/m expansion Correspondence of PW/m and PW/a increments E (2) j 1 j 2 (T ) and E (2) l (T ) Consider the angular parts of ũ l (T 1, 2) and u j1 j 2 (T 1, 2). A useful relation: P k (cos ϑ 12) }{{} = l 4π l l Y m l (δ 1 )Y m (δ l 2 ) c w v 1 v 2 (t 1, t 2, L) ( 1) w+l [t 1, t 2, v 1, v 2] 1/2 { v1 v 2 L t 2 t 1 w The energy increment in the PW/a expansion for l 3: Z t 1 t 2 k (LM L δ 1, δ 2) = cv 1 v (t 2 1, t 2, L) Z v 1 v 2 (LM L δ 1, δ 2). v 1 v 2 E (2) (T ) = α (r l 12) l P l (cos ϑ 12)u 0(T ) r 1 12 u0(t ) = }( t1 v 1 w )( t2 v 2 w l cj 1 j (T )E l j 2 1 j 2 (T ), ) j 1 j 2 E l j 1 j 2 (T ) v l j 1 j 2 (T ) r 1 12 u 0(T ) v l j 1 j 2 (T 1, 2) α ˆB [(r 12) l R ni l i (r 1)R nj l j (r 2)Z j 1 j 2 (LM L δ 1, δ 2) S(SM S 1, 2)].
34 Derivation of formulae for AE coefficients in the PW/m expansion The lth PW/a increment is a sum of constituents of several PW/m increments. For l min 3: E (2) j 1j 2 (T ) (j 1 j 2) l min E (2) l (T ), l=l min Rearranging the terms on rhs we obtain: E (2) j 1j 2 (T ) min(j 1+l i,j 2+l j ) l=max( j 1 l i, j 2 l j ) c l j 1j 2 (T )E l j 1j 2 (T ) (l i and l j defined by T ). (3)
35 Some examples EXAMPLE: Unnatural-parity triplet state (np) 2 3 P The only nonzero indices are: c l l 1 l 1 and cl l+1 l+1. PW/a energy increment: E (2) l = c l l 1 l 1E l l 1 l 1 + c l l+1 l+1e l l+1 l+1 The l-th term of PW/a consists of the (l 1)-th, and (l + 1)-th terms in PW/m. PW/m energy increment: E (2) ll c l 1 ll E l 1 ll + c l+1 ll E l+1 ll.
36 Some examples Extrapolation procedure To derive the formulae for the AE coefficients, extract from Ẽ (2) l (T ) and E l j 1 j 2 (T ), the coefficient of the lowest-k term in the (l + 1/2) k expansion. For the PW/a expansion we can to some extend follow the extrapolation procedure developed by KM. Our approach is different from that of Kutzelnigg and Morgan: We consider directly the pair energies E (2) (T ). KM studied the Hylleraas functional E[Φ] = Φ H (0) E (0) Φ + 2 Re Φ V E (1) Ψ (0), which mplies dealing with kinetic energy terms that mask the simple structure of the AE coefficients.
37 Some examples Comparison of ab initio (a HF 1 ) and numerical (a 1) AE coefficients for Ne (in units of E h ). Pair a HF 1 a 1 1s 2 1 S s2s 1 S s2s 3 S s 2 1 S s2p 1 P s2p 3 P s2p 1 P s2p 3 P p 2 1 S p 2 1 D p 2 3 P
38 Summary Summary Relationship between the second-order He-like 1/Z-PT and Møller-Plesset PT (MP2) theory has been discussed.
39 Summary Summary Relationship between the second-order He-like 1/Z-PT and Møller-Plesset PT (MP2) theory has been discussed. The relevance of using two partial-wave expansions: PW/a and PW/m (directly related to the CI expansion), in asymptotic-expansion (AE) studies is recognized.
40 Summary Summary Relationship between the second-order He-like 1/Z-PT and Møller-Plesset PT (MP2) theory has been discussed. The relevance of using two partial-wave expansions: PW/a and PW/m (directly related to the CI expansion), in asymptotic-expansion (AE) studies is recognized. Equations relating the PW/m and PW/a energy increments are derived.
41 Summary Summary Relationship between the second-order He-like 1/Z-PT and Møller-Plesset PT (MP2) theory has been discussed. The relevance of using two partial-wave expansions: PW/a and PW/m (directly related to the CI expansion), in asymptotic-expansion (AE) studies is recognized. Equations relating the PW/m and PW/a energy increments are derived. Similarity of the formulae for the PW/a AE coefficients of the MP2 pair energies with those for the He-like 1/Z-PT expansion is justified.
42 Summary Summary Relationship between the second-order He-like 1/Z-PT and Møller-Plesset PT (MP2) theory has been discussed. The relevance of using two partial-wave expansions: PW/a and PW/m (directly related to the CI expansion), in asymptotic-expansion (AE) studies is recognized. Equations relating the PW/m and PW/a energy increments are derived. Similarity of the formulae for the PW/a AE coefficients of the MP2 pair energies with those for the He-like 1/Z-PT expansion is justified. Formulae for AE coefficients of the PW/m expansion are derived.
43 Summary Summary Relationship between the second-order He-like 1/Z-PT and Møller-Plesset PT (MP2) theory has been discussed. The relevance of using two partial-wave expansions: PW/a and PW/m (directly related to the CI expansion), in asymptotic-expansion (AE) studies is recognized. Equations relating the PW/m and PW/a energy increments are derived. Similarity of the formulae for the PW/a AE coefficients of the MP2 pair energies with those for the He-like 1/Z-PT expansion is justified. Formulae for AE coefficients of the PW/m expansion are derived. Theoretical PW/m AE coefficients showed close agreement with those extracted from computed energy increments.
44 Summary Summary Relationship between the second-order He-like 1/Z-PT and Møller-Plesset PT (MP2) theory has been discussed. The relevance of using two partial-wave expansions: PW/a and PW/m (directly related to the CI expansion), in asymptotic-expansion (AE) studies is recognized. Equations relating the PW/m and PW/a energy increments are derived. Similarity of the formulae for the PW/a AE coefficients of the MP2 pair energies with those for the He-like 1/Z-PT expansion is justified. Formulae for AE coefficients of the PW/m expansion are derived. Theoretical PW/m AE coefficients showed close agreement with those extracted from computed energy increments. Example of applications of the PW/m coefficients is presented.
45 Summary Summary Relationship between the second-order He-like 1/Z-PT and Møller-Plesset PT (MP2) theory has been discussed. The relevance of using two partial-wave expansions: PW/a and PW/m (directly related to the CI expansion), in asymptotic-expansion (AE) studies is recognized. Equations relating the PW/m and PW/a energy increments are derived. Similarity of the formulae for the PW/a AE coefficients of the MP2 pair energies with those for the He-like 1/Z-PT expansion is justified. Formulae for AE coefficients of the PW/m expansion are derived. Theoretical PW/m AE coefficients showed close agreement with those extracted from computed energy increments. Example of applications of the PW/m coefficients is presented.
46 Summary Summary Relationship between the second-order He-like 1/Z-PT and Møller-Plesset PT (MP2) theory has been discussed. The relevance of using two partial-wave expansions: PW/a and PW/m (directly related to the CI expansion), in asymptotic-expansion (AE) studies is recognized. Equations relating the PW/m and PW/a energy increments are derived. Similarity of the formulae for the PW/a AE coefficients of the MP2 pair energies with those for the He-like 1/Z-PT expansion is justified. Formulae for AE coefficients of the PW/m expansion are derived. Theoretical PW/m AE coefficients showed close agreement with those extracted from computed energy increments. Example of applications of the PW/m coefficients is presented. The present results provide an relationship between the Hartree-Fock wave function and the asymptotic behavior of the MP2 correlation energy.
Chemistry 120A 2nd Midterm. 1. (36 pts) For this question, recall the energy levels of the Hydrogenic Hamiltonian (1-electron):
April 6th, 24 Chemistry 2A 2nd Midterm. (36 pts) For this question, recall the energy levels of the Hydrogenic Hamiltonian (-electron): E n = m e Z 2 e 4 /2 2 n 2 = E Z 2 /n 2, n =, 2, 3,... where Ze is
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 information4πε. me 1,2,3,... 1 n. H atom 4. in a.u. atomic units. energy: 1 a.u. = ev distance 1 a.u. = Å
H atom 4 E a me =, n=,,3,... 8ε 0 0 π me e e 0 hn ε h = = 0.59Å E = me (4 πε ) 4 e 0 n n in a.u. atomic units E = r = Z n nao Z = e = me = 4πε = 0 energy: a.u. = 7. ev distance a.u. = 0.59 Å General results
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 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 informationPRACTICE PROBLEMS Give the electronic configurations and term symbols of the first excited electronic states of the atoms up to Z = 10.
PRACTICE PROBLEMS 2 1. Based on your knowledge of the first few hydrogenic eigenfunctions, deduce general formulas, in terms of n and l, for (i) the number of radial nodes in an atomic orbital (ii) the
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 informationTwo-parameter Study of the 1s2s Excited State of He and Li + - Hund's Rule
Two-parameter Study of the 1ss xcited State of He and Li + - Hund's Rule The trial variational wave functions for the 1s 1 s 1 excited state of helium atom and lithium ion are scaled hydrogen 1s and s
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 information2 Electronic structure theory
Electronic structure theory. Generalities.. Born-Oppenheimer approximation revisited In Sec..3 (lecture 3) the Born-Oppenheimer approximation was introduced (see also, for instance, [Tannor.]). We are
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
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 informationIntroduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule. Vesa Hänninen
Introduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule Vesa Hänninen 1 Introduction In this exercise the equilibrium structure and the electronic energy
More 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 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 information(1.1) In particular, ψ( q 1, m 1 ; ; q N, m N ) 2 is the probability to find the first particle
Chapter 1 Identical particles 1.1 Distinguishable particles The Hilbert space of N has to be a subspace H = N n=1h n. Observables Ân of the n-th particle are self-adjoint operators of the form 1 1 1 1
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 informationChemistry 3502/4502. Exam III. All Hallows Eve/Samhain, ) This is a multiple choice exam. Circle the correct answer.
B Chemistry 3502/4502 Exam III All Hallows Eve/Samhain, 2003 1) This is a multiple choice exam. Circle the correct answer. 2) There is one correct answer to every problem. There is no partial credit. 3)
More information2.4. Quantum Mechanical description of hydrogen atom
2.4. Quantum Mechanical description of hydrogen atom Atomic units Quantity Atomic unit SI Conversion Ang. mom. h [J s] h = 1, 05459 10 34 Js Mass m e [kg] m e = 9, 1094 10 31 kg Charge e [C] e = 1, 6022
More informationPHYS852 Quantum Mechanics II, Spring 2010 HOMEWORK ASSIGNMENT 8: Solutions. Topics covered: hydrogen fine structure
PHYS85 Quantum Mechanics II, Spring HOMEWORK ASSIGNMENT 8: Solutions Topics covered: hydrogen fine structure. [ pts] Let the Hamiltonian H depend on the parameter λ, so that H = H(λ). The eigenstates and
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 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 informationAtomic Structure and Atomic Spectra
Atomic Structure and Atomic Spectra Atomic Structure: Hydrogenic Atom Reading: Atkins, Ch. 10 (7 판 Ch. 13) The principles of quantum mechanics internal structure of atoms 1. Hydrogenic atom: one electron
More informationConsequently, the exact eigenfunctions of the Hamiltonian are also eigenfunctions of the two spin operators
VI. SPIN-ADAPTED CONFIGURATIONS A. Preliminary Considerations We have described the spin of a single electron by the two spin functions α(ω) α and β(ω) β. In this Sect. we will discuss spin in more detail
More information(1/2) M α 2 α, ˆTe = i. 1 r i r j, ˆV NN = α>β
Chemistry 26 Spectroscopy Week # The Born-Oppenheimer Approximation, H + 2. Born-Oppenheimer approximation As for atoms, all information about a molecule is contained in the wave function Ψ, which is the
More informationSolutions Final exam 633
Solutions Final exam 633 S.J. van Enk (Dated: June 9, 2008) (1) [25 points] You have a source that produces pairs of spin-1/2 particles. With probability p they are in the singlet state, ( )/ 2, and with
More informationSelf-consistent Field
Chapter 6 Self-consistent Field A way to solve a system of many electrons is to consider each electron under the electrostatic field generated by all other electrons. The many-body problem is thus reduced
More informationConvergence of many-body wavefunction expansions using a plane wave basis: From the homogeneous electron gas to the solid state
Convergence of many-body wavefunction expansions using a plane wave basis: From the homogeneous electron gas to the solid state TCM Electronic Structure Discussion Group James Shepherd (CUC3, Alavi Group)
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationStudy of Ozone in Tribhuvan University, Kathmandu, Nepal. Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal
Study of Ozone in Tribhuvan University, Kathmandu, Nepal Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal 1 Country of the Mt Everest 2 View of the Mt Everest 3 4 5
More informationLecture 9. Hartree Fock Method and Koopman s Theorem
Lecture 9 Hartree Fock Method and Koopman s Theorem Ψ(N) is approximated as a single slater determinant Φ of N orthogonal One electron spin-orbitals. One electron orbital φ i = φ i (r) χ i (σ) χ i (σ)
More informationSpeed of light c = m/s. x n e a x d x = 1. 2 n+1 a n π a. He Li Ne Na Ar K Ni 58.
Physical Chemistry II Test Name: KEY CHEM 464 Spring 18 Chapters 7-11 Average = 1. / 16 6 questions worth a total of 16 points Planck's constant h = 6.63 1-34 J s Speed of light c = 3. 1 8 m/s ħ = h π
More informationLecture 4 Quantum mechanics in more than one-dimension
Lecture 4 Quantum mechanics in more than one-dimension Background Previously, we have addressed quantum mechanics of 1d systems and explored bound and unbound (scattering) states. Although general concepts
More informationarxiv: v3 [cond-mat.other] 14 Aug 2010
Ground state of two electrons on concentric spheres Pierre-François Loos and Peter M. W. Gill Research School of Chemistry, Australian National University, Canberra, Australian Capital Territory, 000,
More information1 Commutators (10 pts)
Final Exam Solutions 37A Fall 0 I. Siddiqi / E. Dodds Commutators 0 pts) ) Consider the operator  = Ĵx Ĵ y + ĴyĴx where J i represents the total angular momentum in the ith direction. a) Express both
More informationSchrödinger equation for central potentials
Chapter 2 Schrödinger equation for central potentials In this chapter we will extend the concepts and methods introduced in the previous chapter for a one-dimensional problem to a specific and very important
More informationQuantum Theory of Many-Particle Systems, Phys. 540
Quantum Theory of Many-Particle Systems, Phys. 540 Questions about organization Second quantization Questions about last class? Comments? Similar strategy N-particles Consider Two-body operators in Fock
More informationMultielectron Atoms.
Multielectron Atoms. Chem 639. Spectroscopy. Spring 00 S.Smirnov Atomic Units Mass m e 1 9.109 10 31 kg Charge e 1.60 10 19 C Angular momentum 1 1.055 10 34 J s Permittivity 4 0 1 1.113 10 10 C J 1 m 1
More informationQuantum Mechanics Solutions
Quantum Mechanics Solutions (a (i f A and B are Hermitian, since (AB = B A = BA, operator AB is Hermitian if and only if A and B commute So, we know that [A,B] = 0, which means that the Hilbert space H
More informationEnergy levels and radiative rates for Ne-like ions from Cu to Ga
Pramana J. Phys. (2017) 89:79 DOI 10.1007/s12043-017-1469-x Indian Academy of Sciences Energy levels and radiative rates for Ne-like ions from Cu to Ga NARENDRA SINGH and SUNNY AGGARWAL Department of Physics,
More informationFine structure in hydrogen - relativistic effects
LNPhysiqueAtomique016 Fine structure in hydrogen - relativistic effects Electron spin ; relativistic effects In a spectrum from H (or from an alkali), one finds that spectral lines appears in pairs. take
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 informationThe Projector Augmented Wave method
The Projector Augmented Wave method Advantages of PAW. The theory. Approximations. Convergence. 1 The PAW method is... What is PAW? A technique for doing DFT calculations efficiently and accurately. An
More informationSchrödinger equation for central potentials
Chapter 2 Schrödinger equation for central potentials In this chapter we will extend the concepts and methods introduced in the previous chapter ifor a one-dimenional problem to a specific and very important
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 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 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 informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
More informationProblem 1: Spin 1 2. particles (10 points)
Problem 1: Spin 1 particles 1 points 1 Consider a system made up of spin 1/ particles. If one measures the spin of the particles, one can only measure spin up or spin down. The general spin state of a
More informationRydberg States of ThIII and UV
Rydberg States of ThIII and UV W. R. Johnson Department of Physics, Notre Dame University June 25, 2003 Abstract We examine the Rydberg spectrum of the two actinide ions Th 2+ and U 5+. The states of interest
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 informationWave function methods for the electronic Schrödinger equation
Wave function methods for the electronic Schrödinger equation Zürich 2008 DFG Reseach Center Matheon: Mathematics in Key Technologies A7: Numerical Discretization Methods in Quantum Chemistry DFG Priority
More informationLecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions
Lecture Models for heavy-ion collisions (Part III: transport models SS06: Dynamical models for relativistic heavy-ion collisions Quantum mechanical description of the many-body system Dynamics of heavy-ion
More informationWorkshop on: ATOMIC STRUCTURE AND TRANSITIONS: THEORY IN USING SUPERSTRUCTURE PROGRAM
Workshop on: ATOMIC STRUCTURE AND TRANSITIONS: THEORY IN USING SUPERSTRUCTURE PROGRAM PROF. SULTANA N. NAHAR Astronomy, Ohio State U, Columbus, Ohio, USA Email: nahar.1@osu.edu http://www.astronomy.ohio-state.edu/
More informationIntro to Nuclear and Particle Physics (5110)
Intro to Nuclear and Particle Physics (5110) March 13, 009 Nuclear Shell Model continued 3/13/009 1 Atomic Physics Nuclear Physics V = V r f r L r S r Tot Spin-Orbit Interaction ( ) ( ) Spin of e magnetic
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 informationFine Structure Calculations of Atomic Data for Ar XVI
Journal of Modern Physics, 2015, 6, 1609-1630 Published Online September 2015 in SciRes. http://www.scirp.org/journal/jmp http://dx.doi.org/10.4236/jmp.2015.611163 Fine Structure Calculations of Atomic
More informationAb-initio Calculations for Forbidden M1/E2 Decay Rates in Ti XIX ion
EJTP 3, No. 11 (2006) 111 122 Electronic Journal of Theoretical Physics Ab-initio Calculations for Forbidden M1/E2 Decay Rates in Ti XIX ion A. Farrag Physics Department,Faculty of Science, Cairo University,
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 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 information0 belonging to the unperturbed Hamiltonian H 0 are known
Time Independent Perturbation Theory D Perturbation theory is used in two qualitatively different contexts in quantum chemistry. It allows one to estimate (because perturbation theory is usually employed
More informationThe Overhauser Instability
The Overhauser Instability Zoltán Radnai and Richard Needs TCM Group ESDG Talk 14th February 2007 Typeset by FoilTEX Introduction Hartree-Fock theory and Homogeneous Electron Gas Noncollinear spins and
More information1 Rayleigh-Schrödinger Perturbation Theory
1 Rayleigh-Schrödinger Perturbation Theory All perturbative techniques depend upon a few simple assumptions. The first of these is that we have a mathematical expression for a physical quantity for which
More informationLecture 4 Quantum mechanics in more than one-dimension
Lecture 4 Quantum mechanics in more than one-dimension Background Previously, we have addressed quantum mechanics of 1d systems and explored bound and unbound (scattering) states. Although general concepts
More information20 The Hydrogen Atom. Ze2 r R (20.1) H( r, R) = h2 2m 2 r h2 2M 2 R
20 The Hydrogen Atom 1. We want to solve the time independent Schrödinger Equation for the hydrogen atom. 2. There are two particles in the system, an electron and a nucleus, and so we can write the Hamiltonian
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 informationCHM Physical Chemistry II Chapter 9 - Supplementary Material. 1. Constuction of orbitals from the spherical harmonics
CHM 3411 - Physical Chemistry II Chapter 9 - Supplementary Material 1. Constuction of orbitals from the spherical harmonics The wavefunctions that are solutions to the time independent Schrodinger equation
More informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More informationSt Hugh s 2 nd Year: Quantum Mechanics II. Reading. Topics. The following sources are recommended for this tutorial:
St Hugh s 2 nd Year: Quantum Mechanics II Reading The following sources are recommended for this tutorial: The key text (especially here in Oxford) is Molecular Quantum Mechanics, P. W. Atkins and R. S.
More informationHelium and two electron atoms
1 Helium and two electron atoms e 2 r 12 e 1 r 2 r 1 +Ze Autumn 2013 Version: 04.12.2013 2 (1) Coordinate system, Schrödinger Equation 3 slides Evaluation of repulsion term 2 slides Radial Integral - details
More informationMany-Electron Atoms. Thornton and Rex, Ch. 8
Many-Electron Atoms Thornton and Rex, Ch. 8 In principle, can now solve Sch. Eqn for any atom. In practice, -> Complicated! Goal-- To explain properties of elements from principles of quantum theory (without
More informationMagnetism of Atoms and Ions. Wulf Wulfhekel Physikalisches Institut, Karlsruhe Institute of Technology (KIT) Wolfgang Gaede Str. 1, D Karlsruhe
Magnetism of Atoms and Ions Wulf Wulfhekel Physikalisches Institut, Karlsruhe Institute of Technology (KIT) Wolfgang Gaede Str. 1, D-76131 Karlsruhe 1 0. Overview Literature J.M.D. Coey, Magnetism and
More informationr 2 dr h2 α = 8m2 q 4 Substituting we find that variational estimate for the energy is m e q 4 E G = 4
Variational calculations for Hydrogen and Helium Recall the variational principle See Chapter 16 of the textbook The variational theorem states that for a Hermitian operator H with the smallest eigenvalue
More informationChemistry 3502/4502. Final Exam Part I. May 14, 2005
Chemistry 3502/4502 Final Exam Part I May 14, 2005 1. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle (e) The
More informationPotential energy, from Coulomb's law. Potential is spherically symmetric. Therefore, solutions must have form
Lecture 6 Page 1 Atoms L6.P1 Review of hydrogen atom Heavy proton (put at the origin), charge e and much lighter electron, charge -e. Potential energy, from Coulomb's law Potential is spherically symmetric.
More informationDepartment of Physics PRELIMINARY EXAMINATION 2015 Part II. Long Questions
Department of Physics PRELIMINARY EXAMINATION 2015 Part II. Long Questions Friday May 15th, 2014, 14-17h Examiners: Prof. J. Cline, Prof. H. Guo, Prof. G. Gervais (Chair), and Prof. D. Hanna INSTRUCTIONS
More informationOrigin 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 informationPhysics 828 Problem Set 7 Due Wednesday 02/24/2010
Physics 88 Problem Set 7 Due Wednesday /4/ 7)a)Consider the proton to be a uniformly charged sphere of radius f m Determine the correction to the s ground state energy 4 points) This is a standard problem
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 informationLecture 4: Hartree-Fock Theory
Lecture 4: Hartree-Fock Theory One determinant to rule them all, One determinant to find them, One determinant to bring them all and in the darkness bind them Second quantization rehearsal The formalism
More informationElectron Correlation Methods
Electron Correlation Methods HF method: electron-electron interaction is replaced by an average interaction E HF c = E 0 E HF E 0 exact ground state energy E HF HF energy for a given basis set HF E c
More informationChemistry 334 Part 2: Computational Quantum Chemistry
Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical
More 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 informationChapter 1 Basic Concepts: Atoms
Chapter 1 Basic Concepts: Atoms CHEM 511 chapter 1 page 1 of 12 What is inorganic chemistry? The periodic table is made of elements, which are made of...? Particle Symbol Mass in amu Charge 1.0073 +1e
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 informationMagnetism in low dimensions from first principles. Atomic magnetism. Gustav Bihlmayer. Gustav Bihlmayer
IFF 10 p. 1 Magnetism in low dimensions from first principles Atomic magnetism Gustav Bihlmayer Institut für Festkörperforschung, Quantum Theory of Materials Gustav Bihlmayer Institut für Festkörperforschung
More informationCh 125a Problem Set 1
Ch 5a Problem Set Due Monday, Oct 5, 05, am Problem : Bra-ket notation (Dirac notation) Bra-ket notation is a standard and convenient way to describe quantum state vectors For example, φ is an abstract
More informationonly two orbitals, and therefore only two combinations to worry about, but things get
131 Lecture 1 It is fairly easy to write down an antisymmetric wavefunction for helium since there are only two orbitals, and therefore only two combinations to worry about, but things get complicated
More informationJoint ICTP-IAEA Workshop on Nuclear Structure Decay Data: Theory and Evaluation August Introduction to Nuclear Physics - 1
2358-19 Joint ICTP-IAEA Workshop on Nuclear Structure Decay Data: Theory and Evaluation 6-17 August 2012 Introduction to Nuclear Physics - 1 P. Van Isacker GANIL, Grand Accelerateur National d'ions Lourds
More informationThe Hydrogen atom. Chapter The Schrödinger Equation. 2.2 Angular momentum
Chapter 2 The Hydrogen atom In the previous chapter we gave a quick overview of the Bohr model, which is only really valid in the semiclassical limit. cf. section 1.7.) We now begin our task in earnest
More informationNew Frontiers in Nuclear Structure Theory
New Frontiers in Nuclear Structure Theory From Realistic Interactions to the Nuclear Chart Robert Roth Institut für Kernphysik Technical University Darmstadt Overview Motivation Nucleon-Nucleon Interactions
More informationDevelopment of atomic theory
Development of atomic theory The chapter presents the fundamentals needed to explain and atomic & molecular structures in qualitative or semiquantitative terms. Li B B C N O F Ne Sc Ti V Cr Mn Fe Co Ni
More informationChapter 9: Multi- Electron Atoms Ground States and X- ray Excitation
Chapter 9: Multi- Electron Atoms Ground States and X- ray Excitation Up to now we have considered one-electron atoms. Almost all atoms are multiple-electron atoms and their description is more complicated
More informationCONTINUUM STATES IN THE SHELL MODEL
CONTINUUM STATES IN THE SHELL MODEL Andrey Shirokov Skobeltsyn Institute of Nuclear Physics Lomonosov Moscow State University Collaborators: J. Vary, P. Maris (Iowa State University) A. Mazur, I. Mazur
More informationQuantum Mechanics FKA081/FIM400 Final Exam 28 October 2015
Quantum Mechanics FKA081/FIM400 Final Exam 28 October 2015 Next review time for the exam: December 2nd between 14:00-16:00 in my room. (This info is also available on the course homepage.) Examinator:
More information1 r 2 sin 2 θ. This must be the case as we can see by the following argument + L2
PHYS 4 3. The momentum operator in three dimensions is p = i Therefore the momentum-squared operator is [ p 2 = 2 2 = 2 r 2 ) + r 2 r r r 2 sin θ We notice that this can be written as sin θ ) + θ θ r 2
More informationProton Neutron Scattering
March 6, 205 Lecture XVII Proton Neutron Scattering Protons and neutrons are both spin /2. We will need to extend our scattering matrix to dimensions to include all the possible spin combinations fof the
More information" = Y(#,$) % R(r) = 1 4& % " = Y(#,$) % R(r) = Recitation Problems: Week 4. a. 5 B, b. 6. , Ne Mg + 15 P 2+ c. 23 V,
Recitation Problems: Week 4 1. Which of the following combinations of quantum numbers are allowed for an electron in a one-electron atom: n l m l m s 2 2 1! 3 1 0 -! 5 1 2! 4-1 0! 3 2 1 0 2 0 0 -! 7 2-2!
More informationHartree, Hartree-Fock and post-hf methods
Hartree, Hartree-Fock and post-hf methods MSE697 fall 2015 Nicolas Onofrio School of Materials Engineering DLR 428 Purdue University nonofrio@purdue.edu 1 The curse of dimensionality Let s consider a multi
More informationNotes on Density Functional Theory
Notes on Density Functional Theory Rocco Martinazzo E-mail: rocco.martinazzo@unimi.it Contents 1 Introduction 1 Density Functional Theory 7 3 The Kohn-Sham approach 11 1 Introduction We consider here a
More information