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

登錄注冊

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

24小時(shí)熱門版塊排行榜

返回列表

luxia0210

鐵桿木蟲 (著名寫手)

應(yīng)助: 1 (幼兒園)
金幣: 6476.2
紅花: 6
帖子: 1497
在線: 329.7小時(shí)
蟲號(hào): 877190
注冊: 2009-10-19
專業(yè): 材料物理化學(xué)

★
zmrright(金幣+1):謝謝參與

，在賺個(gè)金幣哈

回復(fù)此樓

潛心百年，為求一敗！

11樓2010-05-04 18:50:39

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

acridine

木蟲 (正式寫手)

1ST強(qiáng)帖: 11
應(yīng)助: 13 (小學(xué)生)
金幣: 4901.6
紅花: 13
帖子: 706
在線: 389.5小時(shí)
蟲號(hào): 305478
注冊: 2006-12-10
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★ ★ ★ ★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流
xiaohunhun(金幣+3):xiexie 2010-05-04 20:55:00

樓主的語序有點(diǎn)亂啊，似乎很匆忙。
>1. 我在算LiBH4,但是我的k-point 總是弄不好，我想應(yīng)該是與晶體的晶格參數(shù)有關(guān)，是不是啊
k的設(shè)置不僅僅和晶格參數(shù)有關(guān)，與你計(jì)算的材料、要討論性質(zhì)等等都有關(guān)系。準(zhǔn)確的設(shè)置k-piont就如同樓上各位提到的“收斂性測試”

>2. 我在 run的時(shí)候總是報(bào)錯(cuò)，我想應(yīng)該是k-point 出錯(cuò)了。我改變?yōu)樾〉木秃昧恕�。。但是咋確定呢。。。？
這里給出的信息太少，在給出錯(cuò)誤信息之前，確實(shí)無法確定是否是k點(diǎn)設(shè)置的問題。比如內(nèi)存不足也會(huì)出現(xiàn)這種現(xiàn)象：高精度出錯(cuò)，降低精度才可以計(jì)算下去。

贊一下(2人)

回復(fù)此樓

12樓2010-05-04 19:36:27

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

awmc2008

至尊木蟲 (文壇精英)

學(xué)習(xí)者

應(yīng)助: 62 (初中生)
金幣: 13680.6
紅花: 10
帖子: 12260
在線: 736.8小時(shí)
蟲號(hào): 600677
注冊: 2008-09-12
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流

學(xué)習(xí)學(xué)習(xí)，謝謝高手。

贊一下(1人)

回復(fù)此樓

xuexizhe

13樓2010-05-04 22:48:14

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

zmrright

金蟲 (著名寫手)

應(yīng)助: 49 (小學(xué)生)
金幣: 3459.8
散金: 706
紅花: 9
帖子: 1975
在線: 299.1小時(shí)
蟲號(hào): 828781
注冊: 2009-08-16
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

謝謝，我搜一下。。。

引用回帖:

Originally posted by guifan at 2010-05-03 22:13:50:
忘了提醒一下，可以關(guān)注金屬所的胡青苗教授，直接檢索他的文章就能發(fā)現(xiàn)

謝謝你好。。。

回復(fù)此樓

做好自我，不枉此生

14樓2010-05-04 23:21:44

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

zmrright

金蟲 (著名寫手)

應(yīng)助: 49 (小學(xué)生)
金幣: 3459.8
散金: 706
紅花: 9
帖子: 1975
在線: 299.1小時(shí)
蟲號(hào): 828781
注冊: 2009-08-16
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

謝謝

引用回帖:

Originally posted by acridine at 2010-05-04 19:36:27:
樓主的語序有點(diǎn)亂啊，似乎很匆忙。
>1. 我在算LiBH4,但是我的k-point 總是弄不好，我想應(yīng)該是與晶體的晶格參數(shù)有關(guān)，是不是啊
k的設(shè)置不僅僅和晶格參數(shù)有關(guān)，與你計(jì)算的材料、要討論性質(zhì)等等都有關(guān)系。準(zhǔn)確的 ...

我嘗試一下。。。謝謝

回復(fù)此樓

做好自我，不枉此生

15樓2010-05-04 23:23:22

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

zmrright

金蟲 (著名寫手)

應(yīng)助: 49 (小學(xué)生)
金幣: 3459.8
散金: 706
紅花: 9
帖子: 1975
在線: 299.1小時(shí)
蟲號(hào): 828781
注冊: 2009-08-16
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

剛搜集的，雖然沒太懂，但是大家商量看。。。

★ ★
aylayl08(金幣+2):鼓勵(lì)一下討論 2010-05-07 19:55:52

大家好：

我根據(jù)自己在計(jì)算過程中遇到的收斂性的問題，做了一點(diǎn)關(guān)于CASTEP計(jì)算參數(shù)調(diào)試的總結(jié)。當(dāng)然，這些只是我的個(gè)人心得體會(huì)，大家如果覺得有什么問題請?zhí)岢鲋刚�，或是如果有自己�?dú)到的見解，請跟帖加上來，供大家一起學(xué)習(xí)提高。

如果計(jì)算在Max. cycle全部用完之后仍然不能收斂，要先查看一下.castep文件，如果能量比剛開始計(jì)算時(shí)有明顯下降，而且震動(dòng)幅度也明顯減小，只是還沒有降到收斂標(biāo)準(zhǔn)以下，那說明Max. cycle還不夠，要增加這個(gè)數(shù)值。

但多數(shù)情況下，我們會(huì)發(fā)現(xiàn)，能量并沒有比剛開始計(jì)算時(shí)的能量有明顯下降，而是圍繞著某個(gè)平均值在反復(fù)地震蕩，而且震動(dòng)幅度也沒有減小。那說明是參數(shù)設(shè)定的問題，此時(shí)就算再怎么增加Max. cycle的數(shù)值，計(jì)算最終還是不會(huì)收斂的。

這個(gè)時(shí)候，最好是增加energy cutoff的數(shù)值，降低K point set的數(shù)值，增加smearing width數(shù)值（當(dāng)然，如果通過調(diào)試前兩個(gè)數(shù)值就可以使計(jì)算收斂的話，就最好不要?jiǎng)舆@個(gè)數(shù)值）。

再有，如果計(jì)算是使用的density mixing的話，那要將mix_charge_amp，mix_spin_amp的數(shù)值降低，將mix_energy_cutoff的數(shù)值增加到energy cutoff數(shù)值的3到4倍，增加mix_history_length的數(shù)值。

當(dāng)然，如果調(diào)試上面的參數(shù)仍然發(fā)現(xiàn)收斂很辛苦的話，就將mixing scheme轉(zhuǎn)成All Band/EDFT,這會(huì)解決收斂問題，但是計(jì)算時(shí)間將會(huì)是使用density mixing計(jì)算所用時(shí)間的3倍以上。

總結(jié)的挺好的，但是指出幾點(diǎn)不足，即：

1. 參數(shù)調(diào)試過于經(jīng)驗(yàn)化，如果樓主能結(jié)合對程序原理的認(rèn)識(shí)，或許能夠把每個(gè)參數(shù)調(diào)節(jié)的意義給闡述的更為清晰些；從優(yōu)化算法，體系特性，物理過程等方面著手來進(jìn)行分析，或許可以得到更加全面的認(rèn)識(shí)（如果樓主能夠結(jié)合上述內(nèi)容對CASTEP收斂進(jìn)行詳細(xì)的分析和總結(jié)的話，我會(huì)考慮接受樓主的文章，并且發(fā)表在小木蟲計(jì)算與模擬研究RCS期刊上）。

2. 建議樓主在進(jìn)行平面波DFT優(yōu)化之前，通過其他的辦法，比如經(jīng)典力學(xué)優(yōu)化以及DMol中局域波函數(shù)進(jìn)行優(yōu)化，這樣的話，可以幫助CASTEP在最短的時(shí)間內(nèi)找到合適的最小值，也就是通常意義上說的加快收斂。

great!
其實(shí)在做收斂速度無非是兩個(gè)目的：１，計(jì)算精度要高；２計(jì)算速度要快．
一般平面波程序，在計(jì)算精度上，主要需要測試ecut(倒空間的最大平面波截?cái)�，理論上要求無窮大，因?yàn)槟菢踊覆攀峭陚涞模珜?shí)際上不需要，也做不到)．還有就是ｋ點(diǎn)取樣．理論上要計(jì)算無窮個(gè)ｋ點(diǎn)，但實(shí)際也不可能．我覺得應(yīng)該先測試ecut,因?yàn)槟闶紫纫ＷC每個(gè)ｋ點(diǎn)能量是正確的．然后再測試ｋ點(diǎn)．其實(shí)到了，最后即使ｋ點(diǎn)取的很大的時(shí)候，還是有很大的誤差的，這主要來自于系統(tǒng)誤差．所以有時(shí)候，會(huì)造成震蕩的情況，這主要還是由于ｅｃｕｔ太小，或者是有的ｋ點(diǎn)根本就難以算準(zhǔn)．這個(gè)時(shí)候，如無特殊需要，基本就認(rèn)為它收斂了．除非你想要在０．１ｍｅＶ的誤差．　
但是根據(jù)我的想法，我就不認(rèn)為它能將計(jì)算精度控制在０．１ｍｅＶ．即使收斂測試達(dá)到了．（當(dāng)然,這是個(gè)人看法�。�
第二，就是計(jì)算速度．因?yàn)槟壳埃模疲灾饕€是用自洽的方法，在理論上應(yīng)該是用新的電荷密度等于舊的電荷密度，但實(shí)際上如果這樣做，一是收斂很慢，二是容易形成震蕩情形。一般的做法是rho_in=mix*tho_old+(1-mix)*rho_old。（具體到每個(gè)程序，肯定是做法不同了。）一般認(rèn)為mix=0.3比較合理。這樣收斂的非�？臁！　∷院芏喑绦蛟跍y試各種mix的時(shí)候，就是在測試這個(gè)。
最后，選擇合適的參數(shù)。在你所需要的精度范圍內(nèi)，最快的計(jì)算�！∥覀兊哪康氖且鉀Q問題，有時(shí)候不需要那么精確。

引用：“圍繞著某個(gè)平均值在反復(fù)地震蕩，而且震動(dòng)幅度也沒有減小。
這個(gè)時(shí)候，最好是增加energy cutoff的數(shù)值，降低K point set的數(shù)值，增加smearing width數(shù)值（當(dāng)然，如果通過調(diào)試前兩個(gè)數(shù)值就可以使計(jì)算收斂的話，就最好不要?jiǎng)舆@個(gè)數(shù)值）�！�

我覺得不是減�。�，而應(yīng)當(dāng)增加Ｋ
如果不是金屬體系，調(diào)節(jié)ｓｍｅａｒｉｎｇ沒有意義
如有不妥，還望不吝賜教：Ｐ

贊一下(1人)

回復(fù)此樓

做好自我，不枉此生

16樓2010-05-05 00:03:57

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

wzh上善若水

木蟲 (著名寫手)

應(yīng)助: 0 (幼兒園)
金幣: 7522.5
散金: 103
帖子: 1378
在線: 101.5小時(shí)
蟲號(hào): 980865
注冊: 2010-03-24
性別: GG
專業(yè): 粒子物理學(xué)和場論

★
小木蟲(金幣+0.5):給個(gè)紅包，謝謝回帖交流

需要進(jìn)行進(jìn)行收斂測試，可查看相關(guān)文章學(xué)習(xí)收斂性測試

贊一下(1人)

回復(fù)此樓

知止而后有定；定而后能靜；靜而后能安；安而后能慮；慮而后能得

17樓2010-05-09 20:38:02

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

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

返回列表

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

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

[考研] 本人考085602 化學(xué)工程專碩 +20	不知道叫什么！ 2026-03-15	22/1100	2026-03-21 19:03 by ColorlessPI
[考研] 生物學(xué)調(diào)劑 +3	Surekei 2026-03-21	3/150	2026-03-21 18:31 by 學(xué)員8dgXkO
[考研] 268求調(diào)劑 +9	簡單點(diǎn)0 2026-03-17	9/450	2026-03-21 15:37 by lature00
[考研] 初始318分求調(diào)劑（有工作經(jīng)驗(yàn)） +3	1911236844 2026-03-17	3/150	2026-03-21 02:33 by JourneyLucky
[考研] 材料 336 求調(diào)劑 +3	An@. 2026-03-18	4/200	2026-03-21 01:39 by JourneyLucky
[考研] 324分 085600材料化工求調(diào)劑 +4	llllkkkhh 2026-03-18	4/200	2026-03-21 01:24 by JourneyLucky
[考研] 311求調(diào)劑 +5	冬十三 2026-03-18	5/250	2026-03-21 00:16 by JourneyLucky
[考研] 南京大學(xué)化學(xué)376求調(diào)劑 +3	hisfailed 2026-03-19	6/300	2026-03-20 23:43 by hisfailed
[考研] A區(qū)線材料學(xué)調(diào)劑 +5	周周無極 2026-03-20	5/250	2026-03-20 21:33 by laoshidan
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +3	困困困困坤坤 2026-03-20	3/150	2026-03-20 20:38 by 學(xué)員8dgXkO
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5	1孫悟空 2026-03-17	6/300	2026-03-19 23:46 by zcl123
[考博] 申博26年 +3	八6八68 2026-03-19	3/150	2026-03-19 19:43 by nxgogo
[考研] 286求調(diào)劑 +6	lemonzzn 2026-03-16	10/500	2026-03-19 14:31 by lemonzzn
[考研] 材料考研調(diào)劑 +3	xwt。 2026-03-19	3/150	2026-03-19 11:22 by w沐陽w
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 生物學(xué)071000 329分求調(diào)劑 +3	我愛生物生物愛?/a> 2026-03-17	3/150	2026-03-18 10:12 by macy2011
[考研] 0703化學(xué)336分求調(diào)劑 +6	zbzihdhd 2026-03-15	7/350	2026-03-18 09:53 by zhukairuo
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 302求調(diào)劑 +4	小賈同學(xué)123 2026-03-15	8/400	2026-03-17 10:33 by 小賈同學(xué)123
[考研] 0856求調(diào)劑 +3	劉夢微 2026-03-15	3/150	2026-03-16 10:00 by houyaoxu

24小時(shí)熱門版塊排行榜

謝謝，我搜一下。。。

謝謝

剛搜集的，雖然沒太懂，但是大家商量看。。。

謝謝，我搜一下。。。

剛搜集的，雖然沒太懂，但是大家商量看。。。