Defect formation energy in pyrochlore: The effect of crystal size
Explore the source record for details and available documents.
nuclear (including radiation effects), materials a↗
SEARCH · Search NASA
Results for “projector augmented wave method”
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.
Halogen diffusion on a Ga-stabilized ζ-GaAs(001)–(4 × 2) surface
The atomic and electronic structure of a Ga-stabilized GaAs(001) surface with the ζ(4 × 2) reconstruction and halogens in a number of symmetric sites on the surface are calculated by the plane-wave projector augmented wave method. The energy barriers of halogen-atom diffusion on this surface are calculated, which allow determination of the most preferred paths of their migration. It is shown that there is a low barrier (0.17–0.23 eV) for the diffusion of all halogens under consideration (I, Br, Cl, F) along the surface gallium dimer, whereas the barrier is significantly higher for diffusion between adjacent gallium dimers. In general, the energy barriers for halogen diffusion in both directions ([110] and [1-10]) point to their high surface mobility, despite high binding energies at a number of surface adsorption sites.
36 MATERIALS SCIENCE↗
A Physical Model for Understanding the Activation of MoS 2 Basal-plane Sulfur Atoms for the Hydrogen Evolution Reaction
All the DFT calculations are done with the Vienna Ab Initio Simulation Package (VASP) using the projector augmented wave method. The Bayesian error estimation exchange-correlation functionals (BEEF) with van der Waals interactions are employed. A plane-wave cutoff energy of 400 eV is used together with PAW-PBE potentials where semi core p states are treated as valence. All the calculations allow for spin-polarization. The structures are relaxed until the force is converged to < 0.01 eV/Å. The lattice parameter of MoS 2 unit cell, optimized with this functional, is 3.19 Å. A (4×4) supercell is used to model all the transition metal doped MoS 2 systems studied here, including those with S vacancies. For calculations in the initial dopant structure exploration, the Brillouin zone is sampled with a 3x3x1 Monkhorst-Pack k-point mesh. A 6×6×1 Monkhorst-Pack k-point mesh is used for the H binding energy and density of states calculations. In all calculations, the vacuum layer is set as 15 Å to eliminate periodic interaction perpendicular to the basal plane.
77 NANOSCIENCE AND NANOTECHNOLOGY↗
Carbon diffusion in molten uranium: an ab initio molecular dynamics study
In this work we used ab initio molecular dynamics (AIMD) within the framework of density functional theory (DFT) and the projector-augmented wave (PAW) method to study carbon diffusion in liquid uranium at temperatures above 1600 K. The electronic interactions of carbon and uranium were described using the local density approximation (LDA). The self-diffusion of uranium based on this approach is compared with literature computational and experimental results for liquid uranium. The temperature dependence of carbon and uranium diffusion in the melt was evaluated by fitting the resulting diffusion coefficients to an Arrhenius relationship. We found that the LDA calculated activation energy for carbon was nearly twice that of uranium: 0.55±0.03 eV for carbon compared to 0.32±0.04 eV for uranium. Structural analysis of the liquid uranium-carbon system is also discussed.
Molecular dynamics computer simulation, Nuclear Pr↗
Assessing the Role of Interfacial and Metal Sites in Pt/TiO2-Catalyzed Acetic Acid Hydrodeoxygenation
Thermochemical conversion of biomass to produce drop-in quality biofuels typically involves hydrodeoxygenation (HDO) steps following catalytic fast pyrolysis (CFP) to remove excess oxygen and create a more-stable bio-oil product. HDO reactions are performed by co-feeding the CFP vapor-phase product and H2 gas over a bi-functional catalyst. Noble-metal catalysts supported on reducible metal oxides (e.g., Pt/TiO2) are active and selective toward these HDO reactions. Griffin and co-workers showed that Pt/TiO2 catalysts promote the desired deoxygenation steps for m-cresol HDO while mitigating undesired C-C bond-breaking steps that reduce the overall value/energy density of the biofuel. Such model-compound studies to inform the design of improved catalysts for HDO chemistry are necessary to improve the overall process economics/efficiencies for biofuels production. One important class of bio-derived compounds that has not been studied extensively with respect to HDO chemistry, particularly at the atomic level, is carboxylic acids. Past experimental work indicates that Pt/C and Pt/TiO2 catalysts are selective toward C-C and C-O bond-dissociation products for acetic acid HDO, respectively. To understand the role of Pt/TiO2 active sites in this observed change in selectivity and guide catalyst development, the work in this presentation focuses on modeling the role of Pt-metal and Pt-TiO2-interfacial sites in promoting key C-C bond-breaking, C-O bond-breaking, and (de)hydrogenation steps in acetic acid HDO. Density functional theory (DFT) calculations were performed using the Vienna Ab initio Simulation Package (VASP). The exchange-correlation functional was approximated by the Perdew-Burke-Ernzerhof functional. Dispersion interactions were captured using the D3 method. Projector augmented-wave potentials described electron-ion interactions, and electron wavefunctions were expanded via a planewave basis with an energy cutoff of 500 eV. Activation barriers were calculated using the climbing image nudged elastic band method. Results and Discussion: To discern the role of Pt-metal and Pt-TiO2 interface sites in promoting acetic acid HDO chemistry, Pt(111) slab and anatase TiO2(101)-supported Pt-nanowire (PtNW/OH-TiO2) surface models were constructed, respectively. Because H2 is co-fed in HDO reactions, the anatase support was terminated with OH groups. Interfacial vacancies have been shown to facilitate Ru/TiO2-catalyzed phenol HDO; thus, an interfacial model with an OH vacancy was also considered (PtNW/OHv-TiO2). Pt-TiO2-interface sites stabilize adsorption of all studied acetic acid HDO surface intermediates relative to terrace Pt-metal sites, particularly when an interfacial-OH vacancy is present. Oxygenated species prefer to bind at the OH vacancy through the O atom, suggesting a preference for C-O over C-C bond cleavage at these sites. This hypothesis is supported by net-negative and net-positive average shifts in the reaction energy and activation energy barriers for C-O and C-C bond-breaking steps, respectively, at Pt-TiO2-interface sites relative to Pt-metal sites. Using the calculated energetics, the predicted minimum-energy pathway was determined for each surface. Pt(111) is predicted to follow decarboxylation, producing undesired methane and carbon dioxide. Conversely, PtNW/OH-TiO2 and PtNW/OHv-TiO2 are both predicted to produce desired acetaldehyde and ethane. The interfacial vacancy may also play a key role in facilitating the first C-O bond-breaking step in acetic acid HDO, lowering the barrier by 0.6 eV relative to the defect-free interface model. These results demonstrate the critical role of the Pt-TiO2 interface in the shift in acetic acid HDO selectivity experimentally observed on Pt/C and Pt/TiO2 catalysts. The results herein demonstrate the important role of the Pt-TiO2 interface and interfacial oxygen vacancies in improving the carbon efficiency for HDO reactions in CFP upgrading.
BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗
Relativistic nuclear magnetic resonance J-coupling with ultrasoft pseudopotentials and the zeroth-order regular approximation
We present a method for the first-principles calculation of nuclear magnetic resonance (NMR) J-coupling in extended systems using state-of-the-art ultrasoft pseudopotentials and including scalar-relativistic effects. The use of ultrasoft pseudopotentials is allowed by extending the projector augmented wave (PAW) method of Joyce et al. [J. Chem. Phys. 127, 204107 (2007)]. We benchmark it against existing local-orbital quantum chemical calculations and experiments for small molecules containing light elements, with good agreement. Scalar-relativistic effects are included at the zeroth-order regular approximation level of theory and benchmarked against existing local-orbital quantum chemical calculations and experiments for a number of small molecules containing the heavy row six elements W, Pt, Hg, Tl, and Pb, with good agreement. Finally, {sup 1}J(P-Ag) and {sup 2}J(P-Ag-P) couplings are calculated in some larger molecular crystals and compared against solid-state NMR experiments. Some remarks are also made as to improving the numerical stability of dipole perturbations using PAW.
97 MATHEMATICS AND COMPUTING↗
Materials Design and Discovery: Catalysis and Energy Storage (Mira Early Science Program Final Technical Report): ALCF-2 Early Science Program Technical Report
The investigation and design of new classes of materials for energy and catalysis requires a multi-facetted approach to simulation. Multiple methods are needed to study materials on the length scale 0.1 nm - 10 nm. For simulations where the atomic (and electronic) degrees of freedom are relevant, the methods of choice in the surface science, condensed matter physics, and material science communities are classical molecular dynamics (CMD), Density Functional Theory (DFT), and quantum Monte Carlo (QMC). The original scope for this Early Science Program (ESP) project was to perform fast-accurate DFT calculations on materials for energy and catalysis using the GPAW code on Blue Gene/Q. The types of calculations included significantly reduced time-to-solution on systems sizes accessible on Blue Gene/P (~10, 0000 valence electrons), but also systems which were were at least a factor of two larger (~20, 000 valence electrons). GPAW is a real-space DFT code using the projector augmented wave (PAW) method. DFT calculations on Blue Gene/P were executed on over >100, 000 cores using GPAW; thus it was consider a success on Mira’s predecessor system, Intrepid. One of the co-PIs (NAR), determined that the work necessary to allow the GPAW code for these aforementioned types of calculations could not be accomplish within the time frame of the ESP. This was not simply due to human time required to implement OpenMP parallelism, but also from intrinsic algorithmic limitations in supporting libraries, most notably ScaLAPACK. Additionally, the return-on-investment on $\mathcal{O}(N^3)$ DFT code has become some what tenable at best. NAR argues that what is really needed in preparation for exascale computing and to enable high-fidelity materials research is robust reduced scaling, $\mathcal{O}$(N) or $\mathcal{O}$(N)log(N), DFT approaches. Thus, in agreement with the other co-PIs, the decision was made to purse QMC as a complimentary method on Mira since it would be valuable for the scientific community and could easily leverage the massive parallelism that would be provided by the Blue Gene/Q. We note that the two other atomic-scale methods mentioned here, CMD and DFT, are being pursued by other ESP projects on Mira. The remainder of this report will focus on our progress on QMC.
Benali, Anouar↗
Evolutionary optimization of PAW data-sets for accurate high pressure simulations
We examine the challenge of performing accurate electronic structure calculations at high pressures by comparing the results of all-electron full potential linearized augmented-plane-wave calculations, as implemented in the WIEN2k code, with those of the projector augmented wave (PAW) method, as implemented in Quantum ESPRESSO or Abinit code. In particular, we focus on developing an automated and consistent way of generating transferable PAW data-sets that can closely produce the all electron equation of state defined from zero to arbitrary high pressures. The technique we propose is an evolutionary search procedure that exploits the ATOMPAW code to generate atomic data-sets and the Quantum ESPRESSO software suite for total energy calculations. We demonstrate different aspects of its workability by optimizing PAW basis functions of some elements relatively abundant in planetary interiors. In addition, we introduce a new measure of atomic data-set goodness by considering their performance uniformity over an extended pressure range. - Highlights: • An attempt made to achieve the AE-FLAPW accuracy, at the favorable computational efficiency of the PAW approach for an outspread pressure range. • Provided an automatic and consistent route for developing efficient and transferable PAW data-sets. • Explicated an improved measure for the goodness of data-sets by including uniformity of the data-set performance over an extended pressure range.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Three new silver uranyl diphosphonates: Structures and properties
Explore the source record for details and available documents.
nuclear (including radiation effects), materials a↗
Energy band gap and spectroscopic studies in Mn{sub 1-x}Cu{sub x}WO{sub 4} (0 ≤ x ≤ 0.125)
A study on the effect of nonmagnetic Cu{sup 2+} substitution at Mn{sup 2+} site on the structural and energy band gap of the MnWO{sub 4} is reported. Convenient solid state reaction route has been adopted for the synthesis of Mn{sub 1-x}Cu{sub x}WO{sub 4}. X-ray diffraction (XRD) pattern showed high crystalline quality of the prepared samples. Raman spectroscopic studies were carried out to understand the structural aspects of the doping. 15 Raman active modes were identified out of 18, predicted for wolframite type monoclinic structure of MnWO{sub 4}. UV-visible diffuse reflectance spectra were recorded and analyzed to get energy band gap of the studied system and are found in the range of 2.5 eV to 2.04 eV with a systematic decrease with the increase in Cu{sup 2+} concentration. Energy band gap values are verified by Density Functional Theory calculations based on projector augmented wave (PAW) method. The calculated values are in good agreement with the experimental data.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Noncollinear ground states of solids with a source-free exchange correlation functional
In this paper, we expand upon the source-free (SF) exchange correlation (XC) functional developed by Sangeeta Sharma and coworkers to plane-wave density functional theory (DFT) based on the projector augmented wave (PAW) method. This constraint is implemented by the current authors within the VASP source code, using a fast Poisson solver that capitalizes on the parallel three-dimensional fast Fourier transforms (FFTs) implemented in VASP. Using this modified XC functional, we explore the improved convergence behavior that results from applying this constraint to the GGA-PBE+U+J functional. In the process, we compare the noncollinear magnetic ground state computed by each functional and their SF counterpart for a select number of magnetic materials in order to provide a metric for comparing with experimentally determined magnetic orderings. We observe significantly improved agreement with experimentally measured magnetic ground-state structures after applying the source-free constraint. Furthermore, we explore the importance of considering probability current densities in spin-polarized systems, even under no applied field. We analyze the XC torque as well, in order to provide theoretical and computational analyses of the net XC magnetic torque induced by the source-free constraint. Along these lines, we highlight the importance of properly considering the real-space integral of the source-free local magnetic XC field. Our analyses on probability currents, net torque, and constant terms draw additional links to the rich body of previous research on spin-current density functional theory (SCDFT), and pave the way for future extensions and corrections to the SF corrected XC functional.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Erratum: First-Principles Approach to Calculating Energy Level Alignment at Aqueous Semiconductor Interfaces [Phys. Rev. Lett. 113, 176802 (2014)]
We note, while extending our approach in our Letter to calculate energy level alignment at aqueous TiO 2 interfaces, we discovered two problems with our application of the GW methodology. Briefly, in our approach, we used a DFT + U approach for the underlying structure determination and as an input to the GW calculations. Then, GW calculations were carried out at the full-frequency spectrum-self-consistent level. All calculations were performed with the projector augmented wave (PAW) method as it is implemented in VASP. The specific semiconductors considered were GaN and ZnO. The first problem was caused by improper handling of the U parameter when initializing the follow-up GW calculations, leading to a final calculated band gap for GaN and ZnO that was too large. Second, we used an older set of pseudopotentials for the PAW method. In a paper that was published while our original manuscript was under review, the impact of the completeness of the partial waves representing empty states in the PAW method was analyzed, showing significant issues with the convergence of the GW results for certain materials, particularly ZnO. Specifically, those issues included use of PAW pseudopotentials that deviate too much from norm conservation, an unduly restricted basis set size that constrains the number of empty states available, and use of a small cutoff for the dielectric matrix size that can result in a false convergence with respect to the number of empty states included. Our GW calculations for GaN and ZnO were affected by these problems.
36 MATERIALS SCIENCE↗
Calculated electron paramagnetic resonance g-tensor and hyperfine parameters for zinc vacancy and N related defects in ZnO
Various defects in ZnO, focused on substitutional N O and N 2 in various sites, O-site, interstitial and Zn-site are studied using first-principles calculations with the goal of understanding the electron paramagnetic resonance (EPR) center reported for N 2 in ZnO and substitutional N on the O-site. The g tensors are calculated using the gauge including projector augmented wave (GIPAW) method and compared with experiments. The g-tensor of the free $N^{+}_{2}$ and $N^{–}_{2}$ radicals and their various contributions within the GIPAW theory are analyzed first to provide a baseline reference for the accuracy of the method and for understanding the N 2 behavior in ZnO. Previous controversies on the site location of N 2 in ZnO for this EPR center and on the shallow or deep nature and donor or acceptor nature of this center are resolved. Here, we find that the N 2 on the Zn site is mostly zinc-vacancy like in its spin density and g-tensor, while for the O-site, a model with the N 2 axis lying in-the basal plane and the singly occupied π g -orbital along the c axis provides good agreement with experiment. For the interstitial location, if the N 2 is not strongly interacting with the surroundings, no levels in the gap are found and hence also no possible EPR center. The calculated g-tensors for N O and V Zn are also found to be in good agreement with experiment. The effects of different functionals affecting the localization of the spin density are shown to affect the g-tensor values.
36 MATERIALS SCIENCE↗
Effect of Impurities on the Formation Energy of Point Defects in the γ-TiAl Alloy
The projector augmented-wave method within density functional theory is used to study the effect of 4d elements and Group IIIA and IVA elements on the energetics of point defect formation and elastic moduli of the γ-TiAl alloy. To calculate the point defect formation energy, the grand canonical formalism was used. It is shown that the formation energy of aluminum vacancies decreases by approximately 1.3 eV with increasing its content in the series of alloys Ti{sub 3}Al–TiAl–TiAl{sub 3}, whereas the formation energy of titanium vacancies changes insignificantly. In general, the formation of titanium vacancies in these alloys is preferred to the formation of aluminum vacancies. It has been found that Nb, Mo, Tc, Ru, Rh, and Pd impurities on the aluminum sublattice contribute to an increase in the formation energy of an aluminum vacancies, and Mo and Tc also lead to an increase in the formation energy of titanium vacancies. All 4d impurities, if they substitute for titanium, reduce the formation energy of aluminum vacancies, and Nb and Mo are favorable for increasing the formation energy of titanium vacancies. The influence of impurities on the chemical bond in the γ-TiAl alloy and its elastic moduli and characteristics based on them is discussed.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Effect of Impurities on the Oxygen Adsorption Properties on the NiTi(110) Surface
The effect of 3d–5d elements on the oxygen adsorption energy on the NiTi(110) surface has been studied by the projector augmented-waves method within density functional theory. It is shown that almost all elements, except for a few elements of the end of d periods, lead to an increase in the adsorption energy if they substitute for nickel. On the contrary, the substitutional impurities in the titanium sublattice lower this energy. Based on the analysis of the electronic characteristics of the surface with impurities, it has been found that an increase/decrease in the oxygen adsorption energy on NiTi(110) correlates with a change in the ionic contribution to the mechanism of oxygen bonding with the surface.
36 MATERIALS SCIENCE↗
Ab initio thermodynamic characteristics of the formation of oxygen vacancies, and boron, carbon, and nitrogen impurity centers in anatase
Using the ab initio projector augmented wave (PAW) method, calculations are performed for the electronic energy-band structure of titanium dioxide having the structure of anatase doped with boron, nitrogen, and carbon. Thermodynamic characteristics are determined for the formation of impurity centers, such as the preference energy for the interstitial position, the energy of impurity oxidation, and the energy of oxygen vacancy formation. It is shown that under the conditions of thermodynamic equilibrium the interstitial position of boron atoms is stable, whereas carbon atoms, depending on the oxygen pressure, can occupy both interstitial positions and substitutional positions of oxygen atoms, and nitrogen atoms replace oxygen atoms. It is shown that the presence of oxygen vacancies promotes the thermodynamic stability of carbon and nitrogen atoms. The obtained densities of electronic states correspond to ESR spectroscopy data, which indicates the presence of spin-polarized electrons in the states of the oxygen vacancy.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Adhesion at the Ta(Mo)/NiTi Interface
The atomic and electronic structures of Me(110)/NiTi(110) and TiO{sub 2}(100)/Me(110) interfaces (where Me = Ta, Mo) have been studied by the projector augmented-wave method in the framework of the density functional theory. It is established that the formation of intermediate oxide layers at the metal–oxide interface can lead to decreased adhesion at the interface.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗