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
   +----------------------------------------------------+