计算化学公社

标题: 结构优化的一些小问题 [打印本页]

作者
Author: 萧瑟 时间: 2021-11-22 09:59
标题: 结构优化的一些小问题
从入门到现在也只做到结构优化、晶胞优化这样的，发现了一点小问题，导致我有了一个猜测，就是我在优化一个PBC结构的时候，总感觉是一个小晶胞优化后通过膨胀减小原子之间的斥力要容易一些，但是反过来，一个大晶胞需要通过收缩来增大原子间的引力就很困难。就之前优化一个合金结构，我在Pd晶胞里掺杂Fe 就一直scf不收敛，但是我在Fe晶胞里掺杂Pd，不仅SCF收敛了而且最终也优化完成了。所以不知道我的感觉猜测是不是对的，希望有大佬能指正一下。

作者
Author: sobereva 时间: 2021-11-23 08:23
这是错觉
从优化算法上来说，对同一个体系，从相对于平衡结构更小的结构和更大结构开始优化没有差异（差异仅在收敛限以内）
你的情况两种掺杂方式得到的结构如果存在区别，没有可比性

作者
Author: 萧瑟 时间: 2021-11-23 16:18

sobereva 发表于 2021-11-23 08:23
这是错觉
从优化算法上来说，对同一个体系，从相对于平衡结构更小的结构和更大结构开始优化没有差异（差异 ...

请问sob老师有没有什么比较好的掺杂方法可以避免像我第一种情况那样SCF都不收敛的掺杂结构出现呢？

作者
Author: wzkchem5 时间: 2021-11-23 16:40
其实这只是SCF收敛的问题，和结构优化收敛没关系，因为SCF如果不收敛，即便结构优化理论上再容易，也没法给你算。而Fe和Pd显然不止是尺寸的区别，电子结构也有很大区别。要得出你的结论，起码应该试一下比如说Cu里掺Ag和Ag里掺Cu

作者
Author: sobereva 时间: 2021-11-24 07:31

萧瑟发表于 2021-11-23 16:18
请问sob老师有没有什么比较好的掺杂方法可以避免像我第一种情况那样SCF都不收敛的掺杂结构出现呢？

遇到SCF不收敛，倘若结构本身没有硬伤，首先应该想到的是尝试各种关键词想办法让SCF收敛。仅当实在没辙再考虑换其它结构
怎么掺杂和是否能最终想办法让SCF收敛没有直接联系

作者
Author: 萧瑟 时间: 2021-11-24 09:38

wzkchem5 发表于 2021-11-23 16:40
其实这只是SCF收敛的问题，和结构优化收敛没关系，因为SCF如果不收敛，即便结构优化理论上再容易，也没法给 ...

谢谢回复，有时间计算一下看看这两种情况是不是一样的

作者
Author: 萧瑟 时间: 2021-11-24 09:49

sobereva 发表于 2021-11-24 07:31
遇到SCF不收敛，倘若结构本身没有硬伤，首先应该想到的是尝试各种关键词想办法让SCF收敛。仅当实在没辙再考 ...

结构是从Materials Project上下载的应该是没有问题的（晶胞单独跑一遍的话十几个电子步就收敛了），掺杂的话就是inp文件里改换几个坐标前的元素。有这样的错觉主要是同样的inp文件其他参数一样，我把原子坐标里的Pd改了几个成Fe，然后反过来第二个是把Fe改了几个成Pd ，两种情况原子比例是一样的铁钯七三开的样子，结果就变成了一个收敛一个不收敛了。然后我从MS里看到Pd的晶胞原子间距比Fe的要大一些，所以才有了我之前的一些想法。

作者
Author: wzkchem5 时间: 2021-11-24 15:54

萧瑟发表于 2021-11-24 02:49
结构是从Materials Project上下载的应该是没有问题的（晶胞单独跑一遍的话十几个电子步就收敛了），掺 ...

还是得控制变量才能得到这个结论，Pd和Fe的d电子数也差很多，光由这个计算说明不了尺寸是主要影响因素

作者
Author: 萧瑟 时间: 2021-11-24 16:27

wzkchem5 发表于 2021-11-24 15:54
还是得控制变量才能得到这个结论，Pd和Fe的d电子数也差很多，光由这个计算说明不了尺寸是主要影响因素

所以说是猜测的嘛。毕竟自己学识有限最直观的感觉就是这样的。
当然不是我所想的这样是最好的了，毕竟就没那么玄学了

欢迎光临计算化学公社 (http://bbs.keinsci.com/)