Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic gradient structure”

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
NASA NTRS1981

Metal-semiconductor interfacial reactions - Ni/Si system

X-ray photoelectron spectroscopy and channeling measurements with MeV He-4(+) ions have been used to probe the structure of the interface in the Ni/Si system. It is found that reactions occur where Ni is deposited on Si at 10 to the -10th torr: Si atoms are displaced from lattice sites, the Ni atoms are in an Si-rich environment, and the Ni/Si interface is graded in composition. Composition gradients are present at both interfaces in the Si/Ni2/Si/Ni system. For the Ni-Si system, cooling the substrate to 100 K slows down the reaction rate. The temperature dependence of the interfacial reactivity indicates the kinetic nature of metal-semiconductor interfaces.

Cheung, N. W.↗
DOE OSTI2025

Unraveling the Atomic Mechanism of the Crystalline Phase‐Dependent Structural Features and Special Spectral Design of α‐, β‐, and Ɛ‐Ga₂O₃

Atomic‐scale phase transformations profoundly influence the functional properties of Ga₂O₃ polymorphs. By combining irradiation experiments with microstructure characterization and theoretical approaches, phase‐specific energy‐dissipation pathways in α‐, β‐, and ε‐Ga₂O₃ are uncovered and strategies for targeted property design are outlined. Competing antiphase boundaries (APBs) and twin domain boundaries (TDBs) promote irreversible α→ε interconversion through domain fragmentation. In β‐Ga₂O₃, defect‐induced stress gradients drive two distinct local transformations: surface Ga‐aggregated β→δ that stabilizes transient states, and latent‐track‐confined β→κ phase transition with recoverable distortions via cation reordering. Under electronic excitation, β‐Ga₂O₃ forms nanohillocks via robust GaO₆ octahedra (high density/strong Ga─O bonds), while α/ε‐Ga₂O₃ generates nanopores from tetrahedral Ga looseness (low bonding energy), highlighting phase‐dependent surface dynamics shaped by atomic packing and bonding anisotropy. Defect‐regulated recombination suppresses visible photoluminescence in α/β‐Ga₂O₃, whereas in ε‐Ga₂O₃ bandgap narrowing of ΔE: 0.30 eV is observed, enhancing emission. Linking phase‐dependent defect‐carrier interactions and metastable‐phase engineering in Ga₂O₃ enables property optimization for power‐electronics and optoelectronics devices.

electronic state configuration↗
DOE OSTI2023

Ab Initio Evaluation of Plutonium Dioxide $S(α,β)$ and Thermal Neutron Cross Sections

Plutonium Dioxide (PuO 2 ) is an advanced fuel-cycle nuclear fuel material for thermal and fast reactors, with a melting point at 3261.8 K and thermal conductivity above 3 Wm -1 K -1 for temperatures below 1500 K. Plutonium from spent UO 2 is reprocessed to form chemically stable mixed oxide (MOX) fuel. In contrast to metallic fuels, the low PuO 2 conductivity results in high thermal gradients between fuel center and surface, resulting in structural annealing and efficient fission product migration over burnup. Crystal binding affects the doppler broadening of epithermal resonances for nuclear fuels and scattering behavior of thermal neutrons. Currently, there is no ENDF/B Thermal Scattering Law, i.e., $S(α,β)$, evaluation for PuO 2 , only free atom cross sections for Plutonium and Oxygen, respectively, which neglect lattice contributions to total cross section. Herein, ab initio lattice dynamics (AILD) techniques are employed to calculate the phonon density of states (DOS) using spin-orbit-coupling density functional theory (DFT) to predict the PuO2 paramagnetic ground-state structure. These results highlight the evaluation of $S(α,β)$ for Plutonium and Oxygen, respectively in PuO 2 and consequential generation of thermal neutron scattering cross sections for high fidelity criticality safety analysis and reactor calculations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Extending Interfaces in 3D to Achieve Superior Nanoscale Strength in Ti/Nb Nanolaminates

Tuning the atomic-level structure of nanolaminates enables high strength, increased deformability, and the ability to absorb and mitigate damage due to varied crystalline defects, including dislocations. Here, we present the enhanced strength of Ti/Nb nanolaminates containing thick 3D interfaces (3DIs), relative to their chemically abrupt 2D counterparts. We examine the effects of crystallographic alignment and compositional gradients on mechanical behavior via experiments and phase-field-dislocation dynamics (PFDD) modeling. Mechanical testing reveals that nanolaminates containing thicker 3DIs demonstrate a 28% hardness enhancement compared with sharp-interface nanolaminates. PFDD modeling shows that the critical resolved shear stress (CRSS) increases with the 3DI thickness. Gradual compositional transitions in 3DIs were confirmed via scanning transmission electron microscopy and high-resolution transmission electron microscopy, showing sharp crystallographic transitions and a heightened interface topography. The findings establish a positive function between the 3DI thickness and mechanical robustness for hexagonal-closest-packed-containing composites, emphasizing the role of defect–interface interactions in tailoring the mechanical performance and providing a foundation for future interfacial engineering.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Nucleon-pair coupling scheme in Elliott's SU(3) model

Elliott's SU(3) model is at the basis of the shell-model description of rotational motion in atomic nuclei. Here we demonstrate that SU(3) symmetry can be realized in a truncated shell-model space if constructed in terms of a sufficient number of collective S, D, G,...pairs (i.e., with angular momentum zero, two, four,...) and if the structure of the pairs is optimally determined either by a conjugate-gradient minimization method or from a Hartree-Fock intrinsic state. We illustrate the procedure for six protons and six neutrons in the pf (sdg) shell and exactly reproduce the level energies and electric quadrupole properties of the ground-state rotational band with SDG (SDGI) pairs. The SD-pair approximation without significant renormalization, on the other hand, cannot describe the full SU(3) collectivity. A mapping from Elliott's fermionic SU(3) model to systems with s, d, g,... bosons provides insight into the existence of a decoupled collective subspace in terms of S, D, G,... pairs.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1976

Considerations to achieve directionality for gamma ray lasers

A method of alignment of nuclei for a gamma ray laser and a means of achieving preferential emission of radiation along the crystal axis are studied. Atomic alignment was achieved by materials researchers who made composite structures composed of needle-like single crystals all with the same orientation and all pointing in the same direction contained in a matrix of cobalt or nickel. The proposed method of preferential emission of radiation along the aligned needles is to have a symmetric field gradient at the nucleus and a sequence of excited levels of spin and parity 2(+) and 0(+). The proposed scheme reduces the density of excited states required for lasing and reduces the linewidth due to inhomogenous broadening. Mossbauer absorption experiments intended to test these ideas are outlined.

Jha, S.↗
NASA NTRS1981

Present growth technology of silicon germanium alloys and possible advantages of microgravity growth

The growth technology of Silicon-germanium (Si-Ge) alloys and the possible advantages of growth in microgravity is reviewed. The Si-Ge alloys have a continuous variation of bandgap energy from the germanium bandgap to the silicon bandgap. The unusual two slope behavior of Eg versus composition is due to the differences in the conduction band structure between Si and Ge. Below 17% (atomic), the germanium band structure dominates; and above it, the bands are "silicon like". It is found that the growth of Si-Ge alloys in microgravity is very attractive. In particular, the float zone method, in which a liquid zone of controlled starting composition, used to grow a large amount of useful alloy crystal. Large temperature gradients and relatively flat growth interfaces are necessary to obtain homogeneous crystal growth.

Stafsudd, O. M.↗
NASA NTRS2004

Pickup Ion Velocity Distributions at Titan: Effects of Spatial Gradients

The principle source of pickup ions at Titan is its neutral exosphere, extending well above the ionopause into the magnetosphere of Saturn or the solar wind, depending on the moon's orbital position. Thermal and nonthermal processes in the thermosphere generate the distribution of neutral atoms and molecules in the exosphere. The combination of these processes and the range of mass numbers, 1 to over 28, contribute to an exospheric source structure that produces pickup ions with gyroradii that are much larger or smaller than the corresponding scale heights of their neutral sources. The resulting phase space distributions are dependent on the spatial structure of the exosphere as well as that of the magnetic field and background plasma. When the pickup ion gyroradius is less than the source gas scale height, the pickup ion velocity distribution is characterized by a sharp cutoff near the maximum speed, which is twice that of the ambient plasma times the sine of the angle between the magnetic field and the flow velocity. This was the case for pickup H(sup +) ions identified during the Voyager 1 flyby. In contrast, as the gyroradius becomes much larger than the scale height, the peak of the velocity distribution in the source region recedes from the maximum speed. Iri addition, the amplitude of the distribution near the maximum speed decreases. These more beam like distributions of heavy ions were not observed from Voyager 1 , but should be observable by more sensitive instruments on future spacecraft, including Cassini. The finite gyroradius effects in the pickup ion velocity distributions are studied by including in the analysis the possible range of spatial structures in the neutral exosphere and background plasma.

Hartle, R. E.↗
DOE OSTI2025

Anomalous elastic softening in ferroelectric hafnia under pressure

his study employs first-principles density-functional theory (DFT) calculations to explore the elastic and mechanical properties of ferroelectric hafnia (HfO 2 ) in its polar orthorhombic 𝑃⁢𝑐⁢𝑎⁢2 1 phase under varying hydrostatic pressure conditions up to 30 GPa. Utilizing a plane-wave basis set and Perdew-Burke-Ernzerhof generalized-gradient approximation for solids in our DFT calculations, we investigate both pure and yttrium-substituted HfO 2 . Our findings reveal an anomalous reduction in the 𝐶 33 component of the elastic tensor with increasing pressure, which becomes significant above 15 GPa and signals a potential pressure-driven structural phase transition at higher pressure. The analysis of atomic displacements under pressure sheds light on the unusual mechanical behavior and phase stability of this material. Additionally, we observe a transition from an indirect band gap to a direct band gap with increasing pressure, which could have significant implications for optical applications. Here, the effects of yttrium substitution on the mechanical and electronic properties are further examined, revealing that yttrium substitution softens the elastic response of this material and reduces the electronic band gap. These results enhance our understanding of elastic and mechanical responses of ferroelectric hafnia and its potential for applications in microelectronics, piezoelectric devices, and nonvolatile ferroelectric random-access memories. Further experimental validation is recommended to confirm our predictions and explore the practical implications of the observed phase transitions and electronic behavior of the ferroelectric hafnia under high-pressure conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Microfibril orientation and compositional heterogeneity in fiber and vessel cell walls of poplar xylem studied by AFM-IR and SFG spectroscopy

Understanding the structural organization of cellulose microfibrils (CMFs) within individual plant cell walls is essential for connecting cell wall architecture to its mechanical and physiological functions. However, due to the complex hierarchical structure and nanoscale heterogeneity of cell walls, it remains technically challenging to resolve detailed compositional and orientational information at subcellular levels of individual cell walls. This study investigates the internal 3D structure, chemical composition, and sublayer organization of fiber and vessel cell walls in the xylem tissue of a two-year-old field-grown hybrid poplar tree (Populus alba × P. glandulosa) using photothermal atomic force microscopy coupled with infrared spectroscopy (AFM-IR) and sum frequency generation (SFG) hyperspectral microscopy. AFM-IR provided nanoscale chemical imaging, revealing localized compositional heterogeneity, including variations between adjacent cell walls and transitional layers beyond the traditional S1, S2, and S3 sublayers. SFG microscopy revealed that CMFs in fiber walls are highly aligned along the stem axis, consistent with their role in mechanical support, while vessel cell walls exhibited slightly tilted CMFs, reflecting their function in hydraulic transport. Together, these results offer new insights into cell-type-specific CMF organization and compositional gradients in hybrid poplar xylem. These findings highlight the structural and chemical complexity of secondary cell walls in woody plants and demonstrate the value of AFM-IR and SFG spectroscopy in elucidating plant cell wall architecture.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

High-throughput predictions of metal–organic framework electronic properties: theoretical challenges, graph neural networks, and data exploration

Abstract With the goal of accelerating the design and discovery of metal–organic frameworks (MOFs) for electronic, optoelectronic, and energy storage applications, we present a dataset of predicted electronic structure properties for thousands of MOFs carried out using multiple density functional approximations. Compared to more accurate hybrid functionals, we find that the widely used PBE generalized gradient approximation (GGA) functional severely underpredicts MOF band gaps in a largely systematic manner for semi-conductors and insulators without magnetic character. However, an even larger and less predictable disparity in the band gap prediction is present for MOFs with open-shell 3 d transition metal cations. With regards to partial atomic charges, we find that different density functional approximations predict similar charges overall, although hybrid functionals tend to shift electron density away from the metal centers and onto the ligand environments compared to the GGA point of reference. Much more significant differences in partial atomic charges are observed when comparing different charge partitioning schemes. We conclude by using the dataset of computed MOF properties to train machine-learning models that can rapidly predict MOF band gaps for all four density functional approximations considered in this work, paving the way for future high-throughput screening studies. To encourage exploration and reuse of the theoretical calculations presented in this work, the curated data is made publicly available via an interactive and user-friendly web application on the Materials Project.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fluid structural interaction simulation of the MP-1 plate performance irradiated in the advanced test reactor

Due to the high neutron fluxes they generate and increased heat transfer performance, plate type fuels are used in the U.S. high-performance research reactors. During irradiation, a significant amount of fission energy (i.e., ~200 MeV per atom) is released by the U-235 chain reaction that is carried away by the coolant. Through thermal hydraulic analysis, the coolant’s heat transfer capability is investigated to ensure that the desired fuel temperature can be maintained. At high temperatures, the plate undergoes elastic/plastic deformation, creep, and swelling as a result of both the temperature gradients and the fission gas production within the fuel. These effects are studied via fuel performance analysis. When the plate deformation is small enough that the coolant flow remains relatively unchanged, conducting these two types of analyses independently will suffice. But at high fuel burnups, the swollen plates may encroach into the coolant channels that separate the fuel plates from each other and cause these channels to narrow. Large reductions in channel gap size imperil cooling performance, causing fuel temperatures to rise. Further, if the plate deformation is asymmetric, the fuel centerline will shift toward one side of the channel, causing an uneven reduction in coolability. In addition, a boehmite (oxide) layer will, over time, grow on the plate surface, further obstructing the heat transfer process. To precisely predict fuel temperatures/deformation, a complete coupled analysis that considers coolant flow, heat transfer, oxide growth, elastic/plastic deformation, creep, and swelling is needed; however, this type of analysis method is not available in the literature. To fill the gap, this research developed a fluid structure interaction (FSI) approach to the fuel plate analysis, then successfully applied it to the Mini-Plate (MP)-1 experiment, which was irradiated in the Advanced Test Reactor (ATR) for both one and two cycles. The complete analysis coupled STAR-CCM+, a computational fluid dynamic (CFD) software for calculating flow, with Abaqus, a finite element analysis code for calculating plate deformation. Improvements in the results were found when comparing the fully coupled analysis to the independently conducted analyses but they were not significant due to the miniature size of the plates and the relatively short irradiation time. In the future, the fully coupled approach presented herein will be applied to full-size fuel plates with longer irradiation cycles once additional experiments become available.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Liquid-like solid-state diffusion of lithium ions in super-halide-rich argyrodite

The development of solid electrolytes with high ionic conductivity is essential for advancing safer, high-energy-density solid-state batteries, where lithium site distribution in the sublattice strongly affects ion transport. Here, we report a super-halide-rich argyrodite, Li 5.3 PS 4.3 Cl 1.7 , with remarkable room-temperature ionic conductivity (11.4 ± 0.7 mS cm -1 ) due to population of two additional interstitial lithium sites induced by vacancy redistribution. Prominent lithium density between lithium sites and elevated atomic displacement parameters indicate liquid-like diffusive behavior resembling sublattice melting. Combining electrochemical impedance spectroscopy, pulsed-field gradient NMR, and T 1 relaxation methods, we demonstrate that the augmented conductivity partly arises from a low energy barrier (0.08 eV) at the local scale, attributed to a three-site lithium distribution that drives correlated lithium dynamics. This work advances our understanding of the structure-dynamics interplay in super-halide-rich argyrodites, and highlighting their potential as solid-state battery electrolytes in cells with a coated single-crystal NMC82 cathode that achieve 170 mAh/g capacity at a 0.2 C rate .

25 ENERGY STORAGE↗
DOE OSTI2025

TaTi_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-titanium (Ta-Ti). 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. 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 Hf and Nb. 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 Air Force HPC11 cluster using the VASP 6.5.1. Additional methodology and file structure information is available in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2025

HfNb_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 hafnium-niobium (Hf-Nb). 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. 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 Hf and Nb. 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,029 randomized atomic structures over 31 chemical compositions. Further methodological and structural information is contained in the dataset README.txt file.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analytical nonadiabatic coupling and state-specific energy gradient for the crystal field Hamiltonian describing lanthanide single-ion magnets

Paramagnetic molecules with a metal ion as an electron spin center are promising building blocks for molecular qubits and high-density memory arrays. However, fast spin relaxation and decoherence in these molecules lead to a rapid loss of magnetization and quantum information. Nonadiabatic coupling (NAC), closely related to spin-vibrational coupling, is the main source of spin relaxation and decoherence in paramagnetic molecules at higher temperatures. Predicting these couplings using numerical differentiation requires a large number of computationally intensive ab initio or crystal field electronic structure calculations. To reduce computational cost and improve accuracy, we derive and implement analytical NAC and state-specific energy gradient for the ab initio parametrized crystal field Hamiltonian describing single-ion molecular magnets. Our implementation requires only a single crystal field calculation. In addition, the accurate NACs and state-specific energy gradients can be used to model spin relaxation using sophisticated nonadiabatic molecular dynamics, which avoids the harmonic approximation for molecular vibrations. To test our implementation, we calculate the NAC values for three lanthanide complexes. Finally, the predicted values support the relaxation mechanisms reported in previous studies.

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

NbV_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-vanadium (Nb-V). 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 V. 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. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and 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. and 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. and 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. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2024

NbZr_BCC_SolidSolution_128atoms_VASP6

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

36 MATERIALS SCIENCE↗