计算化学公社

标题: 含有N、H、O、C、I和Cu的一个大分子，赝势混合基组输入文件错误，求大神指点 [打印本页]

作者
Author: Zhanga 时间: 2023-5-23 21:50
标题: 含有N、H、O、C、I和Cu的一个大分子，赝势混合基组输入文件错误，求大神指点
各位老师好，我想要优化一个含有N、H、O、C、I和Cu的一个大分子，想用赝势混合基组，因为计算条件较差，故Cu选用Lanl2DZ,I选用Lanl08(d)，但是在服务器上gif文件提交不上，log文件中显示输入文件错误： A syntax error was detected in input Line . Lanl08(d)请各位老师指点，到底是哪里错误？如何解决？
%nprocshared=32
%mem=4GB
%chk=1-1.chk
#p opt freq b3lyp/genecp em=gd3bj

Title Card Required

0 1
N                -1.49040000 -2.37080000 1.10730000
H                -1.83410000 -1.58510000 1.04690000
C                -2.23780000 -3.33710000 1.57350000
O                -1.84080000 -4.47830000 1.73300000
C                -3.65300000 -2.98350000 1.94680000
C                -4.23840000 -3.60710000 3.01070000
C                -4.39620000 -2.05520000 1.27150000
C                -5.53670000 -3.29240000 3.35150000
H                -3.76480000 -4.23870000 3.50170000
H                -4.02960000 -1.59460000 0.55130000
C                -5.70160000 -1.80610000 1.67110000
N                -6.27000000 -2.40550000 2.70410000
H                -5.91760000 -3.73020000 4.07780000
H                -6.20160000 -1.18730000 1.18940000
Cu             -8.18130000 -2.03190000 3.25950000
I                -8.62090000 -2.41190000 5.75030000
I                -9.21070000 -4.78840000 1.99340000
Cu             -9.65030000 -4.40840000 4.48420000
C                -0.09428684 -2.55192882 0.68426290
H                0.29408345 -1.62075865 0.32790545
H                0.48973422 -2.88712410 1.51580576
H                -0.05046020 -3.27974577 -0.09884739

Cu 0
Lanl2DZ
I 0
Lanl08(d)
****
C O N H 0
6-31G*
****

Cu 0
Lanl2DZ
I 0
Lanl08(d)

作者
Author: sobereva 时间: 2023-5-23 23:50
lanl08系列不是内置的，去BSE拷定义按下文使用
详解Gaussian中混合基组、自定义基组和赝势基组的输入
http://sobereva.com/60

甭给Cu用lanl2DZ那么烂的基组，体系又没多大。改成SDD
看
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）

作者
Author: Zhanga 时间: 2023-5-24 11:08
本帖最后由 Zhanga 于 2023-5-24 16:57 编辑

老师，您好，感谢您的解答。
因为后续工作还有更大的分子，所以我还是想先给Cu用Lanl2DZ基组跑一下，看一下结果如何。后期会考虑SDD。
根据您的解答，我对输入文件进行了修改，服务器上提交成功，但是11分钟之后报错退出。log文件中显示未满足收敛标准。问题一：可能是输入文件依然有问题么？输入文件如下，可否请老师帮我检查一下是否正确呢？
问题二：如果不是输入文件错误，可能是因为Lanl08(d)的基组太大了么？
问题三：SDD基组可否用于碘？
请老师指教，超级感谢。
%nprocshared=32
%mem=4GB
%chk=1-4.chk
#p opt freq b3lyp/genecp em=gd3bj

Title Card Required

0 1
N                -1.49040000 -2.37080000 1.10730000
H                -1.83410000 -1.58510000 1.04690000
C                -2.23780000 -3.33710000 1.57350000
O                -1.84080000 -4.47830000 1.73300000
C                -3.65300000 -2.98350000 1.94680000
C                -4.23840000 -3.60710000 3.01070000
C                -4.39620000 -2.05520000 1.27150000
C                -5.53670000 -3.29240000 3.35150000
H                -3.76480000 -4.23870000 3.50170000
H                -4.02960000 -1.59460000 0.55130000
C                -5.70160000 -1.80610000 1.67110000
N                -6.27000000 -2.40550000 2.70410000
H                -5.91760000 -3.73020000 4.07780000
H                -6.20160000 -1.18730000 1.18940000
Cu             -8.18130000 -2.03190000 3.25950000
I                -8.62090000 -2.41190000 5.75030000
I                -9.21070000 -4.78840000 1.99340000
Cu             -9.65030000 -4.40840000 4.48420000
C                -0.09428684 -2.55192882 0.68426290
H                0.29408345 -1.62075865 0.32790545
H                0.48973422 -2.88712410 1.51580576
H                -0.05046020 -3.27974577 -0.09884739

I    0
S 1 1.00
   0.7242000             1.0000000
S 1 1.00
   0.4653000             1.0000000
S 1 1.00
   0.1336000             1.0000000
P 1 1.00
   1.2900000             1.0000000
P 1 1.00
   0.3180000             1.0000000
P 1 1.00
   0.1053000             1.0000000
P 1 1.00
   0.0308                1.0
D 1 1.00
   0.2940                1.0
****
Cu 0
Lanl2DZ
****
C O N H 0
6-31G*
****

Cu 0
Lanl2DZ

作者
Author: wzkchem5 时间: 2023-5-24 17:46

Zhanga 发表于 2023-5-24 04:08
老师，您好，感谢您的解答。
因为后续工作还有更大的分子，所以我还是想先给Cu用Lanl2DZ基组跑一下，看一 ...

结构不合理，除非你研究的是气相团簇，否则Cu的第一配位层必须完全饱和。
此外这是典型的反铁磁耦合体系，必须用片段初猜或者stable=opt的方法得到开壳层单重态，不能一上来就直接算。

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