[Pw_forum] correction to functional.f90 (4.1): optx identified as a meta-GGA functional

Latévi Max LAWSON DAKU Max.Lawson at unige.ch
Wed Jul 22 18:40:03 CEST 2009


Hi,

Line 399 of functional.f90:

    ismeta     =  (igcx == 7) .or. (igcx == 6 )

should be changed to

    ismeta     =  (igcx == 7)

That is, the OPTX functional was flagged as a meta-GGA.


I noticed this while trying to generate a pseudopotential
with ld1.x  using the OLYP functional (input key "olyp").

By the way, the following comment is written in functional.90.

  !              "olyp"  = "nox+lyp+optx+blyp"!!! UNTESTED !!!

Does this mean, simply, that the olyp shortname is not well tested?
or something else?


By the way, I noticed that the keyword "blyp" is used for designating
both the blyp XC functional and the lyp gradient correction to
correlation. This perhaps explains that the following error message

     from set_dft_value : error #         1
     two conflicting matching values

was issued when I tried to specify the olyp functional as follows:

    'nox-lyp-optx-blyp'

This is a guess based on the presence of Line 261:

    if (matches ('BLYP', dftout) ) call set_dft_value (igcx, 1)



Finally, I've inserted the following lines in the subroutine dft_name(),
Line 609:

  else if (iexch_==0.and.icorr_==3.and.igcx_==6.and.igcc_==3) then
     shortname_ = 'OLYP'


Best regards,
Max


-- 
***********************************************
Latevi Max LAWSON DAKU
Departement de chimie physique 
Universite de Geneve - Sciences II
30, quai Ernest-Ansermet  
CH-1211 Geneve 4 
Switzerland 

Tel: (41) 22/379 6548 ++ Fax: (41) 22/379 6103
***********************************************



More information about the Pw_forum mailing list