应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 如何实现在NVT或者NVE系综下升降温？
81/1返回列表
查看: 2861  |  回复: 7
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

2OLEGEND

新虫 (正式写手)

[交流] 如何实现在NVT或者NVE系综下升降温？

1，在NVT系综下，如何对体系进行降温，并控制降温速率呢？我先前在NVT， T = 3000K条件下计算了500000步，现在需要考虑降温度将为300K，并控制降温速率（100K/ps，timestep为0.1fs），运行270000步，在restart文件中该如何写呢？

fix 4 all nvt 3000K 300K 27000
run 270000

这样写貌似并不能控制降温速率，该结合哪些命令呢，希望各位大神告知！

2，同样的在NVE系综下，跟上述条件一样

fix 4 all nve
fix 5 all temp/berendsen 3000 300 27000
run 270000

这样写是不是只能在规定的时间内降温到目标温度，而没有控制降温速率？该如何控制呢？
回复此楼

» 本帖已获得的红花（最新10朵）

brucefan

» 猜你喜欢
1楼 2018-03-14 15:58:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
2OLEGEND: 金币+10 2018-03-23 11:33:12
你需要的是这个:
fix xxx all temp/rescale 1 3000.0 300.0 0.0 1.0

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

2OLEGEND
2楼2018-03-15 00:18:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2OLEGEND

新虫 (正式写手)
送红花一朵
引用回帖:
2楼: Originally posted by brucefan at 2018-03-15 00:18:42
你需要的是这个:
fix xxx all temp/rescale 1 3000.0 300.0 0.0 1.0

谢谢老师的回答！
再追加一句，使用fix xxx all temp/rescale 1 3000.0 300.0 0.0 1.0命令之后，使用run命令来控制速率。
run 270000，也就是100K/ps的速率降温，是这个意思吗？
赞一下 回复此楼
3楼2018-03-15 10:33:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)
送红花一朵
是的。

[ 发自手机版 http://muchong.com/3g ]
回复此楼
4楼2018-03-15 13:28:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
另外，你还是额外需要用fix nve提供积分操作，因为fix temp/rescale不积分。请仔细阅读这两个命令的手册。

忠告:至少仔细阅读一遍你的脚本中所用到的每一个命令的相关手册部分。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

2OLEGEND
5楼2018-03-15 13:34:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2OLEGEND

新虫 (正式写手)
送红花一朵
引用回帖:
5楼: Originally posted by brucefan at 2018-03-15 13:34:28
另外，你还是额外需要用fix nve提供积分操作，因为fix temp/rescale不积分。请仔细阅读这两个命令的手册。

忠告:至少仔细阅读一遍你的脚本中所用到的每一个命令的相关手册部分。
...

嗯嗯，谢谢老师的回答。我看了一遍命令的含义，有些理解的不够深，所以会有疑问。
手册上说fix temp/rescale不能和fix nvt混合使用，只能和fix nve搭配使用。也就是说我先前的in文件中写了fix nvt，
再读restart文件时，需要将其改为fix nve+fix temp/rescale才能实现温度按速率降温？
谢谢您的回答！
赞一下 回复此楼
6楼2018-03-15 18:10:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

brucefan

专家顾问 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
如果之前有fix nvt ，只要 unfix它，再写后面的组合命令fix nve+fix temp/rescale即可。

[ 发自手机版 http://muchong.com/3g ]
赞一下 回复此楼
7楼2018-03-18 18:49:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gengxin60
8楼2018-03-18 23:25:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 2OLEGEND 的主题更新
81/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:42 by XiYP04It3Evb
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 3xN9SualEUPt 2026-10-24 3/150 2026-10-25 14:23 by XiYP04It3Evb
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 14:07 by XiYP04It3Evb
[论文投稿] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 3xN9SualEUPt 2026-10-24 4/200 2026-10-25 13:54 by XiYP04It3Evb
[考博] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-25 13:28 by XiYP04It3Evb
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 DI1qFHFK0KbE 2026-10-24 4/200 2026-10-25 13:18 by XiYP04It3Evb
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 7xHjyB3bAiwg 2026-10-24 4/200 2026-10-25 13:15 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +6 7xHjyB3bAiwg 2026-10-24 6/300 2026-10-25 13:11 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急 +5 7xHjyB3bAiwg 2026-10-23 10/500 2026-10-25 12:51 by XiYP04It3Evb
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-25 12:49 by XiYP04It3Evb
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 11/550 2026-10-25 12:45 by XiYP04It3Evb
[找工作] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 7xHjyB3bAiwg 2026-10-23 4/200 2026-10-25 12:36 by XiYP04It3Evb
[教师之家] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 GVGRDgNfv7P4 2026-10-23 4/200 2026-10-25 12:10 by XiYP04It3Evb
[公派出国] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 GVGRDgNfv7P4 2026-10-23 9/450 2026-10-25 12:10 by XiYP04It3Evb
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 7xHjyB3bAiwg 2026-10-23 8/400 2026-10-24 22:56 by XcOYG53d3x4p
[教师之家] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 7xHjyB3bAiwg 2026-10-24 3/150 2026-10-24 15:41 by Equinoxhua
[文学芳草园] 说走就走，走到哪便是哪 +5 Amei_123 2026-10-18 10/500 2026-10-23 22:38 by Amei_123
[文学芳草园] 亲密关系 +7 宋昙 2026-10-21 13/650 2026-10-23 17:07 by 小牧11
[考博] 27秋申博自荐 +4 且听风吟时 2026-10-21 4/200 2026-10-23 13:26 by 大大大肥虫
[考博] 郑州大学申博 +5 stars~ 2026-10-19 8/400 2026-10-22 15:06 by stars~
信息提示
关闭
请填处理意见
关闭