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

客户端APP下载

调剂小程序

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (3248)
>虫友互识 (354)
>文献求助 (203)
>休闲灌水 (180)
>导师招生 (154)
>考博 (78)
>博后之家 (50)
>硕博家园 (48)
>考研 (48)
>论文投稿 (40)
>基金申请 (34)
>论文道贺祈福 (31)
>公派出国 (30)
>找工作 (26)
>教师之家 (25)
>招聘信息布告栏 (23)

zhlhhc

铁虫 (小有名气)

应助: 1 (幼儿园)
金币: 43.8
散金: 641
帖子: 177
在线: 86小时
虫号: 2073069
注册: 2012-10-19
专业: 凝聚态物性 II ：电子结构

[求助] 关于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，
这究竟是什么原因呢？

» 猜你喜欢

为什么中国大学教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有3人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急已经有4人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急已经有3人回复
售SCI一区文章，我:8 O5 51O 54,科目齐全,可+急已经有4人回复
“人文社科而论，许多学术研究还没有达到民国时期的水平” 已经有5人回复
过年走亲戚时感受到了所开私家车的鄙视链已经有11人回复
什么是人一生最重要的？已经有4人回复
版面费该交吗已经有3人回复
今年春晚有几个节目很不错，点赞！已经有12人回复
体制内长辈说体制内绝大部分一辈子在底层，如同你们一样大部分普通教师忙且收入低已经有12人回复

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

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人回复

1楼 2015-05-06 15:14:40

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

dgwfly

新虫 (正式写手)

应助: 0 (幼儿园)
金币: 683.9
散金: 258
红花: 9
帖子: 543
在线: 202小时
虫号: 3454332
注册: 2014-10-04
性别: MM
专业: 数论

同求，楼主知道原因了吗

发自小木虫Android客户端

2楼2016-02-05 18:20:45

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

ztfsuiyuer

银虫 (正式写手)

应助: 0 (幼儿园)
金币: 1635
帖子: 311
在线: 74.4小时
虫号: 3424105
注册: 2014-09-17
性别: GG
专业: 理论和计算化学

MD模拟的全名是啥啊

发自小木虫Android客户端

3楼2016-04-09 18:18:33

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

zhlhhc

铁虫 (小有名气)

应助: 1 (幼儿园)
金币: 43.8
散金: 641
帖子: 177
在线: 86小时
虫号: 2073069
注册: 2012-10-19
专业: 凝聚态物性 II ：电子结构

引用回帖:

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 的主题更新