| 查看: 3593 | 回復(fù): 14 | ||
[求助]
救命啊!VASP計(jì)算磁性出錯(cuò) 已有2人參與
|
||
|
小弟在復(fù)現(xiàn)一篇關(guān)于PTO缺陷誘導(dǎo)磁性的文章(Vacancy-driven ferromagnetism in ferroelectric PbTiO3),我是先優(yōu)化PTO單胞,然后組建超胞,2*2*2,主要參數(shù)和這篇APL是一樣的,附我的INCAR Name of system SYSTEM = PbTiO3 Start parameter ISTART = 0 INIWAV = 1 ICHARG = 2 Write flags NWRITE = 2 LWAVE = T LCHARG = T Parallelization NPAR = 4 LPLANE = T LSCALU = F LSCALAPACK = F Electronic Relaxation 1 ENCUT = 500.00 eV EDIFF = 1.0E-06 eV PREC = High NELM = 200 NELMIN = 2 Electronic Relaxation 2 ALGO = Normal LDIAG = T LREAL = .FALSE. Ionic Relaxation #NSW = 150 #IBRION = 2 #ISIF = 3 #EDIFFG = -0.01 #POTIM = 0.5 #ISYM = 2 #IWAVPR = 1 #LCORR = T DOS related values ISMEAR = 0 SIGMA = 0.025 #NEDOS = 3000 # NBANDS = 3958 Magnetic parameter ISPIN = 2 LORBIT = 11 Exchange-correlation functional GGA = PE VOSKOWN = 1 MAGMOM = 8*1 7*1 24*1 我分別用ISIF=2(GGA),ISIF=3(GGA),ISIF=3(GGA+U)去優(yōu)化晶胞,然后自洽計(jì)算磁性,但結(jié)果和文獻(xiàn)的差距很大 比如,我用和文獻(xiàn)一樣的GGA,那么我用ISIF=3計(jì)算的話,當(dāng)初始磁矩為2或3,5的時(shí)候,總磁矩為0.88,文獻(xiàn)里Tii缺陷誘導(dǎo)的磁矩為2.223,我的結(jié)果小很多,但是,自旋電荷密度圖像是差不多的然后我把每個(gè)原子的初始磁矩設(shè)為9,結(jié)果很離譜了,磁矩為3.13,而且缺陷附近的原子磁矩大,遠(yuǎn)的磁矩小,明顯不符合事實(shí)。ISIF=2的話,和ISIF的情況類似,只是磁矩偏小,才0.3,加U也很離譜,我試著把初始磁矩設(shè)為0.5和9兩個(gè)極端,我發(fā)現(xiàn)結(jié)果一樣,磁矩為3.1,但是結(jié)果肯定是錯(cuò)誤的,因?yàn)閷?duì)比文獻(xiàn),磁矩該大的地方小,該小的地方大。希望大神出來解救小弟,我已經(jīng)搗鼓一天了,試了N多種情況,求解答,不勝感激。 |
木蟲 (正式寫手)
|
前輩,我經(jīng)過與一個(gè)師兄的討論,他跟我說,按照洪特定則和泡立不相容原理的電子排布,去給定初始磁矩,比如Pb電子排布4f145d106s26p2,那么初始磁矩為2,氧本身按電子排布沒有磁性,不設(shè)也沒關(guān)系,加入誘導(dǎo)出來的話,最終也會(huì)表現(xiàn)出來,我剛接觸磁性,即得以為前輩的話,是vasp磁矩是從大到小優(yōu)化的,給給的大了會(huì)優(yōu)化到合理的值,你給笑了,可能就出錯(cuò),我現(xiàn)在也不知道,初始磁矩該怎么給,是由特定的規(guī)則,還是我們要不斷的嘗試不同的大小,以得出不同的結(jié)果去對(duì)比。 |
木蟲 (正式寫手)
木蟲 (正式寫手)
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703本科鄭州大學(xué)求調(diào)劑 +3 | nhj_ 2026-03-25 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
[考研] 330一志愿中國海洋大學(xué) 化學(xué)工程 085602 有讀博意愿 求調(diào)劑 +3 | wywy.. 2026-03-27 | 4/200 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +4 | 傾____蕭 2026-03-27 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
[考研] 0703化學(xué)一志愿南京師范大學(xué)303求調(diào)劑 +3 | zzffylgg 2026-03-24 | 3/150 |
|
|
[考研] 351求調(diào)劑 +4 | 麥克阿磊 2026-03-24 | 4/200 |
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17 | RichLi_ 2026-03-25 | 17/850 |
|
|
[考研] 332求調(diào)劑 +6 | 032500 2026-03-25 | 6/300 |
|
|
[考研] 考研一志愿蘇州大學(xué)初始315(英一)求調(diào)劑 +3 | sbdksD 2026-03-24 | 4/200 |
|
|
[考研] 347求調(diào)劑 +4 | L when 2026-03-25 | 4/200 |
|
|
[考研] B區(qū)考研調(diào)劑 +4 | yqdszhdap- 2026-03-22 | 5/250 |
|
|
[考研] 一志愿北化315 求調(diào)劑 +3 | akrrain 2026-03-24 | 3/150 |
|
|
[考研] 一志愿南航材料專317分求調(diào)劑 +5 | 炸呀炸呀炸薯?xiàng)l 2026-03-23 | 5/250 |
|
|
[考研] 300求調(diào)劑,材料科學(xué)英一數(shù)二 +5 | leaflight 2026-03-24 | 5/250 |
|
|
[考研] 一志愿河北工業(yè)大學(xué)0817化工278分求調(diào)劑 +7 | jhybd 2026-03-23 | 12/600 |
|
|
[考研] 308求調(diào)劑 +3 | 墨墨漠 2026-03-21 | 3/150 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-22 | 6/300 |
|