24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 6061  |  回复: 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的回帖

李天帮

木虫 (正式写手)

小小木虫

引用回帖:
7楼: Originally posted by happyzyf at 2015-04-13 09:17:56
手册中说“For SMASS>=0 a canonical ensemble is simulated using the algorithm of Nosé. The Nosé mass controls the frequency of the temperature oscillations during the simulation (see . For SMASS ...

http://muchong.com/bbs/viewthread.php?tid=7106381&fpage=1&target=blank
http://muchong.com/bbs/viewthread.php?tid=7017928&fpage=1&target=blank;楼主这个其实我也是新手一般大家都是先NVT跑一下,然后NVE跑,但是NVT中SMASS选哪个,一般取>0,详细的话估计的看它的原代码.

» 本帖已获得的红花(最新10朵)

事虽难,做则成;路虽远,行则至。
10楼2015-04-14 08:53:06
已阅   回复此楼   关注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的回帖

李天帮

木虫 (正式写手)

小小木虫

引用回帖:
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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 284求调剂 +5 小熊~~ 2026-03-31 5/250 2026-03-31 10:22 by sophie2180
[考研] 085600材料与化工329分求调剂 +5 这是有太阳哇 2026-03-25 5/250 2026-03-31 09:59 by luoyongfeng
[考研] 一志愿郑大材料工程290求调剂 +12 Youth_ 2026-03-30 12/600 2026-03-31 03:34 by 蒙奇奇521
[考研] 285求调剂 +6 AZMK 2026-03-29 9/450 2026-03-30 21:02 by dophin1985
[考研] 346求调剂 一志愿070303有机化学 +7 萝卜炖青菜 2026-03-28 8/400 2026-03-30 19:59 by 无际的草原
[考研] 各位老师好,我的一志愿为北京科技大学085601材料专硕 +10 Koxui 2026-03-28 10/500 2026-03-30 19:33 by 源_2020
[考研] 085601材料工程找调剂 +17 oatmealR 2026-03-29 18/900 2026-03-30 19:21 by Wang200018
[考研] 化工专硕求调剂 +5 question挽风 2026-03-24 5/250 2026-03-30 18:34 by tianyyysss
[考研] 一志愿厦门大学材料工程专硕354找调剂!!! +5 贝呗钡钡 2026-03-30 5/250 2026-03-30 18:16 by 无际的草原
[考研] 材料化工340求调剂 +3 jhx777 2026-03-30 3/150 2026-03-30 17:54 by JourneyLucky
[考研] 334分 一志愿武理 材料求调剂 +16 李李不服输 2026-03-26 16/800 2026-03-30 17:47 by wangjy2002
[考研] 化学308分调剂 +7 你好明天你好 2026-03-30 8/400 2026-03-30 13:02 by 你好明天你好
[考研] 考研调剂 +7 小蜡新笔 2026-03-29 7/350 2026-03-29 19:00 by 学员8dgXkO
[考研] 327求调剂 +6 汲亦昊 2026-03-29 6/300 2026-03-29 13:40 by peike
[考研] 291求调剂 +5 Y-cap 2026-03-29 6/300 2026-03-29 13:18 by mumin1990
[考研] 343求调剂 +6 爱羁绊 2026-03-29 6/300 2026-03-29 12:00 by 无际的草原
[考研] 070300求调剂306分 +4 26要上岸 2026-03-27 4/200 2026-03-28 13:06 by 唐沐儿
[考研] 302求调剂 +4 锦衣卫藤椒 2026-03-25 4/200 2026-03-25 16:29 by 功夫疯狂
[考研] 347求调剂 +4 L when 2026-03-25 4/200 2026-03-25 13:37 by cocolv
[考研] 0854人工智能方向招收调剂 +4 章小鱼567 2026-03-24 4/200 2026-03-25 13:29 by 2177681040
信息提示
请填处理意见