 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 1527 | 回復(fù): 1 | |||
zz1599新蟲 (初入文壇)
|
[交流]
vasp計(jì)算合金的力學(xué)性質(zhì)總是出錯(cuò)
|
|
計(jì)算Th2Zn17結(jié)構(gòu)的Sm2Co17合金的力學(xué)性能。首先優(yōu)化該結(jié)構(gòu)的晶格常數(shù)。該結(jié)構(gòu)為菱方結(jié)構(gòu)，一個(gè)晶胞中有57個(gè)原子，其中6個(gè)Sm原子，51個(gè)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形式的贗勢(shì)文件 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)變對(duì)應(yīng)的體系能量。然后做一次自洽運(yùn)算，將CONTCAR拷為POSCAR，NSW改為0，IBRION改為-1。這時(shí)，不同的應(yīng)變對(duì)應(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.整體的計(jì)算思路是否有問題？ 2.具體在加應(yīng)力或者自洽的過程中，是否是將CONTCAR拷為POSCAR，NSW改為0，IBRION改為-1。有沒有其他的文件需要一起拷貝過來？要不要修改INCAR文件中的其他命令？ 3.為什么加了應(yīng)變之后能量沒有規(guī)律性？而且和自洽得到的能量相差20eV？ |
回復(fù)此樓» 猜你喜歡
誠求調(diào)劑，323分 有專利/科研/實(shí)習(xí)經(jīng)歷
已經(jīng)有5人回復(fù)
-
A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向
已經(jīng)有43人回復(fù)
-
26申博自薦
已經(jīng)有5人回復(fù)
-
售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】
已經(jīng)有7人回復(fù)
-
一志愿南京大學(xué)資源與環(huán)境專業(yè)283分求調(diào)劑
已經(jīng)有4人回復(fù)
-
085701環(huán)境工程295求調(diào)劑，一志愿A區(qū)211
已經(jīng)有7人回復(fù)
-
我嘗試在小木蟲上招計(jì)算機(jī)專業(yè)的調(diào)劑學(xué)生
已經(jīng)有5人回復(fù)
-
接受26屆調(diào)劑生
已經(jīng)有6人回復(fù)
-
一志愿哈爾濱工業(yè)大學(xué)0856材料與化工，前三科206，總分283，求調(diào)劑
已經(jīng)有6人回復(fù)
-
考研一志愿長安大學(xué)材料與化工309分請(qǐng)求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
zz1599
新蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 26.3
- 帖子: 31
- 在線: 50.5小時(shí)
- 蟲號(hào): 3599670
- 注冊(cè): 2014-12-17
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
2樓2019-01-23 08:44:06
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 zz1599 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 招調(diào)劑 +11 | 帆船哥 2026-03-04 | 13/650 |
|
|---|---|---|---|---|
|
|
[考博] 售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】 +3 | ipys00zm58 2026-03-06 | 4/200 |
|
|
|
[考研] 歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組 +19 | lisimayy 2026-03-04 | 31/1550 |
|
|
|
[考研] 0856材料與化工求調(diào)劑！ +5 | 化工考生111 2026-03-04 | 11/550 |
|
|
|
[考研] 一志愿武理085601專碩347分求調(diào)劑 +4 | 啊歐歐歐 2026-03-04 | 5/250 |
|
|
|
[考研] 0856材料求調(diào)劑 +17 | hyf hyf hyf 2026-02-28 | 18/900 |
|
|
[考研] 一志愿中石油，化學(xué)工程，301求調(diào)劑 +3 | 豚饅枳枳 2026-03-04 | 3/150 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑，總分320分，一志愿中南民族大學(xué) +3 | 19286712440 2026-03-04 | 3/150 |
|
|
|
[考研] 求調(diào)劑 +4 | 熬夜的貓頭鷹 2026-03-02 | 4/200 |
|
|
|
[考研] 292求調(diào)劑 +9 | yhk_819 2026-02-28 | 9/450 |
|
|
|
[考研] 291求調(diào)劑 +3 | 咕嚕咕嚕123123 2026-03-02 | 4/200 |
|
|
|
[考博] 26申博 求博導(dǎo) +3 | 愛讀書的小帥 2026-02-28 | 5/250 |
|
|
|
[考研] 267求調(diào)劑 +6 | 釣魚佬as 2026-03-02 | 6/300 |
|
|
|
[考研] 338求調(diào)劑 +5 | 18162027187 2026-03-02 | 6/300 |
|
|
|
[考研] 272求調(diào)劑 +9 | 材紫有化 2026-02-28 | 9/450 |
|
|
|
[考研] 265分求調(diào)劑不調(diào)專業(yè)和學(xué)校有行學(xué)上就 +6 | 禮堂丁真258 2026-02-28 | 9/450 |
|
|
|
[基金申請(qǐng)] 此成果不能導(dǎo)入原因：元數(shù)據(jù)必填信息不完整，可 進(jìn)行補(bǔ)充。 +4 | Kittylucky 2026-03-02 | 5/250 |
|
|
|
[考研] 0856材料求調(diào)劑 +4 | 麻辣魷魚 2026-02-28 | 4/200 |
|
|
|
[基金申請(qǐng)]
剛錄用，沒有期刊號(hào)，但是在線可看的論文可以放為代表作嗎
 10+3
|
arang1 2026-03-01 | 3/150 |
|
|
|
[考研] 317一志愿華南理工電氣工程求調(diào)劑 +6 | Soliloquy_Q 2026-02-28 | 11/550 |
|