24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1931  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华东理工大学,080500学硕,317分,求调剂 +8 s1145 2026-03-31 8/400 2026-03-31 17:37 by 唐沐儿
[考研] 309求调剂 +16 谁不是少年 2026-03-29 16/800 2026-03-31 16:47 by 不吃魚的貓
[考研] 343求调剂 +8 爱羁绊 2026-03-28 8/400 2026-03-31 16:12 by 不吃魚的貓
[基金申请] 面上5B能上会吗? +8 redcom 2026-03-29 8/400 2026-03-31 15:53 by niuailing
[考研] 346求调剂 一志愿070303有机化学 +10 萝卜炖青菜 2026-03-28 11/550 2026-03-31 15:19 by 六翼TMR
[考研] 本科211生物医学工程085409求调剂339分 +7 里子木yy 2026-03-29 7/350 2026-03-31 14:35 by fmesaito
[考研] 一志愿211,335分,0856,求调剂院校和导师 +10 倾____萧 2026-03-27 11/550 2026-03-31 14:32 by fmesaito
[考研] 286求调剂 +5 丢掉懒惰 2026-03-27 8/400 2026-03-31 11:27 by Delta2012
[考研] 总分322求生物学/生化与分子/生物信息学相关调剂 +6 星沉uu 2026-03-26 7/350 2026-03-31 10:19 by GdShizy
[考研] 293分求调剂,外语为俄语 +5 加一一九 2026-03-31 5/250 2026-03-31 09:39 by zhshch
[考研] 370求调剂 +3 080700调剂 2026-03-30 3/150 2026-03-31 01:09 by A_Zhe
[有机交流] 甲基亚磺磺酸钠和甲基磺酸酯反应机理 10+3 kaobao456 2026-03-29 4/200 2026-03-30 23:16 by nBu锂
[考研] 调剂 +4 GK72 2026-03-30 4/200 2026-03-30 20:32 by dick_runner
[考研] 各位老师好,我的一志愿为北京科技大学085601材料专硕 +10 Koxui 2026-03-28 10/500 2026-03-30 19:33 by 源_2020
[考研] 292求调剂 +13 是妍子也是研子 2026-03-30 13/650 2026-03-30 18:01 by 小徐0109
[考研] 317分 一志愿南理工材料工程 本科湖工大 求调剂 +12 芋泥小铃铛 2026-03-28 12/600 2026-03-30 17:06 by wangjy2002
[考研] 349求调剂 +6 李木子啊哈哈 2026-03-25 6/300 2026-03-29 12:47 by 无际的草原
[考研] 356求调剂 +4 gysy?s?a 2026-03-28 4/200 2026-03-29 10:32 by 唐沐儿
[考研] 求调剂 +6 芦lty 2026-03-25 7/350 2026-03-28 13:13 by 唐沐儿
[考研] 网络空间安全0839招调剂 +4 w320357296 2026-03-25 6/300 2026-03-25 17:59 by 255671
信息提示
请填处理意见