作者peter308 (pete)
看板Chemistry
标题[请益]请问板上有人做DFT的分子结构最佳化的吗 …
时间Thu Jul 24 22:12:14 2008
※ [本文转录自 Physics 看板]
作者: peter308 (pete) 看板: Physics
标题: [请益]请问板上有人做DFT的分子结构最佳化的吗????
时间: Thu Jul 24 22:00:45 2008
有人用软体算材料结构最佳化的吗????
我最近用一套软体算纳米金属颗粒(Ag13)的结构优化
发现结构一直无法收敛....
程式会用gradient force的大小判断是否收敛
优化过程gradient force正常的情形下应该是会渐渐变小....
但我会一直在某个值震荡
导致我一直算不出结果
板上有人做这方面相关的吗...???
或是能提供其他我可能找的到答案的讨论区或资源
被这个问题卡了很久.....
感谢分享!!!
--
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 118.161.155.146
--
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 118.161.155.146
1F：推 dkw:你取的起始结构不够稳定...可能会出现这样的现象 07/24 23:33
2F：→ peter308:我取的初始结构是经过优化过的..但用的是empirical位势 07/24 23:43
3F：→ peter308:同样的结构在某些basis set是可以收敛的..这个 07/24 23:44
4F：→ peter308:不收敛的basis为model-core potential basis set 07/24 23:45
※ 编辑: peter308 来自: 118.161.155.146 (07/24 23:47)
5F：推 kkarthur:如果是Gaussian的话 07/25 14:15
6F：→ kkarthur:opt=calcfc scf=xqc 07/25 14:16
7F：推 kancu:真的要先问一下你用的是什麽软体?还有参数你怎麽设? 07/27 07:45
8F：→ kancu:另外~DFT那麽多~你用哪一个? 07/27 07:46
9F：→ kancu:如果真的是用Gaussian的话~如楼上所说~加scf=xqc吧~ 07/27 07:47
10F：→ kancu:xqc不行的~~就改成scf=qc吧~ 07/27 07:54
11F：推 theory:他用的软体是 deMon2k，我个人是没有听过啦.. 07/27 17:15
12F：→ peter308:我用的是deMon2k....真的有点冷门...但还蛮好入手的 08/02 17:35
13F：→ peter308:另外就是guassian的basis library有model core basis 08/02 17:41
14F：→ peter308:可以选吗...???? 我想用guassian来试试....有人知道吗?? 08/02 17:41
15F：→ kancu:model core basis..这个要查一下才知道~~目前是没看过... 08/08 01:16
16F：→ kancu:如果你能知道那个basis set的内容~或许可以试着改成gaussian 08/08 01:18
17F：→ kancu:的格式..那就可以在gaussian用了~~ 08/08 01:19
18F：→ kancu:只是我看你是作材料的计算~用gaussian或许会有点难度吧.. 08/08 01:20
19F：→ kancu:材料的计算，或许可以试试看VASP，满多人用这套在算材料的~ 08/08 01:20
20F：推 kancu:查了一下~目前gaussian是没有MCP的basis set..可是deMon2k有 08/08 13:22
21F：→ kancu:所以你只能自己建gbs~或是改参数试试~~ 08/08 13:24