| 查看: 1535 | 回复: 1 | |||
zz1599新虫 (初入文坛)
|
[交流]
vasp计算合金的力学性质总是出错
|
|
计算Th2Zn17结构的Sm2Co17合金的力学性能。首先优化该结构的晶格常数。该结构为菱方结构，一个晶胞中有57个原子，其中6个Sm原子，51个Co原子。 1.首先是优化晶格常数 INCAR文件 SYSTEM = Sm2Co17R ENCUT = 350 ISTART = 0; ICHARG = 2 ISMEAR = 1; SIGMA = 0.2 NSW = 0; IBRION = 2 NELM = 200 ISIF = 2 POTIM = 0.5 EDIFF = 1E-5 ISPIN = 2 MAGMOM = 6*1 51*3 PREC = Medium NCORE = 4 LREAL = Auto ALGO = Fast POSCAR文件 根据cif文件导出 Sm2Co17R 1.0 8.4020004272 0.0000000000 0.0000000000 -4.2010002136 7.2763458126 0.0000000000 0.0000000000 0.0000000000 12.1719999313 Sm Co 6 51 Direct -0.000000000 -0.000000000 0.333333346 -0.000000000 0.000000000 0.666666614 0.666666670 0.333333340 0.666666693 0.666666670 0.333333340 0.000000000 0.333333340 0.666666681 0.000000000 0.333333340 0.666666681 0.333333346 0.000000000 0.000000000 0.096999998 -0.000000000 -0.000000000 0.902999992 0.666666670 0.333333340 0.430333354 0.666666670 0.333333340 0.236333338 0.333333340 0.666666681 0.763666701 0.333333340 0.666666681 0.569666685 0.500000000 0.000000000 0.500000000 -0.000000003 0.499999994 0.500000000 0.499999997 0.499999994 0.500000000 0.166666755 0.333333340 0.833333346 0.666666639 0.833333335 0.833333346 0.166666667 0.833333335 0.833333346 0.833333397 0.666666681 0.166666752 0.333333319 0.166666752 0.166666752 0.833333400 0.166666687 0.166666752 0.333333333 0.000000000 0.000000000 0.666666610 0.000000000 0.000000000 0.000000004 0.333333340 0.000000000 0.000000003 0.666666615 0.000000000 0.666666613 0.666666615 0.000000000 0.333333337 0.333333340 0.000000000 0.000000004 0.333333340 0.333333346 0.333333337 0.333333340 0.333333346 0.666666674 0.666666681 0.333333346 0.666666667 -0.000000000 0.333333346 0.333333333 -0.000000000 0.333333346 0.000000007 0.666666681 0.333333346 0.666666674 0.666666681 0.666666693 0.000000007 0.666666681 0.666666693 0.333333333 -0.000000000 0.666666693 0.333333337 0.333333340 0.666666693 0.000000004 0.333333340 0.666666693 0.666666667 -0.000000000 0.666666693 0.499999997 0.499999994 0.166666673 0.499999997 0.499999994 0.833333346 0.500000000 -0.000000000 0.166666673 0.500000000 0.000000000 0.833333346 -0.000000003 0.499999994 0.166666673 -0.000000003 0.499999994 0.833333346 0.166666724 0.833333335 0.500000000 0.166666667 0.833333335 0.166666673 0.166666684 0.333333340 0.500000000 0.166666755 0.333333340 0.166666673 0.666666639 0.833333335 0.500000000 0.666666639 0.833333335 0.166666673 0.833333376 0.166666752 0.833333346 0.833333400 0.166666687 0.500000000 0.833333397 0.666666681 0.833333346 0.833333397 0.666666681 0.500000000 0.333333343 0.166666687 0.833333346 0.333333319 0.166666752 0.500000000 KPOINTS文件 A 0 G 4 4 3 0 0 0 POTCAR文件 采用GGA_PBE形式的赝势文件 2.加应力 得到优化之后的晶格常数，然后建立新的文件夹，把优化得到的CONTCAR拷为POSCAR，在POSCAR的三行基矢中，将第一行x轴的基矢改为原来的99.5%、99%、98%（即8.4020004272的99.5%、99%、98%等），相当于给材料施加了0.5%、1%、2%的应力。INCAR文件不变，KPOINTS不变，得到不同应变对应的体系能量。然后做一次自洽运算，将CONTCAR拷为POSCAR，NSW改为0，IBRION改为-1。这时，不同的应变对应的能量完全没有规律，且加了应变的能量和自洽之后的能量相差20eV。以下是部分结果 优化后 -424.27 0.1press/加应力能量-429.46/自洽-405.59 0.5press/加应力能量-429.44/-406.83 1.5press/加应力能量-429.31/-405.86 2press/加应力能量-429.16/406 3press/加应力能量-428.82/-406.08 4press/加应力能量-405.89/-406.17(未收敛) 我想问的是 1.整体的计算思路是否有问题？ 2.具体在加应力或者自洽的过程中，是否是将CONTCAR拷为POSCAR，NSW改为0，IBRION改为-1。有没有其他的文件需要一起拷贝过来？要不要修改INCAR文件中的其他命令？ 3.为什么加了应变之后能量没有规律性？而且和自洽得到的能量相差20eV？ |
回复此楼» 猜你喜欢
281求调剂（0805）
已经有16人回复
-
304求调剂
已经有6人回复
-
材料工程专硕调剂
已经有6人回复
-
一志愿天大材料与化工（085600）总分338
已经有4人回复
-
085700资源与环境308求调剂
已经有3人回复
-
求材料调剂
已经有8人回复
-
294求调剂材料与化工专硕
已经有5人回复
-
一志愿华中科技大学，080502，354分求调剂
已经有4人回复
-
一志愿吉林大学材料学硕321求调剂
已经有6人回复
-
085410人工智能专硕317求调剂（0854都可以）
已经有3人回复
» 本主题相关商家推荐: (我也要在这里推广)
zz1599
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 26.3
- 帖子: 31
- 在线: 50.5小时
- 虫号: 3599670
- 注册: 2014-12-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2019-01-23 08:44:06
