| 查看: 2212 | 回復: 19 | ||
luochuiyi90木蟲 (正式寫手)
|
[求助]
MS勾選自旋優(yōu)化失敗
|
|

木蟲 (小有名氣)
|
我這邊最近估計比較忙,你自己看看我給你回復的研究一下吧?重要的地方我都跟你講過了。 另外你看下我給你的附件,我算了幾個例子在里面,initial spin=1,2,5的,以及還有一個是3倍原胞大小的計算,所有計算的結(jié)論都是每個原胞的net spin是1。你比較下你的參數(shù)有沒有設(shè)置的不對的地方 |
木蟲 (小有名氣)
|
castep的默認設(shè)置是,如果檢測到體系總電子數(shù)是奇數(shù),就會自動進行自旋優(yōu)化。 你這個體系是25個電子所以得到的結(jié)果已經(jīng)是考慮了自旋的結(jié)果了。 如果你設(shè)置initial spin,是參考了別人的結(jié)論,可以看下別人的cell是怎么設(shè)置的,至少目前看來這個單胞得到net spin=1還是比較正常的,畢竟總電子數(shù)是奇數(shù)。 建議做下大的晶胞,比如2個,4個單胞大小的,看看這樣設(shè)置后總電子數(shù)變成偶數(shù)后,還會不會有net spin。 另外后面如果需要考慮磁學性質(zhì)的話,可以加上spin orbital coupling看看,vasp有這個功能,castep好像最新的版本有 |
木蟲 (小有名氣)
捐助貴賓 (職業(yè)作家)
木蟲 (正式寫手)
|
謝謝你的回復。之所以設(shè)置為5,是因為看到有的帖子說參數(shù)可以設(shè)置大點,計算的時候會自行調(diào)整。但是,我發(fā)現(xiàn)設(shè)置為2或者1,對應的 2*Integrated Spin Density 與2*Integrated |Spin Density| 值會發(fā)生變化,比如,設(shè)置為5時: 2*Integrated Spin Density = 5.00000 2*Integrated |Spin Density| = 5.00000 ,但是,設(shè)置為2時: 2*Integrated Spin Density = 2.00000 2*Integrated |Spin Density| = 2.10030 設(shè)置為1時:2*Integrated Spin Density = 1.00000 2*Integrated |Spin Density| = 1.14647,不同的值會有不同的結(jié)果。另外LiNiO2中Ni離子為+3價,核外電子排布應該是3d74s0,所以我個人覺得是不是應該在Modify里面將Ni Spin設(shè)置為3呢?期待您的指教。金幣不夠的話我會追加,謝謝! |

木蟲 (正式寫手)

木蟲 (正式寫手)

木蟲 (小有名氣)
|
按道理來說,確實會即使設(shè)置不同的initial spin也都能優(yōu)化到相同的結(jié)果,你看一看算好的.castep文件的末尾有沒有說優(yōu)化成功。然后分別看下自旋向上和自旋向下的電子分別填充在哪個能級上,castep會分別列出來的,自旋向上和自旋向下的電子數(shù)之差就是優(yōu)化好的net spin。 關(guān)于Ni的外層電子,+3價的我覺得排成3d5,4s2似乎更合理,d軌道半滿,s軌道全滿能量應該比較低。不過Ni的金屬性不強,和O的電負性差別沒那么大,所以Ni和O之間可能是共價性的吧,所以Ni和O之間會公用電子吧,這一點我也不太肯定。 另外推薦個網(wǎng)站 https://www.materialsproject.org/ 我剛才看了,上面有每種LiNiO2晶體構(gòu)型人家算的磁性Final Magnetic Moment,不同的晶體結(jié)構(gòu)算出來的磁性也不同,1到5的都有 |
木蟲 (正式寫手)

木蟲 (小有名氣)
|
對,如果是總電子數(shù)是奇數(shù),即使沒勾選spin polarized,castep算的時候也會自動考慮自旋,這是castep的自動糾錯功能,你要證實的話就檢測一下你的*.castep文件,看看里面會有這么一段 **************************** Electronic Parameters **************************** number of electrons : 348.0 net charge of system : 0.000 net spin of system : 0.000 number of up spins : 174.0 number of down spins : 174.0 treating system as non-spin-polarized 比如我這里是non-spin-polarized,就是算的時候沒自旋。 只要這里寫的是 treating system as spin-polarized就是算的時候考慮了自旋。 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 286分調(diào)劑 +20 | Faune 2026-03-30 | 22/1100 |
|
|---|---|---|---|---|
|
[考研] 材料調(diào)劑 +8 | 一樣YWY 2026-04-02 | 8/400 |
|
|
[考研] 求調(diào)劑推薦 +3 | 南山南@ 2026-04-01 | 3/150 |
|
|
[考研] 0710生物學336分求調(diào)劑 +3 | kiyy 2026-04-01 | 3/150 |
|
|
[基金申請] esi高被引論文是不是能對中標有所加分和幫助呢 +4 | redcom 2026-04-01 | 4/200 |
|
|
[考研] 求調(diào)劑,一志愿 南京航空航天大學 ,080500材料科學與工程學碩,總分289分 +11 | @taotao 2026-03-29 | 11/550 |
|
|
[考研] 367求調(diào)劑 +8 | 芋泥啵啵… 2026-03-28 | 8/400 |
|
|
[考研] 279求調(diào)劑 +6 | 學而思兮知 2026-04-01 | 6/300 |
|
|
[考研] 284求調(diào)劑 +12 | 小熊~~ 2026-03-31 | 12/600 |
|
|
[考研] 330分求調(diào)劑 +11 | qzenlc 2026-03-29 | 11/550 |
|
|
[考研] 土木304求調(diào)劑 +3 | 兔突突突, 2026-03-31 | 3/150 |
|
|
[考研] 349求調(diào)劑 +6 | zwjjjjjj 2026-03-31 | 6/300 |
|
|
[考研] 求調(diào)劑 +9 | akdhjs 2026-03-31 | 10/500 |
|
|
[考研] 266分,求材料冶金能源化工等調(diào)劑 +8 | 哇呼哼呼哼 2026-03-27 | 10/500 |
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +4 | 崔wj 2026-03-31 | 4/200 |
|
|
[考研] 南京大學化學調(diào)劑 +11 | 景隨風 2026-03-29 | 16/800 |
|
|
[考研] 一志愿食品科學與工程083200求調(diào)劑 +4 | XQTJZ 2026-03-30 | 4/200 |
|
|
[考研] 哈爾濱工業(yè)大學材料與化工專碩378求調(diào)劑 +3 | 塔比烏斯 2026-03-30 | 3/150 |
|
|
[考研] 本科雙非材料,跨考一志愿華電085801電氣,283求調(diào)劑,任何專業(yè)都可以 +6 | 芝士雪baoo 2026-03-28 | 8/400 |
|
|
[考研] 298調(diào)劑 +3 | jiyingjie123 2026-03-27 | 3/150 |
|