| 查看: 5455 | 回复: 12 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[求助]
关于用VASP计算弹性常数的两个问题 已有2人参与
|
|||
|
最近在学习用VASP 5.3.5计算弹性常数，设置IBRION = 6, ISIF = 4计算。发现有两个问题想请教大家： 1. 当我用16核计算时，设置NCORE = 4，走完第一个离子步后报错，停止计算：  这是否意味着要直接取消掉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 不知道是否对结果有什么影响？ 谢谢大家了。 |
回复此楼» 猜你喜欢
复试调剂
已经有11人回复
-
311求调剂
已经有13人回复
-
22408调剂求助
已经有9人回复
-
085410 273求调剂
已经有3人回复
-
材料考研调剂
已经有9人回复
-
085600材料与化工329分求调剂
已经有12人回复
-
0854调剂
已经有4人回复
-
调剂求收留
已经有14人回复
-
0831生医工第一轮调剂失败求助
已经有10人回复
-
一志愿东北大学控制工程085406数二英二385，求调剂
已经有8人回复
|
关于问题二的解决方案： The message does not mean an error. The procedure doing the evaluation of symmetry uses a numerical criterion (SYMPREC=1E-5) to determine which atoms are symmetrically equivalent. When this criterion is too soft the programm can find "false" symmetry and reports that "reciprocal lattice and k-lattice belong to different class of lattices." To get rid of this message you can: 1) Switch off the symmetry (ISYM=0) 2) Make the treshold for detecting the symmetry more severe (SYMPREC=1E-8) |
13楼2016-01-27 22:19:55
awdawdawdd
金虫 (小有名气)
- 应助: 24 (小学生)
- 金币: 679
- 散金: 200
- 红花: 6
- 帖子: 111
- 在线: 67小时
- 虫号: 2237233
- 注册: 2013-01-12
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
Vaucanson: 金币+20, ★★★很有帮助, 非常感谢 2015-05-17 19:02:39
感谢参与，应助指数 +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
3楼2015-05-17 19:08:45
awdawdawdd
金虫 (小有名气)
- 应助: 24 (小学生)
- 金币: 679
- 散金: 200
- 红花: 6
- 帖子: 111
- 在线: 67小时
- 虫号: 2237233
- 注册: 2013-01-12
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
Vaucanson: 金币+20, ★★★很有帮助, 多谢回答 2015-05-18 21:07:14
franch: 金币+4, 应助指数+1, 谢谢回帖交流。。。。 2015-05-20 11:16:54
Vaucanson: 金币+20, ★★★很有帮助, 多谢回答 2015-05-18 21:07:14
franch: 金币+4, 应助指数+1, 谢谢回帖交流。。。。 2015-05-20 11:16:54
|
有关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
5