24小时热门版块排行榜    

查看: 2529  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 328求调剂,英语六级551,有科研经历 +6 生物工程调剂 2026-03-17 10/500 2026-03-22 20:22 by edmund7
[考研] 298求调剂 +6 上岸6666@ 2026-03-20 6/300 2026-03-22 20:21 by edmund7
[考研] 一志愿武理材料工程348求调剂 +5  ̄^ ̄゜汗 2026-03-19 7/350 2026-03-22 19:44 by 公瑾逍遥
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-20 3/150 2026-03-22 16:00 by ColorlessPI
[考研] 298求调剂一志愿211 +3 上岸6666@ 2026-03-20 3/150 2026-03-22 15:50 by ColorlessPI
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +3 niko- 2026-03-22 3/150 2026-03-22 15:15 by 杨杨杨紫
[考研] 070300化学求调剂 +3 苑豆豆 2026-03-20 3/150 2026-03-22 14:48 by ColorlessPI
[考研] 考研调剂 +4 来好运来来来 2026-03-21 4/200 2026-03-22 12:15 by 星空星月
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 265求调剂 +9 梁梁校校 2026-03-17 9/450 2026-03-21 02:17 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 296求调剂 +6 www_q 2026-03-18 10/500 2026-03-20 23:56 by JourneyLucky
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 290求调剂 +7 ^O^乜 2026-03-19 7/350 2026-03-20 21:43 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 0703化学调剂 +5 pupcoco 2026-03-17 8/400 2026-03-19 13:58 by houyaoxu
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 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
信息提示
请填处理意见