版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

當(dāng)前主題已經(jīng)存檔。

nkleof

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2093.2
散金: 649
紅花: 5
帖子: 813
在線: 165.9小時
蟲號: 152581
注冊: 2006-01-02
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

[交流] 【求助】晶格優(yōu)化問題已有7人參與

我用vasp優(yōu)化一個體系，所得到的結(jié)果8.36A遠小于實驗值8.51A。
所采用的贗勢是PAWPBE的，過渡金屬元素Fe和Mn分別試過標(biāo)準(zhǔn)勢和pv勢。采用的優(yōu)化方法分別試過先優(yōu)化原胞體積（ISIF=7），后馳豫離子位置（ISIF=2）；以及全優(yōu)化（ISIF=3）。馳豫方式分別試過CG法（IBRION=2）和準(zhǔn)牛頓法（IBRION=1）。為了消除pulay stress，每種方法都以得到的CONTCAR作為POSCAR反復(fù)優(yōu)化了3次。
但是以上所有方法所得到的結(jié)果都是8.36附近，遠小于實驗值，也小于其它文獻所給出的GGA計算值8.49。
請經(jīng)驗豐富的蟲友幫忙想想可能是在什么地方出問題了？

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

VASP

» 猜你喜歡

一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑已經(jīng)有4人回復(fù)
286求調(diào)劑已經(jīng)有9人回復(fù)
286分人工智能專業(yè)請求調(diào)劑愿意跨考！已經(jīng)有8人回復(fù)
資源與環(huán)境調(diào)劑申請(333分) 已經(jīng)有5人回復(fù)
280求調(diào)劑已經(jīng)有12人回復(fù)
269專碩求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑院校信息已經(jīng)有3人回復(fù)
材料學(xué)碩301分求調(diào)劑已經(jīng)有7人回復(fù)
初試 317 已經(jīng)有7人回復(fù)
一志愿211，0703化學(xué)310分求調(diào)劑已經(jīng)有3人回復(fù)

1樓 2010-04-19 20:43:21

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

gfei2008

木蟲 (小有名氣)

應(yīng)助: 6 (幼兒園)
金幣: 1716.7
紅花: 1
帖子: 176
在線: 116.2小時
蟲號: 544222
注冊: 2008-04-12
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★ ★
qli0218(金幣+1):謝謝 2010-04-19 21:14
nkleof(金幣+2):感謝建議！ 2010-04-19 21:14
hometownlove(金幣+2):謝謝交流 2010-04-20 15:24

你在優(yōu)化之前要進行K-point和ENCUT的測試。得到一個可靠的K點和ENCUT，然后晶格常數(shù)的求法也有相關(guān)方面的介紹，就是先求出體系能量與晶格之間的關(guān)系然后通過一個方程可以擬合出一條曲線。這個能量最點對應(yīng)的晶格長度就是晶格常數(shù)。還有，我也用過PAWPBE計算過晶格常數(shù)，是要比試驗值稍微小一點。希望我的回答對你有幫助

贊一下(2人)

回復(fù)此樓

2樓2010-04-19 21:11:59

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

nkleof

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2093.2
散金: 649
紅花: 5
帖子: 813
在線: 165.9小時
蟲號: 152581
注冊: 2006-01-02
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by gfei2008 at 2010-04-19 21:11:59:
你在優(yōu)化之前要進行K-point和ENCUT的測試。得到一個可靠的K點和ENCUT，然后晶格常數(shù)的求法也有相關(guān)方面的介紹，就是先求出體系能量與晶格之間的關(guān)系然后通過一個方程可以擬合出一條曲線。這個能量最點對應(yīng)的晶格長 ...

我用的k網(wǎng)格疏密程度是k點之間距離0.09 /A，PREC=Accurate因此ENCUT是贗勢內(nèi)默認值的1.3倍，所以這兩個優(yōu)化條件我覺得精度還可以，如果有問題請經(jīng)驗豐富的高人再給點建議。
擬合的優(yōu)化正打算做，但是我覺得vasp的自動優(yōu)化功能誤差應(yīng)該不會太大，因此用ISIF=2或3的設(shè)置算出來的值和靜態(tài)計算擬合優(yōu)化出來的值應(yīng)該不會有太大差異，是不是有可能我的什么設(shè)置出錯？
如下是我做優(yōu)化用的INCAR

PREC = Accurate
LREAL = A
IALGO = 48
NPAR = 8
NSIM = 8

NELMIN = 4
NSW = 200
LORBIT = 12
EDIFF = 1E-5
IBRION = 2
ISIF = 3

ISPIN = 2
MAGMOM = 8*-6 16*6 32*0

[ Last edited by nkleof on 2010-4-19 at 21:27 ]

贊一下

回復(fù)此樓

3樓2010-04-19 21:23:24

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zdd2008

金蟲 (職業(yè)作家)

應(yīng)助: 1 (幼兒園)
金幣: 5.4
散金: 2613
紅花: 12
帖子: 3845
在線: 108.2小時
蟲號: 708180
注冊: 2009-02-25
專業(yè): 金屬結(jié)構(gòu)材料

★
nkleof(金幣+1):謝謝討論！ 2010-04-20 11:48
xiaohunhun(金幣+1):xiexie 2010-04-20 11:53

MAGMOM = 8*-6 16*6 32*0
你計算出來的磁矩有多少？

贊一下(1人)

回復(fù)此樓

4樓2010-04-20 10:09:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

calos818

木蟲 (著名寫手)

新人菜鳥

應(yīng)助: 11 (小學(xué)生)
金幣: 3125.7
散金: 77
紅花: 3
帖子: 1201
在線: 730.8小時
蟲號: 698507
注冊: 2009-02-08
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
nkleof(金幣+1):謝謝回帖！ 2010-04-20 11:48
xiaohunhun(金幣+1):xiexie 2010-04-20 11:53

引用回帖:

Originally posted by nkleof at 2010-04-19 20:43:21:
我用vasp優(yōu)化一個體系，所得到的結(jié)果8.36A遠小于實驗值8.51A。
所采用的贗勢是PAWPBE的，過渡金屬元素Fe和Mn分別試過標(biāo)準(zhǔn)勢和pv勢。采用的優(yōu)化方法分別試過先優(yōu)化原胞體積（ISIF=7），后馳豫離子位置（ISIF=2）； ...

1.76%的誤差，可以接受了，但還是幫你頂下，這些問題必須去解決，我也很想知道呵呵

贊一下(1人)

回復(fù)此樓

5樓2010-04-20 10:26:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

nkleof

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2093.2
散金: 649
紅花: 5
帖子: 813
在線: 165.9小時
蟲號: 152581
注冊: 2006-01-02
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
xiaohunhun(金幣+1):xiexie 2010-04-20 11:53

引用回帖:

Originally posted by zdd2008 at 2010-04-20 10:09:34:
MAGMOM = 8*-6 16*6 32*0
你計算出來的磁矩有多少？

8*-3.39 16*3.76 32*0.2
這個應(yīng)該比較貼近Mn和Fe的真實磁矩，磁序?qū)Ц癯?shù)的影響應(yīng)該不會很大吧？

[ Last edited by nkleof on 2010-4-20 at 11:49 ]

贊一下(1人)

回復(fù)此樓

6樓2010-04-20 11:47:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zhang668

木蟲 (正式寫手)

1ST強帖: 3
應(yīng)助: 11 (小學(xué)生)
貴賓: 0.061
金幣: 1735
散金: 1294
紅花: 11
帖子: 843
在線: 324.9小時
蟲號: 733242
注冊: 2009-03-28
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

你這個是尖晶石結(jié)構(gòu)吧。我用USPP優(yōu)化的晶格常數(shù)變化也比較大。

贊一下

回復(fù)此樓

7樓2010-04-20 12:25:08

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

zdd2008

金蟲 (職業(yè)作家)

應(yīng)助: 1 (幼兒園)
金幣: 5.4
散金: 2613
紅花: 12
帖子: 3845
在線: 108.2小時
蟲號: 708180
注冊: 2009-02-25
專業(yè): 金屬結(jié)構(gòu)材料

★
gavinliu7390(金幣+1):謝謝交流！ 2010-04-20 15:17

引用回帖:

Originally posted by nkleof at 2010-04-20 11:47:51:

8*-3.39 16*3.76 32*0.2
這個應(yīng)該比較貼近Mn和Fe的真實磁矩，磁序?qū)Ц癯?shù)的影響應(yīng)該不會很大吧？

[ Last edited by nkleof on 2010-4-20 at 11:49 ]

第一原理研究磁性是一個比較復(fù)雜的東西！

這個順序要與POTCAR的順序?qū)?yīng)��！

[ Last edited by zdd2008 on 2010-4-20 at 14:41 ]

贊一下(1人)

回復(fù)此樓

8樓2010-04-20 14:37:49

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

nkleof

銀蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 2093.2
散金: 649
紅花: 5
帖子: 813
在線: 165.9小時
蟲號: 152581
注冊: 2006-01-02
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

引用回帖:

Originally posted by zdd2008 at 2010-04-20 14:37:49:

第一原理研究磁性是一個比較復(fù)雜的東西！

這個順序要與POTCAR的順序?qū)?yīng)！！

[ Last edited by zdd2008 on 2010-4-20 at 14:41 ]

謝謝提醒，這個順序在計算時經(jīng)過仔細檢查了，應(yīng)該是沒有錯的。
我做的這個是尖晶石氧化物的計算，以前也算過很多體系，沒有出現(xiàn)過這樣的問題，這還是第一次遇到優(yōu)化結(jié)果和實驗值相差這么大的。

贊一下

回復(fù)此樓

9樓2010-04-20 16:05:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

fanchen021

鐵桿木蟲 (正式寫手)

應(yīng)助: 7 (幼兒園)
金幣: 7994.9
散金: 1844
紅花: 4
帖子: 776
在線: 569小時
蟲號: 773188
注冊: 2009-05-18
專業(yè): 理論和計算化學(xué)

nkleof(金幣+2):感謝討論！ 2010-04-20 16:47

K點0.09/A 有點少，一般0.05/A。
Ediff 建議1E-7
IALGO 換成38 試試
以上

贊一下

回復(fù)此樓

10樓2010-04-20 16:31:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 nkleof 的主題更新

返回列表

普通表情龍兔虎貓高級回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 資源與環(huán)境調(diào)劑申請(333分) +5	holy J 2026-03-21	5/250	2026-03-21 22:42 by Catalysis25
[考研] 354求調(diào)劑 +6	Tyoumou 2026-03-18	9/450	2026-03-21 20:47 by lbsjt
[考研] 0703化學(xué)調(diào)劑，六級已過，有科研經(jīng)歷 +14	曦熙兮 2026-03-15	14/700	2026-03-21 19:12 by ColorlessPI
[考研] 0703化學(xué)調(diào)劑 +4	妮妮ninicgb 2026-03-21	4/200	2026-03-21 18:39 by 學(xué)員8dgXkO
[考研] 生物學(xué)調(diào)劑 +3	Surekei 2026-03-21	3/150	2026-03-21 18:31 by 學(xué)員8dgXkO
[考研] 278求調(diào)劑 +9	煙火先于春 2026-03-17	9/450	2026-03-21 17:47 by 學(xué)員8dgXkO
[考研] 311求調(diào)劑 +3	勇敢的小吳 2026-03-20	3/150	2026-03-21 17:40 by ColorlessPI
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +5	脫穎而出 2026-03-16	15/750	2026-03-21 10:16 by 脫穎而出
[考研] 求調(diào)劑 +6	Mqqqqqq 2026-03-19	6/300	2026-03-21 08:04 by JourneyLucky
[考研] 08工科 320總分求調(diào)劑 +6	梨花珞晚風(fēng) 2026-03-17	6/300	2026-03-21 03:40 by JourneyLucky
[考研] 308求調(diào)劑 +3	阿姐阿姐家啊 2026-03-18	3/150	2026-03-20 23:24 by JourneyLucky
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +8	小材化本科 2026-03-18	8/400	2026-03-20 23:16 by JourneyLucky
[考研] 一志愿西北大學(xué) ，070300化學(xué)學(xué)碩，總分287，雙非一本，求調(diào)劑。 +4	晨昏線與星海 2026-03-19	4/200	2026-03-20 22:15 by JourneyLucky
[考研] 材料與化工 322求調(diào)劑 +4	然11 2026-03-19	4/200	2026-03-20 22:12 by luoyongfeng
[考研] 北科281學(xué)碩材料求調(diào)劑 +5	tcxiaoxx 2026-03-20	5/250	2026-03-20 21:35 by laoshidan
[考研] 材料學(xué)碩318求調(diào)劑 +5	February_Feb 2026-03-19	5/250	2026-03-19 23:51 by 23Postgrad
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10	Liwangman 2026-03-15	10/500	2026-03-19 10:25 by 無際的草原
[考研] 344求調(diào)劑 +6	knight344 2026-03-16	7/350	2026-03-18 20:13 by walc
[考研] 0854，計算機類招收調(diào)劑 +3	胡辣湯放糖 2026-03-15	6/300	2026-03-18 12:09 by 上岸上岸……..
[考研] 333求調(diào)劑 +3	文思客 2026-03-16	7/350	2026-03-16 18:21 by 文思客

24小時熱門版塊排行榜

[交流] 【求助】晶格優(yōu)化問題 已有7人參與

» 收錄本帖的淘帖專輯推薦

» 猜你喜歡

[交流] 【求助】晶格優(yōu)化問題已有7人參與