Search NASA⌕ Search

SEARCH · Search NASA

Results for “Frenkel pair”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Recombination radius of a Frenkel pair and capture radius of a self-interstitial atom by vacancy clusters in bcc Fe

The recombination radius of a Frenkel pair is a fundamental parameter for the object kinetic Monte Carlo (OKMC) and mean field rate theory (RT) methods that are used to investigate irradiation damage accumulation in neutron irradiated nuclear materials. The recombination radius in bcc Fe has been studied both experimentally and numerically, however there is no general consensus about its value. The detailed atomistic processes of recombination also remain uncertain. Values from 1:0a₀ to 3:3a₀ have been employed as a recombination radius in previous studies using OKMC and RT. The recombination process of a Frenkel pair is investigated at the atomic level using the self-evolved atomistic kinetic Monte Carlo (SEAKMC) method in this paper. SEAKMC calculations reveal that a self-interstitial atom recombines with a vacancy in a spontaneous reaction from several nearby sites following characteristic pathways. The recombination radius of a Frenkel pair is estimated to be 2.26a₀ by taking the average of the recombination distances from 80 simulation cases. This value agrees well with the experimental estimate. In addition, we apply these procedures to the capture radius of a self-interstitial atom by a vacancy cluster. The capture radius is found to gradually increase with the size of the vacancy cluster. The fitting curve for the capture radius is obtained as a function of the number of vacancies in the cluster.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Frenkel pairs cause elastic softening in zirconia: theory and experiments

Abstract Recent results from molecular dynamics simulations have shown that significant concentrations of Frenkel pairs can be introduced by the proliferation of phonons lying at the edge of the Brillouin zone and when above the Debye temperature. Following the work of Granato (2014 Eur. Phys. J. B 87 18) we extend those calculations to the influence of Frenkels on the elastic modulus. Significant softening is predicted which is confirmed by in situ measurements of the elastic modulus during flash. Frenkel pairs have been proposed to play a central role in the flash phenomena.

Physics↗
DOE OSTI2017

Low energy ion-solid interactions and chemistry effects in a series of pyrochlores

The effect of chemistry on low energy recoil events was investigated at 10 K for each type of atom in pyrochlores using molecular dynamics simulation. Contour plots of the threshold displacement energy (Ed) in Gd2Zr2O7 have been produced along more than 80 directions for each individual species. The Ed surface for each type of atom in Gd2Zr2O7 is highly anisotropic; Ed of Zr exhibits the largest degree of anisotropy, while that of O8b exhibits the smallest. The recommended values of Ed in Gd2Zr2O7 based on the observed minima are 56, 94 and 25 eV, respectively for Gd, Zr and O. The influence of cation radius on Ed in pyrochlores A2B2O7 (with A-site ranging from Lu3+ to La3+ and B-site ranging from Ti4+ to Ce4+) was also investigated along three directions [100], [110] and [111]. The Ed in pyrochlores strongly depended on the atom type, atom mass, knock-on direction, and lattice position. The defects produced after low energy displacement events included cation antisite defects, cation Frenkel pairs, anion Frenkel pairs, various vacancies and interstitials. Ce doping in pyrochlores may affect the radiation response, because it resulted in drastic changes in cation and anion displacement energies and formation of an unusual type of anti-site defect. This work demonstrates links between Ed and amorphization resistance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Ab initio molecular dynamics simulations of ion-solid interactions in zirconate pyrochlores

In this paper, an ab initio molecular dynamics method is employed to study low energy recoil events in zirconate pyrochlores (A 2 Zr 2 O 7 , A = La, Nd and Sm). It shows that both cations and anions in Nd 2 Zr 2 O 7 and Sm 2 Zr 2 O 7 are generally more likely to be displaced than those in La 2 Zr 2 O 7 . The damage end states mainly consist of Frenkel pair defects, and the Frenkel pair formation energies in Nd 2 Zr 2 O 7 and Sm 2 Zr 2 O 7 are lower than those in La 2 Zr 2 O 7 . These results suggest that the order–disorder structural transition more easily occurs in Nd 2 Zr 2 O 7 and Sm 2 Zr 2 O 7 resulting in a defect-fluorite structure, which agrees well with experimental observations. Our calculations indicate that oxygen migration from 48f and 8b to 8a sites is dominant under low energy irradiation. A number of new defects, including four types of cation Frenkel pairs and six types of anion Frenkel pairs, are revealed by ab initio molecular dynamics simulations. The present findings may help to advance the fundamental understanding of the irradiation response behavior of zirconate pyrochlores.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Light–Induced Frenkel Defect Pair Formation Can Lead to Phase–Segregation of Otherwise Miscible Halide Perovskite Alloys

Alloys of ABX 3 halide perovskites (HP) exhibit unique phase behavior compared to traditional III-V and II-VI semiconductor alloys used in solar cells. While the latter typically have good mutual miscibility when their mixed components are size matched, and phase-segregate when size mismatched, HP alloys show good miscibility in the dark but can phase-segregate under light. Quantum mechanical calculations described herein reveal light-induced defect formation and migration hold the key. Specifically, the interaction between a halogen vacancy V X with halogen interstitial X i forming together a Frenkel-pair defect emerges as the enabler for phase-segregation in HP alloys. At a threshold bromine composition in the Br-I alloys, the photogenerated holes in the valence band localize, creating thereby a doubly-charged iodine Frenkel-pair (V I + I i ) 2+ . Faster migration of iodine over bromine interstitial into the vacant iodine V I site leads to the formation of iodine-rich and iodine-depleted regions, establishing phase-segregation. Removal of the mobile defects–the agent of segregation–by dark thermal annealing, supplies the opposing force, leading to reversal of phase-segregation. Furthermore, this atomistic understanding can enable some control of the phase-segregation by selecting substituting elements on the B site–such as replacing some Pb by Sn–that are unable to form stable Frenkel defects.

14 SOLAR ENERGY↗
DOE OSTI2017

Molecular Dynamics Simulation of Fission Fragment Damage in Nuclear Fuel and Surrogate Material

We have performed classical molecular dynamics simulations of swift heavy ion damage, typical of fission fragments, in nuclear fuel (UO2) for energy deposition per unit length of 3.9 keV/nm. We did not observe amorphization. The damage mainly consisted of isolated point defects. Only about 1% of the displacements occur on the uranium sublattice. Oxygen Frenkel pairs are an order of magnitude more numerous than uranium Frenkel pairs in the primary damage state. In contrast, previous results show that the ratio of Frenkel pairs on the two sublattices is close to the stoichiometric ratio in ceria. These differences in the primary damage state may lead to differences in radiation response of UO2 and CeO2.

Molecular dynamics, radiation damage in materials,↗
DOE OSTI2021

The mobility of small, over-pressurized helium bubbles in tungsten at 2000 K

Fusion reactor environments inevitably lead to the formation of high-pressure helium bubbles whose nucleation, growth, and diffusion strongly impact the performance of plasma-facing components. Here this research describes a diffusion mechanism of over-pressurized bubbles via a sequence of Frenkel pair nucleation, self-interstitial migration, and Frenkel pair annihilation. Molecular dynamics was used to simulate the diffusion of small bubbles in tungsten at 2000 K with helium-per-vacancy ratios in the range of 4.5 to 7. The diffusion coefficients are calculated and their dependence on helium content, number of vacancies, and number of attached self-interstitials is characterized. It is found that bubbles are most mobile when the nucleation/annihilation rates of Frenkel pairs are nearly equal and when the bubbles nucleate and annihilate a single self-interstitial. All bubbles experience a peak diffusivity, which can be as high as 10 -11 m 2 /s decreasing with bubble size. The calculated diffusion coefficients provide valuable insight into the mobility of small, high-pressure bubbles, and can be used as input parameters in mesoscale models to improve predictions of plasma-surface interactions. (LA-UR-21-21881)

36 MATERIALS SCIENCE↗
DOE OSTI2014

A Code of Radiation Damage by Gamma - Paper 62

In recent years, the study of radiation damage of materials has seen a thrust because of the expectation of the high rate of radiation damage in fusion reactors, accelerator driven sub critical systems and damage of a control system of space vehicles by the high energy neutrons and gamma showers. Embrittlement of reactor vessels, corrosion of steel and swelling of fuel rods are the well known empirical examples of severe radiation damage. It is also realized that damage by high energy gammas is in no way insignificant, and it is equivalent to the damage by the neutrons. Secondly, value addition of certain materials like enhancement of electrical resistance, mechanical strength, color enhancement of precious stones on irradiation by neutrons and gammas and polymerization in the presence of gamma radiation add to the necessity of a Monte Carlo (MC) code of radiation damage in a wide range of energy. The JAIPU code, for the study of radiation damage by neutrons up to 10 MeV energy, incorporates features like development of secondary atomic cascade and PKA energy corrected for the ionization lose and production of Frenkel pair. MC estimation of both the primary and secondary displacements has been compared with the Frenkel pairs calculated using the NRT model. In the present work, JAIPU code is extended for irradiation of materials by high energy gammas. Here cascade generated by the incident gamma produces electrons by way of photo-electric, Compton and pair production until it is absorbed. The recoiled Frenkel pairs produced in electron + atom collision have been considered for the purpose of atomic cascade. Standard cross section libraries for gamma and electron interaction are used while for atom + atom collision we have made use of the IOTA code. In addition, in the code, the lattice structure of the irradiated material is simulated to know the positions of the photon interactions with atoms and the energy of electron it adds an extra feature of the position and size of the voids and the clusters produced on irradiation. The code has the possibility of calculation of 'absorbed radiation dose' and 'KERMA' and the interatomic distances of the displaced atoms. It may be pointed out that the code uses only the MC simulation of the number of displacements unlike several other approaches where the NRT model with PKA energy is used for gross estimation of the number of the displaced atoms. Thus, the code keeps an account of the positions of displaced atoms for the in-depth modification of the irradiated material. 1 MeV gamma seems to damage Fe up to several millimeters of its depth and this is mostly by the recoiled Compton electrons than photo electrons. The code is in the developmental stage and it has the capability of applications in the developing field of Radiology besides the shielding and other studies of radiation damage. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Effect of the energy of bombarding electrons on the conductivity of n-4H-SiC (CVD) epitaxial layers

The electrical characteristics of epitaxial layers of n-4H-SiC (CVD) irradiated with 0.9 and 3.5MeV electrons are studied. It is shown that the donor removal rate becomes nearly four times higher as the energy of impinging electrons increases by a factor of 4, although the formation cross section of primary radiation defects (Frenkel pairs in the carbon sublattice) responsible for conductivity compensation of the material is almost energy independent in this range. It is assumed that the reason for the observed differences is the influence exerted by primary knocked-out atoms. First, cascade processes start to manifest themselves with increasing energy of primary knocked-out atoms. Second, the average distance between genetically related Frenkel pairs grows, and, as a consequence, the fraction of defects that do not recombine under irradiation becomes larger. The recombination radius of Frenkel pairs in the carbon sublattice is estimated and the possible charge state of the recombining components is assessed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecular dynamics simulations of displacement cascades in LiAlO2 and LiAl5O8 ceramics

Abstract Molecular dynamics was employed to investigate the radiation damage due to collision cascades in LiAlO 2 and LiAl 5 O 8 , the latter being a secondary phase formed in the former during irradiation. Atomic displacement cascades were simulated by initiating primary knock-on atoms (PKA) with energy values = 5, 10 and 15 keV and the damage was quantified by the number of Frenkel pairs formed for each species: Li, Al and O. The primary challenges of modeling an ionic system with and without a core–shell model for oxygen atoms were addressed and new findings on the radiation resistance of these ceramics are presented. The working of a variable timestep function and the kinetics in the background of the simulations have been elaborated to highlight the novelty of the simulation approach. More importantly, the key results indicated that LiAlO 2 experiences much more radiation damage than LiAl 5 O 8 , where the number of Li Frenkel pairs in LiAlO 2 was 3–5 times higher than in LiAl 5 O 8 while the number of Frenkel pairs for Al and O in LiAlO 2 are ~ 2 times higher than in LiAl 5 O 8 . The primary reason is high displacement threshold energies (E d ) in LiAl 5 O 8 for Li cations. The greater E d for Li imparts higher resistance to damage during the collision cascade and thus inhibits amorphization in LiAl 5 O 8 . The presented results suggest that LiAl 5 O 8 is likely to maintain structural integrity better than LiAlO 2 in the irradiation conditions studied in this work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Negative thermal expansion coefficient and amorphization in defective 4H-SiC

This paper presents thermal expansion coefficient (TEC) and amorphization in 4H-SiC containing point defects at different concentrations. We considered vacancy defects, interstitial defects, and Frenkel pair defects and investigated the thermomechanical response of the lattice over a wide range of temperatures using classical molecular dynamics simulations. The results show that 4H-SiC with vacancy defects exhibits a negative TEC above a critical defect density of around 9% (irrespective of the temperature). With interstitial defects, it exhibits a positive TEC (regardless of the defect density), and with Frenkel pair defects it shows a transition from positive TEC to negative TEC for a defect density greater than 8%. The coupling between temperature-induced expansion and defect-introduced stress in the lattice forms the mechanistic basis for the observed variation in TEC. Furthermore, the specific heat decreases rapidly with an increase in defect density at room temperature, with the highest sensitivity of the lattice observed for the Frenkel pair defects followed by interstitial defects and then by vacancy defects. Finally, these findings highlight the critical implications of defects on thermal expansion behavior of 4H-SiC with applications in radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Synchronous parallel spatially resolved stochastic cluster dynamics

In this work, a spatially resolved stochastic cluster dynamics (SRSCD) model for radiation damage accumulation in metals is implemented using a synchronous parallel kinetic Monte Carlo algorithm. The parallel algorithm is shown to significantly increase the size of representative volumes achievable in SRSCD simulations of radiation damage accumulation. Additionally, weak scaling performance of the method is tested in two cases: (1) an idealized case of Frenkel pair diffusion and annihilation, and (2) a characteristic example problem including defect cluster formation and growth in α-Fe. For the latter case, weak scaling is tested using both Frenkel pair and displacement cascade damage. To improve scaling of simulations with cascade damage, an explicit cascade implantation scheme is developed for cases in which fast-moving defects are created in displacement cascades. For the first time, simulation of radiation damage accumulation in nanopolycrystals can be achieved with a three dimensional rendition of the microstructure, allowing demonstration of the effect of grain size on defect accumulation in Frenkel pair-irradiated α-Fe.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Coupling between interfacial charge and mechanical deformation at high temperatures in ceramics. Final report

Basic research, under this project, has led to new fundamental discoveries in the effect of electric field (and temperature) on diffusion kinetics, electroluminescence and electronic conductivity in ceramic materials. For example, yttria stabilized zirconia sinters in just seconds well below 1000°C, while, at the same time, exhibiting a highly non-linear rise in conductivity, and intense electroluminescence. The phenomenon has been shown to occur in nearly all ceramics whether they are ionic or electronic conductors, or insulators. The rise in conductivity leads to Joule heating but in-situ measurements with a platinum standard at APS and Brookhaven (NSLS-II) have shown that the temperatures are much below what would have been needed to sinter in a few seconds. We have explored the hypothesis that huge concentrations of Frenkel pairs of vacancies and interstitials are produced. While these point defects accelerate sintering (by the interstitials being transported to the pores and the vacancies to the grain boundaries), they may also produce electron-hole pairs which can impart high electrical conductivity or recombine to emit photons. Many measurements, including in-situ measurements of lattice expansion are consistent with this hypothesis. However, the energy barrier for generating Frenkels in zirconia for example is much larger than can be expected at temperature and electrical fields in these experiments. Most recent work suggests that short wavelength phonons play a critical role in the genesis of the “flash” phenomenon. For example, the Debye temperature has been shown to be the lower bound for initiating flash. Molecular dynamics simulations have shown that large concentrations of Frenkel pairs can form under the proliferation of “optical” phonons at the edge of the Brillouin zone, that is with a wavelength at or below the lattice parameter. This project has led to the development of new experimental techniques as well as a software hardware interface for process and control and data analysis in real time, with a time scale of less than one second. This approach spells a new paradigm in materials science research. It potential impact on high rate manufacturing is self evident.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Displacement damage caused by gamma-rays and neutrons on Au and Se.

This report documents theoretical calculations of displacement damage produced by gamma rays and neutrons on various materials. The average energy of the gamma rays was 1.24 MeV and 1.0 MeV for the neutrons. The fluence of the gamma rays was 1.2e14 γ/cm 2 , for the neutrons it was 1.0e12 n/cm2. The initial materials of interest were Au and Se. The total doses of the gamma ray exposures were in the 100 kRad range for both elements. An equivalent electron fluence was approximated to be the same as the gamma ray fluence over one gamma ray attenuation length in both materials and at the same 1.24 MeV energy. The maximum recoil energy of the Au and Se for these electrons was calculated relativisticaly to be 29 and 72 eV respectively. The relativisitic McKinley and Feshbach theory for the atomic recoil cross sections produced by the electrons were in the 10s of mbarn range and an upper limit for the concentration of Frenkel pairs for the gamma ray exposures for both elements was in the ppb range. The Robinson Energy Partioning Theory for non-ionizing energy loss (NIEL) of ions in solids was used to calculate the concentration of Frenkel pairs produced by the 1 MeV neutrons, and this concentration was also in the ppb range for both Au and Se. Low damage levels like this can have effects on minority carrier recombination in semiconductors, but are not expected to have any effect on metals like Au, or metalloids such as Se.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Formation of Radiation Defects by Proton Braking in Lightly Doped n- and p-SiC Layers

Mathematical simulation of the cascade of displacements in SiC is used to consider the specific features of Frenkel-pair generation upon the scattering of 8- and 15-MeV protons. The distribution histograms of energies acquired not only by primary knocked-out atoms, but also by recoil atoms generated in displacement cascades, are calculated. An analysis of the histograms considers two energy ranges. In the first range of “low” energies, the spontaneous recombination of genetically related Frenkel pairs is dominant. Recoil atoms in the second range have a higher energy, which enables these atoms to leave the spontaneousrecombination zone and dissociate into isolated components. The compensation of lightly doped n- and p-4H-SiC samples grown by gas-phase epitaxy is experimentally studied under irradiation with 8- and 15-MeV protons. The carrier removal rates are measured. The calculated and experimental data are compared and estimates are obtained for the size of the spontaneous-recombination zone.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Chemical instability leads to unusual chemical-potential-independent defect formation and diffusion in perovskite solar cell material CH 3 NH 3 PbI 3

Methylammonium (MA) lead triiodide (MAPbI 3 ) has recently emerged as a promising solar cell material. But, MAPbI3 is known to have chemical instability, i.e., MAPbI3 is prone to decomposition into MAI and PbI 2 even at moderate temperatures (e.g. 330 K). Here, we show that the chemical instability, as reflected by the calculated negligible enthalpy of formation of MAPbI 3 (with respect to MAI and PbI 2 ), has an unusual and important consequence for defect properties, i.e., defect formation energies in low-carrier-density MAPbI 3 are nearly independent of the chemical potentials of constituent elements and thus can be uniquely determined. This allows straightforward calculations of defect concentrations and the activation energy of ionic conductivity (the sum of the formation energy and the diffusion barrier of the charged mobile defect) in MAPbI 3 . Furthermore, the calculated activation energy for ionic conductivity due to V$+\atop{1}$ diffusion is in excellent agreement with the experimental values, which demonstrates unambiguously that V$+\atop{1}$ is the dominant diffusing defect and is responsible for the observed ion migration and device polarization in MAPbI3 solar cells. The calculated low formation energy of a Frenkel pair (V$+\atop{1}$ -I$-\atop{i}$ and low diffusion barriers of V$+\atop{1}$ and Image I$-\atop{i}$ suggest that the iodine ion migration and the resulting device polarization may occur even in single-crystal devices and grain-boundary-passivated polycrystalline thin film devices (which were previously suggested to be free from ion-migration-induced device polarization), leading to device degradation. Moreover, the device polarization due to the Frenkel pair (which has a relatively low concentration) may take a long time to develop and thus may avoid the appearance of the current–voltage hysteresis at typical scan rates.

14 SOLAR ENERGY↗
DOE OSTI2015

The interstitialcy diffusion in FCC copper: A molecular dynamics study

Damage of materials due to neutron irradiation occurs via energetic cascades caused by energetic primary knock-on atoms (PKA) created by the energetic neutron as it passes through the material. These cascades result in creation of Frenkel Pairs (interstitials and vacancies). The interstitials and vacancies diffuse and recombine to (I) nullify the damage when an interstitial recombines with a vacancy, (II) form interstitial clusters when two or more interstitials recombine, and (III) form vacancy clusters when several vacancies come together. The latter two processes result in change of material properties. Interstitial diffusion has reported time-scales of microseconds and vacancy diffusion has diffusion time-scales of the order of seconds. We have carried out molecular dynamics (MD) simulations of interstitial diffusion in crystal Cu to study the mechanism of diffusion. It is found that interstitialcy diffusion – wherein an interstitial displaces a lattice atom thereby making the lattice atom an interstitial – has time-scales of a few tens of pico-seconds. Therefore we propose that the “interstitialcy diffusion” mechanism could play a major part in the diffusive-recombinations of the Frenkel Pairs created during the cascade.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗