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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【原创】vasp的分子动力学模拟
451/1返回列表
查看: 21656  |  回复: 44
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
【奖励】 本帖被评价29次，作者xbelinda增加金币 25.1 个

xbelinda

银虫 (小有名气)

[资源] 【原创】vasp的分子动力学模拟

vasp做分子动力学的好处，由于vasp是近些年开发的比较成熟的软件，在做电子scf速度方面有较好的优势。
缺点：可选系综太少。
尽管如此，对于大多数有关分子动力学的任务还是可以胜任的。
主要使用的系综是 NVT 和 NVE。
下面我将对主要参数进行介绍！

一般做分子动力学的时候都需要较多原子，一般都超过100个。
当原子数多的时候，k点实际就需要较少了。有的时候用一个k点就行，不过这都需要严格的测试。通常超过200个原子的时候，用一个k点，即Gamma点就可以了。

INCAR：
EDIFF   一般来说，用1E-4 或者1E-5都可以，这个参数只是对第一个离子步的自洽影响大一些，对于长时间的分子动力学的模拟，精度小一点也无所谓，但不能太小。
IBRION=0 分子动力学模拟
IALGO=48 一般用48，对于原子数较多，这个优化方式较好。
NSW=1000   多少个时间步长。
POTIM=3  时间步长,单位fs, 通常1到3.
ISIF=2  计算外界的压力.
NBLOCK= 1  多少个时间步长，写一次CONTCAR，CHG和CHGCAR，PCDAT.
KBLOCK=50 NBLOCK*KBLOCK 个步长写一次 XDATCAR.
ISMEAR=-1  费米迪拉克分布.
SIGMA =0.05 单位:电子伏
NELMIN=8  一般用6到8, 最小的电子scf数.太少的话,收敛的不好.
LREAL=A
APACO=10 径向分布函数距离, 单位是埃.
NPACO=200  径向分布函数插的点数.
LCHARG=F 尽量不写电荷密度,否则CHG文件太大.
TEBEG=300  初始温度.
TEEND=300 终态温度。 不设的话，等于TEBEG.
SMASS   -3  NVE ensemble;-1 用来做模拟退火。大于0 NVT 系综。
回复此楼

» 收录本帖的淘帖专辑推荐

仿真建模与计算 VASP and MS VASP 第一性原理相关文档
MD分子动力学 VASP使用技巧 第一性原理计算辅助工具 分子动力学
VASP相关 我的第一性原理 缺陷计算 田凯
VASP资料 Vasp

» 猜你喜欢

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

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2009-12-20 14:32:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

wangyujia

木虫 (正式写手)

★★★ 三星级,支持鼓励

学习了，但是要指出你的一点小错误。
NBLOCK= 1  指的是多少个离子步（时间步长只是对于MD而言，对于其它的计算，NBLOCK也是起控制作用的）写一次XDATCAR
KBLOCK=50 NBLOCK*KBLOCK 个离子步写一次PCDAT.
CONTCAR是每个离子步之后都会写出来的，但是会用新的把老的覆盖
CHG是在每10个离子步写一次，不会覆盖
CHGCAR是在任务正常结束之后才写的。
赞一下(3人) 回复此楼
15楼2009-12-22 15:46:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

spur

金虫 (正式写手)

★★★ 三星级,支持鼓励

感谢分享。问一下楼主：有办法得到体系在某一个方向上(e.g.  x, y, z)的径向分布函数吗？另外楼主能不能给大家讲讲SMASS=1，2，3的区别？

[ Last edited by spur on 2009-12-21 at 13:00 ]
赞一下(2人) 回复此楼
5楼2009-12-21 12:58:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csfn

荣誉版主 (知名作家)

★★★★★ 五星级,优秀推荐

楼主能否考虑简单介绍下vasp做了md后的主要输出文件（主要结果）分析
赞一下(2人) 回复此楼
11楼2009-12-22 08:10:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yzcluster

金虫 (著名写手)
引用回帖:
Originally posted by csfn at 2009-12-22 08:10:
楼主能否考虑简单介绍下vasp做了md后的主要输出文件（主要结果）分析

偶曾经用CASTEP的MD算了一个晶格的结构相变问题，目前打算把手上那个表面的东西搞完后就做点QMD的东西。现在偶每天除了看点文献外就是马丁的《Electronic Structure》和《Ab Initio Molecular Dynamics: Basic Theory and Advanced Methods》两本书对比看，感触颇多啊！
赞一下(1人) 回复此楼
18楼2009-12-22 22:21:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peng.029

捐助贵宾 (小有名气)

★★★★★ 五星级,优秀推荐

受教了，感谢！不过还有好多具体的问题，比如数据后处理，分析等方面的经验和小工具可以共享吗？ 这里的后处理是不是和经典的分子动力学类似的呢？谢谢了……
赞一下(1人) 回复此楼
30楼2011-06-29 19:20:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

chuzhaonan

铁杆木虫 (著名写手)

★★★★★ 五星级,优秀推荐

这个要认真学习
回复此楼
2楼2009-12-20 16:38:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf

荣誉版主 (职业作家)

★★★★★ 五星级,优秀推荐

感谢分享
分子动力学模拟的参数可选择性也不小，不要完全被某个INCAR所局限，要明白每个参数是什么意思
赞一下 回复此楼
3楼2009-12-20 21:54:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

study163

木虫 (正式写手)

★★★★★ 五星级,优秀推荐

这个要看看，不错
回复此楼
4楼2009-12-20 22:30:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaowandouer

木虫 (正式写手)

★★★★★ 五星级,优秀推荐

楼主太无私了，感动ing。。。。。。。。
收藏起来，以备不时之需
赞一下 回复此楼
6楼2009-12-21 19:27:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf

荣誉版主 (职业作家)
引用回帖:
Originally posted by spur at 2009-12-21 12:58:
感谢分享。问一下楼主：有办法得到体系在某一个方向上(e.g.  x, y, z)的径向分布函数吗？另外楼主能不能给大家讲讲SMASS=1，2，3的区别？

[ Last edited by spur on 2009-12-21 at 13:00 ]

这个还真没有听说过，有一个问题是径向分布函数其本身就是一维的物理量
赞一下 回复此楼
7楼2009-12-21 20:05:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xbelinda

银虫 (小有名气)
引用回帖:
Originally posted by spur at 2009-12-21 12:58:
感谢分享。问一下楼主：有办法得到体系在某一个方向上(e.g.  x, y, z)的径向分布函数吗？另外楼主能不能给大家讲讲SMASS=1，2，3的区别？

[ Last edited by spur on 2009-12-21 at 13:00 ]

SMASS=1,2,3 是没有区别的。都是NVT ensemble。
SMASS只要是大于0就是NVT系综。
另一个问题，我也不会。没用到过。
赞一下 回复此楼
9楼2009-12-21 21:00:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gavinliu7390

木虫 (著名写手)

★★★★★ 五星级,优秀推荐

it is worth to learn.
赞一下 回复此楼
10楼2009-12-21 23:27:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuchenwf

荣誉版主 (职业作家)
引用回帖:
Originally posted by xbelinda at 2009-12-21 21:00:


SMASS=1,2,3 是没有区别的。都是NVT ensemble。
SMASS只要是大于0就是NVT系综。
另一个问题，我也不会。没用到过。

SMASS这个值在MD中，
不同值是对应于振动的频率，也就是热欲的频率和系统的频率要尽量一致。但是在实际中对计算结果影响不是特别大。
赞一下 回复此楼
12楼2009-12-22 08:42:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuan999

专家顾问 (知名作家)

★★★ 三星级,支持鼓励

收藏一下
回复此楼
13楼2009-12-22 09:19:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

qqpie

木虫 (职业作家)

★★★★★ 五星级,优秀推荐

感谢楼主啦~!!@@
赞一下 回复此楼
14楼2009-12-22 13:23:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

freer007

木虫 (正式写手)

★★★★★ 五星级,优秀推荐

谢谢楼主啊。
回复此楼
16楼2009-12-22 16:11:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xbelinda

银虫 (小有名气)
引用回帖:
Originally posted by wangyujia at 2009-12-22 15:46:
学习了，但是要指出你的一点小错误。
NBLOCK= 1  指的是多少个离子步（时间步长只是对于MD而言，对于其它的计算，NBLOCK也是起控制作用的）写一次XDATCAR
KBLOCK=50 NBLOCK*KBLOCK 个离子步写一次PCDAT.
CONTC ...

正解, 当时还是比较忙,有所疏忽!
赞一下 回复此楼
17楼2009-12-22 22:18:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xbelinda

银虫 (小有名气)
引用回帖:
Originally posted by yzcluster at 2009-12-22 22:21:


偶曾经用CASTEP的MD算了一个晶格的结构相变问题，目前打算把手上那个表面的东西搞完后就做点QMD的东西。现在偶每天除了看点文献外就是马丁的《Electronic Structure》和《Ab Initio Molecular Dynamics: Basi ...

看来你在憋大文章啊！
这有像你这样才能做出好工作啊！
赞一下 回复此楼
19楼2009-12-23 01:05:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lei0736

荣誉版主 (职业作家)

★★★★★ 五星级,优秀推荐

马恩靠这个搞了多少高档文章啊 呵呵
赞一下 回复此楼
21楼2010-01-14 21:43:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wjhwoods

铜虫 (正式写手)

en, 学习了。

en, 学习了。
回复此楼
23楼2010-02-05 04:27:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

enola

捐助贵宾 (正式写手)

★★★★★ 五星级,优秀推荐

学习拉，呵呵
回复此楼
24楼2010-06-04 16:17:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

longchang

金虫 (小有名气)

★★★★★ 五星级,优秀推荐

我先顶一个！
回复此楼
25楼2010-06-08 11:00:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

study163

木虫 (正式写手)
不错，总结得很详细
赞一下 回复此楼
26楼2010-06-09 18:48:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)
★
zzy870720z(金币+1):感谢交流，呵呵，细微的差别还是很有用处的 2010-06-10 08:10:57
引用回帖:
Originally posted by xbelinda at 2009-12-21 21:00:15:


SMASS=1,2,3 是没有区别的。都是NVT ensemble。
SMASS只要是大于0就是NVT系综。
另一个问题，我也不会。没用到过。

SMASS =1,2,3 虽说都是NVT 系统，但是，还是有区别的，SMASS 是控制温度震荡的频率的。要选择合适的值，以防止Nose-热浴和离子的运动脱耦，导致非正则的系宗。
赞一下(1人) 回复此楼
27楼2010-06-09 19:53:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

原子修补匠

铁虫 (初入文坛)
真的是非常好的东西，泛泛而谈容易，说仔细的难。
楼主这么细致的帖子，真是无私啊！
赞一下 回复此楼
28楼2010-06-14 15:37:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

juchen

金虫 (小有名气)

★★★★★ 五星级,优秀推荐

学习了~
回复此楼
29楼2011-04-08 14:42:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bird007

荣誉版主 (职业作家)
帖子真精彩!
已经收录到淘贴专辑《VASP》
赞一下 回复此楼
31楼2011-08-31 07:32:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huangll99

木虫 (职业作家)

★★★★★ 五星级,优秀推荐

ding
回复此楼
32楼2012-01-06 21:54:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wujing110

禁虫 (正式写手)
本帖内容被屏蔽
33楼2012-03-15 10:07:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

123左叉

新虫 (小有名气)

★★★★★ 五星级,优秀推荐

请问大侠，VASP跑分子动力学步骤是怎样的呢？是和做其他计算一样吗，结构优化，静态自洽，非静态自洽，求band吗
赞一下 回复此楼
34楼2014-05-06 20:56:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli_2009

新虫 (小有名气)

★★★★★ 五星级,优秀推荐

感谢楼主
回复此楼
35楼2014-07-02 17:18:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

尤图尔m5

铜虫 (小有名气)

★★★★★ 五星级,优秀推荐

有个问题,ISMEAR=-1  费米迪拉克分布.这个是必须要用-1吗,我取的是普通的0,这个会有影响吗?
赞一下 回复此楼
36楼2015-03-08 19:28:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Gphase

禁虫 (小有名气)
本帖内容被屏蔽
37楼2016-05-14 14:37:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gunfawk

新虫 (初入文坛)

★★★★★ 五星级,优秀推荐

引用回帖:
21楼: Originally posted by lei0736 at 2010-01-14 21:43:30
马恩靠这个搞了多少高档文章啊 呵呵

为何hehe，不理解
回复此楼
38楼2017-04-08 20:13:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zhangfrank

金虫 (正式写手)

★★★★★ 五星级,优秀推荐

非常好 顶一下
回复此楼
39楼2017-04-14 14:56:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

2longriver

新虫 (初入文坛)

★★★★★ 五星级,优秀推荐

无私奉献，楼主好人
赞一下 回复此楼
40楼2017-06-03 21:02:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

10133202

铁虫 (小有名气)

★★★★★ 五星级,优秀推荐

引用回帖:
18楼: Originally posted by yzcluster at 2009-12-22 22:21:49
偶曾经用CASTEP的MD算了一个晶格的结构相变问题，目前打算把手上那个表面的东西搞完后就做点QMD的东西。现在偶每天除了看点文献外就是马丁的《Electronic Structure》和《Ab Initio Molecular Dynamics: Basic  ...

亲，您在读这两本书：《Electronic Structure》和《Ab Initio Molecular Dynamics: Basic Theory and Advanced Methods》能否发我一份，非常感谢。邮箱14110720002@fudan.edu.cn.
赞一下 回复此楼
41楼2018-12-07 13:02:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

风度尽量翩

银虫 (正式写手)
您好，我这边摘录一下用公众号转载一下可以吗？给您注明出处
赞一下 回复此楼
42楼2019-09-09 16:33:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

weichengdo

新虫 (小有名气)
可以完整发一下具体的计算流程吗？（脚本的使用和静态计算一样的吗？）
赞一下 回复此楼
43楼2020-06-27 20:36:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

银虫 (著名写手)
引用回帖:
27楼: Originally posted by vasp001 at 2010-06-09 19:53:47
SMASS =1,2,3 虽说都是NVT 系统，但是，还是有区别的，SMASS 是控制温度震荡的频率的。要选择合适的值，以防止Nose-热浴和离子的运动脱耦，导致非正则的系宗。...

如何选择1 2 3?需要测试还是怎么样
赞一下 回复此楼
44楼2021-01-14 16:46:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hailong147

禁虫
本帖内容被屏蔽
45楼2021-11-25 10:19:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
aylayl088楼
2009-12-21 20:15   回复  
 
abcasi20楼
2010-01-14 17:06   回复  
 
zdhlover22楼
2010-01-15 21:26   回复  
 
相关版块跳转 我要订阅楼主 xbelinda 的主题更新
451/1返回列表
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭