计算化学公社

标题: 在gaussian中哪些计算和IRC一样需要初始猜测？ [打印本页]

作者
Author: 元气蛋 时间: 2019-7-27 10:17
标题: 在gaussian中哪些计算和IRC一样需要初始猜测？
各位老师好，在近期的gaussian学习中师姐说计算IRC时，尽量不采用直接输入过渡态优化后out文件中的坐标，而是利用geom=checkpiont guess=read关键词，读取过渡态优化后的chk文件，因为计算IRC时，需要过渡态优化后对力常数的初始猜测，所以我想问一下，在高斯中除了计算IRC还有哪些计算需要初始猜测？

作者
Author: zjxitcc 时间: 2019-7-27 10:41
本帖最后由 zjxitcc 于 2019-7-27 10:44 编辑

几乎所有的量化计算（SCF迭代、结构优化、CC振幅迭代）都需要初始猜测，只是有的计算对初始猜测不那么敏感（如常见闭壳层有机物分子），有的很敏感（过渡金属、自由基）。“不敏感”意为初猜的质量好坏一般不影响结果，只影响过程（即时间长短、收敛快慢）。

你师姐说的不全对，也不全错。
使用irc=rcfc geom=check guess=read，其中guess=read读进上一步收敛的波函数，因此能1圈立即收敛；rcfc表示读进上一步结构优化收敛的二阶导数，不用重新算。这两步都是节省时间用的，一般不影响结果。你也可以写irc=calcfc，且不加guess=read，徒增机时罢了。“计算IRC时，需要过渡态优化后对力常数的初始猜测”这话不妥，TS找到以后做freq任务，得到的二阶导数是精确的，后面IRC任务只不过把它读进来了而已。

作者
Author: 元气蛋 时间: 2019-7-27 10:57

zjxitcc 发表于 2019-7-27 10:41
几乎所有的量化计算（SCF迭代、结构优化、CC振幅迭代）都需要初始猜测，只是有的计算对初始猜测不那么敏感 ...

谢谢老师的指点

作者
Author: sobereva 时间: 2019-7-27 22:04
你师姐根本没有正确理解IRC的基本算法，或者是你理解错她的意思了

guess=read根本和什么“对力常数的初始猜测”毫无关系。而且对于IRC而言，也根本没有什么所谓的“对力常数的初始猜测”。IRC需要提供初始Hessian，Hessian要么直接calcfc计算，要么rcfc从含有当前结构精确Hessian的chk文件中读取。跑IRC根本就不能用近似的Hessian（和opt、opt=TS的情况不同），何来“初始猜测”？（行业内正规叫法是“初猜”）

“尽量不采用直接输入过渡态优化后out文件中的坐标”也根本没法这么说，倒不如说绝大多数情况都是这么来做的。

相关知识仔细看此文，弄清楚IRC该怎么跑、需要什么信息
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

作者
Author: 元气蛋 时间: 2019-8-12 08:38

sobereva 发表于 2019-7-27 22:04
你师姐根本没有正确理解IRC的基本算法，或者是你理解错她的意思了

guess=read根本和什么“对力常数的初 ...

谢谢社长的指点，是我理解有误

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