NWChem5需要修改源代码(NWChem6不用改)。打开文件src/nwdft/input_dft/xc_inp.F,找到下面的行(288行左右,依赖于版本): xfac(2)=1.05151d0改为: if (abs(xfac(2)).lt.1.d-4) xfac(2)=1.05151d0 保存后重新编译。 调用O3LYP的命令为: DFT XC vwn_5 0.19 lyp 0.81 Hfexch 0.1161 slater 0.9262 OPTX 0.8133 END 这是O3LYP原始文献所采用的定义。如果用Gaussian 03程序的定义,需要把OPTX的因子从0.8133改为1.164393477 |