立方相Ag_3PO_4(111)面原子几何及弛豫结构的第一性原理计算

采用平面波超软赝势方法研究了立方相Ag_3PO_4(111)面的表面能和表面原子弛豫结构.首先对Ag_3PO_4(111)面的八种不同原子终止结构的体系总能量进行计算,结果表明B种表面模型被证实为最稳定的(111)面原子几何结构.针对该表面结构,探讨了表面能和原子弛豫与模型中原子层数和真空厚度的关系,当原子层数为24层,真空厚度为0.6 nm时,表面能收敛于1.41 J/m2(LDA-CAPZ)和1.39 J/m2(GGA-PBE).表面原子弛豫后,表面两个三配位的Ag原子均向里移动,超过0.06 nm,而表面次层的O原子则均向外移动约0.0042 nm,导致弛豫后暴露在最表面的是O原子,同时表面原子的核外电子向表面内部发生转移,结构趋于稳定.这些结果为进一步深入研究Ag_3PO_4表面的光催化活性起源提供理论支持.     

锐钛矿型TiO2(101)面原子几何及弛豫结构的第一性原理计算

采用平面波超软赝势方法计算了锐钛矿型TiO2(101)面的表面能和表面原子弛豫结构.首先对TiO2(101)面的6种不同的表面原子终止结构的体系总能量进行了计算,结果表明终止原子为两配位的O原子、次层为五配位的Ti原子的表面结构最为稳定.针对该表面研究了表面能和原子弛豫与模型中原子层数和真空厚度的关系,当原子层数为12层,真空厚度为0.4 nm时,表面能收敛度小于0.01 J/m2.研究发现:表面上两配位的O原子向里移动约0.0012 nm,五配位的Ti原子向里移动约0.0155 nm,弛豫后的Ti-O键长发生了较大变化,核外电荷发生了转移,结构趋于稳定.     

锐钛矿型TiO2(101)面原子几何及弛豫结构的第一性原理计算

采用平面波超软赝势方法计算了锐钛矿型TiO₂(101)面的表面能和表面原子弛豫结构.首先对TiO₂(101)面的6种不同的表面原子终止结构的体系总能量进行了计算,结果表明终止原子为两配位的O原子、次层为五配位的Ti原子的表面结构最为稳定.针对该表面研究了表面能和原子弛豫与模型中原子层数和真空厚度的关系,当原子层数为12层,真空厚度为0.4nm时,表面能收敛度小于0.01J/m².研究发现:表面上两配位的O原子向里移动约0.0012nm,五配 关键词： 第一性原理 2')" href="#">TiO₂ 表面结构 弛豫     

Cu 表面性质的第一性原理分析

采用第一性原理赝势平面波方法, 计算并详细分析了面心立方Cu晶体及其 (100), (110) 和 (111) 这3个低指数表面的原子结构、 表面能量及表面电子态密度. 表面能的计算结果表明, Cu (111) 表面的结构稳定性最好, Cu (100) 表面次之, Cu (110)表面的结构稳定性最差. 3个表面的表面原子弛豫量随着层数的增加而逐渐减弱. Cu (110) 表面的最表层原子相对收缩最大, Cu (100)表面次之, Cu (111) 表面的最表层原子相对收缩最小. 表面原子弛豫不仅引起表面几何结构的变化, 而且使表面层原子的电子态密度峰形相对晶体内部发生变化, 这是表面能产生的主要原因, 而Cu (110)表面相对于Cu (100)与Cu (111)表面具有高表面活性的主要原因则源于其表面层原子电子态密度在高能级处的波峰相对晶体内部显著的升高. 关键词： Cu 晶体 表面结构 表面能 态密度     

c(2×2)O吸附Cu(001)表面结构、电子态与STM图像的研究

用投影子缀加波和CP分子动力学方法研究了贵金属Cu(001)面的表面结构、弛豫以及O原子的c(2×2)吸附状态. 研究结果得出在这种吸附结构中，O原子与衬底Cu原子之间的垂 直距离约为0069nm，Cu—O键长为0.194nm，功函数约为5.29 eV；吸附O原子形成金属性能带结构，由于Cu—O的杂化作用,在费米能以下约6.7 eV处出现了局域的表面态.用Tersoff-Hamann途径计算了该表面的扫描隧道显微镜图像，并讨论了与实验结果之 间的关系. 关键词： Cu(001)-c(2×2)/O 电子态 STM图像     

Ag/Cu(111)与Au/Cu(111)体系异质外延岛演化过程中的应变释放及其对Moiré结构形成的影响

采用嵌入原子方法的原子间相互作用势，利用分子动力学模拟方法研究了Au/Cu(111)和Ag/Cu (111)体系的异质外延结构特征以及外延岛形貌和应变释放的演化过程. 通过对比Au/Cu(111)和Ag/Cu (111)体系的异质外延结构及外延岛演化行为，揭示了导致Ag/Cu (111)体系中异质外延层形成Moiré结构的微观物理机理及其与外延体系的宏观物理特性之间的关系. 研究结果显示，外延岛原子与基体表面原子之间的界面结合强度是形成Moiré结构的重要因素，异质外延体系的界面结合强度取决于二者的合金熔解热. 当异质外延体系的合金熔解热为正值时，界面结合强度较弱，有利于Moiré结构的形成. 同时，外延岛原子之间的相互作用决定着外延岛的面内弛豫行为，对Moiré结构的形成有一定的影响. 外延岛的面内弛豫行为与外延层和基体之间的相对刚度有关，弹性模量较大的外延层具有较强的延展能力，对Moiré结构的形成有利. 此外，Moiré结构的形成与外延岛的尺度有关，主要是外延岛边界原子的钉扎作用对外延岛内原子弛豫行为的约束作用的影响.     

Ag(111)和Au(111)上铋的初始生长行为

半金属铋(Bi)的表面合金具有的Rashba效应,和其具体结构性质有重要关联.本文结合扫描隧道显微镜(STM)和密度泛函理论(DFT),系统地研究了Bi原子在Ag(111)和Au(111)上的不同初始生长行为.在室温Ag(111)上,连续的Ag₂Bi合金薄膜会优先在Ag台阶边缘形成;在570 K Ag(111)上,随着覆盖度增加到0.33分子层(ML),Bi优先取代配位数低的台阶边原子并从单原子随机分布转变为长程有序的Ag₂Bi合金相;随着覆盖度增加,Ag₂Bi通过退合金过程转变成p×3^1/2结构的Bi膜.Bi在室温和570 K的Au(111)上的生长行为一致:在覆盖度低于0.40 ML时,Bi会优先吸附在配位为5的Au原子上,并以单原子和团簇的形式分别分散在Au(111)的密堆积区域和鱼骨纹重构的拐角处;随着覆盖度增加到0.60 ML,无序的Bi会逐渐转变成长程有序的((37)^1/2×(37)^1/2)相;Bi的吸附会导致Au(111)表面应力逐步释放.Bi在Ag(111)和Au(111)上的不同生长行为表明,Bi原子与衬底之间的相互作用起着关键作用.     

α-Al2O3(0001)基片表面结构与能量研究

对α-Al2O3(0001)晶体表层三种不同终止原子结构的计算模型,在三维周期边界条件下的κ空间中,采用超软赝势平面波函数描述多电子体系.应用基于密度泛函理论的局域密度近似,计算了不同表层结构的体系能量,表明最表层终止原子为单层Al的表面结构最稳定.对由10个原子组成的菱形原胞进行了结构优化,得到晶胞参数值(a0＝0.48178nm)与实验报道值误差小于1.3%.进一步计算了超晶胞(2×2)表面弛豫,弛豫后原第2层O原子层成为最表层; 对不同表层O,Al原子最外层电子进行了布居分析,表面电子有更大的概率被定域在O原子的周围,表面明显地表现出O原子的电子表面态.     

GaN(100)表面结构的第一性原理计算

用全势缀加平面波加局域轨道(APW+lo)的方法计算了六方GaN及其非极性(1010)表面的原子及电子结构.计算出的六方GaN晶体结构参数晶格常数和体积弹性模量与实验值符合得很好.用平板超原胞模型来计算GaN(1010)表面的原子与电子结构,结果表明表面顶层原子发生键长收缩并扭转的弛豫特性.表面阳离子向体内移动,趋向于sp2平面构形;而表面阴离子向体外移动,趋向于锥形的p3构形.弛豫后,表面实现由半金属性向半导体性的转变.并且,表面电荷发生大的转移,参与表面键的重新杂化,使得表面原子的离子性减弱共价性增强,认为这就是表面原子键收缩并旋转的原因.     

H原子在Al(111)和Ag(111)面上的吸附

本文用电荷自洽的EHT方法,研究了H原子在Al(111)和Ag(111)面上的吸附,结果指出:在Al(111)面上,H以原子状态吸附在某些对称位置上,它也能渗透到表面层中去,成为填隙原子;H₂分子在表面处发生解离吸附。在Ag(111)表面上,H原子有可能以分子状态吸附,H—H键平行于表面,这与高分辨率电子能量损失谱所得到的实验结果一致;但H₂分子在Ag(111)表面也可能发生解离吸附。 关键词：     

Kinetics of surface segregation in alloys

The kinetics of surface segregation in ordering alloys are studied with the help of a multilayer model in slabs of different thicknesses. The time evolution of the concentrations of atomic layers perpendicular to the (100) planes of an AB-type BCC ordering alloy (CoFe) are calculated. As a result of the competition between the surface segregation of A atoms and ordering in the bulk, a metastable configuration with two anti-phase boundaries (APB) inside the slab, with odd numbers of layers, was obtained if we started from an ordered initial state with B atoms on the free surfaces (B-termination). The effect of the temperature was also investigated for a slab with 41 layers and from the calculation of the free energies of the metastable and stable states it was shown that the metastable state can exist below T₁ = 0.95T_c, where T_c is the critical temperature of the order-disorder transition in bulk. It was obtained that T₁ decreases with decreasing thickness of the slab, while the difference of free energies of the metastable and stable states slightly increases. The effect of the slab thickness is a typical size effect: the time necessary to reach the steady state decreases with decreasing number of layers.     

Relaxation of the (1 1 1) surface of δ-Pu and effects on atomic adsorption: An ab initio study

The computational formalism of the full-potential all-electron linearized augmented plane wave plus local orbitals (FP-LAPW+lo) method has been employed to study the relaxation of the δ-Pu(1 1 1) surface and the consequent effects for atomic adsorption of C, N, and O atoms on this surface. The underlying theoretical principle is the generalized gradient approximation to density functional theory (GGA-DFT) and the surface was modeled by a five-layer slab with a (2×2) surface unit cell. Upon relaxation of the slab, the interlayer separation between the surface and the subsurface layers expanded by 7.1% with respect to the bulk interlayer separation, while the separation between the subsurface and central layers expanded by 0.4%. To study adsorption on the surface, the adatoms were allowed to approach the surface at four high symmetry adsorption sites, namely, the top, bridge, hollow FCC, and hollow HCP sites, the adlayer structure corresponding to a coverage of 0.25 of a monolayer in all cases. The hollow FCC adsorption site was found to be the most stable site for C and N with chemisorption energies of 6.420 and 6.549 eV, respectively, while the hollow HCP adsorption site was found to be the most stable site for O with a chemisorption energy of 7.858 eV. The respective distances of the C, N, and O adatoms from the surface were found to be 1.22, 1.09, and 1.22 Å. The work function and net magnetic moments, respectively, increased and decreased in all cases upon chemisorption compared with the bare δ-Pu(1 1 1) surface. The electronic structure of the interactions between the adsorbates and the substrate is discussed in detail.     

Theoretical study of the initial stage of InN growth on cubic zirconia (111) substrates

An initial stage of InN growth on cubic zirconia (111) substrates has been investigated using first‐principles calculations based on density functional theory (DFT). We have evaluated adsorption energies of indium and nitrogen atoms on cubic zirconia (111) surfaces, and have found that the differences in the adsorption energies of the indium atoms at various adsorption sites were small, indicating that the migration of the indium atoms on zirconia (111) surfaces occurs readily. On the other hand, we have found that the differences in the adsorption energies of the nitrogen atoms at various adsorption sites were large, implying that the nitrogen atoms tend to stay at the stable site with the largest adsorption energy, which was identified as the O–Zr bridge site. These results suggest that the first layer of InN films is the nitrogen layer. In addition, we have found that the energetically favorable arrangement is comprised of InN(0001)//cubic zirconia (111) and InN $ [11\bar 20] $ //cubic zirconia $ [1 \bar 10], $ which is quite consistent with previously obtained experimental data. Furthermore, the hybridization effect between N 2p and O 2p plays a crucial role in determining the interface structure for the growth of InN on cubic zirconia (111) surfaces. (© 2012 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Imaging the surface and the interface atoms of an oxide film on Ag111 by scanning tunneling microscopy: experiment and theory

High-resolution images of the vacuum oxide surface atomic structure of Ag111-p(4x4)-O reveal large terraces of a perfect (4x4) reconstruction. Under certain conditions, the surface Ag atoms at the interface between the reconstructed oxide layer and the underlying Ag111 lattice are also imaged, providing the structural registry. Scanning tunneling microscopy simulations reveal a strong sensitivity to structure and comparison with the experimental images, therefore providing an atom-by-atom model for the entire metal-oxide-vacuum structure.     

Multilayer relaxation of fcc metals (001) surface: A modified embedded atom method study

The (001) surface multilayer relaxation results calculated by the modified embedded atom method (MEAM) show that Ni, Al, Rh and Ir (001) surface are ‘anomalous’ outward relaxation, while Cu, Ag, Au, Pd, Pt and Pb (001) surface are inward relaxation. For the inward relaxation metals, the relaxation between the first two layers increase for the 3d, 4d and 5d metals at the same column in the periodic table, successively. The expansion (contraction) between the first two layers at fcc (001) surfaces is accompanied by the decrease (increase) in the electronic density at the lattice of the first two layers. The surface energies results show that the surface energies decrease for all fcc (001) surfaces due to relaxation, whereas the changes not more than 5%.     

One-dimensional surface reconstruction as an atomic-scale template for the growth of periodically striped Ag films

The role of the In/Si(111)-(4 x 1)-In surface as an atomic-scale geometrical template for the growth of Ag thin films is clarified by scanning tunneling microscopy and low energy electron diffraction. Low-temperature grown Ag films are found to have stripe structures with a transverse periodicity equal to that of indium chains of the In/Si(111)-(4 x 1)-In. The stripes exhibit a structural transformation at the thickness of 6 monolayers (ML); this relaxation allows the stripes to persist up to a thickness as large as 30 ML (approximately = 7 nm) while maintaining their mean periodicity. We attribute this stability to a coincidental matching of the periodicity and the corrugation amplitude between the Ag film and the substrate, which is realized by periodic insertion of stacking faults into a Ag fcc crystal.     

α-Al2O3(0001)基片表面结构与能量研究

对α-Al₂O₃(0001)晶体表层三种不同终止原子结构的计算模型， 在三维周期边界条件下 的κ空间中，采用超软赝势平面波函数描述多电子体系.应用基于密度泛函理论的局域密度 近似，计算了不同表层结构的体系能量，表明最表层终止原子为单层Al的表面结构最稳定. 对由10个原子组成的菱形原胞进行了结构优化，得到晶胞参数值(a₀＝0.48178n m)与实验 报道值误差小于1.3%.进一步计算了超晶胞(2×2)表面弛豫，弛豫后原第2层O原子层成为最 表层； 对不同表层O，Al原子最外层电子进行了布居分析，表面电子有更大的概率被定域在 O原子的周围，表面明显地表现出O原子的电子表面态. 关键词： 2O₃(0001)')" href="#">α-Al₂O₃(0001) 超软赝势 表面结构 表面态     

Surface properties of rutile TiO2(1 1 0) from molecular dynamics and lattice dynamics at 300 K: Variable-charge model results

Geometric structure, atomic vibrations and atomic charges and their thermally induced fluctuations have been calculated as a function of depth in, and thickness of, rutile TiO₂(1 1 0) slabs, within the framework of the variable-charge potential of Swamy and Gale [V. Swamy, J.D. Gale, Phys. Rev. B 62 (2000) 5406] at 300 K. Molecular dynamics simulations and lattice dynamics calculations were performed with a 2D periodic slab model for slab thicknesses between 3 and 11 triple layers (approximately 9-35 Å). Odd-even oscillations with respect to the number of slab layers are found for the surface relaxation for very thin slabs, and for the (slowly converging) rumpling in the middle of the slab. The Ti and O atomic charges in the outermost three atomic layers differ from the rest of the slab (they are less ionic); the thermal vibrations do not alter this picture. The atomic mean-square amplitudes are some 50% larger (more for O, less for Ti) at the surface than in the middle of the slab and decay rather slowly to the bulk values. Comparisons with the results of a rigid-ion potential for titania [M. Matsui, M. Akaogi, Mol. Simul. 6 (1991) 238] are presented for non-electronic properties.