| 查看: 3260 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
babaleo木虫 (正式写手)
|
[交流]
【求助】非零压下弹性常数计算的初始结构
|
||
|
根据侯博士指南，通过系列a值以及ISIF=5 （或4）高精度优化得到六角系Mg 的E-V关系，然后用B-M拟合得到了V0,E0, B,B'. 再然后每个a值的压力便知道了。 请问各位，这些状态方程拟合用的结构可以用来作为非零压下弹性常数初始结构吗？也就是说，这些结构可以拷成侯博士说的 OLDPOS 去做弹性常数计算吗？如果可以，是不是可以把这过程和弹性常数计算一起脚本化？试了一下，发现某些应变下出现能量单调下降现象， 不知如何改进，恳请大虾指教。 |
回复此楼» 猜你喜欢
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有7人回复
-
今天基金会出结果吗？20260819
已经有18人回复
-
我面上完蛋了
已经有14人回复
-
范进中举一文的中心思想
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有27人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
-
没有任何消息-是不是就凉了
已经有8人回复
-
人气不行了
已经有10人回复
-
能否退出参与的面上项目解除限项
已经有24人回复
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
- 专业: 人工晶体
|
补充以下, 下面是我的脚本,请大虾帮看看. #!/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
babaleo
木虫 (正式写手)
- 应助: 18 (小学生)
- 金币: 2391.1
- 散金: 40
- 红花: 1
- 帖子: 375
- 在线: 477.1小时
- 虫号: 981394
- 注册: 2010-03-25
- 专业: 人工晶体
5楼2010-09-25 22:52:49
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
20