24小时热门版块排行榜    

查看: 2527  |  回复: 5

gtssongchi

木虫 (小有名气)

[求助] fix nvt 温度设置为300K,可是结果却显示只有150K,是什么原因呢? 已有2人参与

大家好,lammps初学,请问我在设置中出了什么问题,导致最后热力学输出的温度只有150K
#lammps in file for melt of Cu

echo                        screen

units                        metal
boundary                p p p
atom_style                atomic

lattice                        fcc        3.61        #建立一个晶格参数为3.61埃的fcc晶体
region                        box        block        0 10 0 10 0 10       
create_box                1        box       
create_atoms                1        box       

pair_style                eam
pair_coeff                * * ./potentials/Cu_u3.eam        #格式确保正确


timestep                0.003        #单位是ps
neighbor                2.0        bin
neigh_modify                delay 0 every 5 check yes

velocity                all        create        300  474747 dist gaussian

#variable                N        equal        abs{floor(x[100]+10000)}
#velocity                all        create        N        dist gaussian

#fix                        1        all  nvt  temp  300  300  100.0        # 100是耦合系数,给出驰豫时间100×时间单位
fix                        1         all        nve
fix                        2        all  temp/berendsen 300.0 300.0 100.0


thermo                        100       
#thermo_style                默认
#dump                        1        all        custom  100  dump.lammpstrj  id  type  x  y  z       

run                        10000
回复此楼

» 猜你喜欢

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

cl1994861

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
gtssongchi: 金币+5, ★★★很有帮助 2017-04-20 21:56:48
弛豫时间是不是设置的不合理,官网mannual里给的一般是1000个时间步
2楼2017-04-20 16:23:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gtssongchi

木虫 (小有名气)

引用回帖:
2楼: Originally posted by cl1994861 at 2017-04-20 16:23:07
弛豫时间是不是设置的不合理,官网mannual里给的一般是1000个时间步

我发现问题原因可能是跑得时间步不够,这里之设置了10000步,将时间步改为50000步后,后期温度确实稳定在300K了,只是我很奇怪,为什么lammps刚开始的温度会是150K呢,它从150K慢慢上升到300K,然后稳定。
3楼2017-04-20 17:19:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
gtssongchi: 金币+15, ★★★★★最佳答案 2017-04-20 21:56:31
引用回帖:
2楼: Originally posted by cl1994861 at 2017-04-20 16:23:07
弛豫时间是不是设置的不合理,官网mannual里给的一般是1000个时间步

No, 100 here is too large in the "metal" unit system, where the unit of time is ps. Change it to 0.1.
Also, 10000 steps are too short to reach equilibrium.
4楼2017-04-20 17:22:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gtssongchi

木虫 (小有名气)

引用回帖:
4楼: Originally posted by brucefan at 2017-04-20 17:22:06
No, 100 here is too large in the "metal" unit system, where the unit of time is ps. Change it to 0.1.
Also, 10000 steps are too short to reach equilibrium....

感谢,可能是我理解有误,换成0.1之后,的确解决了这个问题。
5楼2017-04-20 21:12:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hookah

木虫 (正式写手)

除了4楼的提示,你可以考虑去掉:fix                        1         all        nve
你先做的NVT,然后又设置nve系综,NVE系综下T肯定变化 温度当然不是300K
6楼2017-05-07 21:20:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 gtssongchi 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-20 3/150 2026-03-22 16:00 by ColorlessPI
[考研] 招08考数学 +3 laoshidan 2026-03-20 7/350 2026-03-22 12:54 by 在风落中
[考研] 303求调剂 +5 安忆灵 2026-03-22 6/300 2026-03-22 12:46 by 素颜倾城1988
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +10 枫桥ZL 2026-03-18 12/600 2026-03-21 22:02 by peike
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 306求0703调剂一志愿华中师范 +5 纸鱼ly 2026-03-21 5/250 2026-03-21 17:11 by 学员8dgXkO
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 初始318分求调剂(有工作经验) +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 材料与化工 322求调剂 +4 然11 2026-03-19 4/200 2026-03-20 22:12 by luoyongfeng
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 一志愿西南交通 专硕 材料355 本科双非 求调剂 +5 西南交通专材355 2026-03-19 5/250 2026-03-20 21:10 by JourneyLucky
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 一志愿西安交通大学 学硕 354求调剂211或者双一流 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[论文投稿] 有没有大佬发小论文能带我个二作 +3 增锐漏人 2026-03-17 4/200 2026-03-17 09:26 by xs74101122
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
信息提示
请填处理意见