好的老师,谢谢老师啦 |
Weijie 发表于 2023-3-14 10:18 是 |
sobereva 发表于 2023-3-13 21:03 老师,是我的这边关注的问题,我只关注了博文中总电荷,没有细去核算1-12的电荷之和,这是我的疏忽,非常抱歉。 我确实是做的标准两步RESP拟合,那是做一步的ESP就可以用了嘛 |
官网给出的教程进行限定后所得的原子电荷总和也是所构建基元的所有原子电荷总和符合限定,而不是限定的原子范围内的电荷符合限定条件 没有这回事 RESP拟合静电势电荷的原理以及在Multiwfn中的计算 http://sobereva.com/441(http://bbs.keinsci.com/thread-10880-1-1.html) 里面天冬氨酸的例子里,1-12号原子电荷加和就是用户自己设的-1 只可能是你自己弄错了,或者选择选的是两步式RESP电荷拟合,此时没法用电荷约束。博文里该说的都明确说了 在标准两步式RESP电荷计算时也照样可以设置电荷约束、等价性约束,但只是对第一步起效(默认情况下,第一步的时候既没有用任何等价性约束也没有使用电荷约束)。对于此例,载入上述eqvcons.txt后选择选项0做两步式RESP拟合,会发现最终结果里O5和O6的电荷相同了,但是C7、C11的电荷不同,不同甲基上的氢的电荷也不同。这是因为自定义约束对于第二步不生效,而且按照定义,第二步RESP拟合时甲基氢、甲基碳会被重新拟合。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2024-11-27 19:35 , Processed in 0.204868 second(s), 26 queries , Gzip On.