S在Cu(111)表面吸附的第一性原理研究*

基于广义梯度近似的投影缀加平面波(Projector augmented wave) 赝势和具有三维周期性边界条件的超晶胞模型,采用第一原理方法计算并分析了由于S吸附所形成的S/Cu(111)界面体系的吸附结构、吸附能和局域电子结构,考虑了不同覆盖度（1,0.25ML）下S在不同吸附位置的吸附特性. 结果表明：S原子倾向于吸附在高对称的fcc位与hcp位;由于S的负电性而使S/Cu吸附能随覆盖度的减小而增加, 与之相应,S-Cu键长随覆盖度的减小而缩短. DOS图、Bader电荷分析表明杂化主要发生在S的3p态和表面Cu原子的3d态之间,表层近邻的Cu原子向S转移的电子数随覆盖度增加而减小,这表明S与Cu(111)面有强的相互作用.     

Ni(111)表面一氧化碳和氢共吸附的密度泛函理论研究

本文用密度泛函理论(DFT)的总能计算研究了一氧化碳和氢原子在Ni(111)表面上p(2×2)共吸附系统的原子结构和电子态，结果表明CO和H原子分别被吸附于两个对角p(1×1)元胞的hcp和fcc位置.以氢分子和CO分子作为能量参考点，总吸附能为2.81 eV，相应的共吸附表面功函数φ为6.28 eV.计算得到的C—O，C—Ni和H—Ni的键长分别是1.19?, 1.96?和 1.71?，并且CO分子以C原子处于hcp的谷位与金属衬底原子结合.衬底Ni(111)的最外两层的晶面间距在吸附后的相对变化分别是 关键词： Fisher-Tropsch反应 催化作用 Ni(111) p(2×2)/(CO+H) 共吸附     

O在Au(111)表面吸附的密度泛函理论研究

应用密度泛函理论，本文系统地研究了O在Au(111)表面上的吸附能、吸附结构、功函数、电子密度和投影态密度，给出了覆盖度从0.11ML到1.0ML的范围内，O的吸附特性随覆盖度变化的规律.研究发现O的稳定吸附位为3重面心立方(fcc)洞位，O在fcc洞位的吸附能对覆盖度比较敏感，其值随着覆盖度的增加而减小；O诱导Au(111)表面功函数的变化量与覆盖度成近线性关系，原因是Au表面电子向O偏移，形成表面偶极子；O—Au的相互作用形成成键态和反键态，且反键态都被占据，造成O—Au键很弱，O吸附能较小. 关键词： 表面吸附 Au(111)表面 密度泛函理论 电子特性     

Cl原子在γ-TiAl(111)表面吸附的第一性原理研究

基于密度泛函理论,在广义梯度近似下研究了Cl在γ-TiAl(111)表面的吸附.计算结果表明:γ-TiAl(111)表面的面心立方位置(fcc)和六角密排位置(hcp)为Cl吸附的稳定位置,当覆盖度Θ小于一个单层(ML)时,Cl原子倾向于吸附在γ-TiAl(111)表面近邻为多Ti的位置.电子结构分析发现,Cl原子同表面金属原子形成较强的离子键,并且成键具有一定的方向性.当Cl原子和O原子共同在γ-TiAl(111)表面吸附时,二者都趋     

第一性原理研究氧在Ni(111)表面上的吸附能及功函数

采用基于密度泛函理论（DFT）广义梯度近似（GGA）下的第一性原理方法系统地研究了不同覆盖度下O在Ni（111）表面的吸附特性.计算结果表明,O在Ni（111）表面的稳定吸附位为三重面心立方（fcc）洞位,吸附能随着覆盖度的增加而减小,O诱导Ni（111）表面功函数的变化量与覆盖度成近线性关系,并随着覆盖度的增加而增大.同时,通过对电子密度和分波态密度的分析发现：O在Ni（111）表面的吸附使得Ni表面电子向O原子转移,形成表面偶极矩,导致功函数增加;表面Ni原子的3d轨道和O的2p轨道通过耦合、杂化作用形成成键态和反键态,而反键态几乎不被占据,因而O—Ni键相互作用比较强,吸附能较大. 关键词： 表面吸附 密度泛函理论 吸附能 功函数     

H,Cl和F原子钝化Cu_2ZnSnS_4(112)表面态的第一性原理计算

采用基于密度泛函理论的第一性原理计算方法系统研究了Cu_2ZnSnS_4体相的晶格结构、能带、态密度及表面重构与H,Cl和F原子在Cu_2ZnSnS_4(112)表面上的吸附和钝化机理.计算结果表明:表面重构出现在以金属原子Cu-Zn-Sn终止的Cu_2ZnSnS_4(112)表面上,并且表面重构使表面发生自钝化;当单个H,Cl或F原子吸附在S原子终止的Cu_2ZnSnS_4(112)表面上时,相比于桥位(bridge)、六方密排(hcp)位和面心立方(fcc)位点,三种原子均在特定的顶位(top)吸附位点表现出最佳稳定性.当覆盖度为0.5 ML时,无论H,Cl还是F原子占据Cu_2ZnSnS_4(112)表面的2个顶位均具有最低的吸附能.以S原子终止的Cu_2ZnSnS_4(112)表面在费米能级附近的电子态主要由价带顶部Cu-3d轨道和S-3p轨道电子贡献,此即表面态.当H,Cl或F原子在表面的覆盖度达0.5 ML时,费米能级附近的表面态降低,其中H原子钝化表面态的效果最佳,Cl原子的效果次之,F原子的效果最差.表面态降低的主要原因在于吸附原子从S原子获得电子致使表面Cu原子和S原子在费米能级处的态密度峰几乎完全消失.     

铝和银在铱（111）宽范围吸附的稳定性、原子构型及电子特性研究

利用第一性原理密度泛函理论研究了铝和银在铱的111面的宽范围吸附特性。基于密度泛函理论计算了覆盖度在0.11ML到2.00ML的结构稳定性、原子构型及平均结合能。对于铝原子在铱111面的吸附,最稳定的结构是铝原子覆盖度为0.5ML位于密堆六方空位（hcp-hollow）,相应的结合能为-4.68eV;对于亚层铝原子的吸附,最稳定结构是铝原子覆盖度为1.00ML时位于octahedral位置,相应的结合能为-5.28eV。对于覆盖度为2.00ML的满覆盖度混合结构的表层及亚层吸附,最稳定结构是Al位于六方密堆及八方密堆位置,相应的结合能为-4.70eV。这意味着当铝原子以满覆盖度吸附在铱的111面上时,趋向于在铱的111面的亚层形成化学键,而非吸附于表层。相比于铝吸附在铱111面,银的吸附特性呈现出很大的不同,面心位置更为稳定,在覆盖度为0.25ML时其结合能为3.89eV,略微高出密堆六方位置处3.88eV的结合能。     

铝和银在铱（111）宽范围吸附的稳定性、原子构型及电子特性研究

利用第一性原理密度泛函理论研究了铝和银在铱的111面的宽范围吸附特性。基于密度泛函理论计算了覆盖度在0.11ML到2.00ML的结构稳定性、原子构型及平均结合能。对于铝原子在铱111面的吸附,最稳定的结构是铝原子覆盖度为0.5ML位于密堆六方空位（hcp-hollow）,相应的结合能为-4.68eV;对于亚层铝原子的吸附,最稳定结构是铝原子覆盖度为1.00ML时位于octahedral位置,相应的结合能为-5.28eV。对于覆盖度为2.00ML的满覆盖度混合结构的表层及亚层吸附,最稳定结构是Al位于六方密堆及八方密堆位置,相应的结合能为-4.70eV。这意味着当铝原子以满覆盖度吸附在铱的111面上时,趋向于在铱的111面的亚层形成化学键,而非吸附于表层。相比于铝吸附在铱111面,银的吸附特性呈现出很大的不同,面心位置更为稳定,在覆盖度为0.25ML时其结合能为3.89eV,略微高出密堆六方位置处3.88eV的结合能。     

Cu(100) (2×22)R45°-O的表面结构与电子态的密度泛函研究

用第一性原理的总能计算研究了Cu(100))面的表面结构、弛豫以及氧原子的(2×22)吸附状 态.计算给出了Cu(100) (2×22)R45°-O吸附表面的结构参数，并得到了上述结构下氧吸附 的Cu(100)表面氧原子和各层Cu原子的电子态密度.计算得到的吸附表面功函数为4.58 eV ，与清洁Cu(100)表面功函数(～4.53 eV)几乎相同.吸附氧原子与最外层铜原子之间的垂直 距离约为0.02 nm，其能带结构体现出一定的金属性，同时由于Cu-O的杂化作用在费米能以 下约6.4 eV附近出现了局域的表面态.可以认为，在Cu(100) (2×22)R45°的氧吸附表面结 构下，吸附氧原子和衬底之间的结合主要来源于表面最外层铜原子与氧原子的相互作用. 关键词： Cu(100)(2×22)R45°-O表面 缺列再构 表面电子态     

Cs/GaN(0001)吸附体系电子结构和光学性质研究

采用基于第一性原理的密度泛函理论平面波超软赝势方法计算了 1/4ML Cs原子吸附 (2 × 2) GaN(0001) 表面的吸附能、能带结构、电子态密度、电荷布居数、功函数和光学性质. 计算发现, 1/4ML Cs 原子在 GaN(0001) 表面最稳定吸附位为 N 桥位, 吸附后表面仍呈现为金属导电特性, Cs原子吸附GaN(0001)表面后主要与表面 Ga 原子发生作用, Cs6s 态电子向最表面 Ga 原子转移, 引起表面功函数下降. 研究光学性质发现, Cs 原子吸附 GaN(0001) 表面后, 介电函数虚部、吸收谱、反射谱向低能方向移动.     

Coverage-dependent absorption of atomic hydrogen into the sub-surface of Cu(1 1 1) studied by density-functional-theory calculations

With density-functional-theory calculations, we have studied coverage-dependent absorption of H atoms into the sub-surface below a face-centered-cubic (fcc) hollow site of Cu(1 1 1). Both frozen and relaxed surface lattices were considered when the atomic H migrated from the surface to the sub-surface. The potential energy curve for the absorbing H shows that the surface site is in general favored over the sub-surface site, and this trend varies little with the H coverage (0.11-0.67 ML). If the hexagonal-close-packed (hcp) hollow sites immediately vicinal to the absorbing H are pre-adsorbed with other H atoms, the surface adsorption potential is greatly increased, because of the repulsive H-H interaction, to a value near, or even greater than, the sub-surface absorption potential; when two or three H atoms (on the hcp sites) are beside the absorbing H, the energy barrier for the sub-surface absorption is decreased, whereas that for diffusion from the sub-surface to the surface is enhanced. These results indicate that, on an H-saturated Cu(1 1 1) surface (0.67 ML), the sub-surface sites below the fcc sites with two or three neighboring H atoms can trap the sub-surface H.     

Ar adsorptions on Al (111) and Ir (111) surfaces: a first-principles study

We investigate the adsorptions of Ar on Al (111) and Ir (111) surfaces at the four high symmetry sites, i.e., top, bridge, fcc- and hcp-hollow sites at the coverage of 0.25 monolayer (ML) using the density functional theory within the generalized gradient approximation of Perdew, Burke and Ernzerhof functions. The geometric structures, the binding energies, the electronic properties of argon atoms adsorbed on Al (111) and Ir (111) surfaces, the difference in electron density between on the Al (111) surface and on the Ir (111) surface and the total density of states are calculated. Our studies indicate that the most stable adsorption site of Ar on the Al (111) surface is found to be the fcc-hollow site for the (2 × 2) structure. The corresponding binding energy of an argon atom at this site is 0.538 eV/Ar atom at a coverage of 0.25 ML. For the Ar adsorption on Ir (111) surface at the same coverage, the most favourable site is the hcp-hollow site, with a corresponding binding energy of 0.493 eV. The total density of states (TDOS) is analysed for Ar adsorption on Al (111) surface and it is concluded that the adsorption behaviour is dominated by the interaction between 3s, 3p orbits of Ar atom and the 3p orbit of the base Al metal and the formation of sp hybrid orbital. For Ar adsorption on Ir (111) surface, the conclusion is that the main interaction in the process of Ar adsorption on Ir (111) surface comes from the 3s and 3p orbits of argon atom and 5d orbit of Ir atom.     

A first-principles study of surface and subsurface H on and in Ni(1 1 1): diffusional properties and coverage-dependent behavior

Periodic, self-consistent, density functional theory (GGA-PW91) calculations are performed for both surface and subsurface atomic hydrogen on and in Ni(1 1 1). At a low coverage (θ=0.25 ML), the binding energies (BEs) of a hydrogen atom in surface fcc, subsurface octahedral (first layer), and subsurface octahedral (second layer) sites are −2.89, −2.18, and −2.11 eV, respectively. The activation energy barriers for hydrogen diffusion from the surface to the first subsurface layer and from the first to the second subsurface layer are estimated to be 0.88 and 0.52 eV, respectively. In the entire coverage range studied, hydrogen occupies surface fcc and hcp sites and subsurface octahedral sites. In addition, the magnitude of the BE per hydrogen atom and the magnetization of the nickel slabs both decrease as hydrogen coverage increases. Vibrational frequencies of hydrogen at various surface and subsurface sites are calculated and are in reasonable agreement with experimental data. A phase stability calculation with a 2 × 2 surface unit cell shows that a p(2 × 2)-2H overlayer structure (θ=0.5 ML) and a p(1 × 1)-1H structure (θ=1.0 ML) are stable at low hydrogen pressures, in agreement with numerous experimental results. A very large increase in pressure is required to populate subsurface sites. After such an increase occurs, the first subsurface layer is filled completely.     

Extended pattern recognition scheme for self-learning kinetic Monte Carlo simulations

We report the development of a pattern recognition scheme that takes into account both fcc and hcp adsorption sites in performing self-learning kinetic Monte Carlo (SLKMC-II) simulations on the fcc(111) surface. In this scheme, the local environment of every under-coordinated atom in an island is uniquely identified by grouping fcc sites, hcp sites and top-layer substrate atoms around it into hexagonal rings. As the simulation progresses, all possible processes, including those such as shearing, reptation and concerted gliding, which may involve fcc-fcc, hcp-hcp and fcc-hcp moves are automatically found, and their energetics calculated on the fly. In this article we present the results of applying this new pattern recognition scheme to the self-diffusion of 9-atom islands (M(9)) on M(111), where M = Cu, Ag or Ni.     

Adsorption site, core level shifts and charge transfer on the Pd(1 1 1)-I(√3 × √3) surface

We use core level photoelectron spectroscopy and density functional theory (DFT) to investigate the iodine-induced Pd(1 1 1)-I(√3 × √3) structure formed at 1/3 ML coverage. From the calculations we find that iodine adsorbs preferentially in the fcc hollow site. The calculated equilibrium distance is 2.06 Å and the adsorption energy is 68 kcal/mol, compared to 2.45 Å and 54 kcal/mol in the atop position. The adsorption energy difference between fcc and hcp hollows is 1.7 kcal/mol. Calculated Pd 3d surface core level shift on clean Pd(1 l 1) is 0.30 eV to lower binding energy, in excellent agreement with our experimental findings (0.28-0.29 eV). On the Pd(1 1 1)-I(√3 × √3) we find no Pd 3d surface core level shift, neither experimentally nor theoretically. Calculated charge transfer for the fcc site, determined from the Hirshfeld partitioning method, suggests that the iodine atom remains almost neutral upon adsorption.     

Application of low energy ion blocking for adsorption site determination of Na Atoms on a Cu(111) surface

Ion blocking in the low keV energy range is demonstrated to be a sensitive method for probing surface adsorption sites by means of the technique of time-of-flight scattering and recoiling spectroscopy (TOF-SARS). Adsorbed atoms can block the nearly isotropic backscattering of primary ions from surface atoms in the outmost layers of a crystal. The relative adsorption site position can be derived unambiguously by simple geometrical constructs between the adsorbed atom site and the surface atom sites. Classical ion trajectory simulations using the scattering and recoiling imaging code (SARIC) and molecular dynamics (MD) simulations provide the detailed ion trajectories. Herein we present a quantitative analysis of the blocking effects produced by sub-monolayer Na adsorbed on a Cu(111) surface at room temperature. The results show that the Na adsorption site preferences are different at different Na coverages. At a coverage θ = 0.25 monolayer, Na atoms preferentially populate the fcc threefold surface sites with a height of 2.7 ± 0.1 Å above the 1st layer Cu atoms. At a lower coverage of θ = 0.10 monolayer, there is no adsorption site preference for the Na atoms on the Cu(111) surface.