24小时热门版块排行榜    

查看: 3275  |  回复: 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-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 07:58 by C79jjtjAKjEn
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:46 by C79jjtjAKjEn
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:34 by C79jjtjAKjEn
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-14 3/150 2026-09-15 07:14 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-15 06:50 by 5BDX0d0WFp7t
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:38 by 5BDX0d0WFp7t
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 06:01 by 5BDX0d0WFp7t
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 05:50 by 5BDX0d0WFp7t
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:49 by 5BDX0d0WFp7t
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-15 05:34 by 5BDX0d0WFp7t
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:26 by 5BDX0d0WFp7t
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 QUjhNVAcOSff 2026-09-13 8/400 2026-09-15 05:01 by mibUvS8DDCwf
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 4/200 2026-09-14 23:24 by vZfe6xYu34yj
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
信息提示
请填处理意见