ZrNi体系-差分电荷密度出错求助

风飞 · 发表于 Post on 2021-3-20 07:26:55

老师，您好，我在计算Ni/ZrO2(111)体系（ZrO2的111面吸附1个Ni原子）的差分电荷密度时，出现错误：
片段scf计算（Ni-scf,Zr-scf,ZrNi-scf）的输入文件如图所示

片段的电荷密度计算和Δp计算的输入文件，出现如下错误提示：屏幕提示错误，out文件中错误提示，
以及使用VEST进行差分电荷密度计算的错误提示

我能保证片段的坐标是一一对应的
这些命令，我算其他体系也没有问题，请问可能是什么原因呢？

卡开发发 · 发表于 Post on 2021-3-20 09:51:41

本帖最后由卡开发发于 2021-3-20 09:53 编辑

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

风飞 · 发表于 Post on 2021-3-20 11:31:00

本帖最后由风飞于 2021-3-20 11:33 编辑

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

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

卡开发发 · 发表于 Post on 2021-3-20 11:48:58

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

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

风飞 · 发表于 Post on 2021-3-20 12:03:55

本帖最后由风飞于 2021-3-20 12:07 编辑

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

老师，您好，关于Ni的scf计算和电子密度计算文件输入文件如下

ZrNi3-Ni-density.in (186 Bytes, 下载次数 Times of downloads: 1)

ZrNi3-Ni-scf.in (824 Bytes, 下载次数 Times of downloads: 1)
其对应的相关图片如下

卡开发发 · 发表于 Post on 2021-3-20 15:32:49

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

其对应的相关图片如下

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

风飞 · 发表于 Post on 2021-3-20 16:34:11

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

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

卡开发发 · 发表于 Post on 2021-3-20 16:40:45

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

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

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Quantum ESPRESSO] ZrNi体系-差分电荷密度出错求助