| 查看: 1076 | 回复: 2 | |||
[交流]
【求助】求助vasp 分子动力学 已有2人参与
|
|
我想求助vasp做第一性原理分子动力学的脚本文件，因为我想用他做硅纳米线的拉伸和压缩，因此我想大家如果有分子动力学的脚本文件请发我的邮箱 zhangying6971@126.com那样的话我对着mannual 学可能更容易上手一些，我学软件的时间不多了，因此希望各位一起学的能帮助一下，不胜感激。。。。 |
回复此楼» 猜你喜欢
求个博导看看
已经有17人回复
-
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- VASP分子动力学优化的INCAR如下，怎么CONTCAR和POSCAR没有什么变化啊？
已经有13人回复
- 【求助】VASP做分子动力学模拟NVT系综下为啥内能竟然是升高的？
已经有8人回复
- 【任务征集】VASP跑分子动力学数据处理方法
已经有20人回复
- 【求助】关于vasp进行分子动力学模拟的求助
已经有12人回复
- 【求助】请教vasp计算恒温第一原理分子动力学
已经有9人回复
- 【求助】重金求助vasp 算过渡态 80 金币
已经有69人回复
- 【求助】如何利用分子动力学加速DFT计算的收敛
已经有6人回复
- 【原创】vasp的分子动力学模拟
已经有44人回复
beyondstar
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2357.3
- 散金: 601
- 红花: 4
- 帖子: 475
- 在线: 2503.5小时
- 虫号: 96793
- 注册: 2005-11-07
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2010-08-02 17:44:07
hedaors
禁虫 (正式写手)
隐清堂主
- 应助: 11 (小学生)
- 贵宾: 0.225
- 金币: 3948.8
- 散金: 201
- 红花: 8
- 帖子: 855
- 在线: 271.1小时
- 虫号: 449622
- 注册: 2007-11-02
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
★ ★ ★ ★ ★ ★
qasd(金币+3):谢谢 2010-08-03 07:24:47
aylayl08(金币+3):谢谢解答 2010-08-03 14:18:26
qasd(金币+3):谢谢 2010-08-03 07:24:47
aylayl08(金币+3):谢谢解答 2010-08-03 14:18:26
|
这是一个例子，也许对你有提示作用 #------INCAR------ IBRION = 0 # 分子动力学 ISTART = 0 # from scratch IALGO = 48 # 如果发现自恰不收敛的情况可换成 38。 NSW = 5000 # 时间步长，一般在 500 步左右就能达到平衡，但要想充分的话，最好 3000 以上。 NBLOCK = 1; KBLOCK = 50 # 每个 NBLOCK 步数后，对关联函数和 DOS 会被计算，然后离子位置信息写 到相应的 XDATCAR。在 KBLOCK*NBLOCK 步数后，对关联函数和 DOS 会被写到 PCDAT 和 DOSCAR。 ISYM = 0 # 取消对称性。 SMASS = -3 # -3 NVE ensemble, -1 模拟退火. >0 NVT ensemble POTIM = 2.0 # 时间步长。一般 1 到 3. 这个最好跑一下 NVE 系综，看一下总能是否守恒。 ISMEAR = -1; SIGMA = 0.1 # -1 是电子费米达拉克占据，主要用来处理电子含温。 LREAL = A ISIF = 2 # 就算压力。 NPACO = 200,APACO = 10 # 计算对关联分布的。具体可以看 manual，很简单。 NELMIN = 6 # 电子最少的自恰步数，一般是 4 到 8。 NBANDS = 700 # 能带数，这个最好自己设，否则容易出问题。 LCHARG = F # 尽量不写 CHG，太占硬盘资源。 TEBEG = 1000；TEEND=1000 # 初始温度和末态温度，默认是两者相当，所以一般设 TEBEG 就可以。 ##------END--------- PREC 一般也不用设，默认即可。ENCUT 会等于 ENMAX。我能想起来的参数就这些，有问题可以直接问 我。因为分子动力学是一个统计概念，当然原子越多越好的。所以我认为最少也得算 100 个原子。一般 100 个 原子的话，k 点用一个就可以。注意的是在分子动力学里，PSTRESS 是没有的用的。因为没有 NPH 或 NPT 系综。再有分子动力学常用的几个概念就是，对关联函数，这个实验可测。写在 PCDAT 里。速度自关联函数， 这个可以自己看一下公式，很简单，可以自己写程序做。 做分子动力学最好用 vasp5，因为很多参数都不用测了，4.6 好像还得测， 否则慢。还有实际 vasp 在输出温度是需要修正的，这个可以看看 manual，很简单，如果原子数多的话可以不 修正。还有压力也需要修正，但通常影响不大，可以直接用 vasp 的。 |
3楼2010-08-02 22:46:18