GGA - Perdew, Burke and Ernzerhof
Hi,
I am trying to implement the GGA functional of Perdew, Burke and Ernzerhof
(Phys. Rev. Lett. 77, 3865 (1996)). It would appear that the functional
derivative of the gradient correction for the correlation energy (i.e. the
gradient correction to the correlation potential) diverges for a fully
spin polarized system (zeta = +/-1). The paper cites atomization energies
for hydrogen, so there must be some way around this. Has anyone
implemented this functional and solved this problem?
Cheers,
Andrew
+----------------------------------------------------+
Andrew Horsfield e-mail: horsfield "-at-" fecit.co.uk
FECIT, 2 Longwalk Road, Stockley Park, Uxbridge,
Middlesex UB11 1AB, United Kingdom.
phone: +44-(0)181-606-4653 FAX: +44-(0)181-606-4422
+----------------------------------------------------+