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

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

金属钨中氦行为的分子动力学模拟

采用分子动力学方法模拟了氦在金属钨中的扩散聚集行为. 首先,建立了氦与钨原子间相互作用势,短程部分采用ZBL势形式,长程部分采用从头算法数据,实现了两者之间的平滑连接. 通过计算氦在钨中不同间隙位的形成能发现,单个氦原子更易存在于金属钨中的四面体间隙位,这与最新的研究成果是一致的. 在400-1200 K的温度范围内,考察了氦原子在金属钨中的扩散行为,获得了扩散迁移能,其值介于实验值和从头算法结果之间. 最后,研究了氦的聚集行为,从能量的角度考察了氦团簇形成初期的生长机理. 研究发现,在氦团簇形成初期,氦团簇对氦的结合能随着氦团簇的生长有逐渐增大的趋势,说明氦团簇吸收氦的能力逐渐增强. 关键词： 氦扩散 氦团簇 辐照损伤 分子动力学模拟     

bcc Fe的刃型位错中氦-空位团的稳定性

 采用周期性原子阵列方法建立bcc Fe中的刃型位错,利用分子动力学计算了0 K时bcc Fe的位错芯里氦-空位团的稳定性,并与理想Fe晶体里氦-空位团的稳定性进行比较发现,位错的作用导致氦-空位团不稳定。点缺陷（He、空位与自间隙Fe原子）与氦-空位团的结合能与团中氦-空位比例密切相关,当氦与空位数之比在3~6时,结合能趋于稳定。     

金属钛中氦团簇生长行为的分子动力学研究

运用分子动力学方法模拟了常温下金属钛中氦团簇的生长过程.从能量的角度考察了氦团簇的生长机理.研究发现,随着氦团簇的生长,在氦团簇周围逐渐形成位错环缺陷,氦团簇与氦原子的结合能有逐渐下降的趋势,当氦团簇生长到一定尺寸时会通过发射周围缺陷以使得结合能上升,从而增强了进一步吸收氦原子的能力.模拟还发现,随着氦团簇的继续生长,氦团簇的形状由原来的不规则结构逐渐变成了较为规则的棱柱形结构,在随后的生长过程中其生长仅在（001）平面进行,沿［001］轴的厚度几乎不变. 关键词： 氦团簇 缺陷发射 分子动力学模拟     

Al_mNi_n(m+n=2-4)二元合金团簇结构和稳定性的密度泛函研究

在密度泛函理论的B3PW91水平上,优化计算了Al_mNi_n(m+n=2-4)团簇的几何结构、团簇基态的结合能、最高占据轨道与最低空轨道之间的能级间隙和可能的分离能.结果表明,随着团簇尺寸的增加,能级间隙呈现奇偶振荡规律;Al_2Ni_2团簇具有最高的结合能,易分离成Ni原子和Al_2Ni团簇;含Al的合金团簇最易先分离出一个Al原子,而含Ni的合金团簇没有这种趋势.     

压力对金属铌中氦原子凝聚的影响

 采用零温条件下的赝势-平面波方法和有限温度下的Car-Parrinello分子动力学方法,模拟了不同压力环境下氦原子在金属铌中的行为特征,研究了宿主缺陷和氦泡的形成机制。结果表明,闭电子壳层的氦原子在金属铌中具有刚球模型特征,其占据区域为金属自由电子的禁区,从而破坏铌原子之间的金属性键合。在常温条件下,局域高浓度的氦原子优先凝聚于近邻宿主空位缺陷处,从而形成氦泡;完整晶格中高浓度的氦将促使铌原子易位,形成间隙-空位模式的宿主缺陷,氦原子聚集于空位区域。完整宿主在压力（40 GPa）的作用下,晶格参数减小,铌原子之间的相互作用增强,尽管氦原子的存在削弱了铌原子之间的相互作用,位于格点上的铌原子仍难以借助热振动偏离格点形成空位,因而未能形成间隙-空位对和氦泡。     

AlmNin (m+n=2–4)二元合金团簇结构和稳定性的密度泛函研究

在密度泛函理论的B3PW91水平上, 优化计算了AlmNin (m+n=2–4)团簇的几何结构、团簇基态的结合能、最高占据轨道与最低空轨道之间的能级间隙和可能的分离能. 结果表明, 随着团簇尺寸的增加, 能级间隙呈现奇偶振荡规律; Al2Ni2团簇具有最高的结合能, 易分离成Ni原子和Al2Ni团簇; 含Al的合金团簇最易先分离出一个Al原子, 而含Ni的合金团簇没有这种趋势.     

Fe-Cr合金辐照空洞微结构演化的相场法模拟

Fe-Cr合金作为包壳材料在高温高辐照强度等极端环境下服役,产生空位和间隙原子等辐照缺陷,辐照缺陷簇聚诱发空洞、位错环等缺陷团簇,引起辐照肿胀、晶格畸变,导致辐照硬化或软化致使材料失效.理解辐照缺陷簇聚和长大过程的组织演化,能更有效调控组织获得稳定服役性能.本文采用相场法研究Fe-Cr合金中空洞的演化,模型考虑了温度效应对点缺陷的影响以及空位和间隙的产生和复合.选择400—800 K温度区间、0—16 dpa辐照剂量范围的Fe-Cr体系为对象,研究在不同服役温度和辐照剂量下的空位扩散、复合和簇聚形成空洞的过程.在400—800 K温度区间,随着温度的升高,Fe-Cr合金空洞团簇形核率呈现出先升高后下降的趋势.考虑空位与间隙的重新组合受温度的影响可以很好地解释空洞率随温度变化时出现先升高后降低的现象.由于温度的变化将影响Fe-Cr合金中原子离位阀能,从而影响产生空位和间隙原子.同一温度下,空洞半径和空洞的体积分数随辐照剂量的增大而增大.辐照剂量的增大,级联碰撞反应加强,空位与间隙原子大量产生,高温下空位迅速的扩散聚集在Fe-Cr合金中将形成更多数量以及更大尺寸的空洞.     

原子模拟钛中微孔洞的结构及其失效行为

六角金属由于其各向异性等特点,在塑性变形等过程中容易产生形状和构型都相对复杂的点缺陷团簇.这些团簇之间及其与运动位错等缺陷的相互作用直接影响材料的物理和力学性能.然而对相关问题的原子尺度、尤其是空位团簇的演化和微孔洞的形成乃至裂纹形核扩展等的理解还不全面.本文采用激发弛豫算法结合第一原理及原子间作用势,系统考察了钛中的空位团簇构型及不同构型间的相互转变,给出了不同尺寸空位团簇的稳定和亚稳构型、空位团簇合并分解和迁移的激发能垒等关键参数,发现较小的空位团簇形成稳定构型,较大的空位团簇呈现出空间对称分布趋势进而形成微孔洞;采用高通量分子动力学模拟系统研究了不同尺寸的空位团簇在拉应力作用下对变形过程的影响,发现这些空位团簇可以形成层错,并对微裂纹的形核产生影响.     

Al_nCl(n=2-14)团簇结构与电子性质的理论研究

运用卡里普索(CALYPSO)结构预测方法,在杂化密度泛函B3LYP/6-311G+(d)基组水平上,对Al_nCl(n=2-14)团簇的几何结构与电子性质进行优化计算,并讨论了团簇的平均结合能、能隙、二阶能量差分、电离能、亲和能以及电子自然布居和极化率.研究结果表明:Al_nCl(n=2-14)团簇的基态构型由简单平面几何结构向立体结构演化,形成Cl原子戴帽Al_n-1Cl团簇结构;Cl原子的掺杂增大了Al_n团簇的平均结合能;二阶能量差分、能隙、电离能、亲和能的变化表明Al_7Cl是幻数团簇结构;团簇中的电荷总是由Al_原子向Cl原子转移,原子之间的成键作用随着团簇尺寸的增大而增强.     

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.     

Effect of grain boundary structures on the behavior of He defects in Ni: An atomistic study

We investigated the effect of grain boundary structures on the trapping strength of He_N(N is the number of helium atoms) defects in the grain boundaries of nickel. The results suggest that the binding energy of an interstitial helium atom to the grain boundary plane is the strongest among all sites around the plane. The He_N defect is much more stable in nickel bulk than in the grain boundary plane. Besides, the binding energy of an interstitial helium atom to a vacancy is stronger than that to a grain boundary plane. The binding strength between the grain boundary and the He _N defect increases with the defect size. Moreover, the binding strength of the He_N defect to the Σ3(12)[110] grain boundary becomes much weaker than that to other grain boundaries as the defect size increases.     

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.     

First principle study of the behavior of helium in plutonium dioxide

The incorporation and solution of helium in plutonium dioxide have been investigated based on density functional theory. The GGA and GGA + U approximations were used with the projector-augmented-wave method. Several defects that are likely to accommodate the incorporation of helium in PuO₂, such as oxygen vacancy, plutonium vacancy, divacancy and Schotty defects were considered in this work. With GGA approach, the lowest incorporation energy corresponds to neutral trivacancy, followed by divacancy and plutonium vacancy, while the GGA + U scheme gave us that oxygen vacancy is the most favorable incorporation site for He. Both SP-GGA and SP-GGA + U methods obtained a same conclusion that the most favorable solution site for He is oxygen vacancy, interstitial site and plutonium vacancy for under-, perfect and over-stoichiometry, respectively. Additionally, the concentrations of the point defects and the solution energy of He for the different incorporation sites as a function of the stoichiometry were also obtained based on the point-defect model.     

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.     

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.     

Simulation of nanostructure evolution under helium implantation in Fe–(2.5–12.5)at% Cr alloys at a temperature of 343 K

Characterization of helium-implanted Fe–(2.5–12.5)at% Cr alloys with the flux of 7?×?10^–6?dpa/s at a temperature of 343?K has been performed by means of cluster dynamics simulations. We have suggested a model for simulating an Fe–C–Cr system under helium implantation based on a selection of the latest data from atomistic studies and available experiments. Kinetics of carbon-vacancy and helium-vacancy complexes has been studied. Only one parameter is used to guarantee the best reproduction possible of experimental positron annihilation lifetime spectroscopy data for Fe–Cr alloys on dependences of vacancy cluster size on chromium content and irradiation dose via fitting. This is an effective binding energy of self-interstitial atoms to dislocation loops decorated by chromium atoms. It has a “snaky” dependence of chromium content with a minimum of about 9%Cr.     

铝中氦原子行为的密度泛函研究

用密度泛函理论计算了大量He原子存在时He在金属铝中不同位置的能量,并在理论上预测了铝中的氦原子行为。结果表明,铝晶胞内He原子择优占位区是空位,而在整个晶体范围,最有利于容纳He原子的区域是晶界,其次是空位和位错。在fcc-铝的两种间隙位中,He原子优先充填四面体间隙位。间隙He原子的迁移能很小,易于通过迁移在晶内聚集,或被空位、晶界、位错等缺陷束缚。