| 查看: 13790 | 回复: 35 | |||||||||||||||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||||||||||||||||
gavinliu7390木虫 (著名写手)
叶落鹰飞
|
[交流]
【其他】vasp相关技术有奖征集！ 已有26人参与
|
||||||||||||||||||
|
直至目前，国内第一性原理研究已经搞的如火如荼。各类软件也是五花八门，但目前实际流行的基于平面波技术主要就两类。一类就是MS中castep为代表，由于它操作起来简单，易于大众接受，所以很受欢迎。尽管与此，由于门槛低，属于你会我也会，现在已经很难纯靠这类技术发高档次文章了。另一类就是一vasp为代表的，由于必须在linux下操作，相对难一些。但由于它的计算速度快，而且目前很多文章都用这种技术，（当然我不是说castep发不了好文章，prl, nature还是有的。但一般都附加了其它技术）还是受到更多人的喜欢。这一类软件实际有很多，如abinit, quantum expresso(pwscf). 至于其它一些基于轨道例如LMTO，我在这里就不提了。当然 vasp 也不是万能的，很多功能也得用类似软件来完成，例如利用密度泛函微扰计算声子性质就得用abinit或pwscf. 在这里我是希望大家一起集思广益，把一些技术都整理出来，这样一是方便初学者学习，二是方便大家的交流，三是希望认识更多的朋友。 只要大家能够耐心把自己会的相关技术贴出来，即使是思路也好。就有机会得到我发放的5-10个金币。 当然在这里我也会陆续贴出我自己学会的相关技术。（注，因为有的研究小组有规定，不经导师允许，不可以提供绝密技术。所以在这里尽量提供自学或自创技术。） vasp计算能带，态密度我想这都是最基本的技术，在网上一搜也一堆。 计算声子谱，目前主要是用 dircet supercell 方法。工具是vasp+phonon(这个功能还是很强大的，只是需要买)，还有就是 vasp+phon, vasp +fropho , vasp +phonopy, 这都是免费的，我个人比较喜欢最后一个，一是问题少，二是计算的比较快。 详细的看这个帖子，http://muchong.com/bbs/viewthread.php?tid=1710612 我将的第二个技术就是利用vasp计算 ideal strength. 在计算tensile strength的时候，主要是首先在你想要拉伸的方向建立原胞。最好设置 为a为x方向，而且需要的是b和c 与 a垂直。 然后改变啊a的长度，优化b和c。 计算出受力即可。而且要注意每一次改变a都要用前一次的CONTCAR. 其中一个关键问题要改一下vasp代码，重新编译一下。 打开constr_cell_relax.F文件。 在31到35行讲的是如何单独优化x方向。这里你只要将x换成y和z即可。 例如在30行后加入这几行即可， SAVY=FCELL(2,2) SAVZ=FCELL(3,3) FCELL=0 FCELL(2,2)=SAVY FCELL(3,3)=SAVZ 这样即可计算tensile strength了。 好了，再想到什么我还会写。有问题也欢迎。希望大家积极参与！ [ Last edited by gavinliu7390 on 2010-4-3 at 17:33 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 仿真建模与计算 | VASP | 第一性原理相关文档 | 第一性原理计算-tg |
| VASP软件从入门到高深 | vasp计算资料 | 第一性原理计算经验 | ms and vasp学习 |
| 有效质量 | vasp的相关资料 | vasp 学习贴 | 量子化学计算与使用经验 |
| vasp | 我的分享 | soft study | 不错哦 |
» 本帖已获得的红花（最新10朵）
donggaomin» 猜你喜欢
国自然评审意见
已经有8人回复
-
国自然面上复盘~欢迎讨论
已经有14人回复
-
系统查不到
已经有11人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
为什么资助数各大高校都创新高，自己申请怎么就这么难
已经有12人回复
-
有没有仍没收到信息的
已经有3人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有8人回复
-
基金不中，共勉
已经有11人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有9人回复
-
2026年叶企孙基金
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 有关vasp的INCAR的一些问题
已经有5人回复
- VASP+PHON算声子谱
已经有9人回复
- 做第一性原理这么长时间一直在思考这个问题！
已经有12人回复
- 如何计算 lattice mismatch？
已经有6人回复
- 如何用VASP做第一性原理拉伸？
已经有7人回复
- 怎么计算 atomic spin density
已经有4人回复
- Abinit能带和态密度计算中的问题
已经有8人回复
- 怎么用VASP计算声子谱
已经有4人回复
- VASP+phonopy INCAR 参数设定问题？
已经有12人回复
- 用VASP算MD时SMASS参数的设置？
已经有9人回复
- 第一性原理计算用MS好还是VASP好
已经有3人回复
- vasp测试SIGMA
已经有12人回复
- 受动力学因素影响很大的反应。。。。
已经有5人回复
- 动力学计算过程：时间步长、步数一般设置为多少？
已经有10人回复
- ELK 中的原子位置问题
已经有3人回复
- 用phonopy计算声子谱的问题
已经有16人回复
- 【求助】RWIGS和ISMEAR做态密度如何设
已经有7人回复
- 【其他】【诚聘版主】第一性原理版欢迎您的加入
已经有62人回复
- 【求助】castep计算中EOS拟合
已经有4人回复
- 【转帖】动力学蒙特卡洛方法(KMC)及相关讨论
已经有14人回复
- 【原创】elk中的状态方程拟合工具eos
已经有128人回复
- 【求助】询vasp价格
已经有12人回复
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
★ ★ ★ ★ ★
zxzj05(金币+5):谢谢回帖交流 2010-04-04 20:40
zxzj05(金币+5):谢谢回帖交流 2010-04-04 20:40
|
vasp分子动力学模拟。 转自：http://blog.163.com/ly0708@126/b ... 281201001181758620/ 转自：http://muchong.com/bbs/viewthread.php?tid=1744011&fpage=0&view=&highlight=&page=1 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 系综。 NBLOCK= 1 指的是多少个离子步（时间步长只是对于MD而言，对于其它的计算，NBLOCK也是起控制作用的）写一次XDATCAR KBLOCK=50 NBLOCK*KBLOCK 个离子步写一次PCDAT. CONTCAR是每个离子步之后都会写出来的，但是会用新的把老的覆盖 CHG是在每10个离子步写一次，不会覆盖 CHGCAR是在任务正常结束之后才写的。 VASP-md 模拟 INCAR 例子 转自：http://blog.163.com/xiaowei_0905 ... 183592010099379680/ This file determines the kind of job which VASP will perform; single point energy calculation (SPE), geometry optimisation (GO - coarse/fine), molecular dynamics (MD - nve/nvt), spin polarised calculation (mag). Examples can be found in /home/cs/model/vasp_util. Example; INCAR.spe $system = single point energy calc NELMIN = 4 minimum number of electronic SCF cycles EDIFF = 1E-6 stooping criterion for electronic convergence NSW = 0 number of ionic shifts ISMEAR = 0 treatment of partial occupancies of electronic levels Example; INCAR.coarse $system = coarse geom optimisation NELMIN = 4 EDIFF = 1E-2 EDIFFG = -1E-2 stopping criterion for forces Fmax < 0.01 eV/A IBRION = 2 minimisation method, good away from minimum ISIF = 3 optimise coords and cell pars LREAL =.TRUE. do calc in real space - quicker ISTART = 0 start with a random wavefunction NSW = 20 maximum of 20 ionic shifts ISMEAR = 0 LCHARG =.FALSE. don't write CHG and CHGCAR files Example; INCAR.fine $system = geom optimisation NELMIN = 4 EDIFF = 1E-6 EDIFFG = -1E-4 PREC = high increase energy cut-off by 25% IBRION = 1 minimisation method, good close to minimum ISIF = 3 NSW = 50 ISMEAR = 0 LCHARG=.FALSE. Example; INCAR.mag $system = collinear mag structure calc IBRION = 1 ISIF = 3 NPAR = 1 forces mag structure to be written in output file EDIFF = 1E-6 EDIFFG = -1E-3 PREC = high RWIGS = 1.376 0.900 1.233 1.302 radii for spherical integration of spin density, 1 per atom ISPIN = 2 do spin polarised calc MAGMOM = 24*0 5 -5 -5 5 initial mag moments for 28 atoms NSW = 20 Example; INCAR.nve $system = molecular dynamics ALGO = V MAXMIX = 40 IBRION = 0 do molecular dynamics NSW = 6000 number of time steps NBLOCK = 1 store structure every time step POTIM = 3.0 time step 3fs TEBEG = 673 target temperature ISYM = 0 turn off symmetry SMASS = -3 NVE ensemble LREAL =.TRUE. LCHARG =.FALSE. NELMIN = 4 PREC = LOW reduce energy cut-off by 25% for MD ISTART = 0 ISMEAR = 0; SIGMA=0.1 Example; INCAR.nvt $system = molecular dynamics ALGO = V MAXMIX = 40 IBRION = 0 NSW = 6000 NBLOCK = 1 POTIM = 3.0 TEBEG = 673 ISYM = 0 SMASS = 2 NVT ensemble, value determines frequency of coupling to heat bath LREAL =.TRUE. LCHARG =.FALSE. NELMIN = 4 PREC = LOW ISTART = 0 ISMEAR = 0; SIGMA=0.1 Example; INCAR.scale $system = molecular dynamics quench ALGO = V MAXMIX = 40 IBRION = 0 NSW = 50 NBLOCK = 5 rescale temperature every 5 steps POTIM = 3.0 TEBEG = 683 initial temp TEEND = 673 final temp ISYM = 0 SMASS = -1 MD with velocity scaling LREAL =.TRUE. LCHARG =.FALSE. NELMIN = 4 PREC = LOW ISTART = 0 ISMEAR = 0; SIGMA=0.1 |
8楼2010-04-04 10:58:56
freer007
木虫 (正式写手)
- 1ST强帖: 1
- 应助: 1 (幼儿园)
- 金币: 2421.7
- 红花: 1
- 帖子: 676
- 在线: 146小时
- 虫号: 707811
- 注册: 2009-02-24
- 性别: GG
- 专业: 半导体材料
4楼2010-04-03 20:28:15
cuilin0512
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 1859.8
- 散金: 1600
- 红花: 12
- 帖子: 1019
- 在线: 6778.9小时
- 虫号: 641563
- 注册: 2008-10-31
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
5楼2010-04-04 09:18:00
1357246
铁杆木虫 (职业作家)
壮士
- 应助: 4 (幼儿园)
- 贵宾: 0.226
- 金币: 5497.1
- 散金: 4194
- 红花: 30
- 沙发: 5
- 帖子: 3691
- 在线: 307.1小时
- 虫号: 856681
- 注册: 2009-09-25
- 性别: GG
- 专业: 基础物理学
6楼2010-04-04 09:38:15
投票:
30