氢、氮、铍、锂原子在缺陷富勒烯CX（x=58-62）中的吸附与渗透

碳富勒烯包合物由于潜在的应用前景受到广泛关注。我们采用第一性原理方法计算了氢、氮、锂、铍原子在不同碳富勒烯中的吸附和穿越。根据原子在不同碳富勒烯笼中的势能曲线,我们给出了原子穿越碳富勒烯笼的势垒和势阱,归纳出原子穿越碳笼机理分为插入机理、渗透机理以及插入机理和渗透机理的混合机理。     

W_(20)O_(58)(010)表面氢吸附机理的第一性原理研究

采用基于密度泛函理论的第一性原理平面波超软赝势方法,在广义梯度近似下,研究了W_(20)O_(58)晶胞、W_(20)O_(58)(010)表面结构及其氢吸附机理.计算结果表明:W_(20)O_(58)晶体理论带隙宽度为0.8 eV,为间接带隙,具有金属性.W_(20)O_(58)晶体中W—O共振较强,以共价键居多.W_(20)O_(58)(010)表面有WO终止(010)表面和O终止(010)表面,表面结构优化后使得W—O键长和W—O—W键角改变,从而实现表面弛豫.分别计算了H_2分子吸附在WO终止(001)表面和O终止(001)表面的WO-L-O_(1c),WO-V-O_(1c),WO-L-O_(2c),WO-V-O_(2c),O-L-O_(1c)和O-V-O_(1c)六种吸附构型,其中WO-L-O_(1c),WO-V-O_(1c)和WO-L-O_(2c)这三种吸附构型不稳定;而WO-V-O_(2c),O-L-O_(1c)和O-V-O_(1c)这三种吸附构型都很稳定,H_2分子都解离成两个H原子,吸附能均为负值,分别为-1.164,-1.021和-3.11 eV.WO-V-O_(2c)吸附构型的两个H原子分别吸附在O和W原子上;O-L-O_(1c)吸附构型的两个H原子,一个与O原子成键,另一个远离了表面.其中O-V-O_(1c)吸附构型最稳定,两个H原子失去电子,为O原子提供电子.分析其吸附前后的态密度,H的1s轨道电子与O的2p,2s轨道电子相互作用,均形成了一些较强的成键电子峰,两个H原子分别与O_(1c)形成化学键,最终吸附反应生成了一个H_2O分子,同时产生了一个表面氧空位.     

氢分子在Mg2Ni(010)面的吸附与分解

本文采用基于密度泛函理论(DFT)的第一原理赝势平面波(PW-PP)方法,对氢分子在Mg2Ni(010)面的吸附与分解进行了研究,发现氢分子以Horl的方式吸附在表面层Ni原子的顶位时吸附能最高,为0.6769 eV,这表明氢分子最可能以Horl的方式吸附在表面层Ni原子的顶位,此时氢分子跟表面的距离(rd)和氢分子的键长(rH)分别为1.6286 A和0.9174 (A).在分子吸附的基础上计算了氢分子沿着选取的反应路径分解时的反应势垒,发现要使氢分子分解需要0.2778 eV的活化能,而氢分子分解时的吸附能为0.8390 eV,分解后两个氢原子的距离为3.1712(A).在分子吸附和分解吸附时氢原子跟正下方的Ni原子都有较强的相互作用,氢原子所得到的电子主要来自氢分子正下方的Ni原子.     

氢分子在Mg2Ni(010)面的吸附与分解

本文采用基于密度泛函理论(DFT)的第一原理赝势平面波(PW-PP)方法,对氢分子在Mg2Ni(010)面的吸附与分解进行了研究,我们发现氢分子以Hor1的方式吸附在表面层Ni原子的顶位时吸附能最高,为0.6769eV,这表明氢分子最可能以Hor1的方式吸附在表面层Ni原子的顶位,此时氢分子跟表面的距离( )和氢分子的键长( )分别为1.6286Å和0.9174Å. 在分子吸附的基础上计算了氢分子沿着选取的反应路径分解时的反应势垒,发现要使氢分子分解需要0.2778eV的活化能,而氢分子分解时的吸附能为0.8390eV,分解后两个氢原子的距离为3.1712Å. 在分子吸附和分解吸附时氢原子跟正下方的Ni原子都有较强的相互作用,氢原子所得到的电子主要来自氢分子正下方的Ni原子.     

InSb(110)表面S,O原子吸附的第一性原理研究

基于第一性原理方法, 采用广义梯度近似的交换关联势, 对InSb材料(110)表面的硫吸附和氧吸附之后体系性质的差异进行了分析. 讨论了两种吸附下的键长、键角、能带结构和态密度的变化, 从理论上论证了硫吸附比氧吸附对InSb红外探测器表面态的钝化有优越性, 有利于工艺上在钝化时的选择. 关键词： 第一性原理 InSb 硫吸附 氧吸附     

第一性原理研究CO在Cu(110)表面的吸附

本文采用基于密度泛函理论的第一性原理方法,并同时考虑范德华力的作用,计算并分析了CO在Cu(110)表面的吸附情况.结果表明:1)CO在两个表面Cu原子的短桥位位置吸附最强,吸附能为1.28 e V.第二稳定吸附位置为表面Cu原子的顶位,吸附能为1.23 e V.CO在其他两个位置,表面两个Cu的长桥位和表面四个Cu的中心位的吸附要弱一些,约为0.86 e V和0.83 e V.2)在Cu表面吸附的CO的C-O键长有部分拉长,这与较强的吸附能和电荷转移相应.3)电荷分析表明所有吸附的CO整体上从衬底上面获得部分电荷,约为0.2个电荷.     

第一性原理研究CO在Cu(110)表面的吸附

本文采用基于密度泛函理论的第一性原理方法, 并同时考虑范德华力的作用, 计算并分析了CO在Cu(110)表面的吸附情况. 结果表明: 1) CO在两个表面Cu原子的短桥位位置吸附最强, 吸附能为1.28 eV. 第二稳定吸附位置为表面Cu原子的顶位, 吸附能为1.23 eV. CO在其他两个位置, 表面两个Cu的长桥位和表面四个Cu的中心位的吸附要弱一些, 约为0.86 eV 和 0.83 eV. 2) 在Cu表面吸附的CO的C-O键长有部分拉长, 这与较强的吸附能和电荷转移相应. 3) 电荷分析表明所有吸附的CO整体上从衬底上面获得部分电荷, 约为0.2 个电荷.     

双磁性中心内嵌富勒烯Y_2C_2@C_(82)-C_2(1)中的超快自旋动力学行为

自旋翻转和自旋转移是实现基于内嵌富勒体系自旋逻辑功能器件设计的先决条件.本文以双磁性中心内嵌富勒烯Y_2C_2@C8_2-C_2(1)体系为例,采用第一性原理计算方法,结合Λ进程理论模型和自编的遗传算法程序,在该内嵌富勒烯体系中分别实现了亚皮秒时间尺度内的自旋翻转和自旋转移过程.计算结果表明,优化后的内嵌Y_2C_2团簇结构和实验得到的各项数据基本吻合,并且会对外部的C8_2-C_2(1)笼结构产生一定的排斥力,但由于富勒烯笼状结构具有很强的稳定性,所以整个体系仍然保持碳笼结构的完整性.通过对自旋密度分布与激光脉冲作用下自旋期望值演化的具体分析,经由Λ进程的自旋翻转是基于两个Y元素的整体自旋翻转;自旋转移则源自两个磁性中心以及碳笼之间在激光脉冲作用下的自旋密度重新分布.本文结果揭示了Y_2C_2@C8_2-C_2(1)体系中的超快自旋动力学机理,可望为基于实际内嵌富勒烯分子的自旋逻辑功能器件设计提供理论依据.     

NO在Ni_xAg_y(x+y=13)团簇表面吸附分解的第一性原理研究

基于第一性原理计算方法,对Ni_xAg_y(x+y=13)团簇的几何结构进行优化后,研究了NO在此类团簇表面不同位置的吸附分解行为,讨论了团簇的束缚能、束缚能的二阶差分、能隙以及吸附前后键长、吸附能、NO分波态密度的变化情况.结果表明,团簇对称性随着Ni的占比变化而变化,稳定性随Ni占比增加而增强,束缚能的二阶差分随着Ni原子的增多呈现奇偶震荡性,NO@Ni_xAg_y(x+y=13)团簇表面吸附行为主要为化学吸附,吸附后N-O键长的变化在0.028?～0.092?之间.对团簇吸附NO的态密度分析发现,吸附后NO的2π*轨道失去电子,1π轨道得到电子,从而导致吸附能的变化.     

原子分子在Ni(100)表面吸附的第一性原理研究

本文采用密度泛函理论,结合周期性平板模型,通过对原子H、N、O、S和C,分子CO、N2、NH3、NO,以及自由基CH3、CH、CH2、OH在Ni(100)表面吸附的研究,比较了它们的吸附能,稳定吸附位点,吸附结构及扩散能垒等信息.这些吸附质与表面结合能力从小到大依次是N2NH3COCH3NOHOHCH2CNSONCHC.在所有的原子中,O原子倾向于吸附在桥位,而其余的原子则倾向于吸附在空位.除N2之外的分子吸附物(CO、NO、NH3),最佳吸附位点均为四重空位,而N2的最稳定吸附位置为顶位.对于自由基吸附物(CH、CH2、CN、OH)而言,它们倾向于吸附在四重空位,而CH3则稳定吸附在桥位.     

Adsorption and diffusion of H2O molecule on the Be(0001) surface: A density-functional theory study

Using first-principles calculations, we systematically study the adsorption behavior of a single molecular H₂O on the Be(0001) surface. We find that the favored molecular adsorption site is the top site with an adsorption energy of about 0.3 eV, together with the detailed electronic structure analysis, suggesting a weak binding strength of the H₂O/Be(0001) surface. The adsorption interaction is mainly contributed by the overlapping between the s and pz states of the top-layer Be atom and the molecular orbitals 1b₁ and 3a₁ of H₂O. The activation energy for H₂O diffusion on the surface is about 0.3 eV. Meanwhile, our study indicates that no dissociation state exists for the H₂O/Be(0001) surface.     

Density functional investigation of metal encapsulated X@C12Si8 heterofullerene (X=Li, Na, K, Be, Mg, Ca, Al, Ga)

The stability and the possible application of our recently reported SiC heterofullerenes inspire the investigation of their further stabilization through ion encapsulation. The endohedral complexes X@C₁₂Si₈, where X=Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺, Ca²⁺, Al³⁺, and Ga³⁺, are probed at the MPWB1K/6-311G? and B3LYP/6-311G* levels of theory. The optimized geometries show the expanding or contracting capability of C₁₂Si₈ in order to accommodate metal ion guests. The inclusion energies indicate the stability of the complexes compared to the components. Meanwhile, the calculated binding energies show the stabilization of C₁₂Si₈ through the inclusion of Be²⁺, Mg²⁺, Al³⁺, and Ga³⁺. The host-guest interaction that is probed through NBO atomic charges supports the obtained results. This study refers to “metal ion encapsulation” as a strategy for stabilization of SiC heterofullerenes.     

Theoretical study on C100 fullerenes and C96X4 (X=N,P, B,Si)

The geometrical structures and electronic properties of six fullerene isomers of C₁₀₀ were studied at the HF/6-31G^? and B3LYP/6-31G^? levels, respectively. The results of the fully optimized calculations show that three C₁₀₀ isomers 449:D₂, 425:C₁ and 442:C₂ are near isoenergetic isomers. The energies and properties of C₁₀₀ hexaanions were calculated. The C₁₀₀^6? (450:D₅) isomer is predicted to be the most stable isomer at the B3LYP/6-31G^? level, and the C₁₀₀^6? (449:D₂) isomer is 44.1 kcal/mol higher in energy. The heterofullerenes C₉₆X₄ (X=N, P, B, Si) formed from the initial C₁₀₀ (449:D₂) have also been investigated at the B3LYP/6-31G^? level. The HOMO–LUMO gaps and aromaticities show that the replacement of fullerene carbon atoms with four heteroatoms can enhance the electronic stabilization of C₁₀₀ (449:D₂).     

Structural and optical properties of the M@C59X cages (X=N,B and M=Li,Na)

Using B3LYP/6-31G* density functional level of theory, the structural and optical properties of the C₆₀ and M@C₅₉X cages have been investigated. Results indicate that the charge on C atoms and band gap of C₆₀ cage are changed dramatically with the substitution of one B or N atom at one of the C sites and the Li and Na atom encapsulations in the C₆₀ cage. The Mulliken analyses show that the charge is transferred completely between the alkali atoms and the C₅₉X cage. The substitutional and encapsulation doping (SED) reduce the optical gaps of the C₆₀ cage. Also, the oscillator strengths of the absorption peaks are dependent on dopant types.     

Adsorption of SO2 molecule on doped (8, 0) boron nitride nanotube: A first-principles study

Adsorptions of SO₂ on Al-, Ca-, Co-, Cu-, Ge-, Ni-, and Si-doped (8, 0) boron nitride nanotube (BNNT) have been studied using first-principles approach based on density functional theory in order to exploit their potential applications as SO₂ gas sensors. The electronic properties of the BNNT-molecule adsorption adducts are strongly dependent on the dopants. The most stable adsorption geometries, adsorption energies, charge transfers, and density of states of these systems are thoroughly discussed. This work reveals that the sensitivity of (8, 0) BNNT based chemical gas sensors for SO₂ can be drastically improved by introducing appropriate dopant. Si is found to be the best choice among all the dopants.     

Al_3X(X=Zr、Ti、Ce、Er)化合物的第一性原理研究

基于密度泛函理论的第一性原理计算方法,对Al_3X(X=Zr、Ti、Ce、Er)化合物生成焓、结合能、态密度、力学性能进行了计算.生成焓、结合能的计算结果表明:Al_3X(X=Zr、Ti、Ce、Er)化合物中,Al_3Zr最易形成,Al_3Ce则最难形成,计算所得的结合能均为负值,化合物具有热力学稳定性,其中Al_3Zr的热力学稳定性最好,Al_3Er则最差.态密度结果分析:Al_3X(X=Zr、Ti、Ce、Er)化合物在费米能级处态密度不为0,表现出金属特性,通过赝能隙大小判断出它们的共价性具有如下关系:Al_3ZrAl_3TiAl_3CeAl_3Er,这与Al-3s、3p电子分别与Zr-4d,Ti-3d,Ce-5d、4f,Er-4f轨道电子具有较为明显的轨道杂化作用有关.力学性能计算结果表明:Al_3Ti、Al_3Zr抵抗体积变形、剪切变形的能力、刚度较Al_3Er、Al_3Ce大.由于Al_3Zr、Al_3Ti、Al_3Ce、Al_3Er具有一定共价性,共价键为方向性键,不利于塑性变形,因此它们均表现为脆性材料. Al_3Zr、Al_3Ti、Al_3Ce、Al_3Er本征维氏硬度大小排序为:Al_3TiAl_3ZrAl_3ErAl_3Ce.     

Electronic structure and magnetic properties of the N monodoping and (Li, N)-codoped ZnO

Using first-principles calculations based on density functional theory, we investigated systematically the electronic structures and magnetic properties of N monodoping and (Li, N) codoping in ZnO. The results indicate that monodoping of N in ZnO favors a spin-polarized state with a magnetic moment of 0.95 μ_B per supercell and the magnetic moment mainly comes from the unpaired 2p electrons of N and O atoms. In addition, it was found that monodoping of N in ZnO is a weak ferromagnet and it is the spin-polarized O atoms that mediate the ferromagnetic exchange interaction between the two N atoms. Interestingly, by Li substitutional doping at the cation site (Li_Zn), the ferromagnetic stability can be increased significantly and the formation energy can be evidently reduced for the defective system. Therefore, we think that the enhancement of ferromagnetic stability should be attributed to the accessorial holes and the lower formation energy induced by Li_Zn doping.     

Half-metallic ferromagnetism in rocksalt and zinc-blende MS (M=Li, Na and K): A first-principles study

First-principles full-potential linearized augmented plane-wave method is used to investigate the electronic structure and magnetic properties of hypothetical zinc-blende and rocksalt LiS, NaS and KS. We find that all the compounds except rocksalt LiS exhibit half-metallic ferromagnetism with an integer magnetic moment of 1.00 μ_B per formula unit. The ferromagnetism results from the spin-polarization of p states of anion S. Total energies calculations indicate the rocksalt phase is lower in energy than the zinc-blende one. The total energy differences are about 0.38, 0.36 and 0.32 eV per formula unit for LiS, NaS and KS, respectively. Meanwhile, it is shown that rocksalt NaS and KS have the half-metallic gaps of 0.22 and 0.41 eV, respectively, and the half-metallic gaps are 0.03, 0.46 and 0.65 eV for zinc-blende LiS, NaS and KS, respectively. We also find the half-metallicity is robust against the lattice contraction up to 7% and 13% for rocksalt NaS and KS, respectively. Although rocksalt LiS is nonmagnetic and metallic at the equilibrium lattice constant, it shows half-metallic ferromagnetism when the lattice constant is larger than 5.40 Å.