Implementation Strategies for Orbital-dependent Density Functionals

The development of density functional theory (DFT) has been focused primarily on two main pillars: (1) the pursuit of more accurate exchange-correlation (XC) density functionals; (2) the feasibility of computational implementation when dealing with many-body systems. In this context, this work is aimed on using one-dimensional quantum systems as theoretical laboratories to investigate the implementation of orbital functionals (OFs) of density. By definition, OFs are those which depend only implicitly on the density, via an explicit formulation in terms of Kohn-Sham orbitals. Typical examples are the XC functionals arising from the Perdew-Zunger self-interaction correction (PZSIC). Formally, via Kohn-Sham equations, the implementation of OFs must be performed by means of the optimized effective potential method (OEP), which is known by requiring an excessive computational effort even when dealing with few electrons systems. Here, we proceed a systematical investigation aiming to simplify or avoid the OEP procedure, taking as reference the implementation of the PZSIC correction applied to one-dimensional Hubbard chains.     

描述甲烷在Pt(110)-(21)表面的化学吸附解离的密度泛函评估

本文探讨了几种梯度近似(GGA)密度泛函及元梯度近似(meta-GGA)密度泛函在描述甲烷在重构的Pt(110)-(2×1)上的解离化学吸附作用的适用性. 金属的体相和表面结构、甲烷的吸附能量和解离能垒等被用来评估泛函的可靠性. 另外,在从头算分子动力学计算中,采用范德瓦尔斯矫正的GGA函数(optPBE-vdW)和范德瓦尔斯矫正的meta-GGA函数(MS-PBEl-rVV10）计算粘附概率. 计算结果表明,使用这两种泛函能更好地与现有的实验结果吻合,从而为发展甲烷在Pt(110)-(2×1)表面解离的可靠机器学习势能面打下重要基础.     

Prediction of the derivative discontinuity in density functional theory from an electrostatic description of the exchange and correlation potential

We propose an approach to approximate the exchange and correlation (XC) term in density functional theory. The XC potential is considered as an electrostatic potential, generated by a fictitious XC density, which is in turn a functional of the electronic density. We apply the approach to develop a correction scheme that fixes the asymptotic behavior of any approximated XC potential for finite systems. Additionally, the correction procedure gives the value of the XC derivative discontinuity; therefore, it can directly predict the fundamental gap as a ground-state property.     

Localization and delocalization errors in density functional theory and implications for band-gap prediction

The band-gap problem and other systematic failures of approximate exchange-correlation functionals are explained from an analysis of total energy for fractional charges. The deviation from the correct intrinsic linear behavior in finite systems leads to delocalization and localization errors in large and bulk systems. Functionals whose energy is convex for fractional charges such as the local density approximation display an incorrect apparent linearity in the bulk limit, due to the delocalization error. Concave functionals also have an incorrect apparent linearity in the bulk calculation, due to the localization error and imposed symmetry. This resolves an apparent paradox and identifies the physical nature of the error to be addressed to obtain accurate band gaps from density functional theory.     

Hydrogen adsorption and diffusion on the PdTa alloy surface

The hydrogen adsorption on the PdTa alloy surface is studied using a pseudopotential method with a generalized gradient approximation for an exchange-correlation functional. The most preferable hydrogen adsorption sites on two low-index surfaces ((001), (110)) are determined. It is shown that hydrogen adsorption at the bridge site is preferred on the PdTa(001) surface that terminates by one or two tantalum layers and on PdTa(110). The preference of hydrogen adsorption at tantalum-rich sites is caused by partial population of its d shell. During adsorption, the electronic structure of the states involved in interaction with hydrogen is shown to change most substantially, which is accompanied by the corresponding shifts of these states and the appearance of peaks in the densities of states of the metal in the region of the hydrogen valence band. The effect of hydrogen on the electron and structural characteristics of the surfaces is analyzed. The hydrogen diffusion barriers are calculated in the bulk of the alloy and from the surface into the bulk.     

交换关联泛函近似方法对Ru催化剂上甲烷和乙烷形成的影响

密度泛函理论作为多相催化研究中的一个强有力工具,常被用于获得催化过程中关键的热力学及动力学参数,如吸附能、反应焓、活化能垒和速率常数等. 理解密度泛函交换关联近似方法对于揭示催化剂的催化性能及机理至关重要. 本文报道了六种不同的交换关联泛函近似方法,包括PBE、RPBE、BEEF+vdW、optB86b+vdW、SCAN和SCAN+rVV10,对金属Ru(0001)和Ru(1011)表面上甲烷和乙烷形成过程中涉及到的中间体的吸附能、反应能和活化能垒的影响. 当基元反应中反应物和产物与表面的配位数不同时,理论计算的反应能大小强烈依赖于交换关联密度泛函的选择. 对于涉及多个基元步骤的总反应,反应能的计算偏差会逐渐累积,从而导致不同的交换关联泛函近似方法之间的巨大差异. 由于不同泛函对反应涉及到的中间体吸附能之间存在差异,交换关联泛函近似方法的选择显著地影响Ru(0001)表面上甲烷、乙烯和乙烷的选择性. 然而,不同泛函近似方法对于Ru(0001) 和Ru(1011)表面上基元反应的能垒以及结构敏感性影响不大. 本工作不仅揭示了交换关联密度泛函近似方法在理论计算研究催化领域的局限性,也强调了选择合适的交换关联泛函方法对于正确评估催化剂活性和选择性的重要性.     

Oxidation of Co(poly) and Co(1010) surfaces: Structure and temperature effects

The oxidation of Co(poly) and Co(10$\text{1}$0) surfaces has been investigated using AES within the temperature range 300–700 K. It is found that at both surfaces the initial step — dissociative adsorption of oxygen — exhibits identical kinetics, independent of temperature. Approaching the oxide layer formation regime, temperature and surface structure of the substrate as well determine the oxygen uptake. With the aid of sputter profiling through the oxidized surface layers it is seen that at both, single crystal and polycrystal, oxygen is present even far below the actual surface. The diffusion of oxygen into the bulk is found to be faster at the polycrystalline sample.     

Electron loss spectroscopy of clean and oxygen covered GaAs(110) surfaces

Electron energy loss spectra of clean and oxygen covered GaAs(110) surfaces have been measured with a four grid retarding field analyser. Loss spectra of clean cleaved p- and n-type surfaces are slightly different and different states of adsorption for the oxygen on the two surfaces are found. The loss peaks which are common in the spectra obtained from clean surfaces of both types of material have been interpreted in terms of bulk and surface excitations. The data associated with the bulk excitations are in good agreement with previous optical and electron transmission data while loss peaks at 11.5 and 18.5 eV are interpreted as the surface plasma loss and a surface state transition respectively. For n-type material extra loss peaks were observed. In the case of oxygen adsorption on these surfaces new loss peaks were found at 13.5, 17.2 and 28.1 eV in both spectra and are assumed to be characteristic of the oxygen. Further, for n-type material an extra peak occurs at 8.2 eV.     

Density-functional study of the CO chemisorption on bimetallic Pd–Sn(1 1 0) surfaces

We study the ordered PdSn c(2 × 2), (2 × 1), and PdSn₂ (3 × 1) overlayers deposited on Pd(1 1 0) by using first-principles density-functional calculations. It appears that the two PdSn structures are almost degenerate in the energy. Pd–Sn surfaces we consider do not display the marked buckling with Sn atoms displaced towards vacuum that is common for Pt–Sn surfaces. Low-coverage CO chemisorption at these overlayers and on analogous surface structures on Pd₃Sn is considered. It is shown that inclusion of an empirical correction to the CO adsorption energy changes the stable adsorption site from the long-bridge to the top one in most cases. The adsorption energy decreases with the number of Sn atoms in the vicinity of the adsorption site, and this property correlates well with the position of the centre of gravity of the local Pd d-electron band, and also with the variation of the local density of d-electron states at the Fermi level. The centre-of-gravity value is used to assess the core-level shifts for Pd atoms in various geometries. Most of the calculated data compare rather well with the recent measurements on Pd–Sn overlayers at Pd(1 1 0) as well as with other data on related bimetallic systems.     

InN doped with Zn: Bulk and surface investigation from first principles

Structures and stabilities of Zn adsorption and incorporation at InN surfaces are systematically investigated by first-principles calculations. An InN (0001)–(2×2) surface covered by 3/4 monolayer Zn adsorption atoms at the H₃ sites is found to be energetically favorable. The calculated surface energies demonstrate the stability of Zn-incorporated surfaces. Substitutional defects may act as a potential source for the bulk and surface p-type behavior in Zn-doped InN.     

CO adsorption studies on pure and Ni-covered Cu(111) surfaces

The adsorption of CO on pure and Ni-covered Cu(111) surfaces has been studied by means of LEED, TDS, UPS and work function measurements during adsorption and desorption. Different Ni-coverages between 0.1 and 2 monolayers were obtained by Ni-evaporation controlled by a quartz micro balance and by AES. Near room temperature Ni grows in a layer-by-layer mode on Cu(111). The island structure of the surfaces with submonolayer Ni-coverages is clearly demonstrated by TDS und LEED results obtained after CO adsorption. As with surfaces of bulk Cu-Ni alloys CO adsorption on Cu(111) with submonolayer Ni-coverage is dominated by a site effect. Cu-, Ni-, and mixed adsorption sites can be distinguished. The CO induced work function changes for Ni- and Cu-site adsorption show the same sign as observed with the pure metals. Mixed site adsorption has only a minor influence on the work function. A “ligand effect” observed only for the Ni-site adsorption, and only at small Ni-coverages is discussed in detail. Studies on the adsorption kinetics reveal that the Cu-sites may serve as precursor sites for Ni-site adsorption. Detailed UPS studies demonstrate that the CO-induced emission maxima observed on Cu surfaces with submonolayer Ni-coverages can be interpreted as a superposition of the respective adsorption features observed with the pure metals, roughly separated by their work function difference.     

The distribution of charge below the (001) surfaces of gallium arsenide

Whilst the surfaces of several semiconductors have received extensive study, less attention has been paid to what happens below the surface. We have examined the electronic structure of both arsenic and gallium unrelaxed (001) surfaces, our results demonstrating the rather openess of these surfaces for adsorption by impurities. It is also shown that electronic structure up to at least the fourth or fifth layers of atoms into the bulk is affected, and that regions of large charge density are not centred on lattice positions.     

A first principle study on the magnetic properties of Cu2O surfaces

Spin-polarized density functional theory calculations were performed to investigate the magnetism of bulk and Cu₂O surfaces. It is found that bulk Cu₂O, Cu/O-terminated Cu₂O(111) and (110) surfaces have no magnetic moment, while, the O-terminated Cu₂O(100) and polar O-terminated Cu₂O(111) surfaces have magnetism. For low index surfaces with cation and anion vacancy, we only found that the Cu vacancy on the Cu₂O(110) Cu/O-terminated surface can induce magnetism. For atomic and molecular oxygen adsorption on the low index surfaces, we found that atomic and molecular oxygen adsorption on the Cu-terminated Cu₂O(110) surface is much stronger than on the Cu/O-terminated Cu₂O(111) and Cu-terminated Cu₂O(100) surfaces. More interesting, O and O₂ adsorption on the surface of Cu/O terminated Cu₂O(111) and O₂ adsorption on the Cu-terminated Cu₂O(110) surface can induce weak ferromagnetism. In addition, we analysis origin of Cu₂O surfaces with magnetism from density of state, the surface ferromagnetism possibly due to the increased 2p–3d hybridization of surface Cu and O atoms. This is radically different from other systems previously known to exhibit point defect ferromagnetism, warranting a closer look at the phenomenon.     

Quantitative Auger electron spectroscopy analysis of binary alloys with separate matrix correction factors for bulk and surface

A new method which takes into account the separate matrix correction factors for bulk and surface has been tried out for quantitative Auger electron spectroscopy analysis of binary alloys. The calculations use an iteration scheme. It has been applied to the Fe-Cr alloys studied at this Centre and the Auger electron spectroscopy data for the other alloy systems available in literature. The results are now more compatible with the expectation that the surface composition is different from the bulk.     

氢在担载金属表面的吸附研究

应用格林函数方法在紧束缚近似下研究了氢在担载金属表面的吸附性质。采用自洽的Anderson-Newns吸附模型,对氢在Pt/ZnO,Cu/ZnO和Ni/ZnO三种担载式复合体系表面的吸附能△E、吸附态能级E_ad作了计算,并讨论了金属簿层在ZnO衬底上的沉积厚度及金属-衬底相互作用对氢在该类复合体系表面的吸附性质的影响。计算表明,金属-衬底相互作用越强,氢在Pt(Cu,Ni)/ZnO体系表面的吸附能及电荷转移量越小。金属-衬底相互作用抑制了氢在金属表面的吸附。衬底对金属表面吸附性质的影 关键词：     

Interaction of oxygen and CO with thin Ag films at AgTi2 substrates

Clean Ag overlayers with a thickness of about three atomic layers are prepared at polycrystalline AgTi₂ surfaces. Adsorption of oxygen and CO at these overlayers is studied with AES, UPS, and TDS methods. UP spectra reveal that oxygen adsorbs dissociatively on the surface of the Ag overlayers at 120 K, similar as at the surfaces of bulk Ag crystals. At higher temperatures, O atoms diffuse into the overlayer and get trapped at Ti at the interface underneath the Ag layer. Unlike at bulk Ag, oxygen adsorption proceeds with high sticking probability, which is attributed to a reduced activation barrier. CO adsorbs at the Ag overlayers in a weakly chemisorbed state with an adsorption energy of ca. 9 kcal/mol, whereas at bulk Ag true physisorption has been observed. Accordingly, UP spectra of CO adsorbed at 120 K exhibit well separated 4σ, 1π and 5σ emission peaks at energetic positions which are essentially different from the CO/Ag physisorption system. This observation is interpreted as a demonstration of the ligand effect, i.e., the activation of weak CO bonding at Ag by strong bonding Ti.