计算化学公社

 找回密码 Forget password
 注册 Register

模拟氮、氧掺杂石墨烯时,能量极小化报错

查看数: 11055 | 评论数: 36 | 收藏 Add to favorites 3
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2021-12-15 13:58

正文摘要:

本帖最后由 a9471163 于 2021-12-15 14:11 编辑 各位老师,最近有个氮、氧掺杂石墨烯(氧掺杂比例高,达到30%,氮极少,0.6%左右的掺杂比)在gaff力场下用gromacs模拟,拓扑文件是acpype产生的,原子电荷是univer ...

回复 Reply

zhujiabao 发表于 Post on 2025-10-16 10:12:19
a9471163 发表于 2021-12-20 15:22
老师,有重大进展了,我刚才问了同学(因为是跟他合作做这个文章的),他实验那边数据搞错了,本来是氧质 ...

老师您好,想跟您请教一下“他之前一直以为是摩尔比占30%,这大大高估了实验方面所得材料的O含量,事实上这种高O比例的材料可能根本就不存在”。氧化石墨或者氧化石墨烯中,我理解这个氧碳比是正常范围之内的,为什么您会认为这是不存在的呢?我最近也碰到类似的问题了,实在费解
a9471163 发表于 Post on 2022-2-1 23:03:45
喵星大佬 发表于 2021-12-20 07:38
仔细检查你的成键项参数,尤其是涉及环氧的

扭曲到这个程度不可能是静电作用计算的问题,一定是成键参数 ...

老师您好,祝您春节快乐,虎年吉祥,您之前让我认识到oplsaa力场是适合模拟生物炭的,非常感谢您,您之前推断说:扭曲到这个程度不可能是静电作用计算的问题,这个判断总体上也是非常准确的,不过我今天有了新发现。在一个多孔生物炭模拟中(这个模型比之前的更复杂些因为分子的表面缺陷有很多),我发现仅当力场参数合理但原子电荷不合适时,同样模拟报错,并且官能团出现扭曲,所以这个“扭曲到这个程度”也有一部分原因确实是静电作用不合理造成的,新帖在http://bbs.keinsci.com/thread-27631-1-1.html#pid187255,感谢大佬观看和探讨~
a9471163 发表于 Post on 2021-12-31 22:59:29
喵星大佬 发表于 2021-12-22 12:42
opls有1000多个原子类型想必是有环氧相关的

大佬您好,能顺利模拟结束离不开您的帮助,祝您元旦快乐,在数据处理时有1个疑惑,您能否支支招,谢谢您咯~http://bbs.keinsci.com/thread-27222-1-1.html
a9471163 发表于 Post on 2021-12-27 14:56:39
喵星大佬 发表于 2021-12-22 12:42
opls有1000多个原子类型想必是有环氧相关的

大佬,实在是麻烦您了,我的氧化石墨烯在cgenff力场下终于是可以模拟了,多亏了您的帮助,可我还差一个有机小分子(双酚A)的模拟在相同力场下跑不起来,报错归纳到帖子http://bbs.keinsci.com/thread-27141-1-1.html#pid183898里面了,您看能否支支招呢,感谢~
a9471163 发表于 Post on 2021-12-26 23:55:21
本帖最后由 a9471163 于 2021-12-26 23:56 编辑
喵星大佬 发表于 2021-12-22 12:42
opls有1000多个原子类型想必是有环氧相关的

老师您好,想请教一下这种情况您遇到过么,cgenff力场不报错可能量极小化做不了http://bbs.keinsci.com/thread-27132-1-1.html
另外cgenff力场和MMFF力场相比于GAFF而言,参数更多么,能不能很好的描述GO上的环氧基团呢
a9471163 发表于 Post on 2021-12-24 18:20:03
喵星大佬 发表于 2021-12-22 12:42
opls有1000多个原子类型想必是有环氧相关的

大佬您好,我计算RESP2(0.5)作为分子动力学模拟的原子电荷,遇到一个问题即最终的RESP2(0.5)不是整数,在gromacs里加抗衡离子无法中和掉所有电荷,具体描述在帖子http://bbs.keinsci.com/thread-27106-1-1.html#pid183643。想请教一下有解决办法吗,非常感谢
喵星大佬 发表于 Post on 2021-12-22 12:42:31
a9471163 发表于 2021-12-22 10:26
嗯嗯,sob老师刚回复了,说都用resp2就行,我开始算了,这个体系用opls-aa效果很好,配个resp2,初步用op ...

opls有1000多个原子类型想必是有环氧相关的
a9471163 发表于 Post on 2021-12-22 10:26:15
本帖最后由 a9471163 于 2021-12-22 10:28 编辑

嗯嗯,sob老师刚回复了,说都用resp2就行,我开始算了,这个体系用opls-aa效果很好,配个resp2,初步用opls-aa试算(没加电荷的时候),没发现像gaff力场下面的那种环氧基多的区域变形成四面体的问题
喵星大佬 发表于 Post on 2021-12-22 09:55:44
a9471163 发表于 2021-12-22 09:53
谢谢老师,想请教一下对于opls-aa力场,resp2电荷可以跟1.2×CM5在同一体系里混用吗,我的体系有不带电小 ...

当然是不能混用
喵星大佬 发表于 Post on 2021-12-21 17:28:54
a9471163 发表于 2021-12-21 17:13
嗯嗯  我是做吸附的  静电势对静电吸附力影响很大的

那你用resp电荷就好了,你这个体系还不涉及什么包埋原子之类的,也没很多构象,各种ESP电荷应该都差不多
a9471163 发表于 Post on 2021-12-21 17:13:50
喵星大佬 发表于 2021-12-21 16:43
用什么电荷主要看你要模拟的过程和模拟的目的

如果就以结构来说,只看着一个结构的稳定性啊之类的目的 ...

嗯嗯  我是做吸附的  静电势对静电吸附力影响很大的
喵星大佬 发表于 Post on 2021-12-21 16:43:14
a9471163 发表于 2021-12-21 16:08
这是多孔生物炭,布置孔洞模型的。
太详细了大佬,我学到很多您觉得QEq换成RESP,搭配UFF力场 ...

用什么电荷主要看你要模拟的过程和模拟的目的

如果就以结构来说,只看着一个结构的稳定性啊之类的目的,用什么电荷区别倒是不大。这个体系的结构相对刚性,长程静电作用影响有限。所以其实在有些有机力场比如Tripos力场里面直接就没有静电项。但你要是模拟这个结构的吸附过程啊什么的电荷就有影响了,但具体多大不好说
a9471163 发表于 Post on 2021-12-21 16:08:37
本帖最后由 a9471163 于 2021-12-21 16:10 编辑
喵星大佬 发表于 2021-12-21 04:39
看OBGMX的原文就知道了,对于主要近似的就是键角项,总的来说就是用Gromos96形式的键角项代替原版的键角 ...

这是多孔生物炭,布置孔洞模型的。
太详细了大佬,我学到很多您觉得QEq换成RESP,搭配UFF力场以后准确度能提高吗,sob老师说RESP描述静电势比QEq准,但QEq确实是UFF的原配啊
我下一步打算试试oplsaa配RESP2电荷,uff力场配RESP电荷模拟这个体系,还有就是ztop拟合力场参数以后模拟,看看哪种更好

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

GMT+8, 2026-3-12 18:43 , Processed in 0.171688 second(s), 25 queries , Gzip On.

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