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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » WO3的slab模型问题
51/1返回列表
查看: 3157  |  回复: 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

木虫 (小有名气)

无
引用回帖:
6楼: Originally posted by xyj880506 at 2015-07-03 21:25:44
http://pubs.acs.org/doi/abs/10.1021/jp410244c
先优化单胞大小再切slab model，截断能400-500就行，Z方向大小扩大一倍，真空层用15，不用固定底层原子。
NSW 调大
EDIFFG用-0.03到-0.05即可
KPOINTS用3 3 1
...

哇，谢谢谢谢，我仔细瞅瞅
赞一下 回复此楼
可以可以啊
9楼2015-07-03 22:08:55
已阅   回复此楼   关注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的回帖

雪雨星风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的回帖

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的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +10 Kittylucky 2026-08-27 11/550 2026-08-28 10:45 by Kittylucky
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +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
[基金申请] 基金不中，共勉 +10 eulota 2026-08-26 10/500 2026-08-27 17:23 by lqllinqiaoli
[文学芳草园] 梦想 +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 北京莱茵编辑
[基金申请] 我不理解！ +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
[基金申请] 系统查不到 +10 董八千 2026-08-26 10/500 2026-08-26 16:30 by Equinoxhua
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 系统进不去 +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
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 明天应该可查了！？ +6 chengyan1220 2026-08-23 6/300 2026-08-25 19:45 by zfd97
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
信息提示
关闭
请填处理意见
关闭