应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 混合基组的使用
101/1返回列表
查看: 2336  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

heqian9012

新虫 (初入文坛)

[求助] 混合基组的使用

我需要用高斯计算PH3AgCl相互作用的结合能，其中P、H、Cl要使用aug-cc-pvdz基组，Ag要使用aug-cc-pvdz-pp基组，下面是我综合论坛上面的帖子编辑的输入文件，但还是不对，求各位大虾指点指点
%mem=128MW
%nproc=1
# wb97x/gen opt counterpoise=2

Title Card Required

0 1
P                 -1.55000000   -1.59165400    1.54258700 1
H                 -0.32071500   -2.06702800    2.01796100 1
H                 -2.16464200   -1.31424800    2.77074100 1
H                 -2.16464200   -2.81980800    1.26518100 1
Ag                -1.55000000    0.00000000   -0.04906700 1
Cl                -1.55000000    1.58444700   -1.63351400 1
P                  1.55000000    1.59165400    1.54258700 2
H                  0.32071500    2.06702800    2.01796100 2
H                  2.16464200    2.81980800    1.26518100 2
H                  2.16464200    1.31424800    2.77074100 2
Ag                 1.55000000    0.00000000   -0.04906700 2
Cl                 1.55000000   -1.58444700   -1.63351400 2

P H Cl 0 gen
aug-cc-pvdz
****
Ag 0
aug-cc-pvdz-pp
****

Ag 0
aug-cc-pvdz-pp
回复此楼

» 收录本帖的淘帖专辑推荐

量化计算有意义的讨论

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-05-12 21:52:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

思雨G十年

金虫 (正式写手)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
gmy1990: 金币+3 2012-05-13 09:38:47
头文件
# 关键词
空行
标题部分
空行
电荷 自旋多重度
坐标
空一行
元素 0
对应机组
****
元素 0
对应基组
****
空多行
赞一下 回复此楼
要以感性的思想和认识，理性的思考和行动对待科研。
2楼2012-05-12 22:39:31
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gujing136

禁虫 (文坛精英)

【答案】应助回帖

★ ★
感谢参与，应助指数 +1
gmy1990: 金币+2 2012-05-13 19:13:12
heqian,你这里有一个问题，你没有发现么？就是你算的BSSE的，你的多重度还没有完全指定，另外，你在这里问的话，最好把你的out文件的最后20行或者10行贴出来。
赞一下 回复此楼
叫老猫的兔子
3楼2012-05-13 10:53:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heqian9012

新虫 (初入文坛)
引用回帖:
3楼: Originally posted by gujing136 at 2012-05-13 10:53:08:
heqian,你这里有一个问题，你没有发现么？就是你算的BSSE的，你的多重度还没有完全指定，另外，你在这里问的话，最好把你的out文件的最后20行或者10行贴出来。

好麻烦额。。。。。。。。
回复此楼
4楼2012-05-13 19:53:54
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heqian9012

新虫 (初入文坛)
引用回帖:
2楼: Originally posted by 思雨G十年 at 2012-05-12 22:39:31:
头文件
# 关键词
空行
标题部分
空行
电荷 自旋多重度
坐标
空一行
元素 0
对应机组
****
元素 0
对应基组
****
空多行

虽然还是不太懂，还是谢谢。
赞一下 回复此楼
5楼2012-05-13 19:54:38
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gujing136

禁虫 (文坛精英)

【答案】应助回帖

已经解决

因为银离子的基函数的问题，所以需要将银离子的基函数定义才好……

输入文件如下，供大家参考。


#p wb97x/Gen counterpoise=2

Title Card Required

0 1 0 1 0 1
P                 -1.55000000   -1.59165400    1.54258700 1
H                 -0.32071500   -2.06702800    2.01796100 1
H                 -2.16464200   -1.31424800    2.77074100 1
H                 -2.16464200   -2.81980800    1.26518100 1
Ag                -1.55000000    0.00000000   -0.04906700 1
Cl                -1.55000000    1.58444700   -1.63351400 1
P                  1.55000000    1.59165400    1.54258700 2
H                  0.32071500    2.06702800    2.01796100 2
H                  2.16464200    2.81980800    1.26518100 2
H                  2.16464200    1.31424800    2.77074100 2
Ag                 1.55000000    0.00000000   -0.04906700 2
Cl                 1.55000000   -1.58444700   -1.63351400 2

P H Cl 0
aug-cc-pvdz
****
Ag     0
S   7   1.00
    180.0750000              0.0008490        
     21.8987000             -0.0654500        
     13.8670000              0.2977650        
      6.1426300             -0.7531210        
      1.4381400              0.8811750        
      0.6483820              0.4351760        
      0.1288820              0.0147380        
S   7   1.00
    180.0750000             -0.0002030        
     21.8987000              0.0157230        
     13.8670000             -0.0792290        
      6.1426300              0.2226340        
      1.4381400             -0.3491960        
      0.6483820             -0.2559780        
      0.1288820              0.5486660        
S   7   1.00
    180.0750000             -0.0008620        
     21.8987000              0.0525460        
     13.8670000             -0.2081000        
      6.1426300              0.5249720        
      1.4381400             -1.2701730        
      0.6483820              0.5393730        
      0.1288820              1.6533910        
S   1   1.00
      0.0457380              1.0000000        
S   1   1.00
      0.0162000              1.0000000        
P   6   1.00
     11.8751000              0.1162480        
      8.0024500             -0.3072860        
      2.0176600              0.5157360        
      0.9542300              0.5031040        
      0.4231180              0.1420950        
      0.1358850              0.0051530        
P   6   1.00
     11.8751000             -0.0282840        
      8.0024500              0.0783470        
      2.0176600             -0.1567410        
      0.9542300             -0.1886250        
      0.4231180              0.0635650        
      0.1358850              0.5817780        
P   6   1.00
     11.8751000             -0.0528570        
      8.0024500              0.1470390        
      2.0176600             -0.3129730        
      0.9542300             -0.3708280        
      0.4231180              0.4236290        
      0.1358850              0.7536550        
P   1   1.00
      0.0454010              1.0000000        
P   1   1.00
      0.0152000              1.0000000        
D   5   1.00
     26.4320000              0.0034790        
     11.0345000             -0.0138480        
      2.7378700              0.2545990        
      1.1957500              0.4498490        
      0.4820420              0.3757380        
D   5   1.00
     26.4320000             -0.0047330        
     11.0345000              0.0190700        
      2.7378700             -0.4332980        
      1.1957500             -0.4445680        
      0.4820420              0.4914420        
D   1   1.00
      0.1729080              1.0000000        
D   1   1.00
      0.0620000              1.0000000        
F   1   1.00
      1.3299000              1.0000000        
F   1   1.00
      0.4449000              1.0000000
****
赞一下 回复此楼
叫老猫的兔子
6楼2012-05-13 20:00:37
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

heqian9012

新虫 (初入文坛)
引用回帖:
6楼: Originally posted by gujing136 at 2012-05-13 20:00:37:
已经解决

因为银离子的基函数的问题，所以需要将银离子的基函数定义才好……

输入文件如下，供大家参考。


#p wb97x/Gen counterpoise=2

Title Card Required

0 1 0 1 0 1
P                  ...

输入信息就改了方法，方法是b97d
循环到120多。。。

错误信息
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RB97D) =  -4046.92640194     A.U. after  129 cycles
             Convg  =    0.2241D-01             -V/T =  2.9240
KE= 2.103408071264D+03 PE=-1.193043212226D+04 EE= 3.946645829858D+03
Convergence failure -- run terminated.
Error termination via Lnk1e in d:\G09W\l502.exe at Mon May 14 06:36:04 2012.
Job cpu time:  0 days 10 hours  4 minutes 43.0 seconds.
File lengths (MBytes):  RWF=     43 Int=      0 D2E=      0 Chk=      2 Scr=      1
赞一下 回复此楼
7楼2012-05-14 06:41:08
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gujing136

禁虫 (文坛精英)

【答案】应助回帖

引用回帖:
7楼: Originally posted by heqian9012 at 2012-05-14 06:41:08:
输入信息就改了方法，方法是b97d
循环到120多。。。

错误信息
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(RB97D) =  -4046.92640194     A.U. after  129 c ...

收敛失败，一般输入结构的问题……
赞一下 回复此楼
叫老猫的兔子
8楼2012-05-14 07:30:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ahui303

新虫 (初入文坛)
引用回帖:
8楼: Originally posted by gujing136 at 2012-05-14 07:30:40:
收敛失败，一般输入结构的问题……

弱弱的问下，那该怎么去调整结构！！
赞一下 回复此楼
9楼2012-05-14 09:10:26
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gujing136

禁虫 (文坛精英)

【答案】应助回帖

引用回帖:
9楼: Originally posted by ahui303 at 2012-05-14 09:10:26:
弱弱的问下，那该怎么去调整结构！！

收敛失败有很多解决的办法，请谷歌+百度，谢谢···
赞一下 回复此楼
叫老猫的兔子
10楼2012-05-14 10:08:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 heqian9012 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭