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

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

Ab initio study of formation, migration and binding properties of helium-vacancy clusters in aluminum

Ab initio calculations based on density functional theory have been performed to study the dissolution and migration of helium, and the stability of small helium-vacancy clusters He_nV_m (n, m=0-4) in aluminum. The results indicate that the octahedral configuration is more stable than the tetrahedral. Interstitial helium atoms are predicted to have attractive interactions and jump between two octahedral sites via an intermediate tetrahedral site with low migration energy. The binding energies of an interstitial He atom and an isolated vacancy to a He_nV_m cluster are also obtained from the calculated formation energies of the clusters. We find that the di- and tri-vacancy clusters are not stable, but He atoms can increase the stability of vacancy clusters.     

He原子掺杂铝材料的第一性原理研究

采用基于第一性原理的平面波超软赝势方法,结合广义梯度近似(GGA),计算了铝及铝晶胞间隙位置掺入He原子后体系的几何结构、电子结构、总体能量和电荷布居值.计算结果表明:随着氦在金属铝中逐渐形成,铝晶胞体系会发生晶格畸变,但总的趋势是He在铝体系的八面体位置的晶格畸变小于其在四面体位置的晶格畸变.He在铝晶胞八面体和四面体间隙的杂质形成能分别为1.3367 eV和2.4411 eV.由此可知,He在铝晶胞中最稳定位置是八面体间隙位置.同时,文中还从原子尺度层面分析了He原子在铝晶胞中的占位及其键合性质,讨论 关键词： 铝材料 第一性原理 形成能     

Formation energies and swelling of uranium dioxide by point defects

Hybrid density functional theory is used to study the stability and behavior of rare gases in uranium dioxide. Three insertion sites are considered: the octahedral interstitial position and the oxygen and uranium substitution sites. The optimized lattice constant, the volume variation induced by gaseous atom incorporation, and the defect formation energy are studied for each rare gas. Both lattice constants and formation energies increase with increase in radii of the rare gases. The octahedral interstitial position is the most favorable occupation site. The formation energy is found to be negative only for He at an interstitial site.     

Trapping of helium atom by vacancy in tungsten: a density functional theory study

The interaction between helium (He) atom and vacancy defect in tungsten (W) has been investigated by using first-principles simulations. We have obtained that the most stable site for He in tungsten is the substitutional position because He can keep its own electronic structure at this position. In the studied tungsten system, vacancy can act as a trapping center for surrounding He atom with negative trapping energy. The migration behaviors of He atom at tetrahedral interstitial site in W, which can be trapped by vacancy but the final position is almost unchanged comparing with its initial position through structural relaxation, have been predicted and discussed. It is also found that single He atom prefers to go through an octahedral site rather than through a direct path to the vacancy, and the stronger the interaction between He atom and vacancy is, the lower the migration barrier will be.     

Influence of impurity hydrogen on the structure and properties of bulk Li and pressure effects

The structure and properties of a 16-atom body-centered cubic lithium cell with an interstitial hydrogen atom are studied using a pseudopotential-plane-wave method within the density functional theory at 0 K and high pressures. The host lattice is dramatically distorted by the introduction of H. Although the hydrogen atom is stable at the tetragonal site in perfect bcc host lattice, it favors the octahedral site formed by six non-equivalent Li atoms after full relaxation of the cell, showing P4/mmm symmetry within the pressures ranging from 0 to 6 GPa. The lattice ratio (a/c) changes irregularly with external pressure at about 3 GPa. The hydrogen band lies in the bottom of the valence band, separated by a gap from the metallic bands, illustrating the electronegativity of hydrogen. High reflectivity in the low frequency area induced by the impurity hydrogen is observed when only interband transitions are taken account of. A dip in reflectivity due to parallel band transitions is observed at ∼0.4 eV. Another dip at ∼4.3 eV appears when external pressure increases over 4 GPa.     

Atomic positions and diffusion paths of h and he in the α-Ti lattice

The solution energy of H and He in various interstitial and substitution positions in the hcp lattice of α-Ti has been calculated based on the method of electron density functional. The lowest solution energy of He corresponds to the basal octahedral position and that of H corresponds to the octahedral position (next in energy is the tetrahedral position). The calculated vibration frequencies of H in various positions are used for identification of lines in the vibration spectrum obtained by the method of neutron inelastic scattering. Taking into account these spectra, it can be concluded that hydrogen atoms occupy in the hcp lattice of Ti both the octahedral and tetrahedral positions even at 600 K. The available experimental data do not contradict the conclusion that the octahedral position is more preferable in α-Ti. The energy barriers are estimated for various diffusion paths of H and He.     

H,O,N,He原子掺杂Al原子的第一性原理研究

采用密度泛函理论框架下的第一性原理平面波赝势方法,对Al中分别加入H,O,N和He原子后的晶体状态进行了研究.通过晶体结构和形成能分析比较了杂质原子占据不同位置的难易程度及对晶体稳定性的影响,并从态密度、电荷密度和电荷布居的角度,分析了其电子结构.结果表明:H、O和N原子占据金属Al的四面体间隙最稳定,而He原子主要占据金属Al的八面体间隙. O和N原子与Al原子具有强烈的共价作用,H原子与Al原子存在共价作用但相对较弱,而He原子与Al原子的相互作用以范德华力为主.进一步揭示了四种原子在金属Al中不同行为的电子机制.     

Stability of (3)He(2)(4)He(n) and (3)He(3)(4)He(n) L=0 clusters

We have studied the stability of mixed (3)He/(4)He clusters in L=0 states by the diffusion Monte Carlo method, employing the Tang-Toennies-Yiu He-He potential. The clusters (3)He(4)He(N) and (3)He(2)(4)He(N) are stable for N>1. The lighter atoms tend to move to the surface of the cluster. The minimum number of 4He atoms able to bind three 3He atoms in a L=0 state is nine. Two of three fermionic helium atoms stay on the surface, while the third one penetrates into the cluster.     

单轴应变对H在α-Fe中占位及扩散的影响

采用基于密度泛函理论的第一性原理方法,研究了H在不同单轴应变下α-Fe中的间隙占位,计算了H原子的溶解能、态密度、电荷差分密度和电荷布居.结果表明:不同单轴拉压应变作用下,H原子优先占据四面体间隙(Ts)位,且随着压应变减小、拉应变增加,H原子越易溶于α-Fe.压应变使得Ts位的H获得更多的电子,而拉应变减少了这种电荷转移.应用LST/QST过渡态搜索计算垂直应变方向的扩散.八面体间隙位是邻近Ts位H的扩散过渡态.扩散激活能与应变呈线性关系,且随着压应变的增加,扩散激活能降低,扩散更容易.     

Energetics of He and H Atoms in W–Ta Alloys: First-Principle Calculations

Properties of various defects of He and H atoms in W-Ta alloys are investigated based on density functional theory. The tetrahedral interstitial site is the most configured site for self-interstitial He and H in W and W-Ta alloys. Only a single He atom favors a substitutional site in the presence of a nearby vacancy. However, in the coexistence of He and H atoms in the presence of the vacancy, the single H atom favors the tetrahedral interstitial site(TIS) closest to the vacancy, and the He atom takes the vacancy center. The addition of Ta can reduce the formation energy of TIS He or H defects. The substituted Ta affects the charge density distribution in the vicinity of the He atom and decreases the valence electron density of the H atoms. A strong hybridization of the H s states and the nearest W d state s exists in W_(53)He_1 H_1 structure. The sequence of the He p projected DOS at the Fermi energy level is in agreement with the order of the formation energy of the He-H pair in the systems.     

Magnetic interactions influence the properties of helium defects in iron

Density functional theory calculations of He defect properties in iron have shown an unexpected influence of magnetism arising from the defect's electronic structure. In contrast with previous work that neglected such effects, the results indicate that the tetrahedral position is energetically more favorable for the He interstitial than the octahedral site. This may have significant implications for He clustering and bubble nucleation, which will impact material performance in future fusion reactors. These results provide the basis for development of improved atomistic models.     

Vacancy trapping mechanism for multiple hydrogen and helium in beryllium: a first-principles study

The microscopic mechanism for H and He trapping by vacancy defects and bubble formation in a Be host lattice is investigated using first-principles calculations. A single He atom prefers to occupy a vacancy centre while H does not. He can segregate towards the vacancy from the interstitial site much more easily than H. Both H and He exhibit lower diffusion barriers from a remote interstitial to a vacancy with regard to their diffusion barriers inside a perfect Be solid. Up to five H or 12 He atoms can be accommodated into the monovacancy space, and the Be-He interaction is much weaker than Be-H. The physical origin for aggregation of multiple H or He atoms in a vacancy is further discussed. The strong tendency of H and He trapping at vacancies provides an explanation for why H and He bubbles were experimentally observed at vacancy defects in materials. We therefore argue that vacancies provide a primary nucleation site for bubbles of H and He gases inside Be materials.     

氦对LaNi5晶体结构的影响及其扩散特性研究

从密度泛函理论为基础的第一性原理出发,运用全势能线性缀加平面波(FLAPW)方法,对氚衰变后氦在合金中的占位以及LaNi₅He晶体结构进行了理论计算,并系统给出了氦在间隙间的迁移曲线.结果表明,氦原子在十二面体(1b)格位最稳定,并且氦从6m格位向1b格位迁移不存在势垒,而从2d格位向1b格位迁移则需越过1.55eV高的势垒.另外,氦从12n格位穿过12o格位最后到达6m格位也无明显势垒存在,并且处于4h格位之间的氦原子可以自由迁移,而相应12n格位之间的直接迁移则需跨越13.6eV高的势垒.最后还计算给出了氚衰变后合金的态密度、电子密度以及势能分布图,并与相应的LaNi₅H结构作了详细比较. 关键词： 全势能线性缀加平面波(FLAPW) 5He')" href="#">LaNi₅He 态密度 扩散     

低能CU6团簇在CU(001)表面和AU(001)表面沉积的分子动力学模拟研究

关键词：     

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.     

Theory of a surface phase transformation induced by impurity atoms in binary alloys containing interstitial impurity atoms

Mutual influence of the surface segregation of atoms on lattice sites and in interstitials is studied in the self-consistent field approximation. The effect of occurrence of a local surface phase transformation is predicted, which is caused by a local loss of stability of the statistical distribution of atoms of components on the surface due to the presence of interstitial impurity atoms. A theory is constructed for an induced phase transformation for [001] plane of a single crystal of a binary alloy with an fcc lattice whose octahedral interstitials are occupied by interstitial atoms.Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Fizika, No. 9, pp. 50–55, September, 1984.     

Molecular dynamics simulations of point defects in plutonium grain boundaries

A modified analytic embedded atom method (MAEAM) potential is constructed for fcc updelta-Pu. Molecular dynamics (MD) simulations with the potential are performed to investigate the interactions between two symmetrical tilt grain boundaries (GBs) and point defects such as He atom, vacancy and self-interstitial atom (SIA) in Pu. The calculated results show that point defect formation energies are on average lower than those in the lattice but variations from site to site along the GBs are very remarkable. Both substitutional and interstitial He atoms are trapped at GBs. Interstitial He atom is more strongly bound at the GB core than the substitutional He atom. The binding energy of SIA at GB core is higher than those of He atom and vacancy. GB core can bind many He atoms and SIAs due mainly to the fact that it contains many vacancies. Compared with He atom and SIA, the vacancy far from GB core is difficult to diffuse into the core. The GBs can act as sinks and sources of He atoms and SIAs, which may be a reason for the swelling of Pu after a period of self-irradiation because of the higher concentration of vacancy in the bulk.     

Effect of lattice deformations on the local,magnetic dipole field of interstitial sites in H.C.P. crystals

It is shown that lattice deformations around interstitial atoms may give a significant contribution to the local magnetic dipolar field at the interstitial site. An analysis of the dipole field measured for mouns on octahedral interstitial sites in Co shows that 1/4 of the dipole field results from lattice deformations. The double force tensor of the muon is found to be nearly isotropic and of the order of magnitude of 4eV.     

Self-assembly of a colloidal interstitial solid with tunable sublattice doping

We determine the phase diagram of a binary mixture of small and large hard spheres with a size ratio of 0.3 using free-energy calculations in Monte Carlo simulations. We find a stable binary fluid phase, a pure face-centered-cubic (fcc) crystal phase of the small spheres, and binary crystal structures with LS and LS(6) stoichiometries. Surprisingly, we demonstrate theoretically and experimentally the stability of a novel interstitial solid solution in binary hard-sphere mixtures, which is constructed by filling the octahedral holes of an fcc crystal of large spheres with small spheres. We find that the fraction of octahedral holes filled with a small sphere can be completely tuned from 0 to 1. Additionally, we study the hopping of the small spheres between neighboring octahedral holes, and interestingly, we find that the diffusion increases upon increasing the density of small spheres.