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

登录注册

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

24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

liuchenhan

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1536.5
散金: 4
帖子: 170
在线: 108小时
虫号: 1508714
注册: 2011-11-25
专业: 传热传质学

[求助] vasp 分子动力学温度越来越小已有1人参与

本人初学vasp分子动力学模拟，遇到一个比较棘手的问题。平衡态下发现系统的温度越来越低。INCAR文件如下：
system = Gr_MD_300K

ISTART = 0
ICHARGE = 2
ENCUT = 520

NSW = 500
IBRION = 0
POTIM = 1
SMASS = 0

NELMDL = -5
NELM = 200

ISMEAR = -5

PREC = normal
NGX  =72
NGY  =80
NGZ  =120

ISIF = 3
GGA = PE
IALGO = 38

LREAL = Auto
LCHARG = .FALSE.
LWAVE = .FALSE.
NELMIN = 8
TEBEG = 290.625
NBLOCK = 10
KBLOCK = 1

温度数据如下（OSZICAR）：
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1    0.115844102433E+04 0.11584E+04 -0.48542E+04 866 0.120E+03
DAV: 2 -0.902289276473E+02 -0.12487E+04 -0.11843E+04  1195 0.282E+02
DAV: 3 -0.305305726723E+03 -0.21508E+03 -0.21414E+03 996 0.124E+02
DAV: 4 -0.313091051885E+03 -0.77853E+01 -0.77424E+01 986 0.296E+01
DAV: 5 -0.313329284920E+03 -0.23823E+00 -0.23799E+00  1124 0.447E+00 0.212E+01
DAV: 6 -0.300006068430E+03 0.13323E+02 -0.12837E+01 908 0.102E+01 0.112E+01
DAV: 7 -0.294906713533E+03 0.50994E+01 -0.15362E+01 883 0.110E+01 0.197E+00
DAV: 8 -0.294797428997E+03 0.10928E+00 -0.30374E-01 926 0.174E+00 0.322E-01
DAV: 9 -0.294765289986E+03 0.32139E-01 -0.52046E-02 900 0.587E-01 0.264E-01
DAV:  10 -0.294755571584E+03 0.97184E-02 -0.14646E-02 836 0.310E-01 0.143E-01
DAV:  11 -0.294752651723E+03 0.29199E-02 -0.13985E-02 812 0.277E-01 0.120E-01
DAV:  12 -0.294762641447E+03 -0.99897E-02 -0.45852E-03 812 0.189E-01 0.424E-02
DAV:  13 -0.294770431420E+03 -0.77900E-02 -0.15427E-03 725 0.881E-02 0.592E-02
DAV:  14 -0.294773035823E+03 -0.26044E-02 -0.13563E-03 845 0.932E-02 0.218E-02
DAV:  15 -0.294773961754E+03 -0.92593E-03 -0.59405E-05 525 0.226E-02 0.902E-03
DAV:  16 -0.294774644126E+03 -0.68237E-03 -0.53953E-05 490 0.176E-02 0.325E-03
DAV:  17 -0.294774732713E+03 -0.88587E-04 -0.70659E-06 493 0.720E-03
1 T= 291. E= -.29360997E+03 F= -.29477473E+03 E0= -.29476846E+03  EK= 0.11648E+01 SP= 0.00E+00 SK= 0.11E-09
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.294678727663E+03 0.95916E-01 -0.40332E+00 826 0.184E+01 0.308E-01
DAV: 2 -0.294731284934E+03 -0.52557E-01 -0.54279E-01 915 0.187E+00 0.338E-01
DAV: 3 -0.294731721125E+03 -0.43619E-03 -0.69188E-03  1076 0.249E-01 0.331E-01
DAV: 4 -0.294730830045E+03 0.89108E-03 -0.51082E-04 839 0.103E-01 0.296E-01
DAV: 5 -0.294730067707E+03 0.76234E-03 -0.19543E-02 836 0.569E-01 0.248E-01
DAV: 6 -0.294729132040E+03 0.93567E-03 -0.39851E-03 821 0.257E-01 0.140E-01
DAV: 7 -0.294728861627E+03 0.27041E-03 -0.82913E-04 832 0.113E-01 0.792E-02
DAV: 8 -0.294728817390E+03 0.44237E-04 -0.21538E-04 542 0.559E-02
2 T= 280. E= -.29360862E+03 F= -.29472882E+03 E0= -.29472272E+03  EK= 0.11202E+01 SP= 0 2 T= 280. E= -.29360862E+03 F= -.29472882E+03 E0= -.29472272E+03  EK= 0.11202E+01 SP= 0.51E-05 SK= 0.51E-05DAV: 1 -0.294548471642E+03 0.18039E+00 -0.36416E+00 850 0.186E+01 0.271E-01
DAV: 2 -0.294606392170E+03 -0.57921E-01 -0.58813E-01 898 0.223E+00 0.355E-01
DAV: 3 -0.294607403044E+03 -0.10109E-02 -0.78934E-03  1188 0.279E-01 0.370E-01
DAV: 4 -0.294606991865E+03 0.41118E-03 -0.28385E-04 731 0.815E-02 0.338E-01
DAV: 5 -0.294605021689E+03 0.19702E-02 -0.40929E-03 831 0.262E-01 0.252E-01
DAV: 6 -0.294604185544E+03 0.83615E-03 -0.12049E-03 863 0.148E-01 0.195E-01
DAV: 7 -0.294603747802E+03 0.43774E-03 -0.89724E-04 831 0.117E-01 0.117E-01
DAV: 8 -0.294603588157E+03 0.15965E-03 -0.29277E-04 602 0.756E-02 0.659E-02
DAV: 9 -0.294603543452E+03 0.44705E-04 -0.26371E-04 656 0.663E-02
3 T= 249. E= -.29360646E+03 F= -.29460354E+03 E0= -.29459793E+03  EK= 0.99800E+00 SP= - 3 T= 249. E= -.29360646E+03 F= -.29460354E+03 E0= -.29459793E+03  EK= 0.99800E+00 SP= -.11E-02 SK= 0.17E-03DAV: 1 -0.294370573630E+03 0.23301E+00 -0.29607E+00 858 0.186E+01 0.235E-01
DAV: 2 -0.294429865733E+03 -0.59292E-01 -0.59489E-01 882 0.247E+00 0.343E-01
DAV: 3 -0.294430170661E+03 -0.30493E-03 -0.72432E-03  1146 0.284E-01 0.332E-01
DAV: 4 -0.294428883253E+03 0.12874E-02 -0.76441E-04 876 0.125E-01 0.274E-01
DAV: 5 -0.294427108632E+03 0.17746E-02 -0.82661E-03 838 0.367E-01 0.107E-01
DAV: 6 -0.294427101300E+03 0.73316E-05 -0.16730E-03 825 0.155E-01 0.635E-02
DAV: 7 -0.294427077267E+03 0.24033E-04 -0.23777E-04 542 0.605E-02 0.370E-02
DAV: 8 -0.294427063791E+03 0.13475E-04 -0.11320E-04 510 0.460E-02
4 T= 207. E= -.29360341E+03 F= -.29442706E+03 E0= -.29442213E+03  EK= 0.82842E+00 SP= - 4 T= 207. E= -.29360341E+03 F= -.29442706E+03 E0= -.29442213E+03  EK= 0.82842E+00 SP= -.63E-02 SK= 0.15E-02
   N    E                   dE          d eps    ncg    rms       rms(c)
DAV: 1 -0.294181957181E+03 0.24512E+00 -0.22279E+00 868 0.185E+01 0.213E-01
DAV: 2 -0.294239302609E+03 -0.57345E-01 -0.57274E-01 842 0.259E+00 0.324E-01
DAV: 3 -0.294239143143E+03 0.15947E-03 -0.65567E-03  1132 0.285E-01 0.296E-01
DAV: 4 -0.294237097469E+03 0.20457E-02 -0.36171E-03 856 0.257E-01 0.153E-01
DAV: 5 -0.294236778243E+03 0.31923E-03 -0.22958E-03 835 0.175E-01 0.610E-02
DAV: 6 -0.294236732430E+03 0.45813E-04 -0.29561E-04 608 0.746E-02 0.196E-02
DAV: 7 -0.294236742306E+03 -0.98755E-05 -0.97240E-05 489 0.370E-02 0.970E-03
DAV: 8 -0.294236744219E+03 -0.19130E-05 -0.22438E-05 523 0.193E-02
5 T= 162. E= -.29360038E+03 F= -.29423674E+03 E0= -.29423254E+03  EK= 0.64932E+00 SP= -.20E-01 SK= 0.68E-02
   N    E                   dE          d eps    ncg    rms       rms(c)

回复此楼

» 猜你喜欢

0703化学336分求调剂已经有6人回复
304求调剂已经有9人回复
293求调剂已经有13人回复
281求调剂（0805）已经有8人回复
环境领域全国重点实验室招收博士1-2名已经有3人回复
材料专硕306英一数二已经有10人回复
301求调剂已经有6人回复
一志愿天津大学化学工艺专业（081702）315分求调剂已经有7人回复
302求调剂已经有6人回复
26博士申请已经有3人回复

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

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

Vasp相比经典分子动力学的优缺点已经有11人回复
vasp分子动力学报错已经有16人回复
vasp表面和吸附分子的结构分别优化后，如何将它们放到一起已经有11人回复
VASP怎样计算水分子键长键角已经有9人回复
vasp测试SIGMA 已经有12人回复
关于吉布斯函数的问题已经有10人回复
vasp计算单个氧原子和分子的能量已经有20人回复
VASP分子动力学优化的INCAR如下，怎么CONTCAR和POSCAR没有什么变化啊？已经有13人回复
【求助】VASP做分子动力学模拟NVT系综下为啥内能竟然是升高的？已经有8人回复
【任务征集】VASP跑分子动力学数据处理方法已经有20人回复
【求助】vasp 算出的0k下的能量如何扩展到其他温度压力已经有12人回复
【活动】vasp计算分子动力学讨论合集及专贴（筹划中）已经有33人回复
【求助】关于vasp进行分子动力学模拟的求助已经有12人回复
【求助】新手求助vasp计算氧分子结合能的问题已经有12人回复
【求助】请教vasp计算恒温第一原理分子动力学已经有9人回复
【求助】vasp做分子动力学升温或降温模拟问题已经有5人回复
【求助】vasp分子动力学模拟求助已经有15人回复
【原创】vasp的分子动力学模拟已经有44人回复
【讨论】交流讨论下用VASP计算不同温度下材料的热力学性能已经有19人回复

1楼 2014-12-17 19:48:24

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

蓝莓微积分

银虫 (初入文坛)

应助: 4 (幼儿园)
金币: 252
帖子: 47
在线: 89.6小时
虫号: 3368606
注册: 2014-08-16
专业: 凝聚态物性I:结构、力学和

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
ljw4010: 金币+1, 谢谢交流！ 2014-12-19 18:50:43
liuchenhan: 金币+2 2014-12-20 16:35:17

引用回帖:

2楼: Originally posted by liuchenhan at 2014-12-17 19:50:18
温度应该稳定在300K（290.625K，32个原子），然而实际上统计的温度越来越低，动能EK越来越小，不知道什么原因？望大神指教。

你的时间步长加大点试试呢，另外你的程序貌似还没有跑完，跑完后你再看看有没有回到室温。

赞一下

回复此楼

3楼2014-12-19 16:06:01

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

查看全部 6 个回答

liuchenhan

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1536.5
散金: 4
帖子: 170
在线: 108小时
虫号: 1508714
注册: 2011-11-25
专业: 传热传质学

温度应该稳定在300K（290.625K，32个原子），然而实际上统计的温度越来越低，动能EK越来越小，不知道什么原因？望大神指教。

赞一下

回复此楼

2楼2014-12-17 19:50:18

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

liuchenhan

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1536.5
散金: 4
帖子: 170
在线: 108小时
虫号: 1508714
注册: 2011-11-25
专业: 传热传质学

引用回帖:

3楼: Originally posted by 蓝莓微积分 at 2014-12-19 16:06:01
你的时间步长加大点试试呢，另外你的程序貌似还没有跑完，跑完后你再看看有没有回到室温。...

跑了将近200步，还没有。而且温度波动特别大，不知道什么原因？不知道脚本有没有写错。

赞一下

回复此楼

4楼2014-12-20 16:36:09

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

蓝莓微积分

银虫 (初入文坛)

应助: 4 (幼儿园)
金币: 252
帖子: 47
在线: 89.6小时
虫号: 3368606
注册: 2014-08-16
专业: 凝聚态物性I:结构、力学和

【答案】应助回帖

引用回帖:

4楼: Originally posted by liuchenhan at 2014-12-20 16:36:09
跑了将近200步，还没有。而且温度波动特别大，不知道什么原因？不知道脚本有没有写错。...

一般做分子动力学原子数越多温度波动越小，相反原子数越少的话温度可能波动较大，你是不是原子取少了？不知道你模拟的是什么材料呢

赞一下

回复此楼

5楼2014-12-20 20:45:23

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

查看全部 6 个回答

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[考研] 301求调剂 +6	yy要上岸呀 2026-03-17	6/300	2026-03-17 23:58 by 星空星月
[考研] 一志愿天津大学化学工艺专业（081702）315分求调剂 +7	yangfz 2026-03-17	7/350	2026-03-17 23:57 by 星空星月
[基金申请] 被我言中：新模板不强调格式了，假专家开始管格式了 +4	beefly 2026-03-14	4/200	2026-03-17 22:04 by 黄鸟于飞Chao
[考研] 296求调剂 +5	大口吃饭身体健 2026-03-13	5/250	2026-03-17 21:05 by 不惑可乐
[考研] 328求调剂，英语六级551，有科研经历 +3	生物工程调剂 2026-03-16	8/400	2026-03-17 19:03 by Wangjingyue
[考研] 考研化学学硕调剂，一志愿985 +4	张vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 271求调剂 +12	生如夏花… 2026-03-11	14/700	2026-03-17 10:56 by lovewei0727
[考研] 326求调剂 +4	诺贝尔化学奖觊� 2026-03-15	7/350	2026-03-16 17:11 by 诺贝尔化学奖觊�
[考研] 283求调剂 +10	小楼。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 材料与化工一志愿南昌大学327求调剂推荐 +7	Ncdx123456 2026-03-13	8/400	2026-03-16 12:15 by karry wen
[考研] 294求调剂 +3	Zys010410@ 2026-03-13	4/200	2026-03-15 10:59 by zhq0425
[考研] 328求调剂 +3	5201314Lsy！ 2026-03-13	6/300	2026-03-14 15:31 by hyswxzs
[考研] 331求调剂（0703有机化学 +5	ZY-05 2026-03-13	6/300	2026-03-14 10:51 by Jy?
[考研] 招收0805（材料）调剂 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 求调剂（材料与化工327） +4	爱吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] 281求调剂 +9	Koxui 2026-03-12	11/550	2026-03-13 20:50 by Koxui
[硕博家园] 085600 260分求调剂 +3	天空还下雨么 2026-03-13	5/250	2026-03-13 18:46 by 天空还下雨么
[考研] 求调剂 +3	程雨杭 2026-03-12	3/150	2026-03-13 15:06 by JourneyLucky
[考研] 一志愿山大07化学 332分四六级已过本科山东双非求调剂！ +3	不想理你 2026-03-12	3/150	2026-03-13 14:18 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4	shenzxsn 2026-03-11	4/200	2026-03-13 11:55 by ymwdoctor

24小时热门版块排行榜

liuchenhan

[求助] vasp 分子动力学 温度越来越小 已有1人参与

» 猜你喜欢

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

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

蓝莓微积分

【答案】应助回帖

liuchenhan

liuchenhan

蓝莓微积分

【答案】应助回帖

[求助] vasp 分子动力学温度越来越小已有1人参与