计算化学公社

标题: 直线Cs结构虚频无法消除，陷入循环 [打印本页]

作者
Author: R-S 时间: 2025-5-9 20:04
标题: 直线Cs结构虚频无法消除，陷入循环
本帖最后由 R-S 于 2025-5-9 21:10 编辑

在计算一个Cs点群的中性团簇时，有一个破坏直线对称性的虚频，沿虚频调整结构降低对称性后重新优化，会优化到最初的高对称性结构，提示默认点群Cs，但是仍然有这个虚频，等于转了一圈又回去了，这样就陷入到：调整虚频优化会得到更高对称性的结构，但是会有破坏对称性的虚频这样的循环，始终得不到没虚频的结构，甚至连降低对称性的结构都得不到，因为低对称性的结构优化会变成更高点群同时提示有虚频。这种情况该怎么处理呢？计算时遇到了非常多类似的情况，这种原因是什么呢？有什么解决的主旨吗？参考Sober的帖子也做了很多尝试，但是始终存在这样的死循环。
(, 下载次数 Times of downloads: 19) (, 下载次数 Times of downloads: 17)
第一个图就是调整虚频优化后的结果，明显又回到了Cs对称性，但是会有图二的虚频，可以看到这个虚频还很大。

结构和虚频如下：
(, 下载次数 Times of downloads: 16)
虚频就是在碳链部分，碳链的C-C-C角度偏离180° 大概0.几度，所以优化得到的是Cs而不是Cinfv的对称性。但是即使是Cs也还是有一个非常大的虚频，优化用的就是最简单的关键词：#p opt freq 6-311+g(d) bpw91

作者
Author: Stardust0831 时间: 2025-5-9 20:24
虚频原因有很多，当前没有给出原始文件，难以做进一步判断。
建议按照这篇博文里提到的排查一遍《Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法》（http://sobereva.com/278）
同时我希望确认你没有使用类似M06-2X这种极度依赖积分格点精度且容易出虚频的泛函。建议您上传原始文件，或者至少截图一下已有的关键词。

作者
Author: R-S 时间: 2025-5-9 21:00

Stardust0831 发表于 2025-5-9 20:24
虚频原因有很多，当前没有给出原始文件，难以做进一步判断。
建议按照这篇博文里提到的排查一遍《Gaussian ...

已上传，计算用的是bpw91/6-311+g(d)

作者
Author: wzkchem5 时间: 2025-5-9 21:01
此外确定一下计算是否收敛到了正确的波函数，频率计算的波函数和消虚频时收敛到了波函数是否一致

作者
Author: R-S 时间: 2025-5-9 21:05

wzkchem5 发表于 2025-5-9 21:01
此外确定一下计算是否收敛到了正确的波函数，频率计算的波函数和消虚频时收敛到了波函数是否一致

您说的确定波函数是否正确是做波函稳定性检测吗？调整虚频优化得到的波函数是稳定的，也就是得到高对称但有虚频的结构的那一步优化。

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-5-9 21:19
第一眼看还以为是[C9N]-离子，结果末端居然是个裸露的Co原子？电荷和自旋多重度具体是多少，确定没有多参考特征么？关于波函数稳定性检测和几何优化的问题，看http://bbs.keinsci.com/thread-45423-1-1.html二楼。另外，直线型分子得特别注意下优化一开始冗余内坐标的构建。

当前泛函和基组不知道咋选择的，可以看看http://sobereva.com/272和http://sobereva.com/336换一换，比如TPSSh（或PBE0）/def2-TZVP之类的（暂时不确定弥散函数有没有必要）

作者
Author: R-S 时间: 2025-5-9 21:36

Uus/pMeC6H4-/キ发表于 2025-5-9 21:19
第一眼看还以为是[C9N]-离子，结果末端居然是个裸露的Co原子？电荷和自旋多重度具体是多少，确定没有多参考 ...

非常感谢您的建议，我也觉得这个泛函和基组最当前的体系值得商榷，我修改了关键词添加了maxstep=1，并提高积分格点到superfine，这个虚频被消掉了，但是对称性变为了C1，问题是默认高阶的对称性还是Cs，我估计改为Cs对称性还是会有虚频。

作者
Author: wzkchem5 时间: 2025-5-10 08:58

R-S 发表于 2025-5-9 21:05
您说的确定波函数是否正确是做波函稳定性检测吗？调整虚频优化得到的波函数是稳定的，也就是得到高对称但 ...

包括但不限于波函数检测。一个体系的稳定波函数未必是唯一的，光检测波函数稳定性是不够的，必须通过前线轨道、自旋密度、自旋布居等确认你是否收敛到了你想要的那个特定的稳定波函数

作者
Author: R-S 时间: 2025-5-12 09:45

wzkchem5 发表于 2025-5-10 08:58
包括但不限于波函数检测。一个体系的稳定波函数未必是唯一的，光检测波函数稳定性是不够的，必须通过前线 ...

非常感谢您的回复，有个问题想请教，那如果得到的波函数并不是我想要，我如何让他获得一个新的稳定波函数呢？重新优化让他生成吗，还是有什么关键词可以用，谢谢您！

作者
Author: wzkchem5 时间: 2025-5-12 10:27

R-S 发表于 2025-5-12 09:45
非常感谢您的回复，有个问题想请教，那如果得到的波函数并不是我想要，我如何让他获得一个新的稳定波函数 ...

分析一下你希望收敛到的波函数，和当前收敛的波函数，轨道占据数有哪些不同。也就是哪些本该占据的轨道是空轨道，哪些本该不占据的轨道是占据轨道。
如果轨道占据数错误的轨道不多，手动交换相应轨道的顺序，然后再收敛。收敛时可以加一个level shift防止轨道再收敛回原来的错误解。
如果轨道占据数错误的轨道很多，可以从头开始，用片段波函数、flip spin等方法准备合适的初猜，再收敛

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