SEARCH · Engineering Papers
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.
Crystal growth procedures, fabrication techniques, and theoretical analysis were developed in order to make GaAlAs-GaAs solar cell structures which exhibit high performance at air mass 0 illumination and high temperature conditions.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
The University of Chicago instrument on the Heliocentric spacecraft (MEH experiment) will measure the energy spectrum of cosmic electrons in the range 5-400 MeV. In addition, the energy spectra and relative abundances of nuclei from protons to the iron group, with energies ranging from 30 MeV/n to 15 GeV/n, will be determined. Primary scientific objectives involve the study of the long and short term variability of these components as a probe of the structure of the heliosphere. Particles are identified by multiparameter analysis using the pulse height analyzed signals from eight active detectors - silicon solid state, plastic and crystal scintillators are solid and gas Cerenkov counters. Data return is optimized by a three level priority logic scheme.
NASA LeRC has a broad, active cathode technology development program in which both experimental and theoretical studies are being employed to further development of thermionic cathodes for use as electron sources in vacuum devices for communications and other space applications. One important type of thermionic cathode under development is the alkaline-earth oxide-coated (BaO, SrO, CaO) cathode. Significant improvements in the emission characteristics of this cathode have been obtained through modification of the chemical composition and morphology of the oxide coating, with the best result thus far coming from the addition of In2O3 and Sc2O3. Whereas the In2O3 produces a finer, more uniform particle structure, the exact chemical state and role of the Sc2O3 in the emission enhancement is unknown. The purpose of this cooperative agreement is to combine the studies of the surface chemistry and electron emission at NASA LeRC of chemically modified oxide coatings with a study of the thermochemistry and crystal structure using X-ray diffraction equipment and expertise at Clark Atlanta University (CAU). The study at CAU is intended to provide the description and understanding of the structure and thermochemistry needed for further improvement and optimization of the modified coatings. A description of the experimental procedure, preliminary X-ray diffraction test results, together with the design of an ultrahigh vacuum chamber necessary for high temperature thermochemistry studies will be presented.
NASA's Biological Crystal Growth Program (BCG) on the International Space Station (ISS) will consist of two phases. The first phase is during assembly of the ISS and will accommodate generic payloads that currently fly in the orbiter middeck. The second phase is after assembly of the ISS is complete and BCG payloads will occupy part of the Biotechnology Facility aboard the ISS. During both phases of the program, there will be two types of BCG payloads. One type will emphasize the production of crystals for structure determination back on Earth and will have high capacity for screening crystallization conditions. The second type of payload will be designed to study the crystallization process with the primary aim of developing new methods to further optimize the use of the microgravity environment. Beginning immediately, access to the BCG program for Guest Investigators is simplified. Access to all BCG hardware for Guest Investigators will be coordinated through one office at NASA's Marshall Space Flight Center. Details of how to obtain access to microgravity, the hardware available, and the operational aspects of the program will be described.
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)
During the 2001-2002 award period, we performed research on Pt/Ti/bare 6H-SiC and bare 4H-SiC interfaces in order to identify their electronic properties as a function of surface preparation. The overall aim of this work is to optimize the electronic properties of metal contacts to SiC as well as the active SiC material itself as a function of surface preparation and subsequent processing. Initially, this work has involved identifying bare surface, subsurface, and metal induced gap states at the metal-SiC contact and correlating energies and densities of deep levels with Schottky barrier heights. We used low energy electron-excited nanoluminescence (LEEN) spectroscopy, X-ray photoemission spectroscopy (XPS), and Secondary Ion Mass Spectrometry (SIMS) in order to correlate electronic states and energy bands with chemical composition, bonding, and crystal structure. A major development has been the discovery of polytype transformations that occur in 4H-SiC under standard microelectronic process conditions used to fabricate SiC devices. Our results are consistent with the stacking fault generation, defect formation, and consequent degradation of SiC recently reported for state-of-the-art ABB commercial diodes under localized electrical stress. Our results highlight the importance of -optimizing process conditions and material properties - anneal times, temperatures and doping to control such structural changes within epitaxial SiC layers. Thus far, we have established threshold times and temperatures beyond which 4H-SiC exhibits 3C-SiC transformation bands for a subset of dopant concentrations and process conditions. On the basis of this temperature time behavior, we have been able to establish an activation energy of approximately 2.5 eV for polytype transformation and dislocation motion. Work continues to establish the fundamental mechanisms underlying the polytype changes and its dependence on material parameters.
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)
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.