Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solid solution”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Discovery and Characterization of a Metastable Cubic Interstitial Nickel–Carbon System with an Expanded Lattice

Metastable, i.e., kinetically favored but thermodynamically not stable, interstitial solid solutions of carbon in iron are well-understood. Carbon can occupy the interstitial atoms of the host metal, altering its properties. Alloying of the host metal results in the stabilization of the FeC x phases, widening its application. Pure nickel finds niche applications, mainly focusing on catalysis, while nickel alloys are widely applied, e.g., in gas turbines, reactors, and seawater piping. Nickel carbide (Ni 3 C) is the well-known stable Ni–C system displaying a rhombohedral (R3̅c) crystal structure. Some reports describe an elusive cubic Ni–C system, observed during certain catalytic reactions occurring on nickel and formed by the occupation of the interstitials of the metal with carbon: to date, the stabilization and characterization of this phase have not been accomplished. Hereby, we report on the synthesis of a cubic metastable NiC x phase using chemical vapor deposition of methane on supported nickel nanoparticles. The structure was predicted by DFT/ReaxFF, synthesized and monitored with in situ time-resolved synchrotron XRD, and experimentally confirmed by Rietveld refinement and (S)TEM-EELS under ambient conditions. The results show an Fm3̅m phase with a lattice parameter of a = 3.749 ± 0.037 Å at room temperature, with the highest ever reported atomic percentage of carbon occupying the octahedral interstices of 23.1%, resulting in a NiC 0.3 phase. The degree of occupation of the interstitial voids by carbon can be controlled, enabling the tuning of the host metal’s d-spacing and composition, highlighting the applicability of this synthesis route for catalytic nanoparticle preparation.

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

NbZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys niobium-zirconium (Nb-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Nb and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

TaZr_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution binary alloys tantalum-zirconium (Ta-Zr). The electronic structures of alloys have been calculated using Vienna Ab initio Simulation Package (VASP). Within this package the DFT approach is used to reduce many-body Schrodinger equation to set of single particle Kohn-Sham (KS) equations. The generalized electronic exchange-correlation functional is described by generalized gradient approximation with the Perdew-Burke-Ernzerhof parametrization. The electron-ion interactions is described by pseudopotentials developed within the plane-wave basis projector augmented-wave (PAW) approach \cite{PAW}. These pseudopotentials are available at the VASP portal (http://cms.mpi.univie.ac.at/vasp/). Our calculations have been run with the pseudopotentials treating s and p semi-core states as valence in case for the elements Ta and Zr. The electronic densities and potentials are expanded over plane-waves with energy cutoff of 350 eV. 2x2x2 k-mesh and normal precision were used. The alloys were modeled by supercell containing 128 randomly distributed atoms. At initial step the atoms occupy perfect bcc lattice cites. This initial structure was optimized until energy changes less than 1e-6 eV, while forces acting on atoms don't exceed 1e-2 eV/angstrom. The electron-ion interaction is described by PAW pseudopotentials. The calculations have been collected by sampling chemical compositions across the entire compositional range. The chemical compositions have been sampled by progressively changing the number of atoms per constituent by 4. For each chemical composition of binaries and ternaries, the first-principle calculations have been run for 100 randomized arrangements of the constituents on the BCC lattice sites. We collected data for a total of 3,100 randomized atomic structures over 31 chemical compositions. The calculations have been collected on NERSC-Perlmutter and OLCF-Summit using the VASP 6.3.2. The VASP calculations for every atomic structure have been performed in 2 main steps: 1. Starting from an ideal body-centered-cubic (BCC) structure, geometry optimization with low precision has been executed to perform a preliminary optimization of the atomic structure. The output for this calculations is available in the files 0.CONTCAR, 0.OUTCAR, rlx1.out. 2. Using the atomic structure resulting from the preliminary geometry optimization, a second geometry optimization has been performed using normal precision. The output for this calculations is available in the files CONTCAR, OUTCAR, rlx2.out, vaspout.h5, and vasprun.xml. Cases 1-10 have been run without generating the file 'vaspout.h5'. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for binary alloys is AXBY, where A and B refer to the constituents, whereas X and Y are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: FILES contained in each subdirectory with name "case-N" where N ranges between 11 and 80, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. 0.POSCAR: input file that defines the atomic structure of a system 4. 0.CONTCAR: output file that provides the atomic positions and cell parameters after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 5. 0.OUTCAR: output file that contains detailed information about the progress of a calculation after the first geometry optimization has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.out: file with diagnostic information about the execution of the first geometry optimization with precision variable set to PREC=Low in the INCAR file 7. POSCAR: input file that defines the atomic structure of a system after the first geometry optimization has been run at low precision. This represents the input for the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. CONTCAR: output file that provides the atomic positions and cell parameters after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 9. OUTCAR: output file that contains detailed information about the progress of a calculation after the second geometry optimization has been run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.out: file with diagnostic information about the execution of the second geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system FILES contained in each subdirectory with name "case-N" where N ranges between 1 and 10 and between 81 and 100, extremes included: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. KPOINTS: input file that specifies the Bloch vectors (k points) used to sample the Brillouin zone 3. {ID}.POSCAR: input file that defines the atomic structure of a system at the beginning of ID execution of geometry optimization with PREC=LOW 4. {ID}.CONTCAR: output file that provides the atomic positions and cell parameters at the end of ID execution of geometry optimization with PREC=LOW in the INCAR file 5. {ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of geometry optimization that has been run with the precision variable set to PREC=Low in the INCAR file 6. rlx1.{ID}.out: file with diagnostic information about the execution of the ID execution of the geometry optimization with precision variable set to PREC=Low in the INCAR file 7. N{ID}.POSCAR: input file that defines the atomic structure of a system after the geometry optimization run at low precision. This represents the input for the ID execution of the second geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 8. N{ID}.CONTCAR: output file that provides the atomic positions and cell parameters after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 9. N{ID}.OUTCAR: output file that contains detailed information about the progress of a calculation after the ID execution of the geometry optimization run with the precision variable set to PREC=Normal in the INCAR file 10. rlx2.{ID}.out: file with diagnostic information about the ID execution of geometry optimization with precision variable set to PREC=Normal in the INCAR file 11. vaspout.h5: hierarchical HDF5 file containing the inputs and outputs of a VASP calculation. To analyze the data in this file we recommend using py4vasp. This file is only produced if the VASP version used is compiled with HDF5 support 12. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 13. CHGCAR: contains the charge density data of a system. This data is crucial for analyzing electronic structures, calculating electrostatic potential, and studying the distribution of charge in a crystal or molecular system This research is sponsored by the Artificial Intelligence Initiative as part of the Laboratory Directed Research and Development (LDRD) Program of Oak Ridge National Laboratory, managed by UT-Battelle, LLC, for the US Department of Energy under contract DE-AC05-00OR22725. This work used resources of the Oak Ridge Leadership Computing Facility, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC05-00OR22725, under Directorate Discretionary awards MAT025 (Materials Science) and LRN026 (Machine Learning), and INCITE award MAT201. This work also used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231, under award ERCAP0025216. REFERENCES (1) Kresse, G. & Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. & Hafner, J. Ab initio molecular-dynamics simulation of the liquid-metal–amorphous-semiconductor transition in germanium. Phys. Rev. B 49, 14251 (1994) (3) Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. materials science 6, 15–50 (1996) (4) Kresse, G. & Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. review B 54, 11169 (1996) (5) Kresse, G. & Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

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 OSTI2026

Conversion of real-world aluminum scrap streams into high-performance Al–Mg–Si–Cu automotive alloys using shear assisted processing and extrusion

The conversion of post-consumer aluminum (Al) scrap into usable Al alloys without adding primary Al is challenging because of excess impurities. In this work, >99% post-consumer Twitch, used beverage cans (UBCs), and remelt scrap ingots (RSIs) were used as feedstock materials. As-cast and solution heat-treated feedstock billets were extruded using Shear Assisted Processing and Extrusion (ShAPE) at ~510°C, followed by press quenching. To explore the development of the microstructure, texture, and underlying mechanisms and how they contribute to the overall strengthening in as-extruded and artificially aged samples, scanning electron microscopy (SEM) and electron backscatter diffraction (EBSD) were used to collect microstructure and texture data. The enhanced strength and ductility were corroborated with the microstructural features and crystallographic texture. Simple shear $\textrm{A}/\bar{\textrm{A}}$, $\textrm{A}_1^*/\textrm{A} _2^*$ texture components along with weak $\textrm{C}$ and $\textrm{B}/\bar{\textrm{B}}$ texture components were formed during extrusion; the texture was strengthened after heat treatment. The refined second-phase particles helped to retain the deformed microstructure and texture. The contributions of dislocation and precipitate strengthening were maximized when billets were solution-heat-treated prior to extrusion. This is attributed to the formation of effective supersaturated solid solutions during the ShAPE process, which precipitate out during the peak age treatment. Overall, the highest yield strength of 305 MPa, ultimate tensile strength of 350 MPa, and elongation of 12% were achieved in artificially aged samples, which are comparable to those of Al 6082-T6.

Al alloy↗
DOE OSTI2025

Modeling phase transformations in Mn-rich disordered rocksalt cathodes with machine-learning interatomic potentials

Mn-rich disordered rocksalt (DRX) cathode materials exhibit a phase transformation from a disordered to a partially disordered spinel-like structure (δ-phase) during electrochemical cycling. Here, in this computational study, we use charge-informed molecular dynamics with a fine-tuned CHGNet foundation potential to investigate the phase transformation in LixMn 0.8 Ti 0.1 O 1.9 F 0.1 . Our results indicate that transition metal migration occurs and reorders to form the spinel-like ordering in an FCC anion framework. The transformed structure contains a higher concentration of nontransition metal (0-TM) face-sharing channels, which are known to improve Li transport kinetics. Analysis of the Mn valence distribution suggests that the appearance of tetrahedral Mn 2+ is a consequence of spinel-like ordering, rather than the trigger for cation migration as previously suggested. Calculated equilibrium intercalation voltage profiles demonstrate that the δ-phase, unlike the ordered spinel, exhibits solid-solution signatures at low voltage. A higher Li capacity is obtained than in the DRX phase. This study provides atomic insights into solid-state phase transformation and its relation to experimental electrochemistry, highlighting the potential of machine-learning interatomic potentials for understanding complex oxide materials.

Zhong, Peichen [University of California, Berkeley↗
DOE OSTI2025

Synthesis of a series of rare-earth-based multi-anion chalcogenide iodides RE 3 Si 2 Se x S 8− x I (RE = La, Ce, Pr, and Nd) using the flux-assisted boron–chalcogen mixture method

Single crystals and polycrystalline powders of rare earth mixed chalcogenide iodides La 3 Si 2 Se 1.21 S 6.79 I, Ce 3 Si 2 Se 1.39 S 6.61 I, Pr 3 Si 2 Se 1.22 S 6.78 I, and Nd 3 Si 2 Se 1.18 S 6.82 I were prepared using the reactive flux-assisted boron–chalcogen mixture (BCM) method at 850 °C. All compounds crystallized in the monoclinic crystal system, space group C2/c (space group number 15). The series adopts the La 3 Si 2 O 8 Cl structure type, containing isolated SiQ 4 tetrahedra connected by REQ 8 (RE = La, Ce, Pr and Nd) polyhedra; this arrangement creates tunnels that are filled by I atoms. The partial substitution of S by Se was carried out to modulate the optical properties. Phase pure samples and uniform solid solutions were obtained for all compositions as determined using powder X-ray diffraction patterns. Polycrystalline powders were used for physical property measurements, including magnetic susceptibility and UV-Vis diffuse reflectance. The solid-state UV-Vis data for the polycrystalline La 3 Si 2 Se 1.21 S 6.79 I, Ce 3 Si 2 Se 1.39 S 6.61 I, and Pr 3 Si 2 Se 1.22 S 6.78 I samples revealed band gaps of E g = 2.5(1), 2.2(1), and 2.3(1) eV, typical of semiconductors. Magnetic measurements indicated that Ce 3 Si 2 Se 1.39 S 6.61 I and Nd 3 Si 2 Se 1.18 S 6.82 I exhibit paramagnetic behavior with slightly negative Weiss constants θ = −25 and −38. The photoluminescence spectrum of Ce 3 Si 2 Se 1.39 S 6.61 I exhibits a broad emission band around ∼493 nm.

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

Carbonation and Sulfidation of Mg- and Ni-Containing Solutions: Implications for Carbon Mineralization and Critical Element Recovery

Carbonation of alkaline earth metals ( e.g. , magnesium (Mg)) and sulfidation of nickel (Ni) are promising methods to achieve concurrent carbon mineralization and selective Ni recovery. However, the coexistence of alkaline earth metals and Ni from silicate ores or mining wastewater complicates the carbonation and sulfidation owing to cation coprecipitation. To better understand simultaneous metal carbonation and Ni-sulfide formation, we used Mg- and Ni-containing solutions and systematically investigated the Mg and Ni coprecipitates’ phase transformation during sequential/concurrent carbonation and sulfidation. During a single carbonation process, hydromagnesite dehydrated and formed magnesite over time. Nickel bicarbonate formed and became a Mg–Ni carbonate solid solution because of their similar ionic radii. During a single sulfidation process, the pH did not affect Ni-sulfide formation, but it controlled Mg behavior. Specifically, at pH 9.6, brucite formed, while at pH 7.8, Mg 2+ remained in the solution. For the sequential carbonation–sulfidation process, Ni-carbonate formed during carbonation converted to Ni-sulfide because of the low Ni-sulfide K sp . For the sulfidation–carbonation process, Ni-sulfide remained the same even after carbonation and Mg-carbonate precipitates. For the concurrent carbonation and sulfidation process, Mg-carbonate and Ni-sulfide formed simultaneously. Finally, this study develops a scientific foundation of carbonation and sulfidation processes, benefiting coupled CO 2 storage and sulfide-enabled resource recovery.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Systematic discovery of new nano-scale metastable intermetallic eutectic phases in laser rapid solidified Aluminum-Germanium alloy

Laser surface remelting of as-cast Al-Ge eutectic alloy is shown to produce ultrafine lamellar eutectic morphology with interlamellar spacing refined up to ∼60 nm and composed of FCC Al solid solution and unusual Al x Ge y intermetallic phases that do not form during near-equilibrium solidification. The microstructures are characterized and analyzed using a combination of selected area electron diffraction, high-resolution scanning transmission electron microscopy, energy dispersive X-ray spectroscopy to obtain high-resolution elemental maps, and atomistic modeling using density functional theory followed by atomic-scale image simulation. Depending on the local solidification conditions, the crystallography of the Al x Ge y intermetallic phases in the eutectic microstructure is either monoclinic (C 2/c) or monoclinic (P 2 1 ), with high densities of defects in both cases. This is in sharp contrast to the as-cast alloys that showed nominally pure Al and Ge phases with significant solute partitioning and equilibrium FCC and diamond cubic crystal structures, respectively. Corresponding kinetic phase diagrams are proposed to interpret the evolution of nano-lamellar eutectic morphologies with equilibrium Al and metastable Al x Ge y phases, and to explain increased solid solubility in the Al phases manifested by precipitation of ultrafine clusters of Ge. Furthermore, the reasons for the formation of these metastable eutectics under laser rapid solidification are discussed from the perspective of the competitive growth criterion.

Al-Ge eutectic↗
DOE OSTI2025

Strain-associated nanoscale fluctuating lithium transport within single-crystalline LiNi 1/3 Mn 1/3 Co 1/3 O 2 cathode particles

Solid-state lithium diffusion dynamics are critical for the rate capability and longevity of Li-ion batteries. Conventionally, nanoscale lithium diffusion within individual battery particles has been simplified as being primarily driven by concentration gradients, despite the associated processes inducing local lattice expansion, contraction, and strain fields. Using operando scanning transmission soft X-ray microscopy with high spatial resolution and chemical sensitivity to track nanoscale intraparticle lithium transport, and post-cycling Bragg coherent diffraction X-ray imaging to directly reveal three-dimensional intraparticle strain fields, we uncover strain-associated lithium transport dynamics within single-crystalline LiNi 1/3 Mn 1/3 Co 1/3 O 2 (scNMC) particles during cycling. Contrary to the expected thermodynamic solid-solution behavior of scNMC, our observations reveal near-uniform but fluctuating regions of lithium-dense and lithium-dilute areas during cycling. These fluctuations suggest that nanoscale lithium diffusion can proceed counter to concentration gradients. Additionally, we demonstrate that an increased presence of lithium-dilute regions near the surface enhances lithium surface insertion kinetics, emphasizing the importance of controlling surface lithium distribution to improve rate performance. Our study provides insights into nanoscale solid-state ion transport, with potential applications in batteries, solid-state fuel cells, and memristors.

Lee, Danwon [Seoul National Univ. (Korea, Republic↗
DOE OSTI2026

Sodium-Ion Battery Cathode with Dominating Copper and Oxygen Redox Chemistry

Sodium-ion batteries offer low-cost energy storage solutions for the grid and electric vehicles, leveraging the established "rocking-chair" Li-ion design and the natural abundance of sodium. However, SIBs face challenges such as relatively lower voltage and capacity than lithium-ion batteries, as well as dependence on nickel resources. Here, in this work, a new nickel-free cathode material, Na 0.75 Li 0.08 Cu 0.25 Mn 0.66 O 2 , was designed and synthesized. This material has a capacity of ~125 mAh/g and an average discharge voltage of 3.5 V. Notably, more than one-third of the capacity arises from lithium substitution of Cu (~8 mol.%) and high voltage activation to 4.6 V. Multimodal synchrotron x-ray characterization combining spectroscopy, microscopy, and scattering reveal the capacity is primarily from the redox of copper and oxygen, with a minor contribution from the manganese redox. Lithium substitution alters the phase transition mechanism from a two-phase transition in P3-Na 2/3 Cu 1/3 Mn 2/3 O 2 to a solid-solution in Na 0.75 Li 0.08 Cu 0.25 Mn 0.66 O 2 , enhancing the reversibility of this material.

25 ENERGY STORAGE↗
DOE OSTI2025

An energetic link between order and strength in metals: A nanocrystalline strength limit in high-entropy alloys and intermetallic compounds

The metallurgy and materials communities have long understood and exploited fundamental links between chemical and structural ordering in metallic solids to tailor their mechanical properties. We extend these ideas to include prediction of the nanocrystalline strength limit in high-entropy alloys and intermetallic compounds, where a breakdown occurs in the classical Hall-Petch strengthening behavior. The highest reported strength achievable through alloying has rapidly climbed and given rise to new classifications of materials with extraordinary properties, with a notable case being nanocrystalline metals. High-entropy alloys (chemically disordered, concentrated solid solutions) and intermetallic compounds are two boundary cases of how tailored order can be used to manipulate mechanical behavior. Here, we show that the complex electronic-structure mechanisms governing the peak strength of alloys and pure metals can be reduced to a few physically meaningful parameters based on their atomic arrangements and used – with no fitting parameters – to predict the maximum strength of these materials. This includes a generalized energy-based accounting for the degree of structural and chemical ordering that allows for rapid and reasonably accurate prediction of peak strength (validated in the nanocrystalline limit) as a function of temperature. Predictions of maximum strength based on the activation energy (with all materials properties derived from DFT calculations or experiments) for a stress-driven transition to an amorphous state is shown to accurately describe the breakdown in Hall-Petch behavior at the smallest crystallite sizes for pure metals, intermetallic compounds, high-entropy alloys, and metallic glasses. Further, this activation energy is also shown to be directly proportional to interstitial electronic charge density, which is a good predictor of ductility, stiffness (moduli), and phase stability in high-entropy alloys and solid metals generally. The proposed framework suggests the possibility of coupling ordering and intrinsic strength to mechanisms like dislocation nucleation, hydrogen embrittlement, and transport properties, such as through correlations between the activation energies for amorphization with stacking-fault and grain boundary energies. It additionally opens the prospect for greatly accelerated structural materials design and development to address materials challenges limiting more sustainable and efficient use of energy.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Nanoionics Drastically Accelerating Mass Transfer at Elevated Temperatures over 750 °C

Nanoionics were previously considered thermally unstable and infeasible for devices operating above 500 °C. Here, we elucidate the design principle for establishing stable nanoionics from various oxides. We utilized reversible solid oxide cells (SOCs) as the test bed and implemented nanoionics using atomic layer deposition (ALD). We demonstrate a straightforward, interface-controlled, practical approach to render a conformal, ∼15 nm thick ALD film, which initially thermodynamically favors the formation of a solid solution with the substrate into surface nanoionics with single or double layers of nanograins with random crystal orientations. The nanoionics exhibited conductivity estimated to be 7 orders of magnitude higher than that of their bulk-scale counterpart. They demonstrated conformability with uniform grain sizes of ∼15 nm, even after electrochemical operation for ∼500 h at 750 °C and 1000 h at 850 °C. The thermal stability and conductivity of such nanoionics represent a conceptual and technological framework in nanoionics.

atomic layer deposition↗
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 OSTI2025

Experimental evidence of disordered crystalline premixing in sputter-deposited Ni(V)/Al multilayers

The sputter deposition of alternating layers of Ni(V) and Al forms a reactive multilayer known to undergo self-propagating formation reactions when ignited. The sequential deposition process leads to nanometer-scale premixing of reactants at each included interface, which ultimately affects multilayer exothermicity. This work performs the direct measurement of a disordered face-centered cubic (FCC) solid solution premixed phase at the interfaces of Ni(V)/Al multilayers via scanning transmission electron microscopy. The crystallinity of the observed phase differs from previously reported a priori predictions of an amorphous interlayer. The disordered FCC phase retains its symmetry after annealing for 16 h at 135 ± 5 °C, but the lattice parameter shifts consistently with an Al-rich composition. The existence of a crystalline premix in Ni(V)/Al is attributed to the electronic contribution to the entropy of crystallization. The importance of electronic entropy to the phase formation of energetic materials motivates its inclusion when constructing digital twins for atomistic kinetics and ignition sensitivity.

Crystal structure↗
DOE OSTI2024

Tuning Perovskites’ Hydration-Induced Chemical Expansion with Octahedral Tilt Angles

Hydration-induced strains in proton-conducting oxides compromise chemo-mechanical stability when these materials are applied in protonic ceramic electrochemical cells. Here, to develop design principles for zero-strain materials, we systematically studied the hydration coefficients of chemical expansion (CCEs) in perovskite (Sr, Ba)(Ce, Zr, Y)O 3–x solid solutions with in situ dilatometry and thermogravimetric analysis in the range of 430–630 °C. By including and decoupling a wide range of tolerance factors and lattice parameters, we were able to identify a minimum in hydration CCEs (0–0.02) at intermediate tolerance factor values (t ≈ 0.95). Conversely, despite expectations of lower CCEs in larger unit cells, no general trend in CCE versus lattice parameter was found, and opposite trends could be seen for Sr(Ce, Zr, Y)O 3–x versus Ba(Ce, Zr, Y)O 3–x separately. In situ neutron diffraction (ND) enabled atomistic insight. Upon decreasing t, chemical strain anisotropy increased, but this trend did not match the U-shaped dependence of macroscopic CCEs on t. Instead, perovskites with intermediate t, hosting intermediate octahedral tilt angles in the nominally dry state, underwent the largest change in the B–O–B angles during hydration. Accommodating hydration through decreasing B–O–B angles is beneficial because it does not result in large lattice parameter changes. We propose an intermediate tolerance factor as a simple structural descriptor to enable near-zero hydration strains in proton-conducting perovskites.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Electro-chemo-mechanically Driven Ni Exsolution from (Pr,Ce,Ni)O 2−δ : Controlled Nucleation Density and Enhanced Electrode Kinetics

In situ exsolution of metal nanoparticles is a promising strategy to prepare electrocatalysts with enhanced activity and resistance to agglomeration for efficient chemical transformations and energy conversion. Achieving a high nucleation density of nanoparticles under mild conditions and understanding how to tailor the process is important for performance of these electrodes in electrochemical cells. In this work, we demonstrate facile exsolution of Ni nanoparticles using fluorite-structured (Pr,Ce)O 2−δ as the support oxide, driven by electrochemical potential and aided by the metastability of Ni in the solid solution (elastic driving force). We prepare single-phase oriented thin films of (Pr,Ce,Ni)O 2−δ (NPCO) on (Zr,Y)O 2−δ (YSZ) substrates by pulsed laser deposition. With the aid of a high-throughput electrochemical cell that provides a lateral gradient in Nernst voltage, we apply in situ near-ambient pressure synchrotron X-ray photoelectron spectroscopy and ex situ atomic force microscopy to investigate the impact of electrochemical potential on Ni nucleation density. We find that metallic Ni can be successfully exsolved at 550 °C upon cathodic biasing in 20 mTorr O 2 , and its nucleation density increases with increasing electrochemical driving force/decreasing oxygen chemical potential. We further evaluate the electrochemical performance under highly reducing (fuel electrode) conditions by electrochemical impedance spectroscopy. With the exsolved Ni nanoparticles, the surface exchange coefficient of the NPCO is found to be ∼4× higher than for PCO without exsolution. This work confirms mixed conducting fluorites as beneficial host lattices for facile transition-metal exsolution and suggests the possibility for constructing an all ceria-based electrochemical cell with PCO serving as both the cathode and the anode.

36 MATERIALS SCIENCE↗