24小时热门版块排行榜    

查看: 1832  |  回复: 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

专家顾问 (著名写手)

伊卡斯特的小猪

引用回帖:
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的回帖
查看全部 7 个回答

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的回帖

15216105346

金虫 (著名写手)

【答案】应助回帖

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

[ 发自小木虫客户端 ]
5楼2015-08-27 18:06:16
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上项目filecode邪修 +4 西山十月 2026-08-09 6/300 2026-08-10 00:16 by liufeng0619
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +7 Lanmanbaby 2026-08-09 14/700 2026-08-09 23:33 by Tide man
[基金申请] 2026国自然放榜时间 +4 布布和一二 2026-08-08 4/200 2026-08-09 22:24 by otj2008
[基金申请] 关于代码变化问题,想知道的进来 +15 且听虎啸 2026-08-07 22/1100 2026-08-09 22:13 by 布布和一二
[基金申请] 这样的filecode谁见过 +10 布布和一二 2026-08-08 21/1050 2026-08-09 17:21 by 天神眷顾
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 5/250 2026-08-08 21:54 by 13586093586
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +7 a089 2026-08-07 7/350 2026-08-08 18:05 by gltch
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 娱乐 +6 Tide man 2026-08-03 6/300 2026-08-07 22:40 by 铁帽子农民
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
信息提示
请填处理意见