![]() |
Atomic Reference Data |
|
|
(Eq. 1) |
The electron gas parameter, rs, the spin polarization,
, and the ferromagnetic
and paramagnetic exchange energies,
are defined as
|
|
(Eq. 2) |
|
|
(Eq. 3) |
|
|
(Eq. 4) |
and f(
) is given
by
|
|
(Eq. 5) |
where n is the electron number density (implicitly a function of the
spatial coordinates), and
and
its corresponding spin-up and spin-down components
.
The correlation term is related to the function,
![]() |
(Eq. 6) |
where we have x = rs1/2, X(x) = x2+bx+c, Q = (4c-b2)1/2. The parameters x0, b, and c, given in the table below, are used to create three instances of F, using the table below.
| A | x0 | b | c | ||
| Paramagnetic | cP |
0.031 090 7 | -0.104 98 | 3.727 44 | 12.935 2 |
| Ferromagnetic | cF |
0.015 545 35 | -0.325 00 | 7.060 42 | 18.057 8 |
| Spin Stiffness | c |
-1/(6 2) |
-0.004 758 40 | 1.131 07 | 13.004 5 |
The correlation energy is given by
|
|
(Eq. 7) |
The paramagnetic correlation energy is given by
cP(rs) =
F(rs; A, x0, b,
c) with the four parameters taken from the "Paramagnetic" line in the
table. (Similar definitions hold for
cF
and
c.) The polarization
contribution is given by
|
|
(Eq. 8) |
with the function f as defined above, and where we have
|
|
(Eq. 9) |
and
|
|
(Eq. 10) |
is the difference of the ferromagnetic and paramagnetic correlation energies.
The exchange-correlation potential is given by
|
|
(Eq. 11) |
We use this form in all of the codes in this study. To avoid undetected code bugs, the associated subroutine was recoded independently for one of the codes, although the other three codes shared a common subroutine. This subroutine will be provided upon request to charles.clark@nist.gov.