小木虫

当前位置：首页 > 第一原理 >四方ZrO2性质计算遇到的问题

四方ZrO2性质计算遇到的问题

作者 jtmcs: 来源: 小木虫 450 9 举报帖子

各位虫友，有人计算过四方相ZrO2 的性质吗？我用ZrO2 的惯用元胞做的结构优化，之后用优化得出的参数建立的初级元胞，结果用初级元胞做一个测试性的静态计算时发现，元胞中的四个氧原子受到了很大力。我之前的优化结果是没有问题的，构造元胞时也没有问题，求教为何会出现这种情况？返回小木虫查看更多

今日热帖

精华评论

jtmcs

这是我用惯用元胞优化后，OUTCAR中的受力情况，应该没什么问题吧
POSITION                                     TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
   0.00000    0.00000    0.00000       0.000000    0.000000    0.000000
   2.51770    2.51770    0.00000       0.000000    0.000000    0.000000
   0.00000    2.51770    2.55673       0.000000    0.000000    0.000000
   2.51770    0.00000    2.55673       0.000000    0.000000    0.000000
   1.25885    1.25885    1.06437       0.000000    0.000000    -0.002726
   1.25885    1.25885    3.62110       0.000000    0.000000    -0.002726
   3.77655    3.77655    1.06437       0.000000    0.000000    -0.002726
   3.77655    3.77655    3.62110       0.000000    0.000000    -0.002726
   1.25885    3.77655    1.49236       0.000000    0.000000    0.002726
   1.25885    3.77655    4.04909       0.000000    0.000000    0.002726
   3.77655    1.25885    1.49236       0.000000    0.000000    0.002726
   3.77655    1.25885    4.04909       0.000000    0.000000    0.002726
-----------------------------------------------------------------------------------
total drift:                               0.000000    0.000000    0.000000

--------------------------------------------------------------------------------------------------------，
jtmcs

优化得到的CONTCAR如下
ZrO2-tetragonal
1.00000000000000
   5.0353934611750928 0.0000000000000000 0.0000000000000000
   0.0000000000000000 5.0353934611750928 0.0000000000000000
   0.0000000000000000 0.0000000000000000 5.1134555742906809
Zr O
   4    8
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000
  0.5000000000000000  0.5000000000000000  0.0000000000000000
  0.0000000000000000  0.5000000000000000  0.5000000000000000
  0.5000000000000000  0.0000000000000000  0.5000000000000000
  0.2500000000000000  0.2500000000000000  0.2081508361397297
  0.2500000000000000  0.2500000000000000  0.7081508361397297
  0.7500000000000000  0.7500000000000000  0.2081508361397297
  0.7500000000000000  0.7500000000000000  0.7081508361397297
  0.2500000000000000  0.7500000000000000  0.2918491638602703
  0.2500000000000000  0.7500000000000000  0.7918491638602703
  0.7500000000000000  0.2500000000000000  0.2918491638602703
  0.7500000000000000  0.2500000000000000  0.7918491638602703

  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
  0.00000000E+00  0.00000000E+00  0.00000000E+00
jtmcs

这是根据优化的结果，建立的初级元胞的POSCAR
ZrO2-tetragonal-primitive cell
1.0
      3.5605608623       0.0000000000       0.0000000000
      0.0000000000       3.5605608623       0.0000000000
      0.0000000000       0.0000000000       5.1134555743
Zr O
2 4
Direct
   0.000000000       0.000000000       0.500000000
   0.500000000       0.500000000       0.000000000
   0.500000000       0.000000000       0.291849164
   0.500000000       0.000000000       0.791849164
   0.000000000       0.500000000       0.208150836
   0.000000000       0.500000000       0.708150836
hs0218

如果要重复这些个计算使之达到大小模型等价
需要考虑k点密度
模型减半，k点是需要增加的

包括你用了NPAR和KAPR等并行参数都会影响结果
具体的你可以grep k-point OUTCAR grep NABND OUTCAR，看看俩模型对比下
jtmcs

引用回帖:
5楼: Originally posted by hs0218 at 2017-12-01 00:43:22
如果要重复这些个计算使之达到大小模型等价
需要考虑k点密度
模型减半，k点是需要增加的

包括你用了NPAR和KAPR等并行参数都会影响结果
具体的你可以grep k-point OUTCAR grep NABND OUTCAR，看看俩模型对比下

多谢。K点的数量我改过了。后来设置了一下NPAR, 计算结果就奇迹般的对了。但讲道理NPAR影响的应该是计算的速率呀，为啥会对计算的结果产生影响呢
jtmcs

后来经过尝试，发现是NPAR的问题，我就感觉很奇怪，NPAR影响的应该是计算的速率，为啥会对结果产生影响？后来我又进行了进一步的测试，我调用了两个节点，共24个核，NAPR蚕蛹默认值核24 的时候都会出现问题，如果是其他的数值，例如2、4、6、12都没有问题。请教诸位大神这是什么情况
我后来怀疑是节点之间的数据交换出现了问题，座椅又采用了一个节点来计算（12个核），发现还是同样的问题：NPAR=12或默认值都有问题，其他的值如2、4、6等就没有问题，求助各位大侠，这是为什么？（运行环境是天津超算）
hs0218

引用回帖:
7楼: Originally posted by jtmcs at 2017-12-01 12:32:51
后来经过尝试，发现是NPAR的问题，我就感觉很奇怪，NPAR影响的应该是计算的速率，为啥会对结果产生影响？后来我又进行了进一步的测试，我调用了两个节点，共24个核，NAPR蚕蛹默认值核24 的时候都会出现问题，如果是 ...

NBAND不单单是效率，而且是band的分配，这个影响数据的结果，这个早就有人测试过。如果NBAND不能被NPAR整除，就会有问题