| 查看: 2999 | 回复: 3 | |||||||
minmin_0082003金虫 (正式写手)
莫愁莫愁
|
[交流]
【求助】siesta 计算DOS应该注意什么问题 已有3人参与
|
|
输入文件如何设置参数， DOS 与PDOS还有TDOS的主要区别是什么 做DOS是为了解释说明orbital之外，还有其他用途吗？ 请大家指教 |
回复此楼» 收录本帖的淘帖专辑推荐
 第一性原理相关文档 | siesta问题 | 量化计算有意义的讨论 | siesta & vasp |
» 猜你喜欢
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有6人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有7人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有4人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
5个%2F，非固定段22个字母
已经有12人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
-
售T0P一区SCI文章，我:8O5.51.O.54,科目齐全,可＋急
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 看看别人仅用VASP做DOS计算发好文章（两篇PRL）
已经有30人回复
- siesta 计算电子占据数的问题
已经有3人回复
- 【求助】理论计算新手，求帮忙解决带隙及DOS的问题
已经有21人回复
- 【求助】siesta进行ordern计算出错求助
已经有4人回复
- 【求助】DOS和band带隙不一致，是什么问题
已经有20人回复
- [求助]什么结构可以在VNL 中计算能带和DOS
已经有16人回复
- 【求助】DOS计算问题
已经有12人回复
- 【求助】siesta中transiesta模块的收敛问题
已经有13人回复
- 【求助】siesta中关于几何结构优化的问题【完结】
已经有17人回复
- 【求助】siesta怎样计算态密度
已经有8人回复
dalongmao
木虫 (正式写手)
- 应助: 3 (幼儿园)
- 金币: 4358.4
- 散金: 6
- 红花: 2
- 帖子: 613
- 在线: 448.6小时
- 虫号: 596465
- 注册: 2008-09-09
- 专业: 理论和计算化学
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
sunyang1988(金币+2):谢谢交流 2010-11-10 18:52:58
小木虫(金币+0.5):给个红包，谢谢回帖交流
sunyang1988(金币+2):谢谢交流 2010-11-10 18:52:58
|
The block must be a single line with the energies of the range for PDOS projection, (relative to the program's zero, i.e. the same as the eigenvalues printed by the program), the peak width (an energy) for broadening the eigen- values, the number of points in the energy window, and the energy units. %block ProjectedDensityOfStates -20.00 10.00 0.200 500 eV %endblock ProjectedDensityOfStates 手册中说了设置的问题 |
2楼2010-11-10 15:54:04
fengshiquan
铁杆木虫 (正式写手)
- 应助: 10 (幼儿园)
- 金币: 5831.8
- 散金: 80
- 红花: 5
- 帖子: 599
- 在线: 1569.4小时
- 虫号: 884725
- 注册: 2009-10-26
- 专业: 原子和分子物理
3楼2011-03-30 11:46:49
hw45888792
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 5 (幼儿园)
- 金币: 848.5
- 散金: 45
- 红花: 3
- 帖子: 339
- 在线: 282.9小时
- 虫号: 711113
- 注册: 2009-02-28
- 性别: GG
- 专业: 理论和计算化学
★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 首先谢谢你的回复，如果是转载，希望你能给出链接，如果是原创，请注明，我们会给与更多的奖励，谢谢合作 2011-03-31 17:26:16
小木虫(金币+0.5):给个红包，谢谢回帖交流
zzy870720z(金币+2): 首先谢谢你的回复，如果是转载，希望你能给出链接，如果是原创，请注明，我们会给与更多的奖励，谢谢合作 2011-03-31 17:26:16
|
siesta具体设置应该都知道，看看怎么分析，以及有什么意义： MS电荷密度图、能带结构、态密度的分析 如何分析第一原理的计算结果 用第一原理计算软件开展的工作，分析结果主要是从以下三个方面进行定性/定量的讨论： 1、电荷密度图（charge density）； 2、能带结构（Energy Band Structure）； 3、态密度（Density of States，简称DOS）。 电荷密度图是以图的形式出现在文章中，非常直观，因此对于一般的入门级研究人员来讲 不会有任何的疑问。唯一需要注意的就是这种分析的种种衍生形式，比如差分电荷密图（deformation charge density）和二次差分图（difference charge density）等等，加自旋极化的工作还可能有自旋极化电荷密度图（spin-polarized charge density）。所谓“差分” 是指原子组成体系（团簇）之后电荷的重新分布，“二次”是指同一个体系化学成分或者几何构型改变之后电荷的重新分布，因此通过这种差分图可以很直观地看出体系中个原子的成键情况。通过电荷聚集（accumulation）/损失（depletion）的具体空间分布，看成键的极性强弱；通过某格点附近的电荷分布形状判断成键的轨道（这个主要是对d轨道的分析，对于s或者p轨道的形状分析我还没有见过）。分析总电荷密度图的方法类似，不过相对而言，这种图所携带的信息量较小。 能带结构分析现在各个领域的第一原理计算工作中用得非常普遍了。但是因为能带这个概念本身的抽象性，对于能带的分析是让初学者最感头痛的地方。关于能带理论本身，我在这篇文章中不想涉及，这里只考虑已得到的能带，如何能从里面看出有用的信息。首先当然可以看出这个体系是金属、半导体还是绝缘体。判断的标准是看费米能级和导带（也即在高对称点附近近似成开口向上的抛物线形状的能带）是否相交，若相交，则为金属，否则为半导体或者绝缘体。对于本征半导体，还可以看出是直接能隙还是间接能隙：如果导带的最低点和价带的最高点在同一个k点处，则为直接能隙，否则为间接能隙。在具体工作中，情况要复杂得多，而且各种领域中感兴趣的方面彼此相差很大，分析不可能像上述分析一样直观和普适。不过仍然可以总结出一些经验性的规律来。主要有以下几点： 1） 因为目前的计算大多采用超单胞（supercell）的形式，在一个单胞里有几十个原子 以及上百个电子，所以得到的能带图往往在远低于费米能级处非常平坦，也非常密集。原则上讲，这个区域的能带并不具备多大的解说/阅读价值。因此，不要被这种现象吓住，一般的工作中，我们主要关心的还是费米能级附近的能带形状。 2） 能带的宽窄在能带的分析中占据很重要的位置。能带越宽，也即在能带图中的起伏越大，说明处于这个带中的电子有效质量越小、非局域（non-local）的程度越大、组成这条能带的原子轨道扩展性越强。如果形状近似于抛物线形状，一般而言会被冠以类sp带（sp-like band）之名。反之，一条比较窄的能带表明对应于这条能带的本征态主要是由局域于某个格点的原子轨道组成，这条带上的电子局域性非常强，有效质量相对较大。 3） 如果体系为掺杂的非本征半导体，注意与本征半导体的能带结构图进行对比，一般而言在能隙处会出现一条新的、比较窄的能带。这就是通常所谓的杂质态(doping state)，或者按照掺杂半导体的类型称为受主态或者施主态。 4） 关于自旋极化的能带，一般是画出两幅图：majority spin和minority spin。经典的说，分别代表自旋向上和自旋向下的轨道所组成的能带结构。注意它们在费米能级处的差异。如果费米能级与majority spin的能带图相交而处于minority spin的能隙中，则此体系具有明显的自旋极化现象，而该体系也可称之为半金属（half metal）。因为majority spin与费米能级相交的能带主要由杂质原子轨道组成，所以也可以此为出发点讨论杂质的磁性特征。 5） 做界面问题时，衬底材料的能带图显得非常重要，各高对称点之间有可能出现不同的情况。具体地说，在某两点之间，费米能级与能带相交；而在另外的k的区间上，费米能级正好处在导带和价带之间。这样，衬底材料就呈现出各项异性：对于前者，呈现金属性，而对于后者，呈现绝缘性。因此，有的工作是通过某种材料的能带图而选择不同的面作为生长面。具体的分析应该结合试验结果给出。（如果我没记错的话，物理所薛其坤研究员曾经分析过$\beta$-Fe的(100)和(111)面对应的能带。有兴趣的读者可进一步查阅资料。） 原则上讲，态密度可以作为能带结构的一个可视化结果。很多分析和能带的分析结果可以一一对应，很多术语也和能带分析相通。但是因为它更直观，因此在结果讨论中用得比能带分析更广泛一些。简要总结分析要点如下： 1） 在整个能量区间之内分布较为平均、没有局域尖峰的DOS，对应的是类sp带，表明电子的非局域化性质很强。相反，对于一般的过渡金属而言，d轨道的DOS一般是一个很大的尖峰，说明d电子相对比较局域，相应的能带也比较窄。 2） 从DOS图也可分析能隙特性：若费米能级处于DOS值为零的区间中，说明该体系是半导体或绝缘体；若有分波DOS跨过费米能级，则该体系是金属。此外，可以画出分波（PDOS）和局域（LDOS）两种态密度，更加细致的研究在各点处的分波成键情况。 3） 从DOS图中还可引入“赝能隙”（pseudo gap）的概念。也即在费米能级两侧分别有两个尖峰。而两个尖峰之间的DOS并不为零。赝能隙直接反映了该体系成键的共价性的强弱：越宽，说明共价性越强。如果分析的是局域态密度（LDOS），那么赝能隙反映的则是相邻两个原子成键的强弱：赝能隙越宽，说明两个原子成键越强。上述分析的理论基础可从紧束缚理论出发得到解释：实际上，可以认为赝能隙的宽度直接和Hamiltonian矩阵的非对角元相关，彼此间成单调递增的函数关系。 4） 对于自旋极化的体系，与能带分析类似，也应该将majority spin和minority spin分别画出，若费米能级与majority的DOS相交而处于minority的DOS的能隙之中，可以说明该体系的自旋极化。 5） 考虑LDOS，如果相邻原子的LDOS在同一个能量上同时出现了尖峰，则我们将其称之为杂化峰（hybridized peak），这个概念直观地向我们展示了相邻原子之间的作用强弱。 以上是本人基于文献调研所总结的一些关于第一原理工作的结果分析要点。期冀能对刚进入这个领域内的科研工作者有所启发。受本人的水平所限，文章的内容可能会有理论上的不足 甚至错误之处，希望大家指出，共同发展第一原理计算物理的方法和研究内容。 能带结构（翻译） 点击数：743 发布日期：2007-1-27 20:19:00 【收藏】 【评论】 【打印】 【编程爱好者论坛】 【关闭】 Tag:能带 布里渊Brillouin k空间 晶体 周期性 电子结构 平面波 个人原创 在分子中可能的电子能级是分立的、量子化的。但分子变得更大时，这些能级相互就会靠得更近。在晶体里能级之间靠得非常近以致于形成了连续的带子，这些带子的能量具有实际的利用目的。因此，晶体的电子结构可以用其能带结构来描述。 能带的数学描述无限晶体的电子结构用能带图来描述，能带图给出k空间——叫作布里渊（Brillouin）区——中各点的电子轨道的能量。这与角分辨光电子能谱实验结果相一致。 k空间不是一个物理空间，它是对轨道成键性质的一种描述。一个无限长的原子链中，轨道?相位可以是从全成键到全反键（这两个极端情况分别记为k=0和k=π/a）之间的任何状态。其中有时是一条直线有三个成键原子再接着一个反? 的原子的结合方式或者其他什么结合方式。定义了 k空间后，对于某些原子k=0对应于全成键的对称性，而对于其他原子则是全反键对称的，这取决于原子轨道的对称性。 对于三维晶体k空间是三维的，（kx，ky，kz）。k空间中的某些点具有特定的名称。在各维空间中，符号“Γ”指的都是k=0的点，“Μ”指的 都是k=π/a的点。“Χ”、“Y”、“Κ”和“Α”指的是k=0在某些方向上以及k=π/a在其他方向上的点，这取决于晶体的对称性。典型的能带结构图——称为spaghetti图——画出了沿着这些k点所对应的轨道能量，见图34.1。这些符号在参考文献中有更相详细地讨论。 由于轨道展开成了能带，用于形成σ键或σ反键的轨道就展开成更宽的能带，π轨道则形成更窄的能带，而δ轨道则形成最窄的轨道。 计算带隙 有时候研究者只需要知道晶体的带隙。一旦一条完整的能带计算出来，通过观察自然就很容易知道带隙了。但是计算全部能带可能会花费大量的工作，得到许多不必要的信息。估算带隙有一些方法，但并不完全可靠。只在布里渊区的Μ、Κ、Χ和Γ点进行能带结构计算还不足以形成一条能带，因为任何给定的能带的能量极小点和极大点有时会落在这些k点之间。如果计算方法需要较高级别的CPU计算，有时就会进行这样的有限计算。例如，在确定? 否有必要进行高级别的完全计算时，就有可能先进 行这种选点的高级别计算。 有些研究者用分子的计算结果来估计从HOMO到LUMO的带隙。当分子变得更大时，这种带隙会 变得更小，因此就有可能对一些按大小递增的分子 进行量子力学计算，然后通过外推预测无限体系的带隙，这对于通常不是晶体的聚合物很有 用。这些体系也用到一维能带结构，因此必须假定 它们是晶体或者至少是高度的有序的。 计算能带结构 从头算和半经验计算可以得出能量，因而可以用来计算能带结构。但是如果计算一个分子的 能量耗时较长，那么计算布里渊区的一系列点则耗 时更长，要是不需要太精确的结果，可以选用扩展休克尔方法来计算。在能带计算中扩展休 克尔方法有时叫作紧束缚近似。近年来更倾向于使 用从头算或密度泛函（DFT）方法。 就象分子计算那样，从头算需要用基组和一定的方法来计算能量，但计算能带时基组的选择 与计算分子时有些不同。拥有弥散函数的大基组在 相邻的晶胞之间由于存在较大的重叠而发生收缩，这会造成线性相关性，使得方程不能自洽 求解，为此常常用中小基组来解决上述问题。用于 分子计算的原子轨道线性组合（LCAO）方案也可用于晶体的计算，但这并不是唯一的选择。 事实上，以原子为中心的基函数组成布洛赫（Bloch ）函数，布洛赫（Bloch）函数满足体系的平移对称性，但仍然使用LCAO的叫法。 其他有关基组的流行方法时平面波函数方法。之所以提出平面波是因为平面波反映了晶体的 无限平移对称性。最早的平面波计算假定薛定谔方 程在每个原子的附近区域是球对称的（松饼罐头势），但却无法保证电荷守恒。对于离子晶 体松饼罐头计算能给出合理结果，但随着计算技术 和硬件的发展，使人们可以进行更加精确可靠的计算，也就不再采用松饼罐头方法了。还在 使用的一种方法是扩展平面波（APW）方法，是在 Vigner-Seitz晶胞上的晶胞计算。某些类型的问题有许多其他基函数方法。 非常复杂的体系都已经进行了能带结构的计算，然而大多数软件都不够自动化或不够快，不 足以用于临时进行能带计算。计算能带的程序的输 入比大多数计算分子的程序要复杂得多。分子几何构型的输入采用分数坐标，还必须提供原 胞格子矢量和晶体学角度，还可能有必要提供k点的 列表及其简并度。检查各个输入中控制收敛的选项对于计算精度的影响是最保险的措施，软 件附带的手册可能会给出一些推荐值。研究者要想 完成能带计算应当投入大量时间，尤其在学习使用软件阶段。 正如上面所提到的，随着时间推移人们倾向的模拟晶体的计算方法是不断变化的。下面是基 函数方法的列表，按照出现的先后顺序排列： 1. 原子轨道线性组合方法（LCAO） 2. 扩展平面波方法（APW） 3. Korringa、Kohn和Rostoker的格林（Green）函数方法（有时叫作KKR方法） 4. 正交平面波方法（OPW） 5. 赝势方法 6. 各种近似或经验方法 任何基于轨道的方法都可用来计算晶体结构，而趋势是向着更加精确的方法。一些APW和格 函数方法使用了经验参数，因而将它们划到半经验 方法中去。按照使用偏好的顺序，最常用的方法是： 1. 半自洽从头算方法或DFT方法 2. 半经验方法 3. 使用专门的或模拟的势能的方法 描述晶体的电子结构 分子计算中的布居分析方法不能直接应用于能带计算，分析能带结构引入了一系列的方法， 这些方法一般都表示成图线，图线上的数据源于对k 空间中各个点的计算结果。计算大量的点可以得到很好的图线，但为了节省计算时间可以加 大取点间隔，然后用内插法平滑曲线。通常谨慎的 做法是逐次加大取点紧密程度计算几次，看看图线是否有显著变化。 一个重要的问题是，一个给定能级有多少可能的轨道。这可以用态密度图（DOS）来表示，? 图34.2。图中往往用虚线来表示费米（Fermi）能 级。具有半满能带的材料是导体，但如果它们只有少量的未充满的轨道，就可能是不良导体 。有时特别轨道对DOS的贡献会在同一张图上用阴影 区域或虚线画出。 另一个问题是被充满的轨道是成键性的还是反键性的。这可用晶体轨道重叠分布图（COOP） 来表示，如图34.3。一般正的成键区域画在零值线 的右边。 费米能级是填充轨道的最高能级，类似于HOMO能级。如果轨道是半充满的，其能级就会出现 在k空间的点的集合上，称为费米面。 计算晶体性质 晶体计算方面的进展没有分子计算方面的多。经常计算的一个性质是体积弹性模量，它反映 了材料的强度。 在预测热力学条件下会形成什么产物时，可能需要预测哪种晶体结构最稳定，这是一项艰巨 的任务。到目前为止，还没有提出一个完全自动的 的方法试遍由特定的元素集合组成的所有可能的晶体结构。即便这种尝试可以实现，进行计 算所需要消耗的电能也是巨大的。这样的研究经常 用于测试一系列相似的结构，结果无论如何总是正确的。能量最小化也会用到，但须保证起 始结构具有正确的对称性。 缺陷计算 有时并不只对无限体系感兴趣，更关心于晶体中的异类物质，比如晶体吸收的额外的原子。 这时晶体的无限平移对称性并不严格正确。最广泛 应用的模拟方法是Mott-Littleton缺陷方法，这是用来进行晶格局部区域能量最小化的一种 方法。这种方法包含了对晶体中其余部分所受的极 化的连续性描述。 参考文献 描述能带结构的教科书有： P. Y. Yu, M. Cardone, Fundamentals of Semicondutors Springer-Verlag, Berlin(1996 ). R. Ho®mann, Solids and Surfaces A Chemist's View of Bonding in Extended Stru ctures VCH, New York (1988). I. M. Tsidilkovski, Band Structure of Semiconductors Pergamon, Oxford (1982). W. A. Harrison, Solid State Theory Dover, New York (1979). A. A. Levin, Solid State Quantum Chemistry McGraw-Hill, New York (1977). Orbital Theories of Molecules and Solids N. H. March, Ed., Clarendon, Oxford (19 74). J. C. Slater, Quantum Theory of Molecules and Solids McGraw-Hill, New york (1963 ). 综述性文章有： J. K. Burdett, J. Phys. Chem. 100, 13274 (1996). C. R. A. Catlow, J. D. Gale, R. W. Grimes, J. Solid State Chem. 106, 13 (1993). E. Canadell, M.-H. Whangbo, Chem. Rev. 91, 965 (1991). J. H. Harding, Rep. Prog. Phys. 53, 1403 (1990). C. R. A. Catlow, G. D. Price, Nature 347, 243 (1990). JCS Faraday Trans. II Vol. 85, part 5 (1989). R. Ho®mann, Angew. Chem. Int. Ed. Engl. 26, 846 (1987). C. R. A. Catlow, Ann. Rev. Mater. Sci. 16, 517 (1986). J. C. Slater, Adv. Quantum Chem. 1, 35 (1964). 用于作图的程序是YAeHMOP，有两个可执行文件，分别是： G. A. Landrum, W. V. Glassey, Bind 3.0 (1998). G. A. Landrum, Viewkel 2.0 (1998). 都在http://overlap.chem.cornell.edu:8080/yaehmop.html上。 以上文字译自 Computational Chemistry: A Practical Guide for Applying Techniques to Real-World Problems. David C. Young Copyright ( 2001 John Wiley & Sons, Inc. ISBNs: 0-471-33368-9 (Hardback); 0-471-22065-5 (Electronic) -- 哀公问曰：“何为则民服？” 孔子对曰：“举直错诸枉，则民服；举枉错诸直，则民不服。” 摘要：本文总结了对于第一原理计算工作的结果分析的三个重要方面，以及各自的若干要点 关键词：first principles,态密度,能带结构 用第一原理计算软件开展的工作，分析结果主要是从以下三个方面进行定性/定量的讨论： 1、电荷密度图（charge density）； 2、能带结构（Energy Band Structure）； 3、态密度（Density of States，简称DOS）。 电荷密度图是以图的形式出现在文章中，非常直观，因此对于一般的入门级研究人员来讲 不会有任何的疑问。唯一需要注意的就是这种分析的种种衍生形式，比如差分电荷密图（def- ormation charge density）和二次差分图（difference charge density）等等，加自旋极 化的工作还可能有自旋极化电荷密度图（spin-polarized charge density）。所谓“差分” 是指原子组成体系（团簇）之后电荷的重新分布，“二次”是指同一个体系化学成分或者几何 构型改变之后电荷的重新分布，因此通过这种差分图可以很直观地看出体系中个原子的成键情 况。通过电荷聚集（accumulation）/损失（depletion）的具体空间分布，看成键的极性强 弱；通过某格点附近的电荷分布形状判断成键的轨道（这个主要是对d轨道的分析，对于s或者 p轨道的形状分析我还没有见过）。分析总电荷密度图的方法类似，不过相对而言，这种图所 携带的信息量较小。 能带结构分析现在在各个领域的第一原理计算工作中用得非常普遍了。但是因为能带这个 概念本身的抽象性，对于能带的分析是让初学者最感头痛的地方。关于能带理论本身，我在这 篇文章中不想涉及，这里只考虑已得到的能带，如何能从里面看出有用的信息。首先当然可以 看出这个体系是金属、半导体还是绝缘体。判断的标准是看费米能级和导带（也即在高对称点 附近近似成开口向上的抛物线形状的能带）是否相交，若相交，则为金属，否则为半导体或者 绝缘体。对于本征半导体，还可以看出是直接能隙还是间接能隙：如果导带的最低点和价带的 最高点在同一个k点处，则为直接能隙，否则为间接能隙。在具体工作中，情况要复杂得多， 而且各种领域中感兴趣的方面彼此相差很大，分析不可能像上述分析一样直观和普适。不过仍 然可以总结出一些经验性的规律来。主要有以下几点： 1） 因为目前的计算大多采用超单胞（supercell）的形式，在一个单胞里有几十个原子 以及上百个电子，所以得到的能带图往往在远低于费米能级处非常平坦，也非常密 集。原则上讲，这个区域的能带并不具备多大的解说/阅读价值。因此，不要被这种 现象吓住，一般的工作中，我们主要关心的还是费米能级附近的能带形状。 2） 能带的宽窄在能带的分析中占据很重要的位置。能带越宽，也即在能带图中的起伏越 大，说明处于这个带中的电子有效质量越小、非局域（non-local）的程度越大、组 成这条能带的原子轨道扩展性越强。如果形状近似于抛物线形状，一般而言会被冠以 类sp带（sp-like band）之名。反之，一条比较窄的能带表明对应于这条能带的本征 态主要是由局域于某个格点的原子轨道组成，这条带上的电子局域性非常强，有效质 量相对较大。 3） 如果体系为掺杂的非本征半导体，注意与本征半导体的能带结构图进行对比，一般而 言在能隙处会出现一条新的、比较窄的能带。这就是通常所谓的杂质态(doping state)，或者按照掺杂半导体的类型称为受主态或者施主态。 4） 关于自旋极化的能带，一般是画出两幅图：majority spin和minority spin。经典的 说，分别代表自旋向上和自旋向下的轨道所组成的能带结构。注意它们在费米能级处 的差异。如果费米能级与majority spin的能带图相交而处于minority spin的能隙 中，则此体系具有明显的自旋极化现象，而该体系也可称之为半金属（half metal）。因为majority spin与费米能级相交的能带主要由杂质原子轨道组成，所以 也可以此为出发点讨论杂质的磁性特征。 5） 做界面问题时，衬底材料的能带图显得非常重要，各高对称点之间有可能出现不同的 情况。具体地说，在某两点之间，费米能级与能带相交；而在另外的k的区间上，费 米能级正好处在导带和价带之间。这样，衬底材料就呈现出各项异性：对于前者，呈 现金属性，而对于后者，呈现绝缘性。因此，有的工作是通过某种材料的能带图而选 择不同的面作为生长面。具体的分析应该结合试验结果给出。（如果我没记错的话， 物理所薛其坤研究员曾经分析过$\beta$-Fe的(100)和(111)面对应的能带。有兴趣的 读者可进一步查阅资料。） 原则上讲，态密度可以作为能带结构的一个可视化结果。很多分析和能带的分析结果可以 一一对应，很多术语也和能带分析相通。但是因为它更直观，因此在结果讨论中用得比能带分 析更广泛一些。简要总结分析要点如下： 1） 在整个能量区间之内分布较为平均、没有局域尖峰的DOS，对应的是类sp带，表明电 子的非局域化性质很强。相反，对于一般的过渡金属而言，d轨道的DOS一般是一个很 大的尖峰，说明d电子相对比较局域，相应的能带也比较窄。 2） 从DOS图也可分析能隙特性：若费米能级处于DOS值为零的区间中，说明该体系是半导 体或绝缘体；若有分波DOS跨过费米能级，则该体系是金属。此外，可以画出分波 （PDOS）和局域（LDOS）两种态密度，更加细致的研究在各点处的分波成键情况。 3） 从DOS图中还可引入“赝能隙”（pseudogap）的概念。也即在费米能级两侧分别有两 个尖峰。而两个尖峰之间的DOS并不为零。赝能隙直接反映了该体系成键的共价性的 强弱：越宽，说明共价性越强。如果分析的是局域态密度（LDOS），那么赝能隙反映 的则是相邻两个原子成键的强弱：赝能隙越宽，说明两个原子成键越强。上述分析的 理论基础可从紧束缚理论出发得到解释：实际上，可以认为赝能隙的宽度直接和 Hamiltonian矩阵的非对角元相关，彼此间成单调递增的函数关系。 4） 对于自旋极化的体系，与能带分析类似，也应该将majority spin和minority spin分 别画出，若费米能级与majority的DOS相交而处于minority的DOS的能隙之中，可以说 明该体系的自旋极化。 5） 考虑LDOS，如果相邻原子的LDOS在同一个能量上同时出现了尖峰，则我们将其称之为 杂化峰（hybridized peak），这个概念直观地向我们展示了相邻原子之间的作用强 弱。 以上是本人基于文献调研所总结的一些关于第一原理工作的结果分析要点。期冀能对刚进 入这个领域内的科研工作者有所启发。受本人的水平所限，文章的内容可能会有理论上的不足 甚至错误之处，希望大家指出，共同发展第一原理计算物理的方法和研究内容。 能带结构特征分析 [ 2007-3-26 20:52:00 | By: solomon ] ----------------------------------------------------------------------- 0 推荐 对能带结构的分析，除了可以直观的看出材料的导电性能之外，还可以通过以下特征分析： 1. 能带结构图是能量（E）与K矢量的关系图。K点在第一布里渊区内变化，即K∈（0,π/a）.在高对称点，K值对应取相应固定的值。 对于二维能带结构，G＝（0,0）；X＝（π/a，0）；M＝（π/a，π/a）； 对于三维能带结构，如果第一布里渊区为立方体，则G＝（0,0,0）； M＝（π/a，π/a，0）；X＝（π/a，0,0）；R＝（π/a，π/a，π/a） 但为什么有时不是从原点(G 或 τ)开始呢？ 2. 能带数代表晶体中轨道数，但要从能带数推算出轨道数则需要用态密度的概念与算法了。具体等学习态密度时再总结。 3. 分析能带，要找出能带所对应的轨道。首先找出单胞轨道能级图，将各原子的轨道能级标出来，然后考虑成键，具体是如何考虑呢？没看明白。 4. 能带的发散(带宽)说明在该轨道的电子参加成键的能力，发散(带宽)越大，说明处于这个带中的电子有效质量越小、非局域（non-local）的程度越大、组成这条能带的原子轨道扩展性越强，电子参加成键的能力越强。 参考资料： 1. R.霍夫曼，固体与表面 2. carlon (carlon@hit.edu.cn)，能带结构和态密度图的绘制及初步分析(网文) 3. Satchel1979，如何分析第一原理的计算结果，网文(http://www.quantumchemistry.net/ ... hod/200512/82.html) |
4楼2011-03-31 16:41:09
10