| 5 | 1/1 | 返回列表 |
| 查看: 3597 | 回復(fù): 14 | ||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||
[求助]
救命啊!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)閷Ρ任墨I(xiàn),磁矩該大的地方小,該小的地方大。希望大神出來解救小弟,我已經(jīng)搗鼓一天了,試了N多種情況,求解答,不勝感激。 |
木蟲 (正式寫手)
木蟲 (正式寫手)
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 300求調(diào)劑,材料科學(xué)英一數(shù)二 +8 | leaflight 2026-03-24 | 8/400 |
|
|---|---|---|---|---|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研]
|
y7czhao 2026-03-26 | 10/500 |
|
|
[考研] 317求調(diào)劑 +6 | 十閑wx 2026-03-24 | 6/300 |
|
|
[考研] 311求調(diào)劑 +4 | 冬十三 2026-03-24 | 4/200 |
|
|
[考研] 265求調(diào)劑11408 +3 | 劉小鹿lu 2026-03-27 | 3/150 |
|
|
[考研] 348求調(diào)劑 +4 | 小懶蟲不懶了 2026-03-27 | 5/250 |
|
|
[考研] 考研調(diào)劑 +10 | 呼呼?~+123456 2026-03-24 | 10/500 |
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 341求調(diào)劑 +7 | 青檸檬1 2026-03-26 | 7/350 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研]
材料調(diào)劑
5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
[考研] 352求調(diào)劑 +4 | 大米飯! 2026-03-22 | 4/200 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +4 | 李李不服輸 2026-03-26 | 4/200 |
|
|
[考研] 打過很多競賽,085406控制工程300分,求調(diào)劑 +3 | askeladz 2026-03-26 | 3/150 |
|
|
[考研] 材料專碩 335 分求調(diào)劑 +4 | 拒絕冷暴力 2026-03-25 | 4/200 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|