| 5 | 1/1 | 返回列表 |
| 查看: 3595 | 回復: 14 | ||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||
[求助]
救命啊!VASP計算磁性出錯 已有2人參與
|
||
|
小弟在復現(xiàn)一篇關(guān)于PTO缺陷誘導磁性的文章(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)化晶胞,然后自洽計算磁性,但結(jié)果和文獻的差距很大 比如,我用和文獻一樣的GGA,那么我用ISIF=3計算的話,當初始磁矩為2或3,5的時候,總磁矩為0.88,文獻里Tii缺陷誘導的磁矩為2.223,我的結(jié)果小很多,但是,自旋電荷密度圖像是差不多的然后我把每個原子的初始磁矩設(shè)為9,結(jié)果很離譜了,磁矩為3.13,而且缺陷附近的原子磁矩大,遠的磁矩小,明顯不符合事實。ISIF=2的話,和ISIF的情況類似,只是磁矩偏小,才0.3,加U也很離譜,我試著把初始磁矩設(shè)為0.5和9兩個極端,我發(fā)現(xiàn)結(jié)果一樣,磁矩為3.1,但是結(jié)果肯定是錯誤的,因為對比文獻,磁矩該大的地方小,該小的地方大。希望大神出來解救小弟,我已經(jīng)搗鼓一天了,試了N多種情況,求解答,不勝感激。 |
木蟲 (正式寫手)
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 學碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 315求調(diào)劑 +4 | akie... 2026-03-28 | 5/250 |
|
|
[考研] 346求調(diào)劑 一志愿070303有機化學 +3 | 蘿卜燉青菜 2026-03-28 | 3/150 |
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +6 | 腳滑的守法公民 2026-03-27 | 6/300 |
|
|
[考研] 265求調(diào)劑 +8 | 小木蟲085600 2026-03-27 | 8/400 |
|
|
[考研] 考研調(diào)劑 +4 | Sanmu-124 2026-03-26 | 4/200 |
|
|
[考研] 調(diào)劑 +3 | 李嘉圖·S·路 2026-03-27 | 3/150 |
|
|
[碩博家園] 北京林業(yè)大學碩導招生廣告 +6 | kongweilin 2026-03-26 | 8/400 |
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
[考研] 一志愿北化求調(diào)劑 +3 | Jsman 2026-03-22 | 3/150 |
|
|
[考研] 生物學 296 求調(diào)劑 +4 | 朵朵- 2026-03-26 | 6/300 |
|
|
[考研] 中國科學院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學院大學、深圳理工大學聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研] 271求調(diào)劑 +6 | 生如夏花… 2026-03-22 | 6/300 |
|
|
[考研] 352求調(diào)劑 +4 | 大米飯! 2026-03-22 | 4/200 |
|
|
[考研] 085601求調(diào)劑總分293英一數(shù)二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
[考研] 0854電子信息求調(diào)劑 324 +4 | Promise-jyl 2026-03-23 | 4/200 |
|
|
[考研] 一志愿吉林大學材料與化工303分求調(diào)劑 +4 | 為學666 2026-03-24 | 4/200 |
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當牛 2026-03-24 | 4/200 |
|
|
[論文投稿] 急發(fā)核心期刊論文 +3 | 賢達問津 2026-03-23 | 5/250 |
|