24小时热门版块排行榜    

查看: 6395  |  回复: 15
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

happyzyf

铜虫 (小有名气)

[求助] vasp计算分子动力学温度变化不正常 已有1人参与

各位大侠好!我再用vasp的分子动力学做一个模拟退火的过程,感觉结果有问题,OSZICAR里面的温度是从0开始的。INCAR设置如下
NPAR   = 10
LPLANE = .TRUE.
LREAL  = AUTO
ICHARG = 1

ENCUT  = 300
PREC   = Normal
ALGO   = VeryFast
NELMIN = 4
NELM   = 150
EDIFF  = 1E-4
ISMEAR = 0
SIGMA  = 0.1
# ISPIN  = 2

# LWAVE  = .FALSE.
# LCHARG = .FALSE.

IBRION = 0
NSW    = 5000
POTIM  = 2
SMASS  = 3
TEBEG  = 774
TEEND  = 300
NBLOCK = 1
OSZICAR中温度从0开始上升,最高到1400多,然后开始下降,这中间还有震荡,不知道是什么原因。后来我改了INCAR里面的参数,固定体系始末温度都在300K,想让它平衡一段时间看看,设置如下
NPAR   = 10
LPLANE = .TRUE.
LREAL  = AUTO
ICHARG = 1

ENCUT  = 300
PREC   = Normal
ALGO   = VeryFast
NELMIN = 4
NELM   = 150
EDIFF  = 1E-4
ISMEAR = 0
# SIGMA  = 0.1
# ISPIN  = 2

LWAVE  = .FALSE.
LCHARG = .FALSE.

IBRION = 0
NSW    = 3000
POTIM  = 1
SMASS  = 2
TEBEG  = 300
TEEND  = 300
NBLOCK = 1
KBLOCK = 50
OSZICAR里面温度仍然是从0开始上升
多谢~
回复此楼

» 猜你喜欢

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

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

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

李天帮

木虫 (正式写手)

小小木虫

引用回帖:
3楼: Originally posted by happyzyf at 2015-04-08 08:56:40
谢谢回复,我对这三个系综不太了解,对于我第一个INCAR应该设置SMASS=-1是吗?对于第二个INCARSMASS=2有什么问题呢?我看大家通常都是用的大于0的NVT系综。。。...

退火的话就用-1,如果是控制温度用NVT比较好。
事虽难,做则成;路虽远,行则至。
4楼2015-04-10 13:27:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 16 个回答

李天帮

木虫 (正式写手)

小小木虫

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ...
感谢参与,应助指数 +1
ljw4010: 金币+2, 谢谢交流! 2015-04-08 22:54:27
happyzyf: 金币+50, ★★★★★最佳答案, 谢谢 2015-04-12 18:58:54
SMASS参数设置不对啊,-3为NVE系综,-1为退火,>0为NVT系综。
事虽难,做则成;路虽远,行则至。
2楼2015-04-08 08:53:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happyzyf

铜虫 (小有名气)

引用回帖:
2楼: Originally posted by 李天帮 at 2015-04-08 08:53:00
SMASS参数设置不对啊,-3为NVE系综,-1为退火,>0为NVT系综。

谢谢回复,我对这三个系综不太了解,对于我第一个INCAR应该设置SMASS=-1是吗?对于第二个INCARSMASS=2有什么问题呢?我看大家通常都是用的大于0的NVT系综。。。
3楼2015-04-08 08:56:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

happyzyf

铜虫 (小有名气)

引用回帖:
4楼: Originally posted by 李天帮 at 2015-04-10 13:27:11
退火的话就用-1,如果是控制温度用NVT比较好。...

谢谢回复。
NVT中SMASS的取值可以取0-3,具体的数值使用情况能否给个简单介绍呢?
5楼2015-04-12 18:59:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天系统多次维护,明天很可能放榜! +5 zju2000 2026-08-16 5/250 2026-08-16 23:23 by Haru815
[论文投稿] 岩土工程学报什么时候才能终审完呐 25+3 yeager111 2026-08-10 4/200 2026-08-16 20:42 by tfang
[基金申请] 2027广东省杰青 +3 奶牛小黑 2026-08-15 6/300 2026-08-16 20:07 by 奶牛小黑
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +9 Ldrop2023 2026-08-13 9/450 2026-08-16 13:24 by wukun_007
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +7 启萌科技 2026-08-12 26/1300 2026-08-16 12:35 by 启萌科技
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] filecode=后面第一个是大写字母 +6 wangze12014 2026-08-14 7/350 2026-08-15 20:12 by gltch
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +12 Tide man 2026-08-10 13/650 2026-08-15 16:34 by 氺木
[基金申请] 关于Filecode分析方法 +10 majunge000 2026-08-10 13/650 2026-08-15 12:26 by hanpeng972
[基金申请] 是这周出结果还是下周出结果? +4 yuleib84 2026-08-11 4/200 2026-08-14 23:05 by lfy8008
[硕博家园] 读博的好处 +4 lnee 2026-08-11 4/200 2026-08-14 10:20 by ahsoarli
[基金申请] filecode +15 documentary 2026-08-10 17/850 2026-08-14 10:08 by kissu88
[基金申请] FileCode能看出啥? +10 要乐观耀哥 2026-08-10 32/1600 2026-08-14 09:37 by 要乐观耀哥
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 分享一下我之前已中青C的计划书的filecode +4 布布和一二 2026-08-11 5/250 2026-08-13 12:56 by cratir
[基金申请] 结合人工智能,周易传统文化,filecode打分制来了,3分以上希望很大。 +3 Tide man 2026-08-12 4/200 2026-08-13 08:35 by ZJTJZ
[基金申请] 帮忙看看fileCode +7 wwncly 2026-08-10 13/650 2026-08-11 19:36 by 冰心玉壶晴
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
信息提示
请填处理意见