24小时热门版块排行榜    

汕头大学海洋科学、生物学、生物与医药等3个专业接受调剂
查看: 14418  |  回复: 52

fanarsenal

木虫 (正式写手)

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

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

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

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

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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.

» 本帖已获得的红花(最新10朵)

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楼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。这样的结果可靠嘛?
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的关系,里面有啥物理意义吗?请教一下麻烦您了
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
这种情况,优化时也会出来这样的矩阵。
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
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 化工专硕310求调剂 +8 五毛冰淇淋a 2025-04-06 8/400 2025-04-08 10:46 by edmund7
[考研] 085701 环境工程,总分343,求调剂 +4 对山眠 2025-04-05 4/200 2025-04-08 08:36 by 1104338198
[考研] 数二英二,总分276,求化工导师收走 +12 张智翔11 2025-04-07 13/650 2025-04-07 23:21 by tbtin
[考研] 292求调剂 +12 2025417 2025-04-03 12/600 2025-04-07 16:54 by liandandan0
[考研] 304求调剂 +18 妮小汁 2025-04-06 20/1000 2025-04-07 15:43 by 天天爱读研
[考研] 找调剂学校 +5 爱吃土豆的大大� 2025-04-06 7/350 2025-04-07 11:24 by lijunpoly
[考研] 化学0703调剂293分 +8 太阳花向日葵 2025-04-04 8/400 2025-04-06 16:44 by calebby
[考研] 298求调剂 +16 研要上啊 2025-04-05 17/850 2025-04-06 16:32 by dick_runner
[考研] 296求调剂 +11 在结南桃 2025-04-05 12/600 2025-04-06 12:14 by qlm5820
[考研] 324求调剂 +3 无敌海皮皮 2025-04-05 3/150 2025-04-06 10:06 by brbrm
[考研] 一志愿中科大材料与化工338求调剂 +4 今年上岸ya 2025-04-05 4/200 2025-04-06 09:01 by JourneyLucky
[考研] 284求调剂 +4 变态小阿姨 2025-04-05 4/200 2025-04-05 20:04 by calebby
[考研] 材料考研求调剂 +4 wsl111 2025-04-04 4/200 2025-04-05 10:15 by 手札礼记
[考研] 河北大学调剂 +5 努力赚钱666 2025-04-05 6/300 2025-04-05 09:37 by 502874596
[考研] 生物技术与工程086001 318分求调剂 +9 pingjing0. 2025-04-03 9/450 2025-04-04 22:51 by 河大薇薇
[考研] 085500数一英一283求调剂 +3 忧郁的番茄! 2025-04-02 3/150 2025-04-03 23:36 by zhq0425
[考研] 085601找调剂院校 +7 forever.C 2025-04-03 8/400 2025-04-03 15:20 by 川哥v5
[考研] 求调剂 +11 寄于心 2025-04-02 13/650 2025-04-02 23:35 by siwuhen
[考研] 297求调剂 +14 宇不是雨 2025-04-01 15/750 2025-04-02 10:29 by panzt
[考研] 一志愿东南大学材料学硕,初试342分求调剂 +7 好想读研啊啊 2025-04-01 7/350 2025-04-02 09:00 by bingyu0729
信息提示
请填处理意见
婢跺秷鐦拫鍐ㄥ婢跺秷鐦拫鍐ㄥ
婢跺秷鐦拫鍐ㄥ娴滃瞼娣惍锟�閼惧嘲褰囨径宥堢槸 & 鐠嬪啫澧忛幐鍥у础閿涳拷瀵邦喕淇婇幍顐$閹碉拷
閸熷棗濮熼崥鍫滅稊閸熷棗濮熼崥鍫滅稊
閸熷棗濮熼崥鍫滅稊娴滃瞼娣惍锟�閸熷棗濮熼崥鍫滅稊閸溿劏顕�瀵邦喕淇婇幍顐$閹碉拷(婢跺洦鏁為敍姘櫌閸斺€虫値娴o拷)
鐎涳附婀抽獮鑼舵彛鐎涳附婀抽獮鑼舵彛
鐎涳附婀抽獮鑼舵彛娴滃瞼娣惍锟�瀵邦喕淇婇幍顐$閹碉拷娴滃棜袙閺堚偓閺傛澘顒熼張顖氬叡鐠愶拷
鐠佺儤鏋冩潏鍛嚤鐠佺儤鏋冩潏鍛嚤
鐠佺儤鏋冩潏鍛嚤娴滃瞼娣惍锟�瀵邦喕淇婇幍顐$閹碉拷濞h濮為懓浣哥瑎閿涘苯鍘ょ拹鐟版尒鐠囷拷
鐢牐绁幏娑樺珯鐢牐绁幏娑樺珯