应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 四方ZrO2性质计算遇到的问题
101/1返回列表
查看: 1413  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

jtmcs

铁虫 (初入文坛)

[求助] 四方ZrO2性质计算遇到的问题 已有1人参与

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

» 猜你喜欢
1楼 2017-11-30 19:31:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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


--------------------------------------------------------------------------------------------------------
赞一下 回复此楼
2楼2017-11-30 19:37:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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
赞一下 回复此楼
3楼2017-11-30 19:38:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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
赞一下 回复此楼
4楼2017-11-30 19:40:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hs0218

金虫 (正式写手)

DDr.

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
jtmcs: 金币+3, ★有帮助 2017-12-01 13:01:28
如果要重复这些个计算使之达到大小模型等价
需要考虑k点密度
模型减半，k点是需要增加的

包括你用了NPAR和KAPR等并行参数都会影响结果
具体的你可以grep k-point OUTCAR grep NABND OUTCAR，看看俩模型对比下
赞一下(1人) 回复此楼
路漫漫其修远兮
5楼2017-12-01 00:43:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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影响的应该是计算的速率呀，为啥会对计算的结果产生影响呢
赞一下 回复此楼
6楼2017-12-01 12:27:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jtmcs

铁虫 (初入文坛)
后来经过尝试，发现是NPAR的问题，我就感觉很奇怪，NPAR影响的应该是计算的速率，为啥会对结果产生影响？后来我又进行了进一步的测试，我调用了两个节点，共24个核，NAPR蚕蛹默认值核24 的时候都会出现问题，如果是其他的数值，例如2、4、6、12都没有问题。请教诸位大神这是什么情况
我后来怀疑是节点之间的数据交换出现了问题，座椅又采用了一个节点来计算（12个核），发现还是同样的问题：NPAR=12或默认值都有问题，其他的值如2、4、6等就没有问题，求助各位大侠，这是为什么？（运行环境是天津超算）
赞一下 回复此楼
7楼2017-12-01 12:32:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hs0218

金虫 (正式写手)

DDr.

【答案】应助回帖

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

NBAND不单单是效率，而且是band的分配，这个影响数据的结果，这个早就有人测试过。如果NBAND不能被NPAR整除，就会有问题
赞一下 回复此楼
路漫漫其修远兮
8楼2017-12-01 15:34:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jtmcs

铁虫 (初入文坛)
引用回帖:
8楼: Originally posted by hs0218 at 2017-12-01 15:34:19
NBAND不单单是效率，而且是band的分配，这个影响数据的结果，这个早就有人测试过。如果NBAND不能被NPAR整除，就会有问题...

这个我是知道的，NPAR就是每次同时计算多少条能带，是要能被NABNDS整除的。但是我最初说的有问题的结果，NBANDS, NPAR, NCORE等采取的都是默认值，按照默认值应该不会出现不能整除的现象吧
赞一下 回复此楼
9楼2017-12-01 19:52:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hs0218

金虫 (正式写手)

DDr.

【答案】应助回帖

引用回帖:
9楼: Originally posted by jtmcs at 2017-12-01 19:52:29
这个我是知道的，NPAR就是每次同时计算多少条能带，是要能被NABNDS整除的。但是我最初说的有问题的结果，NBANDS, NPAR, NCORE等采取的都是默认值，按照默认值应该不会出现不能整除的现象吧...

你测试下NPAR这个参数就知道了....
赞一下 回复此楼
路漫漫其修远兮
10楼2017-12-02 00:22:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 jtmcs 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭