Predicting the temperature dependence of self-diffusion behavior in Ni-Cr alloys via molecular dynamics
Not Available
Simonnin, Pauline (ORCID:0000000261977435)↗
SEARCH · Engineering Papers
Results for “Vacancy diffusion”
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.
Atomistic modeling of lanthanide diffusion in refractory body-centered cubic molybdenum
Lanthanide fission products can strongly interact with candidate cladding alloys, but their transport properties in refractory metals remain poorly understood. Here, in this work, we investigate the atomic-scale diffusion behavior of La, Ce, Pr, and Nd in body-centered cubic (bcc) molybdenum, a potential candidate for advanced nuclear cladding. Self-consistent mean-field transport modeling is performed to evaluate the fission product transport and vacancy mobility, informed by first-principles and nudged elastic band calculations of vacancy formation energies, migration barriers, and solute–vacancy binding characteristics. Compared with bcc Fe, lanthanide solutes in bcc Mo exhibit slower tracer diffusion due to higher vacancy formation and migration energies. Furthermore, the calculations reveal that the influence of fission products on migration barriers in bcc Mo are not as extensive in range compared to bcc Fe. Among the studied lanthanides, La exhibits the strongest vacancy binding while also being the fastest diffuser in Mo. These findings highlight how refractory bcc alloys can reduce fission product infiltration, offering valuable insight into the development of durable cladding systems for advanced reactors.
36 - MATERIALS SCIENCE↗
Insights into oxygen diffusion in rare earth disilicate environmental barrier coatings
Environmental barrier coatings (EBC) are crucial for the use of SiC-based ceramic matrix composites in high-temperature combustion environments, yet knowledge of oxygen diffusion in these coatings is limited. This study investigates oxygen diffusion dynamics in the β-RE 2 Si 2 O 7 system to minimize oxygen penetration in rare earth disilicates. We analyze defect formation energy under varying oxygen conditions, identifying key diffusion mechanisms. In oxygen-rich environments, the most favorable neutral interstitial oxygen diffuses along the [110] direction. In oxygen-poor conditions, neutral oxygen vacancies rotate around Y and Si atoms, exhibiting a diffusivity of 6.59×10 −22 m 2 /s at 1500 K for β-Yb 2 Si 2 O 7 . Under intermediate oxygen levels, charged interstitial oxygen diffuses via concerted interstitialcy along the [001] direction with a diffusivity of 6.21×10 −17 m 2 /s. Additionally, alloying rare earth Y with Er and Yb increases diffusion barriers, contributing to improved EBC performance in extreme environments. The insights gained provides valuable guidance for designing robust coatings tailored to withstand extreme operational environments.
36 MATERIALS SCIENCE↗
Active and passive defects in tetragonal tungsten bronze relaxor ferroelectrics
Tetragonal tungsten bronze (TTB) based oxides constitute a large family of dielectric materials which are known to exhibit complex distortions producing incommensurately modulated superstructures as well as significant local deviations from their average symmetry. The local deviations produce diffuse scattering in diffraction experiments. The structure as well as the charge dynamics of these materials are anticipated to be sensitive to defects, such as cation or oxygen vacancies. In this work, in an effort to understand how the structural and charge dynamical properties respond to these two types of vacancy defects, we have performed measurements of dielectric susceptibilities and single crystal diffraction experiments of two types of TTB materials with both ‘filled’ (Ba 2 NdFeNb 4 O 15 and Ba 2 PrFeNb 4 O 15 ) and ‘unfilled’ (Sr 0.5 Ba 0.5 Nb 2 O 6 ) cation sublattices. We also perform these measurements before and after oxygen annealing, which alters the oxygen vacancy concentrations. Surprisingly, we find that many of the diffuse scattering features that are present in the unfilled structure are also present in the filled structure, suggesting that the random fields and disorder that are characteristic of the unfilled structure are not responsible for many of the local structural features that are reflected in the diffuse scattering. Furthermore, oxygen annealing clearly affected both color and dielectric properties, consistent with a diminishment of the oxygen vacancy concentration, but had little effect on observed diffuse patterns.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Discrete stochastic model of point defect-dislocation interaction for simulating dislocation climb
Dislocation climb is an important high temperature process in metals plasticity, responsible for the phenomena such as creep, swelling, or hardening. Climb is defined by the ability of dislocations to leave their original glide plane by interacting with point defects. As such, dislocation climb is controlled by point defect diffusion/absorption/emission, all of which involve thermal activation. The existing thermodynamically consistent models for climb are generally formulated in a continuum framework, through the definition of effective defect fluxes and climb propensities in response to thermodynamic driving forces. However, the point-wise discrete nature of vacancies (and/or self-interstitials) confers a highly discrete nature to the climb dynamics, which is also strongly affected by elastic forces. The combination of discreteness, thermal activation, and elasticity makes this process too challenging for direct atomistic methods such as molecular dynamics. Here we develop a kinetic Monte Carlo model that captures vacancy generation and transport kinetics acting in conjuction with the evolving elastic fields provided by discrete dislocation dynamics simulations. The two models are coupled via the applied stresses and stress gradients generated by dislocation structures at vacancy locations. Our simulations reveal two surprising results. First that climb is dominated by vacancy emission even when the background vacancy concentration is much higher than the equilibrium one. And, second, that climb velocities might be much faster than otherwise believed when one uses the classical theories of climb. These effects are due to the locality of vacancy-dislocation processes, which are not captured in classical treatments that assume smooth vacancy fluxes and homogeneous concentrations. We apply the method to study elementary climb processes in body-centered cubic iron and furnish climb mobility functions to be used in parametric dislocation dynamics and/or crystal plasticity simulations. Here, we apply the technique to study non-conservative plastic bypass of spherical precipitates by edge dislocations and point out the differences between our discrete approach and existing continuum formulations.
36 MATERIALS SCIENCE↗
Closed-ampoule diffusion of sulfur into Cd-doped InP substrates - Dependence of S profiles on diffusion temperature and time
In order to optimize the fabrication of n(+)-p InP solar cells made by closed-ampoule diffusion of sulfur into p-InP:Cd substrates, we have investigated the influence of diffusion conditions on sulfur diffusion profiles. We show that S diffusion in InP is dominated by the P vacancy mechanism and is not characterized by a complementary error function as expected for an infinite source diffusion. The S diffusion mechanism in p-InP is qualitatively explained by examining the depth profiles of S, P, and In in the emitter layer and by taking into account the presence and composition of different compounds found to form in the In-P-S-O-Cd system as a result of diffusion.
Faur, Mircea↗
Investigation of Low-Energy Lattice Dynamics and Their Role in Superionic Na Diffusion and Ultralow Thermal Conductivity of Na 3 PSe 4 as a Solid-State Electrolyte
The atomic dynamics of Na 3 PSe 4 were investigated using a combination of neutron scattering experiments and ab initio and machine-learned molecular dynamics simulations to probe the interplay of fast ionic diffusion with atomic vibrations (phonons) of the host lattice. Our results reveal the existence of low-energy vibrational modes, simultaneously involving motions of Na + ions and framework polyanion subunits, and show that these modes become strongly overdamped in the superionic regime as they couple with the Na + hopping process. In particular, the Na+ migration energy landscape is strongly impacted by low-energy phonons derived from a soft acoustic branch of the host lattice, which modulates the diameter of the Na + diffusion channel at the bottleneck. We find that an additional factor for the enhanced Na + conductivity in Na 3 PSe 4 is the presence of Na-vacancies, which also affect the low-frequency dynamics and thermal vibration amplitudes, pointing to an interplay between Na + vacancies and host dynamics, jointly enhancing ionic diffusivity. Finally, we investigate the origin of ultralow thermal conductivities in Na 3 PSe 4 and Na 3 PS 4 using Green–Kubo simulations and find that low-energy acoustic phonon modes of the overall crystal framework provide a dominant contribution to the thermal conductivity.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
The Sluggish Diffusion of Cations in CeO 2 Probed through Molecular Dynamics and Metadynamics Simulations
Cation diffusion in fluorite‐structured CeO 2 , though far slower than anion diffusion, is an important, high‐temperature process because it governs diverse fabrication and degradation phenomena. Herein, cation diffusion is studied by means of classical molecular dynamics and metadynamics simulations. Three different mechanisms are examined: migration involving an isolated cerium vacancy, migration involving a cerium vacancy in a defect associate with an oxygen vacancy, and migration involving a cation divacancy. For each mechanism, defect diffusion coefficients are calculated as a function of temperature, from which the respective activation enthalpy of defect migration is obtained. Through comparisons with experimental cation diffusion data (specifically, of the absolute magnitude of the cation diffusivity as well as its activation enthalpy), it is concluded that cation diffusion takes place predominantly neither by isolated vacancies nor by cation vacancy–oxygen vacancy associates but by cation divacancies.
36 MATERIALS SCIENCE↗
Diffusion Behavior of Oversized Fission Products in bcc Fe Cladding: A First-Principles Study
Fuel-Cladding Chemical Interaction (FCCI) poses significant challenges in nuclear reactors, where fission products from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of oversized fission products, Pr, Nd, Ce, and La, within bcc Fe cladding using density functional theory (DFT), nudged elastic band (NEB) method, and self-consistent mean field (SCMF) theory. Our results reveal significant long-range vacancy binding energies, particularly up to the 6th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The NEB calculations indicate significant high barriers for the dissociation of 1nn vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products was derived in Arrhenius form. The significant trapping effect of vacancies by a very dilute amount of fission products reduces vacancy mobility, leading to an oversaturation of point defects, void nucleation, and swelling. These are critical issues for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. This study provides essential insights for developing advanced cladding materials and design strategies to mitigate FCCI, ultimately enhancing nuclear reactor safety and performance.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Mechanism of Metal Intercalation under Graphene through Small Vacancy Defects
Metal intercalation under graphene has attracted extensive experimental and theoretical research because of its capability to manipulate the electronic structure and properties of graphene. However, the pathways and mechanisms of intercalation are still not well understood. Here, we systematically investigate the intercalation process of metal atoms through graphene vacancies using first-principles calculations. We show that the energy barrier for metal atom penetration through the vacancies in graphene is small as long as the size of the vacancy is larger than a mono-vacancy. However, metal atoms are strongly bonded to the vacancy so that the detachment energy of a metal atom from the vacancy is extremely high. This inhibits the diffusion of the metal atom into the gallery beneath the surface to complete the intercalation process. On the other hand, our calculation results show that the detachment energy of a metal atom from a metal dimer at small vacancy defects is significantly reduced, making intercalation much easier. Therefore, the key step limiting the intercalation process is the detachment of the metal atoms from vacancy defects. This finding from our study provides useful insight into the defect-assisted intercalation mechanism.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
High-throughput approach for investigating interdiffusion in medium- and high-entropy alloys
Interdiffusion experiments are usually time-consuming and tedious since diffusion couples must be annealed at several temperatures for a long time. The efforts required to study interdiffusion in multicomponent alloys increase dramatically as multiple diffusion couples are required to cover broad composition ranges and determine the diffusivities of individual elements in different chemical environments. To circumvent this challenge, we present a high-throughput approach applicable to single-phase and compositionally complex alloys, which are assumed to approximate ideal solid solutions. Here, a simple diffusion-multiple experiment combined with a physically based kinetic model is proposed to efficiently determine the diffusion coefficients of the constituent elements in quaternary CrFeCoNi alloys. Compared with tracer diffusivities reported in the literature, the results, thus, obtained do not differ by more than a factor of 2 and were obtained from a single interdiffusion experiment. In contrast, the diffusivities simulated with commercial mobility and thermodynamic databases are strongly overestimated by a factor ranging from 1 to 16. Therefore, our approach enables high-throughput determination of diffusivities and can help in the design of alloys for high-temperature applications where diffusion plays a key role.
36 MATERIALS SCIENCE↗
Kinetic Monte Carlo simulations of aging in δ -Pu
We have developed a first-passage kinetic Monte Carlo approach for materials aging to investigate the sensitivity of void swelling to model parameters, including helium bubble density and size distribution. In addition to explicitly accounting for the spatial distribution of individual point defects, bubbles, and voids, our approach can simulate total doses equivalent to 100 years of natural aging on statistically representative volumes of materials. This technique enables us to study the effects on swelling and radiation damage evolution due to temperature and dose rate (as altered in artificially aged experiments), differences in effective interaction radii between vacancies and interstitials, and varying defect diffusion activation energies, while providing more detailed information than previous rate-equation based approaches. In conclusion, our results indicate that spatial effects that are not modeled in mean-field rate theories could play a significant role in void swelling initiation and growth for certain regimes of model parameters.
Actinides↗
Oxygen Vacancy Evolution at Li x V 2 O 5 /LiPON Solid State Electrochemical Interfaces Using Depth Resolved Cathodoluminescence Spectroscopy
The formation of oxygen vacancies at buried LiPON/ Li x V 2 O 5 interfaces has been observed on a near-nanometer scale and nondestructively using depth-resolved cathodoluminescence spectroscopy (DRCLS) and interfacial markers. Before electrochemical cycling, as-deposited LiPON/Li x V 2 O 5 exhibits a 1.6 eV defect optical emission, which density functional theory calculations identify as originating from oxygen vacancies. This defect appears first within a few nanometers of the buried LiPON/Li x V 2 O 5 interface without cycling, indicating that spontaneous O diffusion from the Li x V 2 O 5 lattice into LiPON may have caused these interface-localized oxygen vacancy defects. DRCLS measured the intensity and spatial distribution of this oxygen vacancy signal as a function of electrochemical cycling in a LiPON/Li x V 2 O 5 half-cell, showing oxygen vacancy signal increasing and moving deeper into the electrode with increased cycle number. Significant electrochemical irreversibility was also observed, with poor Coulombic efficiency and a 15% drop in capacity over 50 cycles. Theoretical simulations predict that the presence of oxygen vacancies increases the energy barrier for lithium diffusion significantly, indicating that this aggregation of oxygen vacancies could be another battery degradation mechanism accompanying lithiation induced phase changes.
36 MATERIALS SCIENCE↗
Temperature threshold for preferential bubble formation on grain boundaries in tungsten under in-situ helium irradiation
Understanding a material's radiation tolerance requires examining its performance under different irradiation conditions. In this work, we investigate the radiation tolerance in terms of helium bubble damage in tungsten irradiated in-situ with 16 keV helium at 1073 K and 1223 K. Damage evolution represented by helium bubble density, size and total change in volume in the grain matrices and the grain boundaries are quantified as a function of fluence. Preferential large bubble formation and a higher change in volume on the grain boundaries occurred at 1223 K, suggesting faster migration of certain helium-vacancy complexes as confirmed by a diffusion-reaction model.
36 MATERIALS SCIENCE↗
Cluster dynamics simulation of xenon diffusion during irradiation in UO 2
Diffusion of fission gas in UO 2 nuclear fuel impacts several important performance metrics, such as fission gas release, swelling, and thermal conductivity. Current empirical models of fission gas release have significant uncertainty, some of which derives from the bulk diffusion rate and its dependence on, for example, fuel chemistry and irradiation. In this work, we have applied the previously-developed Free Energy Cluster Dynamics (FECD) methodology in the code Centipede to calculate xenon cluster concentrations in UO 2 under intrinsic (high temperature) and irradiation-enhanced (intermediate temperature) conditions in order to develop a model of the xenon diffusion coefficient based on the atomic scale mechanisms responsible for transport. While the diffusion mechanism for xenon in UO 2 is adequately described by the Xe + U 2 O vacancy cluster for intrinsic conditions, a similar process is not capable of capturing measured in-pile fission gas diffusivity at intermediate temperatures. Therefore, a different diffusion mechanism must dominate under this regime. Using calculated atomistic data, we have shown that irradiation-enhanced diffusion at intermediate temperatures occurs via the larger Xe + U 4 O y vacancy clusters, which have lower migration barriers and increase in concentration by several orders of magnitude compared to intrinsic conditions. This mechanism is enabled by the increased uranium vacancy concentration under irradiation due to Frenkel pair production. In addition, the fast migration of uranium interstitials with two attached oxygen interstitials lowers the total uranium interstitial concentration through reactions with sinks. This allows the extended defects, such as Xe + U 4 O y vacancy clusters, to maintain high concentrations by limiting annihilation with attached vacancies. Furthermore, predictions using the Xe + U 4 O y diffusion mechanism are in good agreement with experiment, albeit with some differences in the Arrhenius slope, which we believe may be related to either experimental or model parameter uncertainty. Lastly, an analytical expression suitable for application in fuel performance simulations was derived to capture the predictions of the Centipede simulations.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Diffusion behavior of lanthanide fission products in bcc Fe cladding: A first-principles study
Fuel-cladding chemical interaction poses significant challenges in nuclear reactors, where fission products generated from nuclear fuel interact with Fe-based cladding materials, potentially compromising their structural integrity. This study investigates the diffusion behavior of lanthanide fission products, Lanthanum (La), Cerium (Ce), Praseodymium (Pr), and Neodymium (Nd), within body-centered cubic (bcc) Fe cladding using the density functional theory, nudged elastic band method, and self-consistent mean field theory. Our results reveal significant vacancy binding energies, particularly with the 1st and 2nd nearest neighbors, which diminish beyond the 5th nearest neighbor, with La exhibiting the strongest binding affinity, followed by Nd, Ce, and Pr. The nudged elastic band calculations indicate significant high barriers for the dissociation of 1st nearest neighbor vacancy-solute pairs for all fission products. The tracer diffusion coefficients of these fission products were derived in an Arrhenius form, with a magnetic correction that accounts for the high-temperature paramagnetic state. The significant trapping effect of vacancies caused by a very dilute concentration of fission products reduces vacancy mobility, potentially leading to modifications in point defect supersaturation, void nucleation, and swelling under irradiation. These represent critical challenges for irradiated cladding materials. The tracer diffusion coefficients indicate that Nd diffuses the fastest, followed by La, Ce, and Pr. Furthermore, this study provides essential insights for understanding fission product transport in cladding materials and informs future design strategies to mitigate fuel-cladding chemical interaction, ultimately enhancing nuclear reactor safety and performance.
Diffusion↗
Re-entrance to a ferromagnetic insulator with oxygen-vacancy ordering in the La 0.7 Sr 0.3 MnO 3 /SrTiO 3 superlattice
The charge-density wave is correlated with the properties in transition metal oxides, and regulating this state is important to explore the potential of oxides. Defects could modulate the local environment, in which the effect is subject to the concentration and arrangement, and an artificial arrangement of defect-ordering is highly expected to design the material properties. In this work, the ordering of oxygen vacancies with a preferred location on the (001) La/SrO plane by self-structural optimization was induced in La 0.7 Sr 0.3 MnO 3 /SrTiO 3 superlattices with medium oxygen pressure during SrTiO 3 deposition, and then the ferromagnetic charge-ordering (CO)/charge-density-wave (CDW) state was formed in La 0.7 Sr 0.3 MnO 3 with susceptibility to magnetic field, temperature, and superlattice configuration. In addition, re-entrance to the ferromagnetic insulator CO/CDW could be induced by a high magnetic field. The selective occupancy of oxygen vacancies on the (001) La/SrO plane and the phase fluctuation around the CO/CDW-to-ferromagnetic metal phase transition were critical to the above properties. The interfacial diffusion with atomic-level resolution to control oxygen vacancies could be used in other magnetoelectric systems and corresponding device applications.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗