24小时热门版块排行榜    

查看: 1733  |  回复: 7

skplin

铜虫 (小有名气)

[求助] 分子动力学我设置的恒温为什么降温了 已有1人参与

ISTART=0
ICHARG=2
NPAR=4

PREC=Accurate
ENCUT=300
LREAL=Auto
ISMEAR=-1
SIGMA=0.05
ISPIN=2

IBRION=-1
ISIF=2
NSW=0
POTIM=0.5

APACO=16
NPACO=256

TEBEG=1000
TEEND=1000
SMASS=1

LWAVE=.FALSE.
LCHARG=.TRUE.
哪里错了?????

分子动力学我设置的恒温为什么降温了
$D{68I7D_7G8302HUJ4HKH6.png
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
skplin: 金币+10 2016-05-18 19:52:14
如果体系是优化过的,温度大概会降到你设置的初始值的一半然后慢慢回升到你的heatbath的温度.
看看virial theorem就明白了.
2楼2016-05-10 18:09:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)

【答案】应助回帖

还有IBRION是不是-1吧,你仔细看看manual
3楼2016-05-10 18:10:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skplin

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by luzihen at 2016-05-10 18:09:50
如果体系是优化过的,温度大概会降到你设置的初始值的一半然后慢慢回升到你的heatbath的温度.
看看virial theorem就明白了.

该怎么办
smass如果等于2会有好转吗?
4楼2016-05-18 19:52:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)

引用回帖:
4楼: Originally posted by skplin at 2016-05-18 19:52:51
该怎么办
smass如果等于2会有好转吗?...

读读手册再来问,IBRION先设置成0试试再说。你有没有看过轨迹,原子核动了没有。

6.21.1 IBRION=-1
No update; ions are not moved, but NSW outer loops are performed. In each outer loop the electronic degrees of freedom are
re-optimized (for NSW>0 this obviously does not make much sense, except for test purposes). If no ionic update is required
use NSW=0 instead.
5楼2016-05-20 15:48:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skplin

铜虫 (小有名气)

引用回帖:
5楼: Originally posted by luzihen at 2016-05-20 15:48:58
读读手册再来问,IBRION先设置成0试试再说。你有没有看过轨迹,原子核动了没有。

6.21.1 IBRION=-1
No update; ions are not moved, but NSW outer loops are performed. In each outer loop the electronic d ...

分子动力学怎么看轨迹,,,手册上说分子动力学应该选IBRION=0   然后原子动了
6楼2016-05-21 10:18:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

skplin

铜虫 (小有名气)

引用回帖:
5楼: Originally posted by luzihen at 2016-05-20 15:48:58
读读手册再来问,IBRION先设置成0试试再说。你有没有看过轨迹,原子核动了没有。

6.21.1 IBRION=-1
No update; ions are not moved, but NSW outer loops are performed. In each outer loop the electronic d ...

IBRION=0   nsw= 500   smass= 1   温度都是1000   但是温度变化范围还是很大
7楼2016-05-21 10:20:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

luzihen

金虫 (正式写手)

引用回帖:
7楼: Originally posted by skplin at 2016-05-21 10:20:44
IBRION=0   nsw= 500   smass= 1   温度都是1000   但是温度变化范围还是很大...

什么叫大?
有波动很正常,没有fluctuation就不是NVT系综了
8楼2016-05-21 14:12:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 skplin 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-02 06:20 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-02 06:07 by xzcv0Tehd8lQ
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 xvJS1BYDxvBS 2026-10-31 10/500 2026-11-02 06:00 by xzcv0Tehd8lQ
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 xvJS1BYDxvBS 2026-10-31 11/550 2026-11-02 06:00 by xzcv0Tehd8lQ
[论文投稿] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-02 05:41 by xzcv0Tehd8lQ
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 18/900 2026-11-02 05:37 by xzcv0Tehd8lQ
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 hKe7U1kad5Sg 2026-10-31 9/450 2026-11-02 03:26 by xzcv0Tehd8lQ
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +6 1uaZt2kAZj2S 2026-10-31 9/450 2026-11-02 03:00 by xzcv0Tehd8lQ
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 1uaZt2kAZj2S 2026-10-30 16/800 2026-11-02 02:47 by xzcv0Tehd8lQ
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:27 by xzcv0Tehd8lQ
[找工作] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 hKe7U1kad5Sg 2026-10-31 5/250 2026-11-01 22:27 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 hKe7U1kad5Sg 2026-10-31 6/300 2026-11-01 22:11 by xzcv0Tehd8lQ
[考博] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 xvJS1BYDxvBS 2026-10-31 7/350 2026-11-01 22:01 by xzcv0Tehd8lQ
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 21:45 by xzcv0Tehd8lQ
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-30 8/400 2026-11-01 21:35 by xzcv0Tehd8lQ
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 21:19 by xzcv0Tehd8lQ
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 1uaZt2kAZj2S 2026-10-31 7/350 2026-11-01 12:40 by BCdvBAEpt3zm
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 1uaZt2kAZj2S 2026-10-31 8/400 2026-11-01 12:38 by BCdvBAEpt3zm
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1uaZt2kAZj2S 2026-10-30 17/850 2026-11-01 12:24 by BCdvBAEpt3zm
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 12:14 by BCdvBAEpt3zm
信息提示
请填处理意见