 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 2550 | 回復(fù): 9 | |||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||
yeahgyz銅蟲 (初入文壇)
|
[求助]
黃銅礦幾何結(jié)構(gòu)優(yōu)化 按照文獻(xiàn)參數(shù)不能收斂，且晶胞常數(shù)誤差過大如何解決 已有1人參與
|
||||
|
小弟按照文獻(xiàn)所述參數(shù)對(duì)黃銅礦進(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)該是不能用的。 請(qǐ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ù)此樓» 猜你喜歡
求b區(qū)院校調(diào)劑
已經(jīng)有4人回復(fù)
-
材料與化工考研調(diào)劑
已經(jīng)有9人回復(fù)
-
求調(diào)劑
已經(jīng)有3人回復(fù)
-
347求調(diào)劑
已經(jīng)有3人回復(fù)
-
總分293求調(diào)劑
已經(jīng)有4人回復(fù)
-
284求調(diào)劑
已經(jīng)有15人回復(fù)
-
359求調(diào)劑
已經(jīng)有3人回復(fù)
-
各位老師您好：本人初試372分
已經(jīng)有5人回復(fù)
-
285求調(diào)劑
已經(jīng)有3人回復(fù)
-
282求調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- CASTEP幾何結(jié)構(gòu)優(yōu)化的問題 雖然收斂了 感覺不是很對(duì) 請(qǐng)大家?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ù)對(duì)計(jì)算結(jié)果有什么影響
已經(jīng)有11人回復(fù)
卡開發(fā)發(fā)
專家顧問 (著名寫手)
Ab Initio Amateur
-
專家經(jīng)驗(yàn): +224 - 1ST強(qiáng)帖: 2
- 應(yīng)助: 47 (小學(xué)生)
- 金幣: 3513.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲號(hào): 1369423
- 注冊(cè): 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é)果，這些往往在一些泛函性能評(píng)估的文獻(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)確對(duì)理論計(jì)算其實(shí)沒啥意義。 |
9樓2014-07-19 23:05:18
huxiao.gl
木蟲 (小有名氣)
- 應(yīng)助: 37 (小學(xué)生)
- 金幣: 2939.4
- 紅花: 4
- 帖子: 191
- 在線: 373.4小時(shí)
- 蟲號(hào): 1849937
- 注冊(cè): 2012-06-06
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
yeahgyz2樓2014-07-19 02:48:21
yeahgyz
銅蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 302.5
- 帖子: 48
- 在線: 37.9小時(shí)
- 蟲號(hào): 1478563
- 注冊(cè): 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é)生)
- 金幣: 3513.8
- 散金: 14
- 紅花: 223
- 帖子: 2477
- 在線: 1254小時(shí)
- 蟲號(hào): 1369423
- 注冊(cè): 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的特性，于是不見得是這樣。對(duì)于某些結(jié)構(gòu)，PBEsol也許比PBE好得多。 |
4樓2014-07-19 09:55:00
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +10 | 晨昏線與星海 2026-03-19 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 0703化學(xué)求調(diào)劑 +6 | 奶油草莓. 2026-03-22 | 7/350 |
|
|
|
[考研] 289材料與化工（085600）B區(qū)求調(diào)劑 +4 | 這么名字咋樣 2026-03-22 | 5/250 |
|
|
|
[考研] 0854AI CV方向招收調(diào)劑 +3 | 章小魚567 2026-03-23 | 3/150 |
|
|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-24 | 3/150 |
|
|
|
[考研] 材料專碩找調(diào)劑 +5 | 哈哈哈吼吼吼哈 2026-03-23 | 5/250 |
|
|
|
[考研] 291求調(diào)劑 +3 | HanBeiNingZC 2026-03-24 | 3/150 |
|
|
|
[考研] 材料292調(diào)劑 +8 | 橘頌思美人 2026-03-23 | 8/400 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
|
[考研] 335求調(diào)劑 +4 | yuyu宇 2026-03-23 | 5/250 |
|
|
|
[考研] 一志愿上海交大生物與醫(yī)藥專碩324分，求調(diào)劑 +5 | jiajunX 2026-03-22 | 5/250 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
|
[考研] 298求調(diào)劑 +8 | 上岸6666@ 2026-03-20 | 8/400 |
|
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 313求調(diào)劑 +4 | 肆叁貳壹22 2026-03-19 | 4/200 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
|
[考研] 353求調(diào)劑 +3 | 拉鉤不許變 2026-03-20 | 3/150 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|