| 2 | 1/1 | 返回列表 |
| 查看: 1531 | 回復(fù): 1 | |||
zz1599新蟲 (初入文壇)
|
[交流]
vasp計算合金的力學(xué)性質(zhì)總是出錯
|
|
計算Th2Zn17結(jié)構(gòu)的Sm2Co17合金的力學(xué)性能。首先優(yōu)化該結(jié)構(gòu)的晶格常數(shù)。該結(jié)構(gòu)為菱方結(jié)構(gòu),一個晶胞中有57個原子,其中6個Sm原子,51個Co原子。 1.首先是優(yōu)化晶格常數(shù) 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文件 根據(jù)cif文件導(dǎo)出 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.加應(yīng)力 得到優(yōu)化之后的晶格常數(shù),然后建立新的文件夾,把優(yōu)化得到的CONTCAR拷為POSCAR,在POSCAR的三行基矢中,將第一行x軸的基矢改為原來的99.5%、99%、98%(即8.4020004272的99.5%、99%、98%等),相當(dāng)于給材料施加了0.5%、1%、2%的應(yīng)力。INCAR文件不變,KPOINTS不變,得到不同應(yīng)變對應(yīng)的體系能量。然后做一次自洽運算,將CONTCAR拷為POSCAR,NSW改為0,IBRION改為-1。這時,不同的應(yīng)變對應(yīng)的能量完全沒有規(guī)律,且加了應(yīng)變的能量和自洽之后的能量相差20eV。以下是部分結(jié)果 優(yōu)化后 -424.27 0.1press/加應(yīng)力能量-429.46/自洽-405.59 0.5press/加應(yīng)力能量-429.44/-406.83 1.5press/加應(yīng)力能量-429.31/-405.86 2press/加應(yīng)力能量-429.16/406 3press/加應(yīng)力能量-428.82/-406.08 4press/加應(yīng)力能量-405.89/-406.17(未收斂) 我想問的是 1.整體的計算思路是否有問題? 2.具體在加應(yīng)力或者自洽的過程中,是否是將CONTCAR拷為POSCAR,NSW改為0,IBRION改為-1。有沒有其他的文件需要一起拷貝過來?要不要修改INCAR文件中的其他命令? 3.為什么加了應(yīng)變之后能量沒有規(guī)律性?而且和自洽得到的能量相差20eV? |
新蟲 (初入文壇)
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 標(biāo)題:撿漏預(yù)警|08工科/09農(nóng)學(xué)調(diào)劑!英語要求低,過線即有機會! +9 | 馬超放煙花 2026-03-07 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 304求調(diào)劑(085602一志愿985) +9 | 化工人999 2026-03-09 | 9/450 |
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +4 | mysdl 2026-03-11 | 6/300 |
|
|
[考研] 0703化學(xué)一志愿211 總分320求調(diào)劑 +3 | 瑪卡巴卡啊哈 2026-03-11 | 3/150 |
|
|
[考研] 288求調(diào)劑 +13 | 王曉陽- 2026-03-09 | 18/900 |
|
|
[碩博家園] 深圳大學(xué)碩士招生(2026秋,傳感器方向,僅錄取第一志愿) +4 | xujiaoszu 2026-03-11 | 6/300 |
|
|
[考研] 大連大學(xué)化學(xué)專業(yè)研究生調(diào)劑 +3 | 琪久. 2026-03-10 | 8/400 |
|
|
[考研] 哈工大材料324求調(diào)劑 +6 | 閆旭東 2026-03-10 | 8/400 |
|
|
[考研] 085602化工求調(diào)劑 +7 | 董boxing 2026-03-10 | 7/350 |
|
|
[考研] 一志愿天津大學(xué),英一數(shù)二305分求調(diào)劑,四六級已過 +7 | 小小番的茄 2026-03-09 | 7/350 |
|
|
[考研] 一志愿天大化工(085600)調(diào)劑總分338 +5 | 蔡大美女 2026-03-09 | 5/250 |
|
|
[考研] 085600材料與化工 326 求調(diào)劑 +4 | 熱愛生活ing 2026-03-09 | 4/200 |
|
|
[考研] 數(shù)二英二309分請求調(diào)劑 +3 | dtdxzxx 2026-03-09 | 4/200 |
|
|
[考研] 288求調(diào)劑085600材料與化工 +13 | Daunrin 2026-03-07 | 15/750 |
|
|
[考研] 288求調(diào)劑(0703)一志愿東北大學(xué) +5 | 好好- 2026-03-07 | 5/250 |
|
|
[考研] 346分材料求調(diào)劑 +5 | snow_反季節(jié)版 2026-03-07 | 5/250 |
|
|
[考研] 材料專碩323求調(diào)劑 +5 | 李白26 2026-03-07 | 5/250 |
|
|
[考研] 08工科求調(diào)劑 +3 | 隆LLL 2026-03-06 | 4/200 |
|
|
[考研] 不限學(xué)校專業(yè)的調(diào)劑同學(xué)看過來 +5 | 啊擺啊擺 2026-03-05 | 9/450 |
|
|
[考研] 287求調(diào)劑 +3 | 看看我. 2026-03-05 | 6/300 |
|