24小时热门版块排行榜    

查看: 797  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 3/150 2026-11-01 04:04 by YXM1FDCSE8Bn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 04:04 by YXM1FDCSE8Bn
[找工作] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 03:51 by YXM1FDCSE8Bn
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 03:46 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 03:38 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 20/1000 2026-11-01 03:33 by YXM1FDCSE8Bn
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 X5bHxOULxFqb 2026-10-30 25/1250 2026-11-01 03:30 by YXM1FDCSE8Bn
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 X5bHxOULxFqb 2026-10-30 20/1000 2026-11-01 03:30 by YXM1FDCSE8Bn
[考研] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 12/600 2026-11-01 02:19 by YXM1FDCSE8Bn
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-31 4/200 2026-11-01 00:41 by YXM1FDCSE8Bn
[考博] 2027博士申请-临床检验诊断学、病原生物学、微生物方向 50+3 yyh@ 2026-10-28 4/200 2026-10-31 12:47 by 可淡不可忘
[基金申请] 论文的挂名怎么评判贡献?基金申请中,排名第一的挂名文章 +3 biily 2026-10-25 3/150 2026-10-30 18:33 by 旁观2020
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
[考博] 申博求助 20+3 Katrinae 2026-10-26 4/200 2026-10-27 09:56 by 北京莱茵编辑
信息提示
请填处理意见