Ni,Al及其合金的一种嵌入原子势

嵌入原子势在金属材料的结构及其物性的计算机模拟方面仍然有着重要的作用.针对面心立方结构的镍、铝及其合金,我们拟合了一种简单形式的嵌入原子势.势参数的拟合使用了相应材料的晶格常数,结合能、空位形成能以及三个弹性常数C_(11),C_(12),C_(44).除了一些高频部分的合理偏差之外,使用拟合的势参数得到的声子谱与实验结果符合良好.此外,得到的状态方程也与第一原理的理论结果很好地符合,说明了此嵌入原子势的可靠性.     

第一性原理研究Li_2NH的晶格动力学和热力学性质

采用第一性原理的赝势平面波方法系统地研究了Li_2NH的电子结构、晶格动力学和热力学性质.计算得到的晶格常数与先前的理论和实验结果符合得很好.运用线性响应理论计算了整个布里渊区高对称方向上的声子色散曲线和相应的声子态密度,发现Li_2NH(Pnma)声子色散曲线没有虚频,动力学性能相对最稳定,计算结果和先前实验及理论数据符合得很好.最后,利用得到的声子态密度进一步预测了Li_2NH的热力学性质,包括晶格振动对Helmholtz自由能、内能、熵和热容的贡献,计算结果在一定程度上可为Li-N-H储氢体系的应用提供理论指导.     

高压下铜的声子谱和弹性系数

 以Foiles、Baskes和Daw提出的嵌入原子势（EAM）为基础,并结合Vinet-Rose形式的普适物态方程,建立了金属铜原子间的EAM势。在此基础上,用EAM势和第一原理微扰理论方法计算了高压下铜的声子谱和弹性系数。通过适当地调整嵌入原子势参数,这两种方法几乎能够给出完全一致的结果。而且,晶格振动声子谱和比热的理论值与实验值符合得也非常好。     

UO2声子密度的分子动力学模拟

对UO₂在300～3000 K的分子动力学模拟结果的讨论,计算得到了体系的晶格参数均方位移和UO₂中铀和氧的速度自相关函数得到的声子密度．在模拟计算中,采用了两种原子相互作用势,其中由Basak势计算优化的晶格参数与报道的实验数据符合的非常好;两种势函数得到的均方位移和Bredig态都比较好．最后,讨论了UO₂的声子密度,随着温度的逐渐升高,声子密度的强度是逐渐增大的,并且,振动曲线是随着温度的增加逐渐向低频部分平移并变得逐渐平稳;氧原子的     

赝势对计算石墨烯声子谱线的第一性原理研究

石墨烯是理论与实验方面研究的热点,而探究其声子谱线结构又为研究力学、热力学等提供基础.本文采用基于密度泛函理论的第一性原理,运用不同的交换关联和赝势方法,计算了石墨烯以及石墨的声子谱线.对比研究发现:在声子谱低频率阶段,不同的赝势计算的结果差别很小;而在声子谱的高频率阶段,不同赝势计算的结果差别显著.相对于GGA交换关联,LDA交换关联计算的高频光学支有所软化,计算结果与实验值更加接近.相对于US赝势方法,PAW赝势方法计算的结果与实验值更加接近.综合比较,PAW-LDA赝势的计算结果与实验值最为接近.     

LiFeAs超导体中磁性与声子软化

运用第一性原理密度泛函理论研究了铁基超导体LiFeAs的电子结构和声子谱.计算得到的LiFeAs基态具有涨落的条型反铁磁构型.通过比较LiFeAs在非磁态与条形反铁磁态下的声子态密度,发现,LiFeAs中各向异性自旋互作用的竞争产生了不稳定的自旋密度波和部分晶格位置弛豫,导致Fe和As原子振动模式的软化,从而提高电声子耦合强度.因此,自旋-声子互作用对非常规超导电性有重要贡献. 关键词： 铁基超导体 反铁磁序 超导电性 电声子耦合     

铁基氟化物超导体SrFe_(1-x)Co_xAsF(x=0,0.125)声子特性的第一性原理计算研究

采用基于第一性原理的平面波赝势方法,计算了铁基氟化物及其钴掺杂超导体SrFe1-xCoxAsF(x=0,0.125)在四方非磁态与正交条纹反铁磁态下的声子谱(声子色散曲线、声子态密度)及电-声子耦合常数.计算发现:条纹反铁磁相互作用下的自旋-声子耦合效应强于电-声子耦合作用使声子谱的宽度减小;自旋效应使声子的有效质量增加导致条纹反铁磁态下Fe原子与As原子的耦合振动频率减小.另外,掺杂和自旋效应是提高电-声子耦合常数的两个有效方法,但计算所得超导转变温度远小于实验测量值,表明铁基超导电性非简单的电-声子耦合配对机理.     

NaZn13型Fe基化合物的结构和热力学性质研究

利用Chen-Mbius晶格反演获得的原子间相互作用势,对NaZn13型Fe基金属间化合物进行原子级模拟研究.计算结果表明,Si原子和Co原子均优先占据96i晶位,Si原子和Co原子替代Fe原子后晶体平均结合能降低.随着Co含量的增加,LaFe13-x-yCoySix和NdFe13-x-yCoySix的晶格参数逐渐降低.声子态密度中,稀土原子主要激发低频模,Si原子主要激发高频模.LaFe11.5-yCoySi1.5化合物的德拜温度随Co含量的增加而增高.     

5种金属高压熔化曲线的理论计算

利用从头算法确定了基于对势的平均场模型中的未知势参数, 从而确定了势函数, 然后利用得到的势函数构建平均场, 进一步计算了5种金属材料 Al, Cu, Ni, Na, K在熔点处原子振动的自由体积. 计算的结果表明, 熔化曲线上原子振动的自由体积与原子晶胞体积之比是常数这一假设并不普遍适用; 对该假设修正后可以比较准确地计算材料的熔化曲线. 关键词： 高压熔化 自由体积理论     

NaZn13型Fe基化合物的结构和热力学性质研究

利用Chen-Mbius晶格反演获得的原子间相互作用势,对NaZn₁₃型Fe基金属间化合物进行原子级模拟研究.计算结果表明,Si原子和Co原子均优先占据96i晶位,Si原子和Co原子替代Fe原子后晶体平均结合能降低.随着Co含量的增加,LaFe_13-x-yCo_ySi_x和NdFe_13-x-yCo_ySi_x的晶格参数逐渐降低.声子态密度中,稀土原子主要激发低频模,Si原子主要激发高频模.LaFe_11.5-yCo_ySi_1.5化合物的德拜温度随Co含量的增加而增高. 关键词： 晶格反演 原子间相互作用势 热力学性质 磁致冷材料     

FCC金属结构稳定性和声子谱的MAEAM研究

应用改进分析型嵌入原子方法(MAEAM)计算fcc金属的结构稳定性和声子谱.考虑更远邻原子之间的相互作用,提出新的两体势函数,并采用新的截尾处理函数和加强光滑连接条件.通过拟合单空位迁移能、结合能、三个独立弹性常数及晶体平衡条件,确定了模型参数.在体积不变条件下,计算金属不同结构的能量,说明Ag、Al、Au、Cu、Ir、Ni、Pd、Pt和Rh的fcc结构比较稳定.它们的能量-体积曲线与Rose方程结果基本-致,进-步说明了体积变化时的结构稳定性.同时发现[100]、[110]和[111]三个方向声子谱的计算结果符合实验值和其它计算结果.     

Cu、Ag和Au掺杂AlN的电磁性质的第一性原理研究

采用基于密度泛函理论（DFT）的平面波超软赝势法,研究了Cu、Ag、Au掺杂AlN的晶格常数、磁矩、能带结构和态密度。电子结构表明,Cu、Ag、Au的掺杂使在带隙中引入了由杂质原子的d态与近邻N原子的2p态杂化而成的杂质带,都为p型掺杂,增强了体系的导电性。Cu掺杂AlN具有半金属铁磁性,半金属能隙为0.442eV,理论上可实现100%的自旋载流子注入;Ag掺杂AlN具有很弱的半金属铁磁性;而Au掺杂AlN不具有半金属铁磁性。因此,与Ag、Au相比,Cu更适合用来制作AlN基稀磁半导体。     

Melting curves of FCC-metals by cell-theory

The cell model is developed to account for many body interaction effects, as occurring in the embedded atom method (EAM) potentials, and crystal lattice features to determine the free energy of FCC metals. The well known smearing approximation, which is generally used in conjunction with pair potentials, is also developed for EAM potentials. The free energy so obtained is used to determine melting curves of FCC metals. For this purpose, the liquid phase free energy is calculated using the corrected rigid spheres model. The advantage of our scheme, which is verified with data on Lennard-Jones solids, is that both free energies are based on the same interaction potential. Results match well with the available experimental/theoretical data for Al and Cu. A good agreement is also found for the transition metals, Ni and Pt, for which molecular dynamics as well as theories like Lindemann's law do not give accurate results. It is also found that smearing approximation, which neglects interactions beyond nearest neighbors, also provides good estimates for free energy if many body effects are accounted with EAM potential.     

Improvement of modified analytic embedded atom method potentials for noble metals and Cu

The modified analytic embedded atom model (EAM) potentials considering farther neighbor atoms are improved for the noble metals (Ag, Au, Pt, Pd, Rh) and Cu. We not only adopt an end processing function and an enhanced smooth continuous condition for the pair potential, but also adjust the model parameters of multi-body potential by fitting a cohesive energy, a mono-vacancy formation energy, the Rose equation curve for the cohesive energy as a function of lattice parameter, a structure energy difference, elastic parameters and an equilibrium condition of crystal. The calculation results of structure energy differences misfit the experiment data for the noble metals and Cu in the unimproved EAM, because anyone of these differences have not been considered in the calculation of its model parameters. After the modification, the model showed better simulation results for the noble metals and Cu.     

Reconstructed (1 1 0) surfaces of FCC transition metals

The energies of the ideal, missing row (MR) and missing column (MC) (1 1 0) surfaces have been calculated by using modified embedded atom method (MEAM) for seven face centered cubic (FCC) transition metals Au, Pt, Ag, Pd, Rh, Cu and Ni. The results, that the MC reconstruction can not be formed for all metals, while the MR reconstruction can be formed naturally for Au and Pt, inductively for Ag, Pd, Rh and Cu and difficultly for Ni, are better than EAM calculated results in comparing with experimental results. In addition to the surface energy explanation, the results are also related to the surface topography and valence electron structure.     

ADATOM, VACANCY AND SPUTTERING YIELDS OF ENERGETIC Pt ATOMS IMPACTING ON Pt(100) BY MOLECULAR DYNAMICS SIMULATION

We have studied the influence of incident atoms with low energy on the Pt(100) surface by molecular dynamics simulation. The interaction potential obtained by the embedded atom method (EAM) was used in the simulation. The incident energy changes from 0.1eV to 200eV, and the target temperature ranges from 100 to 500 K. The target scales are 6×6×4 and 8×8×4 fcc cells for lower and higher incident energies, respectively. The adatom, sputtering, vacancy and backscattering yields are calculated. It was found that there is a sputtering threshold for the incident energy. When the incident energy is higher than the sputtering threshold, the sputtering yield increases with the increase of incident energy, and the sputtering shows a symmetrical pattern. We found that the adatom and vacancy yields increase as the incident energy increases. The vacancy yields are much higher than those obtained by Monte Carlo simulation. The dependence of the adatom and sputtering yields on the incident energy and the relative atomistic mechanisms are discussed.     

Liquid/vapor surface tension of metals: embedded atom method with charge gradient corrections

Molecular dynamics simulations for three embedded atom method (EAM) function sets are used to determine the liquid/vapor surface tension gamma for Al, Ni, Cu, Ag, and Au. The three EAM models differ in both the functional forms employed and the fitting procedure used. All the EAM potentials underestimate gamma but one of the models performs consistently better than the others. We show that including a correction to the local charge density associated with gradients in the density together with exploiting the invariance of the EAM potentials to appropriate transformations in the charge density can lead to improved values for gamma, as well as for solid free surface energies, within existing EAM function sets.     

辐照损伤模拟中的Zr 原子势函数

以金属Zr 的一种嵌入原子形式(EAM) 的势函数为基础,通过引入一个调制函数的办法,在不同范围内修改了EAM 势函数的对势部分和原子电子密度分布部分,然后采用分子动力学方法计算点缺陷(间隙原子(SIA) 和空位) 形成能和初级离位原子(PKA) 的阈值能,从而探讨这些物理量对势函数的不同部分的敏感程度。计算结果表明:势函数的对势部分长程范围内的形式对缺陷的形成影响较小,其短程范围的形式对SIA 的形成比对空位的形成影响程度更大;对于PKA 的阈值能,其敏感区域来自于势函数的对势部分和原子电子密度分布的短程范围部分,但在不同晶向上的PKA 的阈值能对势函数的敏感程度有所不同。这些研究结果对于在研究Zr 金属的辐照损伤中势函数的选择或构建有指导意义。For the purpose of detecting the sensitive parts of an embedded atom method(EAM) potential which is considered to be used in molecular dynamics simulation of radiation effects of Zirconium, we introduce a modulation function to modify the pairwise potential and the atomic electron-density distribution of the EAM potential. Based on the modified potential function, the formation energies of the self-interstitial atom (SIA) and the vacancy atom are calculated as well as the displacement threshold energy of primary knock-on atom (PKA). The results indicate that the short range part of the pairwise potential has more greater influence on the SIAs formation than the vacancy formation. The defect formation energies are also very sensitive to the behavior of the atomic electron-density function in the range which is close to the cutoff distance. The displacement threshold is sensitive to the short range behaviors of both the airwise potential and the atomic electron-density function, however, the sensitivity is strongly dependent on the crystal-direction.     

Chen’s lattice inversion embedded-atom method for Ni-Al alloy

The structural properties,the enthalpies of formation,and the mechanical properties of some Ni-Al intermetallic compounds(NiAl,Ni3Al,NiAl3,Ni5Al3,Ni3Al4) are studied by using Chen’s lattice inversion embedded-atom method(CLI-EAM).Our calculated lattice parameters and cohesive energies of Ni-Al compounds are consistent with the experimental and the other EAM results.The results of enthalpy of formation indicate a strong chemical interaction between Ni and Al in the intermetallic compounds.Through analyzing the alloy elastic constants,we find that all the Ni-Al intermetallic compounds discussed are mechanically stable.The bulk moduli of the compounds increase with the increasing Ni concentration.Our results also suggest that NiAl,Ni3Al,NiAl3,and Ni5Al3 are ductile materials with lower ratios of shear modulus to bulk modulus;while Ni3Al4 is brittle with a higher ratio.     

Development of lattice inversion modified embedded atom method and its applications

The modified embedded atom method (MEAM) has been widely used in describing the physical properties of elemental crystals, alloys and compounds with multiple lattice structures. We report here the development of a reliable procedure to reduce the complexity of the MEAM formalism by removing the many-body screening function. In the proposed formulation, the interatomic pair potential is obtained by applying Chen-Möbius lattice inversion up to fifth nearest neighbors, so that the cohesive energy curve can be reproduced faithfully. The newly developed model (Lattice Inversion MEAM, LI-MEAM), which can be viewed as a direct extension of the embedded atom method (EAM), no longer requires the computation of many-body screen functions and has fewer adjustable parameters than MEAM. As an illustration, we optimized the potential parameters of body centered cubic iron (bcc-Fe). The values of the calculated physical properties agree well with experimental results. We further investigated the size-dependent melting behavior of bcc-Fe nanoparticles (NPs) with particle size ranging from 725-atom (∼25 Å) to 22899-atom (∼80 Å) using replica exchange molecular dynamics (REMD) simulations. Our simulations show advantages of LI-MEAM in modeling of the melting process and quantitatively reveals that the liquid skin melting (LSM) process of bcc-Fe NPs.