1c Correlation energy and potential as determined via the local 2c density functional of Vosko, Wilk and Nusair (Can. J. of Phys. 58, 3c 1200 (1980)), functional I, the RPA parameterization. 4* 5* $Id$ 6* 7#define FNL_VWN_1_RPA 8#include "xc_vwn_a.fh" 9c 10c Compute the function deltaEc(rs,zeta) function and its derivatives 11c wrt rs and zeta for the spin-unrestricted case - the rest has the 12c same form for all VWN functionals and is handled in the header files. 13c 14 dec_rsz = dec_rs1*fz 15 d1dec_rsz(1) = d1dec_rs1*fz 16 d1dec_rsz(2) = dec_rs1*d1fz 17#ifdef SECOND_DERIV 18 d2dec_rsz(1) = d2dec_rs1*fz 19 d2dec_rsz(2) = d1dec_rs1*d1fz 20 d2dec_rsz(3) = dec_rs1*d2fz 21#endif 22#ifdef THIRD_DERIV 23 d3dec_rsz(1) = d3dec_rs1*fz 24 d3dec_rsz(2) = d2dec_rs1*d1fz 25 d3dec_rsz(3) = d1dec_rs1*d2fz 26 d3dec_rsz(4) = dec_rs1*d3fz 27#endif 28c 29#include "xc_vwn_b.fh" 30c 31#ifndef SECOND_DERIV 32#define SECOND_DERIV 33c 34c Compile source again for the 2nd derivative case 35c 36#include "xc_vwn_1_rpa.F" 37#endif 38#ifndef THIRD_DERIV 39#define THIRD_DERIV 40c 41c Compile source again for the 3rd derivative case 42c 43#include "xc_vwn_1_rpa.F" 44#endif 45