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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 【讨论】[讨论]奇怪的解离吸附
51/1返回列表
查看: 585  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

dazao

金虫 (正式写手)

[交流] 【讨论】[讨论]奇怪的解离吸附

我在做CO在γ-U（0 0 1 ）面吸附的时候 CO垂直于表面 C原子在下 且是在中心位吸附 slab 3层 vacuum 12 埃   吸附完成后 CO解离 奇怪的是O原子跑到下一个周期supercell的下表面去了   
我把INCAR 和 POSCAR 贴给大家 K是5X5X1 PBE势
   
   ISTART = 0; ICHARG = 2
   ENCUT  = 370
   ISMEAR =  2  ; SIGMA = 0.26
   LREAL=FLASE
NELM=200
NSW=100
IBRION=1
EDIFFG=-0.05
IDIPOL=3
LDIPOL=TRUE
  NGX=44
  NGY=44
  NGZ=98

POSCAT
gama u (1 0 0) CO vert2
1
       6.850000000000000       0.000000000000000       0.000000000000000
       0.000000000000000       6.850000000000000       0.000000000000000
       0.000000000000000       0.000000000000000      15.425000000000001
1 1 12
s
d
     0.2500000000000000   0.2500000000000000   0.2617999911308289 T T T
     0.2500000000000000   0.2500000000000000   0.3359000086784363 T T T
     0.0000000000000000   0.0000000000000000  -0.0000000000000000 F F F
     0.0000000000000000   0.0000000000000000   0.2241546518394000 T T T
     0.0000000000000000   0.4999999999999999  -0.0000000000000000 F F F
     0.0000000000000000   0.4999999999999999   0.2241546518394000 T T T
     0.4999999999999999   0.0000000000000000  -0.0000000000000000 F F F
     0.4999999999999999   0.0000000000000000   0.2241546518394000 T T T
     0.4999999999999999   0.5000000000000000   0.0000000000000000 F F F
     0.4999999999999999   0.5000000000000000   0.2241546518394000 T T T
     0.2500000000000000   0.2500000000000000   0.1120773259197000 F F F
     0.2500000000000000   0.7500000000000000   0.1120773259197000 F F F
     0.7500000000000000   0.2500000000000000   0.1120773259197000 F F F
     0.7500000000000000   0.7499999999999999   0.1120773259197000 F F F

弛豫后的CONTCAR


第一个是C原子 第二个是铀原子
弛豫后的CONTCAR
gama u (1 0 0) CO vert2                 
1.00000000000000000
     6.8499999999999996     .0000000000000000     .0000000000000000
      .0000000000000000    6.8499999999999996     .0000000000000000
      .0000000000000000     .0000000000000000   15.4250000000000007
   1   1  12
Selective dynamics
Direct
   .2500000000000000   .2500000000000000   .2662342893262475   T   T   T
   .2500000000000000   .2500000000000000   .9539376694716014   T   T   T
   .0000000000000000   .0000000000000000   .0000000000000000   F   F   F
   .0134388848036315   .0134388848036318   .2166569367679411   T   T   T
   .0000000000000000   .5000000000000000   .0000000000000000   F   F   F
   .0134388848036309   .4865611151963711   .2166569367679391   T   T   T
   .5000000000000000   .0000000000000000   .0000000000000000   F   F   F
   .4865611151963711   .0134388848036314   .2166569367679408   T   T   T
   .5000000000000000   .5000000000000000   .0000000000000000   F   F   F
   .4865611151963711   .4865611151963711   .2166569367679393   T   T   T
   .2500000000000000   .2500000000000000   .1120773259196994   F   F   F
   .2500000000000000   .7500000000000000   .1120773259196994   F   F   F
   .7500000000000000   .2500000000000000   .1120773259196994   F   F   F
   .7500000000000000   .7500000000000000   .1120773259196994   F   F   F

大家注意O的位置

[ Last edited by wuchenwf on 2009-6-23 at 19:50 ]
回复此楼

» 猜你喜欢
1楼2009-04-04 13:29:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dazao

金虫 (正式写手)
引用回帖:
Originally posted by stractor at 2009-4-10 23:08:
你用三层slab,而弛豫时只允许表面那一层动，觉得你的表面不合理，建议增加1-2层slab,而弛豫3-4层原子。你的O原子向上动了很多，对于U, 12A的真空层是否可以？

[ Last edited by stractor on 2009-4-10 at 23: ...

12A完全可以没问题 测试过的
slab的取法是参照文献选的
赞一下 回复此楼
10楼2009-04-13 08:52:07
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答

hakuna

木虫 (知名作家)
★ ★ ★
fegg7502(金币+3,VIP+0):鼓励交流！ 4-5 16:55
多半是你初始结构不合理，检查一下你的初始结构中原子之间的距离是否合理，如果某个距离给的初始值太小，优化一下会散掉。
你也可查看一下优化过程，看看看是在哪一步C-O散掉的，前一步发生了什么事情，这样就知道什么原因了。
查一下C到表面的距离以及C-O键长。

[ Last edited by hakuna on 2009-4-4 at 19:35 ]
赞一下(10人) 回复此楼
2楼2009-04-04 19:32:20
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dazao

金虫 (正式写手)
我又看了一遍OUTCAR  发现中间过程有O嵌入U层中
师姐刚回来 我请教了 还没解决
大家还有什么解决办法没？

我现在把OC交换了位置 还在计算中
赞一下 回复此楼
3楼2009-04-06 09:33:57
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dazao

金虫 (正式写手)
★ ★
zxzj05(金币+2,VIP+0):3Q！ 期待更多讨论！ 4-7 09:59
现在看可能是偶极修正位置对吸附分子造成了影响。
我没有设置偶极修正的位置（DIPOL缺省）缺省值应该是SUPERCELL的中心
这个位置又和初始O原子的位置很接近 弛豫过程中偶极修正的电场可能给氧原子强烈的作用
把真空层放大到20埃 就没有解离吸附了（为了节约时间K取1x1x1）
又做了个手动设置DIPOL的试验 还在等结果  出了结果告诉大家
赞一下(10人) 回复此楼
4楼2009-04-06 13:59:12
 已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 10 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭