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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Tuning Electronic Properties of 2D Materials Using Metal Adsorbates: Cu at WTe 2 Edges

Two-dimensional materials exhibit properties promising for novel applications. Topologically protected states at their edges can be harnessed for use in quantum devices. We use ab initio simulations to examine properties of edges in 1T'-WTe 2 monolayers, known to exhibit topological order, and their interactions with Cu atoms. Comparison of (010)-oriented edges that have the same composition but different terminations shows that, as the number of Cu atoms increases, their thermodynamically-preferred arrangement depends on the details of the edge structure. Cu atoms aggregate into a cluster at the most stable edge; while the cluster is non-magnetic, it spin-polarizes the W atoms along the edge, which removes the topological protection. At the metastable edge, Cu atoms form a chain incorporated into the WTe 2 lattice; the topological state is preserved in spite of the dramatic edge restructuring. This suggests that exploiting interactions of metal species with metastable edge terminations can provide a path towards non-invasive interfaces.

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

Rate Controlling in Low-Temperature Standard NH 3 -SCR: Implications from Operando EPR Spectroscopy and Reaction Kinetics

Here, a series of seven Cu/SSZ-13 catalysts with Si/Al = 6.7 are used to elucidate key rate-controlling factors during low-temperature standard ammonia-selective catalytic reduction (NH 3 -SCR), via a combination of SCR kinetics and operando electron paramagnetic resonance (EPR) spectroscopy. Strong Cu-loading-dependent kinetics, with Cu atomic efficiency increasing nearly by an order of magnitude, is found when per chabazite cage occupancy for Cu ion increases from ~0.04 to ~0.3. This is due mainly to the release of intercage Cu transfer constraints that facilitates the redox chemistry, as evidenced from detailed Arrhenius analysis. Operando EPR spectroscopy studies reveal strong connectivity between Cu-ion dynamics and SCR kinetics, based on which it is concluded that under low-temperature steady-state SCR, kinetically most relevant Cu species are those with the highest intercage mobility. Transient binuclear Cu species are mechanistically relevant species, but their splitting and cohabitation are indispensable for low-temperature kinetics.

catalytic reactions↗
DOE OSTI2022

Ab initio calculation of X-ray and related core-level spectroscopies: Green's function approaches

X-Ray and related spectroscopies are powerful probes of atomic, vibrational, and electronic structure. In order to unlock the full potential of such experimental techniques, accurate and efficient theoretical and computational approaches are essential. Here we review the status of a variety of first-principles and nearly first principles techniques for X-ray spectroscopies such as X-ray absorption, X-ray emission, and X-ray photoemission, with a focus on Green's function based methods. In particular, we describe the current state of multiple scattering Green's function techniques available in the FEFF10 code and cumulant Green's function techniques for including the effects of many-body electronic excitations. Furthermore, illustrative examples are shown for a variety of materials and compared with other theoretical and experimental results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

First-principles investigation of structure-property relationships in stable and metastable MXenes

Understanding the structure–property relationships in layered transition-metal carbides or nitrides, known as MXenes, is of critical importance for their rational design, synthesis, and application. However, the vast chemical and structural diversity of MXenes, stemming from their wide range of M and X elements, surface terminations, and different atomic coordination environments, makes it challenging to clearly understand these structure–property relationships. Here, in this work, we perform first-principles density functional theory (DFT) calculations and molecular dynamics (MD) simulations to comprehensively investigate the stability and a variety of physical properties of MXenes with different coordination environments. Using Ti- and Mo-based carbide MXenes as model systems, energetic calculations reveal that Ti-based MXenes are most stable in octahedral coordination, whereas Mo-based MXenes preferentially adopt prismatic coordination. This fundamental difference in preferred atomic coordination gives rise to markedly distinct properties between these two systems as a function of the fraction of octahedral and prismatic sites. For instance, the in-plane stiffness of Ti-based MXenes increases as octahedral coordination becomes dominant, but it decreases in the Mo-based MXenes under the same conditions. Additional stability analyses based on mechanical, lattice-dynamical, and temperature-dependent thermodynamic properties demonstrate that many metastable MXenes not only satisfy the strict stability criteria but can also undergo phase transitions among different structures and even become stabilized at elevated temperatures. Although surface terminations, such as F and O atoms, do not alter the energetic ordering or the overall stiffness trends among stable and metastable MXenes, they influence other material properties. For instance, O termination can induce semiconducting behavior in both stable and metastable Ti 2 ⁢CO 2 MXenes. This study significantly advances the fundamental understanding of structure–property relationships in MXenes and provides valuable guidance for developing coordination-based design principles to precisely engineer MXenes with improved properties.

Oyeniran, Noah [University of Alabama, Tuscaloosa,↗
DOE OSTI2022

Ferroelectric HfO 2 and the importance of strain

Ferroelectric oxides based on HfO 2 show tremendous promise for the next generation of memory and logic devices. The ferroelectric polymorph is one of several that can be derived from the high symmetry cubic fluorite structure of HfO 2 . A single grain of HfO 2 may consist of a coherent mixture of multiple orientational and translational variants of different polymorphs. Here, we use symmetry-adapted strain-order parameters to elucidate the relationship between the different HfO 2 polymorphs and their symmetrically equivalent variants. We use first-principles electronic structure methods to identify minimum energy pathways and map them in subspaces of the symmetry-adapted strain order parameters. We next investigate the atomic structure of domain boundaries that separate coexisting variants of ferroelectric HfO 2 . Further, we rely on Gibbsian excess quantities and a precise specification of mechanical boundary conditions to describe the thermodynamic properties of domain boundaries. Our first-principles calculations show that the O and Hf shuffle arrangement within a domain boundary is closely related to the intermediate shuffle patterns of the homogeneous pathways between ferroelectric variants. Furthermore, the preferred structure within a boundary is very sensitive to local strain constraints imposed by the adjacent ferroelectric variants, leading to highly anisotropic domain boundary energies.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Prediction of defect properties in concentrated solid solutions using a Langmuir-like model

The alleged existence of sluggish diffusion in high-entropy alloys has drawn controversy. In high-entropy alloys and, in general, in all solids, transport properties are controlled by point defect concentration, which must be known before performing atomistic simulations to compute transport coefficients. In this work, we present a general Langmuir-like model for defect concentration in an arbitrarily complex solid solution and apply this model to generate expressions for concentrations of vacancies and small interstitial atoms. We then calculate the vacancy concentration as a function of temperature in the equiatomic CoNiCrFeMn and FeAl alloys with modified embedded-atom-method potentials for various chemical orderings, showing there is no clear correlation between vacancy thermodynamics and chemical ordering in the CoNiCrFeMn alloy, but clear systematic patterns for FeAl. We believe this is due to the high stability of disordered, random, and ordered intermetallic phases, respectively, in the CoNiCrFeMn and FeAl systems. Finally, this work provides future avenues to the prediction of thermal interstitials and vacancies in solid solutions, which is necessary for models of nonequilibrium behavior of solid solutions.

composition↗
DOE OSTI2023

Simulation of the dynamics of gas mixtures during plasma processing in the C75 Cavity

Plasma processing using a mixture of noble gas and oxygen is a technique that is currently being used to reduce field emission and multipacting in accelerating cavities. Plasma is created inside the cavity when the gas mixture is exposed to an electromagnetic field that is generated by applying RF power through the fundamental power or higher-order mode couplers. Oxygen ions and atomic oxygen are created in the plasma which breaks down the hydrocarbons on the surface of the cavity and the residuals from this process are removed as part of the process gas flow. Removal of hydrocarbons from the surface increases the work function and reduces the secondary emission coefficient. This work describes the initial results of plasma simulation, which provides insight into the ignition process, distribution of different species, and interactions of free oxygen and oxygen ions with the cavity surfaces. The simulations have been done with an Ar/?2 plasma using COMSOL® multiphysics. These simulations help in understanding the dynamics and control of plasma inside the cavity and the exploration of different gas mixtures.

Dhakal, Pashupati↗
DOE OSTI2020

MOCVD Growth and Characterization of ZnSnN 2

Since the first few reports of its synthesis in 2011 and 20121-3, interest in ZnSnN 2 for diverse applications, including photovoltaics1,3,4, visible light photocatalysis5, and optoelectronics, has been growing. Novel nitride-based LED device designs promise high gains in efficiency at emission wavelengths into the green, amber and even red.6-8 Their realization requires the precise insertion of an ultrathin layer of ZnGeN 2 or ZnSnN 2 into the active region of an InGaAs-GaAs LED structure. While thin films of ZnSnN 2 have been synthesized by RF sputtering1,9 and by MBE 2,10,11, and thin films of ZnGeN 2 , ZnSiN 2 and their alloys by MOCVD, the synthesis of ZnSnN 2 thin films by MOCVD has not yet been reported. It is a challenging goal, given the narrow temperature window between efficient activation of ammonia and the thermal decomposition of ZnSnN 2 at near-atmospheric pressures. We report here the results of studies of ZnSnN 2 films grown on r- and c-plane sapphire substrates, GaN templates on sapphire, and In 0.15 Ga 0.85 N templates on GaN-on-sapphire, at temperatures from 500 °C to 560 °C, using tetramethyltin (TMSn), diethylzinc (DEZn), germane and ammonia. Chamber pressures were varied between 200 and 500 torr. In addition, the ratios of the cation precursors and the total cation-to-anion precursor ratios were varied. Single crystal growth, as evidenced by 2θ-ω x-ray diffraction scans, was achieved on all of these substrates. Wurtzitic (0002) ZnSnN 2 diffraction were observed at 2θ ~ 32.90 for the samples grown on c-sapphire, GaN (0001) and InGaN (0001) templates, and wurtzitic (11-20) peak at 2θ ~ 54.80 for the films grown on r-sapphire. Higher pressures and higher temperatures promoted better crystallinity and more continuous films, as revealed by scanning electron microscopy (SEM) images and 2θ-ω x-ray diffraction peak widths and intensities. Crystal quality was observed to be highly dependent on growth temperature. For example, in one series a 10 °C change in growth temperature resulting in a factor of 15 increase in diffraction peak intensity. SEM images show that films grown on the GaN templates had better uniformity in nucleation than those grown on sapphire. The presence of satellite peaks around the Bragg peak in the XRD 2θ-ω scan for a film grown on a GaN template, and the absence of satellite peaks for the films grown simultaneously on c- and r-sapphire, also indicate superior quality for the film grown on GaN. HAADF-STEM images of films grown on these templates showed smooth interfaces at the film-substrate boundary and show growth rates as high as 167 nm/h, considerably higher than reported previously for MBE growth.9,10 Electron diffraction patterns show single crystallinity. Rocking curves show FWHM as low as 0.100. RMS surfaces roughnesses, measured by atomic force microscopy, are of the order of 3 nm for 1µm x 1µm areas. Current work is focusing on optimizing growth on InGaN templates. These templates are promising for their relatively lower lattice mismatch with ZnSnN 2 . Their use is a next step toward producing LED device structures.

Jayatunga, Benthara Hewage Dinushi↗
DOE OSTI2020

Intermediate-Term Thermal Aging Effect Evaluation for Grade 92 and 316L at The LWR Relevant Temperature

Life extension of the existing nuclear reactors imposes accumulated damages, such as higher fluences and longer period of corrosion, to structural materials, which would result in significant challenges to the traditional reactor materials such as type 304 and 316 stainless steels. Advanced alloys with superior radiation resistance will increase safety margins, design flexibility, and economics for not only the life extension of the existing fleet but also new builds with advanced reactor designs. The Electric Power Research Institute (EPRI) teamed up with Department of Energy (DOE) on the Advanced Radiation Resistant Materials (ARRM) program, aiming to develop and test degradation resistant alloys from current commercial alloy specifications by 2021 to a new advanced alloy with superior degradation resistance in light water reactor (LWR)-relevant environments by 2024. Based on a comprehensive microstructure and property screening performed in Phase-1 of the ARRM program, a total of five alloys, together with 316L and X-750 as references, were down-selected for Phase-2 neutron irradiation studies. Because thermal aging could exert a synergistic effect on neutron irradiation due to the low neutron damage rate on the order of 10 –7 displacements per atom per second (dpa/s), Grade 92 (one of the five down-selected alloys) and two heats of 316L were selected in this task to study the effect of aging at 350°C for 12.5–12.7 kh on microstructure and mechanical properties. The aging time is approximately corresponding to 5 dpa neutron irradiation, which can be served as a reference for the 5-dpa-irradiated samples of the alloys to help understand the independent neutron damage influence on microstructure and mechanical properties. Optical microscopy and scanning electron microscopy were used for microstructural characterization. Hardness, tensile, Charpy impact toughness, and fracture toughness in the ductile regime of the aged samples were examined. The aging of Grade 92 led to the formation of many Laves phase in sizes of ~100–200 nm and resulted in some reduction in hardness and yield/ultimate tensile strength with some increases in uniform and total plastic elongations, which may have helped the enhancement of Charpy impact toughness, e.g., ~4 J increase in upper-shelf energy and 20.5°C reduction in ductile-brittle transition temperature compared with the unaged condition. The fracture toughness of the aged Grade 92 showed decent toughness of ~303 MPa√m (K Jq ) with ~87 tearing modulus at 22°C, which reduced to ~241 MPa√m (~20% reduction) with ~75 tearing modulus at 300°C. The two heats of 316L are differentiated by their amounts of δ-ferrite, i.e., ~1 vol% in heat T1103 and ~4 vol% in heat N5B8. The two heats of 316L were subjected to 15% cold work (CW) as this condition is often used in nuclear reactors to help trapping radiation-induced defects and thus postpone the steady state swelling stage. The aging did not result in any noticeable microstructure changes, except for possible segregations at grain boundaries that need to be further investigated. The aging resulted in slight reduction in hardness of 316L-T1103 and slight increase in hardness of 316L-N5B8, which is consistent with their tensile testing results. The 15%CW might have introduced some inhomogeneity, resulting in large standard deviations in hardness, which were reduced after the aging. The aged 316L followed the same trend as the unaged condition for the strength and elongation results from the tensile testing. Minimum elongations appeared within ~300–600°C, above which the elongations seem to increase. Unlike the negligible or minor changes in microstructure, hardness and tensile results of the two heats of 316L, their fracture toughness showed noticeable difference. The aged 316L-T1103 showed good toughness of ~350 MPa√m with ~80 tearing modulus at 22°C, which reduced to ~271 MPa√m (~23% reduction) with ~83 tearing modulus at 300°C. In contrast, the aged 316L-N5B8 showed lower toughness of ~244 MPa√m with ~53 tearing modulus at 22°C, which reduced to ~180 MPa√m (~26% reduction) with ~31 tearing modulus at 300°C. The preliminary results indicate that the presence of high volume of δ-ferrite would noticeable impair fracture toughness although it may not influence hardness and tensile properties.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Multipole Expansion of Atomic Electron Density Fluctuation Interactions in the Density-Functional Tight-Binding Method

The accuracy of the density-functional tight-binding (DFTB) method in describing noncovalent interactions is limited due to its reliance on monopole-based spherical charge densities. In this study, we present a multipole-extended second-order DFTB (mDFTB2) method that takes into account atomic dipole and quadrupole interactions. Furthermore, we combine the multipole expansion with the monopole-based third-order contribution, resulting in the mDFTB3 method. To assess the accuracy of mDFTB2 and mDFTB3, we evaluate their performance in describing noncovalent interactions, proton transfer barriers, and dipole moments. Here, our benchmark results show promising improvements even when using the existing electronic parameters optimized for the original DFTB3 model. Both mDFTB2 and mDFTB3 outperform their monopole-based counterparts, DFTB2 and DFTB3, in terms of accuracy. While mDFTB2 and mDFTB3 perform comparably for neutral and positively charged systems, mDFTB3 exhibits superior performance over mDFTB2 when dealing with negatively charged systems and proton transfers. Overall, the incorporation of the multipole expansion significantly enhances the accuracy of the DFTB method in describing noncovalent interactions and proton transfers.

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

Signatures of linearized gravity in atom interferometers: A simplified computational framework

We develop a general framework for calculating the leading-order, general relativistic contributions to the gravitational phase shift in single-photon atom interferometers within the context of linearized gravity. We show that the atom gradiometer observable, which only depends on the atom interferometer propagation phase, can be written in terms of three distinct contributions: the Doppler phase shift, which accounts for the tidal displacement of atoms along the baseline, the Shapiro phase shift, which accounts for the delay in the arrival time of photons at atom-light interaction points, and the Einstein phase shift, which accounts for the gravitational redshift measured by the atoms. For specific atom gradiometer configurations, we derive the signal and response functions for two physically motivated scenarios: (i) transient gravitational waves in the transverse-traceless gauge and, for the first time, in the proper detector frame, and (ii) transient massive objects sourcing weak and slow-varying Newtonian potentials. Here, we find that the Doppler contribution of realistic Newtonian noise sources (e.g., a freight truck or a piece of space debris) at proposed atom gradiometer experiments, such as AION, MAGIS, and AEDGE, can exceed the shot noise level and thus affect physics searches if not properly subtracted. Note: 44 pages including appendices, 4 figures; v2: different formatting, fixed typos, 28 pages including appendices, 5 figures, agrees with published version

Atom interferometry↗
DOE OSTI2023

Spin-mediated Mott excitons

Motivated by recent experiments on Mott insulators, in both iridates and ultracold atoms, we theoretically study the effects of magnetic order on the Mott-Hubbard excitons. In particular, we focus on spin-mediated doublon-holon pairing in Hubbard materials. Here, we use several complementary theoretical techniques: Mean-field theory to describe the spin degrees of freedom, the self-consistent Born approximation to characterize individual charge excitations across the Hubbard gap, and the Bethe-Salpeter equation to identify bound states of doublons and holons. The binding energy of the Mott exciton is found to increase with increasing the Néel order parameter, whereas the exciton mass decreases. We observe that these trends rely significantly on the retardation of the effective interaction, and require consideration of multiple effects from changing the magnetic order. Our results are consistent with the key qualitative trends observed in recent experiments on iridates. Moreover, the findings could have direct implications on ultracold atom Mott insulators where the Hubbard model is the exact description of the system and the microscopic degrees of freedom can be directly accessed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Finding simplicity: unsupervised discovery of features, patterns, and order parameters via shift-invariant variational autoencoders *

Abstract Recent advances in scanning tunneling and transmission electron microscopies (STM and STEM) have allowed routine generation of large volumes of imaging data containing information on the structure and functionality of materials. The experimental data sets contain signatures of long-range phenomena such as physical order parameter fields, polarization, and strain gradients in STEM, or standing electronic waves and carrier-mediated exchange interactions in STM, all superimposed onto scanning system distortions and gradual changes of contrast due to drift and/or mis-tilt effects. Correspondingly, while the human eye can readily identify certain patterns in the images such as lattice periodicities, repeating structural elements, or microstructures, their automatic extraction and classification are highly non-trivial and universal pathways to accomplish such analyses are absent. We pose that the most distinctive elements of the patterns observed in STM and (S)TEM images are similarity and (almost-) periodicity, behaviors stemming directly from the parsimony of elementary atomic structures, superimposed on the gradual changes reflective of order parameter distributions. However, the discovery of these elements via global Fourier methods is non-trivial due to variability and lack of ideal discrete translation symmetry. To address this problem, we explore the shift-invariant variational autoencoders (shift-VAEs) that allow disentangling characteristic repeating features in the images, their variations, and shifts that inevitably occur when randomly sampling the image space. Shift-VAEs balance the uncertainty in the position of the object of interest with the uncertainty in shape reconstruction. This approach is illustrated for model 1D data, and further extended to synthetic and experimental STM and STEM 2D data. We further introduce an approach for training shift-VAEs that allows finding the latent variables that comport to known physical behavior. In this specific case, the condition is that the latent variable maps should be smooth on the length scale of the atomic lattice (as expected for physical order parameters), but other conditions can be imposed. The opportunities and limitations of the shift VAE analysis for pattern discovery are elucidated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Multiscale modeling and cinematic visualization of photosynthetic energy conversion processes from electronic to cell scales

We report conversion of sunlight into chemical energy, namely photosynthesis, is the primary energy source of life on Earth. A visualization depicting this process, based on multiscale computational models from electronic to cell scales, is presented in the form of an excerpt from the fulldome show Birth of Planet Earth. This accessible visual narrative shows a lay audience, including children, how the energy of sunlight is captured, converted, and stored through a chain of proteins to power living cells. The visualization is the result of a multi-year collaboration among biophysicists, visualization scientists, and artists, which, in turn, is based on a decade-long experimental-computational collaboration on structural and functional modeling that produced an atomic detail description of a bacterial bioenergetic organelle, the chromatophore. Software advancements necessitated by this project have led to significant performance and feature advances, including hardware-accelerated cinematic ray tracing and instanced visualizations for efficient cell-scale modeling. The energy conversion steps depicted feature an integration of function from electronic to cell levels, spanning nearly 12 orders of magnitude in time scales. This atomic detail description uniquely enables a modern retelling of one of humanity’s earliest stories—the interplay between light and life.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

All-Order Full-Coulomb Quantum Spectral Line-Shape Calculations

Understanding how atoms interact with hot dense matter is essential for astrophysical and laboratory plasmas. Interactions in high-density plasmas broaden spectral lines, providing a rare window into interactions that govern, for example, radiation transport in stars. However, up to now, spectral line-shape theories employed at least one of three common approximations: second-order Taylor treatment of broadening operator, dipole-only interactions between atom and plasma, and classical treatment of perturbing electrons. In this Letter, we remove all three approximations simultaneously for the first time and test the importance for two applications: neutral hydrogen and highly ionized magnesium and oxygen. Finally, we found 15%–50% change in the spectral line widths, which are sufficient to impact applications including white-dwarf mass determination, stellar-opacity research, and laboratory plasma diagnostics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Itinerant Magnetism in Hydride-Synthesized CaCo 12 B 6

A new compound in the underexplored Ca–Co–B phase space has been discovered, validating high-throughput computations from the Open Quantum Materials Database, which predicted thermodynamic stability for CaCo 12 B 6 in the SrNi 12 B 6 structure type. The synthetic effects of different boron precursors and the advantages of using CaH 2 instead of Ca metal were demonstrated by the short synthesis duration and high purity of CaCo 12 B 6 , in contrast with traditional synthesis routes. Powder X-ray diffraction (PXRD) confirmed that CaCo 12 B 6 shares the SrNi 12 B 6 structure ( R $\bar{3}$m (#166), a = 9.469(4) Å, c = 7.468(2) Å, Z = 3) and is water- and air-stable. High-temperature in situ PXRD indicates that CaCo 12 B 6 is stable below 1050 K under vacuum in a silica capillary. CaCo 12 B 6 decomposes between 693 and 773 K during spark-plasma sintering. Density functional theory calculations indicate that CaCo 12 B 6 is metallic with a ferromagnetic ground state. X-ray absorption near-edge spectroscopy and Bader charge analysis indicate that Co atoms in CaCo 12 B 6 lack ionic character. Magnetometry reveals room-temperature paramagnetism with μ eff = 1.7(1)μ B per Co atom and a Weiss constant of +190(10)K. Ferromagnetic ordering occurs below 172(1)K, resulting in a saturation moment of 0.46 μ B per Co atom. Our findings demonstrate that the hydride route is a viable strategy for discovery of new ternary alkaline-earth-transition metal borides analogous to rare-earth-containing counterparts.

diffraction↗
DOE OSTI2021

Multisublattice cluster expansion study of short-range ordering in iron-substituted strontium titanate

Owing to the challenges in obtaining realistic atomic configurations in large chemical phase spaces, it is not straightforward to describe structure–property relations in materials exhibiting configurational disorder. One example is iron-substituted strontium titanate (SrTi 1–x Fe x O 3–d , STF), a promising perovskite-derivative cathode material in solid oxide fuel cells that exhibits full solid solubility 0 ≤ x ≤ 1 and a tendency to exhibit short-range order. Here we demonstrate a multisublattice cluster expansion (CE) framework and apply it to STF across the full composition range. The CE approach is distinct from more traditional CE formulations in that clusters are defined explicitly by the chemical species distributed among multiple sublattices, rather than via cluster functions of occupation variables with decoration. The modified CE approach makes it easy to distinguish meaningful chemical interactions that are harder to extract from conventional CE, since for the latter chemical identity in a cluster is expressed as a product of site occupations. The least absolute shrinkage and selection operator (LASSO) is implemented as a regression analysis tool to select key clusters and avoid overfitting. We demonstrate this formulation on STF, and show that it can accurately predict configurational energies in comparison to conventional CE. From the key clusters, we identify that short-range ordering between substitutional Fe and oxygen vacancies (V O ) results in the formation of Fe–VO strings. In addition, we consider the stability of STF through CE-based Monte Carlo (MC) simulations and confirm the presence of superstructures that were previously observed in transmission electron microscopy. In this work, analysis of atomic configurations from MC samples reveals variations in the oxidation state of Fe atoms, which can be explained by the ordering tendency of Fe and V O . The cluster description and selection formalism described here may be applied to other disordered multisublattice systems for accurate and efficient material modeling.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Structural and chemical disorder enhance point defect diffusion and atomic transport in Ni 3 Al-based γ' phase

The ordered L1 2 γ' Ni 3 Al intermetallic inclusions are routinely used as strengthening constituents in Ni-based superalloys and recent high-entropy alloys. As there are different metallic elements relative to the alloy matrix, the composition of γ' phase inclusions usually suffers some variations. In this work, we studied the elemental partition and its influence on atomic transport properties in a model (NiCo) 3 Al γ' phase with chemical and structural disorder. Combining the Monte Carlo and molecular static techniques, we first determined the elemental distribution in such compositionally-complexed γ' phase. The results suggest that Co tends to segregate, which profoundly influences defect energetics and diffusion properties. By modeling thermally activated diffusion by molecular dynamics, we found that chemical disorder induced by Co can enhance vacancy diffusion while significantly suppress diffusion via interstitial atom migration, as compared to the stoichiometric γ' Ni 3 Al. We further reveal that these observations are closely related to the defect energetics in these systems. Furthermore, our results explain the effects of disorder, long- and short-range orders on the diffusion properties in multicomponent intermetallic alloys, which helps to understand the stability and, thus, γ'-phase-assisted strengthening at different conditions including annealing and irradiation.

36 MATERIALS SCIENCE↗