24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1932  |  回复: 8

453695608

银虫 (正式写手)

[求助] 怎么设置每5步温度升高一度 已有1人参与

怎么样在lammps中设置在5步之后温度升高1度,用的是nvt系宗
回复此楼

» 收录本帖的淘帖专辑推荐

分子动力学

» 猜你喜欢

学习更多的知识,为了未来
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

453695608

银虫 (正式写手)

怎么没有人回答啊
学习更多的知识,为了未来
2楼2015-10-13 21:32:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alundilong

木虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
请参看fix nvt参数temp的说明:
temp values = Tstart Tstop Tdamp
    Tstart,Tstop = external temperature at start/end of run

只要你设置起始,终止温度和合适的Tdamp,那么温度基本就是线性变化。
注意:Tdamp设置要合适,先试一试Tdamp=10个时间步,然后输出温度变化,看看是否合适。
注意:Tdamp不介意太小!
3楼2015-10-14 02:02:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alundilong

木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
453695608: 金币+5, ★★★很有帮助 2015-10-14 09:16:01
引用回帖:
3楼: Originally posted by alundilong at 2015-10-13 12:02:46
请参看fix nvt参数temp的说明:
temp values = Tstart Tstop Tdamp
    Tstart,Tstop = external temperature at start/end of run

只要你设置起始,终止温度和合适的Tdamp,那么温度基本就是线性变化。
注意 ...

注意:Tdamp不建议太小!
4楼2015-10-14 02:04:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

453695608

银虫 (正式写手)

引用回帖:
4楼: Originally posted by alundilong at 2015-10-14 02:04:34
注意:Tdamp不建议太小!...

我想问问,怎么看输出温度啊,要用什么命令
学习更多的知识,为了未来
5楼2015-10-14 09:16:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alundilong

木虫 (正式写手)

【答案】应助回帖

引用回帖:
5楼: Originally posted by 453695608 at 2015-10-13 19:16:03
我想问问,怎么看输出温度啊,要用什么命令...

thermo_style默认输出的温度是基于绝对温度的
compute temp/com你可以用

thermo_style 支持 c_ID, f_ID
对吧?答案显而易见。
6楼2015-10-15 02:17:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

453695608

银虫 (正式写手)

引用回帖:
6楼: Originally posted by alundilong at 2015-10-15 02:17:25
thermo_style默认输出的温度是基于绝对温度的
compute temp/com你可以用

thermo_style 支持 c_ID, f_ID
对吧?答案显而易见。...

fix   3  all nvt temp   200  1200  0.2   
compute 4   all temp/com                                            
thermo_style custom step temp etotal  c_4  f_3                     
dump  2  all  custom  1000  dump.atom    c_1  c_2                  
                                                                    
                                                                    
run                1000000                                                      
                                                                    
                                                                    
Memory usage per processor = 34.0236 Mbytes                        
Step Temp TotEng 4 3                                                
       0            0   -16249.713            0            0        
      50    3.0810635   -16324.981    3.0810635      1612628        
     100    1.0841471     -16363.4    1.0841471      1623777        
     150    8.5668952    -16385.57    8.5668952    1633514.5        
     200  0.068472081   -16396.003  0.068472081      1645359        
     250    161.67359   -16367.088    161.67359      1645578        
     300    46.066746   -16473.096    46.066746    1646098.6        
     350    57.825688   -16509.886    57.825688    1646550.6        
     400    60.454739   -16516.565    60.454739    1646967.1        
     450    67.670996   -16539.283    67.670996    1647399.1        
     500     91.97674   -16534.579     91.97674    1647803.7        
     550    202.90058   -16507.853    202.90058    1648184.4        
     600    446.21263   -16405.259    446.21263    1648504.7        
     650    213.77345   -16518.171    213.77345    1649017.4        
     700    86.530449   -16573.078    86.530449    1649489.2        
     750    108.84435   -16579.516    108.84435    1649902.9        
     800    347.17049   -16483.769    347.17049    1650218.3        
结果中的4   3,是不是就是温度啊,但是感觉3,这个根本不像温度啊
学习更多的知识,为了未来
7楼2015-10-15 20:31:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

alundilong

木虫 (正式写手)

【答案】应助回帖

引用回帖:
7楼: Originally posted by 453695608 at 2015-10-15 06:31:29
fix   3  all nvt temp   200  1200  0.2   
compute 4   all temp/com                                            
thermo_style custom step temp etotal  c_4  f_3                     
dump  2  all   ...

4是温度
3肯定不是,看手册便知(请勤看手册)
8楼2015-10-16 03:37:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

453695608

银虫 (正式写手)

引用回帖:
8楼: Originally posted by alundilong at 2015-10-16 03:37:38
4是温度
3肯定不是,看手册便知(请勤看手册)...

我想问一下为什么温度变化这么不规律啊,nvt不是性线关系吗?
学习更多的知识,为了未来
9楼2015-10-16 15:19:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 453695608 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 311(085601)求调剂 +7 liziyeyeye 2026-03-28 7/350 2026-03-31 19:37 by 星空星月
[考研] 求调剂:085600材料与化工,考材科基,总分319 +16 678lucky 2026-03-31 20/1000 2026-03-31 18:48 by JourneyLucky
[考研] 张芳铭-中国农业大学-环境工程专硕-298 +9 手机用户 2026-03-26 9/450 2026-03-31 18:09 by 544594351
[考研] 311求调剂一志愿合肥工业大学 +11 秋二十二 2026-03-30 11/550 2026-03-31 18:09 by 253863592
[考研] 285求调剂 +6 AZMK 2026-03-30 11/550 2026-03-31 17:33 by AZMK
[考研] 材料求调剂 一志愿哈工大总分298分,前三科223分 +11 dongfang59 2026-03-27 11/550 2026-03-31 16:51 by Wang200018
[考研] 08工科,295,接受跨专业调剂 +4 lmnlzy 2026-03-31 4/200 2026-03-31 15:25 by 西京学院招办
[考研] 367求调剂 +7 芋泥啵啵… 2026-03-28 7/350 2026-03-31 14:55 by 不吃魚的貓
[考研] 一志愿中国科学院大学265求调剂 +6 恬淡ye 2026-03-31 7/350 2026-03-31 11:28 by 恬淡ye
[考研] 085601一志愿西北工业大学初试346 +4 085601初试346 2026-03-30 4/200 2026-03-31 07:47 by jp9609
[考研] 085701环境工程求调剂 +11 多久上课 2026-03-27 12/600 2026-03-30 21:21 by 研究僧导导
[考研] 359求调剂 +5 王了个楠 2026-03-25 5/250 2026-03-30 19:36 by 源_2020
[考研] 289求调剂 +16 新时代材料 2026-03-27 16/800 2026-03-30 19:04 by Wang200018
[考研] 283求调剂(080500) +14 A child 2026-03-27 14/700 2026-03-30 12:06 by 探123
[考研] 296求调剂 +10 彼岸t 2026-03-29 10/500 2026-03-30 10:50 by 探123
[考研] 材料与化工304求B区调剂 +4 邱gl 2026-03-26 7/350 2026-03-30 08:39 by 探123
[考研] 085701环境工程,267求调剂 +16 minht 2026-03-26 16/800 2026-03-28 12:16 by zllcz
[考研] 298调剂 +3 jiyingjie123 2026-03-27 3/150 2026-03-27 11:57 by wxiongid
[考研] 材料求调剂 +5 .m.. 2026-03-25 5/250 2026-03-27 11:08 by 不吃魚的貓
[考研] 化学调剂一志愿上海交通大学336分-本科上海211 +4 小鱼爱有机 2026-03-25 4/200 2026-03-26 10:19 by aa331100
信息提示
请填处理意见