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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP，aimd计算无法达到设定温度，请指点一二
31/1返回列表
查看: 1717  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【悬赏金币】回答本帖问题，作者wfbgc将赠送您 50 个金币

wfbgc

木虫 (小有名气)

[求助] VASP，aimd计算无法达到设定温度，请指点一二 已有2人参与

小弟我用VASP进行aimd计算，NPT系综，设定600K，始终无法达到设定测试，相差200K左右。恳求各位大牛指点一二，跪谢！！
以下为INCAR文件：
##############################Global Parameters###############################

SYSTEM = aimd
ISTART = 0             # (Read existing wavefunction; if there)
ICHARG = 2
LREAL  = Auto          # (Projection operators: automatic)

ENCUT = 390            #  (Cut-off energy for plane wave basis set, in eV)
PREC = Normal          #  Accurate        #  (Precision level)
#SYMPREC = 1.0E-6

ISYM = -1

LWAVE = .F.            #  (Write WAVECAR or not)
LCHARG = .F.           #  (Write CHGCAR or not)


#############################Electronic Relaxation#############################   

ISMEAR = 0            #  (Gaussian smearing; metals:1)
SIGMA  =  0.1         #  (Smearing value in eV; metals:0.2)
NELM = 100             #  (Max electronic SCF steps)
NELMIN =  5           #  (Min electronic SCF steps)
EDIFF = 1.0E-4        # (SCF energy convergence; in eV)

IVDW = 12

############################Ionic Relaxation########################################


IBRION =  0           # (Algorithm: 0-MD; 1-Quasi-New; 2-CG)
POTIM = 1             # (分子动力学步长，一般为1 fs)
NSW = 500            # (设置动力学步数)
TEBEG = 600           # (设置初始温度)
#TEEND = 600           # (设置结束温度)


PSTRESS = 0


MDALGO = 3            # (1 for Andersen, 2 for Nose-Hoover, 3 for Langevin, et al.)
#ANDERSEN_PROB = 1

ISIF   =  3           # (Stress/relaxation: 2-Ions, 3-Shape/Ions/V, 4-Shape/Ions)
EDIFFG = -0.02        # (Ionic convergence; eV/AA)


############################附加关键词########################################

NBLOCK = 1            # After NBLOCK ionic steps the pair-correlation function and the DOS are calculated and the ionic configuration is written to the XDATCAR-file

KBLOCK = 20         # After KBLOCK*NBLOCK ionic steps the averaged pair-correlation function and DOS are written to the files PCDAT and DOSCAR.

SMASS = 1

###########################################

ALGO = FAST
#IALGO = 38
#LMAXMIX = 4


以下为部分log文件：

    327 T=   371. E= -.11023993E+04 F= -.11090206E+04 E0= -.11090206E+04  EK= 0.66213E+01 SP= 0.00E+00 SK= 0.00E+00
bond charge predicted
prediction of wavefunctions
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.107294984080E+04    0.31255E-01   -0.77548E-03  2040   0.273E+00    0.730E-02
RMM:   2    -0.107295130339E+04   -0.14626E-02   -0.16162E-02  1924   0.163E-01    0.744E-02
RMM:   3    -0.107295169989E+04   -0.39650E-03   -0.48419E-03  2219   0.720E-02    0.551E-02
RMM:   4    -0.107295174864E+04   -0.48749E-04   -0.15441E-03  1986   0.429E-02    0.229E-02
RMM:   5    -0.107295175757E+04   -0.89292E-05   -0.30400E-04  1410   0.277E-02
    328 T=   369. E= -.11023992E+04 F= -.11089849E+04 E0= -.11089849E+04  EK= 0.65857E+01 SP= 0.00E+00 SK= 0.00E+00
bond charge predicted
prediction of wavefunctions
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.107293080270E+04    0.20946E-01   -0.68796E-03  1816   0.271E+00    0.710E-02
RMM:   2    -0.107293215553E+04   -0.13528E-02   -0.15078E-02  1920   0.165E-01    0.611E-02
RMM:   3    -0.107293256296E+04   -0.40743E-03   -0.48554E-03  2178   0.727E-02    0.417E-02
RMM:   4    -0.107293264790E+04   -0.84938E-04   -0.15190E-03  1937   0.415E-02    0.200E-02
RMM:   5    -0.107293265800E+04   -0.10108E-04   -0.34766E-04  1431   0.297E-02
    329 T=   368. E= -.11023992E+04 F= -.11089595E+04 E0= -.11089595E+04  EK= 0.65603E+01 SP= 0.00E+00 SK= 0.00E+00
bond charge predicted
prediction of wavefunctions
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.107292292933E+04    0.97186E-02   -0.10482E-02  2072   0.270E+00    0.808E-02
RMM:   2    -0.107292473009E+04   -0.18008E-02   -0.19828E-02  1929   0.182E-01    0.601E-02
RMM:   3    -0.107292515743E+04   -0.42734E-03   -0.52410E-03  2158   0.778E-02    0.285E-02
RMM:   4    -0.107292526876E+04   -0.11133E-03   -0.15188E-03  2054   0.414E-02    0.152E-02
RMM:   5    -0.107292527059E+04   -0.18337E-05   -0.15091E-04  1210   0.274E-02
    330 T=   367. E= -.11023991E+04 F= -.11089458E+04 E0= -.11089458E+04  EK= 0.65466E+01 SP= 0.00E+00 SK= 0.00E+00
bond charge predicted
prediction of wavefunctions
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.107292836787E+04   -0.30991E-02   -0.79190E-03  2048   0.270E+00    0.793E-02
RMM:   2    -0.107293001206E+04   -0.16442E-02   -0.18159E-02  1915   0.174E-01    0.871E-02
RMM:   3    -0.107293045447E+04   -0.44241E-03   -0.54278E-03  2197   0.765E-02    0.702E-02
RMM:   4    -0.107293049320E+04   -0.38726E-04   -0.19311E-03  1963   0.477E-02    0.297E-02
RMM:   5    -0.107293053023E+04   -0.37025E-04   -0.68301E-04  1687   0.320E-02
    331 T=   367. E= -.11023992E+04 F= -.11089447E+04 E0= -.11089447E+04  EK= 0.65455E+01 SP= 0.00E+00 SK= 0.00E+00
bond charge predicted
prediction of wavefunctions
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.107294676754E+04   -0.16274E-01   -0.46392E-03  1712   0.269E+00    0.780E-02
RMM:   2    -0.107294803313E+04   -0.12656E-02   -0.14703E-02  1850   0.164E-01    0.593E-02
RMM:   3    -0.107294841847E+04   -0.38534E-03   -0.46336E-03  2181   0.723E-02    0.456E-02
RMM:   4    -0.107294848955E+04   -0.71085E-04   -0.13568E-03  1953   0.391E-02    0.219E-02
RMM:   5    -0.107294852395E+04   -0.34396E-04   -0.61234E-04  1662   0.301E-02
    332 T=   368. E= -.11023993E+04 F= -.11089564E+04 E0= -.11089564E+04  EK= 0.65571E+01 SP= 0.00E+00 SK= 0.00E+00
bond charge predicted
prediction of wavefunctions
       N       E                     dE             d eps       ncg     rms          rms(c)
DAV:   1    -0.107297702579E+04   -0.28536E-01   -0.62680E-03  1904   0.270E+00    0.823E-02
RMM:   2    -0.107297822101E+04   -0.11952E-02   -0.13954E-02  1892   0.159E-01    0.601E-02
RMM:   3    -0.107297856207E+04   -0.34105E-03   -0.43796E-03  2126   0.705E-02    0.254E-02
RMM:   4    -0.107297865386E+04   -0.91794E-04   -0.12724E-03  1937   0.383E-02    0.142E-02
RMM:   5    -0.107297865594E+04   -0.20798E-05   -0.16222E-04  1141   0.267E-02
    333 T=   369. E= -.11023993E+04 F= -.11089803E+04 E0= -.11089803E+04  EK= 0.65809E+01 SP= 0.00E+00 SK= 0.00E+00
bond charge predicted
prediction of wavefunctions
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)
1楼 2023-06-24 17:02:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)

【答案】应助回帖

就模拟了0.3ps？显然热浴未平衡。
赞一下 回复此楼
2楼2023-08-15 00:59:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

软木塞

新虫 (小有名气)

【答案】应助回帖

要么换个热浴方法，要么搞久一点，1ps起步，5ps合适，10ps基本足够
赞一下 回复此楼
3楼2023-08-24 13:12:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wfbgc 的主题更新
31/1返回列表
不应助 确定回帖应助 (注意：应助才可能被奖励，但不允许灌水，必须填写15个字符以上)
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 读博申请 +5 感dd 2026-03-10 7/350 2026-03-11 17:02 by QGZDSYS
[考研] 一志愿天津大学材料与化工275求调剂 10+5 穿只靴子 2026-03-07 20/1000 2026-03-11 16:25 by 醉在风里
[考研] 311求调剂 +6 zchqwer 2026-03-10 6/300 2026-03-11 15:54 by 30660438
[考研] 081200计算机科学与技术数一英一306求调剂 +3 intankt 2026-03-08 3/150 2026-03-11 10:50 by dandan413
[考研] 求调剂 +4 鹤遨予卿 2026-03-09 4/200 2026-03-11 10:47 by @飒飒飒飒
[考研] 大连大学化学专业研究生调剂 +3 琪久. 2026-03-10 8/400 2026-03-11 10:02 by 琪久.
[基金申请] 第二个面上项目的获批几率如何？ +5 wawm9296 2026-03-04 7/350 2026-03-11 09:59 by 最爱编纂辞典
[考研] 298求调剂 +3 Vv呀！ 2026-03-10 3/150 2026-03-10 22:40 by 剑诗杜康
[考研] 一志愿211，0860总分286食品工程专业求调剂 +8 林林Winnie 2026-03-05 8/400 2026-03-10 18:44 by 清风月
[考研] 308求调剂 +4 是Lupa啊 2026-03-08 7/350 2026-03-10 18:25 by 清风月
[考研] 一志愿天津大学，英一数二305分求调剂，四六级已过 +7 小小番的茄 2026-03-09 7/350 2026-03-10 16:48 by ztnimte
[考研] 一志愿湖师大化学289求调剂 +5 XMCMM3.14159 2026-03-10 5/250 2026-03-10 15:41 by shangxh
[考研] 材料专硕调剂 +7 慕辰123 2026-03-05 10/500 2026-03-10 14:22 by peike
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
[考研] 334求调剂 +5 Trying] 2026-03-07 5/250 2026-03-09 14:47 by ruiyingmiao
[考研] 296求调剂 +4 Xinyu Wu311 2026-03-09 4/200 2026-03-09 11:51 by 勇敢太监王公公
[考研] 调剂 +3 13853210211 2026-03-08 3/150 2026-03-08 19:34 by ms629
[考研] 0817化学工程与技术312分求调剂 +7 T123 tt 2026-03-04 7/350 2026-03-08 11:18 by 我的船我的海
[考研] 求调剂，学校研究所都可以，材料与化工267分 +6 wmx1 2026-03-05 6/300 2026-03-05 20:01 by houyaoxu
[考研] 316求调剂 +3 林小星发大财 2026-03-04 3/150 2026-03-05 07:49 by bxbo
信息提示
关闭
请填处理意见
关闭