| 查看: 10915 | 回复: 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 下载信息：
 这个是用CNEB方法算的CO在Au@CeO2(111) 上吸附并被CeO2表面O氧化的反应路径。前面计算方法中交代了用20images. 图上列出的IS和FS，以及中间的Splillover, MS,Oxidation几个阶段都太有故事了吧。仅仅通过线性插值能得出这样的过程吗？ 这篇文章中用PBE，别的文献也有用PW91的。对于含f电子的稀土金属氧化物和贵金属体系，选择赝势有什么考虑呢？
同时问问，为什么这篇文章可以发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。为什么可以这样选取呢？
 请仔细阅读图下面的说明文字。先列出两个相关联的问题： （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朵）
awmc2008» 猜你喜欢
0703化学求调剂，各位老师看看我！！！
已经有3人回复
-
339求调剂
已经有5人回复
-
287求调剂
已经有8人回复
-
一志愿北京化工大学 070300 学硕 336分 求调剂
已经有8人回复
-
085602 化工专硕 338分 求调剂
已经有5人回复
-
0856，材料与化工321分求调剂
已经有4人回复
-
一志愿北京化工大学材料与化工 264分各科过A区国家线
已经有10人回复
-
266求调剂
已经有8人回复
-
265求调剂
已经有3人回复
-
286求调剂
已经有6人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 想合成一物质，文献的方法为什么这么复杂？有什么玄机？
已经有6人回复
- 关于CO2的插入反应
已经有7人回复
- DFT计算过程中断，如何接着算
已经有3人回复
- neb优化结束 计算频率 INCAR应该怎么改？
已经有10人回复
- NEB计算过渡态和反应路径
已经有9人回复
- 博士第一篇SCI，给延期的我打了强心针
已经有677人回复
- 等效位置的态密度分析
已经有9人回复
- 哪本书对您学习掌握DFT帮助最大？
已经有38人回复
- 文献回顾分析与META有什么异同？
已经有0人回复
- DFT计算 电容
已经有7人回复
- vasp可以固定单个原子的磁矩吗？
已经有5人回复
- 【求助】总态密度与分波态密度
已经有9人回复
- 【求助】两种计算电荷密度的方法得到的结果相差很大
已经有9人回复
- 【求助】物理吸附的态密度图
已经有5人回复
- 【求助】MS中局域态密度的计算
已经有3人回复
- 【求助】反应路径能量图怎么画？
已经有7人回复
- 【求助】表面优化分析
已经有9人回复
- 【求助】用DFT如何计算磁性物质的顺磁态
已经有3人回复
- 【求助】金属掺杂态密度分析
已经有12人回复
- 【求助成功】CNEB后,还需要做结构优化?
已经有5人回复
- 【求助】如何算adsorption energy ?
已经有5人回复
- 【求助】DFT 计算气体
已经有4人回复
- 【求助】如何用origin做pdos图!
已经有9人回复
- 【求助】计算出来的带隙值比文献上计算出来的小
已经有9人回复
- 【求助】VASP算布局数分析
已经有16人回复
- 【讨论】DFT方法属不属于从头算？
已经有13人回复
- 【资源】一篇介绍gaussian计算ECD谱原理的的文献
已经有7人回复
- 【求助】vasp的NEB计算,输出需要哪些参数?
已经有7人回复
» 抢金币啦！回帖就可以得到:
-
基金祈福
+3/1688
-
数一英一硕士调剂（双一流高校）
+1/188
-
老虫子的注册小木虫18周年大礼包
+5/150
-
[b]【材料表征】成分含量丨微观形貌丨分子结构丨材料加工丨预存服务[/b]
+1/72
-
艾里卡特 (Alicat) - 燃料电池测试系统与质量流量和压力仪表
+3/55
-
重庆交通大学院士团队招收材料、能源动力专业硕士研究生
+2/52
-
南京林业大学木质纤维素功能材料国际联合创新中心招收2026级硕士生（调剂）
+2/48
-
沈阳农业大学生物与医药专硕招收调剂生！！！！！！！！！！
+1/37
-
招收2026级环境科学方向（专业代码083001）硕士研究生1-2名
+1/37
-
2026设施园艺作物环境调控与绿色生产团队招收相关专业调剂研究生
+1/34
-
澳大利亚昆士兰大学化工学院招收全奖博士生2名
+1/26
-
南京林业大学-国家级青年人才团队 招2026级博士、硕士（合成化学、植物化学方向）
+1/20
-
中南大学喻桂朋教授团队新增2026年博士研究生
+1/11
-
韩国庆熙大学 Yu Taekyung 教授招2026年博士研究生
+1/11
-
shu湘潭大学化工学院招收化工、能源动力、工程热物理相关研究生调剂
+1/8
-
2026年中国科学院城市环境所与福建农林大学联合培养硕士研究生项目介绍
+1/6
-
湖南农业大学 环境科学与工程0830+资源与环境专硕0857 调剂名额
+1/4
-
湘潭大学教师招聘
+1/1
-
福建师范大学环境与资源学院招收2026级硕士生
+1/1
-
【博士招生】上海交通大学机动学院-招2026年秋“申请-考核制”博士生1名-3D打印方向
+1/1
19楼2011-06-27 22:17:26
★ ★
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
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
3楼2010-12-06 12:44:23
4楼2010-12-06 12:56:19
