Search NASA⌕ Search

SEARCH · Search NASA

Results for “crystal structure optimization”

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

Synthesizing intermetallic thermoelectric materials using metal flux growth method

With zero energy poverty by 2050 goals, exploring innovative energy technologies is essential. Thermoelectric materials (TEMs) have the potential to harness waste heat and convert it to electricity (e.g., for spacecraft power, sensors, and cryogenic cooling). Two limitations of thermoelectric compounds are their low efficiency and limited pool of materials. Optimizing the efficiency of TEMs is challenging due to the interdependent relationship between electrical and thermal conductivity. A deeper understanding of electronic properties and lattice thermal conductivity is necessary to “understand, predict and ultimately control” the thermoelectric properties of TEMs. Understanding the influence of light-atom-doping and intrinsic defects on the electronic properties of the f-element containing TEMs are important questions that need to be answered to design and synthesize TEMs with optimal parameters. The purpose of this work is to understand how doping f-element based TEMs with light elements, such as carbon and boron, would impact the crystal structure, bond lengths, band gaps, and thermoelectric characteristics. High-quality f¬-element containing intermetallic single crystals are grown in metal flux media. In molten metal flux growth method, one or more low-melting metals such as aluminum, gallium, tin, and bismuth are used in excess as the reaction medium. As these are solution-state reactions, we can use lower temperatures than solid-state reactions and isolate the kinetic products instead of the thermodynamically stable compound. Then these are characterized using X-ray diffraction, scanning electron microscopy, energy dispersive spectroscopy, magnetization, electrical and thermal conductivities, and thermoelectric power.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2025

Local orthorhombic phase in zirconium oxide nanocrystals: insights from X-ray pair distribution function analysis

Zirconium dioxide (ZrO 2 ) and hafnium dioxide (HfO 2 ) have emerged as promising alternatives to conventional ferroelectric materials. Understanding the crystal phases of these oxides under different conditions is crucial for optimizing their properties. There are several theories for the (anti)ferroelectric properties; however, comprehensive analysis, particularly at the local structure level, is lacking. In this study, we investigate the local structure of ZrO 2 nanocrystals using X-ray pair distribution function (PDF) analysis, revealing an unexpected local orthorhombic distortion irrespective of crystallite size. This finding suggests the potential existence of an intermediate orthorhombic phase during the microscopic switching pathway observed in previous studies. Additionally, we explore the influence of crystallite size and surface effects on the PDF. These results contribute to a deeper understanding of the structural dynamics in ZrO 2 and offer insights for the design of next-generation ferroelectric materials.

Pokratath, Rohan (ORCID:0000000268383939)↗
DOE OSTI2025

Bio-inspired multiscale design for perovskite solar cells

Metal halide perovskite semiconductors have attractive light-harvesting and charge-carrier transport properties for photovoltaics. Perovskite solar cells (PSCs) and modules have demonstrated their commercial promise with high power conversion efficiencies, but still face stability challenges. In this Review, we explore how biomaterials offer design inspiration for the development of durable and efficient PSCs at three different scales. At the molecular level, bio-inspired molecular interactions are harnessed towards crystallization control and degradation prevention, which offers an enhancement in long-term maximum-power-point tracking stability. At the microstructural level, self-healing and strength-enhancing strategies, utilizing dynamic bonds and interfacial reinforcement, can help PSCs to recover from physical damage and maintain high performance. At the device level, macroscopic functionalities, such as moth-eye-inspired structures tailored to different layers, can collectively enable antireflection, radiative cooling and self-cleaning to optimize light management, heat dissipation and encapsulation in PSCs. Bio-inspired PSC research can combine improved efficiency and lifetime, with abundant, biocompatible alternatives to conventional stabilizers. Future efforts should focus on screening bioinspired molecules to optimize film crystallization and stability, developing self-healing mechanisms triggered by operational stress, designing cost-efficient biomicrostructures, and integrating multifunctional encapsulation to enhance the efficiency and lifespan of PSCs.

Duan, Tianwei↗
DOE OSTI2026

Laser-Induced Enhancement of Acoustic Mode Quality Factor Revealed by Correlated Single-Particle Optical and Electron Microscopy

Acoustic modes in plasmonic nanostructures provide fundamental insights into their optomechanical behavior at the nanoscale, enabling emerging applications in plasmon-enhanced optomechanics, ultrasensitive sensing, and nanoscale energy transduction. Here we explore the modulation of acoustic phonon dynamics in lithographically fabricated gold nanodisks via laser-induced photothermal annealing. Using a correlated approach that utilizes both single-particle transient extinction spectroscopy and advanced electron microscopy, we directly link nanoscale structural transformations to changes in mechanical properties as probed through the coherence of the excited acoustic modes. Specifically, ultrafast pump–probe microscopy reveals an enhancement in the acoustic mode quality factor of annealed gold nanodisks, indicative of reduced damping and improved vibrational coherence. Structural characterization via scanning electron microscopy and electron backscatter diffraction confirms that photoinduced annealing results in smoother surface morphology and overall enhanced crystallinity. The improved crystalline order reduces defect and crystal boundary scattering, which we suggest as the reason underlying the lower quality factor before annealing. Furthermore, these findings demonstrate that targeted structural engineering at the nanoscale offers a powerful strategy for optimizing the optomechanical performance of plasmonic nanostructures, with broad implications for the design of next-generation nanophotonic and optomechanical systems.

Annealing (metallurgy)↗
DOE OSTI2024

Virtual Substrates for Wide Bandgap AlyX1-yN Growth

Lattice-matched substrates are critical for growth of high quality, compositionally-targeted ternary AlyX1-yN compounds with properties suitable for a wide variety of next-generation opto- and power electronic applications. [DOI:10.1149/2.0111702jss] We have identified the (111) plane of transition metal carbides and nitrides as lattice matched "virtual" substrate layers that have additional benefits of electrical conductivity and appropriate coefficients of thermal expansion for nitride layer at both growth and operating temperatures. [arXiv:2208.11769 2022] In this work, (111)-oriented TaC and ZrN are grown by RF sputtering and optimized as substrate layers for AlyGa1-yN and AlyGd1-yN, respectively. TaC is demonstrated as a full proof of concept. (111)-stabilized TaC layers with a rock-salt crystal structure are annealed to improve surface crystal quality and surface morphology and then used as a template for growth of an Al0.7Ga0.3N layer by molecular beam epitaxy. X-ray diffraction (XRD) demonstrates epitaxial registry of the grown layer to the substrate. High resolution transmission electron microscopy (TEM) is used to investigate interface behavior, showing regions of abrupt interface transitions and nitrogen polar termination of the Al0.7Ga0.3N. We also grew (111)-oriented ZrN on Ti-Zr-N graded buffer layers as a lattice-matched virtual substrate for AlyGd1-yN. Structure, morphology, and strain are analyzed using XRD and atomic force microscopy. Structure and quality of initial Al0.8Gd0.2N thin films are investigated for both in-situ and ex-situ growths.

AlyX1-yN↗
DOE OSTI2024

Crystal structure of perfluorononanoic acid, C 9 HF 17 O 2

The crystal structure of perfluorononanoic acid (PFNA) was solved via parallel tempering using synchrotron powder diffraction data obtained from the Brockhouse X-ray Diffraction and Scattering (BXDS) Wiggler Lower Energy (WLE) beamline at the Canadian Light Source. PFNA crystallizes in monoclinic space group P2 1 /c (#14) with lattice parameters a = 26.172(1) Å, b = 5.6345(2) Å, c = 10.9501(4) Å, and β = 98.752(2)°. The crystal structure is composed of dimers, with pairs of PFNA molecules connected by hydrogen bonds via the carboxylic acid functional groups. The Rietveld-refined structure was compared to a density functional theory-optimized structure, and the root-mean-square Cartesian difference was larger than normally observed for correct powder structures. The powder data likely exhibited evidence of disorder which was not successfully modeled.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2026

Open-Air Manufacturing of Efficient Large-Area Perovskite Solar Cells to Meet Stability and Cost Targets

This SETO-funded project (DE-EE0009516) developed a fully open-air manufacturing pathway for perovskite solar cells and modules, eliminating vacuum-based processing steps while improving device stability. Open-air ultrasonic spray deposition, solution combustion synthesis, rapid thermal processing, and atmospheric plasma curing were developed and optimized for all functional device layers — perovskite absorber, hole transport layers (HTLs), electron transport layers (ETLs), transparent conducting oxides (TCOs), and rear electrodes — without reliance on vacuum. Sprayed organic ETLs preserved perovskite crystal structure and enabled power conversion efficiencies exceeding 20%, and open-air processed TCOs and rear electrodes achieved optoelectronic performance comparable to sputtered alternatives at substantially lower manufacturing cost. Device stability was assessed through accelerated aging protocols aligned with IEC and ISOS standards, combined with high-throughput optical diagnostics that identified key degradation mechanisms and antagonistic interactions among heat, humidity, and illumination stressors. A comprehensive technoeconomic analysis quantified the cost and throughput advantages of open-air manufacturing over conventional vacuum-based processes, identifying vacuum deposition as the dominant cost driver and demonstrating a clear pathway toward the $0.02/kWh target. The project produced five peer-reviewed journal publications, eleven conference presentations, two patents/applications, and five Ph.D. dissertations.

14 SOLAR ENERGY↗
DOE OSTI2026

Multinuclear Solid-State NMR and NMR Crystallography of Solid Forms of Creatine and Creatinine

Creatine is a performance-enhancing supplement with two widely available commercial solid forms, namely, creatine monohydrate (creatine·H 2 O) and creatine HCl, the latter of which does not have a reported crystal structure. Moreover, commercial formulations of creatine may contain creatinine, an undesired impurity phase resulting from the self-cyclization of creatine during manufacturing. Therefore, reliable methods for characterizing the different solid forms of creatine and detecting the presence of creatinine are essential. Herein, we address these challenges using 13 C, 15 N, and 35 Cl solid-state NMR (SSNMR) spectroscopy to obtain distinct spectral fingerprints for creatine·H 2 O and creatine HCl, along with creatinine and creatinine HCl. The acquisition of these SSNMR spectra offers a robust approach for both the rapid characterization of each solid form and the detection of the impurity phases. Additionally, quadrupolar NMR crystallography-guided crystal structure prediction (QNMRX-CSP) was applied for the de novo crystal structure determination of creatine HCl, which was validated by the subsequently determined single-crystal X-ray diffraction (SCXRD) structure. Finally, to investigate the relationship between NMR parameters and structural features, 13 C and 15 N chemical shifts and 35 Cl electric field gradient (EFG) tensors were computed from geometry-optimized structures of the four solid forms by using dispersion-corrected DFT-D2* methods. Finally, this integrative approach offers a powerful framework for advancing the structural understanding and quality control of creatine-based supplements and next-generation formulations, as well as a wide range of other solid pharmaceuticals and nutraceuticals.

NMR↗
DOE OSTI2024

Force Field X: A computational microscope to study genetic variation and organic crystals using theory and experiment

Force Field X (FFX) is an open-source software package for atomic resolution modeling of genetic variants and organic crystals that leverages advanced potential energy functions and experimental data. FFX currently consists of nine modular packages with novel algorithms that include global optimization via a many-body expansion, acid–base chemistry using polarizable constant-pH molecular dynamics, estimation of free energy differences, generalized Kirkwood implicit solvent models, and many more. Applications of FFX focus on the use and development of a crystal structure prediction pipeline, biomolecular structure refinement against experimental datasets, and estimation of the thermodynamic effects of genetic variants on both proteins and nucleic acids. The use of Parallel Java and OpenMM combines to offer shared memory, message passing, and graphics processing unit parallelization for high performance simulations. Overall, the FFX platform serves as a computational microscope to study systems ranging from organic crystals to solvated biomolecular systems.

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

Hierarchal structures tuned electrocaloric and electromechanical performance in PVDF-based tetrapolymers

Ferroelectrics with multifunctionalities are gaining increased interest in self-actuated electrocaloric effect (ECE) refrigerators. However, achieving high ECE and electromechanical (EM) coupling concomitantly for maximum heat transfer remains challenging. Here we present the structure-property relationship for poly(vinylidene fluoride-co-trifluoroethylene-co-chlorofluoroethylene-co-double bond), P(VDF-TrFE-CFE-DB), tetrapolymer, which exhibited a high ECE entropy change of 66.5 J Kg −1 K −1 and EM strain of −6.1%. We show that thermal treatment can be a key factor influencing multifunctional properties. High-temperature annealing incorporates DB and CFE units into crystalline grains to form extended-chain crystals, enabling CFE units to induce relaxor behavior and DB units to induce large structural changes at low electric fields. This synergy leads to an enhancement in both ECE and EM performances. Furthermore, at an optimized temperature of 50 °C, the annealed films exhibit giant cross-energy coupling, achieving ECE and EM performances of 100.8 J Kg −1 K −1 and −7.6%. This study provides insights into developing new ferroelectric polymers with electroactive multifunctionalities.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effects of Cooling Rates on Self-Reinforced Polyethylene Composites via Multiscale Experimental Characterizations

This study investigates cooling rate effects on self-reinforced polyethylene composites (PE-SRCs), focusing on macroscopic thermal and mechanical properties, and micro- to sub-micron morphologies. Unlike conventional fiber-reinforced composites (e.g., carbon and glass fiber based), cooling rate effects on SRCs remain less studied. PE-SRCs were fabricated from ultrahigh molecular weight polyethylene (UHMWPE) fabric and high-density polyethylene (HDPE) matrix with cooling rates from 0.76°C/min to 425°C/min. Tensile properties, including modulus (~9.5 GPa), strength (~360 MPa), and ultimate elongation (~6%), showed no significant cooling rate dependence. However, thermal analysis revealed melting temperature of HDPE matrix decreased from ~130°C to ~126.5°C with total crystallinity reducing from 84% to 75%. Wide-angle X-ray scattering confirmed decreasing crystallinity but at lower values (67%–59%), as thermal analysis might overestimate crystallinity. Small-angle X-ray scattering showed long period decreased from 24.3 to 17.4 nm, attributed to thinner lamellae and closer packing, suggesting limited crystal growth and faster nucleation. While cooling rates alter melting point and crystalline structures of HDPE matrix, these changes do not translate to substantial impact on macroscale mechanical properties of PE-SRCs. This insensitivity might be related to low interfacial shear strength at PE/PE interfaces, enabling process optimization while maintaining mechanical performance.

Cooling rate↗
DOE OSTI2026

Carbon doping in GeTe drives differences in local structure and properties

Advances in low-power, energy-efficient information storage and computing require understanding and controlling the atomic and nanoscale structures of functional materials, such as phase-change materials. Phase-change memory technology enables nonvolatile, low-power memory in devices by storing information through reversible changes in a phase-change material's atomic structure (i.e., transformations between amorphous and crystalline phases) that have corresponding changes in properties, including electronic resistivity and optical reflectivity. Here, we apply complementary X-ray absorption spectroscopy and X-ray pair distribution function analyses to experimentally identify the local- and medium-range atomic structure differences of GeTe and C-doped GeTe thin films. Upon controlled heating, composition- and temperature-dependent atomic structure evolution in GeTe and C-doped GeTe films shows differences in bonding behavior and local structure that directly influence crystallization onset temperature. We find that the introduction of C interrupts Ge–Ge bonds in amorphous GeTe, altering the as-deposited structure to be more similar to the distorted rocksalt structure of crystalline α–GeTe. The change alters the response of the amorphous atomic structure to heating and also lowers the crystallization onset temperature, from 230 °C in GeTe to 220 °C in the C-doped film. The combined insights from both X-ray techniques provide understanding of structural transformations that enables the development and optimization of next-generation memory and computing materials.

36 MATERIALS SCIENCE↗
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

Development and Evaluation of Embrittlement Resistant Alloys for Advanced LWR Cladding

The primary challenge preventing the deployment of reactor designs that leverage wrought FeCrAl as an advanced light-water reactor (LWR) cladding is irradiation hardening and embrittlement. Wrought FeCrAl alloys experience a loss of post-irradiation ductility and fracture toughness under low-temperature neutron irradiation (< 350°C) resulting from the combined effects of dislocation loop formation and the irradiation-enhanced precipitation of Cr-rich alpha-prime precipitates throughout the microstructure. Over the past decade, significant improvements in our understanding regarding the effect of Cr and Al content have been enabled through neutron irradiations, but even optimized wrought alloys such as C26M remain vulnerable to brittle failure during storage, transportation, and handling following irradiation. This report summarizes multiple irradiation campaigns initiated over the past 7 years and provides recent insights into the effect of (1) minor alloying composition, (2) alloy processing strategy, (3) crystal structure, and (4) the use of advanced oxide dispersion strengthened (ODS) alloys. The first irradiation campaign irradiated body-centered cubic FeCrAl alloy C26M with and without elements such as Mo and Y. It also included face-centered cubic alumina-forming austenitic (AFA) alloys that have a steeper neutronic penalty but that have far superior high-temperature strength than wrought FeCrAl. Finally, this campaign compared conventional wrought C26M with the same material produced using modern powder metallurgical hotisostatic pressing (PM-HIP). The second irradiation campaign, initiated in 2018, irradiated wrought and FeCrAl-ODS materials up to 50 dpa to assess their irradiation resistance to high-doses. The results of this work indicate that wrought C26M suffers significant degradation following neutron irradiation at LWR-relevant temperatures regardless of Mo and Y content. The PM-HIP variant appears to be more resilient to irradiation-induced ductility loss in comparison with wrought variants but still loses almost all ductility by 8 dpa. Notably, the AFAs retained at least 5% total elongation after irradiation at the same dose and temperature condition, although additional deterioration may be expected at higher fluence levels. FeCrAl-ODS materials show the greatest promise as a transformative longer-term accident-tolerant fuel (LT-ATF) cladding material. The extruded 106ZY10C alloy retained 10% total elongation after 16 dpa irradiation and retained 8% total elongation after 50 dpa. Additional challenges with FeCrAl-ODS alloys also remain, including optimization of end cap joining methodologies, enhancement of fracture toughness, scaling production to prove economic viability, and optimizing post-pilger heat treatments to maximize ductility and irradiation resistance. This series of irradiations demonstrates the efficacy of the LWR testbed available within the United States for the rapid irradiation and down-selection of LT-ATF candidates. As the country rapidly accelerates its timeline for the deployment of advanced reactor concepts, the effective utilization of the High Flux Isotope Reactor for separate-effects style irradiations should continue to be prioritized to answer the final questions pertaining to LT-ATF candidates necessary for the deployment of advanced boiling water reactors and small modular reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Sm 2 Ru 3 Sn 5 : A Noncentrosymmetric Cubic Member of the Ln 2 M 3 X 5 Family

An optimized synthetic method is presented for Sm 2 Ru 3 Sn 5 and investigate its physical properties and electronic structure. Sm 2 Ru 3 Sn 5 is prepared by arc-melting stoichiometric ratios of the elements and is confirmed by single crystal and powder X-ray diffraction. An antiferromagnetic transition is observed at T N = 3.8 K. A modified Curie-Weiss fit to the data in the range 50–150 K yields a Curie-Weiss temperature: θ CW = −36.6 K and an effective magnetic moment: μ eff = 0.83 μ B , in agreement with a Sm 3+ oxidation state. Field-dependent magnetization up to H = 7 T at 2 K shows a maximum response of 0.06 μ B , which is significantly lower than the expected Sm 3+ saturation moment (0.71 μ B ). Resistivity measurements indicate metallic behavior, and analysis of the magnetic entropy from the heat capacity reveals a doublet ground state due to crystal electric field splitting. The electronic structure and density of states are calculated with density function theory and further supported by the local density approximation with dynamical mean-field theory. Finally, the experimental and computational results highlight localized Sm 3+ moments and suggest a possible interplay between Ruddelman–Kitel–Kasuya–Yosida and Kondo interactions, positioning Sm 2 Ru 3 Sn 5 as a promising material for studying topology and complex physical phenomena.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Agentic framework for programmatic crystal structure generation using a fine-tuned worker–supervisor large language model

Platinum group metals (PGMs) underpin many catalytic technologies but face severe supply constraints, motivating the search for alternative materials and computational methods to accelerate discovery. While atomistic simulation tools such as Pymatgen and ASE have streamlined structure manipulation, they require detailed inputs, limiting accessibility for experimentalists and slowing early-stage exploration. Here, in this study, we present an AI-driven agentic framework that orchestrates worker–supervisor large language models (LLMs). The worker translates natural-language prompts of varying abstraction into valid crystallographic structures using a compact LLM fine-tuned with low-rank adaptation on a curated text–code–CIF dataset, emphasizing energy-efficient training. Benchmarking against the baseline CodeGen-350M-mono model shows that fine-tuning reduces hallucination rates from 100% to as low as 5% and improves structural match accuracy to up to 82% for fully specified inputs. Accuracy declines with decreasing prompt detail but remains nontrivial even when only stoichiometry and space group are provided, underscoring the LLM’s capacity for crystallographic inference. The supervisor Claude LLM evaluates the outputs and triggers iterative refinement through the worker’s built-in structure manipulation capabilities (e.g., supercell scaling, strain, vacancy, and substitution operations). We further demonstrate use cases for technologically relevant catalysts, including IrO 2 , pyrochlore Pb 2 Ir 2 O 7 , Ni 2 FeO 4 , and Ni 3 Mo, where the framework generates physically consistent structures that can be refined via geometry optimization. This work introduces a low-energy, language-driven pathway for integrating human and machine intelligence in materials design, paving the way for AI-assisted synthesis planning and high-throughput screening of complex oxides.

AI agent↗
DOE OSTI2025

Phase-field predictions of the influence of cooling rates during AM on the Evolution of Microstructures in Nickel-Based Single Crystal Superalloys

Additive manufacturing of single crystals made of Ni-based superalloys offers major cost savings for gas turbine engines with the inclusion of internal cooling channels. However, the lack of understanding of the effect of transient thermal conditions on solidification grain structure during additive manufacturing hinders the potential for process control to maintain the single crystal quality. The use of high-fidelity simulations through high performance computing to predict the evolution of the solidification microstructure will enhance the abilities to tailor the microstructures through process optimization. Phase field simulations are used to determine the effect local thermal conditions and defects on the stability of the solidification morphology, specifically with respect to the onset of columnar-to-equiaxed transition that results in the loss of the single crystal. The results are expected to be instrumental for developing future surrogate models to speed up the integration of design and manufacturing of turbine blades under the harsh in-service conditions.

36 MATERIALS SCIENCE↗