Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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

Effect of 3D-printed surface textures on wear mechanism in 3-body abrasion of soil

This study systematically investigates the enhancement of wear resistance in 3D printed surface textures through both experimental and theoretical approaches. Three distinct surface morphologies (Smooth Surface, Surface with uniformly distributed Pits, and Surface with uniformly distributed Bumps) were fabricated using High-Impact Polystyrene, where the meso-scale textures were precisely controlled through the 3D printing process. Wear behavior was evaluated using a 3-body wear tester in an abrasive particle environment, analyzing the influence of surface textures under various operating conditions. Systematic wear tests revealed that optimally designed surface textures achieved a remarkable 77 % reduction in wear compared to the worst-performing sample. The wear mechanisms were comprehensively characterized through weight loss measurements, Scanning Electron Microscopy (SEM), and Energy Dispersive Spectroscopy (EDS) analyses, elucidating the surface morphology changes and their interaction with wear particles. Notably, the study identified how the geometric characteristics of surface textures influence the movement of wear particles and the distribution of contact stresses. Discrete element method simulations corroborated the experimental findings, providing theoretical validation for the enhanced wear resistance of the optimal structure. The high correlation between simulated wear patterns and experimental results validates the reliability of the proposed design methodology. In conclusion, these results demonstrate that 3D printed surface texturing offers a cost-effective and scalable approach to significantly improve wear resistance in engineering applications, presenting a practical alternative to conventional, high-cost surface engineering methods.

3-body abrasion↗
DOE OSTI2024

Scalable Surface Micro-Texturing of LLZO Solid Electrolytes for Battery Applications

A challenge for lithium lanthanum zirconate (LLZO)-based solid-state batteries is to increase the critical current density (CCD) to enable high current cycling. A promising strategy is to modify the LLZO surface morphology to provide a larger contact area with the Li metal. Here, a surface-textured thin LLZO electrolyte was prepared through an easily scalable process. The texturing process is a simple pressing of green LLZO tapes between micro-textured substrates. A variety of textures can be produced, depending on the type of substrate, and texturing can be on either one side or both sides. For this work, after pressing and sintering, several micro-patterns are formed on thin LLZO (~118 μm thick). The properties of the various samples were characterized to investigate the impact of surface texturing, and the most promising ones were selected for electrochemical testing in symmetrical lithium cells and full cells. Li symmetric cells using a coarse ridge-textured LLZO exhibit ~2.5 times increased CCD compared to planar non-textured LLZO, and a solid-state full cell shows stable cycling and improved rate performance. Finally, we believe this process offers a favorable trade-off of processing complexity vs structural optimization to maximize CCD.

25 ENERGY STORAGE↗
DOE OSTI2025

Exploring the Structural Origins of Optically Efficient One-Dimensional Lead Halide Perovskite Nanostructures

Metal halide perovskites have excellent optoelectronic properties. This study aims to determine how the optoelectronic properties of a model perovskite, cesium lead bromide (CsPbBr3), change with length and thickness in one dimension (1D). By examining the photophysics of CsPbBr3 quantum dots (QDs), nanowires (NWs), and nanorods (NRs), we observe the influence of confinement, exciton diffusion, and trapping on their optical properties. Our findings reveal that exciton diffusion to trap states limits the photoluminescence quantum yield (PLQY) of 1D CsPbBr3 in the weakly confined regime (8-14 nm) and explains their long-lived exciton dynamics, while enhanced radiative rates contribute to achieving near-unity PLQY in the strongly confined regime (<7 nm). Consequently, blue-emitting, 2.4 nm-thick CsPbBr3 NRs were 3.6X more emissive than the conventional CsPbBr3 QDs. This study underscores how structural optimization can improve the optoelectronic performance of CsPbBr3 and provides insight into the complex interplay of radiative and nonradiative processes in 1D ionic semiconductors.

Oddo, Alexander M↗
DOE OSTI2025

A Conjugated Oligoelectrolyte Exhibiting Room Temperature Spin-Correlated Radical Pair Character for Biological Sensing

We report a water-soluble conjugated oligoelectrolyte (COE) composed of carbazole-benzophenone, COE-CbzBP, that exhibits photogenerated spin-correlated radical pair (SCRP) behavior sensitive to static electric fields from DNA but not from lipid bilayers. The SCRP forms from a thermally activated, spin-polarized state enabled by partial π-conjugation disruption at the donor–acceptor (carbazole-benzophenone) nitrogen–carbon (N–C) junction, which facilitates a twisted intramolecular charge-transfer (TICT) geometry. This state minimizes the singlet–triplet energy gap (ΔE ST = 0.12 eV), radical–pair exchange coupling (J RP ∼ ΔE ST /2), and charge separation free energy (ΔG CS ) in both DNA (−0.19 eV) and lipid bilayers (−0.55 eV). Room-temperature continuous-wave electron paramagnetic resonance (CW-EPR) reveals a photogenerated spin-polarized singlet for COE-CbzBP that splits upon DNA association, consistent with modulation of J RP and hyperfine coupling (A x ), presumably via electric field-spin coupling. No spin-polarized signal was observed under dark, cryogenic conditions, or in liposomes, but was quenched by the spin trap 4-POBN. Transient absorption and spectroelectrochemistry confirmed magnetic-field sensitive long-lived excited-state absorption features attributed to charge-separated states 3 [Cbz •+ -BP •– ]*, which were lengthened by DNA, and quenched in lipid bilayers and 4-POBN. Quantum chemical simulations show that planar geometries (lipid-like) increase ΔE ST by 0.31 eV compared to TICT-optimized structures. This geometry-dependent modulation explains the absence of SCRP signatures in rigid environments, underscoring the importance of TICT states, minimized ΔE ST , and favorable ΔG CS for achieving room-temperature SCRP generation. These findings establish design principles for TICT-enabled molecules exhibiting qubit-like behavior that operate under ambient and biologically relevant conditions, with direct implications for quantum information science (QIS).

Aromatic compounds↗
DOE OSTI2024

Synergizing superwetting and architected electrodes for high-rate water splitting

Water splitting is one of the most promising technologies for generating green hydrogen. To meet industrial demand, it is essential to boost the operation current density to industrial levels, typically in the hundreds of mA cm -2 . However, operating at these high current densities presents significant challenges, with bubble formation being one of the most critical issues. Efficient bubble management is crucial as it directly impacts the performance and stability of the water splitting process. Superwetting electrodes, which can enhance aerophobicity, are particularly favorable for facilitating bubble detachment and transport. By reducing bubble contact time and minimizing the size of detached bubbles, these electrodes help prevent blockage and maintain high catalytic efficiency. Here, in this review, we aim to provide an overview of recent advancements in tackling bubble-related issues through the design and implementation of superwetting electrodes, including surface modification techniques and structural optimizations. We will also share our insights into the principles and mechanisms behind the design of superwetting electrodes, highlighting the key factors that influence their performance. Our review aims to guide future research directions and provides a solid foundation for developing more efficient and durable superwetting electrodes for high-rate water splitting.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Electrochemical reduction of ammonia-captured CO 2 to CO over a nickel single-atom catalyst

Carbon reactive capture and conversion offers a sustainable route to valuable chemicals and fuels while aiding Green House Gas (GHG) reduction. Direct electrochemical conversion of capture solutions like bicarbonate avoids the energy demands of conventional CO 2 regeneration. Ammonium bicarbonate (NH 4 HCO 3 ) is particularly attractive due to its low decomposition temperature and ability to supply in situ CO 2 from dilute sources without requiring purified CO 2 . Meanwhile, single-atom catalysts (SACs) with nitrogen-coordinated metal sites further enhance CO 2 reduction efficiency using Earth-abundant materials. In this study, we demonstrate a nickel single-atom catalyst (Ni-SAC)-based electrolyzer that utilizes NH 4 HCO 3 as the CO 2 source, achieving significantly improved CO production performance compared to the conventional silver cathodes used in the CO 2 reduction reaction (CO 2 RR) to produce CO. The Ni-SAC cathode exhibited a Faradaic efficiency of 60.1% for CO production at −200 mA cm −2 , while the silver cathode achieved a Faradaic efficiency of only 2%, likely due to ammonium-induced poisoning. Furthermore, the integration of a customized microporous layer onto the electrode significantly increased the Faradaic efficiency from 64% to 83% at −100 mA cm −2 , emphasizing the crucial role of electrode structure optimization in enhancing CO selectivity. These findings demonstrate a sustainable and economically viable strategy for green CO production directly from CO 2 capture solutions.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2026

Octave-spanning terahertz quarter-wave plates based on over-coupled Fabry–Pérot resonances in reflective metal–dielectric–metal metasurfaces

Compact devices for broadband polarization control in the terahertz (THz) regime remain challenging to realize due to the intrinsic phase dispersion associated with birefringent materials and resonant structures. Here, we demonstrate achromatic THz quarter-wave plates with near-unity efficiency and operating under a wide range of incidence angles, based on over-coupled metal–dielectric–metal reflective metasurfaces. The devices operate as single-port anisotropic Fabry–Pérot cavities, in which the phase dispersion of the over-coupled resonances is engineered to maintain an approximately constant relative phase delay between orthogonal field components over a broad frequency range. Four complementary metasurface designs operating at an incidence angle of 45° collectively cover the 0.25–3 THz frequency range accessible to a typical THz time-domain spectroscopy system. Each device exhibits an approximately octave-wide operational bandwidth, with axial ratios below 3 dB and polarization-conversion efficiencies exceeding 80% across most of the operating band. Systematic structural optimization suppresses coupling to diffraction and guided-wave modes, thereby further extending the usable bandwidth while preserving the required phase relationship between orthogonal polarizations. The demonstrated metasurfaces are compatible with wafer-scale fabrication processes, and experimental results show excellent agreement with simulations. These findings establish over-coupled reflective metasurfaces as a robust, scalable, and versatile platform for broadband THz polarization control.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Quantum Mechanical MP2 Study of the Electronic Effect of Nonplanarity on the Carbon Pyramidalization of Fullerene C 60

Among C 60 ’s diverse functionalities, its potential application in CO 2 sequestration has gained increasing interest. However, the processes involved are sensitive to the molecule’s electronic structure, aspects of which remain debated and require greater precision. To address this, we performed structural optimization of fullerene C 60 using the QM MP2/6–31G* method. The nonplanarity of the optimized icosahedron is characterized by two types of dihedral angles: 138° and 143°. The 120 dihedrals of 138° occur between two hexagons intersecting at C–C bonds of 1.42 Å, while the 60 dihedrals of 143° are observed between hexagons and pentagons at C–C bonds of 1.47 Å. NBO analysis reveals less pyramidal sp 1.78 hybridization for carbons at the 1.42 Å bonds and more pyramidal sp 2.13 hybridization for the 1.47 Å bonds. Electrostatic potential charges range from -0.04 a.u. to 0.04 a.u. on the carbon atoms. Second-order perturbation analysis indicates that delocalization interactions in the C–C bonds of 1.42 Å (143.70 kcal/mol) and 1.47 Å (34.98 kcal/mol) are 22% and 38% higher, respectively, than those in benzene. MP2/Def2SVP calculations yield a correlation energy of 13.49 kcal/mol per electron for C 60 , slightly higher than the 11.68 kcal/mol for benzene. However, the results from HOMO-LUMO calculations should be interpreted with caution. This study may assist in the rational design of fullerene C 60 derivatives for CO 2 reduction systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

High-Speed and Hysteresis-Free Near-Infrared Optical Hydrogen Sensor Based on Ti/Pd Bilayer Thin Films

Palladium (Pd) and titanium (Ti) exhibit opposite dielectric responses upon hydrogenation, with stronger effects observed in the near-infrared (NIR) region. Leveraging this contrast, we investigated Ti/Pd bilayer thin films as a platform for NIR hydrogen sensing—particularly at telecommunication-relevant wavelengths, where such devices have remained largely unexplored. Ti/Pd bilayers coated with Teflon AF (TAF) and fabricated via sequential electron-beam and thermal evaporation were characterized using optical transmission measurements under repeated hydrogenation cycles. The Ti (5 nm)/Pd (x = 2.5 nm)/TAF (30 nm) architecture showed a 2.7-fold enhancement in the hydrogen-induced optical contrast at 1550 nm compared to Pd/TAF reference films, attributed to the hydrogen ion exchange between the Ti and Pd layers. The optimized structure, with a Pd thickness of x = 1.9 nm, exhibited hysteresis-free sensing behavior, a rapid response time (t90 < 0.35 s at 4% H2), and a detection limit below 10 ppm. It also demonstrated excellent selectivity with negligible cross-sensitivity to CO2, CH4, and CO, as well as high durability, showing less than 6% signal degradation over 135 hydrogenation cycles. These findings establish a scalable, room-temperature NIR hydrogen sensing platform with strong potential for deployment in automotive, environmental, and industrial applications.

Chemistry↗
DOE OSTI2025

Modeling commercial-scale CO 2 storage in the gas hydrate stability zone with PFLOTRAN v6.0

Abstract. Safe and secure carbon dioxide (CO2) storage is likely to be critical for mitigating some of the most dangerous effects of climate change. In the last decade, there has been a significant increase in activity associated with reservoir characterization and site selection for large-scale CO2 storage projects across the globe. These prospective storage sites tend to be selected for their optimal structural, petrophysical, and geochemical trapping potential. However, it has also been suggested that storing CO2 in reservoirs within the CO2 hydrate stability zone (GHSZ), characterized by high pressures and low temperatures (e.g., Arctic or marine environments), could provide a natural thermodynamic barrier to gas leakage. Evaluating the prospect of commercial-scale, long-term storage of CO2 in the GHSZ requires reservoir-scale modeling capabilities designed to account for the unique physics and thermodynamics associated with these systems. We have developed the HYDRATE flow mode and the accompanying fully implicit parallel well model in the massively parallel subsurface flow and reactive transport simulator PFLOTRAN to model CO2 injection into the marine GHSZ. We have applied these capabilities to a set of CO2 injection scenarios designed to reveal the challenges and opportunities for commercial-scale CO2 storage in the GHSZ.

carbon storage↗
DOE OSTI2024

New systems in MOOSE

The Multiphysics Object-Oriented Simulation Environment (MOOSE) serves as a common library of classes between applications developed for advanced reactor analysis, fusion device engineering, spent fuel cask analysis, geochemistry studies, among other fields. These applications drive the development of the framework to meet their needs. Systems in MOOSE group capabilities that share a common purpose and generally common code. They can be leveraged by all downstream applications, providing extensive code re-use and shared maintenance. They facilitate the discovery by new users of the classes meeting at least partially their needs, and offer the same opportunities for customization as other systems. The addition of a new system to MOOSE opens new ways of solving or discretizing nonlinear problems, of performing distributed postprocessing, and a plethora of other needs. While new systems can be introduced in downstream applications rather than at the framework level, the framework team monitors common needs across the community and often triggers their addition. Documentation, training material, development needs can be centralized, limiting duplicated work across the community. The last three years have seen a large expansion in the capabilities of MOOSE. The supporting role of the framework in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has created numerous feature requests to support neutronics, thermal hydraulics, computational fluid dynamics and thermo-mechanics simulations in the Griffin, SAM, Pronghorn and Bison applications respectively. Similarly, laboratory-directed research and development (LDRD) projects in additive manufacturing, high-Reynolds flow simulations, structure optimization also necessitate an expansion of the framework capabilities. This summary reports on the new systems created in MOOSE, their design, their capabilities and some of the relevant interfaces.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Efficient derivative computation for unsteady fatigue-constrained nonlinear aero-structural wind turbine blade optimization

Gradient-based optimization offers significant efficiency advantages for wind turbine blade design, but its application has often been limited by the cost and accuracy of finite-difference derivative calculations, especially when fatigue constraints are considered. In this work, we systematically compare and evaluate four differentiation techniques, namely algorithmic differentiation, implicit differentiation, sparsity exploitation, and parallelization, to determine their effectiveness in computing accurate gradients through time-domain aero-structural simulations. By integrating these techniques with unsteady nonlinear aerodynamic and structural models, we develop software designed for accurate gradient computation. We show that combining these techniques addresses memory and runtime challenges associated with long simulations required by design load cases. Specifically, the most effective combination reduces derivative computation wall time by over an order of magnitude compared to finite differencing while maintaining superior accuracy. We demonstrate this approach in a proof-of-concept aero-structural optimization of a wind turbine blade that improves the cost of energy by 12.78 %. This comparative study establishes a viable approach for fatigue-aware blade design that balances computational efficiency with modeling accuracy.

17 WIND ENERGY↗
DOE OSTI2026

A Flexible Quasi-Static Mooring Design Optimization Method for Floating Structures

This paper presents a flexible and efficient design method for optimizing the mooring systems of floating structures. Mooring system optimization is challenging because of the strong nonlinearity of mooring system behavior and the many technical constraints that must be satisfied. Furthermore, different mooring configurations can have very different design spaces. While some successful examples of mooring design optimization exist in the literature, developing an optimization approach that can work across various mooring design problems is a larger challenge. We present such a method based on a flexible parameterization that allows a wide variety of mooring designs to be described by a list of variables, a quasi-static mooring model that provides efficient evaluation of a mooring design without directly considering mooring system dynamics, and an optimization framework that generates, evaluates, and adjusts the mooring design while considering user-specified constraints such as offset limits, strength safety factors, and seabed contact limits. We demonstrate the design optimization framework on four mooring design problems, each for a different type of mooring system. We compare the use of different design modes to simplify the optimization problem, showing that they can reduce the computation time by up to 75%. We also compare different optimization algorithms and find that the resulting computational speed can vary by up to 51 times. We perform a sensitivity study on one design and find that the local sensitivity of anchoring radius to water depth has a positive correlation of 0.29, but the global sensitivity shows large nonlinearities. Lastly, we perform a coupled dynamic analysis on one of the optimized designs and find that the predicted mean platform motions and mooring line tensions are within 1% of dynamic results and the extreme motions and tensions are within 14%. Lastly, we show that a DEA-Chain-Polyester mooring configuration is cost-optimal for the given design problem of the demonstrations, which aligns with general industry practice.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Integrated process-structure–property–performance optimization of cold-sprayed zinc coatings on AZ91 magnesium alloy

High-pressure die cast (HPDC) AZ91 magnesium alloy is widely used in automotive components such as transmission housings and brackets for its excellent strength-to-weight ratio. Zinc-based cold spray coatings can be applied selectively to vulnerable areas to enhance corrosion resistance, minimize galvanic coupling with dissimilar metals, and eliminate the need for full-surface oxide coatings, making the process more efficient and targeted. A comprehensive evaluation of 16 combinations of nitrogen carrier gas temperatures and pressures led to the identification of an optimal range of process parameters, yielding Zn coatings with porosity <0.5 % by area, wear rates reduced by a factor of two compared to uncoated AZ91, and adhesion strengths up to 35 MPa. The enhanced mechanical performance of the coating is attributed to the low porosity and the formation of a metallurgical bond at the coating-substrate interface. Corrosion studies using macroscale potentiodynamic polarization (PDP) and electrochemical impedance spectroscopy (EIS) revealed a significant decrease in corrosion rate and a shift to more noble corrosion potentials (ZCP) for coated substrates. Furthermore, the Zn cold-sprayed samples exhibited significantly lower corrosion-induced evolved hydrogen content compared to the base AZ91 substrate and AZ91 coated with industrial coatings, demonstrating that the Zn layer effectively protects the substrate from the corrosive environment. Overall, cold spray Zn coatings significantly improve the mechanical and corrosion performance of AZ91 Mg alloys, addressing key material challenges and enabling their broader use in automotive applications.

Corrosion performance↗
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↗