应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » VASP5.2 直接计算弹性常数
531/6返回列表123456下一页
查看: 14016  |  回复: 52
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

fanarsenal

木虫 (正式写手)

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

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

» 收录本帖的淘帖专辑推荐

VASP VASP 第一性原理计算经验 纳米及第一性原理方向
VASP相关 第一性原理 仿真建模与计算

» 本帖已获得的红花（最新10朵）

mitong

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-05-20 20:53:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

doudouding

木虫 (小有名气)
引用回帖:
17楼: Originally posted by mitong at 2012-09-30 20:59:59
请问VASP应力应变法得到的弹性常数矩阵中的 XX          YY          ZZ          XY          YZ          ZX  是否对应着能量法计算出的1 2 3 4 5 6 也就是说应力应变法得到的矩阵对角线上的数是否和能量法的相对 ...

XX -> 1  YY->2  ZZ->3  XY->6  YZ->4  ZX->5
注意不是按照123456的顺序排列的。可以试试计算六角的情况，应该满足c66=(c11-c12)/2.
赞一下(4人) 回复此楼

» 本帖已获得的红花（最新10朵）

Vaucanson
22楼2012-12-04 11:13:06
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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人) 回复此楼
宁愿骑着自行车笑，不愿坐着宝马哭。。。。。
2楼2012-05-22 09:55:48
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ufo991

铁虫 (正式写手)
引用回帖:
4楼: Originally posted by condensed at 2012-05-24 01:23:37
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
...

NFREE=2是不是也可以？还有在计算时候它会自动找对称性，在计算时候出现K点不匹配VERY BAD NEWS! internal error in subroutine IBZKPT:Reciprocal lattice and k-lattice belong to different class of lattices. Often results are still useful...  96，这个对结果会有影响吗？这种方法算弹性常数，准确率高不，我试着计算了一下Al的体弹模量结果为56GPa,实验值是76Gpa。这样的结果可靠嘛？
赞一下(2人) 回复此楼
6楼2012-05-29 10:52:13
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

随便注册

银虫 (初入文坛)
引用回帖:
11楼: Originally posted by liliangfang at 2012-08-27 21:01:09
请问这结果怎样处理恩，谢谢...

TOTAL ELASTIC MODULI  --- 1
ELASTIC MODULI CONTR FROM IONIC RELAXATION --- 2
SYMMETRIZED SYMMETRIZED --- 3
是不是读取第一个？而且是不是1=2+3的关系，里面有啥物理意义吗？请教一下麻烦您了
赞一下(1人) 回复此楼
25楼2013-03-29 08:57:00
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

janet2012222

金虫 (小有名气)
引用回帖:
4楼: Originally posted by condensed at 2012-05-24 01:23:37
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
...

请问，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
这种情况，优化时也会出来这样的矩阵。
赞一下(1人) 回复此楼
26楼2013-05-30 21:44:40
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

李玥玥

铜虫 (正式写手)
引用回帖:
47楼: Originally posted by 不要烦恼 at 2017-04-14 10:29:22
应该是YZ->5,ZX->4吧...

是yz->4,  xz->5
赞一下(1人) 回复此楼
49楼2018-03-09 14:57:14
 已阅   回复此楼   关注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的回帖

fanarsenal

木虫 (正式写手)
谢谢诸位……
回复此楼
5楼2012-05-25 16:57:52
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lbing9002

木虫 (正式写手)
学习了
回复此楼
穷则独善其身，达则兼济天下。
7楼2012-05-29 21:47:27
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

木虫 (著名写手)
引用回帖:
4楼: Originally posted by condensed at 2012-05-24 01:23:37
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
...

你好，请问在这步之前是不是要先进行结构优化啊，谢谢！
赞一下 回复此楼
未来就是现在
8楼2012-08-24 08:54:04
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

future_wl

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

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

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

你好，请问NFREE这个参数应该如何设置，比如我有一个16原子的体系，这个值应该取多少呢？多谢！
赞一下 回复此楼
未来就是现在
9楼2012-08-24 14:34:34
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fanarsenal

木虫 (正式写手)
引用回帖:
8楼: Originally posted by future_wl at 2012-08-24 08:54:04
你好，请问在这步之前是不是要先进行结构优化啊，谢谢！...

肯定是进行结构优化之后的结果再去计算其弹性常数啊……
赞一下 回复此楼
10楼2012-08-25 16:21:55
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fanarsenal 的主题更新
531/6返回列表123456下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[微米和纳米] 【怎么确定一种物质呀】xrd测不出峰的材料 50+3 小裴想上学 2024-12-21 5/250 2024-12-23 14:23 by lyxy0217
[有机交流] 装柱干法上样 +5 刘松垚 2024-12-20 5/250 2024-12-23 13:14 by 风永不止步
[教师之家] 咨询一下，是不是教授之间比较忌讳谈论文造假，因为可能很多人都是造假 +8 akslis2024 2024-12-21 9/450 2024-12-23 12:39 by jurkat.1640
[教师之家] 为啥很多高校老师不愿意当副院长？ +12 意得辑_editage 2024-12-18 16/800 2024-12-23 11:29 by jurkat.1640
[论文投稿] 三个手稿 10+4 s爱吃芒果 2024-12-22 8/400 2024-12-23 11:15 by bobvan
[找工作] 化学博士 +6 实名心动国民萌� 2024-12-17 8/400 2024-12-23 09:33 by zhonghu
[硕博家园] 四线城市公办高校无编or二线城市民办高校 +14 joshhartford 2024-12-20 32/1600 2024-12-23 08:49 by 卓尔me
[考博] 定向就业，毕业证学位证由工作单位领取？ +3 河西夜郎 2024-12-22 4/200 2024-12-22 20:30 by changlu999
[考博] 中山医学院生物信息申博 +3 wcy98117 2024-12-22 3/150 2024-12-22 13:42 by k1f01@4s
[论文投稿] 投稿意见求助，没弄清回答的方向 8+4 moonlig 2024-12-18 4/200 2024-12-22 08:51 by steven_198377
[硕博家园] 博一 +10 幽皮皮皮 2024-12-17 16/800 2024-12-22 07:37 by Ermito
[教师之家] 今年的科研绩效只奖励Q1了 +12 holypower 2024-12-20 17/850 2024-12-21 16:13 by willbuilder
[教师之家] 学校排名靠前发展好，最大的受益者是校领导还是普通老师？ +6 akslis2024 2024-12-17 7/350 2024-12-21 11:22 by 凌晨一点393
[有机交流] 工艺需求，付费 +4 想当然灬 2024-12-17 6/300 2024-12-21 06:23 by a893069456
[教师之家] 删了哦 (EPI+-1)(金币-50) +5 原因在哪里 2024-12-18 6/300 2024-12-19 12:25 by keikofans
[考博] 有未就业的博士吗 +4 大发财树 2024-12-17 4/200 2024-12-18 14:13 by 她在学习
[有机交流] 求这个反应的反应条件 +6 戒骄戒躁+ 2024-12-17 6/300 2024-12-18 11:43 by diaodeng
[有机交流] 细丝状的晶体可以测单晶吗？ +5 好气好气 2024-12-16 5/250 2024-12-17 22:53 by 请让我滚去学习
[职场人生] 网络游戏有意义吗 +3 神明.斯派克. 2024-12-16 3/150 2024-12-17 22:13 by zyqchem
[公派出国] 意大利帕多瓦大学招收CSC博士生--光电催化材料在新能源及可持续绿色化学方向 +3 huangyuxing2 2024-12-16 3/150 2024-12-17 18:25 by zhizhan127
信息提示
关闭
请填处理意见
关闭