Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “quantum programming”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

New Laboratory and Astrophysical Probes of Physics Beyond the Standard Model (Final Technical Report)

The first of these directions centers on new ideas in dark matter detection. As experimental scrutiny narrows the window for new physics at the weak scale, it has become apparent that the solution to the dark matter mystery may not reside there. Theoretical developments such as the hidden sector/valley paradigm have, at the same time, shown that compelling theories of dark matter reside below the weak scale. Searching for dark matter at lower mass scales requires looking beyond the current paradigm of dark matter direct detection based on nuclear recoils. It is the goal of this project to provide new ideas to guide the experimental program looking for light dark matter. I am actively proposing ideas for detecting dark matter as light as a meV, well beyond the current focus of GeV-TeV scale dark matter. Some of the ideas I have proposed, such as superconducting, superfluid, and polar material targets, are actively being developed into experiments. My group will provide the relevant calculations to determine dark matter reach, providing crucial input to experiment on which targets should be developed. We will also provide an effective field theory framework to understand which types of experiments are most sensitive to each interaction type. The second of these directions is on probes of dark matter substructure as a mean to constrain, or observe, models of dark matter. The standard ΛCDM paradigm assumes that the dark matter density perturbations are adiabatic, scale invariant, and produced during inflation. However, many standard models of dark matter will produce modifications of this prediction, such as axion models with symmetry broken below the inflationary scale. Remarkably, we have limited direct measurements of the dark matter clumpiness at mass scales below dwarf galaxies. This allows for the possibility of observing modifications from vanilla ΛCDM due to particle dynamics. In many cases, simulations of dark matter structure in the presence of non-scale invariant fluctu- ations are understudied or completely lacking. I plan to develop both theoretical predictions for small scale structure in models with additional matter power on small scales, as well as observa- tional probes of small scale halos or clumps. One idea I have been recently focused on is pulsar timing, though my group will pursue a variety of lensing and astrometric probes. Lastly, the most risky direction involves spacetime fluctuations from quantum gravity. I showed that if one assumes that metric fluctuations in the Minkowski vacuum, at a surface separating a region in and out of causal contact (which we call a “horizon”), are determined by standard thermodynamic considerations at horizon, these metric fluctuations are large enough to observe in an interferometer with similar sensitivity to LIGO. In a follow-up paper we showed that these assumption holds for the vacuum in AdS/CFT. In future work, I plan to connect this work to recent soft graviton results, to frame this work in a concrete Randall-Sundrum model, and to work out concrete phenomenological predictions from the model. This work unquestionably takes a less trodden path, and works well as part of a well-rounded portfolio of less and more risky ideas. Particle physics is currently at a juncture which requires bold exploration of qualitatively new directions. This proposal outlines some of my plans over the coming years in these directions, leaving room also for surprises.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Direct Higgs-top CP-phase measurement with $t\overline{t}h$ at the 14 TeV LHC and 100 TeV FCC

The study of the Higgs boson’s properties is a cornerstone of the LHC and future collider programs. In this paper, we examine the potential to directly probe the Higgs-top interaction strength and CP-structure in the $t\overline{t}h$ channel with the Higgs boson decaying to bottom-quark pairs and top-quarks in the di-leptonic mode. We adopt the BDRS algorithm to tag the boosted Higgs and exploit the M 2 -assisted reconstruction to compute observables sensitive to the CP-phase at the $t\bar{t}$ rest frame, where the new physics sensitivity can be enhanced. Performing a side-band analysis at the LHC to control the continuum $t\bar{t}b\bar{b}$ background, we find that the Higgs-top strength and CP-phase can be probed up to δ κt ≲ 20% and |a| ≲ 36° at 95% CL, respectively. We also derive that a similar analysis at a 100 TeV future collider could further improve the precision to δ κt ≲ 1% and |a| ≲ 1.5° where the CP-odd observables play a crucial role, boosting the sensitivity on the CP-phase

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Electromagnetic signals of inelastic dark matter scattering

Light dark sectors in thermal contact with the Standard Model can naturally produce the observed relic dark matter abundance and are the targets of a broad experimental search program. A key light dark sector model is the pseudo-Dirac fermion with a dark photon mediator. The dynamics of the fermionic excited states are often neglected. We consider scenarios in which a nontrivial abundance of excited states is produced and their subsequent de-excitation yields interesting electromagnetic signals in direct detection experiments. We study three mechanisms of populating the excited state: a primordial excited fraction, a component up-scattered in the Sun, and a component up-scattered in the Earth. We find that the fractional abundance of primordial excited states is generically depleted to exponentially small fractions in the early universe. Nonetheless, this abundance can produce observable signals in current dark matter searches. MeV-scale dark matter with thermal cross sections and higher can be probed by down-scattering following excitation in the Sun. Up-scatters of GeV-scale dark matter in the Earth can give rise to signals in current and upcoming terrestrial experiments and X-ray observations. We comment on the possible relevance of these scenarios to the recent excess in XENON1T.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

LICS Modeling [Slides]

We propose to design a new diagnostic based on laser inverse Compton scattering (LICS) to study the dynamics of runaway electron formation during killer-pellet triggered disruptions in DIII-D, and their subsequent loss. Obtaining a better understanding of transient runaway electrons (and how to suppress them) is one of the highest priorities of the US fusion program. To enable this, we need better measurements of the runaway electron population during and after disruptions. We intend to measure point-wise time-resolved spatial profile and energy distributions of the runaway electrons while they are in the core of the tokamak plasma. While LICS measurements of relativistic e-beams have been successfully made on accelerators, they have never been done in a tokamak. We bring two innovations to this diagnostic concept which should enable 10 times better signal-to-noise ratios than previously envisioned. In particular, we will marry bright short pulse (10’s of picosecond) laser technologies with the availability of short-pulse gated x-ray imagers developed at Los Alamos for the National Ignition Facility (NIF), to reject noise generated by the tokamak and plasma environment. We will generate a Conceptual Physics Basis, a Preliminary Design, and a Final Design package, with associated DIII-D reviews, along with prototype component testing for the diagnostic over the course of the two year duration of this proposal.

42 ENGINEERING↗
DOE OSTI2020

SIERRA Multimechanics Module: Aria Thermal Theory Manual (V.4.56)

Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process flows via the incompressible Navier-Stokes equations specialized to a low Reynolds number ( Re < 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton’s method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic h -adaptivity and dynamic load balancing are some of Aria’s more advanced capabilities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

SIERRA Multimechanics Module: Aria User Manual - Version 5.8

Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process flows via the incompressible Navier-Stokes equations specialized to a low Reynolds number (Re < 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton’s method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic ℎ-adaptivity and dynamic load balancing are some of Aria’s more advanced capabilities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

SIERRA Multimechanics Module: Aria Thermal Theory Manual - Version 5.8

Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process ows via the incompressible Navier-Stokes equations specialized to a low Reynolds number (Re < 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton’s method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic ℎ-adaptivity and dynamic load balancing are some of Aria’s more advanced capabilities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Precision inclusive Higgs physics at e+e- colliders with tracking detectors and without calorimetry

A primary goal of a future e + e - collider program will be the precision measurement of Higgs boson properties. For practical reasons it is of interest to determine the minimal set of detector specifications required to reach this and other scientific goals. Here we investigate the precision obtainable for the e + e - Zh μ + μ - X inclusive cross section and the Higgs boson mass using the di-muon recoil method, considering a detector that has only an inner tracking system within a solenoidal magnetic field, surrounded by many nuclear interaction lengths of absorbing material, and an outer muon identification system. We find that the sensitivity achievable in these measurements with such a tracking detector is only marginally reduced compared to that expected for a general purpose detector with additional electromagnetic and hadronic calorimeter systems. The difference results mainly from multi-photon backgrounds that are not as easily rejected with tracking detectors. We also comment on the prospects for an analogous measurement of the e + e - → Zh → e + e - X inclusive cross section. Finally, we study searches for light scalars utilizing the di-muon recoil method, estimating the projected reach with a tracking or general purpose detector.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

CO adsorption on MnO(100): Experimental benchmarks compared to DFT

In this study, CO adsorption on the MnO(100) surface was studied using temperature programmed desorption (TPD) and density functional theory (DFT). TPD results show that CO is weakly-bound on MnO(100), with an experimental adsorption energy of -35.6 ± 2.1 kJ/mol at terrace sites in the low coverage limit. PBE simulations suggest that CO adsorption causes an implausible (2 × 2) surface reconstruction. PBE+U simulations show no signs of surface reconstruction, and provide an accurate estimate of the adsorption energy (-36.4 kJ/mol) when combined with the DFT-D3 method with Becke-Jonson damping to correct for van der Waals interactions. This simulation also shows that CO adsorbs C-down onto the Mn 2+ terrace site in a tilted geometry, which is also observed experimentally and computationally on the similarly-structured NiO(100) transition metal oxide surface. TPD results for large doses show a plateauing of the coverage at about 0.4 monolayers of CO at 85 K, with a defect coverage equivalent to 0.08 monolayers. Adsorption associated with defect sites is indicated by a high-temperature desorption tail which is not satisfactorily explained by DFT simulations of simple step or oxygen vacancy defects.

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

Broadband Solenoidal Haloscope for Terahertz Axion Detection

We introduce the Broadband Reflector Experiment for Axion Detection (BREAD) conceptual design and science program. This haloscope plans to search for bosonic dark matter across the [ 10 - 3 , 1 ] eV ([0.24, 240] THz) mass range. BREAD proposes a cylindrical metal barrel to convert dark matter into photons, which a novel parabolic reflector design focuses onto a photosensor. This unique geometry enables enclosure in standard cryostats and high-field solenoids, overcoming limitations of current dish antennas. A pilot 0.7 m 2 barrel experiment planned at Fermilab is projected to surpass existing dark photon coupling constraints by over a decade with one-day runtime. Axion sensitivity requires < 10 - 20 W / Hz sensor noise equivalent power with a 10 T solenoid and 10 m 2 barrel. We project BREAD sensitivity for various sensor technologies and discuss future prospects.

Liu, Jesse↗
DOE OSTI2021

WIYN Open Cluster Study. LXXXV. Li in NGC 2243: Implications for Stellar and Galactic Evolution

High-dispersion spectra in the Li λ6708 region have been obtained and analyzed in the old, metal-deficient cluster NGC 2243. From Hydra spectra for 29 astrometric and radial velocity members, we derive rotational velocities, as well as [Fe/H], [Ca/H], [Si/H], and [Ni/H] based on 17, 1, 1, and 3 lines, respectively. Using ROBOSPECT, an automatic equivalent width measurement program, we derive [Fe/H] = −0.54 ± 0.11 (MAD), for an internal precision for the cluster [Fe/H] below 0.03 dex. Given the more restricted line set, comparable values for [Ca/H], [Si/H], and [Ni/H] are −0.48 ± 0.19, −0.44 ± 0.11, and −0.61 ± 0.06, respectively. With E(B − V) = 0.055, appropriate isochrones imply (m − M) = 13.2 ± 0.1 and an age of 3.6 ± 0.2 Gyr. Using available VLT spectra and published Li abundances, we construct an Li sample of over 100 stars extending from the tip of the giant branch to 0.5 mag below the Li dip. The Li dip is well populated and, when combined with results for NGC 6819 and Hyades/Praesepe, implies a mass/metallicity slope of 0.4 M {sub ⊙}/dex for the high-mass edge of the Li dip. The A(Li) distribution among giants reflects the degree of Li variation among the turnoff stars above the Li dip, itself a function of stellar mass and metallicity and strongly anticorrelated with a v {sub rot} distribution that dramatically narrows with age. Potential implications of these patterns for the interpretation of Li among dwarf and giant field populations, especially selection biases tied to age and metallicity, are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The Present and Future of QCD: QCD Town Meeting White Paper – An Input to the 2023 NSAC Long Range Plan

It is currently understood that there are four fundamental forces in nature: gravitational, electromagnetic, weak and strong forces. The strong force governs the interactions between quarks and gluons, elementary particles whose interactions give rise to the vast majority of visible mass in the universe. The mathematical description of the strong force is provided by the non-Abelian gauge theory Quantum Chromodynamics (QCD). While QCD is an exquisite theory, constructing the nucleons and nuclei from quarks, and furthermore explaining the behavior of quarks and gluons at all energies, remain to be complex and challenging problems. Such challenges, along with the desire to understand all visible matter at the most fundamental level, position the study of QCD as a central thrust of research in nuclear science. Experimental insight into the strong force can be gained using large particle accelerator facilities, which are necessary to probe the very short distance scales over which quarks and gluons interact. The Long Range Plans (LRPs) exercise of 1989 and 1996 led directly to the construction of two world-class facilities: the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Lab (JLab) that is focused on studying how the structure of hadrons emerges from QCD (cold QCD research), and the Relativistic Heavy Ion Collider (RHIC) at Brookhaven National Lab (BNL) that aims at the discovery and study of a new state of matter, the quark-gluon plasma (QGP), at extremely high temperatures (hot QCD research). These past investments have produced major advances. Nucleons and nuclei are being studied with increasing precision with a unified description of the partonic structure utilizing multi-dimensional imaging. Significant progress has been made, paving the way towards a complete picture of how quarks and gluons give rise to the mass, spin, and momentum of the nucleon. In hot QCD, the QGP is created in the collisions of nuclei at RHIC and the Large Hadron Collider (LHC) and is observed to behave like a fluid with very low specific shear viscosity; the current goals are to understand how the fluid behavior emerges from QCD and to characterize the temperature (and chemical potential) dependence of the properties of the QGP. As this White Paper is written, current experimental programs at CEBAF, RHIC and the LHC continue to provide exciting near term opportunities to capitalize on the investments in experimental equipment and accelerator operations. Most importantly, the QCD community looks forward to the construction of the Electron Ion Collider (EIC) as a major new facility to push forward QCD research in the next decades, with significant focus on exploring the properties of gluons, the mediators of the strong force.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Tidally Excited Modes and δ Scuti Pulsations in the Eclipsing Triple Star IM Persei

IM Per is a triple star system whose eclipsing pair masses and radii are accurate to within 1%. The Transiting Exoplanet Survey Satellite (TESS) light curve of the program target exhibits partial eclipses and multiple oscillations with millimagnitude-level amplitudes. It is found that the oscillations affect eclipse timing measurements. Binary modeling of the high-quality TESS data indicates that the eclipsing components of the triple system are twin dwarfs with parameters of M {sub 2}/M {sub 1} = 0.995, R {sub 2}/R {sub 1} = 0.901, and Δ(T{sub eff,1}–T{sub eff,2}) = 12 K in an eccentric (e = 0.049), detached configuration. The third light of l {sub 3} = 0.054 may mostly come from a G-type tertiary companion. Our predicted parallax of 1.52 ± 0.09 mas is concurrent with the Gaia measurement of 1.52 ± 0.05 mas. Multifrequency analysis of the outside-eclipse residuals reveals 22 significant pulsation signals: 4 in the gravity-mode region (0.03−2.22 day{sup −1}) and 18 in the pressure-mode region (9.19−25.12 day{sup −1}). Of the low frequencies, f {sub 11} and f {sub 14} are orbital harmonics that can be identified as tidally excited modes. The pulsation periods and constants for the high frequencies, and the position in the Cepheid instability strip demonstrate that the pulsating component of IM Per is a δ Sct variable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Vidyut3d: A GPU accelerated fluid solver for non-equilibrium plasmas on adaptive grids

We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure three-electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate ~ 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Nanoscale engineering and dynamic stabilization of mesoscopic spin textures

Thermalization, while ubiquitous in physics, has traditionally been viewed as an obstacle to be mitigated. In contrast, we demonstrate here the use of thermalization in the generation, control, and readout of “shell-like” spin textures with interacting 13 C nuclear spins in diamond, wherein spins are polarized oppositely on either side of a critical radius. The textures span several nanometers and encompass many hundred spins; they are created and interrogated without manipulating the nuclear spins individually. Long-time stabilization is achieved via prethermalization to a Floquet-engineered Hamiltonian under the electronic gradient field: The texture is therefore metastable and robust against spin diffusion. This enables the state to endure over multiple minutes before it decays. Our work on spin-state engineering paves the way for applications in quantum simulation and nanoscale imaging.

42 ENGINEERING↗
DOE OSTI2020

The Belle II Physics Book

The Belle II Theory Interface Platform (B2TiP) was created as a physics prospects working group of the Belle II collaboration in June 2014. It offered a platform where theorists and experimentalists could work together to elucidate the potential impacts of the Belle II program, which includes a wide scope of physics topics: B physics, charm, τ, quarkonium physics, electroweak precision measurements, and dark sector searches.

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↗