计算化学公社

标题: 高斯进行柔性扫描和单独做限制性优化的区别？ [打印本页]

作者
Author: JoeyssaWang 时间: 2020-7-16 10:59
标题: 高斯进行柔性扫描和单独做限制性优化的区别？
用Gaussian做柔性扫描和单独做限制性优化opt+freq的区别是什么？
（1）柔性扫描的方法：关键词opt=modredundant
（2）限制性优化：取出柔性扫描的一个结构，单独做opt+freq，是用0/-1的方法限制原子坐标
问题1：这两者计算的力常数不同，做freq得到的力常数应该是对Hessian矩阵做对角化得到，那柔性扫描的力常数是怎么得到的呢？

问题2：这么做限制性优化，opt+freq得到的热力学参数是否准确？

感谢回复！

作者
Author: zjxitcc 时间: 2020-7-16 11:08
问题中有矛盾点和歧义点。
（1）“单独做限制性优化”含义模糊，说清楚到底用的哪种技巧做的限制性优化，技巧有很多，还有-1/0冻结坐标的这种也可叫限制性优化；
（2）力常数与限制性优化是两个问题，单独做限制性优化、不做freq计算并没有力常数。另外，通常应该去掉坐标限制再做freq计算（除非你就是想要带坐标限制的力常数）

作者
Author: JoeyssaWang 时间: 2020-7-16 11:17

zjxitcc 发表于 2020-7-16 11:08
问题中有矛盾点和歧义点。
（1）“单独做限制性优化”含义模糊，说清楚到底用的哪种技巧做的限制性优化， ...

不好意思，我准确说一下问题：
（1）单独做限制性优化，就是用 0 -1的方法限制的原子；
（2）单独做限制性优化时，同时做的opt+freq, 我需要的是在限制坐标下的力常数

作者
Author: sobereva 时间: 2020-7-16 18:30
“柔性扫描的力常数”这描述完全是莫名其妙，柔性扫描和振动分析又没直接关系

做opt freq的时候，如果冻结了某些变量，则freq时Hessian矩阵中相应的矩阵元会被清零。

限制性优化显然之后做freq得到的热力学量在原理上不准确，除非你限制的是与你研究的问题无关的区域。

作者
Author: JoeyssaWang 时间: 2020-7-16 19:52
本帖最后由 JoeyssaWang 于 2020-7-16 20:04 编辑

sobereva 发表于 2020-7-16 18:30
“柔性扫描的力常数”这描述完全是莫名其妙，柔性扫描和振动分析又没直接关系

做opt freq的时候，如果冻 ...

谢谢老师回复，我确实搞错了，我其实看的是柔性扫描过程，每一个原子的受力，就是输出文件里面Forces (Hartrees/Bohr)下面的坐标参数，我误当成了“力常数”。和振动分析得到的力常数搞混了，不好意思。
柔性扫描过程中（拉伸两个原子的距离），每一步scf都有Forces (Hartrees/Bohr)，这时所限制的原子受力就不为零，其他原子为零，这该怎么理解呢？是不是只有在第一个出发点的时候时极小点，这时的E对R的一阶导数为零，然后偏移了这个位置就都不为零了。

作者
Author: sobereva 时间: 2020-7-17 09:52

JoeyssaWang 发表于 2020-7-16 19:52
谢谢老师回复，我确实搞错了，我其实看的是柔性扫描过程，每一个原子的受力，就是输出文件里面Forces (Ha ...

是

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