24小时热门版块排行榜    

Znn3bq.jpeg
查看: 6118  |  回复: 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的回帖

vasp001

金虫 (正式写手)

引用回帖:
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 ...

我也曾被这个问题广泛困扰过。
一定要精通MaterialsStudios,VASP.
8楼2015-04-13 12:09:09
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 271求调剂 +37 2261744733 2026-04-11 39/1950 2026-04-17 10:11 by 黑科技矿业
[考研] 22专硕求调剂 +10 haoyun上岸 2026-04-11 12/600 2026-04-16 22:21 by 猪会飞
[基金申请] RY:中国产出的科学垃圾论文,绝对数量和比例都世界第一 +7 zju2000 2026-04-14 18/900 2026-04-16 11:36 by 欢乐颂叶蓁
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 药学305求调剂 +7 玛卡巴卡boom 2026-04-11 7/350 2026-04-15 13:21 by 西北望—风沙
[考研] 材料工程281还有调剂机会吗 +43 xaw. 2026-04-11 44/2200 2026-04-15 12:46 by 西北望—风沙
[考研] 211本科材料化工求调剂 +19 YHLAH 2026-04-11 23/1150 2026-04-14 22:25 by fenglj492
[考研] 272分材料子求调剂 +41 Loy0361 2026-04-10 54/2700 2026-04-14 18:00 by lhj2009
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[考研] 调剂 +12 月@163.com 2026-04-11 12/600 2026-04-14 15:37 by zs92450
[考研] 求调剂,985材料与化工348分 +9 涵竹刘 2026-04-11 14/700 2026-04-13 22:26 by 涵竹刘
[考研] 2026硕士调剂_能动_河南农业大学 +4 河南农业大学-能 2026-04-12 4/200 2026-04-13 22:01 by bljnqdcc
[考研] 一志愿中南大学 0855 机械 286 求调剂 +11 不会吃肉 2026-04-12 11/550 2026-04-13 21:59 by bljnqdcc
[考研] B区0809 ,数一英一,290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 一志愿985机械学硕380求调剂 +5 关关雎鸠10 2026-04-11 5/250 2026-04-11 10:10 by 知念。A
[考研] 281求调剂 +11 觉得好的吧 2026-04-10 11/550 2026-04-11 09:35 by 逆水乘风
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 本9 一志愿西工大085601 324求调剂 +5 wysyjs25 2026-04-10 5/250 2026-04-10 16:57 by luoyongfeng
信息提示
请填处理意见