计算化学公社

标题: Gaussian优化获得过渡态，验证（IRC和freq计算）时出现问题 [打印本页]

作者
Author: xrpeng 时间: 2024-8-30 16:48
标题: Gaussian优化获得过渡态，验证（IRC和freq计算）时出现问题
本帖最后由 xrpeng 于 2024-8-30 16:47 编辑

求助各位！
我想得到一个小分子解离的IRC，首先在SCSMP2/def2-TZVPP理论水平下优化（MP2 def2TZVPP OPT(TS NoEigenTest calcall Maxcyc=200) 6D 10F IOp(3/125=0333312000)得到了有一个虚频的过渡态，计算频率（MP2 def2TZVPP Freq Nosymm 6D 10F IOp(3/125=0333312000)所得虚频与优化文件中的一致。用这个结构Forword和Reverse跑IRC发现能得到反应物，但产物与预期不一致，而拿这个所谓的产物P0（未完全解离）继续优化却得到了我想要的产物P（完全解离）。
我测试了一系列步长（1-9,20、30、40、50）都无法跨越P0这个结构，我所有IRC的步数都设置的足够大。
后来担心是理论水平的问题，使用DFT优化过渡态，出现的问题是我使用了calcall关键词优化（B3LYP em=GD3 def2TZVPP OPT(TS NoEigenTest calcall Maxcyc=200) 6D 10F）得到了一个虚频的过渡态（-314.6279），但频率计算时却有两个虚频。我现在不知道如何确定是我找的过渡态有问题还是计算的关键词设置有问题

作者
Author: sobereva 时间: 2024-8-31 05:58
用6D 10F又浪费时间又没任何好处
SCS-MP2用于几何优化没有什么好处，又贵而且还未必有合适的泛函靠谱。用得动MP2的情况没理由不用双杂化泛函，哪怕是很老的B2PLYP。SCS-MP2已经没有明显实际意义了，只有历史意义，启发了后来研究者把SCS的思想引入双杂化泛函里。

单纯说你的第一个问题，优化的时候把步长上限改小，以优化的结果判断未完全解离的结构在当前级别下是否存在对应的极小点。

作者
Author: xrpeng 时间: 2024-9-6 10:10

sobereva 发表于 2024-8-31 05:58
用6D 10F又浪费时间又没任何好处
SCS-MP2用于几何优化没有什么好处，又贵而且还未必有合适的泛函靠谱。用 ...

谢谢社长提醒！其他关键词的设置主要是考虑到后续的分析，我还是想拯救一下
对于未完全解离的P0我后来进行了一些优化尝试（1楼提到的优化关键词设置为opt(Maxcyc=200, calcall)）：
01.com:# MP2 def2TZVPP opt(GDIIS,Maxstep=5, notrust,calcall,maxcyc=100)  6D 10F IOp(3/125=0333312000)
02.com:# MP2 def2TZVPP opt(MaxStep=5,NoTrustUpdate,Cartesian)  6D 10F IOp(3/125=0333312000)
03.com:# MP2 def2TZVPP opt(Maxstep=5, notrust,calcall)  6D 10F IOp(3/125=0333312000)
04.com:# MP2 def2TZVPP opt(GDIIS,Maxstep=5, notrust,calcall)  6D 10F IOp(3/125=0333312000)
04_continue.com:# MP2 def2TZVPP opt(GDIIS,Maxstep=5, notrust,calcall)  6D 10F IOp(3/125=0333312000)
05.com:# MP2 def2TZVPP opt(Maxstep=5, notrust,calcall,maxcyc=200)  6D 10F IOp(3/125=0333312000)
上述计算中02、03、04、04c都存在相同的报错
Optimization stopped.
-- Number of steps exceeded,  NStep=  22
-- Flag reset to prevent archiving.
01、05的任务都正常结束，频率计算没有虚频，最终的产物结构是完全解离的。只是产物最后的两个解离单体摆放位置和角度不一致，但距离都挺远的，能量相差0.0x kcal/mol，事实上在优化过程中单体还发生了翻转
请问社长上述测试是否说明IRC得到的结构远非极小值点，我如何设置IRC关键词尽可能靠近解离产物呢
或者如果我不需要分析整个IRC，我已经获得的这段IRC是否能拿来进行进一步分析呢

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