计算化学公社

标题: IRC两端结构跑产物和反应物结构的参数设置求助 [打印本页]

作者
Author: gjf 时间: 2023-4-27 15:15
标题: IRC两端结构跑产物和反应物结构的参数设置求助
各位老师，我想根据IRC两端结构，优化出反应物和产物的结构，跑IRC时的参数设置为#p IRC=(rcfc,LQA,maxpoints=300) ub3lyp/genecp em=gd3 guess=read。在我对两端结构进行结构优化时对参数设置有了点疑问，请各位老师看一下，哪一种参数设置比较合适。谢谢各位老师
（1）#p opt=(calcfc,cartesian,gdiis,) freq ub3lyp/genecp em=gd3
（2）#p opt=(calcfc,cartesian,gdiis,maxstep=20,NoTrust) freq ub3lyp/genecp em=gd3

作者
Author: wzkchem5 时间: 2023-4-27 20:15
都可以，很难说哪种更好（但第一种写法你多写了一个逗号）。如果你有计算类似体系的经验，选了其中一种发现收敛太慢，那就用另一种；否则随便选其中的一种

作者
Author: sobereva 时间: 2023-4-28 01:34
cartesian跟GDIIS根本不兼容，此时用cartesian是无效的
对于分子体系，本来如果没有特殊情况就强烈不建议加cartesian，只会让收敛更慢

作者
Author: gjf 时间: 2023-4-28 08:45

wzkchem5 发表于 2023-4-27 20:15
都可以，很难说哪种更好（但第一种写法你多写了一个逗号）。如果你有计算类似体系的经验，选了其中一种发现 ...

老师，两种速度都差不多，那我就选择第二种是看到一个帖子下回复，说加入maxstep=20，notrust可以增加稳妥性

作者
Author: gjf 时间: 2023-4-28 08:48

sobereva 发表于 2023-4-28 01:34
cartesian跟GDIIS根本不兼容，此时用cartesian是无效的
对于分子体系，本来如果没有特殊情况就强烈不建议 ...

老师，我计算的体系存在原子冻结，是在笛卡尔坐标控制冻结的，所以在计算的时候加入了cartesian

作者
Author: sobereva 时间: 2023-4-28 11:52

gjf 发表于 2023-4-28 08:48
老师，我计算的体系存在原子冻结，是在笛卡尔坐标控制冻结的，所以在计算的时候加入了cartesian

那就甭用GDIIS

作者
Author: wzkchem5 时间: 2023-4-28 16:01

gjf 发表于 2023-4-28 01:45
老师，两种速度都差不多，那我就选择第二种是看到一个帖子下回复，说加入maxstep=20，notrust可以增加稳 ...

增加稳妥性不假，但是当收敛本来就很稳妥时，maxstep=20，notrust可能使收敛变慢。所以要在收敛速度和收敛概率之间做一个权衡，要么收敛概率大而需要的步数多，要么需要的平均步数较少但是收敛概率低

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