Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion Monte Carlo”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Role of Short-Range Order on Diffusion Coefficients in the Li–Mg Alloy

Li–Mg alloys are important because of their beneficial role in fostering uniform plating and stripping of lithium in all-solid-state batteries. The alloy Li x Mg 1–x forms a solid solution on the BCC crystal structure when the lithium content is greater than x ≈ 0.3. The activation barriers of lithium and magnesium exchanges with a vacancy, crucial for substitutional diffusion, are predicted to be exceptionally low and almost identical, with negligible dependence on the alloy composition. The equilibrium vacancy concentration at room temperature is predicted to be very low, and it also remains almost constant with no dependence on Mg content in the alloy (for x Li ≥ 0.5). Nevertheless, both experiments and kinetic Monte Carlo simulations indicate that the tracer diffusion coefficients decrease by almost an order of magnitude with the addition of Mg to the alloy. In this contribution, the crucial role that chemical short-range order plays in affecting the diffusion coefficients is studied. Chemical short-range order is found to increase the effective activation barrier for lithium and magnesium diffusion by making successive atomic hops with vacancies correlated.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Migration energy barriers and diffusion anisotropy of point defects on tungsten surfaces

We report tungsten is one of the most promising candidates for plasma-facing materials in future fusion devices, owing to its high performance under extreme irradiation conditions. However, irradiation-induced surface morphology varies significantly, depending on the irradiation type, fluence, and flux. Therefore, it is critical to examine the dynamics of point defects on tungsten (W) surfaces to understand how irradiation affects surface morphology. In this study, we employ the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) method to search for potential migration paths of point defects on W (1 0 0), (1 1 0) and (1 1 1) surfaces. The first-principles calculations are then used to accurately determine the migration energy barriers. The obtained paths and barriers are incorporated into a kinetic Monte Carlo (KMC) model to determine trajectories and diffusivities, which are described by diffusion tensors to demonstrate their anisotropic features. Multiple diffusion mechanisms have been identified on different surfaces, with various anisotropy factors. Particularly, point defects on the W (1 1 0) surface have the highest diffusivities, with anisotropy factors independent of temperatures. In comparison, the anisotropy of the W (1 0 0) surface decreases as temperature increases, while the W (1 1 1) surface is isotropic for point defects. This study provides insights into defect transport properties on different surfaces, which are essential for understanding the early stages of irradiation-induced microstructural evolutions and surface morphology of tungsten.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Toward a Better Understanding of Ni Coarsening in Solid Oxide Cells: NiH on Ni (111) Examined Using a Combined Theoretical Approach

The coarsening of the Ni particles in the hydrogen electrode of solid oxide cells (SOCs) is an important degradation mechanism. Here, in this paper, density-functional theory and kinetic Monte Carlo methods are used to explore our recent hypothesis that the surface diffusion of NiH may cause faster Ni coarsening in electrolysis cell mode under an overpotential. Using both methods, the diffusion constant or diffusivity of NiH on Ni (111) is determined as the product of the surface coverage and single-molecule diffusivity for the first time considering all possible diffusion paths. It is then determined versus overpotential at the triple-phase boundary of the hydrogen electrode assuming a typical operating temperature of the SOC. Under a significant overpotential, the diffusivity of NiH is found to be sufficiently large to support the above hypothesis that NiH may promote Ni coarsening. However, based on the adsorption configurations identified, the dissociation and reformation of NiH on Ni (111) could occur. Thus, more work is needed to develop a model of Ni coarsening considering both molecular and dissociated forms of NiH.

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

Hydroxyl radical yields in the heavy ion radiolysis of water

The yields of hydroxyl radicals in the radiolysis of water with protons and carbon ions have been examined using experimental scavenging techniques coupled with Monte Carlo track simulations. Combined with previous results using helium ions, this set of data gives valuable information on the potential for radiation induced damage to biological systems for light ions with very different track structures. Carbon dioxide production from aerated formic acid solutions in concentrations ranging from 10 -3 to 1 M was used as a probe of hydroxyl radical yields from about 8 ns to 8 µs. Numerical interpolation of the results at slightly different ion energies in combination with data from gamma radiolysis allows for a systematic analysis of both track average and track segment yields. As expected, considerable track chemistry is found to occur on the nanosecond to microsecond time scales. Monte Carlo track simulations employing stochastic diffusion-kinetic calculations of product yields are found to reproduce experimental observations satisfactorily. The track simulations are used to extract hydroxyl radical kinetics in pure water at neutral conditions.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2023

Local inversion-symmetry breaking in a bismuthate high-Tc superconductor

Abstract The doped perovskite BaBiO 3 exhibits a maximum superconducting transition temperature ( T c ) of 34 K and was the first high- T c oxide to be discovered, yet pivotal questions regarding the nature of both the metallic and superconducting states remain unresolved. Although it is generally thought that superconductivity in the bismuthates is of the conventional s -wave type, the pairing mechanism is still debated, with strong electron-phonon coupling and bismuth valence or bond disproportionation possibly playing a role. Here we use diffuse x-ray scattering and Monte Carlo modeling to study the local structure of Ba 1-x K x BiO 3 across its insulator-metal boundary. We find no evidence for either long- or short-range disproportionation, which resolves a major conundrum, as disproportionation and the related polaronic effects are likely not relevant for the metallic and superconducting states. Instead, we uncover nanoscale structural correlations that break inversion symmetry, with far-reaching implications for the electronic physics. This unexpected finding furthermore establishes that the bismuthates belong to the broader classes of materials with hidden spin-orbit coupling and a tendency towards inversion-breaking displacements.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cluster spin glass correlations and dynamics in Zn 0.5⁢ Mn 0.5⁢ Te

Here, we present a combined magnetometry, muon spin-relaxation (𝜇⁢SR), and neutron-scattering study of the insulating spin glass Zn 0.5 ⁢Mn 0.5 ⁢Te, for which magnetic Mn 2+ and nonmagnetic Zn 2+ ions are randomly distributed on a face-centered cubic lattice. The magnetometry and 𝜇⁢SR results confirm a spin freezing transition around 𝑇 𝑓 ≈ 23 K, with the spin-fluctuation rate decreasing gradually and somewhat inhomogeneously through the sample volume as the temperature decreases toward 𝑇 𝑓 . Characteristic spin-correlation times well above 𝑇 𝑓 are on the order of 10 −10 s, much slower than typically observed in canonical spin glasses but in line with expectations for a cluster spin glass. Using magnetic pair distribution function (mPDF) analysis and reverse Monte Carlo (RMC) modeling of the magnetic diffuse neutron-scattering data, we show that the spin-glass ground state consists of clusters of spins exhibiting short-range-ordered type-III antiferromagnetic correlations with a locally ordered moment of 3.1⁢(1)⁢𝜇 B between nearest-neighbor spins. The type-III correlations decay exponentially as a function of spin separation distance with a correlation length of approximately 5 Å. The diffuse magnetic scattering and corresponding mPDF show no significant changes across 𝑇 𝑓 , indicating that the dynamically fluctuating short-range spin correlations in the paramagnetic state retain the same basic type-III configuration that characterizes the spin-glass state; the only change apparent from the neutron-scattering data is a gradual reduction of the correlation length and locally ordered moment with increasing temperature. Taken together, these results paint a unique and detailed picture of the local magnetic structure and dynamics in Zn 0.5 ⁢Mn 0.5⁢ Te and provide strong evidence that this material is best described as a cluster spin glass. In addition, this work showcases a statistical method for extracting diffuse scattering signals from neutron powder diffraction data, which we developed to facilitate the mPDF and RMC analysis of the neutron data. This method has the potential to be broadly useful for neutron powder diffraction experiments on a variety of materials with short-range atomic or magnetic order.

magnetism↗
DOE OSTI2021

Magnetic Diffuse Scattering of YBaCo 4 O 7 and LuBaCo 4 O 7 on Kagome and Triangular Lattices

In this work, the magnetic diffuse scattering of YBaCo 4 O 7 and LuBaCo 4 O 7 with alternating kagome and triangular lattices has been studied through neutron scattering measurements on single-crystal samples. In YBaCo 4 O 7 , the line-shape-magnetic-diffuse scattering is observed around the magnetic transition temperature, which is consistent with the Monte Carlo simulation. In LuBaCo 4 O 7 , the line-shape-magnetic-diffuse scattering has not been detected. At low temperature, the propagation vectors of the observed magnetic reflections are different between YBaCo 4 O 7 and LuBaCo 4 O 7 . The structural distortion in RBaCo 4 O 7 is suggested to play a critical role for the magnetic properties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Surface coverage dynamics for reversible dissociative adsorption on finite linear lattices

Dissociative adsorption onto a surface introduces dynamic correlations between neighboring sites not found in non-dissociative absorption. We study surface coverage dynamics where reversible dissociative adsorption of dimers occurs on a finite linear lattice. We derive analytic expressions for the equilibrium surface coverage as a function of the number of reactive sites, N, and the ratio of the adsorption and desorption rates. Using these results, we characterize the finite size effect on the equilibrium surface coverage. For comparable N’s, the finite size effect is significantly larger when N is even than when N is odd. Moreover, as N increases, the size effect decays more slowly in the even case than in the odd case. The finite-size effect becomes significant when adsorption and desorption rates are considerably different. These finite-size effects are related to the number of accessible configurations in a finite system where the odd-even dependence arises from the limited number of accessible configurations in the even case. We confirm our analytical results with kinetic Monte Carlo simulations. We also analyze the surface-diffusion case where adsorbed atoms can hop into neighboring sites. As expected, the odd-even dependence disappears because more configurations are accessible in the even case due to surface diffusion.

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

CEFR simulation using diffusion code system RAST-F

In this study, the RAST-F code system, which is based on the two-step approach that couples a multi-group cross-section generation Monte-Carlo (MC) code and a multi-group nodal diffusion solver, was used for the neutronic simulation of the CEFR start-up experiments. The numerical solution of the RAST-F system was verified against the full core MC solution MCS at all control rods fully inserted and withdrawn states. The RAST-F solution of the selected experimental simulations was compared against the measurement data. A good agreement between RAST-F and MCS solutions was observed with less than 120 pcm discrepancies and 1.2% root-mean-square error in terms of k{sub eff} and power distribution, respectively. Meanwhile, the RAST-F result agreed well with the experimental data within two-sigma of experimental uncertainty. The good agreement of these results indicates that RAST-F can be used to neutronic steady-state simulations for small core-size SFR, which was challenged to deterministic code systems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

An Ensemble Score Filter for Tracking High-Dimensional Nonlinear Dynamical Systems

We propose an ensemble score filter (EnSF) for solving high-dimensional nonlinear filtering problems with superior accuracy. A major drawback of existing filtering methods, e.g., particle filters or ensemble Kalman filters, is the low accuracy in handling high-dimensional and highly nonlinear problems. EnSF attacks this challenge by exploiting the score-based diffusion model, defined in a pseudo-temporal domain, to characterizing the evolution of the filtering density. EnSF stores the information of the recursively updated filtering density function in the score function, in stead of storing the information in a set of finite Monte Carlo samples (used in particle filters and ensemble Kalman filters). Unlike existing diffusion models that train neural networks to approximate the score function, we develop a training-free score estimation that uses mini-batch-based Monte Carlo estimator to directly approximate the score function at any pseudo-spatial-temporal location, which provides sufficient accuracy in solving high-dimensional nonlinear problems as well as saves tremendous amount of time spent on training neural networks. Another essential aspect of EnSF is its analytical update step, gradually incorporating data information into the score function, which is crucial in mitigating the degeneracy issue faced when dealing with very high-dimensional nonlinear filtering problems. High-dimensional Lorenz systems are used to demonstrate the performance of our method. EnSF provides surprisingly impressive performance in reliably tracking extremely high-dimensional Lorenz systems (up to 1,000,000 dimension) with highly nonlinear observation processes, which is a well-known challenging problem for existing filtering methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Shape Stability of Truncated Octahedral fcc Metal Nanocrystals

Metallic nanocrystals (NCs) can be synthesized with tailored nonequilibrium shapes to enhance desired properties, e.g., octahedral fcc metal NCs optimize catalytic activity associated with {111} facets. However, maintenance of optimized properties requires stability against thermal reshaping. Thus, we analyze the reshaping of truncated fcc metal octahedra mediated by surface diffusion using a stochastic atomistic-level model with energetic input parameters for Pd. Furthermore, the model describes NC thermodynamics by an effective nearest-neighbor interaction and includes a realistic treatment of diffusive hopping for undercoordinated surface atoms. Kinetic Monte Carlo simulation reveals that the effective barrier, E eff , for the initial stage of reshaping is strongly tied to the degree of truncation of the vertices in the synthesized initial octahedral shapes. This feature is elucidated via exact analytic determination of the energy variation along the optimal reshaping pathway at low-temperature (T), which involves transfer of atoms from truncated {100} vertex facets to form new layers on {111} side facets. Deviations from predictions of the low-T analysis due to entropic effects are more prominent for higher T and larger NC sizes.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A physics-based model of cladding wastage layer formation rate

In this report, lower length scale simulations to inform an engineering-scale model of fuel-cladding chemical interaction (FCCI) conducted under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in FY22 are described. An overall strategy for the implementation of the BISON model is described, and the past and current lower length scale work on FCCI formation is put into the context of this overall strategy. Density functional theory and kinetic Monte Carlo simulations are used to determine the diffusion coefficient of neodymium in iron. Density functional theory calculations are also used to parameterize a parabolic approximation for the dependence of free energy on composition for the intermetallic compound Fe 17 Nd 2 . A phase-field model of the Fe-Nd system was developed that includes the random solid solution body-centered cubic phase and the intermetallic Fe 17 Nd 2 . The model was used to simulate growth of the intermetallic layer in a diffusion couple and the results were compared to experiment. The model’s prediction of the parabolic growth constant is within 40 % of the experimental value. The model is also used to simulate the formation of a grain boundary attack layer that is observed in experiment. The simulations show enhanced Nd concentration along the grain boundaries, but do not clearly show formation of a Fe 17 Nd 2 layer

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Phase-Field Modeling of Materials Interfaces and Nanostructures

Nanostructured materials offer unique properties for a wide range of energy applications. Among various known processing methods, liquid metal dealloying (LMD)—the selective dissolution of a base alloy element into a metallic melt—has emerged as a powerful technique to produce a new class of nano-/meso-scale open porous and bicontinuous composite structures with ultra-high interfacial area. LMD has recently been complemented by the advent of vapor phase dealloying (VPD), a novel technique that exploits the selective evaporation of one element from a parent alloy containing elements with very different vapor pressures, thereby enabling to fabricate open nanoporous structures of various elements from less-noble metals to inorganic elements regardless of their chemical activity without requiring high LMD temperatures and chemical etching. Together LMD and VPD have greatly expanded the scope of dealloying techniques, limited by traditional electrochemical means to noble metals, and boosted the design of new functional and structural materials that combine a wide variety of elements. Topologically-connected open porous structures with ultra-high surface area allow mass transport within the structure while preserving structural integrity, enabling them to serve as catalysts, fuel cells, supercapacitors, or high-capacity battery materials. Bicontinuous composite structures in turn can display high strength and high ductility or superior radiation-damage resistance due to the ultra-high interface area between interpenetrating solid phases. This research program makes use of state-of-the-art computational methods to understand at a basic level the self-organizing dealloying process with main focus on dealloying kinetics and interfacial pattern formation at the dealloying front controlling initial structure size, topology, and phase compositions. Phase-field simulation studies of LMD focus on solid solutions, line compounds, and intermetallic systems that can form ternary composites by nucleation and growth of a new phase. Studies of VPD employ phase-field modeling and a hybrid method combining a kinetic Monte Carlo (KMC) model of evaporation and surface diffusion with molecular dynamics for vapor-phase transport inside nanopores. Simulations explore mechanisms of interface- and diffusion-controlled dealloying kinetics, both observed in VPD but not fundamentally understood. In addition, our newly developed multi-physics phase-field approach of large-volume-change phase transformations is being used to model novel 3D anode geometries based on dealloyed nanoporous structures including novel sandwiched graphene/Si/silica for highrate Li ion battery. Those studies are aimed at elucidating geometric design principles that improve mechanical stability. We expect this research to enhance the capability to tailor nano-/mesoscale structures for a wide range of energy-related materials applications and to yield further advances in computational methodologies that benefit a broad materials research community.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Tunable chemical complexity to control atomic diffusion in alloys

Abstract In this paper we report a new fundamental understanding of chemically-biased diffusion in Ni–Fe random alloys that is tuned/controlled by the intrinsic quantifiable chemical complexity. Development of radiation-tolerant alloys has been a long-standing challenge. Here we show how intrinsic chemical complexity can be utilized to guide the atomic diffusion and suppress radiation damage. The influence of chemical complexity is shown by the example of interstitial atom (IA) diffusion that is the most important defect in radiation effects. We use μs-scale molecular dynamics to reveal sluggish diffusion and percolation of IAs in concentrated Ni–Fe alloys. We develop a mean field diffusion model to take into account the effect of migrating defect energy properties on diffusion percolation, which is verified by a new kinetic Monte Carlo approach addressing detailed processes. We demonstrate that the local variations in the ground state energy of IA configurations in alloys, reflecting the chemical difference between alloying components, drives the percolation effects for atomic diffusion. Percolation, chemically-biased and sluggish diffusion are phenomena that are directly related to the chemical complexity intrinsically to multicomponent alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Towards robust autonomous impedance spectroscopy analysis: A calibrated hierarchical Bayesian approach for electrochemical impedance spectroscopy (EIS) inversion

Distribution-based analyses, such as the distribution of relaxation times (DRT) and the distribution of diffusion times (DDT), present model-free alternatives to equivalent circuit modeling for analysis of electrochemical impedance spectroscopy (EIS) data. However, reconstructing such distributions from noisy impedance data is an ill-posed problem that must be solved with specialized inversion algorithms, requiring careful control and tuning. Furthermore, most inversion algorithms developed to date can only solve problems of limited complexity. Herein, we present a new hierarchical Bayesian method for EIS inversion, leveraging efficient algorithms for optimization and Hamiltonian Monte Carlo (HMC) sampling to solve models of arbitrary complexity. We overcome the challenge of ad-hoc parameter tuning by encoding intrinsic characteristics of the DRT and DDT into flexible prior distributions and “pre-calibrating” the model to simulated data. This approach is versatile, highly robust to noise, and provides quantitative estimates of both the error structure of the data and the uncertainty in the recovered distributions. The model is validated with simulated data to demonstrate accurate recovery of the DRT and the DDT. The method also shows promise for simultaneous recovery of multiple distributions, raising the intriguing possibility of semi-autonomous EIS analysis and ad-hoc model construction. Finally, the practical utility of the method is illustrated with experimental data. Throughout, we draw comparisons to several recently published EIS inversion methodologies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Prediction of vacancy defect diffusion paths in high entropy alloys via machine learning on molecular dynamics data

Identifying the diffusion path of point defects is a critical step in understanding their evolution and the mechanisms of related phenomena. Defect diffusion occurs at small length and time scales, with impacts on material properties that may continue to evolve over ns to μs, ms, and the continuum scale (s, min, etc., and cm, m, etc.). The time scale accessible to molecular dynamics (MD) simulations is limited by small step sizes, typically in the fs range. Thus, surrogate models of MD simulations through machine learning (ML)-based algorithms are of great interest, especially for complex systems such as high entropy alloys (HEAs). In this work, dynamics governing vacancy migration in HEA were approximated with graph convolutional network (GCN) models as ansatzes for kinetic Monte Carlo (KMC) rate catalogs. Network design considered that diffusion in crystalline solids generally depends on interactions between defects and their immediate neighbor atoms. Graphs represented the vacancy surroundings, MD-generated trajectories provided training and comparison datasets, and unsupervised GCN models approximated interatomic dynamics governing vacancy migration in HEAs as ansatzes for KMC. A proof-of-concept model trained on MD data for the Fe, Ni, Cr, Co, and Cu HEA environment was used with two different neighbor interactions to assess the feasibility of training a GCN to predict vacancy defect transition rates in the HEA environment. The resulting setup rapidly generated MD-formatted synthetic trajectories based on dynamics learned from the MD training set, with a time acceleration of roughly two orders of magnitude and a similar diffusion coefficient to MD observations. Additionally, Nudged Elastic Band (NEB) calculations were performed on randomly generated FeNiCrCoCu HEA structures to determine vacancy migration barriers across nearest-neighbor sites. Transition probabilities for each jump, categorized by atomic type, were extracted from these calculations. NEB-based and GCN-based approaches led to similar outcomes.

Reimer, C↗
DOE OSTI2026

Multiscale study of helium diffusion in Ni-Cr alloys: Short-range trapping versus long-range channeling

Ni-Cr alloys are widely employed as structural materials in fast nuclear reactors but are vulnerable to high-temperature helium (He) embrittlement (HTHE) under fast neutron irradiation. A comprehensive understanding of He diffusion in Ni-Cr alloys, which governs the kinetics of HTHE, is therefore essential for developing resilient materials and preventing failure. In this work, we reveal the underlying mechanisms of He diffusion in pure Ni and Ni-Cr alloys by integrating density functional theory (DFT) with atomic kinetic Monte Carlo (AKMC) simulations. Our findings uncover a non-monotonic dependence of He diffusivity on Cr concentration, contradicting the monotonic trends predicted by DFT-parameterized theories. At low Cr concentrations, He diffusion is dominated by short-range trapping, characterized by multiple trapping sites and a distinct mechanism within the first nearest neighbor of Cr, differing from that in pure Ni. At high Cr concentrations, these local traps become interconnected, forming long-range fast diffusion channels that enhance He mobility. The competition between localized trapping and extended channeling results in a diffusivity that first decreases, then increases with rising Cr content. These atomic-scale insights offer critical guidance for the design of radiation-tolerant Ni-based alloys. Moreover, the combined DFT-AKMC methodology and the concept of random walker diffusion through interconnected energy basins present a broadly applicable framework for studying transport phenomena in disordered systems.

36 - MATERIALS SCIENCE↗
DOE OSTI2021

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↗