计算化学公社

标题: 一批Pt配合物的结构优化问题求助，不是没法收敛就全是虚频 [打印本页]

作者
Author: 11111111m 时间: 2024-9-13 10:03
标题: 一批Pt配合物的结构优化问题求助，不是没法收敛就全是虚频
求助各位老师，我在优化一批pt配合物，基组选用LANL2DZ/6-13G(d,p)，泛函分别试了PBE0和B3LYP，PCM溶剂环境，但发现优化出来的结构不是没法收敛就全是虚频，尝试过calcall仍一直不收敛。感觉结构本身好像没有什么问题，麻烦各位老师帮忙看一下，输出文件太大上传不了

。

作者
Author: Loading0760 时间: 2024-9-13 10:19
溶剂化不好收敛的话，可以试试先去掉溶剂化，结构优化，然后读取真空状态下收敛的波函数作为输入，增加溶剂，进行计算。
过渡金属建议用lanl2TZ，而不是LANL2DZ。
你这个是第一步就开始震荡了，没有输出的日志文件，我从中得到的信息不够多。参考社长的文章http://sobereva.com/61。
顺便吐槽一句你的输入文件，96个核，才10GB的内存。

作者
Author: 11111111m 时间: 2024-9-13 10:23

Loading0760 发表于 2024-9-13 10:19
溶剂化不好收敛的话，可以试试先去掉溶剂化，结构优化，然后读取真空状态下收敛的波函数作为输入，增加溶剂 ...

尝试去过溶剂化，也是一样不收敛

作者
Author: sobereva 时间: 2024-9-13 15:14
仔细看此文
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164
应当首先考虑到减小步长上限，上文明确说了

当前情况应当用PBE0-D3(BJ)
区区一个Pt的情况，Pt绝对不该用过时的Lanl2DZ，仔细看下文
谈谈赝势基组的选用
http://sobereva.com/373（http://bbs.keinsci.com/thread-5625-1-1.html）

二氯甲烷这么低极性的溶剂环境做opt freq没绝对必要带溶剂模型

作者
Author: 11111111m 时间: 2024-9-14 10:32
本帖最后由 11111111m 于 2024-9-14 10:34 编辑

sobereva 发表于 2024-9-13 15:14
仔细看此文
量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

sob老师，我用了LanL2TZ，减少了步长上限，但优化出来的结构仍然全是虚频，好奇怪.
#p opt=maxstep=5 freq pbe1pbe/genecp em=GD3BJ，基组是LanL2TZ/6-13g(d,p)

作者
Author: sobereva 时间: 2024-9-15 00:42

11111111m 发表于 2024-9-14 10:32
sob老师，我用了LanL2TZ，减少了步长上限，但优化出来的结构仍然全是虚频，好奇怪.
#p opt=maxstep=5 fr ...

Gaussian中几何优化收敛后Freq时出现NO或虚频的原因和解决方法
http://sobereva.com/278（http://bbs.keinsci.com/thread-633-1-1.html）

那叫6-31G(d,p)

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