24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 10898  |  回复: 30
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

[交流] 【文献学习】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方法算的COAu@CeO2(111) 上吸附并被CeO2表面O氧化的反应路径。前面计算方法中交代了用20images.

图上列出的ISFS,以及中间的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 ]
回复此楼

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

仿真建模与计算 VASP 第一原理资料汇编 vasp
第一性原理和电化学 first principle 第一性 软件 程序拟合等
vasp 应该这么读文献 计算化学 MS软件
收藏

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

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

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

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

查看全部散金贴

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

liyong.gan

木虫 (正式写手)


★ ★
cenwanglai(金币+3):谢谢回复! 2010-12-06 13:51:37
zzy870720z(金币+2):鼓励参加文献讨论,这个经验需要长时间积累啊,呵呵 2010-12-06 15:54:21
引用回帖:
Originally posted by cenwanglai at 2010-12-06 12:56:19:


那针对我所列的这个过渡态计算,初始结构与最终结构差异非常大。如何分段来设计反应路径呢?

一般给定端点,是通过线性插值来找中间点的。如果实现指定若干中间点,再在这些点之间插值,各个image之间的“ ...

那么就要根据经验看你的反应路径中有没有必须经过的某个态,找出来优化,把其当着中间态。
5楼2010-12-06 13:15:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 31 个回答
★ ★
zzy870720z(金币+2):谢谢讨论 2010-12-06 12:04:44
cenwanglai(金币+3):呵呵,谢谢回复。这边讨论着。我写邮件问问文章作者。 2010-12-06 12:45:25
cenwanglai(金币+1):问了文章通讯作者Fabris,中间过程确实是算出来的。后面我附上他的回信。 2010-12-16 11:20:47
我想这些中间过程不是“设计”出来的,而是确实算出来的。image选的多一些,也许确实可以得出这些结构。
对于重金属体系,关键的我觉得这些PBE和PW91其实都是一回事,除非采用杂化泛函。主要是还要考虑强关联效应。
2楼2010-12-06 11:46:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liyong.gan

木虫 (正式写手)


★ ★
cenwanglai(金币+3):谢谢回复。 2010-12-06 12:52:30
zzy870720z(金币+2):鼓励参加文献讨论 2010-12-06 15:53:36
对于有f电子的Ce元素,我看到有文献报道用的+u计算;
对于含有20个image的路径,算是一个比较复杂的路径了,一般不是只设了一个IS和FS,而是把整个过程分为好几个段,然后一段段连接起来的。

[ Last edited by liyong.gan on 2010-12-6 at 12:47 ]
3楼2010-12-06 12:44:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youzhizhe(金币+1):鼓励讨论。 2010-12-06 19:51:58
引用回帖:
Originally posted by liyong.gan at 2010-12-06 12:44:23:
对于有f电子的Ce元素,我看到有文献报道用的+u计算;
对于含有20个image的路径,算是一个比较复杂的路径了,一般不是只设了一个IS和FS,而是把整个过程分为好几个段,然后一段段连接起来的。

[ Last edited  ...

那针对我所列的这个过渡态计算,初始结构与最终结构差异非常大。如何分段来设计反应路径呢?

一般给定端点,是通过线性插值来找中间点的。如果实现指定若干中间点,再在这些点之间插值,各个image之间的“距离”是不相等的。做CNEB计算的时候,需要把事先设计的中间点与插值点练成一条chain,然后在一个任务中跑吗?

[ Last edited by cenwanglai on 2010-12-6 at 13:49 ]
4楼2010-12-06 12:56:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿太原理工安全工程300分,求调剂 +3 0857求调剂. 2026-03-24 3/150 2026-03-27 01:18 by 求知若渴lz
[考研] 349求调剂 +5 杰斯塔里斯 2026-03-21 5/250 2026-03-27 00:31 by wxiongid
[考研] 材料学硕333求调剂 +7 北道巷 2026-03-24 7/350 2026-03-27 00:30 by wxiongid
[考研] 化学工程085602 305分求调剂 +17 RichLi_ 2026-03-25 17/850 2026-03-26 19:44 by plmuchong
[考研] 275求调剂 +9 Micky11223 2026-03-25 11/550 2026-03-26 15:54 by 不吃魚的貓
[考研] 一志愿陕师大生物学071000,298分,求调剂 +4 SYA! 2026-03-23 4/200 2026-03-26 15:27 by caiaijun80
[考研] 332求调剂 +6 032500 2026-03-25 6/300 2026-03-25 22:45 by 418490947
[考研] 材料与化工328分调剂 +6 。,。,。,。i 2026-03-23 6/300 2026-03-25 22:30 by 418490947
[考研] 299求调剂 +7 shxchem 2026-03-20 9/450 2026-03-25 10:41 by lbsjt
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 一志愿华东理工大学081700,初试分数271 +5 kotoko_ik 2026-03-23 6/300 2026-03-24 10:29 by 学术搬砖er
[考研] 一志愿吉大化学322求调剂 +4 17501029541 2026-03-23 6/300 2026-03-24 10:21 by 戴围脖的小蚊子
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[考研] 335求调剂 +4 yuyu宇 2026-03-23 5/250 2026-03-23 23:49 by Txy@872106
[考研] 284求调剂 +3 yanzhixue111 2026-03-23 6/300 2026-03-23 22:58 by pswait
[考研] 361求调剂 +3 Glack 2026-03-22 3/150 2026-03-23 22:03 by fuyu_
[考研] 333求调剂 +3 ALULU4408 2026-03-23 3/150 2026-03-23 19:04 by macy2011
[考研] 接收2026硕士调剂(学硕+专硕) +4 allen-yin 2026-03-23 6/300 2026-03-23 15:04 by 汪!?!
[考研] 材料与化工085600,总分304,本科有两篇sci参与,求调剂 +4 幸运的酱酱 2026-03-22 5/250 2026-03-22 20:15 by edmund7
[考研] 0805材料320求调剂 +3 深海物语 2026-03-20 3/150 2026-03-21 15:46 by 无际的草原
信息提示
请填处理意见