Molecular dynamics simulations of displacement cascades in Feben 10%Cr systems

Molecular dynamics simulations of the displacement cascades in Fe-10%Cr systems are used to simulate the primary knocked-on atom events of the irradiation damage at temperatures 300, 600, and 750 K with primary knocked-on atom energies between 1 and 15 keV. The results indicate that the vacancies produced by the cascade are all in the central region of the displacement cascade. During the cascade, all recoil Fe and Cr atoms combine with each other to form Fe-Cr or Fe-Fe interstitial dumbbells as well as interstitial clusters. The number and the size of interstitial clusters increase with the energy of the primary knocked-on atom and the temperature. A few large clusters consist of a large number of Fe interstitials with a few Cr atoms, the rest are Fe-Cr clusters with small and medium sizes. The interstitial dumbbells of Fe-Fe and Fe-Cr are in the lan111ran and lan110ran series directions, respectively.     

过渡金属掺杂氧化锌团簇的物性研究

本文采用第一性原理密度泛函理论研究了过渡金属（TM）原子Cr和Fe单掺杂和双掺杂(ZnO)12团簇的结构和磁性质。我们考虑了替代掺杂和间隙掺杂。结果表明Cr 和 Fe间隙掺杂团簇结构最稳定。团簇磁矩主要来自TM原子3d态的贡献,4s 和4p 态也贡献了一小部分磁矩。由于轨道杂化,相邻的Zn和O原子上也产生少量自旋。最近邻TM原子间的磁性耦合,主要由两个TM原子之间的直接短程铁磁耦合和TM和O原子之间通过p-d杂化产生的反铁磁耦合这两种相互作用的竞争来决定。不同TM原子掺杂团簇的总磁矩与TM原子种类以及掺杂位置有关,说明在(ZnO)12团簇中掺杂不同TM原子在可调磁矩的磁性材料的领域有潜在应用价值。     

Fe掺杂ZnTe团簇结构和磁性质

本文采用第一性原理密度泛函理论系统的研究了Fe原子单掺杂和双掺杂(ZnTe)12团簇的结构和磁性质。我们考虑了替代掺杂和间隙掺杂。不管是单掺杂还是双掺杂,间隙掺杂团簇都是最稳定结构。团簇磁矩主要来自Fe-3d态的贡献,4s 和4p 态也贡献了一小部分磁矩。由于轨道杂化,相邻的Zn和Te原子上也产生少量自旋。最重要的是,我们指出间隙双掺杂团簇是铁磁耦合,在纳米量子器件有潜在的应用价值。     

Fe掺杂ZnTe团簇结构和磁性质

本文采用第一性原理密度泛函理论系统的研究了Fe原子单掺杂和双掺杂（ ZnTe）12团簇的结构和磁性质。我们考虑了替代掺杂和间隙掺杂。不管是单掺杂还是双掺杂,间隙掺杂团簇都是最稳定结构。团簇磁矩主要来自Fe－3d态的贡献,4s和4p态也贡献了一小部分磁矩。由于轨道杂化,相邻的Zn和Te原子上也产生少量自旋。最重要的是,我们指出间隙双掺杂团簇是铁磁耦合,在纳米量子器件有潜在的应用价值。     

Ion mixing in multilayer films and the doping of the near-surface layers of polycrystalline substrates under irradiation by ion beams with a wide energy spectrum

Various factors influencing ion mixing efficiency are considered. It is established that the energy transfer from ions and primary knocked-on atom films to subsequent displacement cascades underlies the penetration of atoms from multilayer films into a polycrystalline substrate. The penetration of implanted atoms to great depths is in this case determined by the defect density, the radiation induced migration of the implanted atoms, and their interaction with the atoms of the matrix. All of these factors can be described in terms of the isotropic mixing model. It is shown that when doped by atoms from multilayer films, gradient layers form that are determined by the range of radiation in the substrate of the atoms to be implanted and their migration under irradiation by an ion beam with a wide range of energies.     

Fe基体中包含Cu团簇的Fe-Cu二元体系在升温过程中结构变化的原子尺度计算

采用基于嵌入原子方法的分子动力学方法模拟了具有体心立方晶格结构的Fe基体中包含小尺寸Cu纳米粒子的Fe-Cu二元体系在升温过程中的原子堆积结构变化.进行了Cu原子均方位移、Cu原子对分布函数和原子的径向密度分布函数的计算,并对纯Cu原子区、Fe-Cu界面区和纯Fe基体区的分区域原子堆积结构进行了分析.结果表明,Fe基体内Cu团簇的尺寸及其在Fe基体内所能占据区域的大小,对不同温度下的Cu团簇内原子堆积结构及Fe基体的原子堆积结构具有影响.升温过程中不同尺寸受基体约束Cu团簇对Fe基体结构改变的影响表现出很大差异.对于Fe_(bulk)-Cu_(135)体系,基体的应变临近Fe-Cu界面区,同时在团簇中间的基体区域出现大量空位缺陷和应变集中区;对于Fe_(bulk)-Cu_(141)体系,随温度升高,基体中出现的应变区域表现为小尺寸、数量多向大尺寸、小数量的变化.     

Radiation-induced glide motion of interstitial clusters in concentrated alloys

We propose a mechanism for glide motion, i.e. one-dimensional (1D) migration, of interstitial clusters in concentrated alloys driven by high-energy particle irradiation. Interstitial clusters are fundamentally mobile on their respective 1D migration tracks, but in concentrated random alloys they are stationary at the position where the fluctuating formation energy achieves a local minimum. Irradiation changes the microscopic distribution of solute atoms through atomic displacement and recovery of the produced Frenkel pairs, which causes cluster 1D migration into a new stable position. In molecular dynamics simulations of interstitial clusters up to 217i in Fe–Cu alloys, stepwise 1D migration was observed under interatomic mixing or shrinkage of the cluster: a single 1D migration was induced by two exchanges per atom or cluster radius change by two interatomic distances. The 1D migration distance ranged up to several nanometers. We compared the frequency and distance of 1D migration with those for in situ observation using high-voltage electron microscopy, allowing for the extremely large rate of interatomic mixing and cluster shrinkage in the present simulation.     

Application of a solvation shell analysis to the diffusivity of an interstitial species in binary substitutional solutions

A model is presented for the effect of substitutional solute elements on the diffusion coefficient of an interstitial species for low concentrations of both solutes. The model is based on the theory of absolute reaction rates and a Wagner random, noninteractive model for the substitutional solution, with Maxwell-Boltzmann statistics for the distribution of the interstitial atoms. Simplifying assumptions are made concerning the configuration of activated complexes, treating them as formed by the interstitial atom and the two nearest lattice atoms. The calculations are compared with experimental data reviewed for the diffusivity of carbon in Fe-Cr, Fe-Mn, Fe-Co and Fe-Ni alloys.     

Molecular-Dynamics Approach to the Magnetic Structure of Competing Magnetic Alloys: Fe-Cr Alloys

A molecular dynamics (MD) approach which determines automatically the complex magnetic structures in itinerant electron systems is applied to Fe-Cr alloys with use of 250 atoms in a MD unit cell (5×5×5 bcc lattice). It is demonstrated that the Fe-Cr alloys show various complex magnetic structures due to competing interactions: the collinear ferromagnetism (F) of matrix Fe with antiparallel Cr moments beyond 80 at.% Fe, the coexistence of non-collinear structure of Cr and collinear F of Fe between 50 and 75 at.% Fe, the coexistence of broken antiferromagnetism (AF) of Cr and the F of Fe between 25 and 45 at.% Fe, the coexistence of F of Fe and antiferromagnetic long-range order of Cr around 20 at.% Fe, the AF of Cr matrix with non-collinear Fe moments (spin-glass like structure) between 5 and 15 at.% Fe, and the AF below 5 at.% Fe. In the concentration region between 5 and 20 at.% Fe, ferromagnetic Fe pairs which are stabilized with different amplitudes of local moments are found. The magnetic phase diagram and calculated magnetic moments are shown to be consistent with the neutron, Mössbauer, and photoemission experiments.     

Clusters and magnetic anchoring points in (Ga,Fe)N condensed magnetic semiconductors

The stability and magnetic properties of Fe clusters in the (Ga,Fe)N magnetic semiconductor is investigated by using first-principles density functional theory and local spin density+Hubbard U theoretical methods. The present results reveal the existence of ferrimagnetic clusters formed by three or four peripheral Fe atoms neighboring a central Fe atom acting as a robust magnetic anchoring point. These clusters have magnetic moments 2 or 3 times that of a single Fe atom and, when connected by sharing peripheral Fe atoms, can form stable, ordered magnetic regions where all of the central atoms are ferromagnetically coupled. The formation of these ferrimagnetic clusters is proposed here to be at the origin of the ferromagnetic behavior observed in (Ga,Fe)N samples showing chemical phase separation.     

嵌入原子法计算金属钚中点缺陷的能量

钚因放射性衰变而出现老化效应.钚中点缺陷的性质和行为是理解钚老化效应的一个基础和前提.运用分子动力学模拟技术,计算了金属钚中点缺陷和点缺陷团簇的形成能和结合能.其中钚-钚、钚-氦和氦-氦相互作用势分别采用嵌入原子多体势、Morse对势和Lennard-Jones对势.计算结果表明,单个自间隙原子易以〈100〉哑铃状形态存在;间隙氦原子在理想晶格的八面体间隙位置相对较为稳定;氦原子与空位的结合能较大,在钚的自辐照过程中两者易于结合并形成氦-空位团簇;氦-空位团簇的形成能随氦原子数的增加而增大,当氦与空位的数     

SIA energetic and structural morphology of α-iron

In this paper, the stability of various atomic configurations containing a self-interstitial atom (SIA) in a model representing α-iron is investigated. From the energy panorama maps of the SIA located in possible non-equivalent interstitial sites, six relatively stable self-interstitial sites are found, whose structures and formation energies have been described and calculated using the modified analytical embedded atom method and molecular dynamics. The simulation results indicate that the [110] dumbbell interstitial is the energetically most favorable configuration, which is in good agreement with the experimental and ab initio results, and the distances between two displaced atoms that compose the [100], [110] and [111] direction dumbbells have been computed to be 0.68a, 0.65a and 0.29a, respectively, not all being about 0.75a apart. The relaxed displacements up to the fifth-nearest-neighbor atoms around the SIA in O interstitial position are also calculated.     

The Role of Intermixing in the Phase Determination of the Interlayer Exchange Coupling in Multilayers: Application to Fe-Cr Superlattices

The phase of the short-range interlayer exchange coupling oscillations in Fe n Cr m superlattices is expressed in terms of the intermixing at the Fe-Cr as well as at the Cr-Fe interfaces. The interdiffusion is modeled through stochastic algorithms, which presupposes floating of the atoms on the surface of the sample during the epitaxial growth. It automatically leads to the different chemical and magnetic structure of Fe on Cr and Cr on Fe interfaces. Self-consistent calculations of the magnetic moments are performed on the basis of the periodic Anderson model. Although short-range (2 monolayers (ML)) oscillations of exchange coupling were detected for all considered structures, its amplitude and phase strongly depend on the interface alloying. Introduction of the same intermixing at both interfaces does not change the phase as compared to the ideal superlattices with sharp interfaces. However different interdiffusion leads to the ~ -phase shift in accordance with experimental results for the Fe-Cr-Fe trilayers grown on an Fe whiskers. Distribution of magnetic moments on Fe atoms contains several distinct peaks but their position and relative area weakly depend on the alloying and interlayer exchange coupling in the superlattice. On the contrary, magnetic moments on Cr atoms are very sensitive to the intermixing and their behavior determines the strength and the phase of exchange coupling oscillations.     

The multiple hydrogen location near a α-Fe vacancy

The interaction between 10 hydrogen atoms and a α-Fe structure having a vacancy (V) has been studied using a cluster model and a semi-empirical theoretical method. The energy of the system was calculated by the atom superposition and electron delocalization molecular orbital method. The electronic structure was studied using the concept of density of states and crystal orbital overlap population curves.For the study of a sequential absorption, the hydrogen atoms were positioned in their energy minima configurations, near to the tetrahedral sites neighboring the vacancy, except the last H atom that was located far from the vacancy. The energy difference for H agglomeration was also computed. The vacancy-H_n complexes become less stable than VH species for more than three hydrogen's.The changes in the electronic structure of Fe atoms near the vacancy were also analyzed. The interactions mainly involve Fe 3d and 4s atomic orbitals. The contribution of Fe p orbitals is much less important. The Fe-Fe bond weakened as new Fe-H and H-H pairs were formed. The effect of H atoms is limited to its first Fe neighbors. The detrimental effect of H atoms on the Fe-Fe bonds can be related to the mechanism for embrittlement in α-Fe.     

H2在AlnCr(n=1-7)团簇上吸附和解离的密度泛函研究

采用用密度泛函理论中的B3LYP方法研究了H2在AlnCr(n=1-7)团簇上的吸附和解离。结果表明：AlnCr团簇结构与Aln+1团簇结构相似;物理吸附是H2是以侧向的形式吸附在Cr原子上,H-H键长略微增长,H2的振动频率发生了红移;除了n=5外,其它AlnCrH2团簇的最稳定结构均是AlnCr团簇的最稳定结构与两个氢原子成键而成;AlnCr团簇向H原子转移了电荷;AlnCrH2团簇的平均结合能,垂直电离势和能隙均大于AlnCr团簇的,即AlnCrH2团簇比AlnCr团簇更稳定;Al7Cr对H2的化学吸附表现出较强的惰性,而AlnCrH2(n=1,2,6)则表现出较强的化学活性;由化学反应路径跟踪可知,通过改变AlnCr团簇中Al原子的个数可以调节H2的物理化学吸附行为。     

H_2在Al_nCr(n=1-7)团簇上吸附和解离的密度泛函研究

采用密度泛函理论中的B3LYP方法研究了H2在AlnCr(n=1-7)团簇上的吸附和解离.结果表明:AlnCr团簇结构与Aln+1团簇结构相似;物理吸附是H2以侧向的形式吸附在Cr原子上,H-H键长略微增长,H2的振动频率发生了红移;除了n=5外,其它AlnCr H2团簇的最稳定结构均是AlnCr团簇的最稳定结构与两个氢原子成键而成;AlnCr团簇向H原子转移了电荷;AlnCr H2团簇的平均结合能,垂直电离势和能隙均大于AlnCr团簇的,即AlnCr H2团簇比AlnCr团簇更稳定;Al7Cr对H2的化学吸附表现出较强的惰性,而AlnCr H2(n=1,2,6)则表现出较强的化学活性;由化学反应路径跟踪可知,通过改变AlnCr团簇中Al原子的个数可以调节H2的物理化学吸附行为.     

非晶态合金Fe80P20局域结构和性质的量子化学研究

利用原子簇模型Fe_4P简化了非晶态合金Fe_(80)P_(20)的局域结构,设计了四方锥、三角双锥、四面体及平面五边形等十几种构型,对其二、四重态分别进行密度泛函(DFT)优化计算,经过频率验证,获得五种稳定构型。从所得优化构型的键长和键级,可以发现原子簇Fe4P较好地反映了非晶态合金Fe_(80)P_(20)的局域结构。考察了各构型间的过渡转化情况,发现二重态构型的稳定性要好于四重态。分析各构型的能量、成键及电子转移情况,发现与P原子成键的Fe原子个数对这些性质影响较为明显。与P原子成键的Fe原子个数越多,体系的能量就越低,越容易存在；P原子的得电子能力随着与其成键Fe原子个数增多而减少,甚至会将失去自身电子转移到金属原子上。同时通过3d轨道布居数,讨论了原子簇的空穴数及磁学性质。     

铁铬合金团簇点阵中与成分相关的磁性耦合

本文采用气相团簇束流沉积法制备了不同铁铬比的铁铬合金密集团簇点阵,研究了团簇点阵中复杂的多相结构和各种耦合效应.当合金中铬含量较大时,在团簇中能够观测到一种晶格失配的类四方结构,这种结构的出现导致了不同铁铬比的合金团簇中交换偏置效应的不同.随着铬的含量增加,合金团簇的交换偏置场减小,而团簇间的偶极相互作用增强.在合金团簇的铁磁-超顺磁转变温度以上,能明显观察到在类四方结构团簇中有更大的残存矫顽力与剩磁.在场冷条件下,对矫顽力和偏置场随温度变化的研究表明含有更多类四方结构的合金团簇有更好的热稳定性和更大的各向异性.     

α-铁中位移级联过程中缺陷的生成和He-空位复合物的形成

运用分子动力学方法研究了不同He浓度和不同级联能下含He的α-Fe低温时的位移级联过程.模拟的主要撞击原子(PKA)的能量(E_P)(即级联能)从500eV变化到5keV，辐照温 度为100K，He的浓度从1%变化到5%.比较了不同He浓度下的Fe-He混合物和纯α-Fe的位移级 联过程，发现在含有He的α-Fe中总的Fe和He的Frenkel对与纯的α-Fe的Frenkel对的数目相 当.当He的浓度较低时，含有He的α-Fe中的Frenkel对比α-Fe的要低，随He的浓度增加，有 He的α-Fe中的Frenkel对比α-Fe中的高，这主要与He在金属中的性质有关.本研究证明了位 移级联过程可以直接导致He泡的成核.对不同级联能和不同He浓度下Fe的位移级联过程的模 拟，发现在同样的级联能下，随着He浓度的增加，He成团的几率增大；在同样的He浓度下, 随着级联能的增加,He成团的几率也同样增加，并分析了级联下He泡的形成机制. 关键词： Fe 位移级联 He-空位复合物 分子动力学     

Mössbauer spectroscopy of primary knocked-on atoms

In many cases behavior of irradiated alloys is determined by the formation and evolution of the primary damaged zones in the alloys. There are several mechanisms responsible for these processes. One of them is the atomic separation of the alloy species in the cascade zone which can be directly examined by Mössbauers spectroscopy of the primary knocked-on atoms. In this report we summarize some results of measurements on amorphous alloys and intermetallic compounds.