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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 关于用VASP计算弹性常数的两个问题
51/1返回列表
查看: 5450  |  回复: 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

金虫 (正式写手)
引用回帖:
4楼: Originally posted by awdawdawdd at 2015-05-17 19:24:54
有关parallelisation，主要是关系到multi-core multi-node的运算的效率
NCORE和NPAR是差不多的，INCAR内声明一个即可。NCORE默认值为1。parallelisation只会影响运算效率，一般不影响运算的结果。

但IBRION=5， ...

你好，我最近使用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.
9楼2015-06-14 16:04:56
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 广东省 085601 329分求调剂 +13 Eddieddd 2026-04-10 13/650 2026-04-10 23:27 by 314126402
[考研] 291 求调剂 +29 化工2026届毕业� 2026-04-09 29/1450 2026-04-10 22:55 by dick_runner
[考研] 材料专业344求调剂 +16 hualkop 2026-04-10 21/1050 2026-04-10 17:28 by laoshidan
[考研] 301求调剂 +4 XYPLR 2026-04-05 6/300 2026-04-10 09:17 by Delta2012
[考研] 生物与医药273求调剂 +18 荔题南墙 2026-04-05 19/950 2026-04-10 08:14 by kangsm
[考研] 284求调剂 +7 archer.. 2026-04-09 8/400 2026-04-09 21:30 by bljnqdcc
[考研] 296求调剂 +5 汪！？！ 2026-04-09 5/250 2026-04-09 17:47 by 柠檬不酸zy
[考研] 085600材料与化工专硕329 求调剂 +24 额cc 2026-04-06 25/1250 2026-04-09 16:01 by wp06
[考研] 334求调剂 +16 Riot2025 2026-04-08 17/850 2026-04-09 09:28 by wdyheheeh
[考研] 求调剂 +7 chenxrlkx 2026-04-05 9/450 2026-04-09 09:04 by wj165256
[考研] 307分材料专业求调剂 +12 Hll胡 2026-04-05 12/600 2026-04-08 16:33 by luoyongfeng
[考研] 285求调剂 +7 恶法大二的气味� 2026-04-05 10/500 2026-04-08 14:34 by zhq0425
[考研] 259求调剂 +5 就爱吃土豆呀呀 2026-04-07 5/250 2026-04-07 22:40 by JourneyLucky
[考研] 312求调剂 +18 gtw1 2026-04-06 20/1000 2026-04-07 18:16 by 蓝云思雨
[考研] 一志愿苏州大学材料工程（085601）专硕有科研经历三项国奖两个实用型专利一项省级立项 +11 大火山小火山 2026-04-05 11/550 2026-04-06 22:55 by yunlongyang
[考研] 求调剂 +10 chenxrlkx 2026-04-05 10/500 2026-04-06 11:31 by 猪会飞
[考研] 复试调剂 +5 asdasdassda 2026-04-05 5/250 2026-04-06 09:32 by dongzh2009
[考研] 377求调剂 +6 by.ovo 2026-04-05 6/300 2026-04-05 22:18 by dongzh2009
[考研] 一志愿青科085500，初试295分，公共课213分 +3 遇到的人愿望都� 2026-04-05 3/150 2026-04-05 18:45 by 蓝云思雨
[考研] 325求调剂 +4 春风不借意 2026-04-04 4/200 2026-04-04 14:46 by 湘农储能材料
信息提示
关闭
请填处理意见
关闭