Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic ordering”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
NASA NTRS1961

Stability of Titanium Nitride and Titanium Carbide When Exposed to Hydrogen Atoms from 298 to 1950 K

Titanium nitride and titanium carbide deposited on tungsten wires were exposed to hydrogen atoms (10(exp -4) atm pressure) produced by the action of microwave radiation on molecular hydrogen. The results of these experiments in the temperature range 298 to 1950 K indicate that no appreciable reaction takes place between atomic hydrogen and TiN or TiC. The formation of reaction products (NH3, CH4, C2H2) should be favored at lower temperatures. However, because of the high catalytic activity of Ti for H atom recombination, the rate of such reactions with H atoms is controlled by the rate of evaporation of Ti from the surface, this rate being low at temperatures below 1200 K. In order to interpret the stability of TiN and TiC in H atoms more fully, the stability of TiN and TiC in vacuum and H2 gas was also studied. The thermodynamic computations conform in order of magnitude to the experimentally found rates of decomposition of TiN and TiC in vacuum and are also consistent with the fact that no appreciable reaction is found with these compounds in molecular H2 at a pressure of 10(exp -3) atmosphere in the temperature range 2980 to 2060 K. When TiN or TiC was heated in atomic H or molecular H2, no reaction products other than those obtained from the simple decomposition of the nitride and carbide were observed. The gaseous products were analyzed in a mass spectrometer.

Philipp, Warren H.↗
DOE OSTI2024

First-principles and cluster expansion study of the effect of magnetism on short-range order in Fe–Ni–Cr austenitic stainless steels

Short-range order (SRO), the regular and predictable arrangement of atoms over short distances, alters the mechanical properties of technologically relevant structural materials such as medium/high entropy alloys and austenitic stainless steels. In this study, we present a generalized spin cluster expansion (CE) model and show that magnetism is a primary factor influencing the level of SRO present in austenitic Fe-Ni-Cr alloys. The spin CE consists of a chemical cluster expansion combined with an Ising model for Fe-Ni-Cr austenitic alloys. It explicitly accounts for local magnetic exchange interactions, thereby capturing the effects of finite temperature magnetism on SRO. Model parameters are obtained by fitting to a first-principles data set comprising both chemically and magnetically diverse FCC configurations. The magnitude of the magnetic exchange interactions are found to be comparable to the chemical interactions. Compared to a conventional implicit magnetism CE built from only magnetic ground state configurations, the spin CE shows improved performance on several experimental benchmarks over a broad spectrum of compositions, particularly at higher temperatures due to the explicit treatment of magnetic disorder. We find that SRO is strongly influenced by alloy Cr content, since Cr atoms prefer to align antiferromagnetically with nearest neighbors but become magnetically frustrated with increasing Cr concentration. Using the spin CE, we predict that increasing the Cr concentration in typical austenitic stainless steels promotes the formation of SRO and increases order-disorder transition temperatures. Furthermore, this study underscores the significance of considering magnetic interactions explicitly when exploring the thermodynamic properties of complex transition metal alloys. It also highlights guidelines for customizing SRO through adjustments of alloy composition.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Impact of chemical ordering on thermodynamic properties of point defects and Xe substitutional in U-10Mo

The accurate knowledge of defect energetics is critical to understanding the aging and irradiation behavior of U-10Mo nuclear fuel, which is selected as the fuel type for conversion of the United States High-Performance Research Reactors (HPRRs). Furthermore, using hybrid molecular dynamics and Monte Carlo (MDMC) simulation, we studied the impact of chemical ordering on the formation energies of vacancies, interstitials, and the solution energy of the Xe substitutionals. Instead of forming a random solid solution (RSS), substantial short-range-order (SRO) develops in U-10Mo, particularly at low temperatures. Mo atoms are found to repel each other and prefer U-rich local atomic environments within the 1st nearest neighbor (1NN) cutoff. Compared to the case of a RSS, the state with equilibrated Mo ordering shifts the distributions of vacancy and interstitial formation energies due to the dependence of defect energies on the local atomic environment, without a clear effect on Xe solution energy. In the operation temperature range (100–250 °C) of U-Mo fuels, neglecting SRO can lead to an inaccurate estimate of thermal equilibrium point defect concentrations by over an order of magnitude and incorrectly predict the preference among different types of dumbbells, highlighting the critical importance of accounting for the impact of chemical ordering for accurate atomistic calculations of defect properties.

Molecular dynamics and Monte Carlo (MDMC)↗
DOE OSTI2025

Thermodynamics and kinetics of core-shell versus appendage co-precipitation morphologies

What determines precipitate morphologies in co-precipitating alloy systems? We focus on alloys of two precipitating phases, with the fast-precipitating phase acting as heterogeneous nucleation sites for a second phase manifesting slower kinetics. Kinetic lattice Monte Carlo simulations show that the interplay between interfacial and ordering energies, plus active diffusion paths, strongly affect the selection of core-shell verses appendage morphologies. We study a FeCuMnNiSi alloy using the combination of atom probe tomography and simulations, and show that the ordering energy reduction of the MnNiSi phase heterogeneously nucleated on a pre-existing copper-rich precipitate exceeds the energy penalty of a predominantly Fe/Cu interface, leading to initial appendage, rather than core-shell, formation. Diffusion of Mn, Ni and Si around and through the Cu core towards the ordered phase results in subsequent appendage growth. We further show that in cases with higher primary precipitate interface energies and/or suppressed ordering, the coreshell morphology is favored.

Shu, Shipeng↗
DOE OSTI2025

Thermodynamic assessment of the quaternary WTaCrV refractory high entropy alloy as a means to guide experimental approaches

The deployment of fusion energy poses challenges for materials in plasma facing components to withstand high temperatures and thermal gradients, particle implantation and neutron damage. The current material of choice is tungsten, although property degradation limits its consideration in future fusion reactors. Hence, materials with better resistance to harsh environments need to be developed for fusion energy to become a reality. High entropy alloys are being explored as potential candidates with some compositions showing good radiation resistance to defect cluster formation. One of these materials is the WTaCrV system, although only one composition has been tested under ion irradiation. In this work, we study the thermodynamic properties of the entire quaternary alloy composition range. Coupling first principles calculations, cluster expansion approaches, and Monte Carlo methods, we access the free energy functionals, short-range ordering as a function of temperature, and atomic configurations that can be compared to experimental observations. We use this data to inform experiments into compositions with higher propensity to form solid solutions, instead of phase separating. With this formalism we have developed thermodynamic database (TDB) files that can be used to plot quaternary phase diagrams.

Cluster Expansion↗
NASA NTRS2003

Raman-atomic force microscopy of the ommatidial surfaces of Dipteran compound eyes

The ommatidial lens surfaces of the compound eyes in several species of files (Insecta: Diptera) and a related order (Mecoptera) were analyzed using a recently developed Raman-atomic force microscope. We demonstrate in this work that the atomic force microscope (AFM) is a potentially useful instrument for gathering phylogenetic data and that the newly developed Raman-AFM may extend this application by revealing nanometer-scale surface chemistry. This is the first demonstration of apertureless near-field Raman spectroscopy on an intact biological surface. For Chrysopilus testaceipes Bigot (Rhagionidae), this reveals unique cerebral cortex-like surface ridges with periodic variation in height and surface chemistry. Most other Brachyceran flies, and the "Nematoceran" Sylvicola fenestralis (Scopoli) (Anisopodidae), displayed the same morphology, while other taxa displayed various other characteristics, such as a nodule-like (Tipula (Triplicitipula) sp. (Tipulidae)) or coalescing nodule-like (Tabanus punctifer Osten Sacken (Tabanidae)) morphology, a smooth morphology with distinct pits and grooves (Dilophus orbatus (Say) (Bibionidae)), or an entirely smooth surface (Bittacus chlorostigma MacLachlan (Mecoptera: Bittacidae)). The variation in submicrometer structure and surface chemistry provides a new information source of potential phylogenetic importance, suggesting the Raman-atomic force microscope could provide a new tool useful to systematic and evolutionary inquiry.

Eye/ultrastructure↗
DOE OSTI2022

Experimentally informed structure optimization of amorphous TiO 2 films grown by atomic layer deposition

Amorphous titanium dioxide TiO 2 (a-TiO 2 ) has been widely studied, particularly as a protective coating layer on semiconductors to prevent corrosion and promote electron–hole conduction in photoelectrochemical reactions. The stability and longevity of a-TiO 2 is strongly affected by the thickness and structural heterogeneity, implying that understanding the structure properties of a-TiO 2 is crucial for improving the performance. This study characterized the structural and electronic properties of a-TiO 2 thin films (~17 nm) grown on Si by atomic layer deposition (ALD). Fluctuation spectra V(k) and angular correlation functions were determined with 4-dimensional scanning transmission electron microscopy (4D-STEM), which revealed the distinctive medium-range ordering in the a-TiO 2 film. A realistic atomic model of a-TiO 2 was established guided by the medium-range ordering and the previously reported short-range ordering of a-TiO 2 film, as well as the interatomic potential. The structure was optimized by the StructOpt code using a genetic algorithm that simultaneously minimizes energy and maximizes the match to experimental short- and medium-range ordering. The StructOpt a-TiO 2 model presents improved agreements with the medium-range ordering and the k-space location of the dominant 2-fold angular correlations compared with a traditional melt-quenched model. The electronic structure of the StructOpt a-TiO 2 model was studied by ab initio calculations and compared to the crystalline phases and experimental results. Finally, this work uncovered the medium-range ordering in a-TiO 2 thin films and provided a realistic a-TiO 2 structure model for further investigation of structure–property relationships and materials design. In addition, the improved multi-objective optimization package StructOpt was provided for structure determination of complex materials guided by experiments and simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Iron Carbonyl Complexes of [2.2.2]Hericene as a Rigid Tris(1,3-diene) Ligand

Hericene is an unusual hexaolefin consisting of three 1,3-diene units located on a rigid bicyclo [2.2.2]octane framework that restricts the geometrical relationships of metal atoms bonded to these olefinic units. In order to explore possible effects of this rigidity limiting metal–metal interaction in polynuclear derivatives possibly stabilizing coordinatively unsaturated species, the structures and energetics of the hericene iron carbonyl complexes (hericene)Fem(CO)n (m = 1, n = 3; m = 2, n = 6, 5; m = 3, n = 9, 8) have been investigated by density functional theory. The lowest-energy (hericene)Fem(CO)3m (m = 1, 2, 3) structures have the cavities of the hericene ligand filled with a single Fe(CO)3 moiety bonded to a 1,3-diolefin unit. Such species have been synthesized by the reaction of Fe2(CO)9 with hericene. For the (hericene)Fe2(CO)5 system, the lowest energy structures are singlet structures with an Fe(CO)3 unit bonded to a 1,3-diene unit in one hericene cavity and an Fe(CO)2 unit in another hericene cavity bonded to three C=C double bonds from two 1,3-diene units. Higher energy (hericene)Fe2(CO)5 structures include a structure in which a single hericene cavity contains a Fe2(CO)4(µ-CO) moiety with each iron atom bonded to a 1,3-diene unit. In addition, both singlet and triplet (hericene)Fe2(CO)5 structures are found in which an Fe(CO)3 moiety and an Fe(CO)2 moiety located in separate hericene cavities are each bonded to a 1,3-diene unit. The lowest-energy (hericene)Fe3(CO)8 structures have two hericene cavities containing Fe(CO)3 moieties fully bonded to 1,3-diene units and a third hericene cavity containing an Fe(CO)2 moiety fully bonded to a 1,3-diene unit.

Luo, Jinfeng↗
DOE OSTI2024

The derivation of CRSS in pure Ti and Ti-Al alloys

The work focuses on the determination of the critical resolved shear stress (CRSS) in titanium (Ti) and titanium-aluminum (Ti-Al) alloys, influenced by an array of factors such as non-symmetric fault energies and minimum energy paths, dislocation core-widths, short-range order (SRO) effects which alter the local atomic environment, and tension-compression (T-C) asymmetry affected by intermittent slip motion. To address these multifaceted complexities, an advanced theory has been developed, offering an in-depth understanding of the mechanisms underlying slip behavior. The active slip systems in these materials are basal, prismatic, and pyramidal planes, with the latter involving both $\langle$a$\rangle$ and $\langle$c + a$\rangle$ dislocations. Each slip system is characterized by distinct Wigner-Seitz cell configurations for misfit energy calculations, varying partial dislocation separation distances, and unique dislocation trajectories—all critical to precise CRSS calculations. The theoretical CRSS results were validated against a comprehensive range of experimental data, demonstrating a strong agreement and underscoring the model's efficacy.

Critical stress↗
DOE OSTI2025

AUTOIGNITION DELAY TIMES FOR REFORMATE GAS MIXTURES FROM METHANE GAS ENGINES

Methane slip is a prominent issue in natural gas reciprocating engines that are used in transportation and marine applications. The incomplete combustion that results in methane slip can be resolved with the introduction of hydrogen within the combustion mixture to improve methane oxidation and further enable combustion within the engine crevices where methane has previously remained unreacted. Steam methane reforming (SMR) is a common method used to produce hydrogen and can be used to design an onboard device to reduce methane slip from reciprocating engines. The development of this reformer device requires the validation of high-fidelity chemical kinetic models at the low temperatures of the crevice volumes of these engines. In this work, auto-ignition data is obtained using a shock tube at lean (φ—0.714 or λ—1.4) and stoichiometric (φ, λ = 1) equivalence ratios spanning a temperature range of 1042–1234 K at the 80-bar operating pressure of the test engine. Blends of methane, hydrogen, and reformate products from the SMR reaction are shock-heated in synthetic air, with the ignition delay time measured using an OH* chemiluminescence detector at 310 nm and a CH* detector at 430 nm. The experimental results are compared to several state-of-the-art chemical kinetic mechanisms from the literature. In general, most of the mechanisms show very good agreement with experiments at higher temperatures, with simulation results showing little deviation from experiments at lower temperatures. A sensitivity analysis was conducted, and the results reveal that the reaction H2 + CH3O2 = H + CH3O2H has a very significant role in determining low-temperature ignition delay times (IDTs) of SMR mixtures. These findings provide valuable insights into the chemical kinetics governing methane reformate combustion and contribute to the optimization of onboard reformer designs aimed at mitigating methane slip in natural gas-fueled engines.

Fraze, Matthew↗
DOE OSTI2023

Manipulation of Scattering Spectra with Topology of Light and Matter

Structured lights, including beams carrying spin and orbital angular momenta, radially and azimuthally polarized vector beams, as well as spatiotemporal optical vortices, have attracted significant interest due to their unique amplitude, phase front, polarization, and temporal structures, enabling a variety of applications in optical and quantum communications, micromanipulation, and super-resolution imaging. In parallel, structured optical materials, metamaterials, and metasurfaces consisting of engineered unit cells—meta-atoms, opened new avenues for manipulating the flow of light and optical sensing. While several studies explored structured light effects on the individual meta-atoms, their shapes are largely limited to simple spherical geometries. However, the synergy of the structured light and complex-shaped meta-atoms has not been fully explored. Here, in this paper, the role of the helical wavefront of Laguerre–Gaussian beams in the excitation and suppression of higher-order resonant modes inside all-dielectric meta-atoms of various shapes, aspect ratios, and orientations, is demonstrated and the excitation of various multipolar moments that are not accessible via unstructured light illumination is predicted. The presented study elucidates the role of the complex phase distribution of the incident light in shape-dependent resonant scattering, which is of utmost importance in a wide spectrum of applications ranging from remote sensing to spectroscopy.

42 ENGINEERING↗
DOE OSTI2021

Predictability as a probe of manifest and latent physics: The case of atomic scale structural, chemical, and polarization behaviors in multiferroic Sm-doped BiFeO 3

The predictability of a certain effect or phenomenon is often equated with the knowledge of relevant physical laws, typically understood as a functional or numerically derived relationship between the observations and known states of the system. Correspondingly, observations inconsistent with prior knowledge can be used to derive new knowledge on the nature of the system or indicate the presence of yet unknown mechanisms. In this work, we explore the applicability of Gaussian processes (GP) to establish predictability and uncertainty of local behaviors from multimodal observations, providing an alternative to this classical paradigm. Using atomic resolution scanning transmission electron microscopy (STEM) of multiferroic Sm-doped BiFeO 3 across a broad composition range, we directly visualize the atomic structure and structural, physical, and chemical order parameter fields for the material. GP regression is used to establish the predictability of the local polarization field from different groups of parameters, including the adjacent polarization values and several combinations of physical and chemical descriptors, including lattice parameters, column intensities, etc. We observe that certain elements of microstructure, including charged and uncharged domain walls and interfaces with the substrate, are best predicted with specific combinations of descriptors, and this predictability and associated uncertainties are consistent across the composition series. The associated generative physical mechanisms are discussed. It is also found that certain parameter combinations tend to predict the orthorhombic phase in the cases where rhombohedral phase is observed, suggesting a potential role of clamping and confinement phenomena in phase equilibrium in Sm-BiFeO 3 system close to morphotropic phase boundary. We argue that predictability and uncertainty in observational data offer a new pathway to probe the physics of condensed matter systems from multimodal local observations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Toward Automated Sampling of Polymorph Nucleation and Free Energies with the SGOOP and Metadynamics

Understanding the driving forces behind the nucleation of different polymorphs is of great importance for material sciences and the pharmaceutical industry. This includes understanding the reaction coordinate that governs the nucleation process and correctly calculating the relative free energies of different polymorphs. Here, we demonstrate, for the prototypical case of urea nucleation from the melt, how one can learn such a one-dimensional reaction coordinate as a function of prespecified order parameters and use it to perform efficient biased all-atom molecular dynamics simulations. The reaction coordinate is learnt as a function of the generic thermodynamic and structural order parameters using the “spectral gap optimization of order parameters (SGOOP)” approach [Tiwary, P. and Berne, B. J. Proc. Natl. Acad. Sci. U.S.A. (2016)] and is biased using well-tempered metadynamics simulations. The reaction coordinate gives insights into the role played by different structural and thermodynamics order parameters, and the biased simulations obtain accurate relative free energies for different polymorphs. This includes an accurate prediction of the approximate pressure at which urea undergoes a phase transition and one of the metastable polymorphs becomes the most stable conformation. Here, we believe the ideas demonstrated in this work will facilitate efficient sampling of nucleation in complex, generic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2017

Radar Detectability Studies of Slow and Small Zodiacal Dust Cloud Particles. III. The Role of Sodium and the Head Echo Size on the Probability of Detection

We present a path forward on a long-standing issue concerning the flux of small and slow meteoroids, which are believed to be the dominant portion of the incoming meteoric mass flux into the Earth's atmosphere. Such a flux, which is predicted by dynamical dust models of the Zodiacal Cloud, is not evident in ground-based radar observations. For decades this was attributed to the fact that the radars used for meteor observations lack the sensitivity to detect this population, due to the small amount of ionization produced by slow-velocity meteors. Such a hypothesis has been challenged by the introduction of meteor head echo (HE) observations with High Power and Large Aperture radars, in particular the Arecibo 430 MHz radar. Janches et al. developed a probabilistic approach to estimate the detectability of meteors by these radars and initially showed that, with the current knowledge of ablation and ionization, such particles should dominate the detected rates by one to two orders of magnitude compared to the actual observations. In this paper, we include results in our model from recently published laboratory measurements, which showed that (1) the ablation of Na is less intense covering a wider altitude range; and (2) the ionization probability, Beta ip, for Na atoms in the air is up to two orders of magnitude smaller for low speeds than originally believed. By applying these results and using a somewhat smaller size of the HE radar target we offer a solution that reconciles these observations with model predictions.

meteorites↗
DOE OSTI2022

CO oxidation on MgAl 2 O 4 supported Ir n : activation of lattice oxygen in the subnanometer regime and emergence of nuclearity-activity volcano

CO oxidation on Pt group metals is affected by the metal size and reducibility of the oxide support. In this study, we report that Ir supported on MgAl 2 O 4 , traditionally considered non-reducible, exhibits properties similar to reducible oxides when the Ir size is in the subnanometer regime. To show this effect, we synthesized subnanometer Ir clusters and compared their properties to single atoms and nanoparticles (1–1.5 nm). The CO oxidation activity is highest on Ir 0.6–0.8nm while showing distinctly different reaction orders in CO and O 2 (0, +0.4), than single atoms (1, 0) and nanoparticles (–1, +1). Microcalorimetry, in situ X-ray absorption, and infrared spectroscopies show that the CO-saturated Ir 0.6–0.8nm clusters could adsorb and activate O 2 despite binding CO more strongly than nanoparticles. Density functional theory calculations on CO saturated Ir 4 clusters suggest that the increased activity is due to the ability to activate O 2 on oxygen vacancies at the Ir–MgAl 2 O 4 interface. The findings show the important effect of the metal nuclearity on the support and catalyst properties and can guide future design of CO oxidation catalysts.

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

The effect of local chemical ordering on Frank-Read source activation in a refractory multi-principal element alloy

In this work, we investigate the operation of Frank-Read (FR) sources in a refractory multi-principal element alloy (MPEA). Simulations of discrete dislocation motion in MPEAs is enabled by the development of a phase field dislocation dynamics model that treats the atomic-scale fluctuations in lattice energies across the glide plane, present due to local ordering in the chemical composition within the nominally random MPEA atomic structure. We consider, through simulation, a range of length scales over which ordering occurs, varying from short-range lengths, a few times dislocation core width, to long-range lengths, an order of magnitude longer than the core. Characteristic of this body-centered cubic MPEA, the simulations also include screw/edge character dependence in glide resistance, as informed by atomic scale simulation. The critical stresses to activate the source for the same source size are found to be statistically distributed, as a direct consequence of the underlying variation in lattice energy. Furthermore, analysis of critical state for activating edge and screw FR sources in the MPEA reveals that FR source operation occurs via a two-step mechanism, involving athermal kink-pair formation, unlike the conventional FR source operation in a material with no composition fluctuations. This mechanism lowers the average critical stress required to activate the FR source and causes the statistical dispersion in critical stress to depend on the range of composition ordering. More importantly, it leads to a more severe dependence of source strength on FR source length than predicted by line tension alone.

36 MATERIALS SCIENCE↗
NASA NTRS1991

H2/O2 three-body rates at high temperatures

Hydrogen atoms are produced in the presence of excess O2, and the first-order decay are studied as a function of temperature and pressure in order to obtain the rate coefficient for the three-body reaction between H-atoms and O2. Attention is focused on the kinetic scheme employed as well as the reaction cell and photolysis and probe laser system. A two-photon laser-induced fluorescence technique is employed to detect H-atoms without optical-thickness or O2-absorption problems. Results confirm measurements reported previously for the H + O2 + N2 reaction at 300 K and extend these measurements to higher temperatures. Preliminary data indicate non-Arrehenius-type behavior of this reaction rate coefficient as a function of temperature. Measurements of the rate coefficient for H + O2 + Ar reaction at 300 K give a rate coefficient of 2.1 +/- 0.1 x 10 to the -32nd cm exp 6/molecule sec.

Marinelli, William J.↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗