 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 2/1 | 返回列表 |
| 查看: 4465 | 回復: 16 | |||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | |||
[交流]
【求助】如何解決自旋極化計算時的難以收斂問題？ 已有9人參與
|
|||
|
小弟想用castep計算兩種異質界面(Fe/Ni)結合的電子結構和能量學信息， 構建模型的過程如下： 1 .導入程序自帶的Fe的晶胞，幾何優(yōu)化,切出(100)面，構建真空層后幾何優(yōu)化，得出Fe(100)的構型。晶胞參數(shù)a=b=2.775埃 2 .導入程序自帶的Ni的晶胞，幾何優(yōu)化,切出(100)面，構建真空層后幾何優(yōu)化，得出Ni(100)的構型。晶胞參數(shù)c=d=2.513埃 3. 根據(jù)優(yōu)化后的Ni(100)為底，向真空層添加Fe原子，界面間距選擇為優(yōu)化后的兩種結構的層間距的平均值。Fe的層間距采用步驟1的數(shù)據(jù)。以此搭建Fe(100)/Ni(100)結構,晶胞邊距選擇為（a+c)/2，入下圖：  現(xiàn)在的問題是： （1）以上的幾何優(yōu)化都是在加spin polarized下進行的，對這個自己建的結構加自旋極化（自旋設為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），就可以收斂。但這樣的結構不可信的，本來就是磁性體系，況且加自旋極化計算的能量都比不加的低。 是不是我的初始結構不太合理�。炕蚴怯嬎銋�(shù)的問題？希望大家?guī)蛶兔Π桑詈媚芙Y合我的情況給些加速收斂具體的措施。急等！ |
回復此樓» 本帖已獲得的紅花（最新10朵）
歇斯底琳» 猜你喜歡
281求調劑（0805）
已經(jīng)有8人回復
-
環(huán)境領域全國重點實驗室招收博士1-2名
已經(jīng)有3人回復
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復
-
301求調劑
已經(jīng)有6人回復
-
一志愿天津大學化學工藝專業(yè)（081702）315分求調劑
已經(jīng)有7人回復
-
302求調劑
已經(jīng)有6人回復
-
26博士申請
已經(jīng)有3人回復
-
268求調劑
已經(jīng)有3人回復
-
311求調劑
已經(jīng)有10人回復
-
被我言中：新模板不強調格式了，假專家開始管格式了
已經(jīng)有4人回復
16樓2017-09-18 11:49:49
acridine
木蟲 (正式寫手)
- 1ST強帖: 11
- 應助: 13 (小學生)
- 金幣: 4901.6
- 紅花: 13
- 帖子: 706
- 在線: 389.5小時
- 蟲號: 305478
- 注冊: 2006-12-10
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
★ ★ ★
zzy870720z(金幣+3):謝謝指導，很不錯的注意，學習了 2010-06-05 08:32:23
forkn(金幣+6):謝謝你的建議，我試試看 2010-06-06 10:53:22
zzy870720z(金幣+3):謝謝指導，很不錯的注意，學習了 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自洽的收斂情況還是相當漂亮的，這些參數(shù)的設置上不會有太大的問題。那么問題就應該是來自結構的優(yōu)化上。 其實俺一看到真空層+優(yōu)化晶格就很暈，個人認為最好不要這么干。這就是相當于對電腦“下套”，電腦又不象我們這么聰明，它也不知道啥叫真空層，它會費力不討好的不停優(yōu)化c軸長度，而你應該根本就不關心這個方向的晶格長度吧？所以最好不要優(yōu)化，也許不優(yōu)化能很快收斂呢！ 另外如果非極化計算很容易收斂，那么以非極化優(yōu)化出來的結構作為初結構，再進一步的使用自旋極化來優(yōu)化，也許是個不錯的主意。 另外你模型的原子層數(shù)可能太小了，模型的意義可能是個問題，但這和你問的問題沒有什么關系，屬于題外話了。 |
2樓2010-06-05 00:19:47
163.com
金蟲 (著名寫手)
實習版主
- 應助: 2 (幼兒園)
- 貴賓: 0.016
- 金幣: 1489.8
- 散金: 1347
- 紅花: 32
- 帖子: 1215
- 在線: 87小時
- 蟲號: 923817
- 注冊: 2009-12-09
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結構
3樓2010-06-05 06:29:45
4樓2010-06-06 11:31:51
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 考研求調劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 085601專碩，總分342求調劑，地區(qū)不限 +4 | share_joy 2026-03-16 | 4/200 |
|
|
|
[考研] 求調劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學。調劑到任何專業(yè)都可以 +4 | 小丁想進步 2026-03-11 | 5/250 |
|
|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習 +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|
|
|
[考研] 考研調劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
|
[考研] 藥學383 求調劑 +3 | 藥學chy 2026-03-15 | 4/200 |
|
|
|
[考研] 333求調劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
|
[考研] 0854控制工程 359求調劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 0703化學調劑，求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
|
[考研] 0703化學調劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 326求調劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 327求調劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
|
[考研] 材料371求調劑 +9 | 鱷魚? 2026-03-11 | 11/550 |
|
|
|
[考研] 材料專碩288分求調劑 一志愿211 +4 | 在家想你 2026-03-11 | 4/200 |
|
|
|
[考研] 308求調劑 +5 | 是Lupa啊 2026-03-11 | 5/250 |
|
|
|
[考研] 材料專碩350 求調劑 +4 | 王金科 2026-03-12 | 4/200 |
|
|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調劑！ +3 | 不想理你 2026-03-12 | 3/150 |
|
|
|
[考研] 333求調劑 +3 | 152697 2026-03-12 | 4/200 |
|
|
|
[考研] 081200-11408-276學碩求調劑 +3 | 崔wj 2026-03-12 | 4/200 |
|