24小时热门版块排行榜    

查看: 3218  |  回复: 19

lgxyz

铜虫 (正式写手)

[求助] 氧空位CeO2(111)的优化

各位,我最近在优化具有一个氧空位的CeO2(111)负载一个Au原子的结构(Au位于氧空位上),计算吸附能。发现有一个问题:第一次在优化clean的具有一个氧空位的CeO2(111)时是直接从优化过的CeO2体相结构切出(111)面,然后除去一个氧原子,能量为-182.853254eV;第二次是从优化好的Au/Ceo2(111)结构上除去Au原子,再优化,能量为-184.311185eV;二者结构上只是键长键角稍不一样,能量却相差比较大;按道理来讲,应该取能量更低的结构,但是一般我们都是按第一种办法来优化clean的表面。
具有两个氧空位的情况,这两种优化结果能量相差很小。
请教大家是怎样看待这个问题。非常感谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主


gzqdyouxia(金币+1): 谢谢交流 2011-09-23 11:52:06
引用回帖:
1楼: Originally posted by lgxyz at 2011-09-21 16:39:30:
各位,我最近在优化具有一个氧空位的CeO2(111)负载一个Au原子的结构(Au位于氧空位上),计算吸附能。发现有一个问题:第一次在优化clean的具有一个氧空位的CeO2(111)时是直接从优化过的CeO2体相结构切出(11 ...

总能差别会有这么大?怪哉。第一种情况这算氧空位形成能多大?算算DOS和spin polarized charge density,证实这种差别的存在和合理性。

[ Last edited by cenwanglai on 2011-9-21 at 16:50 ]
2楼2011-09-21 16:49:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)

引用回帖:
2楼: Originally posted by cenwanglai at 2011-09-21 16:49:29:
总能差别会有这么大?怪哉。第一种情况这算氧空位形成能多大?算算DOS和spin polarized charge density,证实这种差别的存在和合理性。
[ Last edited by cenwanglai on 2011-9-21 at 16:50 ]

第一种氧空位的形成能有3.95eV,第二种有2.49eV
3楼2011-09-21 16:53:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)

引用回帖:
2楼: Originally posted by cenwanglai at 2011-09-21 16:49:29:
总能差别会有这么大?怪哉。第一种情况这算氧空位形成能多大?算算DOS和spin polarized charge density,证实这种差别的存在和合理性。
[ Last edited by cenwanglai on 2011-9-21 at 16:50 ]

还有一种情况是同一氧空位吸附不同的金属,优化后,除去金属原子,再来优化具有一个氧空位clean CeO2(111),得到的能量也不一样,相差大概0.13eV,不知道是不是由于我建的slab只有6层原子,厚度不够,导致能量浮动。
4楼2011-09-21 17:02:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主

★ ★
gzqdyouxia(金币+2): 谢谢提示 2011-09-23 12:07:46
引用回帖:
4楼: Originally posted by lgxyz at 2011-09-21 17:02:04:
还有一种情况是同一氧空位吸附不同的金属,优化后,除去金属原子,再来优化具有一个氧空位clean CeO2(111),得到的能量也不一样,相差大概0.13eV,不知道是不是由于我建的slab只有6层原子,厚度不够,导致能量 ...

slab有多大?六层原子,垂直真空层方向多大?3*3?第一种情况氧空位形成能偏大,如果是2*2的话似乎差不多。
5楼2011-09-21 17:17:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)

引用回帖:
5楼: Originally posted by cenwanglai at 2011-09-21 17:17:36:
slab有多大?六层原子,垂直真空层方向多大?3*3?第一种情况氧空位形成能偏大,如果是2*2的话似乎差不多。

2x2 的超胞 6层原子(两层O-Ce-O重复)真空层12A
6楼2011-09-21 17:21:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dwma

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★
cenwanglai(金币+2): 有可能是这样的。需要认为调节空穴周围的几何对称性拉降低能量。 2011-09-21 19:25:24
lgxyz(金币+5): 2011-09-22 15:49:46
lgxyz(金币+15): 2011-09-22 15:50:40
From my viewpoint, for the first case, the sixfold (or threefold) symmetry has been maintained during the optimization. However, for the second case, the adsorption of Au has broken that symmetry, thus you can get the real ground state of the reduced Ceria. The situation is same like Jahn-Teller effect. You can also get the real ground state by optimizing the clean surface if you displace one of three Ce cations nearest the O vacancy before the optimization. If you have problems, please contact me.
7楼2011-09-21 19:19:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kjsun

木虫 (正式写手)

★ ★
贺仪(金币+2): 多谢指教! 2011-09-22 09:50:02
引用回帖:
7楼: Originally posted by dwma at 2011-09-21 19:19:07:
From my viewpoint, for the first case, the sixfold (or threefold) symmetry has been maintained during the optimization. However, for the second case, the adsorption of Au has broken that symmetry,  ...

7楼正解,
我也遇到了同样的问题。
对称性的原因
如果是用vasp算的话
设置 ISYM=0
8楼2011-09-22 09:32:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lgxyz

铜虫 (正式写手)

引用回帖:
7楼: Originally posted by dwma at 2011-09-21 19:19:07:
From my viewpoint, for the first case, the sixfold (or threefold) symmetry has been maintained during the optimization. However, for the second case, the adsorption of Au has broken that symmetry,  ...

加了ISYM=0,结果第一种建模方法算的能量为-184.172738eV,与第二种方法的能量(-184.311393eV)相比还是有点差别,但明显比采用默认对称性来算的结果要好一些。请你再给些建议!非常感谢!
9楼2011-09-23 09:34:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kjsun

木虫 (正式写手)

★ ★
gzqdyouxia(金币+2): 谢谢指教 2011-09-23 12:07:25
引用回帖:
9楼: Originally posted by lgxyz at 2011-09-23 09:34:52:
加了ISYM=0,结果第一种建模方法算的能量为-184.172738eV,与第二种方法的能量(-184.311393eV)相比还是有点差别,但明显比采用默认对称性来算的结果要好一些。请你再给些建议!非常感谢!

是这样的,
CeO2上的每一个氧空位会产生2个多余的电子,
这两个电子会分布在周围的Ce上,形成Ce3+ (原来是Ce4+)
两个电子落位的不同,就会形成不同的结构(仔细观察:Ce3+ -O 的键长比Ce4+-O略长0.15A左右),也就对应的不同的能量。
你这两次运算,应该是碰到了电子的不同的落位了。
10楼2011-09-23 11:54:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lgxyz 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 PbwVEMK5haaE 2026-08-25 4/200 2026-08-26 03:20 by cNXvBfCpiZOM
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 Qj98d49dHFb9 2026-08-25 4/200 2026-08-26 03:08 by cNXvBfCpiZOM
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +17 木水思豆 2026-08-25 20/1000 2026-08-26 01:26 by chengyan1220
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 有没有大神帮我看看基金代码 28+4 1234567wang 2026-08-24 10/500 2026-08-25 19:15 by lfy8008
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +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 南海小哥
[教师之家] 导师吐槽:我怎么摊上了这么个极品研究生! +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 范进中举一文的中心思想 +7 炎黄贵胄 2026-08-22 8/400 2026-08-25 08:48 by ZJTJZ
[基金申请] filecode,4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤?
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办? +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗? +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天,20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
请填处理意见