Search NASA⌕ Search

SEARCH · Search NASA

Results for “controlled binary interactions”

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

Skyrmions in an oblique field

Here we study skyrmion states induced by the oblique field in a square magnetic island. To describe stable states and dynamics of the skyrmion, we employ the lattice model that accounts for all relevant interactions and uses parameters in the range of practical interest. Depending on the orientation of the field, two or four spatially separated energy minima emerge in the oblique field. The depth of the energy minima can be controlled by the strength and orientation of the magnetic field. We study the dynamics of the skyrmion and show that it can be moved between any two states by the field gradient. Possible applications for developing binary and quaternary skyrmion-based memory are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Symmetry breaking forms split-off flat bands in quantum oxides controlling metal versus insulator phases

The crystal structure used as input to electronic structure calculation of conventionally bonded solids consists primarily of the standard crystallographic degrees of freedom. Quantum solids sometimes have additional microscopic degrees of freedom (m-DOF) nested within the crystallographic structure. These might include local motifs including positional (Peierls dimers, deformed octahedra, Jahn-Teller distortions), magnetic (local moment configurations) and dipolar (local ferroelectric configurations). Such motifs can be observed experimentally via local probes that avoid averaging, and theoretically as distinct total energy lowering features relative to more simplified “average crystallographic structures”. Here we examine the ability of electronic structure methods independent of strong correlation physics to explain the broad phenomenology of metal-insulator selectivity in quantum oxides by energy-lowering symmetry breaking. We do this by avoiding the restriction of considering only (i) electron-electron interactions in a fixed unresponsive lattice—as done in Mott strong correlation explanations—allowing, however, (ii) local positional and magnetic motifs that coexist within a crystallographic landscape. We find in a broad range of quantum oxides with different symmetry breaking modes the formation of split-off flat bands that can also control metallic versus insulating characteristics. The split-off band effect is common to magnetic and non-magnetic materials, including both binary and ternary oxides. One finds that when such local symmetry breaking motifs are considered in mean-field-like (e.g. density functional theory) electronic structure calculations of quantum oxides, they explain many of the observed trends in metal vs insulator phases discussed otherwise in prior literature via strong correlation in symmetry unbroken view.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Effects of Perfluorinated Alkyl Substances (PFAS) on Amphibian Body and Liver Conditions: Is Lipid Metabolism Being Perturbed throughout Metamorphosis?

Per- and polyfluoroalkyl substances (PFAS) may interact with peroxisome proliferator activated receptors (PPARs) and alter lipid homeostasis. Using Xenopus laevis, we investigated the effect of PFAS on (a) lipid homeostasis and whether this correlated to changes in body and hepatic condition; (b) the expression of hepatic genes regulated by PPAR; and (c) the hepatic lipidome. We chronically exposed tadpoles to 0.5 µg/L of either PFOS, PFHxS, PFOA, PFHxA, a binary mixture of PFOS and PFHxS (0.5 µg/L of each), or a control, from NF stage 52 through metamorphic climax. Growth, development, and survival were not affected, but we detected a sex-specific decrease in body condition at NF 66 (6.8%) and in hepatic condition (16.6%) across metamorphic climax for male tadpoles exposed to PFOS. We observed weak evidence for the transient downregulation of apolipoprotein-V (apoa5) at NF 62 in tadpoles exposed to PFHxA. Acyl-CoA oxidase 1 (acox1) was downregulated only in males exposed to PFHxS (Ln(Fold Change) = −0.54). We detected PFAS-specific downregulation of structural glycerophospholipids, while semi-quantitative profiling detected the upregulation in numerous glycerophospholipids, sphingomyelins, and diglycerides. Overall, our findings indicate that PFAS can induce sex-specific effects that change across larval development and metamorphosis. We demonstrate that PFAS alter lipid metabolism at environmentally relevant concentrations through divergent mechanisms that may not be related to PPARs, with an absence of effects on body condition, demonstrating the need for more molecular studies to elucidate mechanisms of PFAS-induced lipid dysregulation in amphibians and in other taxa.

Bushong, Anna (ORCID:0000000249630578)↗
DOE OSTI2023

Interactions between nonfullerene acceptors lead to unstable ternary organic photovoltaic cells

For organic photovoltaic (OPV) devices to achieve consistent performance and long operational lifetimes, organic semiconductors must be processed with precise control over their purity, composition, and structure. This is particularly important for high volume solar cell manufacturing where control of materials quality has a direct impact on yield and cost. Ternary-blend OPVs containing two acceptor–donor–acceptor (A–D–A)-type nonfullerene acceptors (NFAs) and a donor have proven to be an effective strategy to improve solar spectral coverage and reduce energy losses beyond that of binary-blend OPVs. Here, we show that the purity of such a ternary is compromised during blending to form a homogeneously mixed bulk heterojunction thin film. We find that the impurities originate from end-capping C=C/C=C exchange reactions of A–D–A-type NFAs, and that their presence influences both device reproducibility and long-term reliability. The end-capping exchange results in generation of up to four impurity constituents with strong dipolar character that interfere with the photoinduced charge transfer process, leading to reduced charge generation efficiency, morphological instabilities, and an increased vulnerability to photodegradation. As a consequence, the OPV efficiency falls to less than 65% of its initial value within 265 h when exposed to up to 10 suns intensity illumination. We propose potential molecular design strategies critical to enhancing the reproducibility as well as reliability of ternary OPVs by avoiding end-capping reactions.

14 SOLAR ENERGY↗
DOE OSTI2023

Nonadditive Interactions Unlock Small-Particle Mobility in Binary Colloidal Monolayers

We examine the organization and dynamics of binary colloidal monolayers composed of micron-scale silica particles interspersed with smaller-diameter silica particles that serve as minority component impurities. These binary monolayers are prepared at the surface of ionic liquid droplets over a range of size ratios (σ = 0.16–0.66) and are studied with low-dose minimally perturbative scanning electron microscopy (SEM). The high resolution of SEM imaging provides direct tracking of all particle coordinates over time, enabling a complete description of the microscopic state. In these bidisperse size mixtures, particle interactions are nonadditive because interfacial pinning to the droplet surface causes the equators of differently sized particles to lie in separate planes. By varying the size ratio, we control the extent of nonadditivity in order to achieve phase behavior inaccessible to additive 2D systems. Across the range of size ratios, we tune the system from a mobile small-particle phase (σ < 0.24) to an interstitial solid (0.24 < σ < 0.33) and furthermore to a disordered glass (σ > 0.33). These distinct phase regimes are classified through measurements of hexagonal ordering of the large-particle host lattice and the lattice’s capacity for small-particle transport. Altogether, we explain these structural and dynamic trends by considering the combined influence of interparticle interactions and the colloidal packing geometry. Our measurements are reproduced in molecular dynamics simulations of 2D nonadditive disks, suggesting an efficient method for describing confined systems with reduced dimensionality representations.

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

Determining Key Factors for the Open-Loop Control of Molecular Fragmentation Using Shaped Strong Fields

Pulse shaping has long been employed for tailoring femtosecond laser pulses to study and control the fragmentation of polyatomic molecules. In many cases, a physical explanation connecting the properties of the field to the observed control is difficult to ascertain. We utilized 80 bit binary spectral phase functions to parametrize and map the search space, gaining insight into which pulse parameters most impact the ion yield and fragmentation pattern for the relatively large triethylamine [N(C 2 H 5 ) 3 ] molecule. Pulse structures used to control the m/z 86 branching ratio beyond a simple intensity dependence are identified and compared to pump–probe results. All of these findings are explained in terms of control via a dissociative Rydberg state in the neutral molecule. This methodology may be used to discover new control mechanisms and shed light onto which pulse parameters most influence the interaction between strong field lasers and matter.

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

Arsenic stabilization performance of a novel starch-modified Fe-Mn binary oxide colloid

Arsenic (As) is an environmentally hazardous contaminant and can be a serious threat to human health. In China, remediation needs for a large number of As-contaminated sites renders a strong demand for efficient remedial reagents and cost effective approaches. In this study, a novel starch-modified Fe-Mn binary oxide (SFM), an amorphous colloidal material, has been synthesized as a remedial reagent and its As stabilization performance has been evaluated. A set of laboratory batch experiments were carried out by directly adding SFM of different dosages into three contaminated soils. Results demonstrated that SFM could transform As in soil from non-specifically and specifically sorbed fractions to the more stable form bounded to amorphous iron hydrous oxides, thus reducing the As concentration in TCLP leachates by 42.0-93.2%. Results from adsorption tests and microscopic analysis indicated that the interactions between SFM and As are mainly controlled by adsorption, oxidation, and precipitation processes. SFM has abundant surface hydroxyl groups, with excellent adsorption properties for both As(V) and As(III), with the maximum adsorption capacities of 160.63 and 284.64 mg/g respectively at pH 7.0. The adsorption process closely fitted pseudo second-order kinetics and Freundlich isotherm model. SFM could increase soil Eh and oxidize As(III) to As(V), which facilitated the As stabilization in soilColloidal iron-based material directly used for stabilization in As contaminated soils is reported here for the first time. Starch modification improves both the reactivity and mobility of SFM in soil. Our findings propose an efficient and convenient reagent for As remediation in soil.

Yan, Xiulan↗
DOE OSTI2020

Negative thermal expansion and the role of hybridization in perovskite-type PbTiO 3 -Bi(Cu 0.5 Ti 0.5) O 3

PbTiO 3 -BiMeO 3 ferroelectrics have attracted much attention due to not only their extremely large polarization and piezoelectricity but also their controllable thermal expansion. In this study, we present a new binary system of (1 - x)PbTiO 3 -xBi(Cu 0.5 Ti 0.5 )O 3 possessing large tetragonality and enhanced spontaneous polarization. Interestingly, a combination of near zero thermal expansion and negative thermal expansion over a temperature range of RT to T C (~550 °C) has been observed in tetragonal 0.85PbTiO 3 -0.15Bi(Cu 0.5 Ti 0.5 )O 3 . The comprehensive structure analysis addressed that the abnormal thermal expansion is coupled to the spontaneous volume ferroelectrostriction. Moreover, the chemical bonding investigation by a maximum entropy method further attests that the covalent hybridization of Pb/Bi–O bonds plays a vital role in the improved ferroelectricity and interacts with thermal vibration to dominate the crossover from positive to negative thermal expansion. fINALLY, The present studies would give direct evidence to verify the theoretical calculations and could be extended to reveal the nature of thermal expansion in ferroelectrics.

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

Oxygen–Chlorine Chemisorption Scaling for Seawater Electrolysis on Transition Metals: The Role of Redox

To clarify what controls species oxidation selectivity in seawater electrolysis, density functional theory (DFT) is used to identify chemisorption enthalpy trends and scaling relations for the simplest relevant adsorbates (O, Cl, and H) on relevant surfaces of 3d transition metals, as well as Pd and Pt, in face-centered-cubic and, if different, their ground-state crystal structures. Approximations are tested for electron exchange-correlation (XC) and van der Waals interactions to assess their ability to reproduce experimental adsorption enthalpies of H and O on Pt(111). The vdW-uncorrected generalized gradient approximation to XC of Perdew, Burke, and Ernzerhof (PBE) agrees most closely with experiments. Using DFT-PBE thereafter, it is determined that the O chemisorption enthalpy on this wide range of transition-metal surfaces is proportional to the sum of first and second atomic ionization energies, akin to a Born–Haber cycle for a redox reaction, indicating that metal redox activity controls O chemisorption strength. Then it is shown that the O and Cl chemisorption enthalpies are strongly correlated, suggesting that the transition metals considered will oxidize unselectively water and Cl – . This strong correlation appears also for crystal reduction potentials of binary oxides and chlorides, indicating a fundamental challenge for future seawater electrode materials design.

08 HYDROGEN↗
DOE OSTI2021

New Compounds and Phase Selection of Nickel Sulfides via Oxidation State Control in Molten Hydroxides

Molten salts are promising reaction media candidates for the discovery of novel materials; however, they offer little control over oxidation state compared to aqueous solutions. Here, we demonstrated that when two hydroxides are mixed, their melts become fluxes with tunable solubility, which are surprisingly powerful solvents for ternary chalcogenides and offer effective paths for crystal growth to new compounds. We report that precise control of the oxidation state of Ni is achievable in mixed molten LiOH/KOH to grow single crystals of all known ternary K-Ni-S compounds. It is also possible to access several new phases, including a new polytope of β-K 2 Ni 3 S 4 , as well as low-valence KNi 4 S 2 and K 4 Ni 9 S 11 . KNi 4 S 2 is a two-dimensional low-valence nickel-rich sulfide, and β-K 2 Ni 3 S 4 has a hexagonal lattice. Moreover, using KNi 4 S 2 as a template, we obtained a new layered binary Ni 2 S by topotactic deintercalation of K. The new binary Ni 2 S has a van der Waals gap and can function as a new host layer for intercalation chemistry, as demonstrated by the intercalation of LiOH between its layers. The oxidation states of low-valence KNi 4 S 2 and Ni 2 S were studied using X-ray absorption spectroscopy and X-ray photoelectron spectroscopy. Density functional theory calculations showed large antibonding interactions at the Fermi level for both KNi 4 S 2 and Ni 2 S 2 corresponding to the flat-bands with large Ni-d$_{x^2}$$_{y^2}$ character.

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

ACES: Infrastructure As Code. Model Optimization and Performance Capability

Infrastructure as Code (IaC) refers to managing infrastructure (networks, physical/virtual machines, storage, and connection topology) in a descriptive model/language, rather than configuring it manually or using interactive configuration tools. Just like source code can be compiled to generate the same binary code, IaC enables generating the same environment every time it is applied. IaC is a key DevOps practice and is generally used in conjunction with continuous integration (CI) and continuous delivery (CD). In CI, all code changes are merged into a mainline branch and validated multiple times a day as developers check in their changes to the source code. In CD on the other hand, code changes are automatically packaged for a new release-to-production on a regular basis. This typically enables teams to deliver software changes much more quickly and often. Developing and managing the ACES platform using (IaC) is vital for the robust deployment and continued sustainment of this foundational computing capability. IaC and DevOps practices will help us solve many of the common challenges often encountered in developing and maintaining compute infrastructure. First, it will make the provisioning, deployment, and maintenance of the compute infrastructure across multiple environments much more efficient. Second, these processes help make the overall system much more stable by continuously testing new changes as they are introduced to the system. Third, it allows us to be much more confident of the security controls in place since they can be tested as part of the CI process and all new changes can be audited and tracked. Finally, IaC enables the ACES Platform to be adaptable to the emerging technologies due to its ability to spin up different test beds to evaluate and incorporate these technologies. This document addresses common infrastructure-management challenges, describes what happens if they are not addressed, and highlights the value of utilizing IaC to tackle them. Finally, we will provide a high-level overview of the IaC and DevOps practices being utilized by the ACES Platform team.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Autoignition study of iso -cetane/tetralin blends at low temperature

iso-Cetane and tetralin are the two fuel components commonly considered in literature diesel surrogates, and are hydrocarbon classes representative of iso-alkanes and naphthoaromatics, respectively. Since both surrogate components are involved in the key reactions controlling the autoignition reactivity of diesel surrogates, autoignition studies of iso-cetane, tetralin, and their blending behavior at low temperatures can be helpful in developing/refining chemical kinetic models of diesel surrogates, in order to better predict diesel ignition response under low temperature combustion conditions. In this investigation, autoignition experiments of iso-cetane, tetralin, and their binary mixtures at varying blending ratios have been conducted in a rapid compression machine (RCM) in the temperature/pressure range of 630–930 K and 10–20 bar. For the neat fuel components, the newly-acquired RCM data have been compared with the literature experimental results at some overlapping conditions. In addition, a chemical kinetic model describing the oxidation of iso-cetane and tetralin compiled/updated from the diesel surrogate model developed by Lawrence Livermore National Laboratory is used to simulate the experimental datasets. It is shown that this model exhibits good agreement with the experimental results of tetralin and also predicts the total ignition delay times of iso-cetane at lower temperatures reasonably well. However, the current model fails to fully capture the increase of total ignition delay time with increasing amount of tetralin in the binary blends observed in the present RCM experiments. Therefore, chemical kinetic analyses of the current model are conducted to help identify possible reasons leading to the discrepancies between simulations and experiments. Furthermore, the chemical kinetic interactions between iso-cetane and tetralin are explored by including additional cross-reactions among iso-cetane, tetralin, and their fuel fragments in the current model to assess their effects on model predictions. As a result, based on the chemical kinetic analysis results, future direction for model refinements is discussed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Unusual electronic transport in (1 - x )Cu 2 Se-( x )CuInSe 2 hierarchical composites

The ability to control the relative density of electronic point defects as well as their energy distribution in semiconductors could afford a systematic modulation of their electronic, optical, and optoelectronic properties. Using a model binary hybrid system Cu 2 Se–CuInSe 2 , we have investigated the correlation between phase composition, microstructure, and electronic transport behavior in the synthesized composites. We found that both Cu 2 Se and CuInSe 2 phases coexist at multiple length scales, ranging from sub-ten nanometer to several micrometers, leading to the formation of a hybrid hierarchical microstructure. Astonishingly, the electronic phase diagram of the (1 - x)Cu 2 Se–(x)CuInSe 2 (15% ≤ x ≤ 100%) hierarchical composites remarkably deviates from the trend normally expected for composites between a heavily doped semiconductor (Cu 2 Se) and a poorly conducting phase (CuInSe 2 ). A sudden 3-fold increase in the electrical conductivity and carrier concentration along with a marginal increase in the carrier mobility is observed for composites at the vicinity of equimolar composition (48% ≤ x ≤ 52%). The carrier concentration increases from ~1.5 × 10 20 cm -3 for the composites with x ≤ 45% to 5.0 × 10 20 cm -3 for x = 50%, and remains constant at 4.5 × 10 20 cm -3 with x value in the range of 52% < x ≤ 90%, then quickly drops to 8 × 10 18 cm -3 for pristine CuInSe 2 phase (x = 100%). The atypical electronic behavior was rationalized in the light of the formation of an inter-band (IB) within the band gap, which arises from the hybridization of native electronic point defects from both Cu 2 Se and CuInSe 2 phases in the resulting hierarchical composites. The result points to a new strategy to modulate the electronic structure of semiconductor composites to maximize interaction and coupling between two fundamentally contrasting properties enabling access to electronic hybrid systems with potential applications as interactive and stimuli-responsive multifunctional materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

From Structure to Function: Zn/Mn-Modified Maghemite as an Advanced Nanoplatform for Magnetic Hyperthermia and Radionuclide Therapy

The development of nanoplatforms capable of efficient heat generation and stable radionuclide delivery is essential for effective bimodal cancer therapy. Here, in this study, binary (Fe–M) and ternary (Fe–M–M′) metal oxide nanoparticles were synthesized via a polyol method optimized to produce flower-like γ-Fe 2 O 3 (maghemite) structures, with M and M′ representing Zn and/or Mn. Comprehensive structural and magnetic characterization was conducted to explain the relationship between composition, defect structure, and hyperthermic performance. The analyses revealed that cation substitution induced an Fe-site vacancy, primarily at octahedral positions, leading to local structural distortions, as confirmed by powder X-ray diffraction and pair distribution function analysis. The optimized composition, with Zn/Mn/Fe = 0.040:0.182:1, exhibited the highest concentration of vacancies and structural disorder. These vacancies altered the bonding environment, enhancing magnetic interactions at tetrahedral sites while weakening those at the octahedral positions. The resulting multicore nanoflowers (20–63 nm; core size 13–18 nm) displayed strong heating performance, with intrinsic loss power ranging from 0.34 to 5.77 nHm 2 kg –1 . The optimized sample achieved a temperature increase of 30 °C within 2 min and a specific absorption rate of 369 W g –1 . This composition was further coated with citrate (CA) and successfully radiolabeled with 177 Lu, achieving a radiolabeling yield of 92.7% and excellent stability, thus forming a robust nanoplatform for combined magnetic hyperthermia and radionuclide therapy. Biological evaluation of the optimized S5 composition revealed selective cytotoxicity toward HeLa and LS174 cells, while toxicity was significantly lower to A549, A375, and normal MRC-5 cells. Citrate coating of S5 nanoparticles (S5@CA) drastically reduced their cytotoxicity across all tested cell lines (IC 50 > 200 μg mL –1 ), confirming their enhanced biocompatibility for therapeutic applications. In HeLa cells subjected to magnetic hyperthermia, the viability decreased to approximately 84% after 30 min and 61% after 60 min of treatment, demonstrating the sustained hyperthermic effect at a controlled working temperature of 48 °C. These results underscore the effectiveness of cation substitution and vacancy engineering in tailoring the functional properties of maghemite-based nanomaterials for advanced multimodal cancer therapies.

36 MATERIALS SCIENCE↗
DOE OSTI2023

NbTa_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-tantalum (Nb-Ta). 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. For Ta, p semi-core states as valence were chosen. 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 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 ternary alloys is AXBYCZ, where A, B, and C refer to the constituents, and X, Y, and Z are positive integers that represent the number of atoms for each constituent and their values still sum up to 128. Each atomic structure associated with a specific chemical composition has its own sub-directory within the directory of the corresponding chemical composition. The sub-directories for each atomic structure for each chemical composition are named 'case-*', where * is a positive integer that spans all the values from 1 through 100, extremes included. The files contained in each sub-directory 'case-*' for each atomic structure are as follows: 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. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and 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. and 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. and 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. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2023

TaV_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-vanadium (Tanu-V). 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. 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 element V. For Ta, p semi-core states as valence were chosen. 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. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and 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. and Furthmuller, 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. and Furthmuller, 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. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
DOE OSTI2023

NbV_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-vanadium (Nb-V). 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 V. 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. Subdirectories with name case-N, where N ranges between 1 and 10 (extremes included) contain all the files listed above except 'vaspout.h5'. Subdirectories with name case-N, where N ranges between 41 and 60 (extremes included), contain a duplicate copy of the files listed above except for KPOINTS. The names of the duplicate files end with -bis, and correspond to a second VASP calculation that has converged to a different optimized geometry. 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. and Hafner, J. Ab initio molecular dynamics for liquid metals. Phys. review B 47, 558 (1993). (2) Kresse, G. and 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. and 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. and 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. and Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. review b 59, 1758 (1999)

36 MATERIALS SCIENCE↗
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↗