zjxitcc 发表于 2020-8-17 11:13 谢谢，我们当时是用Turbomole算的，先随便收敛到一个broken symmetry解，然后做轨道局域化，改变局域轨道的占据数，再做SCF迭代，效果也不错。 |
本帖最后由 zjxitcc 于 2020-8-17 11:15 编辑 wzkchem5 发表于 2020-8-17 11:06 全局搜索确实是几乎不可能的。不过团簇体系可以借助“片段组合波函数” http://sobereva.com/82 http://gaussian.com/afc/ 构建初猜，可以十分有效的找出 大多数低能的稳定波函数。我做过不少单核和双核过渡金属，结果靠谱  用这个一下子找出三四个稳定的SCF解出来不是难事 |
zjxitcc 发表于 2020-8-17 01:07 而且波函数即使稳定性检测过了，也只是局部稳定，不一定是全局最优的波函数。找到全局最优的波函数，哪怕对于HF而言，也是NP complete的。我师弟就遇到过一个铁酶的体系，存在局部稳定但是不是全局最优的解，把他坑惨了，是在扫势能曲线的时候波函数自己从一个解跳到另一个解才发现的，最后整个重算 |
wzkchem5 发表于 2020-8-16 14:44 反铁磁就算是UKS配合自旋极化的初猜，结果也未必稳定 说明当前缺乏有效的工具 |
zjxitcc 发表于 2020-8-16 12:44 像这种铁磁性体系还好（自旋平行的条件下遇到波函数不稳定的概率小），要是反铁磁耦合的体系，就很容易被坑了，比如有triplet instability的情况下，仍然从头到尾都用RKS算。 |
本帖最后由 zjxitcc 于 2020-8-16 12:45 编辑 wangxubo 发表于 2020-8-16 08:48 一方面是因为12#说的原因；另一方面，即使加上周期性使波函数变简单了，复杂体系的波函数稳定性仍然是要检测的，但是由于在周期性计算这个领域没有现成的工具（检验波函数稳定性），大家都“掩耳盗铃”地不检测，也就罢了，没有人会去刁难。反过来在团簇领域就不一样了，能检验波函数稳定性的软件还是又不少的。 这也反映出，周期性计算这个领域急需检验波函数稳定性的有效工具。 |
wangxubo 发表于 2020-8-16 08:48 因为用了周期性边界条件的情况下，少了很多边界原子，所以波函数变简单了 |
zjxitcc 发表于 2020-7-20 17:12 你这个说法没道理啊，既然gto方法即使收敛也会收敛到local minima的话，难道用平面波就会收敛到更好的波函数上吗？还是说你想表达其他的意思？ |
beefly 发表于 2020-7-21 17:13 谢谢老师指点 |
beefly 发表于 2020-7-21 17:13 这体系75铁原子，应该是300单电子 |
| 这种体系通常每个Fe原子都是高自旋，整个体系的总自旋自己去算。默认的低自旋组态对应能量超高的激发态，甚至根本没有物理意义，即使收敛了也没用 |
sobereva 发表于 2020-7-21 00:15 谢谢老师，这就学学QE |
zjxitcc 发表于 2020-7-20 17:12 谢谢老师 |
| d、f族金属团簇是最难搞的一类体系。如果你是想用团簇模型代替周期性体系来研究金属表面问题，就别折腾了。QE又免费上手又容易。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-22 20:54 , Processed in 0.171825 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites