Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “B. Metallic material”

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

SiC or B4C-B/low strategic element content composite

Advancements in materials technologies are needed to provide the aerospace industry with alternate material options in the event of future strategic metal shortages and to optimize performance of engines. Composite materials are promising candidates for such an application. The potential of SiC and B4C-B filament reinforced low strategic element iron-base alloy content composites or use at 760 to 870 C is considered. A limiting factor towards developing this material for high temperature use has been the reaction between the filament and matrix material during fabrication and service which degrades filament strength. A low temperature fabrication process to limit filament/matrix reaction is being developed which involves the use of hollow cathode sputtering to coat the filaments with iron-base alloys of various compositions. An investigation is being conducted to determine the interfacial reaction effects of SiC and B4C-B filaments with iron-base alloys to develop an understanding of filament/matrix alloy compatibility for 760 to 870 C service.

Petrasek, D. W.↗

Examples of X-Ray Characterization Techniques in Energy Storage Research

Lithium-ion batteries have revolutionized the portable electronics and transportation sectors. Their performance is often critically dependent on the crystal structures of the anode and cathode electrode materials, which must enable the transport and reversible storage of lithium ions into and out of the lattice. Because lithium is a low-Z element, characterization of materials for lithium-ion batteries can be particularly challenging. Regardless, X-ray techniques enable analysis of material structures to better understand how battery materials perform and degrade, particularly when combined with other materials characterization and electrochemical characterization techniques. While X-ray techniques are most often used in battery research for phase identification of crystal structures, X-ray characterization techniques are also used for a wide variety of other purposes. I will discuss several examples from my research with various collaborators on several projects that highlight the impact that X-ray characterization techniques can have on battery research. The first example will focus on low-temperature microwave-assisted solvothermal synthesis of vanadium-doped LiFePO4 cathode materials for lithium-ion batteries. (1,2) Through a combination of electrochemical and materials characterization, we determined that low temperature synthesis resulted in metastable phases that enabled incorporation of higher dopant levels than resulting from high-temperature synthesis of thermodynamically stable phases. Rietveld refinement of X-ray diffraction data enabled understanding of how lattice parameters changed with doping levels and synthesis temperature. X-ray absorption near edge spectroscopy enabled understanding of the vanadium and iron oxidation states to confirm how vacancies in the structure caused by doping were charge compensated. This was important to understand because the literature suggests doping can improve LiFePO4 electrical conductivity, which improves battery charge and discharge rates. The second example will focus on understanding residual strain in lithium metal anodes. Lithium-ion batteries typically use graphite anodes, but the charge-storage capacity can be theoretically improved ~10x by using lithium metal as the anode material instead. However, lithium anodes suffer from growth of high-aspect-ratio features, such as dendrites, that can pierce nanoporous polymer separators and lead to short circuits and fires. External pressure is commonly applied to cells to enable better morphological control. We hypothesized that applied pressure may promote strain and possibly work hardening during electrochemical cycling, which motivated us to look for evidence of residual strain in lithium metal cycled under applied pressure using X-ray diffraction and sin2(..psi..) analysis. We found that lithium electrodeposited under high pressure exhibited in-plane compressive strain and that that lithium electrodeposited under low pressure did not. (3) The residual strain that accompanies electrodeposition under high pressure may lead to work hardening, which may explain how a soft metal like lithium can puncture separators and why higher pressure does not always decrease short circuits. (4-6) References: 1) Harrison, K. L.; Manthiram, A. Microwave-Assisted Solvothermal Synthesis and Characterization of Metastable LiFe1- x (VO) x PO4 Cathodes. Inorganic chemistry 2011, 50(8), 3613-3620. 2) Harrison, K. L.; Bridges, C. A.; Paranthaman, M. P.; Segre, C. U.; Katsoudas, J.; Maroni, V. A.; Idrobo, J. C.; Goodenough, J. B.; Manthiram, A. Temperature Dependence of Aliovalent-Vanadium Doping in LiFePO4 Cathodes. Chemistry of Materials 2013, 25(5), 768-781. 3) Rodriguez, M. A.; Harrison, K. L.; Goriparti, S.; Griego, J. J.; Boyce, B. L.; Perdue, B. R. Use of a Be-Dome Holder for Texture and Strain Characterization of Li Metal Thin Films via Sin2 (..psi..) Methodology. Powder Diffraction 2020, 35(2), 89-97. 4) Jungjohann, K. L.; Gannon, R. N.; Goriparti, S.; Randolph, S. J.; Merrill, L. C.; Johnson, D. C.; Zavadil, K. R.; Harris, S. J.; Harrison, K. L. Cryogenic Laser Ablation Reveals Short-Circuit Mechanism in Lithium Metal Batteries. ACS Energy Letters 2021, 6(6), 2138-2144. 5) Harrison, K. L.; Merrill, L. C.; Long, D. M.; Randolph, S. J.; Goriparti, S.; Christian, J.; Warren, B.; Roberts, S. A.; Harris, S. J.; Perry, D. L. Cryogenic Electron Microscopy Reveals That Applied Pressure Promotes Short Circuits in Li Batteries. Iscience 2021, 24(12). 6) Harrison, K. L.; Goriparti, S.; Merrill, L. C.; Long, D. M.; Warren, B.; Roberts, S. A.; Perdue, B. R.; Casias, Z.; Cuillier, P.; Boyce, B. L. Effects of Applied Interfacial Pressure on Li-Metal Cycling Performance and Morphology in 4 M LiFSI in DME. ACS Applied Materials & Interfaces 2021, 13(27), 31668-31679.

batteries↗

Infrared nanospectroscopy characterization of metal oxide photoresists

Implementation of extreme ultraviolet (EUV) lithography in high-volume semiconductor manufacturing requires a reliable and scalable EUV resist platform. A mechanistic understanding of the pros and cons of different EUV resist materials is critically important. However, most material characterization methods with nanometer resolution use an x-ray photon or electron beam as the probe, which often cause damage to the photoresist film during measurement. Here, we illustrated the use of non-destructive infrared nanospectroscopy [or nano-Fourier-Transform infrared spectroscopy (nano-FTIR)] to obtain spatially resolved composition information in patterned photoresist films. Further, clear evidence of exposure-induced chemical modification was observed at a spatial resolution down to 40 nm, well below the diffraction limit of infrared light. With improvements, such a nano-FTIR technique with nanoscale spatial resolution, chemical sensitivity, and minimal radiation damage can be a promising candidate for the fundamental study of material properties relevant to EUV lithography.

36 MATERIALS SCIENCE↗

Ingeniously enhanced ferromagnetism in chemically-reduced 2D Ti 3 C 2 T X MXene

Chemical reduction is a facile and cost-effective technique for the modulation of the physical and chemical properties of nanomaterials. In this work, we demonstrate an enhancement of the magnetic behavior of Ti 3 C 2 T x MXene after its chemical reduction via l-ascorbic acid treatment. Small ferromagnetic loops have been observed below 50 K for Ti 3 C 2 T x prepared by hydrofluoric acid (HF) etching of Al layers from Ti 3 AlC 2 . Such a ferromagnetic ordering of spins was significantly enhanced via a chemical reduction of Ti 3 C 2 T x with l-ascorbic acid. Ferromagnetic hysteresis loops were observed for reduced Ti 3 C 2 T x (r-Ti 3 C 2 T x ) up to 150 K indicating a significant upshift of the paramagnetic to the ferromagnetic transition temperature, pushing towards room temperature. The enhancement of ferromagnetism and upshift of the ferromagnetic transition temperature could be attributed to the localized unpaired electron in Ti-3d orbital of the r-Ti 3 C 2 T x crystal and increase in the number of unsaturated Ti atoms upon l-ascorbic acid treatment. Chemical reduction via l-ascorbic acid treatment shows a promising pathway towards the modulation and enhancement of magnetism in various MXene materials for the development of 2D metallic soft ferromagnets and spintronic devices.

36 MATERIALS SCIENCE↗

Kinetics of the Carbothermic Reduction of Manganese Oxide from Slag

Experiments were performed using a range of test conditions to elucidate the rate controlling step during the reaction of liquid iron-carbon droplets and slags containing manganese oxide. Four conditions were tested in the system: initial MnO content in the slag (5, 10, and 15 wt pct), initial carbon content of the metal (1, 2.5, 4.3 wt pct), initial droplet mass (0.5, 1.0, and 1.5 g), and reaction temperature (1823 K [1550 °C], 1873 K [1600 °C], and 1923 K [1650 °C]). Data were collected using the Constant Volume Pressure Increase (CVPI) technique which tracked the continuous pressure increase in the sealed furnace over time. Samples were quenched at the end of each experiment and chemistry was checked using LECO Carbon Analysis and ICP (Inductively Coupled Plasma) for manganese. The rate of reaction can be broken into a faster initial period related to internal CO formation, and a slower second reaction controlled by a complex mechanism involving transport of oxygen from slag to metal via CO{sub 2} and decomposition of the CO{sub 2} at the gas–metal interface.

36 MATERIALS SCIENCE↗

Unconventional polaronic ground state in superconducting LiTi 2 O 4

Geometrically frustrated lattices can display a range of correlated phenomena, ranging from spin frustration and charge order to dispersionless flat bands due to quantum interference. One particularly compelling family of such materials is the half-valence spinel LiB 2 O 4 materials. On the B-site frustrated pyrochlore sublattice, the interplay of correlated metallic behavior and charge frustration leads to a superconducting state in LiTi 2 O 4 and heavy fermion behavior in LiV 2 O 4 . To date, however, LiTi 2 O 4 has primarily been understood as a conventional BCS superconductor despite a lattice structure that could host more exotic ground states. Here, we present a multimodal investigation of LiTi 2 O 4 , combining ARPES, RIXS, proximate magnetic probes, and ab-initio many-body theoretical calculations. Our data reveals a novel mobile polaronic ground state with spectroscopic signatures that underlie co-dominant electron-phonon coupling and electron-electron correlations also found in the lightly doped cuprates. The cooperation between the two interaction scales distinguishes LiTi 2 O 4 from other superconducting titanates, suggesting an unconventional origin to superconductivity in LiTi 2 O 4 . Our work deepens our understanding of the rare interplay of electron-electron correlations and electron-phonon coupling in unconventional superconducting systems. In particular, our work identifies the geometrically frustrated, mixed-valence spinel family as an under-explored platform for discovering unconventional, correlated ground states.

36 MATERIALS SCIENCE↗

Nanoscale Mg–B via Surfactant Ball Milling of MgB 2 : Morphology, Composition, and Improved Hydrogen Storage Properties

Metal borides have attracted the attention of researchers due to their useful physical properties and unique ability to form high hydrogen-capacity metal borohydrides. We demonstrate improved hydrogen storage properties of a nanoscale Mg–B material made by surfactant ball milling MgB 2 in a mixture of heptane, oleic acid, and oleylamine. Transmission electron microscopy data show that Mg–B nanoplatelets are produced with sizes ranging from 5 to 50 nm, which agglomerate upon ethanol washing to produce an agglomerated nanoscale Mg–B material of micron-sized particles with some surfactant still remaining. X-ray diffraction measurements reveal a two-component material where 32% of the solid is a strained crystalline solid maintaining the hexagonal structure with the remainder being amorphous. Fourier transform infrared shows that the oleate binds in a “bridge-bonding” fashion preferentially to magnesium rather than boron, which is confirmed by density functional theory calculations. The Mg–B nanoscale material is deficient in boron relative to bulk MgB 2 with a Mg–B ratio of ~1:0.75. The nanoscale MgB 0.75 material has a disrupted B–B ring network as indicated by X-ray absorption measurements. Hydrogenation experiments at 700 bar and 280 °C show that it partially hydrogenates at temperatures 100 °C below the threshold for bulk MgB 2 hydrogenation. In addition, upon heating to 200 °C, the H–H bond-breaking ability increases ~10-fold according to hydrogen–deuterium exchange experiments due to desorption of oleate at the surface. This behavior would make the nanoscale Mg–B material useful as an additive where rapid H–H bond breaking is needed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Noncollinear spin structure with weak ferromagnetism in NbMnP

NbMnP is a metallic material, which consists of the zigzag chains of Mn moments along the b axis. The magnetic susceptibility as well as the resistivity shows an anomaly at 233 K, which indicates an antiferromagnetic phase transition. Our neutron powder diffraction experiment reveals that the magnetic structure is a Q=0 structure and noncollinear with an easy plane anisotropy perpendicular to the b axis, and the a- and c-axis magnetic components align antiferromagnetically and ferromagnetically along the zigzag chain direction, respectively. The ordered moment is 1.2μ B , which is reduced probably due to the itineracy of the Mn moments. A localized picture model suggests that the Q=0 magnetic structure is formed by the frustration among several exchange couplings. A weak ferromagnetic component is also present in the antiferromagnetic phase, which is considered to be caused by a Dzyaloshinskii-Moriya interaction.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Method of joining metallic and composite components

A method is provided for joining a metallic member to a structure made of a composite matrix material. One or more surfaces of a portion of the metallic member that is to be joined to the composite matrix structure is provided with a plurality of outwardly projecting studs. The surface including the studs is brought into engagement with a portion of an uncured composite matrix material so that fibers of the composite matrix material intertwine with the studs, and the metallic member and composite structure form an assembly. The assembly is then companion cured so as to join the metallic member to the composite matrix material structure.

Semmes, Edmund B.↗

Surface treatment for lithium battery electrode materials

Electrode materials for electrochemical cells and batteries and methods of producing such materials are disclosed herein. A method of preparing an active lithium metal oxide material suitable for use in an electrode for a lithium electrochemical cell comprises the steps of: (a) contacting the lithium metal oxide material with an aqueous acidic solution containing one or more metal cations; and (b) heating the so-contacted lithium metal oxide from step (a) to dryness at a temperature below 200° C. The metal cations in the aqueous acidic solution comprise one or more metal cations selected from the group consisting of an alkaline earth metal ion, a transition metal ion, and a main group metal ion.

25 ENERGY STORAGE↗

Lithium battery electrodes

Electrode materials for electrochemical cells and batteries and methods of producing such materials are disclosed herein. The electrode materials comprise an active lithium metal oxide material prepared by: (a) contacting the lithium metal oxide material with an aqueous acidic solution containing one or more metal cations; and (b) heating the so-contacted lithium metal oxide from step (a) to dryness at a temperature below 200° C. The metal cations in the aqueous acidic solution comprise one or more metal cations selected from the group consisting of an alkaline earth metal ion, a transition metal ion, and a main group metal ion.

25 ENERGY STORAGE↗

Opposed-Flow Flame Spread in a Narrow Channel Apparatus over Thin PMMA Sheets

Flame spread tests have been conducted over polymethylmethacrylate (PMMA) samples in San Diego State University's Narrow Channel Apparatus (SDSU NCA). The Narrow Channel Apparatus (NCA) has the ability to suppress buoyant flow in horizontally spreading flames, and is currently being investigated as a possible replacement or complement to NASA's current material flammability test standard for non-metallic solids, NASA-STD-(I)-6001B Test 1. The buoyant suppression achieved with a NCA allows for tests to be conducted in a simulated microgravity atmosphere-a characteristic that Test 1 lacks since flames present in Test 1 are buoyantly driven. The SDSU NCA allows for flame spread tests to be conducted with varying opposed flow oxidizer velocities, oxygen percent by volume, and total pressure. Also, since the test sample is placed symmetrically between two confining plates so that there is a gap above and below the sample, this gap can be adjusted. This gap height adjustment allows for a compromise between heat loss from the flame to the confining boundaries and buoyancy suppression achieved by those boundaries. This article explores the effect gap height has on the flame spread rate for 75 μm thick PMMA at 1 atm pressure and 21% oxygen concentration by volume in the SDSU NCA. Flame spread results from the SDSU NCA for thin cellulose fuels have previously been compared to results from tests in actual microgravity at various test conditions with the same sample materials and were found to be in good agreement. This article also presents results from the SDSU NCA for PMMA at 1 atm pressure, opposed oxidizer velocity ranging from 3 to 35 cm/s, oxygen concentration by volume at 21%, 30 %, and 50% and fuel thicknesses of 50 and 75 μm. These results are compared to results obtained in actual microgravity for PMMA obtained at the 4.5s drop tower of MGLAB in Gifu, Japan, and the 5.2s drop tower at NASA's Zero-Gravity Research Facility in Cleveland, OH. This comparison confirms that at 1 atm pressure, the SDSU NCA successfully simulates microgravity for not only thin cellulose fuels, but also for thin PMMA sheets as well. This further supports the idea that the NCA is a viable option to complement or replace NASA's Test 1 for material flammability testing. Tests with thick fuels will be conducted in the future to further characterize the SDSU NCA.

flame spread↗

Tensile behavior of Fe-40Al alloys with B and Zr additions

Both Fe-40Al and Fe-40Al-0.1Zr with and without B were produced by the hot extrusion of powdered metal. Tensile properties were determined from room temperature to 1100 K in air. All of the materials possessed some ductility at room temperature, and addition of B caused an increase in ductility and a change in fracture mode from intergranular to transgranular cleavage. At high temperatures, failure was caused primarily by the formation of grain boundary cavities in all of the alloys. The effect of Zr addition was unclear because of the complexity of the various microstructures. Comparison of air and vacuum testing at high temperatures revealed that an apparent oxidation assisted mechanism reduced high temperature ductility in these alloys, especially at 900 K.

Gaydosh, D. J.↗

Computational Modeling Demonstrates the Joint Clearance Effect on the Tensile Strength of Solder Joints.

Soldered joints can be made with a wide range of base materials and filler metals that allow the assembly to meet its performance and reliability requirements. Structural solder joints have, as their foremost requirement, to provide mechanical attachment between base material structures. The joint is typically subjected to one, or a combination of, three loading configurations: (a) tensile or compressive force, (b) shear force, or (c) peel force. Solder filler metals and in particular, the so-called “soft solders” based on tin (Sn), lead (Pb), and indium (In), generally have a bulk strength that is less than that of the base materials. Finally, deformation occurs largely in the solder when the joint is subjected to an applied force.

36 MATERIALS SCIENCE↗

Surface Modified Fly Ash for Value Added Products (SuMo Fly Ash)

The objective of this project is to develop a technology to encapsulate coal fly ash particles in sulfurized vegetable oil, enhancing physical and mechanical properties of the fly ash as a filler material when applied in commercial products. This project significantly advances the knowledge base and technology for synthesizing coated fly ash particles for application in different polymer matrices to increase cross-linking, compatibility, air-entrainment and to decrease the leaching potential of metals of concern. Specifically, the project focuses on (a) collection, and characterization of fly ash material from coal power plants; (b) development of Sulfurized Vegetable Oil (SVO) modified (SuMo) fly ash and their detailed characterization; (c) evaluation of the mechanical properties of the SuMo fly ash incorporated plastic and rubber composites for potential replacement of CaCO 3 or carbon black as filler materials with SuMo fly ash in plastic and rubber composites; (d) leaching potential of SuMo fly ash and SuMo fly ash-incorporated plastic and rubber composites to determine the rate of release of trace metals. The major findings of this project are: i) Sulfurized Vegetable Oil coated fly ash (SuMo fly ash) was successfully prepared with a particle size of ≤ 45 micron which exhibited hydrophobicity of contact angle > 90°; ii) the coating reduces leaching of metals (e.g., B, Cr) from fly ash when exposed to water; iii) incorporation of SuMo fly ash increases thermal stability and yield strength of plastics; iv) SuMo coating helps disperse fly ash particles into cured rubbers, natural rubber; (v) SuMo fly ash-incorporated plastic/cured rubbers compounds protects against leaching of toxic elements.

20 FOSSIL-FUELED POWER PLANTS↗

Electronic structures and magnetic performance related to spintronics of Sr{sub 0.875}Ti{sub 0.125}S

In this study, we have used the first-principle approaches of density functional theory to determine the electronic and magnetic performances related to spintronics applications of titanium (Ti)-substituted SrS such as Sr{sub 1−x}Ti{sub x}S at concentration x = 0.125. The results of electronic and magnetic properties by using TB–mBJ potential revealed that Sr{sub 0.875}Ti{sub 0.125}S compound is half-metallic ferromagnetic with a half-metallic gap and a spin polarization of 100%. The ferromagnetic nature of Sr{sub 0.875}Ti{sub 0.125}S material is confirmed by the ferromagnetic coupling between magnetic moments of Ti, Sr and S atoms. The total magnetic moment of integer Bohr magneton of 2 μ{sub B} and the large half-metallic gap of 1.245 suggest that this compound is a true half-metallic ferromagnet.

36 MATERIALS SCIENCE↗

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↗

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↗