24小时热门版块排行榜    

查看: 666  |  回复: 2

李路lulu

金虫 (小有名气)

[求助] 在高斯候总如何用IOP来确定交换泛函的比例?我用的是PBE0方法 已有1人参与

希望有详细的输入文件和注解!金币好说!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tawfischer

木虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
李路lulu: 金币+20, ★★★★★最佳答案 2014-03-30 11:37:00
怎么又是你?以下是高斯手册里的内容:
The 1996 pure functional of Perdew, Burke and Ernzerhof [Perdew96a, Perdew97], as made into a hybrid by Adamo [Adamo99a]. The keyword is PBE1PBE. This functional uses 25% exchange and 75% correlation weighting, and is known in the literature as PBE0.
如果你改了,那肯定不能叫PBE0了。
下面是如何自定义了,同样来自高斯手册:
User-Defined Models. Gaussian 09 can use any model of the general form:

    P2EXHF + P1(P4EXSlater + P3ΔExnon-local) + P6EClocal + P5ΔECnon-local
The only available local exchange method is Slater (S), which should be used when only local exchange is desired. Any combinable non-local exchange functional and combinable correlation functional may be used (as listed previously).

The values of the six parameters are specified with various non-standard options to the program:

IOp(3/76=mmmmmnnnnn) sets P1 to mmmmm/10000 and P2 to nnnnn/10000. P1 is usually set to either 1.0 or 0.0, depending on whether an exchange functional is desired or not, and any scaling is accomplished using P3 and P4.

IOp(3/77=mmmmmnnnnn) sets P3 to mmmmm/10000 and P4 to nnnnn/10000.

IOp(3/78=mmmmmnnnnn) sets P5 to mmmmm/10000 and P6 to nnnnn/10000.

For example, IOp(3/76=1000005000) sets P1 to 1.0 and P2 to 0.5. Note that all values must be expressed using five digits, adding any necessary leading zeros.

Here is a route section specifying the functional corresponding to the B3LYP keyword:

#P BLYP IOp(3/76=1000002000) IOp(3/77=0720008000) IOp(3/78=0810010000)
The output file displays the values that are in use:

IExCor=  402 DFT=T Ex=B+HF Corr=LYP ExCW=0 ScaHFX=  0.200000
ScaDFX=  0.800000  0.720000  1.000000  0.810000
where the value of ScaHFX is P2, and the sequence of values given for ScaDFX are P4, P3, P6 and P5.
上面给的是如何自定义B3LYP的方法,PBE0自然要从纯泛函的形式PBEPBE出发自定义。剩下的自己好好领悟,不要再问了。不能让别人手把手一步一步教你,应该学会独立思考。可能我说的不太好听,但希望你能理解。
2楼2014-03-29 13:57:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李路lulu

金虫 (小有名气)

引用回帖:
2楼: Originally posted by tawfischer at 2014-03-29 13:57:23
怎么又是你?以下是高斯手册里的内容:
The 1996 pure functional of Perdew, Burke and Ernzerhof , as made into a hybrid by Adamo . The keyword is PBE1PBE. This functional uses 25% exchange and 75% corr ...

好吧,因为在文献里有看到PBE里HF 的杂化泛函是可以从0%50%调的,谢谢你了
3楼2014-03-30 11:36:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 李路lulu 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311求调剂 +3 lin0039 2026-03-26 3/150 2026-03-26 20:39 by 不吃魚的貓
[考研] 071000生物学求调剂,初试成绩343 +4 小小甜面团 2026-03-25 4/200 2026-03-26 20:26 by wxywyq
[考研] 材料与化工304求B区调剂 +3 邱gl 2026-03-26 6/300 2026-03-26 18:03 by 邱gl
[考研] 334分 一志愿武理 材料求调剂 +4 李李不服输 2026-03-26 4/200 2026-03-26 16:00 by 不吃魚的貓
[考研] 303求调剂 +6 蓝山月 2026-03-25 6/300 2026-03-25 22:47 by 418490947
[考研] 求调剂 +3 李李不服输 2026-03-25 3/150 2026-03-25 13:03 by cmz0325
[考研] 085600材料与化工调剂 +9 A-哆啦Z梦 2026-03-23 15/750 2026-03-25 11:18 by Ainin_
[考研] 340求调剂 +5 话梅糖111 2026-03-24 5/250 2026-03-25 06:53 by ilovexiaobin
[考研] 生物学学硕求调剂 +7 小羊睡着了? 2026-03-23 10/500 2026-03-25 02:24 by 清风拂扬。 m
[考研] 307求调剂 +5 超级伊昂大王 2026-03-24 5/250 2026-03-24 15:46 by 星空星月
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 一志愿华东理工大学081700,初试分数271 +5 kotoko_ik 2026-03-23 6/300 2026-03-24 10:29 by 学术搬砖er
[考研] 一志愿重庆大学085700资源与环境,总分308求调剂 +7 墨墨漠 2026-03-23 8/400 2026-03-23 20:36 by Creta
[考研] 070300,一志愿北航320求调剂 +3 Jerry0216 2026-03-22 5/250 2026-03-23 09:16 by 。。堂堂
[考研] 336求调剂 +5 rmc8866 2026-03-21 5/250 2026-03-21 17:24 by 学员8dgXkO
[考研] 一志愿重庆大学085700资源与环境总分308求调剂 +7 墨墨漠 2026-03-20 7/350 2026-03-21 16:36 by barlinike
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 求调剂 +3 @taotao 2026-03-20 3/150 2026-03-20 19:35 by JourneyLucky
信息提示
请填处理意见