| 查看: 7958 | 回复: 25 | ||||||
| 【有奖交流】积极回复本帖子，参与交流，就有机会分得作者 magiczerg 的 3 个金币 | ||||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||||
magiczerg金虫 (小有名气)
|
[交流]
【求助】侯老师的vasp入门的晶格优化部分没看懂，问一下
|
|||||
|
复杂的情况,以六角结构Mg的晶格常数为例 1.先进行一次体积保持不变的粒子驰豫的计算（ISIF=5） 2再把优化得到的CONTCAR拷贝成POSCAR，进行一次静态的计算 3对a的值取10个左右的点，每个点重复上面两步，得到静态计算下volume-Etot的关系。 这三部的脚步可以通过以下脚步run_cell来进行，内容如下： #!/bin/sh rm WAVECAR for i in 5.20 5.23 5.26 5.29 5.32 5.35 5.38 5.41 5.44 5.47 do cat > POSCAR << ! ureadopef $i 0 -1 0 0.8660254 0.5 0 0 0 1.6230530 2 Direct 0.666666667 0.333333333 0.75 0.333333333 0.666666667 0.25 ! cp INCAR.relax INCAR echo "a=$i angestrom"; mpiexec -np 8 vasp cp CONTCAR POSCAR cp INCAR.static INCAR echo "a= $i angstrom ";mpiexec -np 8 vasp E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'` V=`grep "volume" OUTCAR | tail -1 | awk '{printf "%12.4f \n", $5}'` echo $V $E >>EtVo.dat done 4.采用状态方程拟合得到的平衡态下的体积，体弹性模量。 5.在该体积下，重复上面第2步和第3步，得到平衡状态下的其他晶胞参数。这一步也就是：在的到了E(V)曲线后，通过状态方程拟合得到平衡状态下的体积，计算出上面脚本中变量$i的值（弱弱的问一句，晶胞形状都改变了，怎么算a的值啊？就算还是按六角结构来算，a和c的关系不知道，怎么算啊，还是按照初始时a/c的比值来算吗？ （$i 0 -1 0 0.8660254 0.5 0 0 0 1.6230530）） 并改变$i的循环值，再运行run_cell计算一次，得到其他的结构参数c和位置u。 （弱弱的再问：改变$i的循环值不还是在算a吗？算a和算c有什么区别吗？另外算位置u是不是固定晶胞形状和体积只改变原子位置就行了，改一下INCAR就可以了，没必要在运行runcell了吧？） [ Last edited by magiczerg on 2010-8-14 at 20:19 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 仿真建模与计算 | 计算化学-辅助工具 | 科研 |
» 猜你喜欢
你的学校职称晋升认CSCD吗？
已经有3人回复
-
请问大家，我有希望进浙大普博么
已经有8人回复
-
applied catalysis B 投稿模版
已经有5人回复
-
希望有师兄师姐帮忙引荐申博
已经有11人回复
-
毕业论文的数据能否再发小论文
已经有7人回复
-
东南大学电气学院电力电子方向招收硕士博士
已经有7人回复
-
开心！国庆前发个小财~
已经有15人回复
-
QB怎么还没消息
已经有57人回复
-
求推荐医学类杂志
已经有4人回复
-
25有机合成申博
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 急求这一步的机理及详细讲解。谢谢！
已经有8人回复
- E-V曲线拟合后怎么找到能量最低点对应的V值
已经有4人回复
- 晶化原理及容器注意事项等
已经有5人回复
- 请教如何用VASP计算氧化物的晶格常数？以alpha-Al2O3为例
已经有7人回复
- 求助 法语入门。。
已经有5人回复
- 关于六角晶系的弹性常数计算
已经有16人回复
- 几何优化后的晶格参数在哪里查看呢？
已经有4人回复
- VASP 计算E-V曲线求助。。
已经有10人回复
- TED有个视频讲的打印人的器官，没看懂？
已经有20人回复
- 一次型NiAl结构优化，晶格矢量变化非常大，求助~
已经有5人回复
- 求助VASP计算弹性常数的问题
已经有5人回复
- 请教一下那款免费软件的结构优化功能比较好用？
已经有7人回复
- 求助vasp_晶格参数扫描
已经有5人回复
- 真不敢相信，江西师范侯豪情老师将静电纺丝产业化！！！！
已经有28人回复
- VASP的源代码是公开的吗？还有其他可以读懂的程序吗？
已经有12人回复
- 弱弱问一下，characters在这句话里指字母还是单词？
已经有11人回复
- 武汉药企还是湖南衡阳高职老师？
已经有6人回复
- 【求助】vasp优化时选择了isif=3 怎么限制晶格的一个维度呢？
已经有16人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】谁能给解释一下VASP手册中wrap around errors什么意思？
已经有3人回复
- 【求助】晶格常数大了以后 VASP的计算值准确吗
已经有5人回复
- 【求助】vasp优化晶格参数a
已经有13人回复
- 【求助】简单问题----关于vasp晶格优化
已经有7人回复
- 【求助成功】晶格优化的结果很诡异
已经有13人回复
- 【求助】拟合晶格常数
已经有14人回复
- 【分享】教你看懂化验单（图片不能显示问题已解决）【已搜索无重复】
已经有23人回复
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25682.7
- 散金: 602
- 红花: 166
- 帖子: 3778
- 在线: 873.7小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
7楼2010-08-16 18:11:45
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25682.7
- 散金: 602
- 红花: 166
- 帖子: 3778
- 在线: 873.7小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
小木虫: 金币+0.5, 给个红包，谢谢回帖
|
vasp只能算振动频率，其他的比如平动和转动频率不能计算。 计算频率的方法有两种，一种是直接的，基于finite dispacement，另外一种是基于密度泛函微扰论(或线性响应的方法)。 对于第一种很简单。看一下手册的这里就可以：http://cms.mpi.univie.ac.at/vasp/vasp/IBRION_5_IBRION_6.html |
14楼2013-12-27 09:36:02
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25682.7
- 散金: 602
- 红花: 166
- 帖子: 3778
- 在线: 873.7小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
20楼2014-06-24 19:44:53
10