 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 4459 | 回復(fù): 16 | ||||||
[交流]
【求助】如何解決自旋極化計算時的難以收斂問題？ 已有9人參與
|
||||||
|
小弟想用castep計算兩種異質(zhì)界面(Fe/Ni)結(jié)合的電子結(jié)構(gòu)和能量學(xué)信息， 構(gòu)建模型的過程如下： 1 .導(dǎo)入程序自帶的Fe的晶胞，幾何優(yōu)化,切出(100)面，構(gòu)建真空層后幾何優(yōu)化，得出Fe(100)的構(gòu)型。晶胞參數(shù)a=b=2.775埃 2 .導(dǎo)入程序自帶的Ni的晶胞，幾何優(yōu)化,切出(100)面，構(gòu)建真空層后幾何優(yōu)化，得出Ni(100)的構(gòu)型。晶胞參數(shù)c=d=2.513埃 3. 根據(jù)優(yōu)化后的Ni(100)為底，向真空層添加Fe原子，界面間距選擇為優(yōu)化后的兩種結(jié)構(gòu)的層間距的平均值。Fe的層間距采用步驟1的數(shù)據(jù)。以此搭建Fe(100)/Ni(100)結(jié)構(gòu),晶胞邊距選擇為（a+c)/2，入下圖：  現(xiàn)在的問題是： （1）以上的幾何優(yōu)化都是在加spin polarized下進行的，對這個自己建的結(jié)構(gòu)加自旋極化（自旋設(shè)為low和hight都試過）進行幾何優(yōu)化時，一直不收斂，我已經(jīng)嘗試了增加energy cut off (350ev), 增加k point set (7*7*2），empty bonds, 降低SCF 精度，降低smearing等，都不收斂，總是提示： *Warning* max. SCF cycles performed but system has not reached the groundstate. Warning: electronic minimisation did not converge when finding ground state 收斂曲線：  能量曲線：  （2）我計算過不加自旋的優(yōu)化，很容易收斂，energy cut off (300ev), k point set (3*3*2），就可以收斂。但這樣的結(jié)構(gòu)不可信的，本來就是磁性體系，況且加自旋極化計算的能量都比不加的低。 是不是我的初始結(jié)構(gòu)不太合理��？或是計算參數(shù)的問題？希望大家?guī)蛶兔Π�，最好能結(jié)合我的情況給些加速收斂具體的措施。急等！ |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 第一原理資料匯編 | siesta問題 | 心理互助 |
» 本帖已獲得的紅花（最新10朵）
歇斯底琳» 猜你喜歡
311求調(diào)劑
已經(jīng)有9人回復(fù)
-
302求調(diào)劑
已經(jīng)有4人回復(fù)
-
被我言中：新模板不強調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復(fù)
-
303求調(diào)劑
已經(jīng)有5人回復(fù)
-
考研求調(diào)劑
已經(jīng)有4人回復(fù)
-
085601專碩，總分342求調(diào)劑，地區(qū)不限
已經(jīng)有4人回復(fù)
-
085601材料工程專碩求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有6人回復(fù)
-
能源材料化學(xué)課題組招收碩士研究生8-10名
已經(jīng)有6人回復(fù)
-
296求調(diào)劑
已經(jīng)有5人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- castep 自旋極化問題
已經(jīng)有10人回復(fù)
- 為什么我加自旋極化算出來的態(tài)密度用gnuplot畫出來只有自旋向上態(tài)而沒有自旋向下的態(tài)
已經(jīng)有18人回復(fù)
- 關(guān)于自旋極化的設(shè)置
已經(jīng)有4人回復(fù)
- 摻雜TiO2 自旋極化
已經(jīng)有4人回復(fù)
- 【求助】自旋極化出錯
已經(jīng)有9人回復(fù)
- 【求助】關(guān)于擬合得到體彈性模量&自旋極化&LDA+U等問題，望大家指點指點
已經(jīng)有4人回復(fù)
- 【求助】VASP怎么考慮自旋極化？
已經(jīng)有14人回復(fù)
- 【求助】自旋極化計算設(shè)置的問題
已經(jīng)有3人回復(fù)
- 【求助】什么樣的圖形能說明原子發(fā)生了自旋極化...？
已經(jīng)有3人回復(fù)
- 【求助】自旋極化計算
已經(jīng)有4人回復(fù)
- 【前沿】電子自旋極化首次在室溫下實現(xiàn)
已經(jīng)有4人回復(fù)
acridine
木蟲 (正式寫手)
- 1ST強帖: 11
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 4901.6
- 紅花: 13
- 帖子: 706
- 在線: 389.5小時
- 蟲號: 305478
- 注冊: 2006-12-10
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★ ★ ★
zzy870720z(金幣+3):謝謝指導(dǎo)，很不錯的注意，學(xué)習(xí)了 2010-06-05 08:32:23
forkn(金幣+6):謝謝你的建議，我試試看 2010-06-06 10:53:22
zzy870720z(金幣+3):謝謝指導(dǎo)，很不錯的注意，學(xué)習(xí)了 2010-06-05 08:32:23
forkn(金幣+6):謝謝你的建議，我試試看 2010-06-06 10:53:22
|
你提到的energy cut off , k point set等參數(shù)，確實可以影響收斂，但他們最主要還是影響電子scf的收斂。從你的優(yōu)化過程看，已經(jīng)走了20多步了，所以說scf自洽的收斂情況還是相當(dāng)漂亮的，這些參數(shù)的設(shè)置上不會有太大的問題。那么問題就應(yīng)該是來自結(jié)構(gòu)的優(yōu)化上。 其實俺一看到真空層+優(yōu)化晶格就很暈，個人認(rèn)為最好不要這么干。這就是相當(dāng)于對電腦“下套”，電腦又不象我們這么聰明，它也不知道啥叫真空層，它會費力不討好的不停優(yōu)化c軸長度，而你應(yīng)該根本就不關(guān)心這個方向的晶格長度吧？所以最好不要優(yōu)化，也許不優(yōu)化能很快收斂呢！ 另外如果非極化計算很容易收斂，那么以非極化優(yōu)化出來的結(jié)構(gòu)作為初結(jié)構(gòu)，再進一步的使用自旋極化來優(yōu)化，也許是個不錯的主意。 另外你模型的原子層數(shù)可能太小了，模型的意義可能是個問題，但這和你問的問題沒有什么關(guān)系，屬于題外話了。 |
2樓2010-06-05 00:19:47
163.com
金蟲 (著名寫手)
實習(xí)版主
- 應(yīng)助: 2 (幼兒園)
- 貴賓: 0.016
- 金幣: 1489.8
- 散金: 1347
- 紅花: 32
- 帖子: 1215
- 在線: 87小時
- 蟲號: 923817
- 注冊: 2009-12-09
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
3樓2010-06-05 06:29:45
★ ★ ★
xiaohunhun(金幣+3):謝謝 2010-06-08 23:49:18
xiaohunhun(金幣+3):謝謝 2010-06-08 23:49:18
|
引用cenwanglai的測試結(jié)果： “對于含有未配對電子的結(jié)構(gòu),計算能量時是否考慮自旋對結(jié)果影響不可忽略,而且最好在結(jié)構(gòu)優(yōu)化時就考慮自旋最好;但是作為計算時間成本的考慮,在結(jié)構(gòu)優(yōu)化時不考慮自旋,而在后續(xù)能量計算時考慮自旋可能對能量精度影響不大;但是是否在結(jié)構(gòu)優(yōu)化時考慮自旋對結(jié)構(gòu)有一定影響.” 詳見：http://www.gaoyang168.com/bbs/viewthread.php?tid=1532270&fpage=0&view=&highlight=&page=2 本人覺得還是結(jié)構(gòu)優(yōu)化時就加上自旋比較好。至于先不加自旋優(yōu)化，再此結(jié)構(gòu)上再加自旋優(yōu)化的方法，我正在試，希望能收斂吧,呵呵。 |
8樓2010-06-07 09:57:47
zhenfeili
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 73.6
- 散金: 1
- 帖子: 144
- 在線: 57.8小時
- 蟲號: 922736
- 注冊: 2009-12-08
- 專業(yè): 生物電子學(xué)與生物信息處理
9樓2010-06-07 21:03:52
4樓2010-06-06 11:31:51
5樓2010-06-06 11:48:59
6樓2010-06-06 17:51:07
zhenfeili
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 73.6
- 散金: 1
- 帖子: 144
- 在線: 57.8小時
- 蟲號: 922736
- 注冊: 2009-12-08
- 專業(yè): 生物電子學(xué)與生物信息處理
7樓2010-06-06 21:33:06
10樓2010-06-08 10:03:47
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 forkn 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 311求調(diào)劑 +8 | 冬十三 2026-03-15 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
|
[考研] 312求調(diào)劑 +4 | 陌宸希 2026-03-16 | 5/250 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|
|
|
[考研] 085600材料與化工 +4 | 安全上岸！ 2026-03-16 | 4/200 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|
|
|
[考研] 285化工學(xué)碩求調(diào)劑（081700） +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
|
[考研] 機械專碩325，尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|
|
|
[基金申請]
有必要更換申報口嗎
 20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
|
[考研] 295求調(diào)劑 +3 | 小匕仔汁 2026-03-12 | 3/150 |
|
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123！ 2026-03-11 | 4/200 |
|
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|