计算化学公社

标题: gaussian+pysoc可以计算固态下的旋轨耦合吗 [打印本页]

作者
Author: zj123_45 时间: 2020-12-18 19:41
标题: gaussian+pysoc可以计算固态下的旋轨耦合吗
老师好，我想计算一下oniom模型下的单三重态之间的自旋轨道耦合，基于优化后的s0结构输入文件： (, 下载次数 Times of downloads: 20) 得到.log和.rwf文件，
但是pysoc计算后出错：
(, 下载次数 Times of downloads: 11)
Gaussian+pysoc可以用来计算oniom模型下的旋轨耦合吗，orca我们组的服务器太老没法用...所以考虑用pysoc
这是.py文件 (, 下载次数 Times of downloads: 1)

作者
Author: 让你变成回忆 时间: 2020-12-19 09:43
估计ONIOM的输出pysoc是没法利用的。
可以考虑用charge替代MM部分，这样MM与QM的静电相互作用会考虑进去，实际上，对于算单点的目的，MM部分的范德华力不会对QM部分产生影响，所以利用charge背景电荷来替代MM是可以的。

作者
Author: sobereva 时间: 2020-12-19 09:46
估计不能
ONIOM的输出信息和一般情况有极为显著的不同，诸如chk里都没轨道波函数
就算用簇模型算固态，也完全没必要用ONIOM，把周围区域原子用背景电荷表示就完事了。
而且你用QEq这种对静电势描述巨垃圾的电荷表现MM区域对QM区域的静电作用是完全不能接受的。另外，你都没写embedcharge，都没用电子嵌入，相当于周围部分对QM部分波函数根本都没起到极化效果。
你现在等于完全是瞎算。

作者
Author: 让你变成回忆 时间: 2020-12-19 09:55
本帖最后由让你变成回忆于 2020-12-19 10:04 编辑

sobereva 发表于 2020-12-19 09:46
估计不能
ONIOM的输出信息和一般情况有极为显著的不同，诸如chk里都没轨道波函数
就算用簇模型算固态，也 ...

ONIOM在国内都快被玩坏了。
我之前见过有的输入文件，关键词部分没有UFF=QEq，在元素定义部分也没有指定原子电荷。
感觉sob老师有必要专门写一个关于ONIOM的博文。

作者
Author: yygong 时间: 2020-12-20 21:50

让你变成回忆发表于 2020-12-19 09:55
ONIOM在国内都快被玩坏了。
我之前见过有的输入文件，关键词部分没有UFF=QEq，在元素定义部分也没有指定 ...

确实如此，甚至连输入文件都不愿（敢）给人看。

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