计算化学公社

标题: g16对结构优化后的含过渡金属碳面吸附苯分子结构做频率分析时出现l1002报错 [打印本页]

作者
Author: 周瑾根 时间: 2023-12-7 13:57
标题: g16对结构优化后的含过渡金属碳面吸附苯分子结构做频率分析时出现l1002报错
Linux虚拟机下Gaussian16A03.对结构优化后的14环含Fe-N-C结构的碳平面吸附苯分子做频率分析时出现如下报错：
No inversion for dimension 144.
N<=N0 in RedDim.
Error termination via Lnk1e in /home/zjg/g16/l1002.exe at Thu Dec 7 06:51:37 2023.

gjf文件开头关键词为:%nprocshared=6
%mem=25000MB
%chk=14FeN4benfreq.chk
#p freq genecp em=gd3bj
iop(3/174=1000000,3/175=2238200,3/177=452900,3/178=4655000) tpssh

赝势基组分别为Fe：LANL2TZ(f)、C H N：Def2-SVP

查阅了常见报错贴没找到解决办法，恳请各位老师指导，感谢各位前辈🙏

作者
Author: sobereva 时间: 2023-12-7 14:22
尝试Gaussian更新版本，或者换成SDD结合6-31G*试试

作者
Author: zjxitcc 时间: 2023-12-7 14:23
我有个疑问，为啥你不读取opt最后一帧的轨道？这种开壳层过渡金属体系容易有多个SCF解、容易收敛到不稳定波函数上，你现在从零开始算，除了浪费机时外还对SCF十分不利。

作者
Author: 周瑾根 时间: 2023-12-7 14:35

zjxitcc 发表于 2023-12-7 14:23
我有个疑问，为啥你不读取opt最后一帧的轨道？这种开壳层过渡金属体系容易有多个SCF解、容易收敛到不稳定波 ...

抱歉老师我不是很明白您的意思，我是将opt和freq分开计算，这里的freq任务我是选择opt任务后的最后一帧构型来进行计算的，相同的泛函和赝势基组下opt是正常结束了。

作者
Author: 周瑾根 时间: 2023-12-7 14:41

sobereva 发表于 2023-12-7 14:22
尝试Gaussian更新版本，或者换成SDD结合6-31G*试试

好的我尝试一下，感谢sob老师！

作者
Author: zjxitcc 时间: 2023-12-7 14:47

周瑾根发表于 2023-12-7 14:35
抱歉老师我不是很明白您的意思，我是将opt和freq分开计算，这里的freq任务我是选择opt任务后的最后一帧构 ...

你只是取了opt最后一帧结构，并没有读取opt最后一帧的轨道，如果读进来轨道，SCF会1圈收敛。但你扔掉了轨道，选择了重新开始SCF迭代。

作者
Author: 周瑾根 时间: 2023-12-7 15:58

zjxitcc 发表于 2023-12-7 14:47
你只是取了opt最后一帧结构，并没有读取opt最后一帧的轨道，如果读进来轨道，SCF会1圈收敛。但你扔掉了轨 ...

原来如此，我确实不知道还有这种差异，感谢老师提醒，请问老师我应该添加什么命令来让高斯读取最后一帧的轨道呢？

作者
Author: sobereva 时间: 2023-12-8 06:00

周瑾根发表于 2023-12-7 15:58
原来如此，我确实不知道还有这种差异，感谢老师提醒，请问老师我应该添加什么命令来让高斯读取最后一帧的 ...

guess=read，%chk指定成优化任务的chk文件

作者
Author: 周瑾根 时间: 2023-12-8 09:36

sobereva 发表于 2023-12-8 06:00
guess=read，%chk指定成优化任务的chk文件

感谢sob老师的回复，昨天使用SDD结合6-31G*的计算出现了新的l1002错误：
OrtVc1:  Ph=3 IOff=    0 NNewB=    0 IAdd=    5 IPass=20 DotMx2= 1.69+280
OrtVc1 failed #3.
Error termination via Lnk1e in /home/zjg/g16/l1002.exe at Thu Dec  7 21:51:06 2023.

以前也遇到过，不过重新续算之后就没有出现了，这次我加上oldchk和guess=read，g16正在计算中，新的输入文件开头如下（赝势基组仍为SDD、6-31G*）：
%nprocshared=6
%mem=25000MB
%oldchk=/home/zjg/Gauss/chek/14FeN4benopt.chk
%chk=14FeN4benfreq.chk
#p freq genecp em=gd3bj guess=read  iop(3/174=1000000,3/175=2238200,3/177=452900,3/178=4655000) tpssh

请问sob老师这是否正确书写了？以及这个OrtVc1 failed是什么问题导致的错误？

作者
Author: sobereva 时间: 2023-12-8 10:49

周瑾根发表于 2023-12-8 09:36
感谢sob老师的回复，昨天使用SDD结合6-31G*的计算出现了新的l1002错误：
OrtVc1: Ph=3 IOff= 0 NNe ...

bug

作者
Author: 周瑾根 时间: 2023-12-8 11:02

sobereva 发表于 2023-12-8 10:49
bug

了然，感谢老师的回复

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