面心立方单晶镍纳米线稳定性及磁性的第一性原理计算

本文采用第一性原理方法,研究了轴向为低指数晶向的面心立方(fcc)单晶镍纳米线的稳定性和磁性.计算表明,[110] 是fcc镍纳米线最容易出现的取向,[111] 取向次之,而 [001] 取向则很难出现,这一结果与实验事实符合.镍纳米线按照原子位置和磁性强弱的不同,可以分成简单的芯-壳结构,在纳米线芯部,原子的磁矩大小与块体基本一致.在纳米线表面,镍原子的磁矩比芯部原子有所增加.表面原子磁矩与轴向的取向相关,[110] 为轴向的纳米线表面原子磁矩最低,而[001] 为轴向的纳米线表面原子磁矩最高. 关键词： 镍 纳米线 第一性原理 原子磁矩     

氦团簇及氦间隙原子在钨中的稳定性研究

首先采用分子动力学方法研究了在钨中预存氦-空位团簇（He_nV₂₂）后氦原子结合能与氦-空位比的关系。研究发现:当氦-空位比小于4.5时,氦原子结合能随氦-空位比呈线性减小趋势;当氦-空位比大于4.5时,氦原子的结合能随氦-空位比出现剧烈振荡的现象,这种现象是由于钨中预存氦-空位团簇随机挤出位错环使体系能量骤降所导致的。与此同时,氦-空位团簇周围出现了一些处于亚稳态的fcc结构和hcp结构的钨。为了研究氦团簇周围压强对钨基体相变的影响,本文利用第一性原理对钨的三种结构进行了高压相变计算,发现静水压力不能使钨的三种结构互相转变。另外,通过对bcc钨和fcc钨中四面体间隙氦原子和八面体间隙氦原子电荷密度差的计算,发现bcc钨中四面体间隙氦原子的稳定性高于八面体间隙氦原子的稳定性,而在fcc钨中四面体间隙氦原子的稳定性弱于八面体间隙氦原子的稳定性。     

嵌入原子势计算FCC金属的声子谱与表面驰豫

给出一个简单解析的嵌入原子势模型。该模型被用于计算ＦＣＣ金属的声子谱,表面驰豫,结构稳定性和总能曲线,计算结果与实验及更高水平的计算结果进行了比较。     

基于密度泛函理论的5~10元瓜环结构参数与前线轨道的计算

在气相环境中,使用密度泛函理论（DFT）优化由n个苷脲单元组成的瓜环[n]（CB[n]）（n=5~10）,并使用密度泛函（DFT）概念指数和Multiwfn软件包计算和分析CB[n]的结构参数、前线轨道能量和化学稳定性.结果表明：α-N结构最稳定,α=O、γ-γ、γ-H和β-H（2）结构的化学稳定性较差;以CB[6]为界,主要二面角的变化呈现中心对称的形式;随着苷脲单元n的增加,CB[n]的端口直径、空腔直径和圆外径线性增大;前线轨道E_HOMO、E_LUMO值及E_LUMO-E_HOMO值逐渐降低,化学活性逐渐增强、稳定性逐渐减弱;端口O原子是最大的亲电活性位点,LUMO的电子云分布主要与H原子有关,且次甲基H原子对LUMO电子云的贡献最大;次甲基C原子、亚甲基C原子和指向CB[n]端口方向的亚甲基H原子对LUMO的电子云分布具有抑制作用,抑制能力的大小为指向CB[n]端口方向的亚甲基H原子>亚甲基C原子>次甲基C原子.为研究瓜环的超分子组装提供理论依据.     

Co_(2~-)基Heusler合金Co_2FeAl_(1–x)Si_x(x=0.25,x=0.5,x=0.75)的结构、电子结构及热电特性的第一性原理研究

运用基于密度泛函理论的第一性原理方法,对Co_2FeAl_(1–x)Si_x(x=0.25, 0.5, 0.75)系列Heusler合金的电子结构、四方畸变、弹性常数,声子谱以及热电特性进行了计算研究.结果显示, Co_2FeAl_(1–x)Si_x系列合金的电子结构均为半金属特性,向下自旋态(半导体性)均呈现良好的热电特性,并且随着硅原子浓度的增加功率因子随之增加.计算的声子谱不存在虚频,均满足动力学稳定性条件,弹性常数均满足玻恩稳定性条件,机械稳定性均良好.随着晶格常数c/a的比值变化,体系的能量最低点均出现在c/a=1处,即结构稳定性不随畸变度c/a的变化而变化,说明不存在马氏体相变.此系列合金薄膜的电子结构呈现较高的自旋极化率,在替代浓度x=0.75时自旋极化率达到100%,且当x=0.75时薄膜在畸变度c/a=1.2时存在马氏体相变.随着晶格畸变度的改变,总磁矩也发生变化,且主要由Fe和Co两种过渡金属原子的磁矩变化所决定.     

XAFS和XRD研究高能球磨对Fe70Cu30合金结构的影响

利用XRD和XAFS方法研究机械合金化Fe70Cu30二元金属合金随球磨时间的结构变化.XRD结果表明,球磨2 h后,部分金属Fe与Cu生成Fe-Cu合金;球磨20h后,金属Fe与Cu已完全合金化生成Fe-Cu合金,并只在2θ=44°处出现一个宽化的弱衍射峰,认为是在球磨20h后的Fe70Cu30合金中共存着fcc和bcc结构的Fe-Cu合金相.XAFS结果进一步表明,在球磨的初始阶段(2h),fcc结构的Cu颗粒的晶格产生较大的畸变,其无序度σ(σ＝σT+σS)为0.0190nm.球磨5h后,部分fcc结构的Cu原子进入了无序度相对较小的bcc结构的α-Fe相,导致Cu原子的平均无序度σ降为0.0108 nm.球磨10h后,样品中很大比例的Fe原子处于fcc结构的Fe-Cu合金相,其无序度为σ=0.0119 nm;而大部分Cu原子依然保持fcc结构,无序度为σ=0.0110 nm.这是由于扩散到bcc结构α-Fe相的Cu原子超过某一浓度后(约30%-40%),Cu原子能诱导其产生fcc结构相变所致.球磨时间增加到20h,样品中Cu原子和Fe原子在fcc和bcc相的比例与球磨10h基本相同,生成的Fe-Cu合金混合相的组成和结构分别近似于bcc结构的Fe80Cu20和fcc结构的Fe60Cu40.     

长脉冲激光场中H原子的阈上电离的理论研究

利用基于密度泛函理论下的窗算符技术求解含时薛定谔方程方法计算了H原子在长脉冲强激光场中能量谱,并与文献[8]中的动量空间的投影技术方法的计算结果进行了比较.在文献[8]给出的最大光电子能量范围内,我们的计算结果与文献[8]中的计算结果符合得很好,但是我们的计算结果能够给出比文献[8]中能量更高的光电子能谱,光电子能量谱的变化趋势也与实际的物理过程相符.同时我们也给出了最大光电子能量与激光强度关系的经验公式,利用我们的经验公式计算的最大光电子能量与经典计算预测的最大光电子能量符合得很好,这表明们的计算方法可靠,能够给出高精度的计算结果 .因此利用基于密度泛函理论下的窗算符方法能很好地用来研究H原子在长脉冲强激光场中的阈上电离的性质.     

过渡族BCC金属晶格动力学的改进分析型EAM模型

应用改进分析型EAM多体势,运用晶格动力学理论,具体计算了7种过渡族bcc金属(Cr,Fe,W,Mo,Ta,V,Nb)的[100],[110]和[111]三个晶向声子谱和比热.将计算结果与实验值进行了比较,较Johnson势有很大提高;并从方向性键合角度对符合情况进行了分析.     

TiO2分子在GaN(0001)表面吸附的理论研究

采用基于密度泛函理论的平面波超软赝势法,计算了TiO₂分子在GaN(0001)表面的吸附成键过程、吸附能量和吸附位置. 计算结果表明不同初始位置的TiO₂分子吸附后,Ti在fcc或hcp位置,两个O原子分别与表面两个Ga原子成键,Ga-O化学键表现出共价键特征,化学结合能达到7.932-7.943eV,O-O连线与GaN[1120]方向平行,与实验观测(100)[001] TiO₂//(0001)[1120]GaN一致. 通过动力学过程计算分析,TiO₂分子吸附过程经历了物理吸附、化学吸附与稳定态形成的过程,稳定吸附结构和优化结果一致. 关键词： GaN(0001)表面 2分子')" href="#">TiO₂分子 密度泛函理论 吸附     

FCC金属空位的MAEAM理论研究

应用涉及更远邻原子的改进分析型嵌入原子方法(MAEAM)计算了面心立方(fcc)金属(Ag,Al,Au,Cu,Ir,Ni,Pd,Pt,Rh)的空位性能.在MAEAM计算中,考虑了远邻原子相互作用和单空住迁移能,对两体势进行了坚挺处理,并采用新的截尾函数和加强光滑连接条件对两体势作了截尾处理.同时为了更好的符合面心立方晶体的结合能、弹性常数和平衡条件,调整了多体势的模型常数.未弛豫空住性能计算中考虑了两体势的截尾距离和电子密度分布函数的截尾距离之间近邻原子的作用以及双空位迁移途径周围的原子非对称分布.结果与其它方法计算结果基本一致,但更加接近实验值.对双空住迁移能的计算结果有利地说明了fcc金属双空住5种迁移途径的扩散机制.     

FCC金属空位的MAEAM理论研究

应用涉及更远邻原子的改进分析型嵌入原子方法（MAEAM）计算了面心立方（fcc）金属（Ag,Al,Au,Cu,Ir,Ni,Pd,Pt,Rh）的空位性能。在MAEAM计算中,考虑了远邻原子相互作用和单空位迁移能,对两体势进行了坚挺处理,并采用新的截尾函数和加强光滑连接条件对两体势作了截尾处理。同时为了更好的符合面心立方晶体的结合能、弹性常数和平衡条件,调整了多体势的模型常数。未弛豫空位性能计算中考虑了两体势的截尾距离和电子密度分布函数的截尾距离之间近邻原子的作用以及双空位迁移途径周围的原子非对称分布。结果与其它方法计算结果基本一致,但更加接近实验值。对双空位迁移能的计算结果有利地说明了fcc金属双空位5种迁移途径的扩散机制。     

HIGH pressure structural investigations of zirconium using LMTO method

Abstract

The structural energy differences have been calculated for zirconium as a function of pressure at zero temperature using the Andersen force theorem and the linear muffin tin orbital method. The structures included are the following: α (hcp), the room temperature room pressure phase, ω- a three atom simple hexagonal, bcc and fcc. Our calculations show that the bcc structure would become energetically most favourable above 11 GPa. This results is in agreement with well known correlation between the crystal structure and the d-electron population in transition metals at normal volume. The diamond anvil cell based high pressure x-ray diffraction experiments are in progress to verify this result.     

Theoretical study of structural energy, phonon spectra, and elastic constants of Rh and Ir

The transition metal pair potential (TMPP) is used to study band structure energy of Rh and Ir. Both metals are found to be most stable in fcc structure down to atomic volume 0.5V ₀. The pressure at 0.5V ₀ is found to be 5.235 Mbar and 9.216 Mbar in Rh and Ir, respectively. The TMPP is also used to study other properties of these metals like cohesive energy, phonon frequencies at observed volume. The bulk moduli and elastic constants of these metals at observed volume are calculated by including the volume contribution.     

Lattice mechanical properties of some fcc f-shell metals

A pseudopotential depending on an effective core radius is proposed to study the binding energy, equation of state, ion-ion interaction, phonon dispersion curves (q-space and r-space analysis), mode Grüneisen parameters and dynamical elastic constants of some fcc f-shell metals La, Yb, Ce and Th. The contribution of the s-like electrons is calculated in the second-order perturbation theory for the potential while d and f-like electron is taken into account by introducing repulsive short-range Born-Mayer term. The parameter of the potential is evaluated by zero pressure condition. An excellent agreement between theoretical investigations and experimental findings is achieved which confirms the present formalism     

On the mechanical stability of bcc transition elements and alloys

A volume independent and a volume dependent lattice energy function involving short-range interatomic potentials were able to be fitted to the elastic constants, cohesive energy, lattice parameter and for the latter function to the vacancy formation energy and bcc-fcc lattice stability energy, as well, for fcc metals and bcc alkali metals, but not to the cohesive energy and C' elastic constant of bcc transition metals. The assumption that directional, but partial, covalent bonds exist between nearest-neighbors in the bcc transition elements provides an explanation for the latter results and in addition explains the identical dependence of C' and the bcc-fcc lattice stability upon N_d, where N_d is the average number of d electrons, for the bcc transition metals and alloys. Both the mechanical and thermodynamic stability of the bcc structure for transition metals and all transition metal alloys disappears for 5 ? N_d > 2 and <?1.     

Lattice Dynamics of Some Fcc F-Shell Metals

A model pseudopotential depending on an effective core radius but otherwise parameter free is used to study the binding energy, equation of state, ion-ion interaction, phonon dispersion curves (q-space and r-space analysis), phonon density of states, Debye temperature, mode Grüneisen parameters, dynamical elastic constants, Debye-Waller factor, mean-square displacement, Debye-Waller temperature parameter and propagation velocities of elastic waves of some fcc f-shell metals La, Yb, Ce, and Th. The contribution of the s-like electrons is calculated in the second-order perturbation theory for the model potential while d- and f-like electron is taken into account by introduction of repulsive short-range Born-Mayer term. Very recently proposed screening function due to Sarkar et al. has been used to obtain the screened form factor. The parameter of the potential is evaluated by zero pressure condition. Which is independent of any fitting procedure. An excellent agreement between theoretical investigations and experimental findings prove the ability of the potential for d- and f-shell metals exclusively.     

The influence of anharmonicities to the lattice structure of rare gas crystals

The zero-point energy and the leading anharmonic terms of the free energy retained in second-order perturbation theory are calculated numerically at temperature zero for the two closest-packed structures, the face-centered cubic (fcc) and the hexagonal closest-packed (hcp). In agreement with experiment we get the fcc lattices energetically favoured relative to the hep lattices.     

Roton softening in the solid hydrogens

We present a theory of the roton excitations in high-density J = 0 solid hydrogen. We have calculated dispersion relations and density of states as a function of the density of the solids for both hcp and fcc structures. Increasing the density of the rotationally symmetric fcc solid leads to roton softening. When the X₅ roton energy becomes imaginary symmetry breaking of the ground state occurs, resulting in an orientationally ordered solid.     

Calculation of the surface energy of fcc metals with modified embedded-atom method

The surface energies for 38 surfaces of fcc metals Cu, Ag, Au, Ni, Pd, Pt, A1, Pb, Rh and Ir have been calculated by using the modified embedded-atom method. The results show that, for Cu, Ag, Ni, A1, Pb and Ir, the average values of the surface energies are very close to the polycrystalline experimental data. For all fcc metals, as predicted, the close-packed (111) surface has the lowest surface energy. The surface energies for the other surfaces increase linearly with increasing angle between the surfaces (hkl) and (111). This can be used to estimate the relative values of the surface energy.     

Point-defect properties in HCP rare earth metals with analytic modified embedded atom potentials

The analytic embedded atom method (EAM) type many-body potentials of hcp rare earth metals (Dy, Er, Gd, Ho, Nd, Pr, and Tb) have been constructed. The hcp lattice is shown to be energetically most stable when compared with the fcc and bcc structure, and the hcp lattice with ideal c/a. The mechanical stability of the corresponding hcp lattice with respect to large change of density and c/a ratio is examined. The phonon spectra, stacking fault and surface energy are calculated. The activation energy for vacancy diffusion in these metals has been calculated and the most possible diffusion paths are predicted. Finally, the self-interstitial atom (SIA) formation energy and volume have been evaluated for eight possible sites. This calculation suggests that the crowdion and basal split are the most stable configurations. The SIA formation energy increases linearly with the increase of the melting temperature.Received: 26 March 2003, Published online: 9 September 2003PACS: 34.20.Cf Interatomic potentials and forces - 66.30.Fq Self-diffusion in metals, semimetals, and alloys - 61.72.Ji Point defects (vacancies, interstitials, color centers, etc.) and defect clusters - 61.72.Bb Theories and models of crystal defects