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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于用VASP计算弹性常数的两个问题
51/1返回列表
查看: 5452  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

Vaucanson

金虫 (小有名气)

[求助] 关于用VASP计算弹性常数的两个问题 已有2人参与

最近在学习用VASP 5.3.5计算弹性常数，设置IBRION = 6, ISIF = 4计算。发现有两个问题想请教大家：
1. 当我用16核计算时，设置NCORE = 4，走完第一个离子步后报错，停止计算：
关于用VASP计算弹性常数的两个问题
这是否意味着要直接取消掉NCORE一类的相关参数？取消之后是warning，提示我设置NCORE。

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

不知道是否对结果有什么影响？

谢谢大家了。
回复此楼

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

弹性计算

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2015-05-17 16:23:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vasp001

金虫 (正式写手)
引用回帖:
2楼: Originally posted by awdawdawdd at 2015-05-17 18:49:02
IBRION=6的时候 有关parallel的开关只能选择默认值（NPAR NCORE之类的）

我觉得弹性常数最好还是手动算，尤其是对于对称性不高的情况。通过施加strain，每一个类型的strain多弄几个点，在ISIF=2下固定cell，rela ...

你好，我最近使用VASP计算D3h晶体的弹性常数。然后输出的结果不是很懂。
不知道如何从VASP的OUTCAR中获得 Cij, 请多指教。

TOTAL ELASTIC MODULI (kBar)
Direction    XX          YY          ZZ          XY          YZ          ZX
--------------------------------------------------------------------------------
XX         614.9821    162.6485     -1.1175     -0.5210      0.0000      0.0000
YY         162.6485    615.4451     -1.0843      0.9220      0.0000      0.0000
ZZ          -1.1175     -1.0843     -1.7237      0.0287      0.0000      0.0000
XY          -0.5210      0.9220      0.0287    227.1156      0.0000      0.0000
YZ           0.0000      0.0000      0.0000      0.0000     -0.4616      0.0064
ZX           0.0000      0.0000      0.0000      0.0000      0.0064     -0.4690

另外如何根据晶体对称性来确定Cij中独立变量，
如何从VASP的OUTCAR文件中获得Cij数值？


还请不吝指教。

另外你有即时沟通工具吗？ 有的话，方便提供一下吗？
赞一下 回复此楼
一定要精通MaterialsStudios,VASP.
8楼2015-06-14 16:04:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

awdawdawdd

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Vaucanson: 金币+20, ★★★很有帮助, 非常感谢 2015-05-17 19:02:39
IBRION=6的时候 有关parallel的开关只能选择默认值（NPAR NCORE之类的）

我觉得弹性常数最好还是手动算，尤其是对于对称性不高的情况。通过施加strain，每一个类型的strain多弄几个点，在ISIF=2下固定cell，relax所有原子，从OUTCAR中获得stress，然后线性拟合一下就可以获得弹性常数。这样一般都可以获得比较准确的值。利用stress-strain法的话，你的相关参数都要调高一点（ENCUT KPOINTS EDIFFG之类的），否则OUTCAR里面的stress不准确。

或者，也可以利用energy-strain的方法计算，参数要求没有stress-strain那么高，但是有多少个独立的弹性常数就要对cell进行多少组变形，计算量很大。

ref：
stress-strain法：
http://journals.aps.org/prb/abstract/10.1103/PhysRevB.65.104104
http://scitation.aip.org/content ... 0/10.1063/1.2711762
energy-strain法：
http://scitation.aip.org/content ... /9/10.1063/1.368733 (D9的第二行分母我觉得应该是减号，否则不volume conserving)
赞一下 回复此楼
2楼2015-05-17 18:49:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Vaucanson

金虫 (小有名气)
引用回帖:
2楼: Originally posted by awdawdawdd at 2015-05-17 18:49:02
IBRION=6的时候 有关parallel的开关只能选择默认值（NPAR NCORE之类的）

我觉得弹性常数最好还是手动算，尤其是对于对称性不高的情况。通过施加strain，每一个类型的strain多弄几个点，在ISIF=2下固定cell，rela ...

非常感谢你的回答，您说的手动算确实是一个选择，我现在尝试的是面心立方的铝的弹性常数的计算，像这样高对称的结构如果用VASP自带的功能计算而出现我说的提示不知道是否要紧？因为它会继续计算。
另外，不太明白您说的“parallel的开关只能选择默认值”，好像NCORE不存在默认值吧？是NCORE =1的意思？
最后，非常感谢你附的参考文献。
赞一下 回复此楼
3楼2015-05-17 19:08:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

awdawdawdd

金虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
Vaucanson: 金币+20, ★★★很有帮助, 多谢回答 2015-05-18 21:07:14
franch: 金币+4, 应助指数+1, 谢谢回帖交流。。。。 2015-05-20 11:16:54
引用回帖:
3楼: Originally posted by Vaucanson at 2015-05-17 19:08:45
非常感谢你的回答，您说的手动算确实是一个选择，我现在尝试的是面心立方的铝的弹性常数的计算，像这样高对称的结构如果用VASP自带的功能计算而出现我说的提示不知道是否要紧？因为它会继续计算。
另外，不太明白 ...

有关parallelisation，主要是关系到multi-core multi-node的运算的效率
NCORE和NPAR是差不多的，INCAR内声明一个即可。NCORE默认值为1。parallelisation只会影响运算效率，一般不影响运算的结果。

但IBRION=5，6貌似都只能用默认值算(也就是你不在INCAR声明NPAR或者NCORE，或者把这些注释掉)，然后你用了默认值，输出文件里面肯定会提示你NPAR或者NCORE使用默认值计算的效率不高，此时忽略警告即可（因为你只能用默认值）。

还有，IBRION=6默认应该只取两个点（NFREE=2），加上strain=0的情况总共3个点，个人觉得不是很准确，因为拟合直线的话一般多取一些比较好（例如NFREE=4），然后调整好合适的strain的大小（POTIM默认为0.015Angst，根据你的体系调整）

ref：
parallelisation  http://cms.mpi.univie.ac.at/vasp ... PLANE_KPAR_tag.html
NCORE http://cms.mpi.univie.ac.at/wiki/index.php/NCORE
IBRION=6、NFREE、POTIM： http://cms.mpi.univie.ac.at/vasp/vasp/IBRION_5_IBRION_6.html
赞一下 回复此楼
4楼2015-05-17 19:24:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 22408调剂求助 +7 毂12 2026-04-09 9/450 2026-04-11 09:23 by 哦哦123
[考研] 085506-求调剂-285分 +3 雷欧飞踢 2026-04-08 3/150 2026-04-11 08:37 by zhq0425
[考研] 一志愿211，化学学硕，310分，本科重点双非，求调剂 +17 努力奋斗112 2026-04-06 20/1000 2026-04-11 00:31 by wangjihu
[基金申请] 山东省基金2026 +3 jerry681 2026-04-08 4/200 2026-04-10 16:02 by jerry681
[考研] 机械专硕270求调剂，接受跨专业 +10 老师看看我吧aba 2026-04-09 11/550 2026-04-10 10:41 by cal0306
[考研] 调剂申请086000一志愿西北农林科技大学生物与医药320分-本科齐鲁工业大学 +3 美美女士 2026-04-09 3/150 2026-04-10 10:31 by liuhuiying09
[考研] 求调剂 +3 电气300求调剂不 2026-04-08 3/150 2026-04-10 10:19 by LHGeng
[考研] 求 +6 化工专硕323分 2026-04-04 6/300 2026-04-10 10:04 by may_新宇
[考研] 一志愿中科大070300化学，314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[考研] 调剂 +19 2261744733 2026-04-08 19/950 2026-04-09 19:11 by vgtyfty
[考研] 312求调剂 +3 李鸿飞飞 2026-04-06 3/150 2026-04-09 17:32 by wp06
[考研] 一志愿郑州大学 22408 305分求调剂 +3 安小满zzz 2026-04-08 3/150 2026-04-09 16:16 by wp06
[考研] 求调剂 +8 吃口冰激凌 2026-04-07 8/400 2026-04-09 08:03 by 5268321
[考研] 266调剂 +8 daya sun 2026-04-07 9/450 2026-04-08 20:27 by yutian743
[考研] 338求调剂 +8 wxygxsaaaaa 2026-04-06 8/400 2026-04-08 06:58 by 无际的草原
[考研] 388求调剂 +6 四川王涛 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[考研] 一志愿华中农业大学0710（A）初试329分 求调剂 +5 一名26考研生 2026-04-04 5/250 2026-04-07 08:54 by 18828373951
[考研] 材料调剂 +5 小刘同学吖吖 2026-04-06 5/250 2026-04-06 18:34 by sherry_1901
[考研] 324求调剂 +3 k可乐 2026-04-05 4/200 2026-04-06 09:54 by 蓝云思雨
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
信息提示
关闭
请填处理意见
关闭