24小时热门版块排行榜    

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

木虫 (小有名气)

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

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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 323求调剂 +6 洼小桶 2026-03-18 6/300 2026-03-23 00:29 by king123!
[考研] 328求调剂,英语六级551,有科研经历 +6 生物工程调剂 2026-03-17 10/500 2026-03-22 20:22 by edmund7
[考研] 287求调剂 +8 晨昏线与星海 2026-03-19 9/450 2026-03-22 17:01 by i_cooler
[考研] 293求调剂 +12 zjl的号 2026-03-16 17/850 2026-03-22 16:51 by i_cooler
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
[考研] 085600材料与化工306 +4 z1z2z3879 2026-03-21 4/200 2026-03-21 23:44 by ms629
[考研] 一志愿南大,0703化学,分数336,求调剂 +3 收到VS 2026-03-21 3/150 2026-03-21 18:42 by 学员8dgXkO
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 299求调剂 +4 某某某某位 2026-03-21 4/200 2026-03-21 16:30 by barlinike
[考研] 求调剂 +3 .m.. 2026-03-21 4/200 2026-03-21 16:25 by barlinike
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 一志愿西南交大,求调剂 +5 材化逐梦人 2026-03-18 5/250 2026-03-21 00:26 by JourneyLucky
[考研] 274求调剂 +10 S.H1 2026-03-18 10/500 2026-03-20 23:51 by JourneyLucky
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 085601材料工程专硕求调剂 +10 慕寒mio 2026-03-16 10/500 2026-03-19 15:26 by 丁丁*
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
信息提示
请填处理意见