计算化学公社

标题: ZrNi体系-差分电荷密度出错求助 [打印本页]

作者
Author: 风飞 时间: 2021-3-20 07:26
标题: ZrNi体系-差分电荷密度出错求助
老师，您好，我在计算Ni/ZrO2(111)体系（ZrO2的111面吸附1个Ni原子）的差分电荷密度时，出现错误：
片段scf计算（Ni-scf,Zr-scf,ZrNi-scf）的输入文件如图所示
(, 下载次数 Times of downloads: 46)
片段的电荷密度计算和Δp计算的输入文件，出现如下错误提示：屏幕提示错误，out文件中错误提示，
以及使用VEST进行差分电荷密度计算的错误提示
(, 下载次数 Times of downloads: 41)
我能保证片段的坐标是一一对应的
这些命令，我算其他体系也没有问题，请问可能是什么原因呢？

作者
Author: 卡开发发 时间: 2021-3-20 09:51
本帖最后由卡开发发于 2021-3-20 09:53 编辑

三个晶胞参数都一样？可以把三个cube文件格点数据之前的字段发来看看就知道问题在哪了。

作者
Author: 风飞 时间: 2021-3-20 11:31
本帖最后由风飞于 2021-3-20 11:33 编辑

卡开发发发表于 2021-3-20 09:51
三个晶胞参数都一样？可以把三个cube文件格点数据之前的字段发来看看就知道问题在哪了。

老师，您好，我确认了Ni-scf,Zr-scf,ZrNi-scf 三个文件中的晶胞参数是一样的，cube文件中的数据格点如下，的确不一样，请问是什么原因呢？ (, 下载次数 Times of downloads: 38)

作者
Author: 卡开发发 时间: 2021-3-20 11:48

风飞发表于 2021-3-20 11:31
老师，您好，我确认了Ni-scf,Zr-scf,ZrNi-scf 三个文件中的晶胞参数是一样的，cube文件中的数据格点如下 ...

Ni的那个有问题，把Ni涉及到的in文件都发来看看。

作者
Author: 风飞 时间: 2021-3-20 12:03
本帖最后由风飞于 2021-3-20 12:07 编辑

卡开发发发表于 2021-3-20 11:48
Ni的那个有问题，把Ni涉及到的in文件都发来看看。

老师，您好，关于Ni的scf计算和电子密度计算文件输入文件如下 (, 下载次数 Times of downloads: 1)
(, 下载次数 Times of downloads: 1)
其对应的相关图片如下
(, 下载次数 Times of downloads: 62)

作者
Author: 卡开发发 时间: 2021-3-20 15:32

风飞发表于 2021-3-20 12:03
老师，您好，关于Ni的scf计算和电子密度计算文件输入文件如下

其对应的相关图片如下

原子自洽的时候用nosym=.true.

作者
Author: 风飞 时间: 2021-3-20 16:34

卡开发发发表于 2021-3-20 15:32
原子自洽的时候用nosym=.true.

嗯，老师，您好，我加了nosym=.true. 就好了，
我看了我的其他体系，对于Ni原子的scf计算没有加nosym 也成功了，请问是什么原因呢？

作者
Author: 卡开发发 时间: 2021-3-20 16:40

风飞发表于 2021-3-20 16:34
嗯，老师，您好，我加了nosym=.true. 就好了，
我看了我的其他体系，对于Ni原子的scf计算没有加nosym ...

可能格点比较巧合，正好其他的体系判断不出这种不合理的“对称性”，目前我也很难判断原因。

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