| 5 | 1/1 | 返回列表 |
| 查看: 3601 | 回復(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)閷?duì)比文獻(xiàn),磁矩該大的地方小,該小的地方大。希望大神出來解救小弟,我已經(jīng)搗鼓一天了,試了N多種情況,求解答,不勝感激。 |
木蟲 (正式寫手)
木蟲 (正式寫手)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +3 | 彼岸t 2026-03-29 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 298求調(diào)劑 +3 | 種圣賜 2026-03-29 | 3/150 |
|
|
[考研] 356求調(diào)劑 +4 | gysy?s?a 2026-03-28 | 4/200 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 生物學(xué)學(xué)碩,一志愿湖南大學(xué),初試成績(jī)338 +6 | YYYYYNNNNN 2026-03-26 | 7/350 |
|
|
[考研] 291求調(diào)劑 +15 | hhhhxn.. 2026-03-23 | 21/1050 |
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢(mèng) 2026-03-23 | 16/800 |
|
|
[考研] 314求調(diào)劑 +3 | 溪云珂 2026-03-26 | 3/150 |
|
|
[考研] 考研調(diào)劑 +9 | 小蠟新筆 2026-03-26 | 9/450 |
|
|
[考研] 一志愿吉大071010,316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|
|
[考研] 0703化學(xué)338求調(diào)劑! +6 | Zuhui0306 2026-03-26 | 7/350 |
|
|
[考研] 機(jī)械學(xué)碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 一志愿 南京郵電大學(xué) 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 材料與化工328分調(diào)劑 +6 | 。,。,。,。i 2026-03-23 | 6/300 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +3 | 邱gl 2026-03-25 | 3/150 |
|
|
[考研] 302求調(diào)劑 +4 | 錦衣衛(wèi)藤椒 2026-03-25 | 4/200 |
|
|
[有機(jī)交流]
20+3
|
FENGSHUJEI 2026-03-23 | 5/250 |
|
|
[考研] 344求調(diào)劑 +3 | desto 2026-03-24 | 3/150 |
|
|
[考研] 284求調(diào)劑 +3 | yanzhixue111 2026-03-23 | 6/300 |
|