24小时热门版块排行榜    

查看: 1259  |  回复: 0

52qing

新虫 (初入文坛)

[交流] nvt系综下 分子动力模拟 温度没有波动?

运行vasp分子动力模拟,发现温度一直保持在1600 K。按理说温度应该有个波动。这是为什么,请指教。

这是INCAR:
ISTART = 0 ; ICHARG=2
ENCUT  =  350
ISMEAR = 1; SIGMA = 0.1;
EDIFF=0.0001
#EDIFFG=0.0001
PREC =Accurate
NSW=2000
IBRION=0
ISIF=2
POTIM=1
LREAL= Auto
LWAVE=.FALSE.
IALGO=48
NBLOCK=1 #(xdatcar)
KNLOCK=50 #(pcdat)
NELMIN=6
APACO=8
NPACO=256
LCHARG=F
TEBEG= 1600  #initial temperature
TEEND= 1600  #end temperature
SMASS=0 #-3,NVE;0=<,NVT;

观察OSZICAR:
                    17 T=  1600. E= -.38408436E+03 F= -.41035015E+03 E0= -.41035833E+03  EK= 0.26266E+02 SP= 0.00E+00 SK= 0.29E-30
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1    -0.400443322025E+03    0.99068E+01   -0.24980E-02   435   0.150E-01    0.361E-02
RMM:   2    -0.400443365162E+03   -0.43137E-04   -0.85710E-04   377   0.329E-02    0.194E-02
RMM:   3    -0.400443359729E+03    0.54329E-05   -0.35693E-05   299   0.118E-02    0.133E-02
RMM:   4    -0.400443357240E+03    0.24887E-05   -0.98096E-06   270   0.654E-03    0.678E-03
RMM:   5    -0.400443357291E+03   -0.51237E-07   -0.30765E-06   201   0.414E-03    0.202E-03
RMM:   6    -0.400443357330E+03   -0.39056E-07   -0.70952E-07   155   0.309E-03
   18 T=  1600. E= -.37417757E+03 F= -.40044336E+03 E0= -.40045174E+03  EK= 0.26266E+02 SP= 0.00E+00 SK= 0.32E-30
       N       E                     dE             d eps       ncg     rms          rms(c)
RMM:   1    -0.389579708957E+03    0.10864E+02   -0.19781E-02   436   0.142E-01    0.323E-02
RMM:   2    -0.389579754317E+03   -0.45360E-04   -0.69683E-04   364   0.305E-02    0.175E-02
RMM:   3    -0.389579755743E+03   -0.14266E-05   -0.29629E-05   299   0.112E-02    0.121E-02
RMM:   4    -0.389579753288E+03    0.24558E-05   -0.76749E-06   246   0.604E-03    0.514E-03
RMM:   5    -0.389579753461E+03   -0.17318E-06   -0.21662E-06   192   0.401E-03    0.125E-03
RMM:   6    -0.389579753512E+03   -0.51539E-07   -0.35787E-07   149   0.309E-03
   19 T=  1600. E= -.36331396E+03 F= -.38957975E+03 E0= -.38958786E+03  EK= 0.26266E+02 SP= 0.00E+00 SK= 0.36E-30
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 52qing 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +15 低垂的野花 2026-07-31 20/1000 2026-08-01 14:10 by gltch
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 18/900 2026-08-01 13:56 by gltch
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 8/400 2026-08-01 13:53 by y4XuShGyllvT
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:51 by y4XuShGyllvT
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:51 by y4XuShGyllvT
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 7/350 2026-08-01 13:49 by y4XuShGyllvT
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:49 by y4XuShGyllvT
[基金申请] 小木虫看见有人已经知道结果了 +16 1234567wang 2026-07-25 20/1000 2026-08-01 08:57 by type2type
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +11 winnerche 2026-07-29 15/750 2026-07-30 23:08 by jnhyjjm
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 你们的时间戳变了吗 +3 archvillain 2026-07-30 4/200 2026-07-30 18:53 by levinzhwen
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见