| 查看: 1646 | 回复: 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？ |
回复此楼» 猜你喜欢
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有9人回复
-
时间戳今天，20号变了
已经有5人回复
-
2026国自然函评费到账
已经有9人回复
-
放榜前的不淡定
已经有14人回复
-
建议基金发布提前给出明确的时间点
已经有13人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
» 本主题相关商家推荐: (我也要在这里推广)
zz1599
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 26.3
- 帖子: 31
- 在线: 50.7小时
- 虫号: 3599670
- 注册: 2014-12-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2019-01-23 08:44:06
