Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atomic structure”

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 1,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Photoelectric characteristics of silicon carbide–silicon structures grown by the atomic substitution method in a silicon crystal lattice

Data obtained in an experimental study of the photoelectric characteristics of silicon–silicon carbide structures grown by the atomic substitution method on silicon (100) and (111) substrates are presented. It is found that the maximum sunlight conversion efficiency of a silicon–silicon carbide (silicon carbide–silicon) heterojunction is 5.4%. The theory of dilatation dipole formation upon synthesis by the atomic substitution method is used to account for the mechanism of electrical barrier formation at the silicon–silicon carbide interface.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Crystal Structure of an i‐Motif from the HRAS Oncogene Promoter

Abstract An i‐motif is a non‐canonical DNA structure implicated in gene regulation and linked to cancers. The C‐rich strand of the HRAS oncogene, 5′‐CG CCC GTG CCC TGCG CCC GCAA CCC GA‐3′ (herein referred to as iHRAS), forms an i‐motif in vitro but its exact structure was unknown. HRAS is a member of the RAS proto‐oncogene family. About 19 % of US cancer patients carry mutations in RAS genes. We solved the structure of iHRAS at 1.77 Å resolution. The structure reveals that iHRAS folds into a double hairpin. The two double hairpins associate in an antiparallel fashion, forming an i‐motif dimer capped by two loops on each end and linked by a connecting region. Six C−C + base pairs form each i‐motif core, and the core regions are extended by a G−G base pair and a cytosine stacking. Extensive canonical and non‐canonical base pairing and stacking stabilizes the connecting region and loops. The iHRAS structure is the first atomic resolution structure of an i‐motif from a human oncogene. This structure sheds light on i‐motifs folding and function in the cell.

Li, Kevin S.↗
DOE OSTI2023

Crystal Structure of an i‐Motif from the HRAS Oncogene Promoter

Abstract An i‐motif is a non‐canonical DNA structure implicated in gene regulation and linked to cancers. The C‐rich strand of the HRAS oncogene, 5′‐CG CCC GTG CCC TGCG CCC GCAA CCC GA‐3′ (herein referred to as iHRAS), forms an i‐motif in vitro but its exact structure was unknown. HRAS is a member of the RAS proto‐oncogene family. About 19 % of US cancer patients carry mutations in RAS genes. We solved the structure of iHRAS at 1.77 Å resolution. The structure reveals that iHRAS folds into a double hairpin. The two double hairpins associate in an antiparallel fashion, forming an i‐motif dimer capped by two loops on each end and linked by a connecting region. Six C−C + base pairs form each i‐motif core, and the core regions are extended by a G−G base pair and a cytosine stacking. Extensive canonical and non‐canonical base pairing and stacking stabilizes the connecting region and loops. The iHRAS structure is the first atomic resolution structure of an i‐motif from a human oncogene. This structure sheds light on i‐motifs folding and function in the cell.

Li, Kevin S.↗
DOE OSTI2021

Electron diffraction radial distribution function analysis of amorphous boron carbide synthesized by ion beam irradiation and chemical vapor deposition

Amorphous boron carbide (a-BxC) networks consist of light elements, and their low atomic scattering factors makes structural analysis by x-ray diffraction difficult. Electron diffraction has an advantage of detecting the light elements, because of the strong interaction between the matter and electrons. In this work, we prepared a-BxC by ion beam technologies and plasma-enhanced chemical vapor deposition, and characterized their structures via atomic pair-distribution functions derived from electron diffraction intensity profiles. It was found that a pentagonal pyramid is the most favorable cluster in a-B4C generated by ion irradiation, while C—C homonuclear bonds were formed in the deposited a-B x C thin film. X-ray photoemission spectroscopy revealed that the a-BxC thin film possesses more carbon than B 4 C, which is responsible for the formation of the homonuclear bonds.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Mapping the effects of physical and chemical reduction parameters on local atomic distributions within bimetallic nanoparticles

Bimetallic nanoparticles prove advantageous over their monometallic counterparts due to the tunable, hybrid properties that result from combining different atomic species in a controlled way. The favorable optical and catalytic properties resulting from AgAu nanoparticle formation have been widely attributed to the existence of Ag–Au bonds, the maximization of which assumes the formation of a homogeneous alloy. Despite the importance of atomic scale structure in these systems, synthetic studies are typically not paired with structural characterization at the atomic scale. Herein, a comprehensive synthetic exploration of physical and chemical reduction parameters of resulting nanoparticle products is complemented with thorough X-ray characterization to probe how these parameters affect atomic scale alloy distributions within AgAu nanoparticles. Presented evidence shows Ag is substantially underincorporated into nanoparticle constructs compared with solution Ag : Au ratios regardless of precursor :reductant ratio or volume of reductant added. Both Ag and Au exhibit significant local clustering, with Ag distributed preferentially towards the nanoparticle surface. Most significantly, the results of this investigation suggest that reduction parameters alone can affect the local alloy distributions and homogeneity within bimetallic nanoparticles, even when the ratio of metallic precursors remains constant. Overall, this investigation presents the ability to control alloy distributions using kinetics and provides new considerations for optimizing synthetic methods to produce functional bimetallic nanoparticles.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Nuclear deformation as a source of the nonlinearity of the King plot in the Yb + ion

We perform atomic relativistic many-body calculations of the field isotope shifts and calculations of corresponding nuclear parameters for all stable even-even isotopes of the Yb + ion. We demonstrate that if we take nuclear parameters of the Yb isotopes from a range of the state of the art nuclear models, which all predict strong quadrupole nuclear deformation, and then calculate nonlinearity of the King plot caused by the difference in the deformation in different isotopes, the result is consistent with the nonlinearity observed in the experiment. In conclusion, the changes of nuclear rms radius between isotopes extracted from experiment are consistent with those obtained in the nuclear calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Polar catastrophe and the structure of KT a 1 - x N b x O 3 surfaces: Results from elastic and inelastic helium atom scattering

In this paper, the structure and dynamics of cleaved (001) surfaces of potassium tantalates doped with niobium, KTa 1-x Nb x O 3 (KTN), with x ranging from 0% to 30%, were measured by helium atom scattering (HAS). Through HAS time-of-flight (TOF) experiments, a dispersionless branch (Einstein phonon branch) with energy of 13-14meV was observed across the surface Brillouin zone in all samples. When this observation is combined with the results from earlier experimental and theoretical studies on these materials, a consistent picture of the stable surface structure emerges: After cleaving the single-crystal sample, the surface should be composed of equal areas of KO and TaO 2 /NbO 2 terraces. The data, however, suggest that K + and O 2- ions migrate from the bulk to the surface, forming a charged KO lattice that is neutralized primarily by additional K + ions bridging pairs of surface oxygens. This structural and dynamic modification at the (001) surface of KTN appears due to its formally charged KO(-1) and TaO 2 /NbO 2 (+1) layers and avoids a “polar catastrophe.” Finally, this behavior is contrasted with the (001) surface behavior of the fluoride perovskite KMnF 3 with its electrically neutral KF and MnF 2 layers.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Confirmation of a de novo structure prediction for an atomically precise monolayer-coated silver nanoparticle

Fathoming the principles underpinning the structures of monolayer-coated molecular metal nanoparticles remains an enduring challenge. Notwithstanding recent x-ray determinations, coveted veritable de novo structural predictions are scarce. Building on recent syntheses and de novo structure predictions of M 3 Au x Ag 17-x (TBBT) 12 , where M is a countercation, x = 0 or 1, and TBBT is 4-tert-butylbenzenethiol, we report an x-ray–determined structure that authenticates an a priori prediction and, in conjunction with first-principles theoretical analysis, lends force to the underlying forecasting methodology. The predicted and verified Ag(SR) 3 monomer, together with the recently discovered Ag 2 (SR) 5 dimer and Ag 3 (SR) 6 trimer, establishes a family of unique mount motifs for silver thiolate nanoparticles, expanding knowledge beyond the earlier-known Au-S staples in thiol-capped gold nanoclusters. These findings demonstrate key principles underlying ligand-shell anchoring to the metal core, as well as unique T-like benzene dimer and cyclic benzene trimer ligand bundling configurations, opening vistas for rational design of metal and alloy nanoparticles.

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

Electronic Structure of Single-Atom Alloys and Its Impact on The Catalytic Activities

Single-atom alloys (SAAs) are promising materials for heterogeneous catalysis due to their unique structure and electronic properties. SAAs have active sites narrowed down to the single-atom level, which combines the advantages of alloy materials and single-site catalysts. Given the unique structural feature of SAAs, their electronic properties can be more flexibly tailored than for their monometallic counterparts, which can be used to effectively control their catalytic activities. One interesting feature commonly observed for SAAs is the lower density of state (DOS) near the Fermi level than their bulk references. Comparing with results for their monometallic bulk reference, the most noticeable electronic property change in SAAs is the narrowing of the valence band, which gives them free-atom-like character. Moreover, the d-band position of both single atoms and their host metals can show a pronounced shift. These changes of electronic structure in SAAs could largely affect the adsorption behavior of adsorbates during the catalytic processes. Close examination of the relationship between electronic structure and catalytic activity can provide useful guidance for rational design of new catalysts with improved performance.

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

Unraveling the 13C NMR Chemical Shifts in Single-Walled Carbon Nanotubes: Dependence on Diameter and Electronic Structure

The atomic specificity afforded by nuclear magnetic resonance (NMR) spectroscopy could enable detailed mechanistic information about single-walled carbon nanotube (SWCNT) functionalization as well as the noncovalent molecular interactions that dictate ground-state charge transfer and separation by electronic structure and diameter. However, to date, the polydispersity present in as-synthesized SWCNT populations has obscured the dependence of the SWCNT {sup 13}C chemical shift on intrinsic parameters such as diameter and electronic structure, meaning that no information is gleaned for specific SWCNTs with unique chiral indices. In this article, we utilize a combination of {sup 13}C labeling and density gradient ultracentrifugation (DGU) to produce an array of {sup 13}C-labeled SWCNT populations with varying diameter, electronic structure, and chiral angle. We find that the SWCNT isotropic {sup 13}C chemical shift decreases systematically with increasing diameter for semiconducting SWCNTs, in agreement with recent theoretical predictions that have heretofore gone unaddressed. Furthermore, we find that the {sup 13}C chemical shifts for small diameter metallic and semiconducting SWCNTs differ significantly, and that the full-width of the isotropic peak for metallic SWCNTs is much larger than that of semiconducting nanotubes, irrespective of diameter.

14 SOLAR ENERGY↗
NASA NTRS1997

Doping Scheme in Atomic Chain Electronics

Due to the dramatic reduction in MOS size, there appear many unwanted effects. In these small devices, the number of dopant atoms in the channel is not macroscopic and electrons may suffer significantly different scattering from device to device since the spatial distribution of dopant atoms is no longer regarded as continuous. This prohibits integration, while it is impossible to control such dopant positions within atomic scale. A fundamental solution is to create electronics with simple but atomically precise structures, which could be fabricated with recent atom manipulation technology. All the constituent atoms are placed as planned, and then the device characteristics are deviation-free, which is mandatory for integration. Atomic chain electronics belongs to this category. Foreign atom chains or arrays form devices, and they are placed on the atomically flat substrate surface. We can design the band structure and the resultant Fermi energy of these structures by manipulating the lattice constant. Using the tight-binding theory with universal parameters, it has been predicted that isolated Si chains and arrays are metallic, Mg chains are insulating, and Mg arrays have metallic and insulating phases [1]. The transport properties along a metallic chain have been studied, emphasizing the role of the contact to electrodes [2]. For electronic applications, it is essential to establish a method to dope a semiconducting chain, which is to control the Fermi energy position without altering the original band structure. If we replace some of the chain atoms with dopant atoms randomly, the electrons will see random potential along die chain and will be localized strongly in space (Anderson localization). However, if we replace periodically, although the electrons can spread over the chain, there will generally appear new bands and band gaps reflecting the new periodicity of dopant atoms. This will change the original band structure significantly. In order to overcome this dilemma, we may place a dopant atom beside the chain at every N lattice periods (N > 1). Because of the periodic arrangement of pant atoms, we can avoid the unwanted Anderson localization. Moreover, since the dopant atoms do not constitute the chain, the overlap interaction between them is minimized, and the band structure modification can be made smallest. Some tight-binding results will be discussed to demonstrate the present idea.

Toshishige, Yamada↗
DOE OSTI2017

Intermetallic Pd3Pb Nanowire Networks Boost Ethanol Oxidation and Oxygen Reduction Reaction with Significantly Improved Methanol Tolerance

Intermetallic nanocrystals are currently receiving extensive attentions due to their well- defined crystal structures, highly ordered atomic distribution and superior structural stability that endow them with optimized catalytic activities, stabilities and high selectivity as electrocatalysts for fuel cells. Here, for the first time, we reported the facile synthesis of intermetallic Pd3Pb nanowire networks (IM-Pd3Pb NNs) with a one-step wet-chemical strategy at a relatively low temperature (i.e. 170 °C) in 1 h. The as-prepared IM-Pd3Pb NNs exhibited superior bifunctional catalytic performances toward oxygen reduction reaction (ORR) and ethanol oxidation reaction (EtOR) compared to commercial Pt/C and Pd black, respectively. Significantly, IM-Pd3Pb NNs also showed excellent methanol- and CO-tolerant ability as ORR cathode and EtOR anode electrocatalysts, respectively. The electrochemical active surface area and mass activity of IM-Pd3Pb NNs are about 3.4 times and 2 times higher than Pd black toward EtOR, respectively. As the Pt-free bifunctional electrocatalysts, 3D IM-Pd3Pb architectures with exceptional catalytic performances hold great promise in various applica-tions of energy conversion and storage devices, sensors, electronics, optics and so on.

Shi, Qiurong↗
DOE OSTI2016

Structural properties of lithium atom under weakly coupled plasma environment

The Rayleigh-Ritz variational technique with a Hylleraas basis set is being tested for the first time to estimate the structural modifications of a lithium atom embedded in a weakly coupled plasma environment. The Debye-Huckel potential is used to mimic the weakly coupled plasma environment. The wave functions for both the helium-like lithium ion and the lithium atom are expanded in the explicitly correlated Hylleraas type basis set which fully takes care of the electron-electron correlation effect. Due to the continuum lowering under plasma environment, the ionization potential of the system gradually decreases leading to the destabilization of the atom. The excited states destabilize at a lower value of the plasma density. The estimated ionization potential agrees fairly well with the few available theoretical estimates. The variation of one and two particle moments, dielectric susceptibility and magnetic shielding constant, with respect to plasma density is also been discussed in detail.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Using Complementary Methods of Synchrotron Radiation Powder Diffraction and Pair Distribution Function to Refine Crystal Structures with High Quality Parameters—A Review

Determination of the atomic-scale structures of certain fine-grained minerals using single-crystal X-ray diffraction (XRD) has been challenging because they commonly occur as submicron and nanocrystals in the geological environment. Synchrotron powder diffraction and scattering techniques are useful complementary methods for studying this type of minerals. In this review, we discussed three example studies investigated by combined methods of synchrotron radiation XRD and pair distribution function (PDF) techniques: (1) low-temperature cristobalite; (2) kaolinite; and (3) vernadite. Powder XRD is useful to determine the average structure including unit-cell parameters, fractional atomic coordinates, occupancies and isotropic atomic displacement parameters. X-ray/Neutron PDF methods are sensitive to study the local structure with anisotropic atomic displacement parameters (ADP). The results and case studies suggest that the crystal structure and high-quality ADP values can be obtained using the combined methods. The method can be useful to characterize crystals and minerals that are not suitable for single-crystal XRD.

58 GEOSCIENCES↗
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↗