| 查看: 1278 | 回复: 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？ |
回复此楼» 猜你喜欢
前些天开会有个人见到人就搞关系，一查此人全是MDPI/Hindawi论文，鄙视!
已经有11人回复
-
化学B02口青基 代表作都是什么水平的？向大佬求助
已经有16人回复
-
EI学报，一审返修后，为啥不再送审，直接终审中?
已经有7人回复
-
实验室太吵闹，无法安静学习，怎么办？
已经有16人回复
-
《控制理论与应用》期刊收版面费吗？
已经有5人回复
-
Scientific reports 投稿
已经有6人回复
-
电化学传感器
已经有5人回复
-
职能部门工作人员态度不好是普遍的吗？怎么让他们态度好一些？
已经有7人回复
-
为什么美国要封锁华为而不封锁丘C同呢？
已经有10人回复
-
每天学术时间不能保证，能保证的只有：
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
zz1599
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 26.3
- 帖子: 31
- 在线: 50.2小时
- 虫号: 3599670
- 注册: 2014-12-17
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2019-01-23 08:44:06
5