| 查看: 5621 | 回复: 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 不知道是否对结果有什么影响？ 谢谢大家了。 |
回复此楼» 收录本帖的淘帖专辑推荐
 弹性计算 |
» 猜你喜欢
博士申请
已经有6人回复
-
评审有感
已经有11人回复
-
面上本子正文33页，违规吗？会被低分嘛？
已经有10人回复
-
国自然上会要求
已经有9人回复
-
今年审到国自然15份，谈谈感受
已经有16人回复
-
上海大学实验技术岗位非升即走
已经有8人回复
-
考博自荐
已经有6人回复
-
青C资助名额大幅增加！
已经有16人回复
-
重磅！青年科学基金项目（C类）资助增幅预计超过50%
已经有10人回复
-
我在等一个没有答案的答案
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于BCC-Fe的弹性常数的能量-应变计算
已经有47人回复
- 采用VASP计算得到的禁带宽度为什么比实验值低
已经有11人回复
- 救命啊！VASP计算磁性出错
已经有14人回复
- 感觉VASP计算缺陷存在很大问题
已经有17人回复
- VASP算DFT+U求助
已经有7人回复
- VASP做结构优化时用ISIF=3只想对晶格常数中的一个维度放开经行优化
已经有7人回复
- 弹性常数的计算
已经有6人回复
- VASP计算TiO2表面能问题 求助
已经有8人回复
- vasp计算弹性时的错误（使用defvector方法）
已经有25人回复
- vasp自旋极化计算时，初始设置自旋能否不被优化
已经有6人回复
- vasp5.2如何计算DOS和band structure？
已经有10人回复
- 计算弹性常数时的问题
已经有10人回复
- 弹性常数的计算
已经有12人回复
- VASP5.2 直接计算弹性常数
已经有52人回复
- vasp计算中，如何选择正确的K点
已经有7人回复
- 【求助】vasp 计算 INCAR的设置
已经有9人回复
- 【求助】一定压强下的弹性常数计算方法
已经有25人回复
- 【求助】算VASP弹性常数时的能量-应变出错
已经有18人回复
- 【讨论】用vasp计算过弹性常数的请进！！更新中！！
已经有52人回复
11楼2015-06-14 20:20:17
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
Li小喵
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 91
- 帖子: 39
- 在线: 29.3小时
- 虫号: 3480917
- 注册: 2014-10-17
- 性别: MM
- 专业: 凝聚态物性 II ：电子结构
5楼2015-05-19 15:15:49
huangsunchao
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1437
- 散金: 144
- 帖子: 70
- 在线: 75.1小时
- 虫号: 1306932
- 注册: 2011-05-26
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
6楼2015-06-12 09:27:32
7楼2015-06-13 09:46:51
vasp001
金虫 (正式写手)
- 应助: 69 (初中生)
- 金币: 1513.7
- 散金: 134
- 红花: 24
- 帖子: 702
- 在线: 326小时
- 虫号: 727265
- 注册: 2009-03-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
你好，我最近使用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数值？ 还请不吝指教。 另外你有即时沟通工具吗？ 有的话，方便提供一下吗？ |
8楼2015-06-14 16:04:43
vasp001
金虫 (正式写手)
- 应助: 69 (初中生)
- 金币: 1513.7
- 散金: 134
- 红花: 24
- 帖子: 702
- 在线: 326小时
- 虫号: 727265
- 注册: 2009-03-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
你好，我最近使用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数值？ 还请不吝指教。 另外你有即时沟通工具吗？ 有的话，方便提供一下吗？ |
9楼2015-06-14 16:04:56
vasp001
金虫 (正式写手)
- 应助: 69 (初中生)
- 金币: 1513.7
- 散金: 134
- 红花: 24
- 帖子: 702
- 在线: 326小时
- 虫号: 727265
- 注册: 2009-03-20
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
|
你好，我最近使用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数值？ 还请不吝指教。 另外你有即时沟通工具吗？ 有的话，方便提供一下吗？ |
10楼2015-06-14 16:05:03
30