计算化学公社

标题: 做电场作用下的分子优化和频率分析总是出现错误 [打印本页]

作者
Author: ddddnight 时间: 2021-4-14 18:03
标题: 做电场作用下的分子优化和频率分析总是出现错误
基本基组的设置如下：
%nprocshared=4
%mem=1111MB
%chk=G:\oil gaussian computation\C-O-C(EF).chk
# opt freq m062x/6-311+g(d,p) Nosymm Field=read SCF=QC

out文件最后出错如下：
Mulliken charges with hydrogens summed into heavy atoms:
            1
   1  C 9.013123
   5  C 8.228843
   8  C 7.745332
11  C 20.874124
14  C  -15.954833
17  C  -38.399315
21  O 8.492727
Electronic spatial extent (au):  <R**2>=       14115.1065
Charge=             0.0000 electrons
Dipole moment (field-independent basis, Debye):
X=          -137.4560 Y=             1.0211 Z=       -1114.1707  Tot=          1122.6182
Quadrupole moment (field-independent basis, Debye-Ang):
XX=          -221.5484 YY=          -543.7749 ZZ=       -11600.2543
XY=          -207.8923 XZ=          143.4816 YZ=       -1320.5827
Traceless Quadrupole moment (field-independent basis, Debye-Ang):
XX=          3900.3108 YY=          3578.0843 ZZ=       -7478.3951
XY=          -207.8923 XZ=          143.4816 YZ=       -1320.5827
Octapole moment (field-independent basis, Debye-Ang**2):
  XXX=          206.7945  YYY=       -2386.2105  ZZZ=       -107252.0399  XYY=          -59.1879
  XXY=          -416.1467  XXZ=       -2966.7481  XZZ=          4608.9835  YZZ=       -13638.5834
  YYZ=       -5911.5846  XYZ=          -160.6552
Hexadecapole moment (field-independent basis, Debye-Ang**3):
XXXX=       -3537.5846 YYYY=       -12158.3333 ZZZZ=       -967195.0136 XXXY=          270.1639
XXXZ=          3794.7861 YYYX=          254.9390 YYYZ=       -19868.2802 ZZZX=       56411.7212
ZZZY=       -124881.0623 XXYY=       -2321.2718 XXZZ=       -27158.4777 YYZZ=       -50798.0341
XXYZ=       -4343.7791 YYXZ=          1746.6517 ZZXY=          3036.2550
N-N= 3.567498093180D+06 E-N=-2.459538750037D+07  KE= 1.045898382421D+02
Calling FoFJK, ICntrl=    2127 FMM=F ISym2X=0 I1Cent= 0 IOpClX= 0 NMat=1 NMatS=1 NMatT=0.
dEdEF cannot handle higher than dipole.
Error termination via Lnk1e in F:\G09W\l716.exe at Wed Apr 14 17:49:01 2021.
Job cpu time:    0 days  1 hours  7 minutes 23.0 seconds.
File lengths (MBytes):  RWF=    30 Int=    0 D2E=    0 Chk=    2 Scr=    1

作者
Author: sobereva 时间: 2021-4-14 19:25
8成是输入文件末尾的外电场的设置有问题。没输入文件没法判断
别随便用scf=qc

作者
Author: ddddnight 时间: 2021-4-14 20:33

sobereva 发表于 2021-4-14 19:25
8成是输入文件末尾的外电场的设置有问题。没输入文件没法判断
别随便用scf=qc

老师，您看，这是完整的输入文件
%nprocshared=4
%mem=1111MB
%chk=G:\oil gaussian computation\C-O-C(EF).chk
# opt freq m062x/6-311+g(d,p) Nosymm Field=read SCF=QC

Title Card Required

0 1
C                -1.16837663 0.99474172 0.08973719
H                -0.81046075 -0.01362105 0.08888385
H                -0.81270513 1.50019390 -0.78371337
H                -2.23837572 0.99341609 0.09018955
C                -0.65541026 1.72046490 1.34743007
H                0.41458882 1.72179054 1.34697771
H                -1.01332615 2.72882768 1.34828340
C                -1.16731130 0.99299168 2.60454582
H                -0.80939542 -0.01537110 2.60369249
H                -2.23731038 0.99166605 2.60499818
C                -1.16632206 0.99136664 4.93972527
H                -0.80840618 -0.01699614 4.93887193
H                -2.23632115 0.99004100 4.94017763
C                -0.65335569 1.71708981 6.19741815
H                -1.01127157 2.72545259 6.19827148
H                0.41664339 1.71841545 6.19696579
C                -1.16525673 0.98961659 7.45453391
H                -0.80734085 -0.01874619 7.45368057
H                -0.80884503 1.49385282 8.32838545
H                -2.23525581 0.98829096 7.45498626
O                -0.69098538 1.66687749 3.77240349

1 2 1.0 3 1.0 4 1.0 5 1.0
2
3
4
5 6 1.0 7 1.0 8 1.0
6
7
8 9 1.0 10 1.0 21 1.0
9
10
11 12 1.0 13 1.0 14 1.0 21 1.0
12
13
14 15 1.0 16 1.0 17 1.0
15
16
17 18 1.0 19 1.0 20 1.0
18
19
20
21

0.000 0.000 0.005

作者
Author: sobereva 时间: 2021-4-15 06:44
你得把原子连接关系部分删了啊，你又没写geom=connectivity，程序明显会把这些内容当成读入的电场信息

作者
Author: ddddnight 时间: 2021-4-15 10:01

sobereva 发表于 2021-4-15 06:44
你得把原子连接关系部分删了啊，你又没写geom=connectivity，程序明显会把这些内容当成读入的电场信息

哦哦哦，好的，谢谢老师

作者
Author: ddddnight 时间: 2021-4-16 16:12
本帖最后由 ddddnight 于 2021-4-16 16:15 编辑

sobereva 发表于 2021-4-15 06:44
你得把原子连接关系部分删了啊，你又没写geom=connectivity，程序明显会把这些内容当成读入的电场信息

老师，IRC计算得到的起点结构一定是反应物吗，为什么我IRC计算出来的起点是我想要的产物结构，终点是反应物结构。是不是要根据自己的构型和需要判断反应物和产物，因为毕竟逆反应就对得上了

作者
Author: sobereva 时间: 2021-4-21 15:35

ddddnight 发表于 2021-4-16 16:12
老师，IRC计算得到的起点结构一定是反应物吗，为什么我IRC计算出来的起点是我想要的产物结构，终点是反应 ...

此文明确说了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

此问题和此帖主题无关，和主题不相关的问题应当在别处问

作者
Author: ddddnight 时间: 2021-4-23 11:02

sobereva 发表于 2021-4-21 15:35
此文明确说了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com ...

好的老师，知道了，谢谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)