24小时热门版块排行榜    

查看: 437  |  回复: 2
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

肖钧江

铜虫 (小有名气)

[求助] 原子温度升不上去

现在遇到一个问题,希望牛人帮忙解答一下
lattice   hcp 3.184
region    box block 0 20 0 10 0 10
create_box 1 box
create_atoms 1 box
首先用上面四条命令创建了8000个hcp原子
1、用fix 1 all npt temp 2000 2000 0.1 iso 1.0 1.0 1,运行出来的dump.atom中的原子不动
但是改为fix 1 all npt temp 2000 2000 0.01 iso 1.0 1.0 0.1运行出来的结果体系完全融化了
2、另外还有一种方式就是先用velocity   all create 1 23435赋予原子1K初始速度
再用fix 1 all npt temp 2000 2000 0.1 iso 1.0 1.0 1命令体系也融化了

有人知道其中的原因么,求指导~~~
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

肖钧江

铜虫 (小有名气)

对了,另外在补充下,log文件里面记录温度是到了2000K,按道理这个温度点下体系早应该融化了
3楼2014-03-25 11:17:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

肖钧江

铜虫 (小有名气)

另外timestep为0.002ps,单位为metal
2楼2014-03-25 10:53:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见