Ab initio asymptotic-expansion coefficients for pair energies in Møller-Plesset perturbation theory for atoms

Size: px
Start display at page:

Download "Ab initio asymptotic-expansion coefficients for pair energies in Møller-Plesset perturbation theory for atoms"

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

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 information

Basis sets for electron correlation

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

More information

4πε. me 1,2,3,... 1 n. H atom 4. in a.u. atomic units. energy: 1 a.u. = ev distance 1 a.u. = Å

4πε. 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 information

OVERVIEW OF QUANTUM CHEMISTRY METHODS

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

More information

Exercise 1: Structure and dipole moment of a small molecule

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

More information

PRACTICE PROBLEMS Give the electronic configurations and term symbols of the first excited electronic states of the atoms up to Z = 10.

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

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

More information

Two-parameter Study of the 1s2s Excited State of He and Li + - Hund's Rule

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

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

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

More information

2 Electronic structure theory

2 Electronic structure theory Electronic structure theory. Generalities.. Born-Oppenheimer approximation revisited In Sec..3 (lecture 3) the Born-Oppenheimer approximation was introduced (see also, for instance, [Tannor.]). We are

More information

The 3 dimensional Schrödinger Equation

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

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

More information

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

Lecture 5: More about one- Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory.

Lecture 5: More about one- Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory. Lecture 5: More about one- determinant wave functions Final words about the Hartree-Fock theory. First step above it by the Møller-Plesset perturbation theory. Items from Lecture 4 Could the Koopmans theorem

More information

Computational Methods. Chem 561

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

More information

(1.1) In particular, ψ( q 1, m 1 ; ; q N, m N ) 2 is the probability to find the first particle

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

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

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

More information

Chemistry 3502/4502. Exam III. All Hallows Eve/Samhain, ) This is a multiple choice exam. Circle the correct answer.

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

2.4. Quantum Mechanical description of hydrogen atom

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

PHYS852 Quantum Mechanics II, Spring 2010 HOMEWORK ASSIGNMENT 8: Solutions. Topics covered: hydrogen fine structure

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

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

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

More information

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

Chemistry 3502/4502. Final Exam Part I. May 14, 2005 Advocacy chit Chemistry 350/450 Final Exam Part I May 4, 005. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle

More information

Atomic Structure and Atomic Spectra

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

Consequently, the exact eigenfunctions of the Hamiltonian are also eigenfunctions of the two spin operators

Consequently, 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 = α>β

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

Solutions Final exam 633

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

Self-consistent Field

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

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

Electron Correlation

Electron Correlation Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation

More information

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

Lecture 9. Hartree Fock Method and Koopman s Theorem

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

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

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

Lecture 4 Quantum mechanics in more than one-dimension

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

arxiv: v3 [cond-mat.other] 14 Aug 2010

arxiv: 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 information

1 Commutators (10 pts)

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

Schrödinger equation for central potentials

Schrö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 information

Quantum Theory of Many-Particle Systems, Phys. 540

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

Multielectron Atoms.

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

Quantum Mechanics Solutions

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

Energy levels and radiative rates for Ne-like ions from Cu to Ga

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

Fine structure in hydrogen - relativistic effects

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

Electron Correlation - Methods beyond Hartree-Fock

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

More information

The Projector Augmented Wave method

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

Schrödinger equation for central potentials

Schrö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 information

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

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

More information

Introduction to Computational Chemistry

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

More information

Preliminary Quantum Questions

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

More information

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

MO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij ) MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation

More information

Problem 1: Spin 1 2. particles (10 points)

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

Rydberg States of ThIII and UV

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

Electric properties of molecules

Electric properties of molecules Electric properties of molecules For a molecule in a uniform electric fielde the Hamiltonian has the form: Ĥ(E) = Ĥ + E ˆµ x where we assume that the field is directed along the x axis and ˆµ x is the

More information

Wave function methods for the electronic Schrödinger equation

Wave function methods for the electronic Schrödinger equation Wave function methods for the electronic Schrödinger equation Zürich 2008 DFG Reseach Center Matheon: Mathematics in Key Technologies A7: Numerical Discretization Methods in Quantum Chemistry DFG Priority

More information

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

Workshop on: ATOMIC STRUCTURE AND TRANSITIONS: THEORY IN USING SUPERSTRUCTURE PROGRAM

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

Intro to Nuclear and Particle Physics (5110)

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

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

Session 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools

More information

Fine Structure Calculations of Atomic Data for Ar XVI

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

Ab-initio Calculations for Forbidden M1/E2 Decay Rates in Ti XIX ion

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

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

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

More information

Highly accurate quantum-chemical calculations

Highly accurate quantum-chemical calculations 1 Highly accurate quantum-chemical calculations T. Helgaker Centre for Theoretical and Computational Chemistry, Department of Chemistry, University of Oslo, Norway A. C. Hennum and T. Ruden, University

More information

0 belonging to the unperturbed Hamiltonian H 0 are known

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

The Overhauser Instability

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

1 Rayleigh-Schrödinger Perturbation Theory

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

Lecture 4 Quantum mechanics in more than one-dimension

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

20 The Hydrogen Atom. Ze2 r R (20.1) H( r, R) = h2 2m 2 r h2 2M 2 R

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

Computational Chemistry I

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

More information

CHM Physical Chemistry II Chapter 9 - Supplementary Material. 1. Constuction of orbitals from the spherical harmonics

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

Density Functional Theory. Martin Lüders Daresbury Laboratory

Density Functional Theory. Martin Lüders Daresbury Laboratory Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei

More information

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

Helium and two electron atoms

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

Many-Electron Atoms. Thornton and Rex, Ch. 8

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

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

r 2 dr h2 α = 8m2 q 4 Substituting we find that variational estimate for the energy is m e q 4 E G = 4

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

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

Chemistry 3502/4502. Final Exam Part I. May 14, 2005 Chemistry 3502/4502 Final Exam Part I May 14, 2005 1. For which of the below systems is = where H is the Hamiltonian operator and T is the kinetic-energy operator? (a) The free particle (e) The

More information

Potential energy, from Coulomb's law. Potential is spherically symmetric. Therefore, solutions must have form

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

Department of Physics PRELIMINARY EXAMINATION 2015 Part II. Long Questions

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

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

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

More information

Physics 828 Problem Set 7 Due Wednesday 02/24/2010

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

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

Quantum mechanics can be used to calculate any property of a molecule. The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is, Chapter : Molecules Quantum mechanics can be used to calculate any property of a molecule The energy E of a wavefunction Ψ evaluated for the Hamiltonian H is, E = Ψ H Ψ Ψ Ψ 1) At first this seems like

More information

Lecture 4: Hartree-Fock Theory

Lecture 4: Hartree-Fock Theory Lecture 4: Hartree-Fock Theory One determinant to rule them all, One determinant to find them, One determinant to bring them all and in the darkness bind them Second quantization rehearsal The formalism

More information

Electron Correlation Methods

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

More information

Chemistry 334 Part 2: Computational Quantum Chemistry

Chemistry 334 Part 2: Computational Quantum Chemistry Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical

More information

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

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

More information

Chapter 1 Basic Concepts: Atoms

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

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

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

More information

Magnetism in low dimensions from first principles. Atomic magnetism. Gustav Bihlmayer. Gustav Bihlmayer

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

Ch 125a Problem Set 1

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

only two orbitals, and therefore only two combinations to worry about, but things get

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

Joint ICTP-IAEA Workshop on Nuclear Structure Decay Data: Theory and Evaluation August Introduction to Nuclear Physics - 1

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

The Hydrogen atom. Chapter The Schrödinger Equation. 2.2 Angular momentum

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

New Frontiers in Nuclear Structure Theory

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

Development of atomic theory

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

Chapter 9: Multi- Electron Atoms Ground States and X- ray Excitation

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

CONTINUUM STATES IN THE SHELL MODEL

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

Quantum Mechanics FKA081/FIM400 Final Exam 28 October 2015

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

1 r 2 sin 2 θ. This must be the case as we can see by the following argument + L2

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

Proton Neutron Scattering

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

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

Hartree, Hartree-Fock and post-hf methods

Hartree, 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 information

Notes on Density Functional Theory

Notes 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