24小时热门版块排行榜    

查看: 2537  |  回复: 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的回帖

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

cl1994861

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
gtssongchi: 金币+5, ★★★很有帮助 2017-04-20 21:56:48
弛豫时间是不是设置的不合理,官网mannual里给的一般是1000个时间步
2楼2017-04-20 16:23:07
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +4 vv迷 2026-03-22 4/200 2026-03-22 23:29 by king123!
[考研] 280分求调剂 一志愿085802 +4 PUMPT 2026-03-22 7/350 2026-03-22 22:13 by 星空星月
[考研] 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
[考研] 284求调剂 +5 Zhao anqi 2026-03-22 5/250 2026-03-22 17:38 by barlinike
[考研] 289求调剂 +7 怀瑾握瑜l 2026-03-20 7/350 2026-03-22 15:57 by ColorlessPI
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 材料工程专硕 348分求调剂 +3 冬辞. 2026-03-17 5/250 2026-03-21 18:47 by 学员8dgXkO
[考研] 265求调剂 +12 梁梁校校 2026-03-19 14/700 2026-03-21 13:38 by lature00
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 求调剂 +3 Ma_xt 2026-03-17 3/150 2026-03-21 02:05 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 求调剂 +3 eation27 2026-03-20 3/150 2026-03-20 19:32 by JourneyLucky
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 材料学硕318求调剂 +5 February_Feb 2026-03-19 5/250 2026-03-19 23:51 by 23Postgrad
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 材料,纺织,生物(0856、0710),化学招生啦 +3 Eember. 2026-03-17 9/450 2026-03-18 10:28 by Eember.
信息提示
请填处理意见