| 查看: 1202 | 回复: 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 |
回复此楼» 猜你喜欢
全日制（定向）博士
已经有5人回复
-
假如你的研究生提出不合理要求
已经有10人回复
-
萌生出自己或许不适合搞科研的想法，现在跑or等等看？
已经有4人回复
-
Materials Today Chemistry审稿周期
已经有4人回复
-
参与限项
已经有3人回复
-
实验室接单子
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
» 本主题相关商家推荐: (我也要在这里推广)
