 |
 |
 |
| 24小時(shí)熱門(mén)版塊排行榜 |
- >論壇更新日志 (1341)
- >考研 (113)
- >導(dǎo)師招生 (103)
- >蟲(chóng)友互識(shí) (51)
- >教師之家 (26)
- >論文道賀祈福 (23)
- >考博 (23)
- >碩博家園 (17)
- >無(wú)機(jī)/物化 (17)
- >論文投稿 (16)
- >文獻(xiàn)求助 (12)
- >休閑灌水 (11)
- >公派出國(guó) (9)
- >基金申請(qǐng) (7)
- >找工作 (5)
- >留學(xué)生活 (4)
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 2/1 | 返回列表 |
| 查看: 4469 | 回復(fù): 16 | |||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||
[交流]
【求助】如何解決自旋極化計(jì)算時(shí)的難以收斂問(wèn)題？ 已有9人參與
|
|||
|
小弟想用castep計(jì)算兩種異質(zhì)界面(Fe/Ni)結(jié)合的電子結(jié)構(gòu)和能量學(xué)信息， 構(gòu)建模型的過(guò)程如下： 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)在的問(wèn)題是： （1）以上的幾何優(yōu)化都是在加spin polarized下進(jìn)行的，對(duì)這個(gè)自己建的結(jié)構(gòu)加自旋極化（自旋設(shè)為low和hight都試過(guò)）進(jìn)行幾何優(yōu)化時(shí)，一直不收斂，我已經(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）我計(jì)算過(guò)不加自旋的優(yōu)化，很容易收斂，energy cut off (300ev), k point set (3*3*2），就可以收斂。但這樣的結(jié)構(gòu)不可信的，本來(lái)就是磁性體系，況且加自旋極化計(jì)算的能量都比不加的低。 是不是我的初始結(jié)構(gòu)不太合理啊？或是計(jì)算參數(shù)的問(wèn)題？希望大家?guī)蛶兔Π�，最好能結(jié)合我的情況給些加速收斂具體的措施。急等！ |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
歇斯底琳» 猜你喜歡
293求調(diào)劑
已經(jīng)有14人回復(fù)
-
0703化學(xué)336分求調(diào)劑
已經(jīng)有6人回復(fù)
-
304求調(diào)劑
已經(jīng)有9人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有8人回復(fù)
-
環(huán)境領(lǐng)域全國(guó)重點(diǎn)實(shí)驗(yàn)室招收博士1-2名
已經(jīng)有3人回復(fù)
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復(fù)
-
301求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復(fù)
-
302求調(diào)劑
已經(jīng)有6人回復(fù)
-
26博士申請(qǐng)
已經(jīng)有3人回復(fù)
297186690
木蟲(chóng) (正式寫(xiě)手)
- 應(yīng)助: 3 (幼兒園)
- 金幣: 1889.6
- 紅花: 1
- 帖子: 364
- 在線: 59.4小時(shí)
- 蟲(chóng)號(hào): 1994420
- 注冊(cè): 2012-09-12
- 專業(yè): 結(jié)構(gòu)化學(xué)
15樓2013-10-31 10:58:40
acridine
木蟲(chóng) (正式寫(xiě)手)
- 1ST強(qiáng)帖: 11
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 4901.6
- 紅花: 13
- 帖子: 706
- 在線: 389.5小時(shí)
- 蟲(chóng)號(hào): 305478
- 注冊(cè): 2006-12-10
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★ ★ ★
zzy870720z(金幣+3):謝謝指導(dǎo)，很不錯(cuò)的注意，學(xué)習(xí)了 2010-06-05 08:32:23
forkn(金幣+6):謝謝你的建議，我試試看 2010-06-06 10:53:22
zzy870720z(金幣+3):謝謝指導(dǎo)，很不錯(cuò)的注意，學(xué)習(xí)了 2010-06-05 08:32:23
forkn(金幣+6):謝謝你的建議，我試試看 2010-06-06 10:53:22
|
你提到的energy cut off , k point set等參數(shù)，確實(shí)可以影響收斂，但他們最主要還是影響電子scf的收斂。從你的優(yōu)化過(guò)程看，已經(jīng)走了20多步了，所以說(shuō)scf自洽的收斂情況還是相當(dāng)漂亮的，這些參數(shù)的設(shè)置上不會(huì)有太大的問(wèn)題。那么問(wèn)題就應(yīng)該是來(lái)自結(jié)構(gòu)的優(yōu)化上。 其實(shí)俺一看到真空層+優(yōu)化晶格就很暈，個(gè)人認(rèn)為最好不要這么干。這就是相當(dāng)于對(duì)電腦“下套”，電腦又不象我們這么聰明，它也不知道啥叫真空層，它會(huì)費(fèi)力不討好的不停優(yōu)化c軸長(zhǎng)度，而你應(yīng)該根本就不關(guān)心這個(gè)方向的晶格長(zhǎng)度吧？所以最好不要優(yōu)化，也許不優(yōu)化能很快收斂呢！ 另外如果非極化計(jì)算很容易收斂，那么以非極化優(yōu)化出來(lái)的結(jié)構(gòu)作為初結(jié)構(gòu)，再進(jìn)一步的使用自旋極化來(lái)優(yōu)化，也許是個(gè)不錯(cuò)的主意。 另外你模型的原子層數(shù)可能太小了，模型的意義可能是個(gè)問(wèn)題，但這和你問(wèn)的問(wèn)題沒(méi)有什么關(guān)系，屬于題外話了。 |
2樓2010-06-05 00:19:47
163.com
金蟲(chóng) (著名寫(xiě)手)
實(shí)習(xí)版主
- 應(yīng)助: 2 (幼兒園)
- 貴賓: 0.016
- 金幣: 1489.8
- 散金: 1347
- 紅花: 32
- 帖子: 1215
- 在線: 87小時(shí)
- 蟲(chóng)號(hào): 923817
- 注冊(cè): 2009-12-09
- 性別: GG
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
★
zzy870720z(金幣+1):謝謝參與，最好能拿出證據(jù)再說(shuō)，你做過(guò)的測(cè)試？ 2010-06-05 08:33:36
forkn(金幣+3):謝謝交流 2010-06-06 11:32:43
zzy870720z(金幣+1):謝謝參與，最好能拿出證據(jù)再說(shuō)，你做過(guò)的測(cè)試？ 2010-06-05 08:33:36
forkn(金幣+3):謝謝交流 2010-06-06 11:32:43
|
你的體系不收斂是體系結(jié)構(gòu)引起的，你的參數(shù)設(shè)置這樣沒(méi)有很大不妥，沒(méi)問(wèn)題； 你如果做過(guò)不加自旋優(yōu)化，可以用之前優(yōu)化模型加 自旋計(jì)算； 第一性原理計(jì)算指出，加不加自旋對(duì)對(duì)體系的結(jié)構(gòu)變化不大。 |
3樓2010-06-05 06:29:45
4樓2010-06-06 11:31:51
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 303求調(diào)劑 +3 | 睿08 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +4 | 錦鯉0909 2026-03-17 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
|
[考研] 293求調(diào)劑 +6 | 世界首富 2026-03-11 | 6/300 |
|
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
|
[考研] 一志愿蘇州大學(xué)材料工程（085601）專碩有科研經(jīng)歷三項(xiàng)國(guó)獎(jiǎng)兩個(gè)實(shí)用型專利一項(xiàng)省級(jí)立項(xiàng) +6 | 大火山小火山 2026-03-16 | 8/400 |
|
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學(xué)獎(jiǎng)覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +4 | ahbd 2026-03-14 | 4/200 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考博] 東華理工大學(xué)化材專業(yè)26屆碩士博士申請(qǐng) +6 | zlingli 2026-03-13 | 6/300 |
|
|
|
[考研] 復(fù)試調(diào)劑 +3 | 呼呼？~+123456 2026-03-14 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
|
[考研] 材料與化工（0856）304求B區(qū)調(diào)劑 +6 | 邱gl 2026-03-12 | 7/350 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
|
[考研] 277求調(diào)劑 +4 | anchor17 2026-03-12 | 4/200 |
|
|
|
[考研] 283求調(diào)劑，材料、化工皆可 +8 | 蘇打水7777 2026-03-11 | 10/500 |
|
|
|
[考博] 26讀博 +4 | Rui135246 2026-03-12 | 10/500 |
|