Supersonic N-Crowdions in a Two-Dimensional Morse Crystal

NASA Astrophysics Data System (ADS)

Dmitriev, S. V.; Korznikova, E. A.; Chetverikov, A. P.

2018-03-01

An interstitial atom placed in a close-packed atomic row of a crystal is called crowdion. Such defects are highly mobile; they can move along the row, transferring mass and energy. We generalize the concept of a classical supersonic crowdion to an N-crowdion in which not one but N atoms move simultaneously with a high velocity. Using molecular dynamics simulations for a close-packed two-dimensional Morse crystal, we show that N-crowdions transfer mass much more efficiently, because they are capable of covering large distances while having a lower total energy than that of a classical 1-crowdion.

Effects of applied strain on nanoscale self-interstitial cluster formation in BCC iron

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gao, Ning; Setyawan, Wahyu; Kurtz, Richard J.

2017-09-01

The effect of applied strains on the configurational evolution of self-interstitial clusters in BCC iron (Fe) is explored with atomistic simulations. A novel cluster configuration is discovered at low temperatures (<600 K), which consists of <110> dumbbells and <111> crowdions in a specific configuration, resulting in an immobile defect. The stability and diffusion of this cluster at higher temperatures is explored. In addition, an anisotropy distribution factor of a particular [hkl] interstitial loop within the family of loops is calculated as a function of strain. The results show that loop anisotropy is governed by the angle between the stress directionmore » and the orientation of the <111> crowdions in the loop, and directly linked to the stress induced preferred nucleation of self-interstitial atoms.« less

Effects of applied strain on nanoscale self-interstitial cluster formation in BCC iron

NASA Astrophysics Data System (ADS)

Gao, Ning; Setyawan, Wahyu; Kurtz, Richard J.; Wang, Zhiguang

2017-09-01

The effect of applied strains on the configurational evolution of self-interstitial clusters in BCC iron (Fe) is explored with atomistic simulations. A novel cluster configuration is discovered at low temperatures (<600 K), which consists of 〈 110 〉 dumbbells and 〈 111 〉 crowdions in a specific configuration, resulting in an immobile defect. The stability and diffusion of this cluster at higher temperatures is explored. In addition, an anisotropy distribution factor of a particular [ hkl ] interstitial loop within the family of 〈 hkl 〉 loops is calculated as a function of strain. The results show that loop anisotropy is governed by the angle between the stress direction and the orientation of the 〈 111 〉 crowdions in the loop, and directly linked to the stress induced preferred nucleation of self-interstitial atoms.

Suppression of radiation-induced point defects by rhenium and osmium interstitials in tungsten

PubMed Central

Suzudo, Tomoaki; Hasegawa, Akira

2016-01-01

Modeling the evolution of radiation-induced defects is important for finding radiation-resistant materials, which would be greatly appreciated in nuclear applications. We apply the density functional theory combined with comprehensive analyses of massive experimental database to indicate a mechanism to mitigate the effect of radiation on W crystals by adding particular solute elements that change the migration property of interstitials. The resultant mechanism is applicable to any body-centered-cubic (BCC) metals whose self-interstitial atoms become a stable crowdion and is expected to provide a general guideline for computational design of radiation-resistant alloys in the field of nuclear applications. PMID:27824134

EFFECT OF STRAIN FIELD ON THRESHOLD DISPLACEMENT ENERGY OF TUNGSTEN STUDIED BY MOLECULAR DYNAMICS SIMULATION

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wang, D.; Gao, Ning; Setyawan, Wahyu

The influence of hydrostatic strain on point defect formation energy and threshold displacement energy (Ed) in body-centered cubic (BCC) tungsten was studied with molecular dynamics simulations. Two different tungsten potentials (Fikar and Juslin) were used. The minimum Ed direction calculated with the Fikar-potential was <100>, but with the Juslin-potential it was <111>. The most stable self-interstitial (SIA) configuration was a <111>-crowdion for both potentials. The stable SIA configuration did not change with applied strain. Varying the strain from compression to tension increased the vacancy formation energy but decreased the SIA formation energy. The SIA formation energy changed more significantly thanmore » for a vacancy such that Ed decreased with applied strain from compression to tension.« less

Low Temperature Diffusion Transformations in Fe-Ni-Ti Alloys During Deformation and Irradiation

NASA Astrophysics Data System (ADS)

Sagaradze, Victor; Shabashov, Valery; Kataeva, Natalya; Kozlov, Kirill; Arbuzov, Vadim; Danilov, Sergey; Ustyugov, Yury

2018-03-01

The deformation-induced dissolution of Ni3Ti intermetallics in the matrix of austenitic alloys of Fe-36Ni-3Ti type was revealed in the course of their cascade-forming neutron irradiation and cold deformation at low temperatures via employment of Mössbauer method. The anomalous deformation-related dissolution of the intermetallics has been explained by the migration of deformation-induced interstitial atoms from the particles into a matrix in the stress field of moving dislocations. When rising the deformation temperature, this process is substituted for by the intermetallics precipitation accelerated by point defects. A calculation of diffusion processes has shown the possibility of the realization of the low-temperature diffusion of interstitial atoms in configurations of the crowdions and dumbbell pairs at 77-173 K. The existence of interstitial atoms in the Fe-36Ni alloy irradiated by electrons or deformed at 77 K was substantiated in the experiments of the electrical resistivity measurements.

QUANTUM CALCULATIONS OF ENERGETICS OF RHENIUM CLUSTERS IN TUNGSTEN

DOE Office of Scientific and Technical Information (OSTI.GOV)

Setyawan, Wahyu; Nandipati, Giridhar; Roche, Kenneth J.

2015-09-22

Density functional theory was employed to explore the energetic properties of clusters up to size 2 of Re in W. While WW<111> is the most stable intrinsic dumbbell, ReW<110> is more stable than ReW<111>. However, when they are trapped by a substitutional Re (Re_s), ReW<111> becomes more stable than ReW<110>. In this case, the most stable configuration forms a ReWRe crowdion with the W atom in between the Re atoms. Simulations of a ReW[111] (dumbbell’s vector is from Re to W) approaching a Re_s along [111] indicate that the binding energy decreases from 0.83 eV at the first nearest neighbormore » (NN1) to 0.10 eV at NN3 and ~0 at NN4. In addition, while ReW<111> and ReW<110> are stable near a Re_s at NN1, the ReW<100> instantaneously rotates toward ReW<111>.« less

Effect of Strain Field on Threshold Displacement Energy of Tungsten Studied by Molecular Dynamics Simulation

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wang, Dong; Gao, Ning; Setyawan, W.

The influence of strain field on defect formation energy and threshold displacement energy (Ed) in body-centered cubic (BCC) tungsten (W) has been studied with molecular dynamics simulations. Two different W potentials (Fikar and Juslin) were compared and the results indicate that the connection distance and selected function linking the short-range and long-range portions of the potentials affects the threshold displacement energy and its direction-specific values. The minimum Ed direction calculated with the Fikar-potential is <100> and with the Juslin-potential is <111>. Nevertheless, the most stable self-interstitial configuration is found to be a <111>-crowdion for both potentials. This stable configuration doesmore » not change with applied strain. Varying the strain from compression to tension increases the vacancy formation energy but decreases the self-interstitial formation energy. The formation energy of a self-interstitial changes more significantly than a vacancy such that Ed decreases with applied hydrostatic strain from compression to tension.« less

Defect structures induced by high-energy displacement cascades in γ uranium

DOE Office of Scientific and Technical Information (OSTI.GOV)

Miao, Yinbin; Beeler, Benjamin; Deo, Chaitanya

Displacement cascade simulations were conducted for the c uranium system based on molecular dynamics. A recently developed modified embedded atom method (MEAM) potential was employed to replicate the atomic interactions while an embedded atom method (EAM) potential was adopted to help characterize the defect structures induced by the displacement cascades. The atomic displacement process was studied by providing primary knock-on atoms (PKAs) with kinetic energies from 1 keV to 50 keV. The influence of the PKA incident direction was examined. The defect structures were analyzed after the systems were fully relaxed. The states of the self-interstitial atoms (SIAs) were categorizedmore » into various types of dumbbells, the crowdion, and the octahedral interstitial. The voids were determined to have a polyhedral shape with {110} facets. The size distribution of the voids was also obtained. The results of this study not only expand the knowledge of the microstructural evolution in irradiated c uranium, but also provide valuable references for the radiation-induced defects in uranium alloy fuels.« less

Interstitial loop transformations in FeCr

DOE PAGES

Béland, Laurent Karim; Osetsky, Yuri N.; Stoller, Roger E.; ...

2015-03-27

Here, we improve the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) algorithm by integrating the Activation Relaxation Technique nouveau (ARTn), a powerful open-ended saddle-point search method, into the algorithm. We use it to investigate the reaction of 37-interstitial 1/2[1 1 1] and 1/2[View the MathML source] loops in FeCr at 10 at.% Cr. They transform into 1/2[1 1 1], 1/2[View the MathML source], [1 0 0] and [0 1 0] 74-interstitial clusters with an overall barrier of 0.85 eV. We find that Cr decoration locally inhibits the rotation of crowdions, which dictates the final loop orientation. Moreover, the final loop orientationmore » depends on the details of the Cr decoration. Generally, a region of a given orientation is favored if Cr near its interface with a region of another orientation is able to inhibit reorientation at this interface more than the Cr present at the other interfaces. Also, we find that substitutional Cr atoms can diffuse from energetically unfavorable to energetically favorable sites within the interlocked 37-interstitial loops conformation with barriers of less than 0.35 eV.« less

Kinetics of self-interstitial migration in bcc and fcc transition metals

NASA Astrophysics Data System (ADS)

Bukkuru, S.; Bhardwaj, U.; Srinivasa Rao, K.; Rao, A. D. P.; Warrier, M.; Valsakumar, M. C.

2018-03-01

Radiation damage is a multi-scale phenomenon. A thorough understanding of diffusivities and the migration energies of defects is a pre-requisite to quantify the after-effects of irradiation. We investigate the thermally activated mobility of self-interstitial atom (SIA) in bcc transition metals Fe, Mo, Nb and fcc transition metals Ag, Cu, Ni, Pt using molecular dynamics (MD) simulations. The self-interstitial diffusion involves various mechanisms such as interstitialcy, dumbbell or crowdion mechanisms. Max-Space Clustering (MSC) method has been employed to identify the interstitial and its configuration over a wide range of temperature. The self-interstitial diffusion is Arrhenius like, however, there is a slight deviation at high temperatures. The migration energies, pre-exponential factors of diffusion and jump-correlation factors, obtained from these simulations can be used as inputs to Monte Carlo simulations of defect transport. The jump-correlation factor shows the degree of preference of rectilinear or rotational jumps. We obtain the average jump-correlation factor of 1.4 for bcc metals and 0.44 for fcc metals. It indicates that rectilinear jumps are preferred in bcc metals and rotational jumps are preferred in fcc metals.

Elastic fields, dipole tensors, and interaction between self-interstitial atom defects in bcc transition metals

NASA Astrophysics Data System (ADS)

Dudarev, S. L.; Ma, Pui-Wai

2018-03-01

Density functional theory (DFT) calculations show that self-interstitial atom (SIA) defects in nonmagnetic body-centered-cubic (bcc) metals adopt strongly anisotropic configurations, elongated in the <111 > direction [S. Han et al., Phys. Rev. B 66, 220101 (2002), 10.1103/PhysRevB.66.220101; D. Nguyen-Manh et al., Phys. Rev. B 73, 020101 (2006), 10.1103/PhysRevB.73.020101; P. M. Derlet et al., Phys. Rev. B 76, 054107 (2007), 10.1103/PhysRevB.76.054107; S. L. Dudarev, Annu. Rev. Mater. Res. 43, 35 (2013), 10.1146/annurev-matsci-071312-121626]. Elastic distortions, associated with such anisotropic atomic structures, appear similar to distortions around small prismatic dislocation loops, although the extent of this similarity has never been quantified. We derive analytical formulas for the dipole tensors of SIA defects, which show that, in addition to the prismatic dislocation looplike character, the elastic field of a SIA defect also has a significant isotropic dilatation component. Using empirical potentials and DFT calculations, we parametrize dipole tensors of <111 > defects for all the nonmagnetic bcc transition metals. This enables a quantitative evaluation of the energy of elastic interaction between the defects, which also shows that in a periodic three-dimensional simple cubic arrangement of crowdions, long-range elastic interactions between a defect and all its images favor a <111 > orientation of the defect.

Analytic modified embedded atom potentials for HCP metals

NASA Astrophysics Data System (ADS)

Hu, Wangyu; Zhang, Bangwei; Huang, Baiyun; Gao, Fei; Bacon, David J.

2001-02-01

Analytic modified embedded atom method (AMEAM) type many-body potentials have been constructed for ten hcp metals: Be, Co, Hf, Mg, Re, Ru, Sc, Ti, Y and Zr. The potentials are parametrized using analytic functions and fitted to the cohesive energy, unrelaxed vacancy formation energy, five independent second-order elastic constants and two equilibrium conditions. Hence, each of the constructed potentials represents a stable hexagonal close-packed lattice with a particular non-ideal c/a ratio. In order to treat the metals with negative Cauchy pressure, a modified term has been added to the total energy. For all the metals considered, the hcp lattice is shown to be energetically most stable when compared with the fcc and bcc structure and the hcp lattice with ideal c/a. The activation energy for vacancy diffusion in these metals has been calculated. They agree well with experimental data available and those calculated by other authors for both monovacancy and divacancy mechanisms and the most possible diffusion paths are predicted. Stacking fault and surface energy have also been calculated and their values are lower than typical experimental data. Finally, the self-interstitial atom (SIA) formation energy and volume have been evaluated for eight possible sites. This calculation suggests that the basal split or crowdion is the most stable configuration for metals with a rather large deviation from the ideal c/a value and the non-basal dumbbell (C or S) is the most stable configuration for metals with c/a near ideal. The relationship between SIA formation energy and melting temperature roughly obeys a linear relation for most metals except Ru and Re.

Annihilating vacancies via dynamic reflection and emission of interstitials in nano-crystal tungsten

NASA Astrophysics Data System (ADS)

Li, Xiangyan; Duan, Guohua; Xu, Yichun; Zhang, Yange; Liu, Wei; Liu, C. S.; Liang, Yunfeng; Chen, Jun-Ling; Luo, G.-N.

2017-11-01

Radiation damage not only seriously degrades the mechanical properties of tungsten (W) but also enhances hydrogen retention in the material. Introducing a large amount of defect sinks, e.g. grain boundaries (GBs) is an effective method for improving radiation-resistance of W. However, the mechanism by which the vacancies are dynamically annihilated at long timescale in nano-crystal W is still not clear. The dynamic picture for eliminating vacancies with single interstitials and small interstitial-clusters has been investigated by combining molecular dynamics, molecular statics and object Kinetic Monte Carlo methods. On one hand, the annihilation of bulk vacancies was enhanced due to the reflection of an interstitial-cluster of parallel ≤ft< 1 1 1 \\right> crowdions by the GB. The interstitial-cluster was observed to be reflected back into the grain interior when approaching a locally dense GB region. Near this region, the energy landscape for the interstitial was featured by a shoulder, different to the decreasing energy landscape of the interstitial near a locally loose region as indicative of the sink role of the GB. The bulk vacancy on the reflection path was annihilated. On the other hand, the dynamic interstitial emission efficiently anneals bulk vacancies. The single interstitial trapped at the GB firstly moved along the GB quickly and clustered to be the di-interstitial therein, reducing its mobility to a value comparable to that that for bulk vacancy diffusion. Then, the bulk vacancy was recombined via the coupled motion of the di-interstitial along the GB, the diffusion of the vacancy towards the GB and the accompanying interstitial emission. These results suggest that GBs play an efficient role in improving radiation-tolerance of nano-crystal W via reflecting highly-mobile interstitials and interstitial-clusters into the bulk and annihilating bulk vacancies, and via complex coupling of in-boundary interstitial diffusion, clustering of the interstitial and vacancy diffusion in the bulk.

Primary damage formation in bcc iron

NASA Astrophysics Data System (ADS)

Stoller, R. E.; Odette, G. R.; Wirth, B. D.

1997-11-01

Primary defect formation in bee iron has been extensively investigated using the methods of molecular dynamics (MD) and Monte Carlo (MC) simulation. This research has employed a modified version of the Finnis-Sinclair interatomic potential. MD was used in the simulation of displacement cascades with energies up to 40 keV and to examine the migration of the interstitial clusters that were observed to form in the cascade simulations. Interstitial cluster binding energies and the stable cluster configurations were determined by structural relaxation and energy minimization using a MC method with simulated annealing. Clusters containing up to 19 interstitials were examined. Taken together with the previous work, these new simulations provide a reasonably complete description of primary defect formation in iron. The results of the displacement cascade simulations have been used to characterize the energy and temperature dependence of primary defect formation in terms of two parameters: (1) the number of surviving point defects and (2) the fraction of the surviving defects that are contained in clusters. The number of surviving point defects is expressed as a fraction of the atomic displacements calculated using the secondary displacement model of Norgett-Robinson-Torrens (NRT). Although the results of the high energy simulations are generally consistent with those obtained at lower energies, two notable exceptions were observed. The first is that extensive subcascade formation at 40 keV leads to a higher defect survival fraction than would be predicted from extrapolation of the results obtained for energies up to 20 keV. The stable defect fraction obtained from the MD simulations is a smoothly decreasing function up to 20 keV. Subcascade formation leads to a slight increase in this ratio at 40 keV, where the value is about the same as at 10 keV. Secondly, the potential for a significant level of in-cascade vacancy clustering was observed. Previous cascade studies employing this potential have reported extensive interstitial clustering, but little evidence of vacancy clustering. Interstitial clusters were found to be strongly bound, with binding energies in excess of 1 eV. The larger clusters exhibited a complex, 3D structure and were composed of <111> crowdions. These clusters were observed to migrate by collective <111> translations with an activation energy on the order of 0.1 eV.