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

24小时热门版块排行榜    

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

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的回帖

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的回帖

xyj880506

银虫 (小有名气)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
雪雨星风17: 金币+20, ★★★很有帮助 2015-07-03 22:09:06
liliangfang: 金币+2, 谢谢交流 2015-07-04 08:12:40
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
SIGMA = 0.1
赞一下 回复此楼
6楼2015-07-03 21:25:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk

至尊木虫 (著名写手)
引用回帖:
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
...

算slab不用固定底层原子?能否说明原因？
赞一下 回复此楼
我们最深的恐惧不是我们能力不够，我们最深的恐惧是能力超越了极限。
7楼2015-07-03 21:31:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

雪雨星风17

木虫 (小有名气)

无
引用回帖:
5楼: Originally posted by hs0218 at 2015-07-03 17:55:09
我用你的POSCAR和你的INCAR计算了下，也是原子力很大。结果和你差不多
我没去查你固定的是哪些原子，是不是基体原子固定，优化靠近真空的原子？

是的是的，我就是这样做的
赞一下 回复此楼
可以可以啊
8楼2015-07-03 22:07:27
已阅   回复此楼   关注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的回帖

xyj880506

银虫 (小有名气)

【答案】应助回帖

★
liliangfang: 金币+1, 谢谢交流 2015-07-04 08:12:49
引用回帖:
7楼: Originally posted by magicmonk at 2015-07-03 21:31:30
算slab不用固定底层原子?能否说明原因？...

以前有一篇文章说不固定底层得到的能带和实验的更接近

[ 发自小木虫客户端 ]
赞一下 回复此楼
10楼2015-07-03 23:24:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 雪雨星风17 的主题更新
131/2返回列表12下一页
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭