SEARCH · Search NASA
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.
Here, we use pair distribution function analysis of synchrotron X-ray total scattering data to inspect the local structure of BaTiO 3 nanocrystals from 253 K < T < 413 K. The local structure consists of ca. r = 35 Å domains in which Ti 4+ displacements have the same coherence length as the overall perovskite structure.
We present a multiwavelength data analysis of IRAS 05463+2652 (hereafter I05463+2652) to study star formation mechanisms. A shell-like structure around I05463+2652 is evident in the Herschel column density map, which is not associated with any ionized emission. Based on the Herschel submillimeter images, several parsec-scale filaments (including two elongated filaments, “s-fl” and “nw-fl” having lengths of ∼6.4 and ∼8.8 pc, respectively) are investigated in the I05463+2652 site. The Herschel temperature map depicts all these features in a temperature range of ∼11–13 K. 39 clumps are identified and have masses between ∼70--945 M{sub ⊙}. The majority of clumps (having M{sub clump}≳300 M{sub ⊙}) are distributed toward the shell-like structure. 175 young stellar objects (YSOs) are selected using the photometric 1–5 μm data and a majority of these YSOs are distributed toward the four areas of high column density (≳5×10{sup 21} cm{sup −2}; A {sub V} ∼ 5.3 mag) in the shell-like structure, where massive clumps and a spatial association with filament(s) are also observed. The knowledge of observed masses per unit length of elongated filaments and critical mass length reveals that they are supercritical. The filament “nw-fl” is fragmented into five clumps (having M{sub clump}∼100--545 M{sub ⊙}) and contains noticeable YSOs, while the other filament “s-fl” is fragmented into two clumps (having M{sub clump}∼170--215 M{sub ⊙}) without YSOs. Together, these observational results favor the role of filaments in the star formation process in I05480+2545. This study also reveals the filament “s-fl,” containing two starless clumps, at an early stage of fragmentation.
Graph Neural Networks (GNN) are indispensable in learning from graph-structured data, yet their rising computational costs, especially on massively connected graphs, pose significant challenges in terms of execution performance. To tackle this, distributed-memory solutions such as partitioning the graph to concurrently train multiple replicas of GNNs are in practice. However, approaches requiring a partitioned graph usually suffer from communication overhead and load imbalance, even under optimal partitioning and communication strategies due to irregularities in the neighborhood minibatch sampling. This paper proposes practical trade-offs for improving the sampling and communication overheads for representation learn- ing on distributed graphs (using popular GraphSAGE architecture) by developing a parameterized prefetch and eviction scheme on top of the state-of-the-art Amazon DistDGL distributed GNN framework, demonstrating about 15–40% improvement in end-to-end training performance on the NERSC Perlmutter supercomputer for various OGB datasets.
Geophysical imaging techniques are a non-invasive way to image the subsurface and understand both subsurface solid (rock/soil) and fluid property distributions and their evolution in time. Inversions of the geophysical data, such as Electrical Resistance Tomography (ERT) data, are solved to estimate the subsurface property distributions, such as conductivity, and many inversion techniques smooth out sharp gradients in rock or fluid property distributions. Sharp gradients in subsurface properties tend to be present in situations with complex subsurface structures, which are common in many subsurface applications. We have successfully demonstrated that it is possible to inform, or constrain, inversions with neural networks trained on synthetic data with complex subsurface structures. Initial results suggest this process may be optimizable to yield property distributions that better represent the true property distributions than the same inversion process without the neural network constraint. Future work would optimize the neural network performance for this application and then apply the synthetic-data trained neural network to real data to understand the utility and performance of this technique for real data sets.
Jurisdictions are increasingly adopting compensation structures for distributed PV and PV-battery systems that price exported energy lower than energy consumed onsite (net billing rates). Standard methods for calculating the bill savings from PV and PV-battery systems were developed for net metering structures, and applying these same methods to net billing structures (such as using Typical Meteorological Year weather with actual year load) introduces bias errors that underestimate PV-battery system bill savings by between 1.5\% and 9\%, depending on the utility rate. We assess the magnitude of these errors and compare them to other sources of uncertainty when estimating the bill savings from PV and PV-battery systems under more complex utility rates.
Topochemical reactions involving ionic exchange have been used to assess a large number of metastable compositions, particularly in layered metal oxides. This method encompasses complex reactions that are poorly explored, yet are of prime importance to understand and control the materials’ properties. In this work, we embark on investigating the reactions involved during the ionic exchange between a layered Na-titanate (lepidocrocite-type structure) and an acidic solution (HCl), leading to a protonic (H 3 O + ) titanate (trititanate structure). The reactions involve an ionic exchange provoking a structural change from the lepidocrocite-type to the trititanate structure as shown by real-space refinements of ex situ pair distribution function data. Mobile Na + ions are exchanged by hydronium ions inducing high proton mobility in the final structure. Moreover, the reaction was followed by ex situ 23 Na and 1 H solid-state MAS NMR which allowed, among other things, confirming that the Na + ions are in the interlayer space and specifying their local environment. Strikingly, the ionic exchange reaction induces progressive exfoliation of the Na-titanate particles leading to 2–5 nm thin elongated crystallites. To further understand the different steps associated with the ionic exchange, the evolution of the electrolytic conductivity, using conductimetric titration, has been monitored upon HCl addition, enabling characterization of the intercalation(H + )/de-intercalation(Na + ) reactions and assessing kinetic parameters. Accordingly, it is hypothesized that the exfoliation of the particles is due to the accumulation of charges at the particle level in relation to the rapid intercalation of protons. This work provides novel insights into ionic exchange reactions involved in layered oxide compounds.
Abstract Large‐scale estimations of flood losses are often based on spatially aggregated inputs. This makes risk assessments vulnerable to aggregation bias, a well‐studied, sometimes substantial outcome in analyses that model fine‐grained spatial phenomena at coarse spatial units. To evaluate this potential in the context of large‐scale flood risk assessments, we use data from a high‐resolution flood hazard model and structure inventory for over 1.3 million properties in Massachusetts and examine how prominent data aggregation approaches affect the magnitude and spatial distribution of flood loss estimates. All considered aggregation approaches rely on aggregate structure inventories but differ in whether flood hazard is also aggregated. We find that aggregating only structure inventories slightly underestimates overall losses (−10% bias), and when flood hazard data is spatially aggregated to even relatively small spatial units (census block), statewide aggregation bias can reach +366%. All aggregation‐based procedures fail to capture the spatial covariation of inputs distributions in the upper tails that disproportionately generate total expected losses. Our findings are robust to several key assumptions, add important context to published risk assessments and highlight opportunities to improve flood loss estimation uncertainty quantification.
Electron counting helped realize the resolution revolution in single-particle cryoEM and is now accelerating the determination of MicroED structures. Its advantages are best demonstrated by new direct electron detectors capable of fast (kilohertz) event-based electron counting (EBEC). This strategy minimizes the inaccuracies introduced by coincidence loss (CL) and promises rapid determination of accurate structures. We used the Direct Electron Apollo camera to leverage EBEC technology for MicroED data collection. Given its ability to count single electrons, the Apollo collects high-quality MicroED data from organic small-molecule crystals illuminated with incident electron beam flux densities as low as 0.01–0.045 e − /Å 2 /s. Under even the lowest flux density (0.01 e − /Å 2 /s) condition, fast EBEC data produced ab initio structures of a salen ligand (268 Da) and biotin (244 Da). Each structure was determined from a 100° wedge of data collected from a single crystal in as few as 50 s, with a delivered fluence of only ∼0.5 e − /Å 2 . Fast EBEC data collected with a fluence of 2.25 or 3.33 e − /Å 2 also facilitated a 1.5 Å structure of thiostrepton (1665 Da). While refinement of these structures appeared unaffected by CL, a CL adjustment applied to EBEC data further improved the distribution of intensities measured from the salen ligand and biotin crystals. However, CL adjustment only marginally improved the refinement of their corresponding structures, signaling the already high counting accuracy of detectors with counting rates in the kilohertz range. Overall, by delivering low-dose structure-worthy data, fast EBEC collection strategies open new possibilities for high-throughput MicroED.
First-principles fusion plasma simulations are both compute and memory intensive, and CGYRO is no exception. The use of many HPC nodes to fit the problem in the available memory thus results in significant communication overhead, which is hard to avoid for any single simulation. That said, most fusion studies are composed of ensembles of simulations, so we developed a new tool, named XGYRO, that executes a whole ensemble of CGYRO simulations as a single HPC job. By treating the ensemble as a unit, XGYRO can alter the global buffer distribution logic and apply optimizations that are not feasible on any single simulation, but only on the ensemble as a whole. The main saving comes from the sharing of the collisional constant tensor structure, since its values are typically identical between parameter-sweep simulations. This data structure dominates the memory consumption of CGYRO simulations, so distributing it among the whole ensemble results in drastic memory savings for each simulation, which in turn results in overall lower communication overhead.
We present DMNetwork, a high-level package included in the PETSc library for the simulation of multiphysics phenomena over large-scale networked systems. The library aims at applications that have networked structures such as those in electrical, gas, and water distribution systems. DMNetwork provides data and topology management, parallelization for multiphysics systems over a network, and hierarchical and composable solvers to exploit the problem structure. DMNetwork eases the simulation development cycle by providing the necessary infrastructure through simple abstractions to define and query the network components. This article presents the design of DMNetwork, illustrates its user interface, and demonstrates its ability to solve multiphysics systems, such as an electric circuit, a network of power grid and water subnetworks, and transient hydraulic systems over large networks with more than 2 billion variables on extreme-scale computers using up to 30,000 processors.
Abstract One of the fundamental goals of particle physics is to gain a microscopic understanding of the strong interaction. Electromagnetic form factors quantify the structure of hadrons in terms of charge and magnetization distributions. While the nucleon structure has been investigated extensively, data on hyperons are still scarce. It has recently been demonstrated that electron-positron annihilations into hyperon-antihyperon pairs provide a powerful tool to investigate their inner structure. We present a method useful for hyperon-antihyperon pairs of different types which exploits the cross section enhancement due to the effect of vacuum polarization at theJ/ψresonance. Using the 10 billionJ/ψevents collected with the BESIII detector, this allows a precise determination of the hyperon structure function. The result is essentially a precise snapshot of the$$\bar{\Lambda }{\Sigma }^{0}\,(\Lambda {\bar{\Sigma }}^{0})$$ Λ ¯ Σ 0 ( Λ Σ ¯ 0 ) transition process, encoded in the transition form factor ratio and phase. Their values are measured to beR = 0.860 ± 0.029(stat.) ± 0.015(syst.),$$\Delta {\Phi }_{\bar{\Lambda }{\Sigma }^{0}}=(1.011\pm 0.094({{\rm{stat.}}})\pm 0.010({{\rm{syst.}}}))\,{{\rm{r}}}ad$$ Δ Φ Λ ¯ Σ 0 = ( 1.011 ± 0.094 ( stat. ) ± 0.010 ( syst. ) ) r a d and$$\Delta {\Phi }_{\Lambda {\bar{\Sigma }}^{0}}=(2.128\pm 0.094({{\rm{stat.}}})\pm 0.010({{\rm{syst.}}}))\,{{\rm{r}}}ad$$ Δ Φ Λ Σ ¯ 0 = ( 2.128 ± 0.094 ( stat. ) ± 0.010 ( syst. ) ) r a d . Furthermore, charge-parity (CP) breaking is investigated in this reaction and found to be consistent with CP symmetry.
Tensor algebra is a powerful tool for computing, but writing optimized codes that operate on sparse tensors can be very complex. This project enables a Tensor Algebra Compiler (TACO) that simplifies this task from man-years to man-days and extends TACO to support complex and large distributed systems. This report details the hypotheses, approaches used, and findings in this project.
We demonstrate that the Discrete Persistent Source Extractor (DisPerSE) can be used with spectroscopic redshifts to define the cosmic web and its distance to galaxies in small-area deep fields. Here we analyze the use of DisPerSE to identify structure in observational data. We apply DisPerSE to the distribution of galaxies in the Cosmic Evolution Survey (COSMOS) field and find the best parameters to identify filaments. We compile a catalog of 11,500 spectroscopic redshifts from the Galaxy and Mass Assembly (GAMA) G10 data release. We analyze two-dimensional slices, extract filaments, and calculate the distance for each galaxy to its nearest filament. We find that redder and more massive galaxies are closer to filaments. To study the growth of galaxies across cosmic time, and environment, we are carrying out an H i survey covering redshifts of z = 0–0.45, the COSMOS H i Large Extragalactic Survey (CHILES). In addition we present the predicted H i mass fraction as a function of distance to filaments for the spectroscopically known galaxies in CHILES. Lastly we discuss the cold gas morphology of a few individual galaxies and their positions with respect to the cosmic web. The identification of the cosmic web, and the ability of CHILES to study the resolved neutral hydrogen morphologies and kinematics of galaxies, will allow future studies of the properties of neutral hydrogen in different cosmic web environments across the redshift range of z = 0.1–0.45.
Hyperstoichiometric UO 2.15 was characterized by neutron total scattering at high temperature in the single-phase UO 2+x region of the U/O phase diagram. The diffraction data confirmed a single-phase fluorite structure at high temperature. Analysis of the short-range data showed that the same structural model does not fit the pair distribution functions well. Instead, structural models containing specific configurations of oxygen defect clusters best represent the local atomic arrangement. In conclusion, prevalent defect clusters previously proposed were fit to the experimental data, and moderately distorted oxygen cuboctahedra hypothesized by recent molecular dynamics simulations fit the data most accurately.
The prospect of creating ferroelectric or high permittivity nanomaterials provides motivation for investigating complex transition metal oxides of the form Ba(Ti, MV)O3, where M = Nb or Ta. Solid state processing typically produces mixtures of crystalline phases, rarely beyond minimally doped Nb/Ta. Using a modified sol-gel method, we prepared single phase nanocrystals of Ba(Ti, M)O3. Compositional and elemental analysis puts the empirical formulas close to BaTi0.5Nb0.5O3−δ and BaTi0.5Ta0.5O3−δ. For both materials, a reversible temperature dependent phase transition (non-centrosymmetric to symmetric) is observed in the Raman spectrum in the region 533–583 K (260–310 °C); for Ba(Ti, Nb)O3, the onset is at 543 K (270 °C); and for Ba(Ti, Ta)O3, the onset is at 533 K (260 °C), which are comparable with 390–393 K (117–120 °C) for bulk BaTiO3. The crystal structure was resolved by examination of the powder x-ray diffraction and atomic pair distribution function (PDF) analysis of synchrotron total scattering data. It was postulated whether the structure adopted at the nanoscale was single or double perovskite. Double perovskites (A2B′B″O6) are characterized by the type and extent of cation ordering, which gives rise to higher symmetry crystal structures. PDF analysis was used to examine all likely candidate structures and to look for evidence of higher symmetry. The feasible phase space that evolves includes the ordered double perovskite structure Ba2(Ti, MV)O6 (M = Nb, Ta) Fm-3m, a disordered cubic structure, as a suitable high temperature analog, Ba(Ti, MV)O3Pm-3m, and an orthorhombic Ba(Ti, MV)O3Amm2, a room temperature structure that presents an unusually high level of lattice displacement, possibly due to octahedral tilting, and indication of a highly polarized crystal.
In the last decade, significant progress has been made in the many-body Hamiltonian formulation of liquid dynamics theory. Earlier analysis of experimental data for a wide variety of elemental liquids had provided the reliable but qualitative description of the atomic motion: vibrations in a representative 3N -dimensional potential energy valley, plus transits, in which atoms cross the intersections between these valleys. Recent comparisons of first-principles theory with experiment for several elemental liquids at melt revealed a highly accurate and versatile theory. In the present work, we report on an extensive quench study of the entire condensed-matter structure-energy distribution for a metallic Na MD system. With these results, all that was learned from experimental data is confirmed, refined in detail, and made more accurate. Here, we show the entire structure-energy distribution, composed of widespread symmetrics and higher lying randoms, we show the increasing dominance of the randoms as N increases, until the symmetrics vanish completely, and we show the random distribution continue its spectacular narrowing as N continues to increase. This behavior certifies our early assignment of the random distribution to the liquid phase and our prediction of macroscopic uniformity of the random structures. Procedures are discussed to identify and calibrate a single random structure to represent the liquid, and the role of the structure energy in liquid thermodynamics is described. A comparison with other liquid dynamics theories is observed in the Introduction, and the relation to the equation-of-state program is noted in the Conclusion.
This report summarizes the work performed under the project "Advanced Data Structures for Improved Cyber Resilience and Awareness in Untrusted Environments." The goal of the project was to design, analyze, and test new data structures for cybersecurity applications. We had two major thrusts: 1) using new/improved write-optimized data structures and/or algorithms to better man- age and analyze high-speed massive-volume cyberstreams, and 2) adding security features to data structures at minimum cost. Write optimization allows data structures to better use secondary memory to store and search a larger amount of useful information. Secondary memory is large compared to main memory, but data movement to and from secondary memory must be carefully managed to run quickly enough to keep up with fast streams. The first thrust included managing cyberstreams in parallel, both multi-threaded and distributed, and improving the benchmarking infrastructure for testing new streaming data structures. We considered both (near) real-time discovery of particular patterns, and improved logging for improved forensics. We considered two kinds of security-feature problem. The first was high-performance history- independent external-memory data structures. These provide certain protections to data if a disk is stolen. We also prove some trade-offs between speed and security in this setting. The second data-security problem is more secure data look-up in secret-shared data bases. This report summarizes the project's major accomplishments, with the background to under- stand these accomplishments. It gathers the abstracts and references for the six refereed publications that have appeared as part of this work. We summarize several accomplishments that will be submitted for publication. We then archive one piece of partial work that is not likely to be published in the near future: validation of history-independent data structure implementations.
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)