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

Similar documents
Chemistry 120A 2nd Midterm. 1. (36 pts) For this question, recall the energy levels of the Hydrogenic Hamiltonian (1-electron):

Basis sets for electron correlation

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

OVERVIEW OF QUANTUM CHEMISTRY METHODS

Exercise 1: Structure and dipole moment of a small molecule

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

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

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

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

2 Electronic structure theory

The 3 dimensional Schrödinger Equation

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

Introduction to computational chemistry Exercise I: Structure and electronic energy of a small molecule. Vesa Hänninen

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

Computational Methods. Chem 561

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

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

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

2.4. Quantum Mechanical description of hydrogen atom

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

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

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

Atomic Structure and Atomic Spectra

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

(1/2) M α 2 α, ˆTe = i. 1 r i r j, ˆV NN = α>β

Solutions Final exam 633

Self-consistent Field

Convergence of many-body wavefunction expansions using a plane wave basis: From the homogeneous electron gas to the solid state

Electron Correlation

Study of Ozone in Tribhuvan University, Kathmandu, Nepal. Prof. S. Gurung Central Department of Physics, Tribhuvan University, Kathmandu, Nepal

Lecture 9. Hartree Fock Method and Koopman s Theorem

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.

Lecture 4 Quantum mechanics in more than one-dimension

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

1 Commutators (10 pts)

Schrödinger equation for central potentials

Quantum Theory of Many-Particle Systems, Phys. 540

Multielectron Atoms.

Quantum Mechanics Solutions

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

Fine structure in hydrogen - relativistic effects

Electron Correlation - Methods beyond Hartree-Fock

The Projector Augmented Wave method

Schrödinger equation for central potentials

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

Introduction to Computational Chemistry

Preliminary Quantum Questions

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

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

Rydberg States of ThIII and UV

Electric properties of molecules

Wave function methods for the electronic Schrödinger equation

Lecture Models for heavy-ion collisions (Part III): transport models. SS2016: Dynamical models for relativistic heavy-ion collisions

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

Intro to Nuclear and Particle Physics (5110)

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

Fine Structure Calculations of Atomic Data for Ar XVI

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

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

Highly accurate quantum-chemical calculations

0 belonging to the unperturbed Hamiltonian H 0 are known

The Overhauser Instability

1 Rayleigh-Schrödinger Perturbation Theory

Lecture 4 Quantum mechanics in more than one-dimension

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

Computational Chemistry I

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

Density Functional Theory. Martin Lüders Daresbury Laboratory

St Hugh s 2 nd Year: Quantum Mechanics II. Reading. Topics. The following sources are recommended for this tutorial:

Helium and two electron atoms

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

Magnetism of Atoms and Ions. Wulf Wulfhekel Physikalisches Institut, Karlsruhe Institute of Technology (KIT) Wolfgang Gaede Str. 1, D Karlsruhe

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

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

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

Department of Physics PRELIMINARY EXAMINATION 2015 Part II. Long Questions

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

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

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

Lecture 4: Hartree-Fock Theory

Electron Correlation Methods

Chemistry 334 Part 2: Computational Quantum Chemistry

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

Chapter 1 Basic Concepts: Atoms

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

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

Ch 125a Problem Set 1

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

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

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

New Frontiers in Nuclear Structure Theory

Development of atomic theory

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

CONTINUUM STATES IN THE SHELL MODEL

Quantum Mechanics FKA081/FIM400 Final Exam 28 October 2015

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

Proton Neutron Scattering

" = 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,

Hartree, Hartree-Fock and post-hf methods

Notes on Density Functional Theory

Transcription:

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 87-100 Toruń, Poland b Universidad de Vigo 36-200 Vigo, Spain 03.06.2006

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

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.

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.

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.

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.

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

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.

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.

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,

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.

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.

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 >

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

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!

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 = 15 128 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=0 5 64 F 0R (7) + 15 256 (F 1 G 1 ) R (5) + 35 384 F 1R (7) + 77 3072 (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 0 0 0 ρ +(r, r)r 5 dr, R (7) = )( l2 l 2 n 0 0 0 0 ρ (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. )

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

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

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

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,

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)

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.

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

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 )

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.

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

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

Disagreement of theoretical and computed AE coefficients Flores-Słupski-Jankowski-Malinowski (JCP 121,12334(2004)) For Zn 2+ 56 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)).

Disagreement of theoretical and computed AE coefficients Flores-Słupski-Jankowski-Malinowski (JCP 121,12334(2004)) For Zn 2+ 56 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) + 15 256 ( 7 3 3 a UT 6 = 15 128 ( 7 3 ) NF1R(5) + [ 4 a US 8 = N 35 384 ( 9 5 ) F1R(7) + 77 ] ) (F1 G1) R(5) + 3072 ( 9 2 ] ) (F2 G2) R(5) +

Disagreement of theoretical and computed AE coefficients Flores-Słupski-Jankowski-Malinowski (JCP 121,12334(2004)) For Zn 2+ 56 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) + 15 256 ( 7 3 3 a UT 6 = 15 128 ( 7 3 ) NF1R(5) + [ 4 a US 8 = N 35 384 ( 9 5 ) F1R(7) + 77 ] ) (F1 G1) R(5) + 3072 ( 9 2 ] ) (F2 G2) R(5) + WHAT IS THE REASON OF THIS DISAGREEMENT?

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]

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

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 0 0 0 )( t2 v 2 w 0 0 0 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)].

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)

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.

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.

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 -0.173114-0.173115 1s2s 1 S -0.019364-0.019370 1s2s 3 S -0.029555-0.029553 2s 2 1 S -0.235741-0.235740 1s2p 1 P -0.030622-0.030638 1s2p 3 P -0.047228-0.047227 2s2p 1 P -0.411775-0.411778 2s2p 3 P -0.308022-0.308017 2p 2 1 S -0.561612-0.561563 2p 2 1 D -0.224645-0.224641 2p 2 3 P -0.819018-0.818899

Summary Summary Relationship between the second-order He-like 1/Z-PT and Møller-Plesset PT (MP2) theory has been discussed.

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.

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.

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.

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.

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.

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.

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.

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.

2024 © sciencedocbox.com Privacy Policy | Terms of Service | Feedback