24小时热门版块排行榜    

查看: 729  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上项目filecode邪修 +3 西山十月 2026-08-09 3/150 2026-08-09 19:39 by Tide man
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +7 Lanmanbaby 2026-08-09 13/650 2026-08-09 19:04 by 天神眷顾
[基金申请] 关于代码变化问题,想知道的进来 +14 且听虎啸 2026-08-07 20/1000 2026-08-09 18:03 by 天神眷顾
[基金申请] 这样的filecode谁见过 +10 布布和一二 2026-08-08 21/1050 2026-08-09 17:21 by 天神眷顾
[基金申请] filecode与中标关系的预测 +5 布布和一二 2026-08-07 5/250 2026-08-09 16:15 by 袁向阳007
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 16:47 by oEVWOejN9taj
[基金申请] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 16:27 by oEVWOejN9taj
[论文投稿] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 15:07 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 KXLV3nuBVBY7 2026-08-07 6/300 2026-08-08 14:01 by oEVWOejN9taj
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见