24小时热门版块排行榜

返回列表

【奖励】本帖被评价16次，作者yjcmwgk增加金币 14.7 个

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

yjcmwgk

禁虫 (文坛精英)

[资源] 【分享】一些特殊的DFT在Gaussian03中如何输入？（不断更新中）

引用回帖:
9楼的问题：特殊的？是否指G03中并不存在的？或者是隐藏的？或者是需要用户自己添加的？若G03中找不到，该到何处寻觅？

回答：我所谓的特殊的指的是G03中不存在的DFT，但是可以通过用户自定义来得到的DFT。
Gaussian03可以使用具有下列一般形式的任何模型
P1*[P3*ΔEx(non-local)+P4*Ex(Slater)]+P2*Ex(HF)+P5*ΔEc(non-local)+P6*Ec(local)
P1~P6的值由IOP(3)来指定
一套指定方式，支持四位有效数字
IOp(3/45=mmmmnnnn)
指定P1为mmmm/1000，P2为nnnn/1000.
IOp(3/46=mmmmnnnn)
指定P3为mmmm/1000，P4为nnnn/1000。
IOp(3/47=mmmmnnnn)
指定P5为mmmm/1000，P6为nnnn/1000。
另一套指定方式，支持五位有效数字
IOp(3/76=MMMMMNNNNN)
指定P1为MMMMM/10000，P2为NNNNN/10000.
IOp(3/77=MMMMMNNNNN)
指定P3为MMMMM/10000，P4为NNNNN/10000.
IOp(3/78=MMMMMNNNNN)
指定P5为MMMMM/10000，P6为NNNNN/10000.
其实依照这个例子，你自己发明一个新泛函也行（你的新发明是否合理，那就有待商榷了）
举个例子
# BLYP IOp(3/45=10000200) IOp(3/46=07200800) IOp(3/47=08101000)
# BLYP IOp(3/76=1000002000) IOp(3/77=0720008000) IOp(3/78=0810010000)
这两个输入路径其实都相当于B3LYP

(1)mpw1k/6-31+G(d,p)
# mpwpw91/6-31+G(d,p) iop(3/76=0572004280)
(2)MPWB1K/6-31+G(d,p)
# mpwb95/6-31+G(d,p) IOp(3/76=0560004400)
(3)MPW1B95/6-31+G(d,p)
# mpwb95/6-31+G(d,p) IOp(3/76=0690003100)
(4)TPSS1KCIS/6-31+G(d,p)
# tpsskcis/6-31+G(d,p) IOp(3/76=0870001300)
(5)MPW3LYP/6-31+G(d,p)
# mpwlyp/6-31+G(d,p) IOp(3/76=1000002000) IOp(3/77=0720008000) IOp(3/78=0810010000)
(6)MPWKCIS1K/6-31+G(d,p)
# MPWKCIS/6-31+G(d,p) IOp(3/76= 0590004100)
(7)MPW1KCIS/6-31+G(d,p)
# MPWKCIS/6-31+G(d,p) IOp(3/76= 0850001500)
(8)PBE1KCIS/6-31+G(d,p)
# pbekcis/6-31+G(d,p) IOp(3/76=0780002200)
(9)PBE1W/6-31+G(d,p)
# pbepbe/6-31+G(d,p) IOp(3/78=0740010000)
(10)MPWLYP1W/6-31+G(d,p)
# mpwv5lyp/6-31+G(d,p) IOp(3/78=0880010000)
(11)PBELYP1W
# pbev5lyp/6-31+G(d,p) IOp(3/78=0540010000)
(12)TPSSLYP1W
# tpssv5lyp/6-31+G(d,p) IOp(3/78=0740010000)
(13)MPWLYP1M/6-31+G(d,p)
# mpwlyp/6-31+G(d,p) IOp(3/76=0950000500)
(14)MOHLYP/6-31+G(d,p)
# OV5LYP/6-31+G(d,p) IOp(3/77=1292010000) IOp(3/78=0500010000)
(15)MOHLYP2/6-31+G(d,p)
# OV5LYP/6-31+G(d,p) IOp(3/77=1849810515) IOp(3/78=0500005000) (for Gaussian03 B&C)
# OV5LYP/6-31+G(d,p) IOp(3/77=1292010000) IOp(3/78=0500005000) (for Gaussian03 D&E)

[ Last edited by yjcmwgk on 2010-1-14 at 09:50 ]

回复此楼