24小时热门版块排行榜

返回列表

cenwanglai

荣誉版主 (知名作家)

[交流] 【文献学习】JACS上一篇纯DFT计算的文献： CI-NEB计算，Leuwdin布局分析，态密度分析

根据第一性原理板块版务讨论精神，尝试开展文献讨论活动。水平有限，欢迎大家批评指导。尽管如此，我还是希望以讨论而不是求助的语气来展开。希望更多虫友参与此项活动：贴出自己阅读的有意思的文献，拿出来讨论一下，供大家学习参考！

【标题】：Reaction Mechanisms for the CO Oxidation on Au/CeO2 Catalysts: Activity of Substitutional Au3+/Au+ cations and Deactivation of Supported Au+ Adatoms
【作者】：Matteo Farnesi Camellone and Stefano Fabris J.AM.CHEM.SOC.2009, 131.10473-10483

【使用软件】：Quantum-ESPRESSO
【所用赝势】：UPP Au [Xe]4f14|5d9.5-6s1-6p0.5 Ce [Kr]4d10|5s2-5p6-4f1-5d1-6s2 O [He]|2s2-2p4
【分析方法】：Leuwdin population analysis ; DOS and PDOS; CI-NEB

下载信息：
邮箱：emuch_dft@163.com
密码：cenwanglai
主题：DFT001_JACs
下载链接：http://eg1a164.mail.163.com/js3/ ... jitFoppdXRaKeHaDovp

内容提要

这篇文章用GGA/GGA+U方法研究了Au单个原子在CeO2（111）表面上的吸附，认为Au在O-O桥位和O顶位上的化学吸附，存在Au向O有电子转移，使得Au带正电荷。这是CO能够在 Au上发生化学吸附的原因。桥位吸附优于顶位吸附。 CO在化学计量CeO2(111)表面上为物理吸附。

桥位吸附的Au吸附CO分子后，前者被CO拉到O顶位，形成Au-C-O表面垂直结构。接着CO在 Au上溢流，与临近顶位O反应生成CO2,同时产生一个空位。Au被拉到O空位上，形成负价Au离子，失去CO吸附能力而发生失活。

文章认为，Au替位掺杂（取代一个次表层Ce原子）活性高于吸附掺杂（Au近吸附于CeO2表面）；而且单个Au原子桥位或者顶位吸附结构在CO催化氧化过程中都存在失活现象。文章认为Au cluster吸附可以限制Au的移动，从而不会占据反应过程中的O空位，氧空位可以从气相复氧，催化过程可以继续。

我的问题(一)

这个是用CNEB方法算的CO在Au@CeO2(111) 上吸附并被CeO2表面O氧化的反应路径。前面计算方法中交代了用20images.

图上列出的IS和FS，以及中间的Splillover, MS,Oxidation几个阶段都太有故事了吧。仅仅通过线性插值能得出这样的过程吗？

这篇文章中用PBE，别的文献也有用PW91的。对于含f电子的稀土金属氧化物和贵金属体系，选择赝势有什么考虑呢？

NEW

2010-12-16
写邮件问了文章通讯作者Fabris，他的回复如下：
The MS image was initially found as a metastable adsorption configuration in a full NEB calculation from the IS to the FS. With this information, we then split the search for the minimum-energy path into two separate calculations, one for the IS-MS and one for the MS-FS, as you correctly guessed.

同时问问，为什么这篇文章可以发JACS呢？好像不是很牛逼呢。或者是JACS也不是很牛？（文章写得确实好，计算是不难，难在作者的谋篇布局，论述说理）

google了一下 Stefano Fabris. 2009年他拿这篇文章到天津大学讲座。这人至今还不是Democritos的终生雇员，一个tenure track。

我的问题(二)

上表是Support Information中的Leuwdin（eu表示o上两点）五种结构的布局分析最相关的三个原子的结果，包括一个Au原子，和两个Ce原子。其它原子布局可以忽略。

Support Information中把Au布局结果分成两组，9.5以下的两个数取平均，认为是Au3+离子；9.5以上三个取平均，认为是Au+离子。Ce也分为两组，小于1的认为是Ce4+，大于1的，认为是Ce3+. 为什么可以这样求平均呢？

实际上元素周期表中Au外层为5d10-6s1。而文献中用5d9.5-6s1-6p0.5。为什么可以这样选取呢？

我的问题(三)待完善ing

请仔细阅读图下面的说明文字。先列出两个相关联的问题：
（1）上图中左边a和d represent the bonding charge integrated in planes perpendicular to the surface and plotted as a function of the height from the surface.
（2）中间图b和e display the bonding charge plotted at the value of 0.06 electrons/A^3. Electron accumulation and depletion are represented by red and blueareas, respectively.
上面两处蓝色字体是什么意思呢？bonding charge是怎么得到的？差分电荷？

上面图片中的文字是在文献中扣出来的。这里交代了如何计算bonding charge. 即在同一个box里分别计算AB,A，B三个部分的电荷密度。这个可不是差分电荷密度呢。火力集中在A,B的界面上。

总结起来，分析电荷密度，研究成键情况，至少有有三种方法处理电荷密度了：差分电荷密度，局域函数，和这里的bonding charge.

注意到其中能来那个计算Eads=Etot-(Esub-Ex)。为什么这样呢？不是应该是Etot-(Esub+Ex)吗？这是不是ESPRESSO软件的特点之一？

文中说Ex是单个原子Au在气相中的能量。但是没有交代怎么算的。那么Etot,Esub和Ex是不是要用同一个box,同一k设置，同一截断能呢？我自己试过，单个Au原子在不同情况下的基态（截断能和k设置合理的情况下）能量差别不大。

[ Last edited by cenwanglai on 2011-6-27 at 17:35 ]

24小时热门版块排行榜

[交流] 【文献学习】JACS上一篇纯DFT计算的文献： CI-NEB计算，Leuwdin布局分析，态密度分析

» 收录本帖的淘帖专辑推荐

» 本帖已获得的红花（最新10朵）

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

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

» 抢金币啦！回帖就可以得到: 查看全部散金贴

» 本帖已获得的红花（最新10朵）

» 本帖已获得的红花（最新10朵）

» 抢金币啦！回帖就可以得到:

查看全部散金贴