| 查看: 2472 | 回复: 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， 这究竟是什么原因呢？ |
回复此楼» 猜你喜欢
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低
已经有13人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有11人回复
-
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？
已经有8人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有9人回复
-
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MD分子动力学NVT系综模拟时间选择
已经有4人回复
- 分子动力学的结果：看系综平均还是看时间演进？
已经有6人回复
- 关于npt系综设定？
已经有5人回复
- MD系综选择的问题？
已经有5人回复
- 分子模拟中MSD的问题讨论
已经有25人回复
- MS discover模块计算问题
已经有6人回复
- NPT系综下判断是否平衡的标志
已经有12人回复
- 关于NVT系统下的压力设定
已经有10人回复
- NPT系综动力学
已经有19人回复
- 分子动力学模拟中的系综问题
已经有5人回复
- 做MD计算时系综NVT中体积不变 但是这个体积怎么规定的 有具体大小吗？
已经有4人回复
- 三方晶体的晶胞体积怎么计算？应该有公式的吧？
已经有13人回复
- 第一原理动力学的系综选择，和各种系综的输入文件Example
已经有13人回复
- nvt系综，系统总能量震荡厉害
已经有9人回复
- 【求助】VASP做分子动力学模拟NVT系综下为啥内能竟然是升高的？
已经有8人回复
- NPT系综的平衡问题
已经有10人回复
- 计算表面时：Vacuum真空层和slab层厚度收敛性测试，保持超晶胞体积不变吗？？
已经有10人回复
- 【活动】vasp计算分子动力学讨论合集及专贴（筹划中）
已经有33人回复
- 【讨论】系综问题
已经有32人回复
- 【讨论】关于md中系统稳定的问题
已经有10人回复
- 【讨论】lammps中怎样检验模拟系统在NVT系综下弛豫是否达到平衡？
已经有6人回复
- 【原创】vasp的分子动力学模拟
已经有44人回复
ztfsuiyuer
银虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 1635
- 帖子: 311
- 在线: 74.4小时
- 虫号: 3424105
- 注册: 2014-09-17
- 性别: GG
- 专业: 理论和计算化学
3楼2016-04-09 18:18:33
2楼2016-02-05 18:20:45
zhlhhc
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 43.8
- 散金: 641
- 帖子: 177
- 在线: 86小时
- 虫号: 2073069
- 注册: 2012-10-19
- 专业: 凝聚态物性 II ：电子结构
4楼2016-04-10 09:08:11
