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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » 【讨论】[讨论]奇怪的解离吸附
51/1返回列表
查看: 657  |  回复: 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的回帖

stractor

金虫 (著名写手)
★
zdhlover(金币+1,VIP+0):谢谢交流，欢迎常来 4-11 01:05
你用三层slab,而弛豫时只允许表面那一层动，觉得你的表面不合理，建议增加1-2层slab,而弛豫3-4层原子。你的O原子向上动了很多，对于U, 12A的真空层是否可以？

[ Last edited by stractor on 2009-4-10 at 23:12 ]
赞一下(10人) 回复此楼
9楼2009-04-10 23:08:55
已阅   回复此楼   关注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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 申博/考博 +4 啃面包的小书虫 2026-04-17 8/400 2026-04-21 16:26 by 啃面包的小书虫
[考研] 一志愿A区211，22408 321求调剂 +7 随心所欲☆ 2026-04-15 8/400 2026-04-21 08:22 by Equinoxhua
[考研] 295分求调剂 +6 ?要上岸? 2026-04-17 6/300 2026-04-21 08:18 by Equinoxhua
[考研] 一志愿中科大材料与化工，353分还有调剂学校吗 +11 否极泰来2026 2026-04-15 13/650 2026-04-20 22:31 by Equinoxhua
[考研] 085600材料与化工调剂 5+3 孜孜不倦2002 2026-04-19 6/300 2026-04-20 21:25 by babero
[论文投稿] 期刊推荐 +3 材料研究生 2026-04-15 5/250 2026-04-20 16:02 by 豆豆7758
[考研] 337求调剂 +3 jyz04 2026-04-18 3/150 2026-04-20 12:24 by 研可安
[考博] 申博 +3 Xyyx. 2026-04-18 3/150 2026-04-20 10:44 by YuY66
[考博] 湖南大学刘巧玲课题组2026年第二批次博士研究生招生信息 +3 南风观火 2026-04-18 5/250 2026-04-20 10:13 by 南风观火
[考研] 求计算机方向调剂 +3 Toffee2 2026-04-16 6/300 2026-04-19 22:37 by ll叶
[考研] 294求调剂 +8 淡然654321 2026-04-17 9/450 2026-04-19 19:51 by Equinoxhua
[考研] 304求调剂 +8 castLight 2026-04-16 8/400 2026-04-19 17:14 by 中豫男
[考研] 求调剂 +10 小聂爱学习 2026-04-16 12/600 2026-04-19 16:51 by 中豫男
[考研] 307中医考研调剂 +9 于以采蘩 2026-04-14 9/450 2026-04-19 08:41 by 烟雨流涯
[考研] 294求调剂 +15 淡然654321 2026-04-15 15/750 2026-04-19 08:20 by cuisz
[考研] 0854求调剂 +23 门路摸摸 2026-04-15 27/1350 2026-04-19 01:59 by 烟雨流涯
[考研] 300求调剂 +12 橙a777 2026-04-15 12/600 2026-04-18 23:51 by 路病情
[考研] 260求调剂 +4 Zyt1314520.. 2026-04-17 5/250 2026-04-18 08:28 by babysonlkd
[考研] 急需调剂 +9 绝不放弃22 2026-04-15 10/500 2026-04-18 08:09 by chixmc
[有机交流] 二苯甲酮酸类衍生物 50+3 小白爱主人 2026-04-17 6/300 2026-04-17 18:47 by kf2781974
信息提示
关闭
请填处理意见
关闭