计算化学公社

 找回密码 Forget password
 注册 Register
Views: 644|回复 Reply: 2
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] Gaussian输入文件中的自定义基组与输出文件不完全一致

[复制链接 Copy URL]

19

帖子

0

威望

241

eV
积分
260

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
使用Gaussian16自定义基组,输出文件中基组的系数不同,请问这是我输入的格式问题吗?

输入文件

%chk=uGAUorb.chk
%mem=90GB
%nprocshared=36
# 6d 10f genecp nosymm upbe1pbe scf=(maxcycle=120,xqc) scfcyc=999 pop=full IOp(3/24=1)

Title Card Required

0 1
Fe                     0.04379120   -0.04849457    0.03454395
O                      1.64539120   -0.04849457    0.03454395

Fe 0
S    3   1.00
      6.4220000             -0.3927882
      1.8260000              0.7712643
      0.7135000              0.4920228
S    4   1.00
      6.4220000              0.1786877
      1.8260000             -0.4194032
      0.7135000             -0.4568185
      0.1021000              1.1035048
S    1   1.00
      0.0363000              1.0000000
P    3   1.00
     19.4800000             -0.0470282
      2.3890000              0.6248841
      0.7795000              0.4722542
P    1   1.00
      0.0740000              1.0000000
P    1   1.00
      0.0220000              1.0000000
D    4   1.00
     37.0800000              0.0329000
     10.1000000              0.1787418
      3.2200000              0.4487657
      0.9628000              0.5876361
D    1   1.00
      0.2262000              1.0000000
****
O 0
6-31G*
****

Fe 0
FE-ECP     2     10
d potential
  3
1    392.6149787            -10.0000000
2     71.1756979            -63.2667518
2     17.7320281            -10.9613338
s-d potential
  5
0    126.0571895              3.0000000
1    138.1264251             18.1729137
2     54.2098858            339.1231164
2      9.2837966            317.1068012
2      8.6289082           -207.3421649
p-d potential
  5
0     83.1759490              5.0000000
1    106.0559938              5.9535930
2     42.8284937            294.2665527
2      8.7701805            154.4244635
2      8.0397818            -95.3164249

输出文件截取

*******************************************************************************************************************************
          ATOMIC CENTER             *        ATOMIC ORBITAL       *                     GAUSSIAN FUNCTIONS                    *
*******************************************************************************************************************************
                                    * FUNCTION    SHELL    SCALE  *                                                           *
ATOM   X-COORD  Y-COORD  Z-COORD   *  NUMBER     TYPE     FACTOR *  EXPONENT    S-COEF      P-COEF      D-COEF      F-COEF   *
*******************************************************************************************************************************
  Fe     0.08275 -0.09164  0.06528                                                                                            *
                                            1   S           1.00                                                               
                                                                  0.182600D+01-.774426D-010.000000D+000.000000D+000.000000D+00
                                                                  0.713500D+00-.130105D+000.000000D+000.000000D+000.000000D+00*
                                                                  0.102100D+000.143371D+000.000000D+000.000000D+000.000000D+00*
                                            2   S           1.00                                                              *
                                                                  0.642200D+01-.112933D+010.000000D+000.000000D+000.000000D+00
                                                                  0.182600D+010.863453D+000.000000D+000.000000D+000.000000D+00*
                                                                  0.713500D+000.272233D+000.000000D+000.000000D+000.000000D+00*
                                            3   S           1.00                                                              *
                                                                  0.363000D-010.592707D-010.000000D+000.000000D+000.000000D+00
                                        4-  6   P           1.00                                                              *
                                                                  0.194800D+020.000000D+00-.274337D+010.000000D+000.000000D+00
                                                                  0.238900D+010.000000D+000.264551D+010.000000D+000.000000D+00*
                                                                  0.779500D+000.000000D+000.493044D+000.000000D+000.000000D+00*
                                        7-  9   P           1.00                                                              *
                                                                  0.740000D-010.000000D+000.550148D-010.000000D+000.000000D+00
                                       10- 12   P           1.00                                                              *
                                                                  0.220000D-010.000000D+000.120773D-010.000000D+000.000000D+00
                                       13- 18   D           1.00                                                              *
                                                                  0.370800D+020.000000D+000.000000D+000.301717D+020.000000D+00
                                                                  0.101000D+020.000000D+000.000000D+000.168344D+020.000000D+00*
                                                                  0.322000D+010.000000D+000.000000D+000.571711D+010.000000D+00*
                                                                  0.962800D+000.000000D+000.000000D+000.905120D+000.000000D+00*
                                       19- 24   D           1.00                                                              *
                                                                  0.226200D+000.000000D+000.000000D+000.122116D+000.000000D+00


3813

帖子

4

威望

8002

eV
积分
11895

Level 6 (一方通行)

MOKIT开发者

2#
发表于 Post on 2023-9-19 11:15:38 | 只看该作者 Only view this author
估计是你没写int=nobasistransform的原因。
顺带一提,传轨道前建议在gjf文件里写上nosymm int=nobasistransform,确保传轨道效果最好。
自动做多参考态计算的程序MOKIT

19

帖子

0

威望

241

eV
积分
260

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-9-19 12:00:14 | 只看该作者 Only view this author
zjxitcc 发表于 2023-9-19 11:15
估计是你没写int=nobasistransform的原因。
顺带一提,传轨道前建议在gjf文件里写上nosymm int=nobasistra ...

谢谢!就是这个问题,加上之后就对了!非常感谢!

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-26 16:52 , Processed in 0.563768 second(s), 21 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list