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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » WO3的slab模型问题
51/1返回列表
查看: 3181  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

雪雨星风17

木虫 (小有名气)

无

[求助] WO3的slab模型问题 已有4人参与

首先从material studio的数据库中优化优化得到了WO3的结构，然后切割得到了001面的slab模型，真空厚度10A, encut选了750，poscar如下：
CONTCAR\(0\0\1)                        
   1.00000000000000     
     7.5412001609999999    0.0000000000000000    0.0000000000000000
     0.0000000000000000    7.7389998436000003    0.0000000000000000
     0.0000000000000000    0.0000000000000000   17.8393001556000002
  24   8
Selective dynamics
Direct
  0.2500000000000000  0.2590700089999984  0.0000000000000000   F   F   F
  0.2500000000000000  0.2590200010000032  0.2204499990000031   F   F   F
  0.2500000000000000  0.0120099999999965  0.1107200009999971   F   F   F
  0.2500000000000000  0.0172809866630946  0.3285616818588990   T   T   T
  0.7500000000000000  0.2409300059999993  0.2204400000000035   F   F   F
  0.7500000000000000  0.2409799989999968  0.0000100000000032   F   F   F
  0.7500000000000000  0.4827085057864299  0.3285591810238420   T   T   T
  0.7500000000000000  0.4878999889999989  0.1107399990000033   F   F   F
  0.7500000000000000  0.7157909681018981  0.4471848534696731   T   T   T
  0.7500000000000000  0.7409800289999993  0.2175499949999988   F   F   F
  0.7500000000000000  0.9863484665249046  0.3264384093378995   T   T   T
  0.7500000000000000  0.9879000190000014  0.1068200019999992   F   F   F
  0.2500000000000000  0.7590699790000031  0.2175599930000018   F   F   F
  0.2500000000000000  0.7842769586426328  0.4472094799095390   T   T   T
  0.2500000000000000  0.5120099780000018  0.1068399999999983   F   F   F
  0.2500000000000000  0.5136721265377703  0.3264131468700368   T   T   T
-0.0574929106060290  0.0308639091850143  0.4732777707551415   T   T   T
  0.0000000000000000  0.5000000000000000  0.2189999969999974   F   F   F
  0.5574929106060290  0.0308639091850143  0.4732777707551415   T   T   T
  0.5000000000000000  0.5000000000000000  0.2189999969999974   F   F   F
  0.0000000000000000  0.0000000000000000  0.2189999969999974   F   F   F
  0.0575100652447704  0.4692103557750934  0.4733084597596518   T   T   T
  0.5000000000000000  0.0000000000000000  0.2189999969999974   F   F   F
  0.4424899347552297  0.4692103557750934  0.4733084597596518   T   T   T
  0.2500000000000000  0.0277400000000014  0.0110600000000005   F   F   F
  0.2500000000000000  0.0277099999999990  0.2315299959999990   F   F   F
  0.7500000000000000  0.4722599979999984  0.2314999999999969   F   F   F
  0.7500000000000000  0.4722900089999982  0.0110799999999998   F   F   F
  0.7500000000000000  0.9433945793964041  0.4283109619600028   T   T   T
  0.7500000000000000  0.9722899790000028  0.2064699980000029   F   F   F
  0.2500000000000000  0.5277400020000016  0.2064999939999979   F   F   F
  0.2500000000000000  0.5566217165151475  0.4283746100846041   T   T   T
在优化outcar结果里面：

POSITION                                       TOTAL-FORCE (eV/Angst)
-----------------------------------------------------------------------------------
      1.88530      2.00494      0.00000         0.000000     -5.098157    -0.081823
      1.88530      2.00456      3.93267         0.000000      0.811526     -0.398017
      1.88530      0.09295      1.97517         0.000000     -0.316509     -4.473482
      1.88530      0.13336      5.86231         0.000000     -0.009223      0.005466
      5.65590      1.86456      3.93250         0.000000     -0.803723     -0.398596
      5.65590      1.86494      0.00018         0.000000      5.094488     -0.083175
      5.65590      3.73604      5.86180         0.000000      0.009235      0.004896
      5.65590      3.77586      1.97552         0.000000      0.321206     -4.475273
      5.65590      5.53847      7.97806         0.000000     -0.008509      0.016186
      5.65590      5.73444      3.88094         0.000000     -0.964166     -0.124007
      5.65590      7.63218      5.82294         0.000000     -0.033566     -0.000981
      5.65590      7.64536      1.90559         0.000000     -0.330897      0.921050
      1.88530      5.87444      3.88112         0.000000      0.953343     -0.125655
      1.88530      6.07055      7.97856         0.000000      0.000961      0.013845
      1.88530      3.96245      1.90595         0.000000      0.333991      0.929429
      1.88530      3.97655      5.82309         0.000000      0.033371      0.000669
      7.10701      0.23848      8.44284        -0.002915     -0.001649      0.007800
      0.00000      3.86950      3.90681        -0.086346      0.042062     -0.338047
      4.20479      0.23848      8.44284         0.002915     -0.001649      0.007800
      3.77060      3.86950      3.90681         0.086346      0.042062     -0.338047
      0.00000      0.00000      3.90681         0.086992     -0.041792     -0.337944
      0.43410      3.63148      8.44358         0.011632      0.005666      0.001942
      3.77060      0.00000      3.90681        -0.086992     -0.041792     -0.337944
      3.33650      3.63148      8.44358        -0.011632      0.005666      0.001942
      1.88530      0.21468      0.19730         0.000000      5.038306      2.986906
      1.88530      0.21445      4.13033         0.000000     -0.467185      2.745592
      5.65590      3.65482      4.12980         0.000000      0.468972      2.755327
      5.65590      3.65505      0.19766         0.000000     -5.038394      2.983682
      5.65590      7.30162      7.64188         0.000000     -0.022045     -0.000543
      5.65590      7.52455      3.68328         0.000000      0.848326     -0.936172
      1.88530      4.08418      3.68382         0.000000     -0.851840     -0.945213
      1.88530      4.30719      7.64203         0.000000      0.021915      0.012387
-----------------------------------------------------------------------------------
    total drift:                                0.000000      0.001047      0.021587

里面有几个原子受力特别大，不停的弛豫也无法降低下来，   1.88530      2.00494      0.00000         0.000000   -5.098157     -0.081823 ；      5.65590      1.86494      0.00018         0.000000      5.094488     -0.083175。这些原子对应的是设置的F F F 原子，求助各位大神，如何解决。
回复此楼

» 收录本帖的淘帖专辑推荐

答疑解惑

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
可以可以啊
1楼 2015-07-02 15:59:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪雨星风17

木虫 (小有名气)

无
引用回帖:
2楼: Originally posted by hs0218 at 2015-07-02 20:50:04
ENCUT为什么选这么高！
还有看不到你的INCAR，也不好分析。
你这个某些原子受力确实大的离谱。

INCAR 选的是候博士的手册上的
SYSTEM = name
ENCUT = 750
ISMEAR = 1; SIGMA = 0.20
ISTART = 0; ICHARG = 2
EDIFF = 1E-5; EDIFFG = -1.0E-3
NSW = 60; IBRION = 2
POTIM = 0.1
PREC= Accurate
ALGO = F
LWAVE = .FALSE.
LCHARG = .FALSE
KPOINTS是7 7 1,跟encut选的高有关系么？
赞一下 回复此楼
可以可以啊
3楼2015-07-03 09:34:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

hs0218

金虫 (正式写手)

DDr.

【答案】应助回帖

★
感谢参与，应助指数 +1
liliangfang: 金币+1, 谢谢交流 2015-07-03 07:54:38
ENCUT为什么选这么高！
还有看不到你的INCAR，也不好分析。
你这个某些原子受力确实大的离谱。
赞一下 回复此楼
路漫漫其修远兮
2楼2015-07-02 20:50:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liutaifeng

新虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
雪雨星风17: 金币+3, ★有帮助 2015-07-03 22:09:28
liliangfang: 金币+1, 谢谢交流 2015-07-04 08:12:26
表面构造是否是非极性面 可以参考一下文献看看 文献上有很多wo3表面的
赞一下 回复此楼
4楼2015-07-03 11:01:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hs0218

金虫 (正式写手)

DDr.

【答案】应助回帖

★ ★ ★ ★
雪雨星风17: 金币+3, ★有帮助 2015-07-03 22:09:39
liliangfang: 金币+1, 谢谢交流 2015-07-04 08:12:32
我用你的POSCAR和你的INCAR计算了下，也是原子力很大。结果和你差不多
我没去查你固定的是哪些原子，是不是基体原子固定，优化靠近真空的原子？
赞一下 回复此楼
路漫漫其修远兮
5楼2015-07-03 17:55:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-29 18:04 by 4FFAWE8HcgUD
[硕博家园] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 10/500 2026-08-29 16:49 by zICmwzsBXjbN
[论文投稿] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 14:13 by jCd0dEvKHShX
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[论文投稿] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 ASdOkHsho7FD 2026-08-28 4/200 2026-08-29 13:51 by jCd0dEvKHShX
[考研] 售一区SCI文章T0P，我:8O.551.O54,科目全,可十急 +4 ASdOkHsho7FD 2026-08-28 6/300 2026-08-29 11:40 by jCd0dEvKHShX
[基金申请] 中青基了要发朋友圈吗？ +3 349506619 2026-08-28 3/150 2026-08-29 10:49 by 木风9012
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +12 Kittylucky 2026-08-27 13/650 2026-08-29 00:04 by superceng
[基金申请] 基金未中，这种答复是模板吗？ +5 zhaosm1982 2026-08-27 6/300 2026-08-27 16:00 by lfy8008
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了，中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
关闭
请填处理意见
关闭