24小时热门版块排行榜    

查看: 684  |  回复: 1
【悬赏金币】回答本帖问题,作者举重若轻I将赠送您 10 个金币

举重若轻I

新虫 (小有名气)

[求助] 仿真一开始原子间距变大,晶体炸开

请问一下各位大神,我使用MBKS势函数,其中bond使用Morse势来表示共价键程度势函数参数如图所示,可是在NPT系综下晶体会出现很奇怪的现象,原子与原子之间会有很大的空隙,似乎是bond导致的原子之间作用力过大?在NVT系综下,原子相互挤压。
       以下是关于势函数的in文件,请各位大神帮忙看看哪里出现的错误,是势函数写的不对吗?万分感谢!!
units     metal
boundary   p p p
dimension  3
atom_style full
neighbor 2.0  nsq

read_data hfo2-4-4-10.data

group    Hf   type 1               
group    O    type 2        

set type 1 charge   2.40               
set type 2 charge  -1.20

pair_style   hybrid/overlay  buck 8.0   coul/long 10.0

pair_coeff   1 2 buck  12372.16   0.2286   81.36        # hf-o
pair_coeff   2 2 buck  1388.773   0.3623   175.0        # o-o
pair_coeff   * * coul/long

bond_style   morse
bond_coeff   1  0.3474    1.6230   2.0480
kspace_style ewald 1.0e-6

仿真一开始原子间距变大,晶体炸开
势函数参数.png


仿真一开始原子间距变大,晶体炸开-1
初始模型.png


仿真一开始原子间距变大,晶体炸开-2
NPT下原子炸开.png


仿真一开始原子间距变大,晶体炸开-3
势函数表达式,Dij为morse势.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

V11dot

新虫 (初入文坛)

您好,请问可以发一下势函数或者文献的doi号吗,急求!谢谢!!
2楼2022-08-03 03:08:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 举重若轻I 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 啥也不想干 +5 lr小机灵鬼 2024-05-10 5/250 2024-05-10 23:21 by tjmamhtk
[基金申请] 河南工业大学过年开学后,就不让用国基金报账,这合法合规吗? 100+4 知己而已 2024-05-09 22/1100 2024-05-10 22:33 by likeac
[考博] 一篇一作AM,寻求西安985博导 +5 LYWwrz 2024-05-06 6/300 2024-05-10 18:26 by FLIPED
[硕博家园] 29 岁去读一个双非学校机械工程的博士 还有前途吗 +11 funaizheng 2024-05-09 13/650 2024-05-10 17:41 by godlikeme88
[基金申请] 函评这一关首先干掉一大半,今年卷上天了 (金币+5) +26 mycaas 2024-05-07 36/1800 2024-05-10 17:22 by 降温器
[考博] 上海大学2024年秋季博士招生 +5 94872300 2024-05-08 8/400 2024-05-10 17:15 by 甜妹儿~
[基金申请] E05送审了么? +4 wly78528 2024-05-08 4/200 2024-05-10 16:51 by 路遥还有谁
[硕博家园] 5000mg/L的盐分换成百分比等于多少呀 +3 大星Dax 2024-05-06 9/450 2024-05-10 13:30 by zailiyang456
[基金申请] 特聘教授,算高级职称吗?青年基金还没结题的情况下还能报面上吗? +19 aasahr 2024-05-05 28/1400 2024-05-10 13:29 by dxcharlary
[硕博家园] 科研不是打打杀杀,科研就是人情世故 +8 杞天大圣 2024-05-08 8/400 2024-05-10 07:44 by dianda502
[论文投稿] 论文投稿这种情况希望大吗? +3 woshichuxuan 2024-05-08 6/300 2024-05-09 17:27 by topedit
[基金申请] 基金函评会不会跨学部的送啊 +5 cls512 2024-05-09 5/250 2024-05-09 17:08 by cauzcl
[硕博家园] 矿大本-中科大硕(推免)-24年科研助理求职-方向不限 +7 Kelaizhang 2024-05-08 13/650 2024-05-09 16:56 by 安德里斯
[教师之家] 请问有没有想要联合培养学生的老师? +3 吴小懒 2024-05-08 4/200 2024-05-09 15:35 by ou0551
[基金申请] E06送审了吗? +3 edge099 2024-05-08 3/150 2024-05-09 06:32 by Tyler_Ong
[访问学者] NUS Bala教授招CSC联培 +4 山高水远来日方 2024-05-07 5/250 2024-05-08 20:33 by lucas1553
[有机交流] 低沸点物质浓缩 100+3 桃桃PXS 2024-05-04 8/400 2024-05-08 19:41 by 谢大奇
[教师之家] 请教下各位,南京航空航天大学天目湖校区如何? +3 pp爱写论文 2024-05-05 5/250 2024-05-07 16:48 by 似水流年Youth
[基金申请] 基金函评开始了吗 +5 fuzy 2024-05-05 5/250 2024-05-06 11:57 by hdzw9071
[基金申请] 活页超过7000字,有影响吗? +4 青竹82 2024-05-04 5/250 2024-05-04 23:28 by 华氏911.
信息提示
请填处理意见