| 查看: 904 | 回复: 4 | |||
| 当前主题已经存档。 | |||
gavinliu7390木虫 (著名写手)
叶落鹰飞
|
[交流]
【求助】vasp跑分子动力学
|
||
|
有没有用vasp跑分子动力学的。可以分享一下输入文件或经验吗？ Thanks in advance! |
回复此楼» 猜你喜欢
一个有机合成实验室都需要哪些设备？
已经有7人回复
-
2026年国自然面上资助率
已经有20人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有22人回复
-
微信指数没变化，科研之友没阅读
已经有18人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有8人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
你们的时间戳变了吗
已经有4人回复
-
产物和副产物价值比较
已经有5人回复
-
时间戳他又来了
已经有17人回复
» 本主题相关商家推荐: (我也要在这里推广)
wuchenwf
荣誉版主 (职业作家)
- 应助: 0 (幼儿园)
- 贵宾: 3.433
- 金币: 19427.4
- 散金: 10
- 红花: 4
- 帖子: 3566
- 在线: 1037小时
- 虫号: 398569
- 注册: 2007-06-10
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
2楼2009-11-05 23:21:43
3楼2009-11-06 00:58:38
wangzg9292
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3933
- 散金: 56
- 帖子: 348
- 在线: 89.6小时
- 虫号: 723844
- 注册: 2009-03-16
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
4楼2009-11-06 13:28:11
★ ★ ★
gavinliu7390(金币+1):谢谢参与
aylayl08(金币+2,VIP+0):谢谢提供信息 11-6 14:41
gavinliu7390(金币+1):谢谢参与
aylayl08(金币+2,VIP+0):谢谢提供信息 11-6 14:41
|
http://cms.mpi.univie.ac.at/vasp/vasp/node233.html Molecular -- Dynamics About which files do you have to worry: param.inc INCAR POSCAR POTCAR KPOINTS Use the makeparam utility to create the param.inc file. For a molecular dynamics PREC= Low is definitely sufficient. The INCAR file might be similar to SYSTEM = Se ENCUT = 150 eV ! energy cutoff (opt) IALGO = 48 ! RMM-DIIS algorithm for electrons LREAL = A ! evaluate projection operators in real space NELMIN = 4 ! do a minimum of four electronic steps BMIX = 2.0 ! mixing parameter MAXMIX = 50 ! keep dielectric function between ionic movements Ionic Relaxation ISYM = 0 ! switch of symmetry NSW = 100 ! number of steps for IOM NBLOCK = 1 ; KBLOCK = 100 SMASS = 2.0 ! Nose mass-parameter (am) POTIM = 3.00 ! time-step for ion-motion TEBEG = 573 ! temperature PC-function APACO = 10.0 ! distance for P.C. Use IALGO=48 (RMM-DIIS for electrons) for large molecular dynamic runs. You should also evaluate the projection operators in real space (LREAL=A), and require at least 4 electronic iterations per ionic step (NELMIN = 4). For surface you might need to increase this value to NELMIN = 8. Special consideration require the parameters BMIX and MAXMIX: it is usually desirable to use optimal mixing parameters for molecular dynamics simulations. This can be done by performing a few static calculations with varying AMIX and BMIX parameters and do determine the one leading to the fastest convergence. However in the latest version of VASP the dielectric function is reused when the ions are updated (an optimal AMIX and BMIX is no longer that important). The dielectric function is reused after ionic updates, if MAXMIX ist set. MAXMIX should be about three times as larger as the number of iterations required to converge the electronic wavefunctions in the first iteration. After doing executing VASP once, it is only necessary to copy CONTCAR to POSCAR and to restart VASP. Usually a shell script is used for this task. An example shell script can be found on the vamp account in the file vamp/scripts/iter. |
5楼2009-11-06 13:29:00