计算化学公社

标题: 采用g16采用def2-tzvp.1的基组优化甜菜碱分子输入行语法错误 [打印本页]

作者
Author: 陈小黑啊啊啊 时间: 2023-3-16 10:27
标题: 采用g16采用def2-tzvp.1的基组优化甜菜碱分子输入行语法错误
采用def2-tzvp.1 基组几何优化甜菜碱，在g16的basis中没有找到此基组，通过sob老师的博客找到了def2-tzvp.1的gbs文件如下附件所示，其log和报错也如下图所示，想求助输入语法应该怎么输入

作者
Author: yaoyuan0711 时间: 2023-3-16 10:38
把基组文件放到某个目录，在Gaussian输入文件最下面空一行，再 @目录就行了

作者
Author: 陈小黑啊啊啊 时间: 2023-3-16 10:46

yaoyuan0711 发表于 2023-3-16 10:38
把基组文件放到某个目录，在Gaussian输入文件最下面空一行，再 @目录就行了

%Nproc=12
%chk=BHC.chk
%mem=32GB
#p opt def2-tzvp.1 geom=connectivity

247

0 2
O                1.95210610 -0.86637287 -0.05329635
O                1.12417571 -2.46280327 1.42788802
N                -0.68813052 0.14206312 -0.22319701
C                -0.37783178 -0.60308383 1.00536561
C                -1.19887346 -0.56833227 -1.40446739
C                -2.06508143 0.34704990 0.24887506
C                -0.48768109 1.59756746 -0.27049473
C                0.93924682 -1.37850231 0.81661198
H                -0.27493080 0.08142824 1.82130680
H                -1.16998136 -1.29121666 1.21483204
H                -1.52532979 -1.54649309 -1.11892130
H                -2.02175813 -0.02469478 -1.81946241
H                -0.42130296 -0.65090623 -2.13485490
H                -2.07187517 1.09335445 1.01561200
H                -2.67778490 0.66860055 -0.56727175
H                -2.44786073 -0.57159369 0.64190964
H                -0.36940356 1.97381830 0.72416340
H                0.38951484 1.81841195 -0.84202731
H                -1.33725745 2.05992181 -0.72804375
H                0.24859635 0.00263306 -0.54434294

1 8 1.0
您好，这是我的gaussian输入的gjf文件，想请问怎么输入

作者
Author: yaoyuan0711 时间: 2023-3-16 15:56
%Nproc=12
%chk=BHC.chk
%mem=32GB
#p opt def2-tzvp.1 geom=connectivity

247

0 2
O                1.95210610 -0.86637287 -0.05329635
O                1.12417571 -2.46280327 1.42788802
N                -0.68813052 0.14206312 -0.22319701
C                -0.37783178 -0.60308383 1.00536561
C                -1.19887346 -0.56833227 -1.40446739
C                -2.06508143 0.34704990 0.24887506
C                -0.48768109 1.59756746 -0.27049473
C                0.93924682 -1.37850231 0.81661198
H                -0.27493080 0.08142824 1.82130680
H                -1.16998136 -1.29121666 1.21483204
H                -1.52532979 -1.54649309 -1.11892130
H                -2.02175813 -0.02469478 -1.81946241
H                -0.42130296 -0.65090623 -2.13485490
H                -2.07187517 1.09335445 1.01561200
H                -2.67778490 0.66860055 -0.56727175
H                -2.44786073 -0.57159369 0.64190964
H                -0.36940356 1.97381830 0.72416340
H                0.38951484 1.81841195 -0.84202731
H                -1.33725745 2.05992181 -0.72804375
H                0.24859635 0.00263306 -0.54434294

@/home/...... 改成你的路径

另外你计算的泛函呢？末尾的1 8 1.0是什么?

作者
Author: 陈小黑啊啊啊 时间: 2023-3-16 16:31

yaoyuan0711 发表于 2023-3-16 15:56
%Nproc=12
%chk=BHC.chk
%mem=32GB

您好，这个是打开gjf，文件后自带的粉红色的字体

作者
Author: 啊不错的飞过海 时间: 2023-3-16 23:50

yaoyuan0711 发表于 2023-3-16 15:56
%Nproc=12
%chk=BHC.chk
%mem=32GB

1 8 1.0大概是Gaussview画分子时候生成的键合关系，做分子力学计算时候会读，QM方法应该是直接忽略。
Gaussian不写方法时候印象是默认HF？总之应该能跑，不过HF配3-ζ基组做优化好像不是什么常见做法、def2-TZVP.1也不是非常常见的基组，不是很理解楼主拿这套组合做甜菜碱这么个普通有机分子优化的目的

作者
Author: mfdsrax2 时间: 2023-3-17 15:40
#行只能写关键字，def2-tzvp.1不是关键字，所以会报错。正确的写法是
#p opt 方法/gen
你没有指定方法，默认是HF，你确定想用HF计算
在原子坐标的后面空一行写
@ def2-tzvp.1的路径
geom=connectivity对计算没有任何用，可以删掉，后面的连通性列表也一起删，免得出错

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