 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 10 | 1/1 | 返回列表 |
| 查看: 2531 | 回復(fù): 9 | ||||
yeahgyz銅蟲 (初入文壇)
|
[求助]
黃銅礦幾何結(jié)構(gòu)優(yōu)化 按照文獻(xiàn)參數(shù)不能收斂，且晶胞常數(shù)誤差過大如何解決 已有1人參與
|
|
小弟按照文獻(xiàn)所述參數(shù)對黃銅礦進(jìn)行了幾何結(jié)構(gòu)的優(yōu)化，但是卻一直不能收斂是怎么回事，而且算得的晶胞常數(shù)與實(shí)驗(yàn)值相差超過了2%，甚至達(dá)到了4%或5%，是什么地方出現(xiàn)了問題嗎？下面吧具體模擬過程及文獻(xiàn)提供給各位大神分析，謝謝！ 文獻(xiàn)中的設(shè)置如下：黃銅礦，理想黃銅礦的晶胞參數(shù)為：a=b=5.289 Å，c=10.423 Å，α=β=γ=90°。采用GG-PBE梯度修正函數(shù)，在進(jìn)行幾何優(yōu)化時(shí)，選取合適的截?cái)嗄芎蚹點(diǎn)，收斂標(biāo)準(zhǔn)為：原子位移為0.005 Å，原子間作用力為0.01 eV/nm，原子間的內(nèi)應(yīng)力為0.2 GPa，體系總能量的變化不大于5.0×10-4 eV/atom。 截?cái)嗄芎蚹點(diǎn)選擇過程及數(shù)據(jù)如圖一和圖二： 我按照245計(jì)算，計(jì)算到98步進(jìn)行不下去了，受不了就停止了。 手頭上有一個(gè)270ev的數(shù)據(jù)，收斂過程如圖三圖四，設(shè)置參數(shù)如圖五圖六。 另外我算了E點(diǎn)為310，350和400的，k點(diǎn)位5x5x2，很快收斂，但是優(yōu)化后的晶胞常數(shù)誤差超過2%，最大的達(dá)到4%，應(yīng)該是不能用的。 請問1.為什么我按照文獻(xiàn)設(shè)置不能收斂？ 2.為什么提高E點(diǎn)可以收斂，但是晶胞常數(shù)誤差過大？是不是K點(diǎn)設(shè)置問題，還是 其他問題？ 謝謝指教！ 360截圖20140718223046593.jpg 360截圖20140718223114573.jpg 360截圖20140718223437784.jpg 360截圖20140718223452372.jpg 360截圖20140718223536488.jpg 360截圖20140718223610451.jpg |
回復(fù)此樓» 猜你喜歡
271材料工程求調(diào)劑
已經(jīng)有5人回復(fù)
-
281求調(diào)劑（0805）
已經(jīng)有16人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
材料工程專碩調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天大材料與化工（085600）總分338
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有3人回復(fù)
-
求材料調(diào)劑
已經(jīng)有8人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有5人回復(fù)
-
一志愿華中科技大學(xué)，080502，354分求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有6人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- CASTEP幾何結(jié)構(gòu)優(yōu)化的問題 雖然收斂了 感覺不是很對 請大家?guī)兔φ乙幌略?/a>
已經(jīng)有11人回復(fù)
- 收斂測試和結(jié)構(gòu)優(yōu)化
已經(jīng)有10人回復(fù)
- 【討論】MS中sorption模塊的MC計(jì)算
已經(jīng)有24人回復(fù)
- 【求助】castep中的參數(shù)對計(jì)算結(jié)果有什么影響
已經(jīng)有11人回復(fù)
huxiao.gl
木蟲 (小有名氣)
- 應(yīng)助: 37 (小學(xué)生)
- 金幣: 2939.4
- 紅花: 4
- 帖子: 191
- 在線: 373.4小時(shí)
- 蟲號: 1849937
- 注冊: 2012-06-06
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
yeahgyz
2樓2014-07-19 02:48:21
yeahgyz
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 302.5
- 帖子: 48
- 在線: 37.9小時(shí)
- 蟲號: 1478563
- 注冊: 2011-11-06
- 專業(yè): 礦物工程與物質(zhì)分離科學(xué)
3樓2014-07-19 09:01:42
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3509.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲號: 1369423
- 注冊: 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
|
測試過程建議采用結(jié)構(gòu)單點(diǎn)，一旦確定參數(shù)，做結(jié)構(gòu)優(yōu)化，并在優(yōu)化結(jié)構(gòu)的單點(diǎn)做參數(shù)的檢驗(yàn)，這樣可以節(jié)約時(shí)間。最后確定截?cái)嗄芎蚹不是隨便設(shè)定一個(gè)，而是由測試確定，如果很小的截?cái)嗄苣軌驖M足測試條件，直接使用也無妨（當(dāng)然我覺得不太可能）。如果測試不能達(dá)到標(biāo)準(zhǔn)，或是標(biāo)準(zhǔn)較低，導(dǎo)致的結(jié)果就是計(jì)算“力”及“應(yīng)力”會(huì)產(chǎn)生偏差，導(dǎo)致結(jié)構(gòu)優(yōu)化失敗。除了這兩個(gè)參數(shù)，空帶數(shù)目是否充分也很重要。 原則上來說，截?cái)嗄芎蚹越充分結(jié)果是越好的，但是由于DFT的特性，于是不見得是這樣。對于某些結(jié)構(gòu)，PBEsol也許比PBE好得多。 |
4樓2014-07-19 09:55:00
huxiao.gl
木蟲 (小有名氣)
- 應(yīng)助: 37 (小學(xué)生)
- 金幣: 2939.4
- 紅花: 4
- 帖子: 191
- 在線: 373.4小時(shí)
- 蟲號: 1849937
- 注冊: 2012-06-06
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
5樓2014-07-19 16:42:34
yeahgyz
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 302.5
- 帖子: 48
- 在線: 37.9小時(shí)
- 蟲號: 1478563
- 注冊: 2011-11-06
- 專業(yè): 礦物工程與物質(zhì)分離科學(xué)
6樓2014-07-19 16:55:59
yeahgyz
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 302.5
- 帖子: 48
- 在線: 37.9小時(shí)
- 蟲號: 1478563
- 注冊: 2011-11-06
- 專業(yè): 礦物工程與物質(zhì)分離科學(xué)
7樓2014-07-19 17:01:39
huxiao.gl
木蟲 (小有名氣)
- 應(yīng)助: 37 (小學(xué)生)
- 金幣: 2939.4
- 紅花: 4
- 帖子: 191
- 在線: 373.4小時(shí)
- 蟲號: 1849937
- 注冊: 2012-06-06
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
8樓2014-07-19 20:53:54
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3509.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲號: 1369423
- 注冊: 2011-08-16
- 性別: GG
- 專業(yè): 金屬材料的磨損與磨蝕
- 管轄: 第一性原理
|
只計(jì)算能量而不解析導(dǎo)數(shù)優(yōu)化結(jié)構(gòu)，即單點(diǎn)能（single point）。 能不能吻合實(shí)驗(yàn)在于XC functional和pseudopotential的選取，怎么選取不妨參考不少測試的文獻(xiàn)。據(jù)我（經(jīng)驗(yàn)）所知，一些離子型的結(jié)構(gòu)PBEsol稍微好些，不過均勻結(jié)構(gòu)（如純金屬）往往LDA就能給出好的結(jié)果，這些往往在一些泛函性能評估的文獻(xiàn)能夠查到。重點(diǎn)看后續(xù)工作要做什么，目前沒有一種泛函能夠保證任何性質(zhì)都和實(shí)驗(yàn)吻合。另外我了解到的情況是含Cu的晶格常數(shù)似乎經(jīng)常不準(zhǔn)確，也許和關(guān)聯(lián)效應(yīng)有關(guān)（這點(diǎn)目前我不確定）。 測試過程并不是一味增加。事實(shí)上，樓上的結(jié)論正確，2%的偏差不算大。如我前面所說，如果泛函的性能不合適，提升精度反倒是會(huì)產(chǎn)生偏離，粗略的精度恰好更接近實(shí)驗(yàn)也只是巧合，這種準(zhǔn)確對理論計(jì)算其實(shí)沒啥意義。 |
9樓2014-07-19 23:05:18
10樓2016-06-01 22:06:46
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 yeahgyz 的主題更新
| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085700資源與環(huán)境308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 085410人工智能專碩317求調(diào)劑（0854都可以） +3 | xbxudjdn 2026-03-18 | 3/150 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +7 | 妮妮ninicgb 2026-03-15 | 11/550 |
|
|
|
[考研] 【同濟(jì)軟件】軟件（085405）考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢人 2026-03-18 | 4/200 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +6 | 困于星晨 2026-03-17 | 6/300 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 307求調(diào)劑 +3 | 冷笙123 2026-03-17 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 318求調(diào)劑 +3 | Yanyali 2026-03-15 | 3/150 |
|
|
|
[考研]
 機(jī)械專碩調(diào)劑
+3
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
|
[考研] 0856專碩279求調(diào)劑 +5 | 加油加油！? 2026-03-15 | 5/250 |
|
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
|
[考研] 307求調(diào)劑 +5 | 超級伊昂大王 2026-03-12 | 5/250 |
|
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
|
[考研] 085600材料與化工 309分請求調(diào)劑 +7 | dtdxzxx 2026-03-12 | 8/400 |
|
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|
|
|
[考研] 321求調(diào)劑（食品/專碩） +3 | xc321 2026-03-12 | 6/300 |
|