Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “First principles”

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

Towards a Deeper Fundamental Understanding of (Al,Sc)N Ferroelectric Nitrides

Density functional theory (DFT) calculations, within the virtual crystal alloy approximation, are performed, along with the development of a Landau-type model employing a symmetry-allowed analytical expression of the internal energy and having parameters determined from first principles, to investigate properties and energetics of Al1-xScxN ferroelectric nitrides in their hexagonal forms. These DFT computations and this model predict the existence of two different types of minima, namely, the fourfold-coordinated wurtzite (WZ) polar structure and a five-fold coordinated paraelectric hexagonal phase (denoted as H5), for any Sc composition up to 40%. The H5 minimum progressively becomes the lowest-energy state within hexagonal symmetry as the Sc concentration increases from 0 to 0.4. Furthermore, the model points to several key findings. Examples include the crucial role of the coupling between polarization and strains to create the WZ minimum, in addition to polar and elastic energies, and that the origin of the H5 state overcoming the WZ phase as the global minimum within hexagonal symmetry when increasing the Sc composition mostly lies in the compositional dependency of only two parameters-one linked to the polarization and another one being purely elastic in nature. Other examples are that forcing Al1-xScxN systems to have no or a weak change in lattice parameters when heating them allows us to reproduce their finite-temperature polar properties well and that a value of the axial ratio close to that of the ideal WZ structure implies a large polarization at low temperatures but not necessarily at high temperatures because of the ordered-disordered character of the temperature-induced formation of the WZ state. Such findings should allow for a better fundamental understanding of (Al,Sc)N ferroelectric nitrides, which may be used to design efficient devices having, e.g., low operating voltages.

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

Exploring New Physics on Cosmological Scales

In recent work with her group, the PI developed a general formalism to compute from first principles the projected mass density (convergence) power spectrum of the substructure in galactic halos under different populations of dark matter sub halos. She constructed a halo model-based formalism, computing the 1-subhalo and the 2-subhalo terms from first principles for the first time. She found that the asymptotic slope of the substructure power spectrum at large wave number reflects the internal density profile of the sub halos, and proposed this as a key observable to discern between different dark matter scenarios.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Effects of pressure on the electronic and magnetic properties of bulk NiI 2

Transition metal dihalides have recently garnered interest in the context of two-dimensional van der Waals magnets as their underlying geometrically frustrated triangular lattice leads to interesting competing exchange interactions. In particular, NiI 2 is a magnetic semiconductor that has been long known for its exotic helimagnetism in the bulk. Recent experiments have shown that the helimagnetic state survives down to the monolayer limit with a layer-dependent magnetic transition temperature that suggests a relevant role of the interlayer coupling. Here, we explore the effects of hydrostatic pressure as a means to enhance this interlayer exchange and ultimately tune the electronic and magnetic response of NiI 2 . We study first the evolution of the structural parameters as a function of external pressure using first-principles calculations combined with x-ray diffraction measurements. We then examine the evolution of the electronic structure and magnetic exchange interactions via first-principles calculations and Monte Carlo simulations. Further, we find that the leading interlayer coupling is an antiferromagnetic second-nearest-neighbor interaction that increases monotonically with pressure. The ratio between isotropic third- and first-nearest-neighbor intralayer exchanges, which controls the magnetic frustration and determines the magnetic propagation vector q of the helimagnetic ground state, is also enhanced by pressure. As a consequence, our Monte Carlo simulations show a monotonic increase in the magnetic transition temperature, indicating that pressure is an effective means to tune the magnetic response of NiI 2 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Microscopic origin of the high thermoelectric figure of merit of n -doped SnSe

Excellent thermoelectric performance in the out-of-layer n-doped SnSe has been observed experimentally [Chang et al., Science 360, 778-783 (2018)]. Though, a first-principles investigation of the dominant scattering mechanisms governing all thermoelectric transport properties is lacking. In the present work, by applying extensive first-principles calculations of electron-phonon coupling associated with parameterized calculation of the scattering by ionized impurities, we investigate the reasons behind the superior figure of merit as well as the enhancement of zT above 600 K in n-doped out-of-layer SnSe, as compared to p-doped SnSe with similar carrier densities. For the n-doped case, the relaxation time is dominated by ionized impurity scattering and increases with temperature, a feature that maintains the power factor at high values at higher temperatures and simultaneously causes the carrier thermal conductivity at zero electric current (κ el ) to decrease faster for higher temperatures, leading to an ultrahigh-zT = 3.1 at 807 K. We rationalize the roles played by κ el and κ 0 (the thermal conductivity due to carrier transport under isoelectrochemical conditions) in the determination of zT. Our results show the ratio between κ 0 and the lattice thermal conductivity indeed corresponds to the upper limit for zT, whereas the difference between calculated zT and the upper limit is proportional to κ el .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Accelerating Multiscale Materials Modeling with Machine Learning

The focus of this project is to accelerate and transform the workflow of multiscale materials modeling by developing an integrated toolchain seamlessly combining DFT, SNAP, LAMMPS, (shown in Figure 1-1) and a machine-learning (ML) model that will more efficiently extract information from a smaller set of first-principles calculations. Our ML model enables us to accelerate first-principles data generation by interpolating existing high fidelity data, and extend the simulation scale by extrapolating high fidelity data (10 2 atoms) to the mesoscale (10 4 atoms). It encodes the underlying physics of atomic interactions on the microscopic scale by adapting a variety of ML techniques such as deep neural networks (DNNs), and graph neural networks (GNNs). We developed a new surrogate model for density functional theory using deep neural networks. The developed ML surrogate is demonstrated in a workflow to generate accurate band energies, total energies, and density of the 298K and 933K Aluminum systems. Furthermore, the models can be used to predict the quantities of interest for systems with more number of atoms than the training data set. We have demonstrated that the ML model can be used to compute the quantities of interest for systems with 100,000 Al atoms. When compared with 2000 Al system the new surrogate model is as accurate as DFT, but three orders of magnitude faster. We also explored optimal experimental design techniques to choose the training data and novel Graph Neural Networks to train on smaller data sets. These are promising methods that need to be explored in the future.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Ab initio ultrafast spin dynamics in solids

Spin relaxation and decoherence is at the heart of spintronics and spin-based quantum information science. Currently, theoretical approaches that can accurately predict spin relaxation of general solids including necessary scattering pathways and are capable of nanosecond to millisecond simulation time are urgently needed. We present a first-principles real-time density-matrix approach based on Lindblad dynamics to simulate ultrafast spin dynamics for general solid-state systems. Through the complete first-principles descriptions of pump, probe, and scattering processes including electron-phonon, electron-impurity, and electron-electron scatterings with self-consistent electronic spin-orbit couplings, our method can directly simulate the ultrafast pump-probe measurements for coupled spin and electron dynamics over nanoseconds at any temperatures and doping levels. We first apply this method to a prototypical system GaAs and obtain excellent agreement with experiments. We found that the relative contributions of different scattering mechanisms and phonon modes differ considerably between spin and carrier relaxation processes. In sharp contrast to previous work based on model Hamiltonians, we point out that the electron-electron scattering is negligible at room temperature but becomes dominant at low temperatures for spin relaxation in n-type GaAs. We further examine ultrafast dynamics in novel spin-valleytronic materials: monolayer and bilayer WSe 2 with realistic defects. We find that spin relaxation is highly sensitive to local symmetry and chemical bonds around defects. For the bilayer WSe 2 , we identify the scattering pathways in ultrafast dynamics and determine relevant dynamical properties, essential to its utilization of unique spin-valley-layer locking effects. In conclusion, our work provides a predictive computational platform for spin dynamics in solids, which has potential for designing new materials ideal for spintronics and quantum information technology.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

On the elastic anisotropy of the entropy-stabilized oxide (Mg, Co, Ni, Cu, Zn)O compound

Herein, we study the elastic properties of the entropy-stabilized oxide (Mg, Co, Ni, Cu, Zn)O using experimental and first principles techniques. Our measurements of the indentation modulus on grains with a wide range of crystallographic orientations of the entropy-stabilized oxide revealed a high degree of elastic isotropy at ambient conditions. First principles calculations predict mild elastic anisotropy for the paramagnetic structure, which decreases when the system is considered to be non-magnetic. When the antiferromagnetic state of CoO, CuO, and NiO is accounted for in the calculations, a slight increase in elastic anisotropy is observed, suggesting a coupling between magnetic ordering and the orientation dependent elastic properties. Furthermore, an examination of the local structure reveals that the isotropy is favored through local ionic distortions of Cu and Zn—due to their tendencies to form tenorite and wurtzite phases. The relationships between the elastic properties of the multicomponent oxide and those of its constituent binary oxides are reviewed. These insights open up new avenues for controlling isotropy for technological applications through tuning composition and structure in the entropy-stabilized oxide or the high-entropy compounds in general.

Entropy↗
DOE OSTI2025

Cost analysis of alternative large-scale high-temperature solid oxide electrolysis hydrogen production facilities

We extend our past cost analysis of gigawatt-scale solid oxide electrolysis (SOE) facilities that produce high purity hydrogen gas from water by estimating construction and operating costs for three new alternative design cases: (1) offsite feed steam generation; (2) near-atmospheric pressure (NAP) stack; and (3) onsite electric boiler feed steam generation. Pressure effects on hydrogen electrode-(cathode-)supported SOE cell (SOEC) stack performance are estimated for the same assumed cell and stack construction and used to determine facility-wide stack capital costs for achieving a fixed H2 production at different pressures. Modular balance of plant (BOP) process equipment capital costs are estimated for each new alternative design case using our past equipment sizing, design, and cost data and scaling relationships. Furthermore, we update BOP equipment sizing and design for the NAP case using Aspen®. Vendor quotes for electric boilers are used to estimate costs for the electric boiler design case. Factory and onsite assembly and installation costs for SOEC stacks and BOP equipment are calculated using our past simplified first-principles approach. First-of-a-kind (FOAK) and N th -of-a-kind (NOAK) production maturity cost estimates are included for all results. The case with NAP stacks offers the lowest facility total capital cost (TCC, ~23% lower than base) while use of small electric boilers requires the highest TCC (~3% higher than base). H 2 production prices decrease from the base of ~$\$2.17$ /kgH 2 to ~$\$1.92$/kgH 2 for 1 GW e DC SIP facilities utilizing NAP stacks supplied by offsites steam situated in large modules and blocks for $\$0.030$/kWh e and $0.009/kWh t prices for electricity and thermal energy, respectively. We report all costs in 2021 US dollars.

Balance of plant (BOP) process equipment↗
DOE OSTI2023

Thermodynamics of diamond formation from hydrocarbon mixtures in planets

Hydrocarbon mixtures are extremely abundant in the Universe, and diamond formation from them can play a crucial role in shaping the interior structure and evolution of planets. With first-principles accuracy, we first estimate the melting line of diamond, and then reveal the nature of chemical bonding in hydrocarbons at extreme conditions. We finally establish the pressure-temperature phase boundary where it is thermodynamically possible for diamond to form from hydrocarbon mixtures with different atomic fractions of carbon. Notably, here we show a depletion zone at pressures above 200 GPa and temperatures below 3000 K-3500 K where diamond formation is thermodynamically favorable regardless of the carbon atomic fraction, due to a phase separation mechanism. The cooler condition of the interior of Neptune compared to Uranus means that the former is much more likely to contain the depletion zone. Our findings can help explain the dichotomy of the two ice giants manifested by the low luminosity of Uranus, and lead to a better understanding of (exo-)planetary formation and evolution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Near-room-temperature ferromagnetic ordering in the pressure-induced collapsed-tetragonal phase in SrC⁢o 2 ⁢P 2

We present high-pressure electrical transport, magnetization, and single-crystal x-ray diffraction data on SrC⁢o 2 ⁢P 2 single crystals, revealing a transition from the uncollapsed tetragonal (ucT) to the collapsed tetragonal (cT) phase >∼10 GPa. This structural transition is accompanied by ferromagnetic (FM) ordering with a Curie temperature up to 260 K. First-principles-based density functional theory calculations show that there is a first-order transition between ucT and cT phases, with an onset ∼10 GPa as well as the appearance of FM ordering in the cT phase. Above ∼30 GPa, the experimental signatures of the magnetic ordering vanish in a first-order-like manner, consistent with the theoretical calculation results, indicating that SrC⁢o 2 ⁢P 2 is another example of the avoidance of quantum criticality in FM intermetallic compounds with almost room-temperature critical temperature. SrC⁢o 2 ⁢P 2 provides clear evidence that the structural, electronic, and magnetic properties associated with the cT transition are strongly entangled and are not only qualitatively captured by our first-principles-based calculations but are quantitatively reproduced as well.

density functional theory↗
DOE OSTI2022

Response to Comment on “Reversible disorder-order transitions in atomic crystal nucleation”

Yu et al. suggested calculating precisely the size ranges of the three parts of our figure 3A, adjusting the free-energy levels in figure 3B, and considering the shape effect in the first-principles calculation. The first and second suggestions raise strong concerns for misinterpretation and overinterpretation of our experiments. Finally, the original calculation is sufficient to support our claim about crystalline-to-disordered transformations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Electronic and magnetic properties of the topological semimetal SmMg 2 Bi 2

Dirac semimetals show nontrivial physical properties and can host exotic quantum states like Weyl semimetals and topological insulators under suitable external conditions. Here, by combining angle-resolved photoemission spectroscopy measurements (ARPES) and first-principle calculations, we demonstrate that the Zintl-phase compound SmMg 2 Bi 2 is in close proximity to a topological Dirac semimetallic state. ARPES results show a Dirac-like band crossing at the zone center near the Fermi level (E F ), which is further confirmed by first-principle calculations. Theoretical studies also reveal that SmMg 2 Bi 2 belongs to a Z 2 topological class and hosts spin-polarized states around the E F . Zintl's theory predicts that the valence state of Sm in this material should be Sm 2+ , however, we detect many Sm-4f multiplet states (flat-bands) whose energy positions and relative intensities suggest the presence of both dominant Sm 2+ and minor Sm 3+ . The small concentration (2.5%) of Sm 3+ in the bulk of a crystal is inferred to arise from Sm vacancies in the crystal. It is also evident that these flat bands and other dispersive states are strongly hybridized when they cross each other. Due to the presence of Sm 3+ ions, the temperature dependence of the magnetic susceptibility χ (T) shows a Curie-Weiss-like contribution in the low-temperature region, in addition to the Van Vleck-like behavior expected for the Sm 2+ ions. The present study will help to better understand the electronic structure, magnetism, and transport properties of related materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Is Ba 3 In 2 O 6 a high- T c superconductor?

It has been suggested that Ba 3 In 2 O 6 might be a high-Tc superconductor. Experimental investigation of the properties of Ba 3 In 2 O 6 was long inhibited by its instability in air. Recently epitaxial Ba 3 In 2 O 6 with a protective capping layer was demonstrated, which finally allows its electronic characterization. The optical bandgap of Ba3In2O6 is determined to be 2.99 eV in-the (001) plane and 2.83 eV along the c-axis direction by spectroscopic ellipsometry. First-principles calculations were carried out, yielding a result in good agreement with the experimental value. Various dopants were explored to induce (super-)conductivity in this otherwise insulating material. Neither A- nor B-site doping proved successful. The underlying reason is predominately the formation of oxygen interstitials as revealed by scanning transmission electron microscopy and first-principles calculations. Additional efforts to induce superconductivity were investigated, including surface alkali doping, optical pumping, and hydrogen reduction. To probe liquid-ion gating, Ba 3 In 2 O 6 was successfully grown epitaxially on an epitaxial SrRuO 3 bottom electrode. So far none of these efforts induced superconductivity in Ba 3 In 2 O 6 , leaving the answer to the initial question of whether Ba 3 In 2 O 6 is a high-T c superconductor to be 'no' thus far.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Topological Circular Dichroism in Chiral Multifold Semimetals

Uncovering the physical contents of the nontrivial topology of quantum states is a critical problem in condensed matter physics. Here, we study the topological circular dichroism in chiral semimetals using linear response theory and first-principles calculations. We show that, when the low-energy spectrum respects emergent SO(3) rotational symmetry, topological circular dichroism is forbidden for Weyl fermions, and thus is unique to chiral multifold fermions. This is a result of the selection rule that is imposed by the emergent symmetry under the combination of particle-hole conjugation and spatial inversion. Using first-principles calculations, we predict that topological circular dichroism occurs in CoSi for photon energy below about 0.2 eV. In conclusion, our Letter demonstrates the existence of a response property of unconventional fermions that is fundamentally different from the response of Dirac and Weyl fermions, motivating further study to uncover other unique responses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Effect of interparticle fields and radiation reaction on beam dynamics

The dynamics of relativistic particles in an intense electromagnetic field can be described by the Landau-Lifshitz (LL) equation, where the radiation reaction (RR) is accounted for via a self-force, and interparticle fields are often neglected as an approximation. However, the inclusion of interparticle fields is necessary to ensure energy-momentum conservation, particularly during coherent emission. Here we present (i) an analytical proof showing that the energy-momentum conservation law of the Hamilton-Rohrlich-Dirac action, which is divergence free and describes a generic system of interacting charges, respects causality and provides physically sensible results; (ii) a simple generalization of the LL equation for many particles evaluated as a function of the total field, i.e., the sum of the external and interparticle fields. By performing first-principles numerical simulations of a neutral, relativistic bunch of electrons and positrons (e − /e + ) colliding with a laser pulse, this theory is shown to satisfy energy-momentum conservation when interparticle fields and RR are simultaneously taken into account; and (iii) the combined effect of interparticle fields and RR primarily affects the tail of the particle energy distribution. Additionally, our first-principles simulations show that the effect of interparticle fields on beam energy loss becomes smaller when most of the radiated energy is incoherent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Fermiology and transport properties of the candidate topological crystalline insulator SrAg 4 ⁢Sb 2

Compared to time-reversal symmetry-protected ℤ 2 topological insulators and Dirac/Weyl semimetals, there are significantly fewer candidates for topological crystalline insulators. SrAg 4 ⁢Sb 2 is predicted to exhibit topological crystalline insulator behavior when considering spin-orbit coupling. In this study, we systematically investigate single crystals of SrAg 4 ⁢Sb 2 using electrical transport and magnetic torque measurements, along with first-principles calculations. Our transport data reveals its compensated semimetal nature with a magnetoresistance up to around 700% at 2 K and 9 T. Analysis of de Haas–van Alphen oscillations uncovers a Fermi surface consisting of three distinct Fermi pockets with light effective masses. Comparison between the three-dimensional fermiology obtained from our oscillation data and the first-principles calculations demonstrates excellent agreement. This confirms the accuracy of the calculations, which indicate a band inversion centered at the 𝑇 point and identify the existence of nontrivial tube and needle hole Fermi pockets at Γ, alongside one trivial diamond electron pocket at the 𝐹 point in the Brillouin zone. Furthermore, symmetry and topology analysis results in two potential sets of topological invariants, suggesting the emergence of two-dimensional gapless Dirac surface states either on the 𝑎⁢𝑏 planes or on both the 𝑎⁢𝑏 planes and mirror planes, protected by crystal symmetries. Therefore, SrAg 4 ⁢Sb 2 emerges as a promising candidate topological crystalline insulator.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Ferroelectric HfO 2 and the importance of strain

Ferroelectric oxides based on HfO 2 show tremendous promise for the next generation of memory and logic devices. The ferroelectric polymorph is one of several that can be derived from the high symmetry cubic fluorite structure of HfO 2 . A single grain of HfO 2 may consist of a coherent mixture of multiple orientational and translational variants of different polymorphs. Here, we use symmetry-adapted strain-order parameters to elucidate the relationship between the different HfO 2 polymorphs and their symmetrically equivalent variants. We use first-principles electronic structure methods to identify minimum energy pathways and map them in subspaces of the symmetry-adapted strain order parameters. We next investigate the atomic structure of domain boundaries that separate coexisting variants of ferroelectric HfO 2 . Further, we rely on Gibbsian excess quantities and a precise specification of mechanical boundary conditions to describe the thermodynamic properties of domain boundaries. Our first-principles calculations show that the O and Hf shuffle arrangement within a domain boundary is closely related to the intermediate shuffle patterns of the homogeneous pathways between ferroelectric variants. Furthermore, the preferred structure within a boundary is very sensitive to local strain constraints imposed by the adjacent ferroelectric variants, leading to highly anisotropic domain boundary energies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

NbTaV_BCC_SolidSolution_128atoms_VASP6

We performed density functional theory (DFT) calculations for body-centered-cubic (BCC) structures with 128 lattices sites of solid solution ternary alloys niobium-tantalum-vanadium (Nb-Ta-V). The first-principle code that has been used to run the calculations is the closed-source Vienna Ab-Initio Simulation Package (VASP). 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 8. 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 10,500 randomized atomic structures over 105 chemical compositions. The calculations have been collected on NERSC-Perlmutter 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. Every chemical composition sampled across the composition range in the dataset has its own directory. The convention used to name the directories for ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z 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: 1. INCAR: input file that contains various parameters and settings for controlling the behavior of the electronic structure calculations 2. 0.POSCAR: input file that defines the atomic structure of a system 3. 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 4. 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 5. 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 6. 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 7. 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 8. OUTCARL: 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 9. 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 10. 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 11. vasprun.xml: contains similar information to OUTCAR, but in an xml format. 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.

36 MATERIALS SCIENCE↗