| 查看: 1626 | 回复: 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？ |
回复此楼» 猜你喜欢
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有17人回复
-
2026年国自然面上资助率
已经有16人回复
-
小木虫看见有人已经知道结果了
已经有20人回复
-
微信指数没变化，科研之友没阅读
已经有18人回复
-
基础研究怎么拉横向，学校到款任务越来越多，难以完成 拉横向，都有哪些途径啊
已经有8人回复
-
HXDI做水性聚氨酯乳液，是不是特别容易出渣
已经有3人回复
-
系统今天又提示维护了，估计离放榜不远了
已经有15人回复
-
申请2027年材料/化学类博士
已经有5人回复
-
你们的时间戳变了吗
已经有4人回复
-
今年的WR进展到哪一步了？
已经有9人回复
» 本主题相关商家推荐: (我也要在这里推广)
zz1599
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 26.3
- 帖子: 31
- 在线: 50.7小时
- 虫号: 3599670
- 注册: 2014-12-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2019-01-23 08:44:06