24小时热门版块排行榜    

查看: 2492  |  回复: 3

zhlhhc

铁虫 (小有名气)

[求助] 关于NVT系综计算MD的一个疑问

ISTART = 0
ICHARG=2
ENCUT  = 400
ISMEAR =0
SIGMA = 0.05
IBRION=0
ISIF= 2
NSW = 1000
NELM=40
NPAR=4
#ISPIN =2
#LORBIT=11
LREAL=Auto
#EDIFF=1E-5
EDIFFG=-0.02
PREC = Normal
ALGO = Very Fast
MAXMIX = 40
NCORE= 4
ISYM = 0
NELMIN =6
NWRITE = 0
LCHARG = .FALSE.
LWAVE = .FALSE.
TEBEG = 300
TEEND = 300
SMASS = 1
NBLOCK = 50
POTIM = 1
这是我的参数设置,初温和终温都是300K,奇怪的是我的初始温度设置的是300K,模拟计算的时候初始温度却是0:
  1 T=     0. E= -.38355582E+03 F= -.38355591E+03 E0= -.38355591E+03  EK= 0.32792E-07 SP= 0.00E+00 SK= 0.94E-04
   2 T=     0. E= -.38355587E+03 F= -.38355597E+03 E0= -.38355597E+03  EK= 0.27653E-06 SP= -.76E-03 SK= 0.85E-03
   3 T=     0. E= -.38355587E+03 F= -.38355597E+03 E0= -.38355597E+03  EK= 0.69827E-06 SP= -.23E-02 SK= 0.24E-02
   4 T=     0. E= -.38355587E+03 F= -.38355597E+03 E0= -.38355597E+03  EK= 0.13173E-05 SP= -.45E-02 SK= 0.46E-02
   5 T=     0. E= -.38355587E+03 F= -.38355597E+03 E0= -.38355597E+03  EK= 0.21672E-05 SP= -.76E-02 SK= 0.76E-02
   6 T=     0. E= -.38355588E+03 F= -.38355597E+03 E0= -.38355597E+03  EK= 0.32043E-05 SP= -.11E-01 SK= 0.11E-01
   7 T=     0. E= -.38355588E+03 F= -.38355597E+03 E0= -.38355597E+03  EK= 0.44692E-05 SP= -.16E-01 SK= 0.16E-01
   8 T=     0. E= -.38355588E+03 F= -.38355597E+03 E0= -.38355597E+03  EK= 0.59776E-05 SP= -.21E-01 SK= 0.21E-01
   9 T=     0. E= -.38355588E+03 F= -.38355598E+03 E0= -.38355598E+03  EK= 0.76788E-05 SP= -.27E-01 SK= 0.27E-01
  10 T=     0. E= -.38355588E+03 F= -.38355598E+03 E0= -.38355598E+03  EK= 0.94657E-05 SP= -.34E-01 SK= 0.34E-01
  11 T=     0. E= -.38355588E+03 F= -.38355598E+03 E0= -.38355598E+03  EK= 0.11228E-04 SP= -.42E-01 SK= 0.42E-01
这说明我设置的初始温度没有起作用啊,这是为什么呢?
更奇怪的是,用同样的参数设置,算其他的体系,初始温度就是300K,
这究竟是什么原因呢?
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

dgwfly

新虫 (正式写手)

2楼2016-02-05 18:20:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ztfsuiyuer

银虫 (正式写手)

3楼2016-04-09 18:18:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhlhhc

铁虫 (小有名气)

引用回帖:
3楼: Originally posted by ztfsuiyuer at 2016-04-09 18:18:33
MD模拟的全名是啥啊

原因是计算MD时的初始POSCAR除原子坐标外,还有一大堆0000000,删掉就正常了,MD模拟的全名是ab initio molecular dynamics
4楼2016-04-10 09:08:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 zhlhhc 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[有机交流] 有机合成求助 20+3 FENGSHUJEI 2026-03-23 5/250 2026-03-24 19:31 by 88817753
[考研] 306求0703调剂一志愿华中师范 +10 纸鱼ly 2026-03-21 11/550 2026-03-24 17:22 by qingfeng258
[考研] 材料学硕,求调剂 6+3 糖葫芦888ll 2026-03-22 7/350 2026-03-24 17:11 by hello七七
[考研] 一志愿北京化工大学 070300 学硕 336分 求调剂 +7 vv迷 2026-03-22 7/350 2026-03-23 23:44 by Txy@872106
[考研] 一志愿国科过程所081700,274求调剂 +3 三水研0水立方 2026-03-23 3/150 2026-03-23 23:11 by MajorWen
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
[考研] 一志愿070300浙大化学358分,求调剂! +4 酥酥鱼.. 2026-03-21 4/200 2026-03-23 08:12 by Iveryant
[考研] 293求调剂 +3 涛涛Wjt 2026-03-22 5/250 2026-03-22 22:21 by jiangpengfei
[考研] 石河子大学(211、双一流)硕博研究生长期招生公告 +3 李子目 2026-03-22 3/150 2026-03-22 21:01 by 怎么释怀
[考研] 一志愿中南化学(0703)总分337求调剂 +9 niko- 2026-03-19 10/500 2026-03-22 16:08 by ColorlessPI
[考研] 260求调剂 +3 朱芷琳 2026-03-20 4/200 2026-03-22 15:12 by 朱芷琳
[考研] 286求调剂 +10 Faune 2026-03-21 10/500 2026-03-21 23:34 by 314126402
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 材料 336 求调剂 +3 An@. 2026-03-18 4/200 2026-03-21 01:39 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
信息提示
请填处理意见