作者peter308 (pete)
看板Physics
標題[請益]請問板上有人做DFT的分子結構最佳化的嗎????
時間Thu Jul 24 22:00:45 2008
有人用軟體算材料結構最佳化的嗎????
我最近用一套軟體算納米金屬顆粒(Ag13)的結構優化
發現結構一直無法收斂....
程式會用gradient force的大小判斷是否收斂
優化過程gradient force正常的情形下應該是會漸漸變小....
但我會一直在某個值震盪
導致我一直算不出結果
板上有人做這方面相關的嗎...???
或是能提供其他我可能找的到答案的討論區或資源
被這個問題卡了很久.....
感謝分享!!!
--
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 118.161.155.146
1F：→ nidoran:你用哪一套??? 07/24 22:34
2F：→ nidoran:收斂的方法好像又有差 07/24 23:46
3F：→ peter308:回樓上大大....我使用的軟體為deMon2k,蠻冷門的...不常見 07/24 23:49
4F：→ peter308:至於你說的收斂方法是一種改良過的牛頓法...叫RFO 07/24 23:50
5F：→ peter308:rational functional optimization.... 07/24 23:51
6F：→ peter308:可以選的收斂選項有三個:RFO,steepest descent 07/24 23:52
7F：→ peter308:potential surface walking...軟體預設是RFO 07/24 23:52
8F：→ pipidog:把mixing scheme改成pulay,history step加多? 07/25 11:16
9F：推 Keelungman:通常減小每一步驟的長度有助穩定收斂 07/25 16:10
10F：→ peter308:其實我有找過資料...不外乎說減少步長,增大cycle數,改變 07/25 17:17
11F：→ peter308:一下初始結構,或是用小一點的基組..等.改變步長我有試過 07/25 17:19
12F：→ peter308:似乎沒什麼大太用處...請問要改到多小呢???預設0.3a.u. 07/25 17:20
13F：→ peter308:我只試過減到0.1a.u.,請問還要再更小嗎????? 07/25 17:21
14F：推 nidoran:振盪的範圍??? 07/26 00:52
15F：推 fauna:運氣不好的話會在不同的electronic state 之間跑來跑去 07/26 02:12
16F：→ fauna:所以收斂不了 07/26 02:12
17F：→ peter308:樓上可以解釋得更詳細一點嗎??????感激!! 07/26 11:38
18F：推 fauna:先改成steepest descent 作幾步之後再換RFO繼續可能比較好 07/29 04:42