24小时热门版块排行榜    

查看: 1828  |  回复: 6

lkd-syl

金虫 (著名写手)

[求助] 动力学跑出这样的结构,是什么原因啊? 已有1人参与

RT,跑铁的超胞动力学出来一个这样的结果,感觉是错了,请各位大神分析一下是什么原因,不胜感激,,,
   输入参数:
md conp conse qok nomod pres spat pot efg
ensemble npt
tau_thermostat  0.1 ps
tau_barostat  0.1 ps

temperature 1500.000000 K
pressure 0.000101 GPa

equilibration 10 ps
timestep 2 fs
production 50 ps
write 100
sample 100


title
GULP calculation from Materials Studio for Fe9826
end


cell
48.728800 48.728800 48.728800 90.000000 90.000000 90.000000

fractional
Fe1 core 0.000000 0.000000 0.000000 0.000000 1.000000 0.0 1 1 1
Fe1 core 0.029412 0.029412 0.029412 0.000000 1.000000 0.0 1 1 1

输入的结果以几张附图所示,,,

动力学跑出这样的结构,是什么原因啊?
结构.jpg


动力学跑出这样的结构,是什么原因啊?-1
constant of motion.jpg


动力学跑出这样的结构,是什么原因啊?-2
pressure.jpg


动力学跑出这样的结构,是什么原因啊?-3
temperature.jpg
回复此楼
Nothing impossible, anything impossible!!
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

pigrass

专家顾问 (著名写手)

伊卡斯特的小猪

50ps后npt系综的温度还有这么大的波动,说明结构还不稳定
如果是觉得原子在盒子“外”比较奇怪,右键display style里改in-cell
#有问题发论坛,不接受私聊#爱晶体,爱MaterialsStudio#MS专辑点击右边#http://emuch.net/taotie/41673.html#物性数据#http://webbook.nist.gov/chemistry/#
2楼2015-08-27 08:39:06
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lkd-syl

金虫 (著名写手)

引用回帖:
2楼: Originally posted by pigrass at 2015-08-26 12:39:06
50ps后npt系综的温度还有这么大的波动,说明结构还不稳定
如果是觉得原子在盒子“外”比较奇怪,右键display style里改in-cell

结构不稳定是因为在这个原子数下50ps的时间不够吗?  另外原来“盒子”内的原子跑完动力学后成为这样,感觉像是整体的移动,这是不是说明已经有比较大的问题了呀?
Nothing impossible, anything impossible!!
3楼2015-08-27 09:47:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

pigrass

专家顾问 (著名写手)

伊卡斯特的小猪

引用回帖:
3楼: Originally posted by lkd-syl at 2015-08-27 09:47:02
结构不稳定是因为在这个原子数下50ps的时间不够吗?  另外原来“盒子”内的原子跑完动力学后成为这样,感觉像是整体的移动,这是不是说明已经有比较大的问题了呀?...

周期性结构没有整体移动的说法,你以为原子“移动到外面"了,可是你把display style设置成incell,你还能看得出移动吗,可能只是显示问题罢了。除非你看到原先松散排列的都团聚到一起,那才是真的移动

平衡的方法建议找找文献,看看别人用的什么势函数,用的什么系综,另外,原子多的话平衡时间肯定长
#有问题发论坛,不接受私聊#爱晶体,爱MaterialsStudio#MS专辑点击右边#http://emuch.net/taotie/41673.html#物性数据#http://webbook.nist.gov/chemistry/#
4楼2015-08-27 11:12:34
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

15216105346

金虫 (著名写手)

【答案】应助回帖

★ ★ ★
感谢参与,应助指数 +1
lkd-syl: 金币+3, 有帮助 2015-08-27 19:29:33
个人感觉,你的势参数有问题!要是用文献的参数注意单位换算问题!我曾经也有过这样的经历

[ 发自小木虫客户端 ]
5楼2015-08-27 18:06:16
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lkd-syl

金虫 (著名写手)

引用回帖:
5楼: Originally posted by 15216105346 at 2015-08-26 22:06:16
个人感觉,你的势参数有问题!要是用文献的参数注意单位换算问题!我曾经也有过这样的经历

请问你的后边是怎么解决的,这个也是查阅文献所引用的,比较粗糙,具体一点就需要拟合新的势函数,,,
Nothing impossible, anything impossible!!
6楼2015-08-27 19:29:16
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

lkd-syl

金虫 (著名写手)

引用回帖:
4楼: Originally posted by pigrass at 2015-08-26 15:12:34
周期性结构没有整体移动的说法,你以为原子“移动到外面"了,可是你把display style设置成incell,你还能看得出移动吗,可能只是显示问题罢了。除非你看到原先松散排列的都团聚到一起,那才是真的移动

平衡 ...

恩,原来我也是感觉看着不舒服而已,我再调整一下试试,,,,
Nothing impossible, anything impossible!!
7楼2015-08-27 19:30:38
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ONGIONG 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 这样的filecode谁见过 +3 布布和一二 2026-08-08 5/250 2026-08-08 12:22 by 仁砚薪传
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +4 a089 2026-08-07 4/200 2026-08-08 12:15 by arzu_hma
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] fileCode有新解读? +8 Tide man 2026-08-08 16/800 2026-08-08 09:59 by Tide man
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 08:59 by 6vVgjDL4CnGu
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 HEQlVqMTIA7d 2026-08-07 4/200 2026-08-08 08:27 by 6vVgjDL4CnGu
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 KXLV3nuBVBY7 2026-08-07 4/200 2026-08-08 08:21 by 6vVgjDL4CnGu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 关于filecode +4 布布和一二 2026-08-07 7/350 2026-08-07 22:55 by zhanghaozhu
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode与中标关系的预测 +3 布布和一二 2026-08-07 3/150 2026-08-07 15:09 by gltch
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] 关于豆爷回答的JTJC与%2F数量 +4 yang182083 2026-08-06 6/300 2026-08-07 12:43 by 虫友是什么虫
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
信息提示
请填处理意见