计算化学公社

 找回密码 Forget password
 注册 Register

求助:晶胞的原子在模拟过程中有几个原子跑散了

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

正文摘要:

两边是黏土,中间是溶液,为什么晶胞有些原子会自己跑到中间,我用一颗黏土和溶液的时候没有出现这种现象,两颗黏土反而出现晶胞有的原子散了跑向溶液

回复 Reply

wuxiaofei 发表于 Post on 2023-6-9 18:06:12

好的老师,我查了一下文献,发现他们没有把金属与 O H的键角考虑进去,我删除这部分参数就可以了
wuxiaofei 发表于 Post on 2023-6-9 18:04:51
萌新杨同学 发表于 2023-6-9 09:35
感觉很可能是你的拓扑文件有问题,看看每种原子的类型有没有搞错

谢谢,现在我把金属与O H之间的键角参数删除就可以了
萌新杨同学 发表于 Post on 2023-6-9 09:35:02
wuxiaofei 发表于 2023-6-8 15:29
那有什么解决办法吗,我也做了位置限制还是有OH跑出去,还是不用管也没什么影响

感觉很可能是你的拓扑文件有问题,看看每种原子的类型有没有搞错
sobereva 发表于 Post on 2023-6-9 06:19:46
wuxiaofei 发表于 2023-6-8 15:30
老师,我也做了位置限制,还是跑

肯定没设对
wuxiaofei 发表于 Post on 2023-6-8 15:30:13
sobereva 发表于 2023-6-8 13:39
干脆把slab加位置限制固定住就完了,省得排查拓扑文件

老师,我也做了位置限制,还是跑
wuxiaofei 发表于 Post on 2023-6-8 15:29:52
萌新杨同学 发表于 2023-6-8 14:47
可能是周期性边界导致你两个粘土片层间的距离太小,部分离子由于不合理的初始位置就排斥飞走了,毕竟clayff ...

那有什么解决办法吗,我也做了位置限制还是有OH跑出去,还是不用管也没什么影响
萌新杨同学 发表于 Post on 2023-6-8 14:47:53
可能是周期性边界导致你两个粘土片层间的距离太小,部分离子由于不合理的初始位置就排斥飞走了,毕竟clayff不考虑Bond
sobereva 发表于 Post on 2023-6-8 13:39:19
wuxiaofei 发表于 2023-6-8 11:24
老师,这是我的初始构型以及部分itp、mdp文件

干脆把slab加位置限制固定住就完了,省得排查拓扑文件
sobereva 发表于 Post on 2023-6-8 07:22:34
应当把盒子显示出来,否则别人都不知道你模拟用的模型到底是什么样
如果模型没问题就是参数有问题。单个slab模拟时没跑出去可能只是巧合
想绝对避免跑出去可以用谐振势描述slab里的键

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

GMT+8, 2024-11-27 14:32 , Processed in 1.000191 second(s), 26 queries , Gzip On.

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