| 查看: 3270 | 回复: 9 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
babaleo木虫 (正式写手)
|
[交流]
【求助】非零压下弹性常数计算的初始结构
|
||
|
根据侯博士指南，通过系列a值以及ISIF=5 （或4）高精度优化得到六角系Mg 的E-V关系，然后用B-M拟合得到了V0,E0, B,B'. 再然后每个a值的压力便知道了。 请问各位，这些状态方程拟合用的结构可以用来作为非零压下弹性常数初始结构吗？也就是说，这些结构可以拷成侯博士说的 OLDPOS 去做弹性常数计算吗？如果可以，是不是可以把这过程和弹性常数计算一起脚本化？试了一下，发现某些应变下出现能量单调下降现象， 不知如何改进，恳请大虾指教。 |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
国自然面上复盘~欢迎讨论
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
基金不中，共勉
已经有5人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
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
