| 查看: 3284 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
babaleo木虫 (正式写手)
|
[交流]
【求助】非零压下弹性常数计算的初始结构
|
||
|
根据侯博士指南，通过系列a值以及ISIF=5 （或4）高精度优化得到六角系Mg 的E-V关系，然后用B-M拟合得到了V0,E0, B,B'. 再然后每个a值的压力便知道了。 请问各位，这些状态方程拟合用的结构可以用来作为非零压下弹性常数初始结构吗？也就是说，这些结构可以拷成侯博士说的 OLDPOS 去做弹性常数计算吗？如果可以，是不是可以把这过程和弹性常数计算一起脚本化？试了一下，发现某些应变下出现能量单调下降现象， 不知如何改进，恳请大虾指教。 |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有3人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有5人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有5人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有7人回复
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
7楼2010-09-27 01:37:13
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
|
补充以下, 下面是我的脚本,请大虾帮看看. #!/bin/sh rm WAVECAR for i in 1.90 1.95 2.00 2.05 2.10 2.15 2.20 2.30 do cat > POSCAR < Mg 1 $i 0.0 -1.0 0.0 0.8660254037844 0.5 0.0 0.0 0.0 1.580 2 Direct 0.6666666666666667 0.3333333333333333 0.750 0.3333333333333333 0.6666666666666667 0.250 ! cp INCAR.relax INCAR mpirun -np 4 vasp cp CONTCAR OLDPOS echo "running C11+C12 ...." ./geo-opt-dd0000 echo "Running C11-C12...." ./geo-opt-00000d echo "Running C33...." ./geo-opt-00d000 echo "Running C44...." ./geo-opt-123dd6 echo "Running bulk...." ./geo-opt-ddd456 done |
2楼2010-09-25 10:34:07
sunyang1988
至尊木虫 (著名写手)
- 1ST强帖: 1
- 应助: 24 (小学生)
- 贵宾: 0.984
- 金币: 9397.7
- 散金: 1593
- 红花: 23
- 沙发: 1
- 帖子: 1744
- 在线: 1255.8小时
- 虫号: 965316
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
4楼2010-09-25 17:55:08
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
5楼2010-09-25 22:52:49
