| 查看: 1496 | 回复: 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？ |
回复此楼» 猜你喜欢
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
假如你的研究生提出不合理要求
已经有6人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
-
26申博（荧光探针方向，有机合成）
已经有4人回复
» 本主题相关商家推荐: (我也要在这里推广)
zz1599
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 26.3
- 帖子: 31
- 在线: 50.5小时
- 虫号: 3599670
- 注册: 2014-12-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2019-01-23 08:44:06