计算化学公社

标题: 直线结构的过渡金属碳团簇优化，虚频始终消不掉 [打印本页]

作者
Author: R-S 时间: 2026-1-27 00:21
标题: 直线结构的过渡金属碳团簇优化，虚频始终消不掉
优化一个过渡金属碳团簇，直线构型，无法收敛使用xqc收敛后，有虚频，尝试了各种方法始终没办法将这个虚频消去。有时候甚至会出现5.6个破坏碳链线性的虚频，但是调整虚频结构后，优化始终不收敛。甚至会出现当前点群C1，显示有虚频默认点群Cinfv的情况。计算后的文件已上传，实在没办法了，不知道是不是方法选择问题，试了几个其他的b3lyp，tpssh，mn15l都是一样的问题。那么这个结构到底是不是一个稳定结构呢？能否直接忽视虚频呢？

作者
Author: wal 时间: 2026-1-27 00:34
确定是直线吗？如果怎么优化都是破坏对称性的虚频，还这么大，那可能确实不是直线的。
另外xqc收敛之后stable了吗，qc收敛到的波函数可能不稳定

作者
Author: UW_0728. 时间: 2026-1-27 05:48
本帖最后由 UW_0728. 于 2026-1-27 05:52 编辑

你的电荷和自旋多重度设置确定没问题吧
不清楚你在这里带弥散的意义是什么，本来加弥散函数就会增大SCF的收敛难度
解决SCF不收敛不要一上来就用xqc
输出文件里看你用的BPW91，搜了下，应该是不推荐的。可以试试PBE0-D3(BJ)

作者
Author: R-S 时间: 2026-1-27 09:40

UW_0728. 发表于 2026-1-27 05:48
你的电荷和自旋多重度设置确定没问题吧
不清楚你在这里带弥散的意义是什么，本来加弥散函数就会增大SCF的 ...

加弥散是因为负离子优化带了弥散为保持结果可对比，那我试一下不加弥散，但是负离子加弥散一样是这样的问题，其他的方法都没作用，只有xqc能优化的下去。

作者
Author: 122332 时间: 2026-1-27 16:58
试试用PBE0，scf不收敛可以试试用精确计算Hessian矩阵calcall

作者
Author: R-S 时间: 2026-1-27 20:59

122332 发表于 2026-1-27 16:58
试试用PBE0，scf不收敛可以试试用精确计算Hessian矩阵calcall

几乎所有的方法都试过了，PBE0也是一样的情况，关键是我不让他高对称性都做不到，即使用C1优化，他都会给你优化到默认Cinfv 然后告诉你有两个虚频。太折磨人了

作者
Author: 122332 时间: 2026-1-27 21:29
你这个虚频两个简并，而且500多波数那个虚频会带着链上很多原子振动，scf=xqc收敛到的波函数应该不是稳定的波函数，这个结构可能也不稳定

作者
Author: bxyao 时间: 2026-1-28 11:18

R-S 发表于 2026-1-27 20:59
几乎所有的方法都试过了，PBE0也是一样的情况，关键是我不让他高对称性都做不到，即使用C1优化，他都会给 ...

你确定这个分子是真实存在的么？这么多C还是直线型的分子么？

作者
Author: Luciferly 时间: 2026-1-28 11:42
负离子你的电荷怎么设置的是0

作者
Author: KazusaT 时间: 2026-1-28 12:09
1. 你的log文件中有

-------------------
Title Card Required
-------------------
Symbolic Z-matrix:
Charge = 0 Multiplicity = 2

建议你再检查一下电荷和自旋多重度的设置
2. 如前面各位老师所说，很可能收敛到了不稳定的波函数，可以尝试先用小基组、不加弥散的基组找到一个稳定的波函数并读入来做后续计算

作者
Author: R-S 时间: 2026-1-31 17:42

bxyao 发表于 2026-1-28 11:18
你确定这个分子是真实存在的么？这么多C还是直线型的分子么？

这就是在计算来确定他的结构，不是知道他的结构再去算的

作者
Author: R-S 时间: 2026-1-31 17:42

Luciferly 发表于 2026-1-28 11:42
负离子你的电荷怎么设置的是0

这个算的是中性的结构，负离子和中性都存在这样的问题

作者
Author: R-S 时间: 2026-1-31 17:45

KazusaT 发表于 2026-1-28 12:09
1. 你的log文件中有

建议你再检查一下电荷和自旋多重度的设置

这个是中性结构，波函数都是稳定的，连降低对称性都做不到，因为你降低了之后就会被优化到默认Cinfv点群，然后同时告诉你有虚频。在尝试其他办法了，只要是没虚频是什么点群算什么点群。

作者
Author: Luciferly 时间: 2026-2-3 12:04

R-S 发表于 2026-1-31 17:42
这个算的是中性的结构，负离子和中性都存在这样的问题

你这个结构很违反化学直觉，这么多连续的积累烯烃应该不会很稳定。而且烯烃的金属配合物很少端基配位，你看看有没有侧基配位的情况。

作者
Author: UW_0728. 时间: 2026-2-3 13:33
本帖最后由 UW_0728. 于 2026-2-3 13:56 编辑

你这个结构很怪异，是否有可能是静态相关强，杂化泛函本身就描述不好？

我试了试PBE纯泛函，这个结构反而能够很好地优化收敛且无虚频（用的关键词是pbepbe/tzvp/fit；保险起见默认加了D3-(BJ)色散校正，尽管对于该体系应该没啥必要）。但是我也不太敢打包票说这里用纯泛函优化就一定可行。

P.S. 奇怪的是，拿PBE纯泛函收敛的结果用PBE0-D3(BJ)优化也很快收敛了，但用初始结构就连SCF都不收敛

(, 下载次数 Times of downloads: 1) (, 下载次数 Times of downloads: 1)

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