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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 四方ZrO2性质计算遇到的问题
101/1返回列表
查看: 1482  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-29 11:39 by jCd0dEvKHShX
[基金申请] 我就是申请一个面上项目而已，这评审意见是按照杰青的条件评的吧？ +6 gouxfjh 2026-08-28 9/450 2026-08-29 09:59 by jklily
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 05:45 by gy1nBQXYQJqL
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +4 gdd2018 2026-08-28 9/450 2026-08-28 14:22 by jklily
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉…… +5 a089 2026-08-26 6/300 2026-08-27 14:05 by jonewore
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 国合里面能看到了 +7 一怀馨秋 2026-08-26 7/350 2026-08-26 11:23 by zhaosm1982
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
信息提示
关闭
请填处理意见
关闭