24小时热门版块排行榜    

查看: 2779  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085601材料工程315分求调剂 +3 yang_0104 2026-03-15 3/150 2026-03-15 10:58 by peike
[考研] 289求调剂 +5 步川酷紫123 2026-03-11 5/250 2026-03-15 00:45 by kruisytel
[考研] 中科大材料与化工319求调剂 +3 孟鑫材料 2026-03-14 3/150 2026-03-14 20:10 by ms629
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 267一志愿南京工业大学0817化工求调剂 +5 SUICHILD 2026-03-12 5/250 2026-03-14 14:53 by jean5056
[考研] 290求调剂 +4 @将就将就看 2026-03-10 8/400 2026-03-14 14:23 by 千千运气
[考研] 331求调剂(0703有机化学 +5 ZY-05 2026-03-13 6/300 2026-03-14 10:51 by Jy?
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 求调剂 +3 清风问长安 2026-03-09 3/150 2026-03-14 02:15 by JourneyLucky
[考研] 求调剂 +6 yfihxh 2026-03-09 6/300 2026-03-14 01:18 by JourneyLucky
[考研] 复试调剂 +9 Copy267 2026-03-10 9/450 2026-03-13 23:45 by userper
[考研] 341求调剂 +4 番茄头--- 2026-03-10 4/200 2026-03-13 23:12 by JourneyLucky
[考研] [0860]321分求调剂,ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 工科,求调剂 +3 我887 2026-03-11 3/150 2026-03-13 21:39 by JourneyLucky
[考研] 26调剂/材料科学与工程/总分295/求收留 +9 2026调剂侠 2026-03-12 9/450 2026-03-13 20:46 by 18595523086
[考研] 289求调剂 +3 李政莹 2026-03-12 3/150 2026-03-13 11:02 by 求调剂zz
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 081200-11408-276学硕求调剂 +3 崔wj 2026-03-12 4/200 2026-03-12 19:33 by 求调剂zz
[考研] 哈工大材料324求调剂 +6 闫旭东 2026-03-10 8/400 2026-03-10 22:49 by 星空星月
信息提示
请填处理意见