计算化学公社

标题: 机械力化学中的 CoGEF方法如何实现 [打印本页]

作者
Author: chuan437 时间: 2022-12-4 23:34
标题: 机械力化学中的 CoGEF方法如何实现
本帖最后由 chuan437 于 2022-12-4 23:35 编辑

现在研究聚合物机械力化学化学有一种方法叫做CoGEF方法。如下图所示，固定这个大分子两端的碳原子。然后不断增加这两个碳原子之间的距离，直到把这个分子最中间的碳碳键拉断。
请问如何用高斯程序实现这样的CoGEF方法。我现在知道势能面扫描只是改变两个成键原子之间的键长距离，对于这种方法是不是需要先把这个大分子分成两个结构固定的基团（以被拉断的键为界限）？
请教各位朋友，谢谢。

作者
Author: Kalinite 时间: 2022-12-5 00:15
本帖最后由 Kalinite 于 2022-12-5 00:18 编辑

大多数文献里CoGEF一般是用spartan，可以在分子力学层面做。
要用Gaussian做的话，可以做柔性扫描，把锚定的两个原子（而不是基团）距离作为扫描变量，步长设置的小一些，计算量应该是比较大的。另外就是Gaussian只能固定步长做。
严格考察施加外力的情况，还得是ORCA。

不过我没具体用过spartan，最近合作者问了我这方面的东西，我个人以为分子力学做是不太合适的，毕竟参数都是针对比较靠近平衡位置进行拟合，对于拉开一定距离，偏离平衡位置的势能面描述的会有问题，对于可能涉及交叉的问题就更大了。

作者
Author: sobereva 时间: 2022-12-5 02:09
“能面扫描只是改变两个成键原子之间的键长距离”完全没这回事。明明任意两个原子间的距离都可以扫描，比如北京科音初级量子化学培训班（http://www.keinsci.com/workshop/KEQC_content.html）里的例子

(, 下载次数 Times of downloads: 8)

Spartan在这方面没有任何优势，而且又收费又冷门
这种问题不适合分子力场。本来精度就远不如常用的DFT，而且大部分常见力场对化学键都没有正确的解离行为。

作者
Author: ClaudiaZ 时间: 3 day ago
同问，如何在高斯中实现用CoGEF，在原有基础上对两个原子额外施加1.2埃伸长长度的机械力

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