请教自设泛函HFE参数的问题

sunnyyg

大家好，我重复文献的G09计算，结果总不对，请大家看看下面这段文献描述怎么设置。
“Among the methods used, TD-LC-ωPBE/6-31+G(d)//mPWPW91α50(*with 50% HFE values*)/6-31G(d) in a solvent effect with ω = 0.20 Bohr-
1 gave the best accuracy of <0.03 eV for the analogues with spatial overlap of 0.42 e ΛHL e 0.52.”

我的计算关键词是这样的：
For GS opt&freq: #p opt mpwpw91/6-31G(d) IOp(3/130=05000) scrf=(cpcm,solvent=ethanol) pop=(nbo,savenbo) geom=connectivity

For Excited State: #p TD-LC-wPBE/6-31G(d) iop(3/107=0200000000) scrf=(cpcm,solvent=thf) geom=connectivity

是不是我把参数设错了，或者漏了什么。结果对不上，差30%。

谢谢大家。

sobereva

杂化泛函调整HF成份不是3/130。opt freq应当是mPWPW91 IOp(3/76=1000005000) IOp(3/77=0500005000)

激发态计算应当是LC-wPBE IOp(3/107=0200000000) IOp(3/108=0200000000)

激发态计算文中用了弥散，你没用弥散。

sunnyyg

sobereva 发表于 2015-11-18 03:56
杂化泛函调整HF成份不是3/130。opt freq应当是mPWPW91 IOp(3/76=1000005000) IOp(3/77=0500005000)

激发 ...

非常感谢，我试试。

sunnyyg

sobereva 发表于 2015-11-18 03:56
杂化泛函调整HF成份不是3/130。opt freq应当是mPWPW91 IOp(3/76=1000005000) IOp(3/77=0500005000)

激发 ...

你好，你给的建议果然算出来了，实在是佩服。IOp(3/76=1000005000) 是指DFT交换为1，HF交换为0.5，为什么DFT交换要设为1呢；IOp(3/77=0500005000)是指局域和非局域交换各0.5，为什么要设置这个呢？

我的问题可能有点蠢或者是偷懒，还是想请你点拨一下，呵呵，谢谢啦。轻拍。

sobereva

sunnyyg 发表于 2015-11-18 22:24
你好，你给的建议果然算出来了，实在是佩服。IOp(3/76=1000005000) 是指DFT交换为1，HF交换为0.5，为什么 ...

若记为IOp(3/76=ab)、IOp(3/77=cd)、IOp(3/78=ef)，则实际形式是
a*[ d*E_X_LDA + c*ΔE_X_GGA ] + b*E_X_HF + f*E_C_LDA + e*ΔE_C_GGA

a项在所有泛函中都是1。
对于PBE1PBE、BHandHLYP等等，c和d都是等同的，所以你自己设的时候也让它等同。