| 查看: 3291 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
babaleo木虫 (正式写手)
|
[交流]
【求助】非零压下弹性常数计算的初始结构
|
||
|
根据侯博士指南，通过系列a值以及ISIF=5 （或4）高精度优化得到六角系Mg 的E-V关系，然后用B-M拟合得到了V0,E0, B,B'. 再然后每个a值的压力便知道了。 请问各位，这些状态方程拟合用的结构可以用来作为非零压下弹性常数初始结构吗？也就是说，这些结构可以拷成侯博士说的 OLDPOS 去做弹性常数计算吗？如果可以，是不是可以把这过程和弹性常数计算一起脚本化？试了一下，发现某些应变下出现能量单调下降现象， 不知如何改进，恳请大虾指教。 |
回复此楼» 猜你喜欢
能否申诉？
已经有8人回复
-
梦想
已经有5人回复
-
中青基了要发朋友圈吗？
已经有7人回复
-
有多少人是今天查系统知道结果的？
已经有19人回复
-
29号明天会评吗
已经有4人回复
-
国社科又开始会评了，不知道这次命运如何
已经有10人回复
-
基金不中，共勉
已经有12人回复
-
面上意见出来了
已经有19人回复
-
为什么到现在没收到通知？
已经有5人回复
-
有没有仍没收到信息的
已经有8人回复
sunyang1988
至尊木虫 (著名写手)
- 1ST强帖: 1
- 应助: 24 (小学生)
- 贵宾: 0.984
- 金币: 9397.7
- 散金: 1593
- 红花: 23
- 沙发: 1
- 帖子: 1744
- 在线: 1255.8小时
- 虫号: 965316
- 注册: 2010-03-09
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
6楼2010-09-26 01:03:58
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
30