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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 四方ZrO2性质计算遇到的问题
101/1返回列表
查看: 1481  |  回复: 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一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:49 by G6APbkg8SA6w
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 09:48 by G6APbkg8SA6w
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:39 by G6APbkg8SA6w
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 ASdOkHsho7FD 2026-08-28 5/250 2026-08-29 09:28 by G6APbkg8SA6w
[考研] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +5 ASdOkHsho7FD 2026-08-28 8/400 2026-08-29 09:17 by G6APbkg8SA6w
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +13 晴天加油 2026-08-26 14/700 2026-08-28 18:36 by huagongfeihu
[基金申请] 基金不中，共勉 +11 eulota 2026-08-26 11/550 2026-08-28 14:22 by 火星超人xi
[基金申请] 基金系统什么内容也没有 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
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[文学芳草园] 梦想 +3 myrtle 2026-08-26 3/150 2026-08-27 10:01 by angelyueyi
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭