| 查看: 3596 | 回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +6 | c297914 2026-03-28 | 6/300 |
|
|
[考研] 317分 一志愿南理工材料工程 本科湖工大 求調(diào)劑 +7 | 芋泥小鈴鐺 2026-03-28 | 7/350 |
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
[考研] 322求調(diào)劑 +5 | 舊吢 2026-03-24 | 5/250 |
|
|
[考研] 086502化學(xué)工程342求調(diào)劑 +6 | 阿姨復(fù)古不過 2026-03-27 | 6/300 |
|
|
[考研] 求調(diào)劑推薦 材料 304 +15 | 荷包蛋hyj 2026-03-26 | 15/750 |
|
|
[考研] 化學(xué)調(diào)劑 +4 | 愛吃番茄的旭 2026-03-24 | 5/250 |
|
|
[考研] 283求調(diào)劑(080500) +4 | A child 2026-03-27 | 4/200 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研]
材料學(xué)碩,求調(diào)劑
6+5
|
糖葫蘆888ll 2026-03-22 | 10/500 |
|
|
[考研] 化學(xué)工程085602 305分求調(diào)劑 +17 | RichLi_ 2026-03-25 | 17/850 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
[考研] 293求調(diào)劑 +7 | 加一一九 2026-03-24 | 7/350 |
|
|
[考研] 一志愿武理085500機(jī)械專業(yè)總分300求調(diào)劑 +3 | an10101 2026-03-24 | 7/350 |
|
|
[考研] 300分,材料,求調(diào)劑,英一數(shù)二 +5 | 超贊的 2026-03-24 | 5/250 |
|
|
[考研] 求調(diào)劑一志愿武漢理工大學(xué)材料工程(085601) +5 | WW.' 2026-03-23 | 7/350 |
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達(dá)問津 2026-03-23 | 5/250 |
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|