计算化学公社

标题: 求助：请问含有羟基自由基的体系做AIMD模拟的计算级别选用 [打印本页]

作者
Author: Lisy 时间: 2025-9-6 13:58
标题: 求助：请问含有羟基自由基的体系做AIMD模拟的计算级别选用
请问含有8个羟基自由基的水溶液体系用CP2K做AIMD模拟，计算级别应该用什么比较好呢？用PBE/DZVP-MOLOPT-SR-GTH可以吗？自旋多重度是多少呢？

作者
Author: sobereva 时间: 2025-9-6 16:24
带上D3校正后可以接受

自旋多重度1，并且开UKS，且通过magnetization对初猜波函数里的各个OH的O设置初始自旋磁矩，四个+1四个-1

作者
Author: Lisy 时间: 2025-9-6 19:37

sobereva 发表于 2025-9-6 16:24
带上D3校正后可以接受

自旋多重度1，并且开UKS，且通过magnetization对初猜波函数里的各个OH的O设置初始 ...

谢谢卢老师，请问是不是选择4个OH中对应O的原子编号，然后将magnetization设置为4，然后输入另外四个原子变化，然后输入-4呀

作者
Author: sobereva 时间: 2025-9-7 14:57

Lisy 发表于 2025-9-6 19:37
谢谢卢老师，请问是不是选择4个OH中对应O的原子编号，然后将magnetization设置为4，然后输入另外四个原子 ...

“选择4个OH中对应O的原子编号”语义不明。把前提描述充分

作者
Author: Lisy 时间: 2025-9-8 11:41

sobereva 发表于 2025-9-7 14:57
“选择4个OH中对应O的原子编号”语义不明。把前提描述充分

就是假如我这八个O的编号是1-8，那我可不可以把1-4归为一类O，总磁矩设置为4，另外5-8归为另外一类O，总磁矩设置为-4

作者
Author: sobereva 时间: 2025-9-8 12:42

Lisy 发表于 2025-9-8 11:41
就是假如我这八个O的编号是1-8，那我可不可以把1-4归为一类O，总磁矩设置为4，另外5-8归为另外一类O，总 ...

1-4用一个原子类型，初始磁矩都设+1，5-8用另一个原子类型，初始磁矩都用-1

根本没法对某一类原子类型里的所有原子设总磁矩（除非你要用CDFT）

作者
Author: Lisy 时间: 2025-9-8 14:49

sobereva 发表于 2025-9-8 12:42
1-4用一个原子类型，初始磁矩都设+1，5-8用另一个原子类型，初始磁矩都用-1

根本没法对某一类原子类型 ...

明白了，谢谢卢老师！

作者
Author: ljh123 时间: 2025-9-16 14:48

sobereva 发表于 2025-9-8 12:42
1-4用一个原子类型，初始磁矩都设+1，5-8用另一个原子类型，初始磁矩都用-1

根本没法对某一类原子类型 ...

老师请教一下，这里8个羟基自由基按照自旋极化单重态来计算是什么原因呢。是因为溶液体系下，这样的能量是最低的吗

作者
Author: sobereva 时间: 2025-9-16 20:52

ljh123 发表于 2025-9-16 14:48
老师请教一下，这里8个羟基自由基按照自旋极化单重态来计算是什么原因呢。是因为溶液体系下，这样的能量 ...

想想为什么诸如丁烷双自由基基态是自旋极化单重态。显然因为高自旋状态能量更高
虽然单电子位点充分分离的情况两种状态能量是简并的，但保不齐动力学过程中出现单电子位点恰好距离较近的情况

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