作者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