计算化学公社

标题: oniom报错l701求助 [打印本页]

作者
Author: 〇66 时间: 2025-7-10 22:31
标题: oniom报错l701求助
(, 下载次数 Times of downloads: 2) (, 下载次数 Times of downloads: 2)
麻烦老师可以帮忙看看这个报错是什么原因吗？谢谢老师！

opt=(calcfc,ts,noeigen) oniom(td/wb97xd/def2svp:uff)

作者
Author: wal 时间: 2025-7-10 23:14
感觉像是ONIOM(TDDFT:力场)的时候使用opt=calcfc这方面的问题
试试能不能直接freq 如果可以的话就不要calcfc了改用rcfc读hessian

作者
Author: sobereva 时间: 2025-7-10 23:46
这个计算本身没什么意义。原子电荷都没有指定，完全没法描述分子间静电作用。而且MM层也没冻结，优化过程会严重变形而没法维持晶体下的状态

作者
Author: 〇66 时间: 2025-7-11 09:29

sobereva 发表于 2025-7-10 23:46
这个计算本身没什么意义。原子电荷都没有指定，完全没法描述分子间静电作用。而且MM层也没冻结，优化过程会 ...

谢谢老师答复，刚开始学QM/MM，可能还不太上道，可以麻烦问一下怎么指定原子电荷以及冻结MM层吗？谢谢老师

作者
Author: 〇66 时间: 2025-7-11 09:40

wal 发表于 2025-7-10 23:14
感觉像是ONIOM(TDDFT:力场)的时候使用opt=calcfc这方面的问题
试试能不能直接freq 如果可以的话就不要cal ...

谢谢！我尝试一下# opt=(calcall,ts,noeigen) freq oniom(td/wb97xd/def2svp:uff)

作者
Author: wal 时间: 2025-7-11 09:45

〇66 发表于 2025-7-11 09:40
谢谢！我尝试一下# opt=(calcall,ts,noeigen) freq oniom(td/wb97xd/def2svp:uff)

如果calcfc不行的话calcall大概率也不行 XD
如果确实出问题的话我的建议是测试一下单独freq 如果可以的话用rcfc读

作者
Author: sobereva 时间: 2025-7-11 22:35

〇66 发表于 2025-7-11 09:29
谢谢老师答复，刚开始学QM/MM，可能还不太上道，可以麻烦问一下怎么指定原子电荷以及冻结MM层吗？谢谢老 ...

参考北京科音中级量子化学培训班（http://www.keinsci.com/KBQC）里全面讲Gaussian的力场计算和ONIOM部分的幻灯片：

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

作者
Author: 〇66 时间: 2025-7-12 09:54

sobereva 发表于 2025-7-11 22:35
参考北京科音中级量子化学培训班（http://www.keinsci.com/KBQC）里全面讲Gaussian的力场计算和ONIOM部 ...

谢谢老师答复，可能还有一点小问题想麻烦问一下

作者
Author: wal 时间: 2025-7-12 10:16

〇66 发表于 2025-7-12 09:54
谢谢老师答复，可能还有一点小问题想麻烦问一下

第一列以'-'作为分隔符，字段1为元素符号，字段2为原子类型(记不清名，反正是给力场用的，让力场知道用什么参数)，字段3为原子电荷。你现在只有两个字段，再加个-，后面写电荷即可设置电荷

冻结设置是在gview的atom group里设置的，有个freeze

作者
Author: 〇66 时间: 2025-7-12 10:21

wal 发表于 2025-7-12 10:16
第一列以'-'作为分隔符，字段1为元素符号，字段2为原子类型(记不清名，反正是给力场用的，让力场知道用什 ...

原子冻结的找到啦非常感谢！另外想问一下，指定原子电荷只能一个一个手动改吗？一般每个原子该指定多少电荷是怎么判断呢

作者
Author: wal 时间: 2025-7-12 10:24

〇66 发表于 2025-7-12 10:21
原子冻结的找到啦非常感谢！另外想问一下，指定原子电荷只能一个一个手动改吗？一般每个原子该指定多少 ...

multiwfn有个计算resp电荷的功能，用那个
转移到gjf里可以直接插入列，嫌麻烦可以考虑我的banetop

作者
Author: 〇66 时间: 2025-7-12 10:30

wal 发表于 2025-7-12 10:24
multiwfn有个计算resp电荷的功能，用那个
转移到gjf里可以直接插入列，嫌麻烦可以考虑我的banetop

好的好的谢谢你太好啦呜呜呜

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