| 查看: 1953 | 回復(fù): 10 | ||||
[求助]
Elk與QE聯(lián)合使用,是否可行?
|
||||
|
最近利用ELK計(jì)算一系列TM-Al體系,涉及空位、元素?fù)诫s等。遇到了計(jì)算量的問(wèn)題。 (1)首先,含缺陷體系屬于超晶胞,用ELK來(lái)進(jìn)行優(yōu)化,幾乎是不可能完成的任務(wù)。于是利用QE,進(jìn)行VC-relax,得到結(jié)構(gòu)用于ELK計(jì)算。這樣是否會(huì)被雜志接受?為了減少計(jì)算量,QE VC-relax時(shí),可否不考慮自旋激化? (2)其次,在ELK計(jì)算中,為了得到結(jié)合能,需要計(jì)算sublattice的能量,但煩人的是,當(dāng)對(duì)總體系計(jì)算時(shí),空帶已經(jīng)足夠了,但對(duì)TM的sublattice總提示空帶不夠,如果增加空帶,這不就不在一個(gè)精度或水平下來(lái)進(jìn)行比較,得到的結(jié)合能就不可靠了,不知理解是否正確? (3)想要得到彈性模量,如果用ELK計(jì)算不同體積下的基態(tài)能量,然后擬合得到,這個(gè)計(jì)算量也相當(dāng)要命。如果只用QE計(jì)算的話,感覺(jué)分量不夠! 各位有何見(jiàn)解?請(qǐng)指教! |
金蟲 (著名寫手)
|
如果是磁性體系,不加自旋會(huì)可能獲得不同的結(jié)果。因?yàn)橛辛俗孕,體系的能量和離子受力都會(huì)不一樣。當(dāng)然,第一步你可以不加自旋優(yōu)化,在此之后,再加自旋。 此外,結(jié)合能肯定都要在同一個(gè)標(biāo)準(zhǔn)尺度下進(jìn)行計(jì)算,這樣才真的有了比較的標(biāo)準(zhǔn),要不沒(méi)法比較。如果你的體系不涉及內(nèi)殼層電子,也沒(méi)有必要用ELK。我不明白你說(shuō)的分量不夠是什么意思,這些都主要看其中的物理實(shí)現(xiàn)是什么。 |
專家顧問(wèn) (著名寫手)
![]() |
專家經(jīng)驗(yàn): +174 |
|
即使是同一個(gè)軟件,使用不同的XC泛函或者贗勢(shì),都可能收斂的模型尺寸有出入,更何況一個(gè)是平面波贗勢(shì)方法,一個(gè)是全勢(shì)線性綴加波方法的呢……最好是用QE做個(gè)大體的結(jié)構(gòu)出來(lái),然后再用ELK仔細(xì)查找,并且可以把QE的工作在此省略。 至于只用QE計(jì)算,是完全可以的,別說(shuō)QE這種功能強(qiáng)悍的,就是VASP那樣功能比較少的軟件,在不掛其他輔助程序的前提下做一個(gè)有好的idea的計(jì)算,照樣也能發(fā)nature級(jí)的文章啊~ |

木蟲 (正式寫手)
star watcher
專家顧問(wèn) (著名寫手)
![]() |
專家經(jīng)驗(yàn): +174 |

|
謝謝!的確如你所言,結(jié)果剛出來(lái),我就是先沒(méi)加自旋優(yōu)化,然后在此基礎(chǔ)上加上自旋,得到的結(jié)果還是有比較明顯的差異。因?yàn)镕eAl金屬間化合物中,B2結(jié)構(gòu)的FeAl實(shí)驗(yàn)基態(tài)是非磁態(tài),但無(wú)論是GGA還是LDA都只能得到鐵磁基態(tài),當(dāng)采用LDA+U后,可以得到非磁基態(tài)。文獻(xiàn)認(rèn)為是Fe的3d電子的t2g和eg劈裂時(shí)在U的引入后所占比重不同。 我也面臨這樣一個(gè)問(wèn)題, 當(dāng)對(duì)FeAl單胞選擇U=6.33 eV可得到非磁基態(tài),但對(duì)3x3x3FeAl時(shí),相同的計(jì)算參數(shù),磁性有所增加,約為0.04。難道還要針對(duì)這樣的超晶胞重新選擇U值?在ELK里就會(huì)計(jì)算量很大。 另外,當(dāng)計(jì)算3x3x3FeAl中心一個(gè)Fe空位時(shí),ELK的輸入文件是否需要把vkloff做一定的偏移? |
專家顧問(wèn) (著名寫手)
![]() |
專家經(jīng)驗(yàn): +174 |

金蟲 (著名寫手)
專家顧問(wèn) (著名寫手)
![]() |
專家經(jīng)驗(yàn): +174 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿天大材料與化工(085600)總分338 +4 | 蔡大美女 2026-03-13 | 4/200 |
|
|
[考研] 一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑 +4 | Ymlll 2026-03-18 | 6/300 |
|
|
[考研] 330求調(diào)劑 +3 | 小材化本科 2026-03-18 | 3/150 |
|
|
[考研] 274求調(diào)劑 +5 | S.H1 2026-03-18 | 5/250 |
|
|
[考研] 085600材料與化工 +5 | 安全上岸! 2026-03-16 | 5/250 |
|
|
[考研] 266求調(diào)劑 +5 | 陽(yáng)陽(yáng)哇塞 2026-03-14 | 9/450 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +6 | zbzihdhd 2026-03-15 | 7/350 |
|
|
[考研] 268求調(diào)劑 +6 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 6/300 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考研] 材料與化工求調(diào)劑 +6 | 為學(xué)666 2026-03-16 | 6/300 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[基金申請(qǐng)]
今年的國(guó)基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 26考研一志愿中國(guó)石油大學(xué)(華東)305分求調(diào)劑 +3 | 嘉年新程 2026-03-15 | 3/150 |
|
|
[考研] 材料專碩350 求調(diào)劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|