| 查看: 11041 | 回复: 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化学招61最终排名62化学求调剂
已经有32人回复
-
化工学硕294分，求导师收留
已经有16人回复
-
考研英一数一338分
已经有7人回复
-
302求调剂
已经有6人回复
-
271求调剂
已经有25人回复
-
266调剂
已经有10人回复
-
生物学调剂
已经有11人回复
-
求调剂
已经有8人回复
-
085801电气专硕272求调剂
已经有10人回复
-
一志愿南科大生物学297分，求调剂推荐
已经有11人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 想合成一物质，文献的方法为什么这么复杂？有什么玄机？
已经有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人回复
» 抢金币啦！回帖就可以得到:
-
国家“千人计划”青年项目获得者团队最后2人07、08等学科研究生调剂
+1/285
-
陕科大环境学院招收0857和0830调剂考生【可跨门类，最后2两小时】
+1/200
-
哈尔滨工程大学核动力装置研究所招收1-2名博士生
+1/88
-
安心播种，静待花开
+1/77
-
五邑大学 化学工程与技术调剂 4月10号九点截止！！！！！
+1/42
-
接收调剂：组内氛围轻松活跃，欢迎加入
+1/40
-
中国科学院超导与储能材料研究组诚聘特别研究助理（博士后）
+1/36
-
招收调剂生（一志愿07开头专业），海洋生物专业4个名额
+1/36
-
西交利物浦大学黄彪院士招收26年全奖博士生1名（工业系统数字孪生方向）
+1/35
-
延安大学化学与化工学院接收化学专业调剂生（工科、理学均可调剂）
+5/30
-
一本院校，化学、环境工程，上车就走
+1/22
-
欢迎调剂到郑州轻工业大学材料与化学工程学院
+1/12
-
陕西科技大学轻工科学与工程学院刘老师接收08专业研究生调剂
+1/8
-
西安工业大学光电学院接收调剂
+1/7
-
烟台大学精准材料高等研究院26年材料方向研究生招生
+1/5
-
江苏海洋大学+环境与化学工程学院+课题组招调剂生+第一批调剂
+2/4
-
河南城建学院生命科学与工程学院（085700资源与环境方向）现有调剂名额，欢迎调剂咨询
+1/2
-
湖南理工大学 化学、化工调剂指标充足，录取率高
+1/2
-
北大口腔招人工智能、生物再造科研助理（优秀者可留校或推荐读博）
+1/1
-
安徽工程大学化学与环境工程学院2026级硕士研究生还有部分调剂名额-资源与环境专硕
+1/1
9楼2010-12-07 09:19:12
★ ★
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
