20.02.2013 Views

Atomic Structure Theory

Atomic Structure Theory

Atomic Structure Theory

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

230 7 Introduction to MBPT<br />

From the Schrödinger equation<br />

one easily shows<br />

(H0 + V )Ψl = ElΨl,<br />

HeffΨ (0)<br />

l<br />

(0)<br />

= ElΨ l , (7.113)<br />

where<br />

Heff = PH0P + PVΩP (7.114)<br />

is an effective Hamiltonian defined within the model space that gives the exact<br />

energy. A perturbation expansion of Heff leads to a perturbation expansion<br />

for the energy.<br />

7.6.2 First-Order Perturbation <strong>Theory</strong><br />

The first-order approximation to Heff is PH0P + PVP. Substituting into<br />

(7.113) leads immediately to the d-dimensional eigenvalue equation<br />

<br />

〈Φl|H0 + V |Φk〉ck = E (1)<br />

l cl. (7.115)<br />

k<br />

Note that the first-order energy above corresponds to the sum of the zerothand<br />

first-order energy for systems with one valence electron. There is no well<br />

defined counterpart of the zeroth-order energy in cases where the model space<br />

consists of non-degenerate states.<br />

In the magnesium example, we let the index k in (7.115) represent the<br />

two-particle state vw(JM) and the index l represent the state xy(JM). The<br />

eigenvalue equation then becomes<br />

<br />

where<br />

v≤w<br />

x≤y<br />

Vxy,vw = ηxyηvw<br />

[(ɛv + ɛw)δvxδwy + Vxy,vw] cvw = Ecxy,<br />

<br />

<br />

(−1)<br />

L<br />

J+L+jy+jv<br />

<br />

jx jy J<br />

TL(xyvw)<br />

jw jv L<br />

+(−1) L+jy+jv<br />

<br />

jx jy J<br />

TL(xywv) . (7.116)<br />

jv jw L<br />

In the relativistic case, TL(rskl) =XL(rskl)+BL(rskl)isthesumofCoulomb<br />

(7.38) and Breit (7.92) interaction matrix elements. In first order, we are<br />

led to a d-dimension two-particle CI calculation, essentially identical to that<br />

described in Section 7.5.1. The basic orbitals here, however, are obtained in<br />

the V (N−2) HF potential of the ionic core.<br />

As an illustration, we consider the magnesiumlike ion P +3 . The diagonal<br />

terms in the first-order Hamiltonian matrix are dominated by sums ɛv + ɛw

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!