计算化学公社

标题: ORCA跑AIMD可不可以分别设两个片段电荷和多重度 [打印本页]

作者
Author: -HMDS- 时间: 2023-12-1 04:14
标题: ORCA跑AIMD可不可以分别设两个片段电荷和多重度
各位大佬，ORCA可不可以输入两个片段进行AIMD计算，比如算水合质子和羟基自由基二聚的动力学，分别设定初始是两部分的电荷多重度，或者从其他文件读波函数作为初始状态的波函数再跑动力学，如果没有这个限制往往电荷会分散在两个片段上，这样计算还有意义吗。

作者
Author: sobereva 时间: 2023-12-1 05:59
没法设。顶多提供一个初猜波函数（可以利用Multiwfn产生）满足你要求的状态，但动力学过程中没法严格维持。
CP2K支持AIMD过程中用下文介绍的CDFT约束片段的带电和单电子状态，但耗时会额外增加很多
谈谈约束性DFT(CDFT)
http://sobereva.com/271（http://bbs.keinsci.com/thread-535-1-1.html）

通常完全没必要试图做人为约束，电子该怎么分布怎么分布，只要用的泛函的离域性误差别太显著就行，诸如wB97XD通常就足够合理。

作者
Author: -HMDS- 时间: 2023-12-1 18:02

sobereva 发表于 2023-12-1 05:59
没法设。顶多提供一个初猜波函数（可以利用Multiwfn产生）满足你要求的状态，但动力学过程中没法严格维持。 ...

好的，多谢sob，主要问题在于我模拟的几个过程都是一个正离子和一个自由基的反应，我用Multiwfn画了一下电荷和自旋的分布，发现结束状态是对的，但是初始状态是反的，就是电荷基本都在应该是自由基那边，而自旋基本都在应该是正电荷那边，这个正离子在质谱里面是能看到的，所以才有这个问题。
泛函用的是wB97系列的，跑之前做过测试，几个过程参考CCSD(T)还是很准的。

作者
Author: sobereva 时间: 2023-12-2 11:13

-HMDS- 发表于 2023-12-1 18:02
好的，多谢sob，主要问题在于我模拟的几个过程都是一个正离子和一个自由基的反应，我用Multiwfn画了一下 ...

可能说明两个物质接近后会有自发的电子转移

可以做个从远到进的柔性扫描，考察此过程中电子结构状态的变化

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