| 查看: 15202 | 回复: 52 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
fanarsenal木虫 (正式写手)
|
[求助]
VASP5.2 直接计算弹性常数
|
||
|
需要直接用应力应变关系得到弹性常数，看Manul上5.2可以做到这一点。可是试了几次都有问题，有用过的能不能给点儿提示。 除了在INCAR里面加上ISIF = 3 IBRION = 6之外，还要怎么设置？ |
回复此楼» 本帖已获得的红花（最新10朵）
mitong» 猜你喜欢
青椒八年已不青，大家都被折磨成啥样了？
已经有4人回复
-
救命帖
已经有8人回复
-
招博士
已经有4人回复
-
限项规定
已经有6人回复
-
西南交通大学国家级人才团队2026年博士研究生招生（考核制）—机械、材料、力学方向
已经有3人回复
-
英文综述是否需要润色及查重
已经有5人回复
-
为什么nbs上溴 没有产物点出现呢
已经有9人回复
-
最失望的一年
已经有18人回复
janet2012222
金虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 416.8
- 散金: 630
- 红花: 3
- 帖子: 294
- 在线: 181.4小时
- 虫号: 1633598
- 注册: 2012-02-22
- 专业: 碳素材料与超硬材料
|
请问，VASP.5.2.12可以直接用这种方法计算弹性常数吗。我的INCAR按照您这种方法设置，却没有出现这样的矩阵。出来的只有Direction XX YY ZZ XY YZ ZX -------------------------------------------------------------------------------------- Alpha Z 168.35185 168.35185 168.35185 Ewald -438.11662 -438.11665 -333.25913 0.00000 0.00000 -0.00000 Hartree 3.07456 3.07456 22.22788 -0.00000 -0.00000 -0.00000 E(xc) -102.35527 -102.35527 -103.34240 0.00000 0.00000 0.00000 Local 62.67338 62.67339 -51.50159 -0.00000 0.00000 0.00000 n-local -42.23117 -45.57686 -37.94608 -3.45887 -0.23611 -0.25461 augment 19.13534 19.13534 18.25597 -0.00000 -0.00000 -0.00000 Kinetic 329.90385 332.13257 317.41678 -7.90017 1.01621 1.72256 Fock 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 这种情况，优化时也会出来这样的矩阵。 |
26楼2013-05-30 21:44:40
2楼2012-05-22 09:55:48
fanarsenal
木虫 (正式写手)
- 应助: 47 (小学生)
- 金币: 2021.7
- 红花: 23
- 帖子: 515
- 在线: 659.4小时
- 虫号: 1069215
- 注册: 2010-08-04
- 专业: 金属结构材料
3楼2012-05-22 17:38:26
condensed
木虫 (著名写手)
- 应助: 21 (小学生)
- 金币: 3223.5
- 散金: 44
- 红花: 3
- 帖子: 1002
- 在线: 208.4小时
- 虫号: 292567
- 注册: 2006-11-04
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★
fanarsenal: 金币+3, ★★★很有帮助 2012-05-25 16:57:46
fanarsenal: 金币+3, ★★★很有帮助 2012-05-25 16:57:46
|
Diamond: INCAR: SYSTEM = diamond PREC = Accurate ENCUT = 600.0 EDIFF = 1e-6 ISMEAR = 0 SIGMA = 0.05 POTIM = 0.100 LCHARG = FALSE LWAVE = FALSE IBRION=6 ISIF=3 NFREE=4 EDIFFG = 1e-5 EDIFFG = -0.001 OUTCAR 92772 TOTAL ELASTIC MODULI (kBar) 92773 Direction XX YY ZZ XY YZ ZX 92774 -------------------------------------------------------------------------------- 92775 XX 10542.0731 1253.9340 1253.8826 -33.2003 0.0000 0.0000 92776 YY 1253.9340 10542.0731 1253.8826 -33.2003 0.0000 0.0000 92777 ZZ 1253.8826 1253.8826 10542.0889 -4.3655 0.0000 0.0000 92778 XY -33.2003 -33.2003 -4.3655 5633.1588 0.0000 0.0000 92779 YZ 0.0000 0.0000 0.0000 0.0000 5633.9957 -11.1251 92780 ZX 0.0000 0.0000 0.0000 0.0000 -11.1251 5633.9957 92781 -------------------------------------------------------------------------------- OK |
4楼2012-05-24 01:23:37