版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

返回列表

52qing

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 8
散金: 15
帖子: 23
在线: 51.1小时
虫号: 1404507
注册: 2011-09-16
专业: 凝聚态物性I:结构、力学和

[交流] 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

回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

1楼 2018-08-31 13:32:56

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 52qing 的主题更新