计算化学公社

标题: 求助晶体环境下内重组能的计算 [打印本页]

作者
Author: jjm 时间: 2022-4-25 21:47
标题: 求助晶体环境下内重组能的计算
各位老师好：
我想要计算晶体环境下内重组能，用oniom实现。
第一步：优化中性状态的结构，给MM区设定了chelpg电荷
看sob老师的帖子，第二步是将中心分子设成带电状态，对a的结构做单点。
请问大家第二步，是否是将中心分子带电状态下的的chelpg电荷拟合出来，同时oniom中的电荷多重度设置为：0 1 1 2 1 2 呢？
谢谢大家

作者
Author: sobereva 时间: 2022-4-26 07:34
以ONIOM(QM:MM)方式计算这种时，一般假定周围环境分子的原子电荷不受中心分子带电状态影响，否则考虑起来很麻烦（而且确切来说还需要考虑迭代来让中心分子和环境分子的电荷分布达到自洽，见http://bbs.keinsci.com/thread-20960-1-1.html），因此第二步不需要重算环境分子的CHELPG电荷

作者
Author: jjm 时间: 2022-4-26 09:03

sobereva 发表于 2022-4-26 07:34
以ONIOM(QM:MM)方式计算这种时，一般假定周围环境分子的原子电荷不受中心分子带电状态影响，否则考虑起来很 ...

非常感谢老师的回复。
老师，不知道我理解的对不对：计算第二步时，MM区的原子电荷就用计算第一步的原子电荷，但是QM区的中心分子需要添加带电状态的 CHELPG 电荷。oniom的的电荷自旋多重度设置为 0 1 1 2 1 2
谢谢老师

作者
Author: sobereva 时间: 2022-4-27 05:17

jjm 发表于 2022-4-26 09:03
非常感谢老师的回复。
老师，不知道我理解的对不对：计算第二步时，MM区的原子电荷就用计算第一步的原子 ...

对

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