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↗
Application of Ab Initio Methods in the Development of Advanced Technical Ceramics
A significant need exists to develop materials not only capable of providing desired electronic and mechanical properties but also survival at extreme temperatures during service and device fabrication, such as co-firing. Computational methods offers an efficient and systematic manner to design new materials and guide their development. As an example computational -based material approaches can be used to determine the suitability of a given materials as a practical thermoelectric for energy harvesting. In this presentation examples are given of applications to relevant technical ceramics such as thermoelectrics, dielectrics, and magnetic systems. The calculation were based on density functional theory and carried out with norm conserving and projector augmented wave (PAW) methods using commercial codes Materials Studio ( Biovia, Inc) with the Cambridge Serial Total Energy Package (CASTEP) and MedeA (Materials Design Inc.) utilizing the Vienna Ab-initio Simulation Package (VASP) as the respective computational engines. This study makes predictions of relevant technical properties of ceramic materials.
thermoelectric↗
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↗