计算化学公社

标题: 求助含23个Cu原子的团簇的计算过程中出现l508报错 [打印本页]

作者
Author: 化学girl 时间: 2026-1-13 09:31
标题: 求助含23个Cu原子的团簇的计算过程中出现l508报错
老师们好，在过渡金属团簇的计算过程中，出现了l508的报错，这是有23个Cu原子的团簇
之前还算了有12个Cu原子的团簇就收敛了，用相同的方法算这个大一点的就不收敛，还有一个25个Cu的，也不收敛

本人属于理论计算的新手小白，Cu23的输入文件和输出文件上传至附件，请各位老师批评指教~

作者
Author: wal 时间: 2026-1-13 09:53
是xqc收敛失败，本质上还是SCF不收敛，参考SCF不收敛解决方案http://sobereva.com/61
题外话，为什么直接做freq？

作者
Author: paramecium86 时间: 2026-1-13 10:35
你如果是要用lanl2dz 那么对于Cu是有赝势的。而你的关键词目前用的是 Gen 也就是你只读取了Cu的电子基组而省略了赝势。这样基组不完整收敛不了很正常。你应该把Gen 换成 GenECP应该至少qc下能收敛。

作者
Author: 化学girl 时间: 2026-1-13 11:21

wal 发表于 2026-1-13 09:53
是xqc收敛失败，本质上还是SCF不收敛，参考SCF不收敛解决方案http://sobereva.com/61
题外话，为什么直接 ...

谢谢解答，也按照论坛里面的步骤尝试过了，但是还是不收敛，才用了qc
因为这个构型是解单晶以后的构型了，处于能量最低的位点，合作的小伙伴说不需要进行结构优化了，所以就直接进行了频率计算。

作者
Author: 化学girl 时间: 2026-1-13 11:22

paramecium86 发表于 2026-1-13 10:35
你如果是要用lanl2dz 那么对于Cu是有赝势的。而你的关键词目前用的是 Gen 也就是你只读取了Cu的电子基组而 ...

谢谢解答，但是不太明白，如果是这个原因的话，那为什么我之前算的Cu12团簇就能够收敛呢？

作者
Author: wal 时间: 2026-1-13 11:40
本帖最后由 wal 于 2026-1-13 11:43 编辑

化学girl 发表于 2026-1-13 11:21
谢谢解答，也按照论坛里面的步骤尝试过了，但是还是不收敛，才用了qc
因为这个构型是解单晶以后的构型了 ...

如果是做波函数分析，那你直接用单晶不做结构优化还说得过去，但做freq必须必须重新优化吧，不然振动分析是无意义的。
我楼下说的没错，我看了你的输入，不止是genecp的问题，你基组写的也有问题，目前你的写法铜根本没给基组定义。你此前的计算如果也是用的这个写法，即使收敛了估计也是有问题的。应该是这样的

坐标...
H 4.73970000 12.56840000 19.84830000
C S H 0
6-31G*
****
Cu 0
LANL2DZ
****
Cu 0
LANL2DZ

复制代码

作者
Author: paramecium86 时间: 2026-1-13 12:47
本帖最后由 paramecium86 于 2026-1-13 12:49 编辑

化学girl 发表于 2026-1-13 11:22
谢谢解答，但是不太明白，如果是这个原因的话，那为什么我之前算的Cu12团簇就能够收敛呢？

能收敛也不代表就对。无论收不收敛你的体系首先得保证用的基组是正确的。只是说你的体系更大的时候基组用错了，不收敛的概率就更大。而且下面混合基组的格式也写法有问题。不同基组之间你多了一个空行。这样就根本也读取不到Cu的基组。

作者
Author: 化学girl 时间: 2026-1-13 16:31

wal 发表于 2026-1-13 11:40
如果是做波函数分析，那你直接用单晶不做结构优化还说得过去，但做freq必须必须重新优化吧，不然振动分析 ...

好的，谢谢您，那请问我这样算完的结构进行前线轨道分析的话，数据是可靠的吗？
我也再修改一下试着再算一算。

作者
Author: 化学girl 时间: 2026-1-13 16:34

paramecium86 发表于 2026-1-13 12:47
能收敛也不代表就对。无论收不收敛你的体系首先得保证用的基组是正确的。只是说你的体系更大的时候基组 ...

嗷嗷，好的，谢谢您~
一个是需要修改写法，请问一下那更大的体系的时候是不是还需要考虑修改基组呀？使用def2-TZVP是不是更合适一点呢？

作者
Author: pikachuupup 时间: 2026-1-13 16:40

化学girl 发表于 2026-1-13 16:34
嗷嗷，好的，谢谢您~
一个是需要修改写法，请问一下那更大的体系的时候是不是还需要考虑修改基组呀？使 ...

横向比较时不应更换基组，如果只进行波函数分析，单点直接都用def2-TZVP显然更方便，也不用设置混合基组

作者
Author: sgwzq 时间: 2026-1-13 16:43
本帖最后由 sgwzq 于 2026-1-13 17:23 编辑

化学girl 发表于 2026-1-13 16:34
嗷嗷，好的，谢谢您~
一个是需要修改写法，请问一下那更大的体系的时候是不是还需要考虑修改基组呀？使 ...

你是说几何优化振动分析还是单点能？前面那两个是“按能量下降的方向找极小点结构”（能量数值不是主要研究对象），后面那个是精确算电子能量。

如果几何优化的话，你可以简单算一个小一点的分子，看看TZ和DZ会不会明显影响结构（键长键角等，再看看小基组下有没有虚频），并且看看耗时差距。
想看耗时的话，在Linux提交计算任务之前，加个time命令，比如原来是

g09 xx xx

复制代码

加个time变成

time g09 xx xx

复制代码

在计算完成时会同步输出耗时

如果耗时差距特别大，几何结构优化相差不多，那么优化就用个DZ就可以了（比如SVP等优化，TZVP算单点）
算单点能的时候再都换成大基组。

如果小基组算出来结构结果完全不可信，那优化就必须用大基组来充分描述原子外电子分布情况。
基组做的事就是，别人预先定义好每个原子外的电子长成啥形。
泛函做的事是别人预先提供一个“把原子周围电子密度（或者说波函数）组合成分子周围电子密度“的算法。
高斯做的事就是按这个算法，把电子密度组合一下。
至于几何优化/振动分析，可以看这个

作者
Author: wal 时间: 2026-1-13 16:45

化学girl 发表于 2026-1-13 16:31
好的，谢谢您，那请问我这样算完的结构进行前线轨道分析的话，数据是可靠的吗？
我也再修改一下试着再算 ...

如果你之前的基组也是这么写的，并不可靠，铜根本没基组。

作者
Author: 化学girl 时间: 2026-1-13 18:44

wal 发表于 2026-1-13 16:45
如果你之前的基组也是这么写的，并不可靠，铜根本没基组。

我重新修改了计算，但是直接就报错了，请问为什么会是这样呀？

作者
Author: 化学girl 时间: 2026-1-13 18:46

sgwzq 发表于 2026-1-13 16:43
你是说几何优化振动分析还是单点能？前面那两个是“按能量下降的方向找极小点结构”（能量数值不是主要 ...

好的，谢谢您，我再好好学习一下~

作者
Author: wal 时间: 2026-1-13 18:53

化学girl 发表于 2026-1-13 18:44
我重新修改了计算，但是直接就报错了，请问为什么会是这样呀？

赝势没有****，你看6l我发的那个

作者
Author: sobereva 时间: 2026-1-14 06:13

化学girl 发表于 2026-1-13 18:44
我重新修改了计算，但是直接就报错了，请问为什么会是这样呀？

认真把此文看了
在Gaussian中自定义范围分离泛函的方法
http://sobereva.com/550（http://bbs.keinsci.com/thread-17256-1-1.html）

都没做优化上来就做freq根本毫无意义

下文说了，纯泛函比杂化泛函普遍更难收敛。建议改用PBE0甭用PBE
解决SCF不收敛问题的方法
http://sobereva.com/61
上文也说了，别习惯性带着xqc，非要靠qc才能收敛的情况很可能收敛到了不稳定波函数，至少你知道什么是波函数稳定性测试的前提下再用

16核跑133原子体系的freq居然才给10GB内存不可思议，认真看下文了解内存分配常识
Gaussian的安装方法及运行时的相关问题
http://sobereva.com/439（http://bbs.keinsci.com/thread-10814-1-1.html）

强烈建议好好看看此文
谈谈学量子化学如何正确地入门
http://sobereva.com/355（http://bbs.keinsci.com/thread-4447-1-1.html）
一定要先有足够的知识基础再做计算

作者
Author: 化学girl 时间: 2026-1-14 08:42

wal 发表于 2026-1-13 18:53
赝势没有****，你看6l我发的那个

好呢，这回明白了，太感谢啦！

作者
Author: 化学girl 时间: 2026-1-14 09:04

sobereva 发表于 2026-1-14 06:13
认真把此文看了
在Gaussian中自定义范围分离泛函的方法
http://sobereva.com/550（http://bbs.keinsci. ...

好的，谢谢老师，，我再好好读读这几篇文章。现在还没有确定全面投入理论计算，就是在自己摸索一下，后期有需要的话也非常想跟着老师系统的学习一下。

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