24小时热门版块排行榜    

查看: 3050  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

何飘同学

新虫 (小有名气)

[求助] 高斯采用genecp 混合基组做opt计算报错!已有1人参与

我最近在用Gaussian软件做Ru金属配合物的opt,采用genecp 混合基组(其中过渡金属Ru 原子采用Aug-cc-pVDZ-PP 相关一致基组,其他的轻原子采用6-311++G(d,p)全电子基组),但是计算报错,提示信息为语法错误,还请各位大牛指点,并给出输入文件的写法,十分感谢!
1.输入文件:
%chk=Ru-c-anion-1-pvdz.chk
%mem=6GB
%nprocshared=8
#p b3lyp/genecp opt freq nosymm scf(maxcycle=1000)

opt Ru-c-anion

-2 1
几何坐标

C H O N P Cl 0
6-311++G(d,p)
****
Ru 0
aug-cc-pvdz-pp
****

Ru 0
aug-cc-pvdz-pp
空行
END

2.输出文件报错部分:
Centers:       1     24
aug-cc-pvdz-pp
QPErr --- A syntax error was detected in the input line.
aug-cc-pvdz-pp
            '
Last state= "Last"
TCursr=      5805 LCursr=        11
Error reading general basis specification.
Error termination via Lnk1e in /public/software/g09D/l301.exe at Thu Dec 27 18:28:56 2018.
Job cpu time:       0 days  0 hours  0 minutes  4.7 seconds.
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

何飘同学

新虫 (小有名气)

引用回帖:
3楼: Originally posted by sophiazbs at 2018-12-27 23:26:45
基组没有读到,把ru换成原子编号试试

你好!感谢您的回复!方便请你写一个输入文件吗?万分感激!
4楼2018-12-28 10:24:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

sophiazbs

银虫 (正式写手)

不如行动
2楼2018-12-27 23:18:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sophiazbs

银虫 (正式写手)

基组没有读到,把ru换成原子编号试试

发自小木虫IOS客户端
不如行动
3楼2018-12-27 23:26:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
zhou2009: 金币+10 2018-12-29 10:18:58
高斯对于Ru 没内置aug-cc-pvdz-pp基组  手动导入就行了  

PS 这是基组库 要找别的基组可以在这里找到。
https://bse.pnl.gov/bse/portal#opennewwindow

把基组那里改成下面这样就行。

C H O N P Cl 0
6-311++G(d,p)
****
Ru     0
S   7   1.00
    159.3720000              0.0006730        
     20.9016000             -0.0371610        
     13.0741000              0.1665170        
      4.7147400             -0.6159540        
      1.1361200              0.8428300        
      0.5154640              0.4477630        
      0.1067670              0.0184450        
S   7   1.00
    159.3720000             -0.0001690        
     20.9016000              0.0093920        
     13.0741000             -0.0477130        
      4.7147400              0.2011130        
      1.1361200             -0.3734210        
      0.5154640             -0.2738590        
      0.1067670              0.6093840        
S   7   1.00
    159.3720000              0.0008950        
     20.9016000             -0.0396350        
     13.0741000              0.1442430        
      4.7147400             -0.5192870        
      1.1361200              1.5854040        
      0.5154640             -0.9082030        
      0.1067670             -1.4720020        
S   1   1.00
      0.0392180              1.0000000        
S   1   1.00
      0.0144000              1.0000000        
P   6   1.00
     12.8738000              0.0299680        
      6.1475000             -0.1879300        
      1.5652400              0.4999900        
      0.7361940              0.4957320        
      0.3263170              0.1475150        
      0.1124070              0.0060880        
P   6   1.00
     12.8738000             -0.0078010        
      6.1475000              0.0536560        
      1.5652400             -0.1719400        
      0.7361940             -0.1991700        
      0.3263170              0.0744110        
      0.1124070              0.5971370        
P   6   1.00
     12.8738000             -0.0135020        
      6.1475000              0.0937330        
      1.5652400             -0.3159100        
      0.7361940             -0.3618000        
      0.3263170              0.4008770        
      0.1124070              0.7537660        
P   1   1.00
      0.0391720              1.0000000        
P   1   1.00
      0.0137000              1.0000000        
D   5   1.00
      6.2126500             -0.0245970        
      2.7853900              0.1273340        
      1.2685100              0.3570810        
      0.5517900              0.4082070        
      0.2247310              0.2812400        
D   5   1.00
      6.2126500              0.0261160        
      2.7853900             -0.1459400        
      1.2685100             -0.4425210        
      0.5517900             -0.1690200        
      0.2247310              0.5052370        
D   1   1.00
      0.0824110              1.0000000        
D   1   1.00
      0.0302000              1.0000000        
F   1   1.00
      0.8500000              1.0000000        
F   1   1.00
      0.2719000              1.0000000        
****

RU     0
RU-ECP     4     28
g-ul potential
  1
2      1.0000000              0.0000000        
s-ul potential
  2
2     11.50059000           209.78649300      
2      5.06857500            30.21430700      
p-ul potential
  4
2     10.53263400            48.75124400      
2     10.19201000            97.49652900      
2      4.73489200             7.86018800      
2      4.50906500            15.32975100      
d-ul potential
  4
2      8.87797700            26.96750600      
2      8.76612200            40.43230300      
2      3.17019600             3.34075800      
2      3.22885100             5.25635200      
f-ul potential
  2
2      7.82024900            -8.84752500      
2      7.83964700           -11.83551800

» 本帖已获得的红花(最新10朵)

5楼2018-12-28 10:48:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见