求助含23个Cu原子的团簇的计算过程中出现l508报错

化学girl

老师们好，在过渡金属团簇的计算过程中，出现了l508的报错，这是有23个Cu原子的团簇
之前还算了有12个Cu原子的团簇就收敛了，用相同的方法算这个大一点的就不收敛，还有一个25个Cu的，也不收敛

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

wal

是xqc收敛失败，本质上还是SCF不收敛，参考SCF不收敛解决方案http://sobereva.com/61
题外话，为什么直接做freq？

paramecium86

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

化学girl

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

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

化学girl

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

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

wal

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

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

1. 坐标...
   
2. H                  4.73970000   12.56840000   19.84830000
   
3. 
   
4. C S H 0
   
5. 6-31G*
   
6. ****
   
7. Cu 0
   
8. LANL2DZ
   
9. ****
   
10. 
    
11. Cu 0
    
12. LANL2DZ
    
13. 
    
14. 
    

复制代码

paramecium86

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

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

化学girl

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

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

化学girl

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

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

pikachuupup

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

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

sgwzq

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

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

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

1. g09 xx xx

复制代码

加个time变成

1. time g09 xx xx

复制代码

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



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

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

wal

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

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

化学girl

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

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

化学girl

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

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

wal

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

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