Search NASA⌕ Search

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.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Time-dependent density functional theory with the orthogonal projector augmented wave method

The projector augmented wave (PAW) method of Blöchl linearly maps smooth pseudo wavefunctions to the highly oscillatory all-electron DFT orbitals. Compared to norm-conserving pseudopotentials (NCPP), PAW has the advantage of lower kinetic energy cutoffs and larger grid spacing at the cost of having to solve for non-orthogonal wavefunctions. We earlier developed orthogonal PAW (OPAW) to allow the use of PAW when orthogonal wavefunctions are required. In OPAW, the pseudo wavefunctions are transformed through the efficient application of powers of the PAW overlap operator with essentially no extra cost compared to NCPP methods. Previously, we applied OPAW to DFT. Here, we take the first step to make OPAW viable for post-DFT methods by implementing it in real-time time-dependent (TD) DFT. Using fourth-order Runge–Kutta for the time-propagation, we compare calculations of absorption spectra for various organic and biological molecules and show that very large grid spacings are sufficient, 0.6–0.7 bohr in OPAW-TDDFT rather than the 0.4–0.5 bohr used in traditional NCPP-TDDFT calculations. This reduces the memory and propagation costs by around a factor of 3. Our method would be directly applicable to any post-DFT methods that require time-dependent propagations such as the GW approximation and the Bethe–Salpeter equation.

Chemistry↗
DOE OSTI2015

Implementation of density functional embedding theory within the projector-augmented-wave method and applications to semiconductor defect states

We report a new implementation of the density functional embedding theory (DFET) in the VASP code, using the projector-augmented-wave (PAW) formalism. Newly developed algorithms allow us to efficiently perform optimized effective potential optimizations within PAW. The new algorithm generates robust and physically correct embedding potentials, as we verified using several test systems including a covalently bound molecule, a metal surface, and bulk semiconductors. We show that with the resulting embedding potential, embedded cluster models can reproduce the electronic structure of point defects in bulk semiconductors, thereby demonstrating the validity of DFET in semiconductors for the first time. Compared to our previous version, the new implementation of DFET within VASP affords use of all features of VASP (e.g., a systematic PAW library, a wide selection of functionals, a more flexible choice of U correction formalisms, and faster computational speed) with DFET. Furthermore, our results are fairly robust with respect to both plane-wave and Gaussian type orbital basis sets in the embedded cluster calculations. This suggests that the density functional embedding method is potentially an accurate and efficient way to study properties of isolated defects in semiconductors.

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

Orbital-free density functional theory implementation with the projector augmented-wave method

We present a computational scheme for orbital-free density functional theory (OFDFT) that simultaneously provides access to all-electron values and preserves the OFDFT linear scaling as a function of the system size. Using the projector augmented-wave method (PAW) in combination with real-space methods, we overcome some obstacles faced by other available implementation schemes. Specifically, the advantages of using the PAW method are twofold. First, PAW reproduces all-electron values offering freedom in adjusting the convergence parameters and the atomic setups allow tuning the numerical accuracy per element. Second, PAW can provide a solution to some of the convergence problems exhibited in other OFDFT implementations based on Kohn-Sham (KS) codes. Using PAW and real-space methods, our orbital-free results agree with the reference all-electron values with a mean absolute error of 10 meV and the number of iterations required by the self-consistent cycle is comparable to the KS method. The comparison of all-electron and pseudopotential bulk modulus and lattice constant reveal an enormous difference, demonstrating that in order to assess the performance of OFDFT functionals it is necessary to use implementations that obtain all-electron values. The proposed combination of methods is the most promising route currently available. We finally show that a parametrized kinetic energy functional can give lattice constants and bulk moduli comparable in accuracy to those obtained by the KS PBE method, exemplified with the case of diamond.

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

GPAW: An open Python package for electronic structure calculations

We review the GPAW open-source Python package for electronic structure calculations. GPAW is based on the projector-augmented wave method and can solve the self-consistent density functional theory (DFT) equations using three different wave-function representations, namely real-space grids, plane waves, and numerical atomic orbitals. The three representations are complementary and mutually independent and can be connected by transformations via the real-space grid. This multi-basis feature renders GPAW highly versatile and unique among similar codes. By virtue of its modular structure, the GPAW code constitutes an ideal platform for the implementation of new features and methodologies. Moreover, it is well integrated with the Atomic Simulation Environment (ASE), providing a flexible and dynamic user interface. In addition to ground-state DFT calculations, GPAW supports many-body GW band structures, optical excitations from the Bethe–Salpeter Equation, variational calculations of excited states in molecules and solids via direct optimization, and real-time propagation of the Kohn–Sham equations within time-dependent DFT. A range of more advanced methods to describe magnetic excitations and non-collinear magnetism in solids are also now available. In addition, GPAW can calculate non-linear optical tensors of solids, charged crystal point defects, and much more. Recently, support for graphics processing unit (GPU) acceleration has been achieved with minor modifications to the GPAW code thanks to the CuPy library. We end the review with an outlook, describing some future plans for GPAW.

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

SU-C-204-03: DFT Calculations of the Stability of DOTA-Based-Radiopharmaceuticals

Purpose: Application of the density function theory (DFT) to investigate the structural stability of complexes applied in cancer therapy consisting of the 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA) chelated to Ac225, Fr221, At217, Bi213, and Gd68 radio-nuclei. Methods: The possibility to deliver a toxic payload directly to tumor cells is a highly desirable aim in targeted alpha particle therapy. The estimation of bond stability between radioactive atoms and the DOTA chelating agent is the key element in understanding the foundations of this delivery process. Thus, we adapted the Vienna Ab-initio Simulation Package (VASP) with the projector-augmented wave method and a plane-wave basis set in order to study the stability and electronic properties of DOTA ligand chelated to radioactive isotopes. In order to count for the relativistic effect of radioactive isotopes we included Spin-Orbit Coupling (SOC) in the DFT calculations. Five DOTA complex structures were represented as unit cells, each containing 58 atoms. The energy optimization was performed for all structures prior to calculations of electronic properties. Binding energies, electron localization functions as well as bond lengths between atoms were estimated. Results: Calculated binding energies for DOTA-radioactive atom systems were −17.792, −5.784, −8.872, −13.305, −18.467 eV for Ac, Fr, At, Bi and Gd complexes respectively. The displacements of isotopes in DOTA cages were estimated from the variations in bond lengths, which were within 2.32–3.75 angstroms. The detailed representation of chemical bonding in all complexes was obtained with the Electron Localization Function (ELF). Conclusion: DOTA-Gd, DOTA-Ac and DOTA-Bi were the most stable structures in the group. Inclusion of SOC had a significant role in the improvement of DFT calculation accuracy for heavy radioactive atoms. Our approach is found to be proper for the investigation of structures with DOTA-based-radiopharmaceuticals and will enhance our understanding of processes occurring at subatomic levels.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Local structure study of Fe dopants in Ni-deficit Ni 3 Al alloys

We found that the local electronic and magnetic structure, hyperfine interactions, and phase composition of polycrystalline Ni–deficient Ni 3-x FexAl (x = 0.18 and 0.36) were investigated by means of 57 Fe Mössbauer spectroscopy. The samples were characterized by X–ray diffraction and magnetization measurements. The ab initio calculations performed with the projector augmented wave method and the calculations of the energies of iron point defects were done to elucidate the electronic structure and site preference of Fe doped Ni 3 Al. Moreover, the value of calculated electric field gradient tensor V zz =1.6 10 21 Vm -2 matches well with the results of Mössbauer spectroscopy and indicates that the Fe atoms occupy Ni sites.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2012

Crystal structure and magnetic properties of '{alpha} Prime Prime -Fe{sub 16}N{sub 2}' containing residual {alpha}-Fe prepared by low-temperature ammonia nitridation

Slight enhancement of saturation magnetization to 219 A m{sup 2} kg{sup -1} was observed from 199 A m{sup 2} kg{sup -1} for the original {alpha}-Fe on the intermediate nitrided mixture of '{alpha} Prime Prime -Fe{sub 16}N{sub 2}' with residual {alpha}-Fe among the low temperature ammonia nitridation products under 5 T magnetic field at room temperature. The value changed not linearly against the yield as had been expected. Crystal structure refinement indicated that the phase similar to {alpha} Prime Prime -Fe{sub 16}N{sub 2} had deviations on its lattice constants and positional parameters, compared to previously reported values for {alpha} Prime Prime -Fe{sub 16}N{sub 2}. Spin-polarized total energy calculations were performed using the projector-augmented wave method as implemented in the Vienna ab-initio simulation package (VASP) to calculate magnetic moment on the refined crystal structure of the intermediate '{alpha} Prime Prime -Fe{sub 16}N{sub 2}'. The calculations supported the observed magnetization enhancement in the intermediate nitridation product. - Graphical abstract: Crystal structural parameters slightly change in the intermediate nitrided '{alpha} Prime Prime -Fe{sub 16}N{sub 2}' from those in {alpha} Prime Prime -Fe{sub 16}N{sub 2} to show the magnetization maxima in the mixture of '{alpha} Prime Prime -Fe{sub 16}N{sub 2}' and the residual {alpha}-F. Highlights: Black-Right-Pointing-Pointer Larger magnetization was observed than the value of Fe{sub 16}N{sub 2} on its intermediate nitrided mixture with residual {alpha}-Fe. Black-Right-Pointing-Pointer The enhancement was related to the crystal structural deviation from Fe{sub 16}N{sub 2} on the intermediate nitride. Black-Right-Pointing-Pointer It was supported by spin-polarized total energy calculation using the deviated structure.

36 MATERIALS SCIENCE↗
DOE OSTI2015

High pressure structural behavior of YGa{sub 2}: A combined experimental and theoretical study

High pressure structural stability studies were carried out on YGa{sub 2} (AlB{sub 2} type structure at NTP, space group P6/mmm) up to a pressure of ~35 GPa using both laboratory based rotating anode and synchrotron X-ray sources. An isostructural transition with reduced c/a ratio, was observed at ~6 GPa and above ~17.5 GPa, the compound transformed to orthorhombic structure. Bulk modulus B{sub 0} for the parent and high pressure phases were estimated using Birch–Murnaghan and modified Birch–Murnaghan equation of state. Electronic structure calculations based on projector augmented wave method confirms the experimentally observed two high pressure structural transitions. The calculations also reveal that the ‘Ga’ networks remains as two dimensional in the high pressure isostructural phase, whereas the orthorhombic phase involves three dimensional networks of ‘Ga’ atoms interconnected by strong covalent bonds. - Graphical abstract: High pressure X-ray diffraction patterns of YGa{sub 2} up to ~35 GPa shows an isostructural phase transition at ~5 GPa and transition to an orthorhombic structure ~14 GPa. - Highlights: • High pressure structural stability studies were carried out on YGa{sub 2} up to 35 GPa. • An isostructural transition with reduced c/a ratio was observed above 6 GPa. • Above 17.5 GPa, the compound transformed to orthorhombic structure. • PAW based electronic structure calculations have been carried out. • Calculations confirm the experimentally observed structural transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Effect of van der Waals interactions on the structural and binding properties of GaSe

The influence of van der Waals interactions on the lattice parameters, band structure, elastic moduli and binding energy of layered GaSe compound has been studied using projector-augmented wave method within density functional theory. We employed the conventional local/semilocal exchange-correlation functionals and recently developed van der Waals functionals which are able to describe dispersion forces. It is found that application of van der Waals density functionals allows to substantially increase the accuracy of calculations of the lattice constants a and c and interlayer distance in GaSe at ambient conditions and under hydrostatic pressure. The pressure dependences of the a-parameter, Ga–Ga, Ga–Se bond lengths and Ga–Ga–Se bond angle are characterized by a relatively low curvature, while c(p) has a distinct downward bowing due to nonlinear shrinking of the interlayer spacing. From the calculated binding energy curves we deduce the interlayer binding energy of GaSe, which is found to be in the range 0.172–0.197 eV/layer (14.2–16.2 meV/Å{sup 2}). - Highlights: • Effects of van der Waals interactions are analyzed using advanced density functionals. • Calculations with vdW-corrected functionals closely agree with experiment. • Interlayer binding energy of GaSe is estimated to be 14.2–16.2 meV/Å{sup 2}.

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

Gold in the layered structures of R 3 Au 7 Sn 3 : From relativity to versatility

A new isotypic series of ternary rare earth element-gold-tetrel intermetallic compounds has been synthesized and their structures and properties have been characterized. R 3 Au 7 Sn 3 (R = Y, La-Nd, Sm, Gd-Tm, Lu) crystallize with the hexagonal Gd 3 Au 7 Sn 3 prototype (Pearson symbol hP26; P6 3/m , a = 8.110-8.372 Å, c = 9.351-9.609 Å, V cell = 532.7-583.3 Å3, Z = 2), an ordered variant of the Cu 10 Sn 3 -type. Their structure is built up by GdPt 2 Sn-type layers, which feature edge-sharing Sn@Au 6 trigonal antiprisms connected by trigonal R3 groups. Additional insertion of gold atoms leads to the formation of new homoatomic Au clusters, Au@Au 6 ; alternatively, the structure can be considered as a superstructural polyhedral packing of the ZrBeSi-type. The magnetization, heat ca-pacity and electrical resistivity have been measured for R 3 Au 7 Sn 3 (R = Ce, Pr, Nd and Tb). All four compounds order antiferromagnetically with the highest T N of 13 K for Tb 3 Au 7 Sn 3 . In Ce 3 Au 7 Sn 3 , which has a T N of 2.9 K, the heat capacity and electrical resistivity data in zero and applied fields indicate the presence of Kondo interactions. The coefficient of the linear term in the electronic heat capacity, γ, derived from the heat capacity data below 0.5 K is 211 mJ/Ce mol K 2 suggesting strong electronic correlations due to the Kondo interaction. The electronic structure calculations based on the projector augmented wave method for particular representatives of the series suggest different tendencies of the localized R-4f AOs to hybridize with the valence states. LMTO-based bonding analysis on the non-magnetic La 3 Au 7 Sn 3 indicates that the integrated crystal orbital Hamilton popu-lations (COHPs) are dominated by the heteroatomic Au–Sn contacts; however, contributions from La–Au and La–Sn separations are significant, both together exceeding 40 % in the overall bonding. Furthermore, homoatomic Au–Au interactions are evident for the Au@Au 6 units but, despite of the high atomic concentration of Au in the compound, they do not dominate the entire bonding picture.

36 MATERIALS SCIENCE↗
DOE OSTI2015

From the ternary Eu(Au/In) 2 and EuAu 4 (Au/In) 2 with remarkable Au/In distributions to a new structure type: The gold-rich Eu 5 Au 16 (Au/In) 6 structure

The ternary Eu(Au/In) 2 (EuAu 0.46 In 1.54 (2) ) (I), EuAu 4 (Au/In) 2 (EuAu 4+x In 2–x with x = 0.75(2) (II), 0.93(2), and 1.03(2)), and Eu 5 Au 16 (Au/In) 6 (Eu 5 Au 17.29 In 4.71(3) ) (III) have been synthesized, and their structures were characterized by single-crystal X-ray diffraction. I and II crystallize with the CeCu 2 -type (Pearson Symbol oI 12; Imma ; Z = 4; a = 4.9018(4) Å; b = 7.8237(5) Å; c = 8.4457(5) Å) and the YbAl 4 Mo 2 -type (tI14; I4/ mmm ; Z = 2; a = 7.1612(7) Å; c = 5.5268(7) Å) and exhibit significant Au/In disorder. I is composed of an Au/In-mixed diamond-related host lattice encapsulating Eu atoms, while the structure of II features ribbons of distorted, squared Au 8 prisms enclosing Eu, Au, and In atoms. Combination of these structural motifs leads to a new structure type as observed for Eu 5 Au 16 (Au/In) 6 (Eu 5 Au 17.29 In 4.71(3) ) (oS108; Cmcm ; Z = 4; a = 7.2283(4) Å; b = 9.0499(6) Å; c = 34.619(2) Å), which formally represents a one-dimensional intergrowth of the series EuAu 2 –“EuAu 4 In 2 ”. The site preferences of the disordered Au/In positions in II were investigated for different hypothetical “EuAu 4 (Au/In) 2 ” models using the projector-augmented wave method and indicate that these structures attempt to optimize the frequencies of the heteroatomic Au–In contacts. Furthermore, a chemical bonding analysis on two “EuAu 5 In” and “EuAu 4 In 2 ” models employed the TB-LMTO-ASA method and reveals that the subtle interplay between the local atomic environments and the bond energies determines the structural and site preferences for these systems.

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

Transition metal dimer on Au(111) surface: A first principle study

The adsorption behaviour of transition metal dimers M{sub 2} (M= Cu, Ag, Au) on the Au(111) surface have been studied using the density functional theory formalism. The projector augmented wave method under the spin polarized version of generalized gradient approximation scheme was employed to calculate the total energy. The results suggest that all dimers prefer to orient in parallel to the surface plane, where two M atoms are adsorbed on two nearby threefold fcc sites. We have investigated the chemical interaction between M atoms and Au surface through electronic density of state analysis. It is found that on deposition, the electronic density of states (EDOS) of M{sub 2} dimer becomes broader in comparison to their gas phase spectrum.

36 MATERIALS SCIENCE↗
DOE OSTI2014

NMR characterization of hydrocarbon adsorption on calcite surfaces: A first principles study

The electronic and coordination environment of minerals surfaces, as calcite, are very difficult to characterize experimentally. This is mainly due to the fact that there are relatively few spectroscopic techniques able to detect Ca{sup 2+}. Since calcite is a major constituent of sedimentary rocks in oil reservoir, a more detailed characterization of the interaction between hydrocarbon molecules and mineral surfaces is highly desirable. Here we perform a first principles study on the adsorption of hydrocarbon molecules on calcite surface (CaCO{sub 3} (101{sup ¯}4)). The simulations were based on Density Functional Theory with Solid State Nuclear Magnetic Resonance (SS-NMR) calculations. The Gauge-Including Projector Augmented Wave method was used to compute mainly SS-NMR parameters for {sup 43}Ca, {sup 13}C, and {sup 17}O in calcite surface. It was possible to assign the peaks in the theoretical NMR spectra for all structures studied. Besides showing different chemical shifts for atoms located on different environments (bulk and surface) for calcite, the results also display changes on the chemical shift, mainly for Ca sites, when the hydrocarbon molecules are present. Even though the interaction of the benzene molecule with the calcite surface is weak, there is a clearly distinguishable displacement of the signal of the Ca sites over which the hydrocarbon molecule is located. A similar effect is also observed for hexane adsorption. Through NMR spectroscopy, we show that aromatic and alkane hydrocarbon molecules adsorbed on carbonate surfaces can be differentiated.

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

Dirac cones in transition metal doped boron nitride

The transition metal (TM) doped zinc blende boron nitride (c-BN) is studied by using the first principle calculation. TM atoms fill in the interstitials in c-BN and form two-dimensional honeycomb lattice. The generalized gradient approximation and projector augmented wave method are used. The calculated density of states and band structures show that d electrons of TM atoms form impurity bands in the gap of c-BN. When the TM-BN system is in ferromagnetic or non-magnetic state, Dirac cones emerge at the K point in Brillouin zone. When TM is Ti and Co, the Dirac cones are spin polarized and very close to the Fermi level, which makes them promising candidates of Dirac half-metal [H. Ishizuka and Y. Motome, Phys. Rev. Lett. 109, 237207 (2012)]. While TM is Ni and Cu, the system is non-magnetic and Dirac cones located above the Fermi level.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗