| 查看: 4145 | 回复: 7 | ||||
天棚元帅木虫 (著名写手)
|
[求助]
请教如何用VASP计算氧化物的晶格常数？以alpha-Al2O3为例
|
|
如题，请教各位虫友， 计算fcc金属（a = b =c），用murg script（Birch Murnaghan equation of state）拟合可以得到优化好的lattice constant （a）。 而金属氧化物，以alpha-Al2O3为例，a = b ≠ c, 如何拟合其晶格常数？（a ？c？） 希望有知道的虫友不吝赐教哦~~~  ~~~最好是能结合我想计算的alpha-Al2O3的晶格常数（a，c）来给予尽可能详细的说明~~~非常感谢~~~------------------------- 我参考了论坛中以下几个帖子，也不是很明白，其中的给出的帮助链接也大都失效。。  。。【原创】elk中的状态方程拟合工具eos http://muchong.com/bbs/viewthread.php?tid=2056902&fpage=1 【求助】E-V曲线晶体状态方程拟合问题 http://muchong.com/bbs/viewthread.php?tid=2755021 【求助】利用VASP算出E-V曲线后，怎么拟合得到体弹性模量？ http://muchong.com/bbs/viewthread.php?tid=2069233 【求助成功】测试晶格常数后用什么方程来拟合？ http://muchong.com/bbs/viewthread.php?tid=2026674&page=1 【求助】侯老师的vasp入门的晶格优化部分没看懂，问一下 http://muchong.com/bbs/viewthread.php?tid=2299924  [ Last edited by 天棚元帅 on 2013-2-17 at 10:17 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 Photochemistry |
» 猜你喜欢
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
麻烦专家们看看评委们的意见（F口面上）
已经有10人回复
-
面上意见出来了
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有9人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有9人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教我的用vasp优化Ni的晶格常数的时候怎么不收敛啊。内含INCAR
已经有10人回复
- 【求助】关于VASP计算材料的晶格热导率
已经有8人回复
- VASP一次性优化锐钛矿TiO2晶格常数时出错！！
已经有6人回复
- 【求助】vasp优化时选择了isif=3 怎么限制晶格的一个维度呢？
已经有16人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】晶格常数大了以后 VASP的计算值准确吗
已经有5人回复
- 【求助】vasp优化晶格参数a
已经有13人回复
- 【求助】侯老师的vasp入门的晶格优化部分没看懂，问一下
已经有25人回复
- 【转帖】如何用VASP计算晶格常数
已经有8人回复
天棚元帅
木虫 (著名写手)
- 应助: 100 (初中生)
- 金币: 3141.3
- 散金: 625
- 红花: 3
- 帖子: 1186
- 在线: 283.9小时
- 虫号: 609050
- 注册: 2008-09-21
- 性别: GG
- 专业: 理论和计算化学
2楼2013-02-17 10:27:30
guohuazhon
木虫 (正式写手)
- 应助: 5 (幼儿园)
- 金币: 12479.2
- 散金: 10
- 红花: 2
- 帖子: 931
- 在线: 92.6小时
- 虫号: 807030
- 注册: 2009-07-11
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
3楼2013-02-17 10:53:30
天棚元帅
木虫 (著名写手)
- 应助: 100 (初中生)
- 金币: 3141.3
- 散金: 625
- 红花: 3
- 帖子: 1186
- 在线: 283.9小时
- 虫号: 609050
- 注册: 2008-09-21
- 性别: GG
- 专业: 理论和计算化学
|
追加问题啊~~ ================================== 【侯老师指南上部分内容摘录】 复杂的情况,以六角结构Mg的晶格常数为例 1.以实验的晶格参数为基础做好POSCAR，等 2. 先进行一次体积保持不变的粒子驰豫的计算（ISIF=5） 3 再把优化得到的CONTCAR拷贝成POSCAR，进行一次静态的计算 4 对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 5.采用状态方程拟合得到的平衡态下的体积，体弹性模量。 6.在该体积下，重复上面第2步和第3步，得到平衡状态下的其他晶胞参数。这一步也就是：在的到了E(V)曲线后，通过状态方程拟合得到平衡状态下的体积，计算出上面脚本中变量$i的值，并改变$i的循环值，再运行run_cell计算一次，得到其他的结构参数c和位置u。 ======================================== 我是准备计算出Al2O3的晶格参数后，切出表面，作表面的吸附及反应等计算模拟。 仔细看完侯老师的指南，还是有些地方没明白， 1）在我这种情况下，在得到平衡体积以及相应的晶格常数a后，另一个晶格常数c直接不就可以用V/a方算出来了么？为什么还要再重头对晶格常数c算一遍呢？ 2）在我的这种情况下，在进行体积保持不变的原子弛豫计算中（侯老师指南），INCAR的ISIF值应该设置为多少呀？4？ 3）如何保证在进行体积保持不变的原子弛豫计算中（侯老师指南），a/b值为1呢？还是不用管这个问题呢？ 麻烦高手来解答下~~非常感谢~~重金酬谢啊~~ |
4楼2013-02-17 14:29:06
5楼2013-02-19 16:27:07
天棚元帅
木虫 (著名写手)
- 应助: 100 (初中生)
- 金币: 3141.3
- 散金: 625
- 红花: 3
- 帖子: 1186
- 在线: 283.9小时
- 虫号: 609050
- 注册: 2008-09-21
- 性别: GG
- 专业: 理论和计算化学
★ ★
franch: 金币+2, 鼓励交流，， 2013-03-11 23:16:59
franch: 金币+2, 鼓励交流，， 2013-03-11 23:16:59
|
问题已于前段时间解决，现来回复下： 解决方案基本采用如下贴中二楼虫友的回复， http://muchong.com/bbs/viewthread.php?tid=2299924 我的计算结果比实验值超了大概1%，还算合理。 |
6楼2013-03-11 11:04:43
Brilliantly
捐助贵宾 (正式写手)
- 应助: 7 (幼儿园)
- 金币: 420
- 散金: 2788
- 红花: 8
- 帖子: 768
- 在线: 312.2小时
- 虫号: 1360339
- 注册: 2011-08-03
- 性别: GG
- 专业: 功能陶瓷
7楼2014-05-06 20:55:53
天棚元帅
木虫 (著名写手)
- 应助: 100 (初中生)
- 金币: 3141.3
- 散金: 625
- 红花: 3
- 帖子: 1186
- 在线: 283.9小时
- 虫号: 609050
- 注册: 2008-09-21
- 性别: GG
- 专业: 理论和计算化学
8楼2014-05-08 04:08:58
