Search NASA⌕ Search

SEARCH · Search NASA

Results for “electronic wave function”

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

Connectivity-Dependent Exciton–Phonon Coupling in Cesium Bismuth Halide Quantum Dots

Metal halide octahedra form the fundamental functional building blocks of metal halide perovskites, dictating their structures, optical properties, electronic structures, and dynamics. Here, in this study, we show that the connectivity of bismuth halide octahedra in Cs 3 Bi 2 Br 9 and Cs 3 Bi 2 I 9 quantum dots (QDs) changes with different halide elements. We use first-principles calculations to reveal the key role of the connectivity of bismuth halide octahedra on the wave function symmetry, Huang-Rhys factor, and exciton-phonon interaction strength. Following QD synthesis via a ligand-mediated transport method, the effect of connectivity is verified with transient absorption spectroscopy, where we contrast Cs 3 Bi 2 Br 9 and Cs 3 Bi 2 I 9 QD exciton dynamics. In photoexcited Cs 3 Bi 2 I 9 QDs, phonons related to the vibrational motions of face-sharing [BiI 6 ] 3- bioctahedra couple strongly to the electronic state and drive rapid carrier relaxation. Equivalent signals are not observed for photoexcited Cs 3 Bi 2 Br 9 QDs, implying a lack of phonon involvement in band-edge absorption and subsequent exciton relaxation. Our findings suggest that structural engineering can effectively tune the exciton-phonon coupling and therefore influence exciton relaxation and recombination in perovskite nanomaterials.

TDDFT↗
DOE OSTI2020

Electron beam injection from a hollow cathode plasma into a downstream reactive environment: Characterization of secondary plasma production and Si 3 N 4 and Si etching

A material etching system was developed by combining beam electron injection from a direct current (DC) hollow-cathode (HC) electron source with the downstream reactive environment of a remote CF 4 /O 2 low temperature plasma. The energy of the injected beam electrons is controlled using an acceleration electrode biased positively relative to the HC Argon discharge. For an acceleration voltage greater than the ionization potential of Ar, the extracted primary electrons can produce a secondary plasma in the process chamber. We characterized the properties of the secondary plasma by performing Langmuir probe measurements of the electron energy probability function (EEPF) 2.5 cm below the extraction ring. The data indicate the existence of two major groups of electrons, including electrons with a primary beam electron energy which varies as the acceleration voltage is varied, along with low energy electrons produced by ionization of the Ar gas atoms in the process chamber by the injected beam electrons. When combining the HC Ar beam electron with a remote CF 4 /O 2 electron cyclotron wave resonance (ECWR) plasma, the EEPF of both the low energy plasma electron and beam electron components decreases. Additionally, we studied surface etching of Si3N4 and poly-crystalline Si (poly-Si) thin films as a function of process parameters, including the acceleration voltage (0-70 V), discharge current of the HC discharge (1-2 A), pressure (2- 100 mTorr), source to substrate distance (2.5 to 5 cm), and feed gas composition (with or without CF 4 /O 2 ). The direction of the incident beam electrons was perpendicular to the surface. Si 3 N 4 and poly-crystalline silicon etching are seen and indicate an electronneutral synergy effect. No/little remote plasma spontaneous etching was observed for the conditions used in this study, and the etching is confined to the substrate area irradiated by the injected beam electrons. The electron etched Si 3 N 4 surface etching rate profile distribution is confined within a ~30 mm diameter circle, which is slightly broader than the area for which poly-Si etching is seen, and coincides closely with the spatial profile of beam electrons as determined by the Langmuir probe measurements. The magnitude of poly-Si etching rate is by a factor of two times smaller than the Si 3 N 4 etching rate. In this paper, we discuss possible explanations of the data and the role of surface charging.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Efficient solutions of fermionic systems using artificial neural networks

In this study, we explore the similarities and differences between variational Monte Carlo techniques that employ conventional and artificial neural network representations of the ground-state wave function for fermionic systems. Our primary focus is on shallow neural network architectures, specifically the restricted Boltzmann machine, and we examine unsupervised learning algorithms that are appropriate for modeling complex many-body correlations. We assess the advantages and drawbacks of conventional and neural network wave functions by applying them to a range of circular quantum dot systems. Our findings, which include results for systems containing up to 90 electrons, emphasize the efficient implementation of these methods on both homogeneous and heterogeneous high-performance computing facilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Disentangling the physics of the attractive Hubbard model as a fully interacting model of fermions via the accessible and symmetry-resolved entanglement entropies

The complicated ways in which electrons interact in many-body systems such as molecules and materials have long been viewed through the lens of local electron correlation and associated correlation functions. However, quantum information science has demonstrated that more global diagnostics of quantum states like the entanglement entropy can provide a complementary and clarifying lens on electronic behavior. One particularly useful measure that can be used to distinguish between quantum and classical sources of entanglement is the accessible entanglement, the entanglement available as a quantum resource for systems subject to conservation laws, such as fixed particle number, due to superselection rules. In this work, we introduce an algorithm and demonstrate how to compute accessible and symmetry-resolved entanglements for interacting fermion systems. This is accomplished by combining an incremental version of the swap algorithm with a recursive auxiliary field quantum Monte Carlo algorithm recently developed by the authors. We apply these tools to study the pairing and charge density waves exhibited in the paradigmatic attractive Hubbard model via entanglement. We find that the particle and spin symmetry-resolved entanglements and their related full probability distribution functions show very clear—and unique—signatures of the underlying electronic behavior even when those features are less pronounced in conventional correlation functions. Altogether, this work provides a systematic means of characterizing the entanglement within quantum systems that can grant a deeper understanding of the complicated electronic behavior that underlies quantum phase transitions and crossovers in many-body systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

RF antenna helicity dependent particle heating in a helicon source

Experiments have demonstrated that ion phenomena, such as the lower hybrid resonance, play an important role in helicon source operation. Damping of the slow branch of the bounded whistler wave at the edge of a helicon source (i.e. the Trivelpiece-Gould mode) has been correlated with the creation of energetic electrons, heating of ions at the plasma edge, and anisotropic ion heating. Here we present ion velocity distribution function measurements, electron density and temperature measurements, and magnetic fluctuation measurements on both sides of an m = |1| helical antenna in a helicon source as a function of the driving frequency, magnetic field strength, and magnetic field orientation relative to the antenna helicity. Significant electron and ion heating (up to two times larger) occurs on the side of the antenna consistent with the launch of the m = +1 mode. The electron and ion heating occurs within one electron skin depth of the plasma edge, where slow wave damping is expected. The source parameters for enhanced particle heating are also consistent with lower hybrid resonance effects, which can only occur for Trivelpiece-Gould wave excitation.

wave modes↗
DOE OSTI2021

Orbital contributions in the element-resolved valence electronic structure of Bi 2 Se 3

In this work, we studied the bulk band structure of a topological insulator (TI) Bi 2 Se 3 and determined the contributions of the Bi and Se orbital states to the valence bands using standing-wave excited hard x-ray photoemission spectroscopy (SW-HAXPES). This SW technique can provide the element-resolved information and extract individual Bi and Se contributions to the Bi2Se3 valence band. Comparisons with density-functional theory calculations (local density approximation and GW) reveal that the Bi 6s, Bi 6p, and Se 4p states are dominant in the Bi 2 Se 3 HAXPES valence band. Furthermore, these findings pave a way for studying the element-resolved band structure and orbital contributions of this class of TIs.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Copolymers of vinylferrocene and vinylimidazole for potential sensing applications

Abstract Since imidazole is easily protonated and well known to coordinated to various transition metals, its incorporation into a polymer with ferrocene may result in interesting electronic interactions between the ferrocenyl and imidazolium moieties. To this end, copolymers of vinylferrocene and vinylimidazole were prepared in various ratios by free radical polymerization. The interactions between the ferrocenyl and imidazolium moieties were studied by cyclic voltammetry, UV–Vis spectroscopy and density functional theory (DFT) studies. Cyclic voltammetry of polymer thin films showed multiple oxidation waves at low pH or in the presence of Pb 2+ or Cd 2+ ions. The UV–Vis spectrum of thin films showed an additional absorbance at 625 nm suggesting electronic interactions. In addition, DFT studies indicated that the morphology of the polymer would place some of the ferrocenyl and imidazolium moieties in close proximity.

Polymer Science↗
DOE OSTI2023

Variable magnetic field electron spectrometer to measure hot electrons in the range of 50–460 keV

Resonance absorption (RA) occurs when a p-polarized electromagnetic wave, obliquely incident on an inhomogeneous plasma, tunnels past its turning point and resonantly excites an electron plasma wave (EPW) at the critical density. This phenomenon is important, for instance, in the direct drive approach to inertial fusion energy and is a particular example of a wider phenomenon in plasma physics known as mode conversion, which is crucial for heating magnetic fusion devices, such as tokamaks, via RF heating. Direct measurement of these RA-generated EPW accelerated hot electrons, with energy in the range of a few tens to a few hundreds of keV, is a challenging task due to the relatively low deflecting magnetic fields needed. The solution described here is a magnetic electron spectrometer (MES) with a continually changing magnetic field, lower at the entrance of the MES and gradually increasing toward the end, that enables the measurement of a wide spectral range of electrons with energies between 50 and 460 keV. Electron spectra taken in a LaserNetUS RA experiment were acquired from plasmas generated by irradiating polymer targets with the combination of an ~300 ps pulse followed by a series of ten high intensity 50–200 fs duration laser pulses from the ALEPH laser at Colorado State University. In conclusion, the high intensity beam is designed as spike trains of uneven duration and delay pulses in order to modify the RA phenomenon.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

237 Np Mössbauer Isomer Shifts: A Lesson About the Balance of Static and Dynamic Electron Correlation in Heavy Element Complexes

A large set of neptunium compounds with different oxidation states (III to VII) was assembled to study the Mössbauer isomer shift by wave function calculations and better understand covalency in f-elements complexes. The contact density approach was used to calculate the isomer shift using complete active space self-consistent field (CASSCF) multiconfiguration wave functions, as well as matrix product states [from Density Matrix Renormalization Group (DMRG) algorithms] for large active spaces. Dynamic correlation effects for the isomer shifts were treated via CASPT2 energy derivatives with respect to the nuclear radius. The CASSCF calculations appear to produce different orbital overlocalization errors for low and high Np oxidation states. For compounds with low Np oxidation numbers, the errors can be attributed to the overlocalization of the 5f orbitals. For the compounds with high Np oxidation numbers, the main errors arise from an overlocalization of ligand orbitals and concomitant to weak donation bonding. Attempts to mitigate the overlocalization errors with large active spaces using DMRG were only partially successful, showing that explicit treatment of dynamic correlation is necessary for accurate predictions of Mössbauer isomer shifts. The CASPT2 calculations perform very satisfactorily. For a subset of Np compounds, both static and dynamic correlation effects were substantial. Furthermore, a rational active space selection based on orbital entanglement diagrams proved beneficial for determining the optimal reference wave function.

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

Variational Discrete Action Theory

In this work, we propose the variational discrete action theory (VDAT) to study the ground state properties of quantum many-body Hamiltonians. VDAT is a variational theory based on the sequential product density matrix (SPD) ansatz, characterized by an integer $\mathscr{N}$, which monotonically approaches the exact solution with increasing $\mathscr{N}$. To evaluate the SPD, we introduce a discrete action and a corresponding integer time Green’s function. We use VDAT to exactly evaluate the SPD in two canonical models of interacting electrons: the Anderson impurity model and the d = ∞ Hubbard model. For the latter, we evaluate $\mathscr{N}$ = 2 – 4, where $\mathscr{N}$ = 2 recovers the Gutzwiller approximation (GA), and we show that $\mathscr{N}$ = 3, which exactly evaluates the Gutzwiller-Baeriswyl wave function, provides a truly minimal yet precise description of Mott physics with a cost similar to that of the GA. VDAT is a flexible theory for studying quantum Hamiltonians, competing both with state-of-the-art methods and simple, efficient approaches all within a single framework.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

The ABCs of phase retrieval: Connecting the acronyms of scanning transmission electron microscopy

High-resolution scanning transmission electron microscopy (S/TEM) is an indispensable tool for characterizing the structure and properties of materials down to the atomic scale. Conventional S/TEM imaging, however, is limited by the phase problem, whereby the phase of the electron exit wave is lost upon detection. Recent advances in diffractive imaging and 4D-STEM have enabled a range of phase-retrieval techniques that computationally reconstruct the missing information encoded in the phase of the transmission function. These approaches offer improved dose efficiency and enhanced sensitivity to weakly scattering signals, extending quantitative imaging to beam-sensitive materials composed of light elements. In this work, we introduce the phase problem in electron microscopy and survey the diverse landscape of phase-retrieval techniques used in the field. Despite their many acronyms and algorithmic variations, these techniques share a common physical and mathematical foundation. We present a unified framework that connects these seemingly distinct methods, from parallax imaging and tilt-corrected bright-field (tcBF-STEM), to aberration-corrected bright-field (acBF-STEM), optimum bright-field (OBF-STEM) and single-sideband (SSB) ptychography, as well as first-moment integrated center of mass techniques (iCOM) and iterative ptychographic algorithms. Based on these insights, we discuss the opportunities and practical limitations of applying these methods across different materials systems, detector designs, and microscope configurations.Graphical abstractRepresentative electron microscopy configurations used for phase retrieval and diffractive imaging in S/TEM: (a) Zernike phase-contrast transmission electron microscopy (TEM), (b) small-convergence-angle four-dimensional scanning transmission electron microscopy (4D-STEM) for nanobeam-based phase reconstruction methods, and (c) large-convergence-angle 4D-STEM for ptychographic and related diffractive imaging techniques reviewed in this work.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multinode Multi-GPU Two-Electron Integrals: Code Generation Using the Regent Language

The computation of two-electron repulsion integrals (ERIs) is often the most expensive step of integral-direct self-consistent field methods. Formally it scales as O(N 4 ), where N is the number of Gaussian basis functions used to represent the molecular wave function. In practice, this scaling can be reduced to O(N 2 ) or less by neglecting small integrals with screening methods. The contributions of the ERIs to the Fock matrix are of Coulomb (J) and exchange (K) type and require separate algorithms to compute matrix elements efficiently. We previously implemented highly efficient GPU-accelerated J-matrix and K-matrix algorithms in the electronic structure code TeraChem. Although these implementations supported the use of multiple GPUs on a node, they did not support the use of multiple nodes. This presents a key bottleneck to cutting-edge ab initio simulations of large systems, e.g., excited state dynamics of photoactive proteins. We present our implementation of multinode multi-GPU J- and K-matrix algorithms in TeraChem using the Regent programming language. Regent directly supports distributed computation in a task-based model and can generate code for a variety of architectures, including NVIDIA GPUs. We demonstrate multinode scaling up to 45 GPUs (3 nodes) and benchmark against hand-coded TeraChem integral code. Finally, we also outline our metaprogrammed Regent implementation, which enables flexible code generation for integrals of different angular momenta.

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

VZr_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 vanadium-zirconium (V-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 V 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.

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

Marginal stability constraint on runaway electron distribution

High-frequency kinetic instabilities of the strongly anisotropic runaway electrons (RE) can enhance the pitch-angle scattering of the runaways significantly. This wave-induced scattering can easily prevail over runaway scattering on high-Z impurities. In a steady state, collisional damping balances the kinetic drive of the unstable waves, keeping the RE distribution function at marginal stability. The marginal stability constraint limits the achievable RE densities and the shape of the RE distribution function. In this study, we consider whistler and compressional Alfvén waves as the primary source of enhanced elastic scattering of the runaways. By balancing the anomalous Doppler resonance drive with the collisional wave damping, we find the RE distribution function in the ultra-relativistic range of the phase space. We also derive an expression for the spectral energy density of the waves. We show that the power needed to compensate for the wave dissipation is negligible compared to the work of the electric field. The latter is in balance with the synchrotron losses of the runaway electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Quasi-particle interference of the van Hove singularity in Sr2RuO4

Abstract The single-layered ruthenate Sr 2 RuO 4 is one of the most enigmatic unconventional superconductors. While for many years it was thought to be the best candidate for a chiral p -wave superconducting ground state, desirable for topological quantum computations, recent experiments suggest a singlet state, ruling out the original p -wave scenario. The superconductivity as well as the properties of the multi-layered compounds of the ruthenate perovskites are strongly influenced by a van Hove singularity in proximity of the Fermi energy. Tiny structural distortions move the van Hove singularity across the Fermi energy with dramatic consequences for the physical properties. Here, we determine the electronic structure of the van Hove singularity in the surface layer of Sr 2 RuO 4 by quasi-particle interference imaging. We trace its dispersion and demonstrate from a model calculation accounting for the full vacuum overlap of the wave functions that its detection is facilitated through the octahedral rotations in the surface layer.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Calculation of spin–orbit couplings using RASCI spinless one-particle density matrices: Theory and applications

Here we present the formalism and implementation for calculations of spin–orbit couplings (SOCs) using the Breit–Pauli Hamiltonian and non-relativistic wave functions described by the restricted active space configuration interaction (RASCI) method with general excitation operators of spin-conserving spin-flipping, ionizing, and electron-attaching types. The implementation is based on the application of the Wigner–Eckart theorem within the spin space, which enables the calculation of the entire SOC matrix based on the explicit calculation of just one transition between the two spin multiplets. Numeric results for a diverse set of atoms and molecules highlight the importance of a balanced treatment of correlation and adequate basis sets and illustrate the overall robust performance of RASCI SOCs. The new implementation is a useful addition to the methodological toolkit for studying spin-forbidden processes and molecular magnetism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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↗