Search NASA⌕ Search

SEARCH · Search NASA

Results for “quantum compilation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Infrared /IR/ spectroscopy in support of atmospheric measurements

In general, the data required for the more abundant species are: (1) linewidths as a function of temperature; (2) line strength measurements obtained at various temperatures when since quantum assignments for the lines are unknown the lower state energy levels are unknown; and (3) line strength measurements of the 'hot' bands and bands of the isotopic species. It is pointed out that an extensive list of line center frequencies and line strengths from 2400 to 3200 per cm of (CH4)4 has recently been added to the Air Force Geophysics Laboratory compilation. Since many of the lines in the list could not be assigned, the laboratory research required measurements of samples at various temperatures in order to experimentally determine the temperature dependence on the line strengths. The 'hot' band here refers to transitions arising from vibrational states other than the ground state. 'Hot' bands are usually near the parent band (arising from the ground state). Very few 'hot' bands are without parent bands.

Toth, R. A.↗
NASA NTRS1995

Low Intensity Low Temperature (LILT) measurements and coefficients on new photovoltaic structures

Past NASA missions to Mars, Jupiter, and the outer planets were powered by radioisotope thermal generators (RTG's). Although these devices proved to be reliable, their high cost and highly toxic radioactive heat source has made them far less desirable for future planetary missions. This has resulted in a renewed search for alternate energy sources, some of them being photovoltaic (PV) and thermophotovoltaic (TPV). Both of these alternate energy sources convert light/thermal energy directly into electricity. In order to create a viable PV and TPV data base for planetary mission planners and cell designers, we have compiled low temperature low intensity (LILT) I-V data on single junction and multi-junction high efficiency solar cells. The cells tested here represent the latest photovoltaic technology. Using this LILT data to calculate dI(sub SC)/dT, dV(sub OC)/dT, dFF/dT, and also as a function of intensity, an accurate prediction of cell performance under the AMO spectrum can be determined. When combined with QUantum efficiency at Low Temperature (QULT) data, one can further enhance the data by adding spectral variations to the measurements. This paper presents an overview of LILT measurements and is only intended to be used as a guideline for material selection and performance predictions. As single junction and multi-junction cell technologies emerge, new test data must be collected. Cell materials included are Si, GaAs/Ge, GainP/GaAs/Ge, InP, InGaAs/InP, InP/InGaAs/InP, and GainP. Temperatures range as low as -175 C and intensities range from 1 sun to .02 suns.

Schelman, David A.↗
NASA NTRS1995

Low Intensity Low Temperature (LILT) Measurements and Coefficients on New Photovoltaic Structures

Past NASA missions to Mars, Jupiter and the outer planets were powered by radioisotope thermal generators (RTGs). Although these devices proved to be reliable, their high cost and highly toxic radioactive heat source has made them far less desirable for future planetary missions. This has resulted in a renewed search for alternate energy sources, some of them being photovoltaics (PV) and thermophotovoltaics (TPV). Both of these alternate energy sources convert light/thermal energy directly into electricity. In order to create a viable PV data base for planetary mission planners and cell designers, we have compiled low intensity low temperature (LILT) I-V data on single junction and multi-junction high efficiency solar cells. The cells tested here represent the latest photovoltaic technology. Using this LILT data to calculate Short Circuit Current (I(sub sc)), Open Circuit Voltage (V(sub os)), and Fill Factor (FF) as a function of temperature and intensity, an accurate prediction of cell performance under the AM0 spectrum can be determined. When combined with QUantum efficiency at Low Temperature (QULT) data, one can further enhance the data by adding spectral variations to the measurements. This paper presents an overview of LILT measurements and is only intended to be used as a guideline for material selection and performance predictions. As single junction and multi-junction cell technologies emerge, new test data must be collected. Cell materials included are Si, GaAs/Ge, GaInP/GaAs/GaAs, InP, InGaAs/InP, InP/InGaAs/InP, and GaInP. Temperatures range down to as low as -180 C and intensities range from 1 sun down to 0.02 suns. The coefficients presented in this paper represent experimental results and are intended to provide the user with approximate numbers.

Scheiman, David A.↗
NASA NTRSDate not supplied

Symmetry and Error Mitigation in Quantum Alternating Operator Ansatz

Quantum Alternating Operator Ansatz (QAOA) provides a platform for probing quantum advantage in optimization problems on NISQ hardware. By exploring symmetries in the system and the structure of unitaries in the algorithm, we proposed using problem-specific mixing operators to limit the quantum evolution within a subspace where the target state lives. As an example, in case of equality constraint, we studied the use of XY operators that preserve total Pauli Z. This QAOA variant, in the noiseless case, significantly outperforms the vanilla version of QAOA where the constraints are encoded as penalty in the cost Hamiltonian and the search is in the whole Hilbert space. In presence of noise, although the advantage still holds, error mitigation is needed for QAOA. I will overview our work on the design of QAOA in presence of symmetry, the implementation and compilation of XY-QAOA, its performance and error mitigation techniques in presence of noise.

Zhihui Wang↗
NASA NTRS1993

Theoretical calculations on the electron absorption spectra of selected Polycyclic Aromatic Hydrocarbons (PAH) and derivatives

As a theoretical component of the joint effort with the laboratory of Dr. Lou Allamandola to search for potential candidates for interstellar organic carbon compound that are responsible for the visible diffuse interstellar absorption bands (DIB's), quantum mechanical calculations were performed on the electron absorption spectra of selected polycyclic aromatic hydrocarbons (PAH) and derivatives. In the completed project, 15 different species of naphthalene, its hydrogen abstraction and addition derivatives, and corresponding cations and anions were studied. Using semiempirical quantum mechanical method INDO/S, the ground electronic state of each species was evaluated with restricted Hartree-Fock scheme and limited configuration interaction. The lowest energy spin state for each species was used for electron absorption calculations. Results indicate that these calculations are accurate enough to reproduce the spectra of naphthalene cation and anion observed in neon matrix. The spectral pattern of the hydrogen abstraction and addition derivatives predicted based on these results indicate that the electron configuration of the pi orbitals of these species is the dominant determinant. A combined list of 19 absorptions calculated from 4500 A to 10,400 A were compiled and suggested as potential candidates that are relevant for the DIB's absorptions. Continued studies on pyrene and derivatives revealed the ground state symmetries and multiplicities of its neutral, anionic, and cationic species. Spectral calculations show that the cation (B(sub 3g)-2) and the anion (A(sub u)-2) are more likely to have low energy absorptions in the regions between 10 kK and 20 kK, similar to naphthalene. These absorptions, together with those to be determined from the hydrogen abstraction and addition derivatives of pyrene, can be used to provide additional candidates and suggest experimental work in the search for interstellar compounds that are responsible for DIB's.

Du, Ping↗
DOE OSTI2026

Quantum Choreography of the Nucleus: Rotations, Vibrations, and Emergent Structure

Abstract Nuclei are complex many-body quantum systems where interactions of the neutrons and protons via the strong, the weak, and the electromagnetic forces lead to the emergence of simple patterns of energy states that have been described by various theoretical approaches. One of the goals of all the theoretical models is the development of a universal theory that can be applied across the entire chart of nuclides. Significant progress has been made by experiments as well as the increasing sophistication of models, but a universal theory has yet to be established. A recent review of nuclei in the Z = 50–82 region of the chart of nuclides has analyzed all the available compiled data from several decades of studies towards a clarification of the low-lying structure of nuclei. Other reviews have reported and explained the emergence of multiple different shapes in nuclei at somewhat higher excitation energies than the ground state. Some have challenged the interpretation of the first excited K π = 0 + band as a vibration of ground state. This work attempts to provide a guide to determining the nature of the first excited K π = 0 + band in nuclei by the combined use of nuclear lifetimes, energy level evolutions, dynamic moments of inertia, and intrinsic quadrupole moments extracted from transition probabilities. The result is that for a subset of the nuclei in this region, the K π = 0 + band is consistent with the traditional 𝛽-vibration description of an oscillation built on the ground state.

CHFB+5DCH↗
NASA NTRS2020

Propionitrile in the Two Lowest Excited Vibrational States in the Laboratory and on Titan

The rotational spectrum of propionitrile (ethyl cyanide, C2H5CN) in the two lowest excited vibrational states, ν(sub 13) = 1 and ν(sub 21) = 1, was studied in detail. The compiled experimental data set consisted of published line lists, lines measured from published spectra up to 645 GHz, and new measurements at frequencies from 8 to 331 GHz. The two states are subject to considerable mutual interaction but application of a- and b-axis Coriolis coupling model to the measured A-symmetry internal rotation substate transitions was found to be successful. A total of over 1900 transitions were fitted for J values up to 76 and Ka up to 24, resulting in determined vibrational separation between the two states ΔE = 6:2070246(14) cm(exp -1), and Coriolis coefficients |ζ(sup a)(sub 13,21)| = 0.378, |ζ(sup b) (sub 13,21)| = 0.127, in good agreement with quantum chemistry calculations. The new line catalog was used to model ALMA spectra of Titan’s atmosphere from 2016, fitting many weak lines not previously modeled. By including the vibrationally excited lines of C2H5CN and a new partition function, we derive a disk-averaged volume mixing ratio estimate in 2016 of 7.8 ± 0.2 ppb assumed constant above 300 km.

Analysis of rotational spectrum↗
NASA NTRS1972

Explicit expressions of the potential and its derivatives at the origin in terms of the scattering data

The quantum mechanical theory of scattering of a particle by a spherically symmetrical potential is presented. As in the inverse scattering problem, the input of the calculation is the scattering and bound-state data, and the output is data on the potential. The results discussed are explicit expressions for the values of the potential and its derivatives at the origin in terms of the scattering and boundstate data. Various methods to obtain these results are outlined. The presentation is aimed at introducing these various approaches. The simplest scattering problem (nonrelativistic S-wave scattering on a holomorphic potential without bound states) is used as the basis for discussion, and technicalities are omitted whenever possible without loss of clarity. A complete compilation is given of the results obtained to date in this field, including the treatment of higher partial waves and the Klein-Gordon and Dirac equations.

Calogero, F.↗
DOE OSTI2026

An Introduction to Variational Calculus

Variational Calculus is an advanced topic in mathematics. This document exists to compile the basics of Variational Calculus in a readable and digestible manner for application by engineers. This document is a living document and will be updated as time moves forward. It is the hope that by reading this document, the basics of Variational Calculus can be readily applied to any system/function of interest. This does not serve as a replacement for textbooks or other learned sources but should act as a companion piece to better serve the reader in times of confusion and distress.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Fault-tolerant resource comparison of qudit and qubit encodings for diagonal quadratic operators

Finite local Hilbert-space truncations arise naturally in quantum simulations of lattice field theories and motivate qudit encodings, but their fault-tolerant advantage over qubit encodings remains unclear. We compare the non-Clifford cost of implementing quadratic diagonal evolutions, exemplified by 𝑈 = 𝑒$^{−𝑖⁢𝑡⁢𝜙^2_𝑥}$ in a uniform field-amplitude discretization of a real scalar field, using either one logical 𝑑-level qudit or 𝑛 𝑏 = ⌈log 2⁡ 𝑑⌉ logical qubits. We analyze two standard settings: product-formula simulation and linear combination of unitaries (LCU) per block encoding, taking the resource metric to be the number of non-Clifford gates after synthesis into a discrete logical gate set. Because tight synthesis bounds for general single-qudit rotations are not known, we express the qudit constructions in terms of embedded two-level SU⁡(2) rotations and derive explicit finite-𝑑 break-even conditions for their synthesis cost; these serve as compiler targets for when qudit encodings can outperform the qubit baseline. Within the constructive models studied here, product-formula implementations would require an exponentially stronger per-primitive synthesis advantage for qudits to win asymptotically, while in the LCU setting the qubit encoding is asymptotically cheaper in 𝑑. Nevertheless, the finite-𝑑 threshold analysis identifies low-dimensional regions in which qudits can yield meaningful constant-factor savings, particularly for LCU-based implementations. As a secondary analysis of the LCU construction, we use an idealized negligible-overhead qubit-qudit code-switching model to give an absolute 𝑇-count comparison and reinterpret the savings as an allowable per-switch overhead budget.

Godwood, Samuel [Univ. of Liverpool (United Kingdo↗
NASA NTRS2020

Disruptive Technologies and Their Putative Impacts Upon Society and Aerospace- Entering The Virtual Age

Developments in technology over the recent decades have been extraordinary. They include the IT, bio, nano, and now quantum and energetics technology arenas and their many combinatorial interactions and impacts. In the main, these are at the frontiers of the small and in a combinational, synergistic feeding frenzy with each other. They fall under the broad category of Disruptive Technologies and have greatly altered society. The outlook for the runout of these and other technology developments augers mid-term to later alterations in components of the human existence theorem, including the requirement to work for our living and our physiological makeup and longevity (Ref 1). The IT revolution began in the 1950s with the development of solid-state electronics. The biologics revolution began later in the 1960s and 1970s with DNA and genomics, and the nano revolution in the 1990s with self-forming nano systems and carbon nanotubes. Quantum technology is now developing rapidly, aided by enabling nano systems, and the energetics revolution is providing ever more efficient and less expensive renewable energy sources. The IT revolution has produced improvements of an astounding eleven orders of magnitude in computing speed since the late 1950s. As we shift from silicon to biological, optical, nano, molecular, and atomic computing, improvements of some 4 orders of magnitude are evidently possible from either optical or DNA computing [Refs 2and 3], then there are combinatorials. Then there is quantum computing, under development worldwide for an increasing number of applications and proffering phenomenal capabilities. The current fastest computers are considerably beyond human brain speed. Machine intelligence is developing well after decades of inadequate machine capability, now no longer the case, and a detour into expert systems. Researchers in machine intelligence are now pursuing deep learning approaches using neural nets, which are proving to be extremely useful. Some believe the frontier of potential human-level machine intelligence may be found in biomimetics and brain-emulation approaches. There is even a possibility of “emergence”—i.e., when the machine intelligence is complex enough that it “wakes up,” as when human intelligence emerged via evolution during the million-plus years of the hunter-gatherer epoch [ Ref 4]. In fact, some posit that human intelligence can be improved upon and is only a cul-de-sac of what is conceivable. The IT revolution has produced massive changes in human society and economics—from the Internet, enabling the rapid expansion of knowledgeability (and even what is knowable), to an increasingly pervasive trend of “tele-everything.” The extraordinary compilation, storage, and availability of truly massive amounts of information could, when combined with AI and under the mantra of “big data,” greatly improve many of our technical and commercial processes and their content including elucidating new heuristic governing laws.

Dennis M. Bushnell↗
DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗
NASA NTRS2007

Comparison of Integrated Radiation Transport Models with TEPC Measurements for the Average Quality Factors in Spaceflights

The purpose of this work is to test our theoretical model for the interpretation of radiation data measured in space. During the space missions astronauts are exposed to the complex field of radiation type and kinetic energies from galactic cosmic rays (GCR), trapped protons, and sometimes solar particle events (SPEs). The tissue equivalent proportional counter (TEPC) is a simple time-dependent approach for radiation monitoring for astronauts on board the International Space Station. Another and a newer approach to Microdosimetry is the use of silicon-on-insulator (SOI) technology launched on the MidSTAR-1 mission in low Earth orbit (LEO). In the radiation protection practice, the average quality factor of a radiation field is defined as a function of linear energy transfer (LET), Q(sub ave)(LET). However, TEPC measures the average quality factor as a function of the lineal energy y, Q(sub ave)(y), defined as the average energy deposition in a volume divided by the average chord length of the volume. Lineal energy, y, deviates from LET due to energy straggling, delta-ray escape or entry, and nuclear fragments produced in the detector volume. Monte Carlo track structure simulation was employed to obtain the response of a TEPC irradiated with charged particle for an equivalent site diameter of 1 micron of wall-less counter. The calculated data of the energy absorption in the wall-less counter were compiled for various y values for several ion types at various discrete projectile energy levels. For the simulation of TEPC response from the mixed radiation environments inside a spacecraft, such as, Space Shuttle and International Space Station, the complete microdosimetric TEPC response, f( y, E, Z), were calculated with the Monte Carlo theoretical results by using the first order Lagrangian interpolation for a monovariate function at a given y value (y = 0.1 keV/micron 5000 keV/micron) at any projectile energy level (E = 0.01 MeV/u to 50,000 MeV/u) of each specific radiation type (Z = 1 to 28). Because the anomalous response has been observed at large event sizes in the experiment due to the escape of energy out of sensitive volume by delta-rays and the entry of delta-rays from the high-density wall into the low-density gas-volume cavity, Monte Carlo simulation was also made for the response of a walled-TEPC with wall thickness 2 mm and density 1 g/cm(exp 3). The radius of cavity was set to 6.35 mm and a gas density 7.874 x 10(exp -5) g/cm(exp 3). The response of the walled- and the wall-less counters were compared. The average quality factor Q(sub ave)(y) for trapped protons on STS-89 demonstrated the good agreement between the model calculations and flight TEPC data as shown. Using an integrated space radiation model (this includes the transport codes HZETRN and BRYNTRN, the quantum nuclear interaction model QMSFRG) and the resultant response distribution functions of walled-TEPC from Monte-Carlo track simulations, we compared model calculations with walled-TEPC measurements from NASA missions in LEO and made predictions for the lunar and the Mars missions. The Q(sub ave)(y) values for the trapped or the solar protons ranged from 1.9-2.5. This over-estimates the Qave(LET) values which ranged from 1.4-1.6. Both quantities increase with shield thickness due to nuclear fragmentation. The Q(sub ave)(LET) for the complete GCR spectra was found to be 3.5-4.5, while flight TEPCs measured 2.9-3.4 for Q(sub ave)(y). The GCR values are decreasing with the shield thickness. Our analysis for a proper interpretation of data supports the use of TEPCs for monitoring space radiation environment.

Kim, Myung-Hee Y.↗
DOE OSTI2024

Determining Exact RANS Operators with the Macroscopic Forcing Method (Final Report)

This report contains a compilation of key results and findings from the ACT project “Determining Exact RANS Operators with the Macroscopic Forcing Method.” The Macroscopic Forcing Method (MFM), a numerical tool for determining closure operators, is used to measure eddy diffusivity moments in Rayleigh-Taylor (RT) instability. It is first applied to low-Atwood 2D RT instability; that work is then extended to 3D RT at different finite Atwood numbers. It is found that nonlocality is important for modeling the mean scalar transport closure operator in RT mixing. Additionally, work is done to improve the statistical convergence of MFM for chaotic problems like RT mixing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Leafweb: Leaf Gas Exchange and Pulse-Amplitude Modulated Fluorometry for C4 Species, June 2026 Release

This dataset contains leaf gas exchange and Pulse-Amplitude Modulated (PAM) fluorometry for 98 C4 species. The C4 photosynthetic pathway employs specialized CO2 concentration mechanisms and Kranz anatomy to enrich CO2 concentration around Rubisco, the enzyme that catalyzes carbon fixation in the Calvin-Benson cycle to suppress photorespiration and increase the use efficiencies of light, nitrogen, and water as compared to the C3 photosynthetic pathways. Large-scale C4 photosynthetic datasets are relatively scarce, which has affected C4 photosynthesis research. To improve C4 photosynthetic data availability, Leafweb organized an effort to systematically collect, compile, standardize, and organize measurements of leaf gas exchange and/or Pulse-Amplitude Modulated (PAM) fluorometry of C4 species. This derived a C4 photosynthetic dataset containing measurements made by independent researchers in multiple countries in various environments (field, garden, or greenhouse). It covers three biochemical subtypes – the nicotinamide adenine dinucleotide phosphate-malic enzyme (NADP-ME), nicotinamide adenine dinucleotide-malic enzyme (NAD-ME), and phosphoenolpyruvate carboxykinase (PEP-CK) subtypes. This dataset is useful for using Artificial Intelligence / Machine Learning and mechanistic models to study C4 photosynthesis and compare across different biochemical subtypes. This dataset contains 3 compressed (*.zip) folders containing 1,892 data files in comma-separate values (*.csv) format. Additional metadata are provided: one data dictionary and a file-level metadata file in comma-separate values (*.csv) format and a user guide in PDF (*.pdf) format.

Zhou, Haoran [Tianjin University, China]↗
DOE OSTI2025

Synthesis of single-qutrit circuits from Clifford+𝑅 gates

Here, we present two deterministic compilation algorithms for single-qutrit unitaries with O ( log 1 / ɛ ) gate depth. Each algorithm selects a nearby approximation to the target unitary and then exactly synthesizes the approximation over the Clifford + R basis. The first algorithm exhaustively searches over the group; while the second algorithm searches only for Householder reflections. The exhaustive search algorithm yields an average R count of 2.193 ( 11 ) + 8.621 ( 7 ) log 10 ( 1 / ɛ ) , albeit with a time complexity of O ( ɛ − 4.4 ) . The Householder search algorithm results in a larger average R count of 3.20 ( 13 ) + 10.77 ( 3 ) log 10 ( 1 / ɛ ) at a reduced time complexity of O ( ɛ − 0.42 ) , greatly extending the reach in ɛ . These costs correspond asymptotically to 35% and 69% more non-Clifford gates compared with synthesizing the same unitary with two qubits. Such initial results are encouraging for using the R gate as the nontransversal gate for qutrit-based computation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

NbZr_BCC_SolidSolution_128atoms_VASP6

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

36 MATERIALS SCIENCE↗
DOE OSTI2024

TaZr_BCC_SolidSolution_128atoms_VASP6

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

36 MATERIALS SCIENCE↗