wal 发表于 2025-7-12 10:24 好的好的 谢谢 你太好啦 呜呜呜 ![]() |
〇66 发表于 2025-7-12 10:21 multiwfn有个计算resp电荷的功能,用那个 转移到gjf里可以直接插入列,嫌麻烦可以考虑我的banetop |
参与人数Participants 1 | eV +5 | 收起 理由Reason |
---|---|---|
| + 5 | 谢谢 |
〇66 发表于 2025-7-12 09:54 第一列以'-'作为分隔符,字段1为元素符号,字段2为原子类型(记不清名,反正是给力场用的,让力场知道用什么参数),字段3为原子电荷。你现在只有两个字段,再加个-,后面写电荷即可设置电荷 冻结设置是在gview的atom group里设置的,有个freeze |
〇66 发表于 2025-7-11 09:29 参考北京科音中级量子化学培训班(http://www.keinsci.com/KBQC)里全面讲Gaussian的力场计算和ONIOM部分的幻灯片: ![]() |
〇66 发表于 2025-7-11 09:40 如果calcfc不行的话calcall大概率也不行 XD 如果确实出问题的话我的建议是测试一下单独freq 如果可以的话用rcfc读 |
参与人数Participants 1 | eV +5 | 收起 理由Reason |
---|---|---|
| + 5 |
wal 发表于 2025-7-10 23:14 谢谢!我尝试一下# opt=(calcall,ts,noeigen) freq oniom(td/wb97xd/def2svp:uff) |
sobereva 发表于 2025-7-10 23:46 谢谢老师答复,刚开始学QM/MM,可能还不太上道,可以麻烦问一下怎么指定原子电荷以及冻结MM层吗?谢谢老师 |
这个计算本身没什么意义。原子电荷都没有指定,完全没法描述分子间静电作用。而且MM层也没冻结,优化过程会严重变形而没法维持晶体下的状态 |
感觉像是ONIOM(TDDFT:力场)的时候使用opt=calcfc这方面的问题 试试能不能直接freq 如果可以的话就不要calcfc了 改用rcfc读hessian |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-16 18:34 , Processed in 0.194794 second(s), 26 queries , Gzip On.