计算化学公社

标题: 求助：cp2k优化单镍晶胞输出文件4个yes，1个No。算优化成功了么？ [打印本页]

作者
Author: ldwf 时间: 2024-10-23 14:13
标题: 求助：cp2k优化单镍晶胞输出文件4个yes，1个No。算优化成功了么？
本帖最后由 ldwf 于 2024-10-23 16:13 编辑

想用Ni111晶面进行氨基和羧基的吸附能计算。
单镍晶胞优化输出文件显示4个yes。
不知道输入文件是否有错误。
请前辈们看看。

作者
Author: kirSTY 时间: 2024-10-23 14:55
1.首先可以确定目前的计算是没有意义的，不要开“IGNORE_CONVERGENCE_FAILURE”,参见“http://sobereva.com/61”的帖子最后“附：对Gaussian中IOp(5/13=1)的大批判”
2.过渡元素计算确实收敛难度高，多看论文
3.其余参见“http://sobereva.com/665”

作者
Author: ldwf 时间: 2024-10-23 14:58

kirSTY 发表于 2024-10-23 14:55
1.首先可以确定目前的计算是没有意义的，不要开“IGNORE_CONVERGENCE_FAILURE”,参见“http://sobereva.com ...

cp2k 2024.1版本的，没加这个，第一步就停了。然后报错提示，让加这个语句。

作者
Author: ldwf 时间: 2024-10-23 15:03

kirSTY 发表于 2024-10-23 14:55
1.首先可以确定目前的计算是没有意义的，不要开“IGNORE_CONVERGENCE_FAILURE”,参见“http://sobereva.com ...

单胞和切面，在优化结构输入文件主要参数上的差别会很大么？

作者
Author: kirSTY 时间: 2024-10-23 15:05
报错的意思应当是“可以通过加入XXX关键词避免提示”，而不是“建议添加XXX关键词”，这二者之间有区别的。
“http://sobereva.com/61”中说了，忽略SCF不收敛必然导致受力计算烂，那么在这个基础上优化出来的所有结果皆是不可信的，所以目前应当解决SCF不收敛，而不是忽略SCF不收敛。

作者
Author: qczcb2 时间: 2024-10-23 15:48
这是切好了表面吗？如果是表面，z方向不需要k点；xy方向k点数目不够，最好改成9*9*1以上；自旋多重度不对，如果镍的磁矩为1的话，铁磁的多重度应为5；反铁磁需要手动设置两个镍物种分别带1和-1的磁矩，设置多重度为1并开UKS；另外，表面切的太薄了，建议再多几层。

作者
Author: ldwf 时间: 2024-10-23 15:59

qczcb2 发表于 2024-10-23 15:48
这是切好了表面吗？如果是表面，z方向不需要k点；xy方向k点数目不够，最好改成9*9*1以上；自旋多重度不对， ...

这个没切面，切面的输入文件在另一个帖子里，优化不成功。
然后这个是优化一下单胞的。这样看来错误很多。

作者
Author: ldwf 时间: 2024-10-23 16:12

qczcb2 发表于 2024-10-23 15:48
这是切好了表面吗？如果是表面，z方向不需要k点；xy方向k点数目不够，最好改成9*9*1以上；自旋多重度不对， ...

我重新编辑了一下帖子，把111晶面优化的输入文件放上了，请前辈空闲时间帮忙看看。

作者
Author: qczcb2 时间: 2024-10-23 21:16
如果你是想算反铁磁的话，
&KIND Ni
   ELEMENT Ni
   BASIS_SET DZVP-MOLOPT-SR-GTH-q18
   POTENTIAL GTH-PBE
加上  MAGNETIZATION  -1.00
&END KIND
&KIND Ni_1
   ELEMENT Ni
   BASIS_SET DZVP-MOLOPT-SR-GTH-q18
   POTENTIAL GTH-PBE
   MAGNETIZATION  1.00
&END KIND
  &END SUBSYS
多重度要改为1；铁磁的话就按multiwfn把所有镍都选上设1的磁矩自动设多重度就好了。现在多重度7莫名其妙，就算按你现在的设法只有一半原子带磁性也应该是1+8=9.另外，应该先把晶胞优化好，你的晶胞明显不对，z方向不可能这么大，如果是用这个晶胞结构切的晶面那不可信

作者
Author: ldwf 时间: 2024-10-23 21:36

qczcb2 发表于 2024-10-23 21:16
如果你是想算反铁磁的话，
&KIND Ni
ELEMENT Ni

晶面结构不是从帖子里面的Ni单胞切出来的，是我直接用MP里面下载的单胞切的。然后优化不成，我就打算试试单胞使用同样的条件。

冻结的原子也要带磁性么，所有原子都带磁性的话，生成的输入文件里面没有Ni_1关键词字样。

作者
Author: qczcb2 时间: 2024-10-23 22:37

ldwf 发表于 2024-10-23 21:36
晶面结构不是从帖子里面的Ni单胞切出来的，是我直接用MP里面下载的单胞切的。然后优化不成，我就打算试试 ...

冻结的原子只是不被移动，但是scf计算也被包括在内，如果有磁性不加的话肯定算不准啊，所以需要加上。如果某种元素全是同一个磁矩，multiwfn不会额外设置kind。额外设置是方便你做反铁磁等情况用的。

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