24小时热门版块排行榜    

查看: 16213  |  回复: 52
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

fanarsenal

木虫 (正式写手)

[求助] VASP5.2 直接计算弹性常数

需要直接用应力应变关系得到弹性常数,看Manul上5.2可以做到这一点。可是试了几次都有问题,有用过的能不能给点儿提示。
除了在INCAR里面加上ISIF = 3 IBRION = 6之外,还要怎么设置?
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

janet2012222

金虫 (小有名气)

引用回帖:
27楼: Originally posted by fanarsenal at 2013-05-31 09:26:54
你找的内容不对,你可以在OUTCAR中找一下kBar...

您好!我看见OUTCAR出来这样的话 INCAR:
POTCAR:    PAW_PBE Re 17Jan2003                  
POTCAR:   PAW_PBE B 06Sep2000                    

-----------------------------------------------------------------------------
|                                                                             |
|           W    W    AA    RRRRR   N    N  II  N    N   GGGG   !!!           |
|           W    W   A  A   R    R  NN   N  II  NN   N  G    G  !!!           |
|           W    W  A    A  R    R  N N  N  II  N N  N  G       !!!           |
|           W WW W  AAAAAA  RRRRR   N  N N  II  N  N N  G  GGG   !            |
|           WW  WW  A    A  R   R   N   NN  II  N   NN  G    G                |
|           W    W  A    A  R    R  N    N  II  N    N   GGGG   !!!           |
|                                                                             |
|      For optimal performance we recommend that you set                      |
|        NPAR = approx SQRT( number of cores)                                 |
|      This will greatly improve the performance of VASP for DFT.             |
|      The default NPAR=number of cores might be grossly inefficient          |
|      on modern multi-core architectures or massively parallel machines.     |
|      Unfortunately you need to use the default for hybrid, GW and RPA       |
|      calculations.                                                          |
|                                                                             |
-----------------------------------------------------------------------------
OUTCAR 最后也出来|      VASP internal routines  have requested a change of the k-point set.    |
|      Unfortunately this is only possible if NPAR=number of nodes.           |
|      Please remove the tag NPAR from the INCAR file and restart the         |
|      calculations.                                                          |
|                                                                             |
我把INCAR里面的NPAR删掉了也会出来这样的话。请问这是怎么回事呢
29楼2013-05-31 10:27:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 53 个回答

zsl_321123

铁虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
liliangfang: 金币+2, 谢谢指教 2012-05-22 10:02:54
fanarsenal: 金币+2, ★★★很有帮助 2012-05-25 16:57:39
vasp5.2计算可以直接计算弹性常数。

关键参数:  将IBRION=6;ISIF=3;NFREE=4;

弹性矩阵会直接写在OUTCAR里。
宁愿骑着自行车笑,不愿坐着宝马哭。。。。。
2楼2012-05-22 09:55:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)

引用回帖:
2楼: Originally posted by zsl_321123 at 2012-05-22 09:55:48:
vasp5.2计算可以直接计算弹性常数。

关键参数:  将IBRION=6;ISIF=3;NFREE=4;

弹性矩阵会直接写在OUTCAR里。

其他的设置会不会跟这个有冲突的说,我增加了这三项,为保险起见还增加了NSW=1,结果却是没法从INCAR中读取ISTART,#掉ISTART之后,又会读不出下面的,不知为什么……
3楼2012-05-22 17:38:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

condensed

木虫 (著名写手)

【答案】应助回帖

★ ★ ★
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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 2JOx3r2CYEgw 2026-08-22 4/200 2026-08-23 03:55 by OEbVnUOu01ol
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 2JOx3r2CYEgw 2026-08-22 7/350 2026-08-23 03:55 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 2JOx3r2CYEgw 2026-08-22 6/300 2026-08-23 03:40 by OEbVnUOu01ol
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QTy3jDtz1uLt 2026-08-21 9/450 2026-08-23 03:04 by OEbVnUOu01ol
[基金申请] 2026国自然函评费到账 +11 羊腰板 2026-08-21 11/550 2026-08-22 22:03 by lsbin3733
[基金申请] 范进中举一文的中心思想 +3 zjuhero 2026-08-22 4/200 2026-08-22 20:02 by mycsru
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 2JOx3r2CYEgw 2026-08-21 4/200 2026-08-22 16:52 by sunzitan
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 今日不放榜?网传国自然预计 8 月 27 日可查结果 +16 医学老男孩 2026-08-20 20/1000 2026-08-21 21:13 by Ldrop2023
[基金申请] 科研孤儿太难了 +17 我4大白菜 2026-08-20 18/900 2026-08-21 20:57 by zhangev
[基金申请] 看来今天不会放榜了? +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 让我中一个面上吧! +11 大萍1987 2026-08-20 13/650 2026-08-21 14:48 by 20120902066
[基金申请] 感觉是下周放榜了 +7 angus9576 2026-08-17 12/600 2026-08-21 13:38 by weiyin
[论文投稿] 投稿咨询 +5 wwm09 2026-08-17 7/350 2026-08-21 10:11 by 期刊论文帮手
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 估计是周四 +3 archvillain 2026-08-18 3/150 2026-08-21 01:48 by jnhyjjm
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
信息提示
请填处理意见