辐照损伤模拟中的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.     

质子对BaTiO3薄膜辐照损伤的计算机模拟

本文先应用分子动力学模拟BaTiO₃体系在初级击出原子(primary knock-on atom, PKA)轰击下缺陷产生和复合的动力学过程, 模拟结果表明：PKA的方向和能量对缺陷数目有重要影响, 并计算了Ba, O和Ti原子的平均位移阈能分别为69 eV, 51 eV和123 eV, 远大于SRIM程序默认的位移阈能25 eV. 然后应用蒙特卡罗软件包SRIM, 模拟质子在BaTiO₃薄膜中的能量损失过程, 比较位移阈能对模拟结果的影响, 分析质子能量和入射角度对空位数量以及分布的影响. 结果表明空位数量随着质子能量增加而增加, 增加的速率随能量的增加是降低的;当入射角度大于60°, 空位数量随入射角增大而明显减少.     

Atomistic simulations to characterize the influence of applied strain and PKA energy on radiation damage evolution in pure aluminum

Knowledge of defects generation, their mobility, growth rate, and spatial distribution is the cornerstone for understanding the surface and structural evolution of a material used under irradiation conditions. In this study, molecular dynamics simulations were used to investigate the coupled effect of primary knock-on atom (PKA) energy and applied strain (uniaxial and hydrostatic) fields on primary radiation damage evolution in pure aluminum. Cascade damage simulations were carried out for PKA energy ranging between 1 and 20 keV and for applied strain values ranging between ?2% and 2% at the fixed temperature of 300 K. Simulation results showed that as the atomic displacement cascade proceeds under uniaxial and hydrostatic strains, the peak and surviving number of Frenkel point defects increases with increasing tension; however, these increments were more prominent under larger volume changing deformations (hydrostatic strain). The percentage fraction of point defects that aggregate into clusters increases under tension conditions; compared to the reference conditions with no strain, these increases are around 13% and 7% for interstitials and vacancies, respectively (under 2% uniaxial strain), and 19% and 11% for interstitials and vacancies, respectively (under 2% hydrostatic strain). Clusters formed of vacancies and interstitials were both larger under tensile strain conditions, with increases in both the average and maximum cluster sizes. The rate of increase/decrease in the number of Frenkel pairs, their clustering, and their size distributions under expansion/compression strain conditions were higher for higher PKA energies. Overall, the present results suggest that strain effects should be considered carefully in radiation damage environments, specifically for conditions of low temperature and high radiation energy. Compressive strain conditions could be beneficial for materials used in nuclear reactor power systems.     

稀土La对bcc-Fe中Cu扩散行为影响的第一性原理研究

采用基于密度泛函理论的第一性原理,研究了稀土La对bcc-Fe中Cu析出行为的影响.计算了α-Fe中La原子和Cu原子与空位之间,以及La原子和Cu原子之间的点缺陷结合能;在此基础上,讨论了α-Fe中La对Cu扩散激活能的关系.结果表明:La原子与空位之间有较强的相互吸引作用,且对近邻Cu原子也有一定的束缚.此外,La的加入使Cu原子近邻的空位形成能显著升高,这表明La,Cu偏聚区形成空位较为困难.与此同时,由于La原子对近邻空位和Cu原子的吸引作用,使Cu原子向近邻空位跳跃的迁移能有所升高.迁移能与空位形成能变化的计算结果显示,La原子的加入能够使α-Fe中Cu的扩散激活能显著升高,从而延缓了铁素体区富铜相的偏聚和析出.     

Radiation loss of electrons under scattering by a Thomas-Fermi atom

A universal theory and calculation results for the bremsstrahlung of electrons on complex atoms are presented. The theory accounts for the dynamic polarization of the core in the energy range from 0.5 to 10 keV, which is characteristic of radiation energy losses in a hot plasma with heavy ions. The treatment is based on the statistical atom model and the quasi-classical approximation of the incident electron. The model accounts for the penetration of the incident electron into the atomic core, which affects the relationship between the polarization and static radiation channels. The contribution of the polarization channel in both the spectral and the total radiation loss of electrons at various frequencies, nucleus charges, and energies of the incident particle is analyzed. It is shown that the contribution of the polarization channel is comparable with that of the static channel (which was calculated elsewhere) in a wide range of parameters. The results obtained are in a reasonable quantitative agreement with the detailed quantum-mechanical calculation carried out for individual atoms.     

B,N协同掺杂金刚石电子结构和光学性质的第一性原理研究

基于金刚石的稳定结构, 在实验研究的基础上, 本文采用基于周期性密度泛函理论计算了B/N单掺杂、共掺杂金刚石的晶体结构, 并就掺杂方式和掺杂后形成能进行了对比研究, 得到了B/N双掺杂的最稳定结构. 在此基础上, 进一步计算了N单掺杂及B/N共掺杂最稳定结构的吸收光谱、电子结构和态密度. 通过与实验结果对比可以看出, 较N单掺杂, B/N共掺杂的吸收光谱发生明显红移, 与实验符合较好. 计算结果表明: N原子单掺杂优先于B原子; 由于原子间的协同作用, B/N近邻共掺杂体系的形成能最低, 为掺杂的最可能结构.     

Energetics of small hydrogen-vacancy clusters in bcc iron

Hydrogen may be trapped in voids in iron, leading to undesirable material properties. In this paper, the energetics of small hydrogen-vacancy clusters in body centered cubic iron are investigated. Results from two interatomic potentials are compared. We use molecular dynamics and Monte Carlo methods to find the minimum energy configurations of voids of up to ten vacancies containing up to 50 hydrogen atoms with ratios of hydrogen to vacancy of 10 or less. The formation energies and binding energies of defects to these clusters are calculated. Our results indicate that the hydrogen stabilizes bubbles by causing vacancies to be more tightly bound to clusters, while neighboring irons are less tightly bound. Hydrogen itself becomes less well bound to clusters as the inventory increases. The more physically relevant potential indicates a maximum supported ratio of hydrogen atoms to vacancies of about 4.     

Atom abstraction in the scattering of state-selected NO+(X1Sigma(+)) on O/Al(111)

The reactive scattering of state-selected NO+(X1Sigma(+)) on oxygen-covered Al(111) is explored at hyperthermal collision energies. Relative ion yields and mean translational energies of scattered NO-2 are presented as a function of oxygen exposure and NO+ collision energy. Above the 9+/-1 eV threshold for reaction, NO-2 products emerge with an average kinetic energy that depends linearly on incident NO+ energy. The formation of NO-2 is assigned to the direct, Eley-Rideal abstraction of an adsorbed O atom by an incident NO molecule.     

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

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

Point defects in L1_0 FePt studied by molecular dynamics simulations based on an analytic bond-order potential

The point defects and their related physical properties in L10 FePt are investigated by molecular dynamics simulations based on an analytic bond-order potential. The calculated results agree well with the experimental value, indicating that the analytic bond-order potential is suitable to describe the structural properties and surface energies of the FePt alloy in the L10 phase. However, the calculated vacancy formation energy of an Fe atom is higher than that of a Pt atom, which disagrees with some other p...     

Atomistic structure of stepped surfaces

Atomistic computer simulation with embedded atom method (EAM) interatomic forces was used to study the structure of surface steps on the {111} unreconstructed surface in fcc metallic materials. The energetics and local atomic relaxation behavior of ledges parallel to the 110 direction were studied using a potential describing lattice properties of Au. The vacancy formation energies in the stepped surfaces was also studied, and it was found that the energy of formation of a vacancy in a terrace is the same as that in the perfect unstepped surface. This value is 30% lower than that of the bulk. The vacancy formation energy in the ledge is reduced by a factor of two with respect to that of the terraces. The structure of the “up ledge” (A step) is different from the “down ledge” (B step). These differences do not significantly affect the energy of the ledges, although they do affect the vacancy formation energies in sites in the second surface layer near the ledge. The implications of the results for the formation of kinks and the general structure of high index surfaces are discussed.     

Energy dependence of site‐specific radiation damage in protein crystals

It is important to consider radiation damage to crystals caused by data collection when solving structures and critical when determining protein function, which can often depend on very subtle structural characteristics. In this study the rate of damage to specific sites in protein crystals cooled at 100 K is found to depend on the energy of the incident X‐ray beam. Several lysozyme crystals were each subjected to 3–26 MGy of cumulative X‐ray exposure by collecting multiple data sets from each crystal at either 9 keV or 14 keV. The integrated electron density surrounding each S atom in the structure was calculated for each data set and the change in electron density was evaluated as a function of dose at the two energies. The rate of electron density decrease per cubic Å per MGy was determined to be greater at 14 keV than at 9 keV for cysteine sulfurs involved in disulphide bridges; no statistically significant differences in the decay rates were found for methionine sulfurs. These preliminary results imply that it might be possible to minimize certain types of specific radiation damage by an appropriate choice of energy. Further experiments studying a variety of photolabile sites over a wider range of energies are needed to confirm this conclusion.     

Molecular dynamics simulation of displacement damage in 6H-SiC

ABSTRACT

As a new generation of low-loss components, 6H-SiC is widely used in optoelectronic devices, electronic devices and other fields, especially in high temperature and strong radiation environment. Radiation will lead to a great large of defects generated in the material, then forming permanent displacement damage which will result in performance down or failure. In this paper, the molecular dynamics method was used to study the spatial distribution of defects and the effect of temperature and PKA energy on the various properties of the point defects. The main contributions were the evolution of defects in 6H-SiC crystal, the influence of PKA energy and temperature on the defect distribution and the number of point defect. In this paper, the spatial distribution of point defects in 6H-SiC crystal after PKA incidence was obtained, the recombination rate under four kinds of energy PKA was calculated, and the effect of temperature and incident PKA energy on the number of defects, the rate of the vacancy cluster and the rate of vacancy defects at steady state were investigated. The results show that the number of defects at steady-state increases linearly with the increase of incident PKA energy, the effect of temperature and energy on defects producing and the rate of clusters are insignificant.     

Hydrogen chemisorption on metal surfaces

The adsorption of H atoms on metal (jellium) surfaces has been investigated using linear response theory within the density functional formalism. The adsorbate is represented initially by a 1S orbital on the H atom, which perturbs the jellium surface and indirectly the adsorbate itself. The interaction energy curves, atomic binding energies, induced dipole moments, chemical shifts associated with the adsorbate, and vibrational excitation energies at the equilibrium internuclear separation have been calculated for a single H atom chemisorbed on metal surfaces. The sum of the atomic binding energy and the ionization potential of the H 1S level may be regarded as the initial state energy in the case of photoemission from the chemisorbed H. The rather satisfactory overall agreement between the theory and the experimental results for binding energies, vibrational excitation energies, and dipole moments suggests that this simple formalism could also have useful applications in more complicated chemisorbed systems.     

外延压应变对BaTiO$lt;sub$gt;3$lt;/sub$gt;铁电体抗辐射性能影响的分子动力学研究

采用基于壳模型的分子动力学模拟方法, 研究了存在外延压应变时BaTiO₃铁电体的辐射位移效应, 以O原子作为初冲原子(primary knock-on atom, PKA), 能量为1 keV, 方向为[001], 分别计算了外延压应变为0, 0.4%, 0.8%, 1.2%, 1.6%, 2.0%时体系的缺陷数量、分布, 以及辐射前后的极化强度, 比较了压应变为2%以及无应变下损伤区域、缺陷离位距离和反向外电场下PKA的迁移距离. 结果表明, 随外延压应变增加体系极化近似线性增加, 辐射后极化降低幅度降低、缺陷产生的数量有所减小, 2% 压应变存在时缺陷原子的离位距离、PKA在反向外电场作用下的迁移距离和损伤区域都小于无应变的情况, 说明外延压应变的存在对辐射造成的晶格损伤具有抑制作用, 对辐射损伤具有改善作用, 可以通过引入外延压应变来调控BaTiO₃的辐射损伤. 关键词： 应变 3')" href="#">BaTiO₃ 辐射损伤 分子动力学模拟     

The kinetic MC modelling of reversible pattern formation in initial stages of thin metallic film growth on crystalline substrates

The results of kinetic MC simulations of the reversible pattern formation during the adsorption of mobile metal atoms on crystalline substrates are discussed. Pattern formation, simulated for submonolayer metal coverage, is characterized in terms of the joint correlation functions for a spatial distribution of adsorbed atoms. A wide range of situations, from the almost irreversible to strongly reversible regimes, is simulated. We demonstrate that the patterns obtained are defined by a key dimensionless parameter: the ratio of the mutual attraction energy between atoms to the substrate temperature. Our ab initio calculations for the nearest Ag-Ag adsorbate atom interaction on an MgO substrate give an attraction energy as large as 1.6 eV, close to that in a free molecule. This is in contrast to the small Ag adhesion and migration energies (0.23 and 0.05 eV, respectively) on a defect-free MgO substrate.     

On possibility of measurement of the electron beam energy using absorption of radiation by electrons in magnetic field

A possibility of precise measurement of the electron beam energy using absorption of radiation by electrons in a homogeneous magnetic field for electrons of high energy in the range up to a few hundred GeV, was considered earlier. In this paper, with the purpose of experimental checking of this method in the range of several tens MeV of electrons energies, a possibility of measurement of absolute energy of the electron beam with a relative accuracy up to 10^?4, is considered. We take into account influence of the laser beam diffraction, of the spread of electrons over energies, and of the length of formation of radiation absorption in the process of electron beam energy measurement. The laser wavelength and the length of the magnet are chosen depending on the length of photon absorption formation. It is found that the kinematical restrictions on the photon absorption process lead to the selection in angles of propagation of photons, which can be absorbed by the beam electrons. It is shown that parameters of the electron beam will noticeably not vary during the measurement of the energy.     

Be原子在Be基底上的沉积过程研究

利用分子动力学方法模拟了Be原子在Be基底上的沉积过程. 模拟了沉积粒子不同入射动能条件下, 沉积薄膜表面形态的差异. 在一定能量范围内, 增加粒子入射动能可以减小薄膜的表面粗糙度. 但是, 过高的入射动能, 不利于减小薄膜表面粗糙度. 通过沉积薄膜中原子配位数以及单个原子势能沿薄膜厚度的分布, 分析沉积原子入射动能对于薄膜及表面结构的影响. 沉积动能较大时, 薄膜的密度较大; 单个原子势能沿薄膜厚度分布较为连续; 同时薄膜中原子应力沿薄膜厚度分布较为连续. 最后, 分析了沉积粒子能量转化的过程、粒子初始动能对基底表面附近粒子局部动能增加的影响.     

Molecular dynamics study of structural damage in amorphous silica induced by swift heavy-ion radiation

In this paper, the radiation defects induced by the swift heavy ions and the recoil atoms in amorphous SiO₂ were studied. The energy of recoil atoms induced by the incident Au ions in SiO₂ was calculated by using Monte Carlo method. Results show that the average energies of recoils reach the maximum (200?eV for Si and 130?eV for O, respectively) when the incident energy of Au ion is 100?MeV. Using Tersoff/zbl potential with the newly built parameters, the defects formation processes in SiO₂ induced by the recoils were studied by using molecular dynamics method. The displacement threshold energies (E_d) for Si and O atoms are found to be 33.5 and 16.3?eV, respectively. Several types of under- and over-coordinated Si and O defects were analyzed. The results demonstrate that Si₃, Si₅, and O₁ are the mainly defects in SiO₂ after radiation. Besides, the size of cylindrical damage region produced by a single recoil atom was calculated. The calculation shows that the depth and the radius are up to 2.0 and 1.4?nm when the energy of recoils is 200?eV. Finally, it is estimated that the Au ion would induce a defected track with a diameter of 4?nm in SiO₂.     

Al掺杂对HfO2俘获层可靠性影响第一性原理研究

采用基于MS(Materials Studio)软件和密度泛函理论的第一性原理方法, 研究了HfO₂ 俘获层的电荷俘获式存储器(Charge Trapping Memory, CTM)中电荷的保持特性以及耐擦写性. 在对单斜晶HfO₂中四配位氧空位(V_O4) 缺陷和V_O4 与Al替位Hf掺杂的共存缺陷体(Al+V_O4)两种超晶胞模型进行优化之后, 分别计算了其相互作用能、形成能、Bader电荷、态密度以及缺陷俘获能. 相互作用能和形成能的计算结果表明共存缺陷体中当两种缺陷之间的距离为2.216 Å时, 结构最稳定、缺陷最容易形成; 俘获能计算结果表明, 共存缺陷体为双性俘获, 且与V_O4缺陷相比, 俘获能显著增大; Bader电荷分析表明共存缺陷体更有利于电荷保持; 态密度的结果说明共存缺陷体对空穴的局域能影响较强; 计算两种模型擦写电子前后的能量变化表明共存缺陷体的耐擦写性明显得到了改善. 因此在HfO₂俘获层中可以通过加入Al杂质来改善存储器的保持特性和耐擦写性. 本文的研究可为改善CTM数据保持特性和耐擦写性提供一定的理论指导.