| 2 | 1/1 | 返回列表 |
| 查看: 1526 | 回復(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)的體系能量。然后做一次自洽運(yùn)算,將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ā)表 | |
|---|---|---|---|---|
|
[考研] 復(fù)試調(diào)劑 +4 | 呼呼?~+123456 2026-03-05 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +3 | 馳光絕景 2026-03-04 | 3/150 |
|
|
[考研] 材料與化工,291,求調(diào)劑 +7 | 咕嚕咕嚕123123 2026-03-05 | 8/400 |
|
|
[考研] 301求調(diào)劑 +4 | 李LJR 2026-03-04 | 4/200 |
|
|
[考研] 085701環(huán)境工程 求調(diào)劑 +7 | xiiiia 2026-03-04 | 7/350 |
|
|
[考研] 22408-273求調(diào)劑-擔(dān)任3個項(xiàng)目的負(fù)責(zé)人-1篇國際期刊論文(一作)1篇核心期刊論文在投。 +3 | 沒想好取什么名 2026-03-03 | 3/150 |
|
|
[考研] 一志愿985材料與化工 326分求調(diào)劑 +3 | Hz795795 2026-03-04 | 3/150 |
|
|
[考研] 成績276,專業(yè)代碼0856求調(diào)劑 +10 | 小陳朵 2026-03-03 | 10/500 |
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|
[考研] 0856材料工程,初試313調(diào)劑 +7 | 賣個關(guān)子吧 2026-03-03 | 7/350 |
|
|
[基金申請] 請問大家,研究風(fēng)險與應(yīng)對措施那里, 大家都怎么寫呢 ? +3 | cauasen 2026-03-02 | 3/150 |
|
|
[考研] 292求調(diào)劑 +3 | sgbl 2026-03-03 | 3/150 |
|
|
[考研] 085602化學(xué)工程350,調(diào)劑,有沒有211的 +5 | 利好利好. 2026-03-02 | 9/450 |
|
|
[考研] 一志愿華中科技大學(xué),化學(xué)專業(yè)344分,求調(diào)劑 +6 | 邢xing1 2026-03-02 | 6/300 |
|
|
[考研]
|
glwshine 2026-03-02 | 5/250 |
|
|
[考研] 0856材料調(diào)劑 +5 | 沿岸有貝殼OUC 2026-03-02 | 5/250 |
|
|
[考研] 275求調(diào)劑 +7 | 明遠(yuǎn)求學(xué) 2026-03-01 | 7/350 |
|
|
[考研] 材料085601調(diào)劑 +5 | 多多子. 2026-03-02 | 5/250 |
|
|
[考研] 材料與化工328求調(diào)劑 +3 | 。,。,。,。i 2026-03-02 | 3/150 |
|