计算化学公社

 找回密码 Forget password
 注册 Register

求助,fdp3-在em过程中出现成键异常导致nvt无法继续进行

查看数: 1869 | 评论数: 10 | 收藏 Add to favorites 0
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2023-3-20 11:14

正文摘要:

问题背景:fdpna3单晶结构如图 这里用到的fdp3-离子是从单晶结构中抠出的;原子电荷和itp,top,gro文件是高斯+multiwfn+sobtop得到的。 1.在em的过程中总是出现成键异常问题,优化结束后总是有如下Energy minimizat ...

回复 Reply

WuQs 发表于 Post on 2023-4-12 21:36:31
已解决
借鉴CHARMM TIP3P修改O3的范德华参数可以解决该问题
sobereva 发表于 Post on 2023-3-29 04:32:31
WuQs 发表于 2023-3-28 20:44
sob老师我有个问题,为什么这种静电吸引会导致H和O完全重合呢,按理说无限接近的时候应该有很大的斥力才 ...

这种情况上面的氢原子对应的原子类型的范德华参数为0(力场的定义)

评分 Rate

参与人数
Participants 1
eV +3 收起 理由
Reason
WuQs + 3 谢谢

查看全部评分 View all ratings

WuQs 发表于 Post on 2023-3-28 20:44:04
sobereva 发表于 2023-3-20 15:17
Gaussian根本没法计算RESP电荷,结合Multiwfn才能计算
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
...

sob老师我有个问题,为什么这种静电吸引会导致H和O完全重合呢,按理说无限接近的时候应该有很大的斥力才对
牧生 发表于 Post on 2023-3-22 12:34:38
本帖最后由 牧生 于 2023-3-22 13:55 编辑
WuQs 发表于 2023-3-22 10:39
但是我有个问题sob老师,如果只是VMD的显示问题,那为什么em.log里会显示在这附近有异常呢;就是我并不需 ...

你试试我这个,已经替换为RESP电荷 f5.chg (1.63 KB, 下载次数 Times of downloads: 1)
f5.itp (28.4 KB, 下载次数 Times of downloads: 1)
f5.mol2 (1.76 KB, 下载次数 Times of downloads: 1)

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
WuQs + 4 虽然但是还是非常感谢

查看全部评分 View all ratings

sobereva 发表于 Post on 2023-3-21 23:58:08
WuQs 发表于 2023-3-21 22:38
在使用# M062X/gen em=GD3 opt SCRF=SMD int=ultrafine 优化结构,结合multiwfn计算出新的原子电荷之后, ...

注意VMD的成键判断规则,看下文
谈谈VMD可视化程序的连接关系的判断和设置问题
http://sobereva.com/534http://bbs.keinsci.com/thread-16396-1-1.html

应当让VMD重新判断成键。从你当前的图来看,原子间距离、角度等都没有异常,重新判断成键后应该看起来很正常。
WuQs 发表于 Post on 2023-3-20 19:43:37
本帖最后由 WuQs 于 2023-3-20 19:45 编辑
sobereva 发表于 2023-3-20 15:17
Gaussian根本没法计算RESP电荷,结合Multiwfn才能计算
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
...

感谢sob老师,我没有说清楚,确实是高斯和multiwfn一起结合算的RESP电荷;也是学着sob老师的博文做的。
自旋度我一直说习惯了,说着说着就忘了本来的名字,我的错,不好意思。
太感谢了,我一直觉得是高斯的命令没用好,但是我没有系统的学习过高斯的使用。
我去试试,如果能解决的话会来更新帖子。
sobereva 发表于 Post on 2023-3-20 15:17:45
Gaussian根本没法计算RESP电荷,结合Multiwfn才能计算
RESP拟合静电势电荷的原理以及在Multiwfn中的计算
http://sobereva.com/441http://bbs.keinsci.com/thread-10880-1-1.html
计算RESP原子电荷的超级懒人脚本(一行命令就算出结果)
http://sobereva.com/476http://bbs.keinsci.com/thread-12858-1-1.html

没有叫"自旋度"的词,那叫自旋多重度

带负电荷这么多的体系得到波函数的时候至少要对带显著负电荷的部分带弥散函数

仔细照着检查http://sobereva.com/soft/Sobtop#FAQ8。里面专门还提到过
比如曾经我用某个程序产生带磺酸基的分子的拓扑文件...

评分 Rate

参与人数
Participants 1
eV +4 收起 理由
Reason
WuQs + 4 谢谢,感动终于有希望了

查看全部评分 View all ratings

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-25 13:54 , Processed in 0.209065 second(s), 26 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list