First-principles local density approximation + U and generalized gradient approximation + U study of plutonium oxides

The electronic structure and properties of PuO2 and Pu2O3 have been studied from first principles by the all-electron projector-augmented-wave method. The local density approximation+U and the generalized gradient approximation+U formalisms have been used to account for the strong on-site Coulomb repulsion among the localized Pu 5f electrons. We discuss how the properties of PuO2 and Pu2O3 are affected by the choice of U as well as the choice of exchange-correlation potential. Also, oxidation reaction of Pu2O3, leading to formation of PuO2, and its dependence on U and exchange-correlation potential have been studied. Our results show that by choosing an appropriate U, it is promising to correctly and consistently describe structural, electronic, and thermodynamic properties of PuO2 and Pu2O3, which enable the modeling of redox process involving Pu-based materials possible.     

Description of core excitations by time-dependent density functional theory with local density approximation, generalized gradient approximation, meta-generalized gradient approximation, and hybrid functionals

Time-dependent density functional theory (TDDFT) is employed to investigate exchange-correlation-functional dependence of the vertical core-excitation energies of several molecules including H, C, N, O, and F atoms. For the local density approximation (LDA), generalized gradient approximation (GGA), and meta-GGA, the calculated X1s-->pi* excitation energies (X = C, N, O, and F) are severely underestimated by more than 13 eV. On the other hand, time-dependent Hartree-Fock (TDHF) overestimates the excitation energies by more than 6 eV. The hybrid functionals perform better than pure TDDFT because HF exchange remedies the underestimation of pure TDDFT. Among these hybrid functionals, the Becke-Half-and-Half-Lee-Yang-Parr (BHHLYP) functional including 50% HF exchange provides the smallest error for core excitations. We have also discovered the systematic trend that the deviations of TDHF and TDDFT with the LDA, GGA, and meta-GGA functionals show a strong atom-dependence. Namely, their deviations become larger for heavier atoms, while the hybrid functionals are significantly less atom-dependent.     

Away from generalized gradient approximation: orbital-dependent exchange-correlation functionals

The local-density approximation of density functional theory (DFT) is remarkably accurate, for instance, for geometries and frequencies, and the generalized gradient approximations have also made bond energies quite reliable. Sometimes, however, one meets with failure in individual cases. One of the possible routes towards better functionals would be the incorporation of orbital dependence (which is an implicit density dependency) in the functionals. We discuss this approach both for energies and for response properties. One possibility is the use of the Hartree-Fock-type exchange energy expression as orbital-dependent functional. We will argue that in spite of the increasing popularity of this approach, it does not offer any advantage over Hartree-Fock for energies. We will advocate not to apply the separation of exchange and correlation, which is so ingrained in quantum chemistry, but to model both simultaneously. For response properties the energies and shapes of the virtual orbitals are crucial. We will discuss the benefits that Kohn-Sham potentials can offer which are derived from either an orbital-dependent energy functional, including the exact-exchange functional, or which can be obtained directly as orbital-dependent functional. We highlight the similarity of the Hartree-Fock and Kohn-Sham occupied orbitals and orbital energies, and the essentially different meanings the virtual orbitals and orbital energies have in these two models. We will show that these differences are beneficial for DFT in the case of localized excitations (in a small molecule or in a fragment), but are detrimental for charge-transfer excitations. Again, orbital dependency, in this case in the exchange-correlation kernel, offers a solution.     

DFT+U study of defects in bulk rutile TiO(2)

We present a systematic study of electronic gap states in defected titania using our implementation of the Hubbard-U approximation in the grid-based projector-augmented wave density functional theory code, GPAW. The defects considered are Ti interstitials, O vacancies, and H dopants in the rutile phase of bulk titanium dioxide. We find that by applying a sufficiently large value for the Hubbard-U parameter of the Ti 3d states, the excess electrons localize spatially at the Ti sites and appear as states in the band gap. At U=2.5?eV, the position in energy of these gap states are in fair agreement with the experimental observations. In calculations with several excess electrons and U=2.5?eV, all of these end up in gap states that are spatially localized around specific Ti atoms, thus effectively creating one Ti(3+) ion per excess electron. An important result of this investigation is that regardless of which structural defect is the origin of the gap states, at U=2.5?eV, these states are found to have their mean energies within a few hundredths of an eV from 0.94 eV below the conduction band minimum.     

Long-range and short-range Coulomb correlation effects as simulated by Hartree–Fock,local density approximation,and generalized gradient approximation exchange functionals

 Exchange functionals used in density functional theory (DFT) are generally considered to simulate long-range electron correlation effects. It is shown that these effects can be traced back to the self-interaction error (SIE) of approximate exchange functionals. An analysis of the SIE with the help of the exchange hole reveals that both short-range (dynamic) and long-range (nondynamic) electron correlation effects are simulated by DFT exchange where the local density approximation (LDA) accounts for stronger effects than the generalized gradient expansion (GGA). This is a result of the fact that the GGA exchange hole describes the exact exchange hole close to the reference electron more accurately than the LDA hole does. The LDA hole is more diffuse, thus leading to an underestimation of exchange and stronger SIE effects, where the magnitude of the SIE energy is primarily due to the contribution of the core orbitals. The GGA exchange hole is more compact, which leads to an exaggeration of exchange in the bond and the nonbonding region and negative SIE contributions. Partitioning of the SIE into intra-/interelectronic and individual orbital contributions makes it possible to test the performance of a given exchange functional in different regions of the molecule. It is shown that Hartree–Fock exchange always covers some long-range effects via interelectronic exchange while self-interaction-corrected DFT is lacking these effects. Received: 25 May 2002 / Accepted: 7 October 2002 / Published online: 21 January 2003 Correspondence to: E. Kraka e-mail: kraka@theoc.gn.se Acknowledgements. This work was supported financially by the Swedish Natural Science Research Council (NFR). Calculations were done on the supercomputers of Nationellt Superdatorcentrum (NSC), Link?ping, Sweden. The authors thank the NSC for a generous allotment of computer time.     

Density functional LCAO calculations for solids: Comparison among Hartree–Fock,DFT local density approximation,and DFT generalized gradient approximation structural properties

In this article, the results of a recently implemented DFT a posteriori and Kohn-Sham (KS ) linear combination of atomic orbital computational scheme for solids are presented. The equilibrium lattice parameters, bulk moduli, and lattice energies are calculated for eight crystallized systems. Local density approximation (LDA ) and generalized gradient approximation (GCA ) functionals and potentials are used. The maps of the Hartree-Fock (HF ) and Ks electronic densities and band structures are depicted. The KS results confirm the trend of the a posteriori scheme. Very good agreement between calculated and experimental lattice energies has been found for GGA potentials. © 1995 John Wiley & Sons, Inc.     

Avoiding singularity problems associated with meta-GGA (generalized gradient approximation) exchange and correlation functionals containing the kinetic energy density

Convergence problems of meta-GGA (generalized gradient approximation) XC (exchange and correlation) functionals containing a self-interaction correction term are traced back to a singularity of the latter that occurs at critical points of the electron density. This is demonstrated for the bond critical point of equilibrium and stretched H2. A simple remedy is suggested that cures meta-XC functionals such as VSXC, TPSS, M05, M06, and their derivatives without extra cost.     

First-principles structures and stabilities of AlN+ (N = 46-62) clusters

We present plausible candidates for the global minimum structures of Al(N)(+) (N = 46-62) cluster ions, determined by pseudopotential density functional theory static calculations under the spin-polarized generalized gradient approximation. Our calculations provide a first important step toward the rationalization of recent calorimetric experiments on the meltinglike transition of Al(N)(+). Most clusters with N > or = 48 clearly adopt fragments of the face-centered-cubic (fcc) crystalline lattice, although with significant distortions and a substantial proportion of defects in some cases. Another important driving force for stabilization comes from (111)-like surfaces, as the clusters often prefer to adopt less compact structures in order to keep the proportion of (100)-like surfaces at a minimum level. Al(46)(+) and Al(47)(+) adopt rather disordered structures instead. We find indications of enhanced stabilities for N = 51, 57, and 61 and of a substantial structural change between Al(55)(+) and Al(56)(+). These features correlate, albeit qualitatively, with the experimental observations.     

Generalized gradient approximation adjusted to transition metals properties: Key roles of exchange and local spin density

Perdew-Burke-Ernzerhof (PBE) and PBE adapted for solids (PBEsol) are exchange-correlation (xc) functionals widely used in density functional theory simulations. Their differences are the exchange, μ, and correlation, β, coefficients, causing PBEsol to lose the Local Spin Density (LSD) response. Here, the μ/β two-dimensional (2D) accuracy landscape is analyzed between PBE and PBEsol xc functional limits for 27 transition metal (TM) bulks, as well as for 81 TM surfaces. Several properties are analyzed, including the shortest interatomic distances, cohesive energies, and bulk moduli for TM bulks, and surface relaxation degree, surface energies, and work functions for TM surfaces. The exploration, comparing the accuracy degree with respect experimental values, reveals that the found xc minimum, called VV, being a PBE variant, represents an improvement of 5% in mean absolute percentage error terms, whereas this improvement reaches ~11% for VVsol, a xc resulting from the restoration of LSD response in PBEsol, and so regarded as its variant.     

First-principles study of the electronic structures of icosahedral TiN (N=13,19,43,55) clusters

We have studied the electronic structures of icosahedral Ti(N) clusters (N=13, 19, 43, and 55) by using a real-space first-principles cluster method with generalized gradient approximation for exchange-correlation potential. The hexagonal close-packed and fcc close-packed clusters have been studied additionally for comparisons. It is found that the icosahedral structures are the most stable ones except for Ti(43), where fcc close-packed structure is favorable in energy. We present and discuss the variation of bond length, the features of the highest occupied molecular orbitals and the lowest unoccupied molecular orbital, the evolution of density of states, and the magnetic moment in detail. The results are in good agreement with the predictions from the collision-induced dissociation and size-selected anion photoelectron spectroscopy experiments.     

Investigations of hydrogen-bonded systems: Local density approximation and gradient corrections

The applicability of the local density approximation (LDA ) and of corresponding gradient corrections (for the exchange and correlation energy) for the treatment of the hydrogen bond is investigated. As test systems, we consider the water dimer and the H₂O…?HX complexes (X = F, Cl, Br): Using an LCAO scheme, their equilibrium geometries and interaction energies are ?alculated and compared with experimental data and with other calculations. We obtain that the LDA gives the geometries in qualitative agreement with other data, whereas the energies are overestimated. The use of the gradient corrections (GC ) according to Becke and Perdew leads to a significant improvement of the geometry, and especially of the interaction energies. The calculations indicate further that LDA + GC should also be able to describe weaker intermolecular interactions than the usual hydrogen bond. Finally, a short discussion of the charge distribution and the dipole moments of the H₂O…?HX complexes is performed. © 1994 John Wiley & Sons, Inc.     

A density-functional study of the structures, binding energies and total spins of Ni-Fe clusters using nonlocal norm-conserving pseudopotentials and the generalized gradient approximation

We report ab initio calculations of the structures, binding energies, and total spins of the clusters Ni(13), Ni(19), Ni(23), Ni(26), Ni(12)Fe, Ni(11)Fe(2), Ni(18)Fe, Ni(17)Fe(2), Ni(22)Fe, Ni(20)Fe(3), and Ni(25)Fe using a density-functional method that employs linear combination of atomic orbitals as basis sets, nonlocal norm-conserving pseudopotentials, and the generalized gradient approximation for exchange and correlation. Our results show that the Fe-doped Ni clusters, which have icosahedral or polyicosahedral ground-state structures similar to those of the corresponding pure Ni clusters, are most stable with the Fe atoms occupying internal positions, as has also been inferred from experimental results on the adsorption of molecular nitrogen on the cluster surfaces. We also rule out the possibility that the experimentally observed difference between the (nonpolyicosahedral) configurations of N(2)-saturated Ni(26) and N(2)-saturated Ni(25)Fe be due to the influence of the Fe atom on the energy of the underlying metal cluster.     

A density-functional study of the possibility of noncollinear magnetism in small Mn clusters using SIESTA and the generalized gradient approximation to exchange and correlation

We investigated the possibility of noncollinear magnetism in small Mn(n) clusters (n=2-6) using the density-functional method SIESTA with the generalized gradient approximation (GGA) to exchange and correlation. The lowest-energy states identified were collinear, with the atomic spin magnetic moments pointing in the same direction, for Mn(2) and Mn(3), and noncollinear for Mn(4), Mn(5) and, most decidedly, Mn(6). These SIESTA/GGA results, which are compared with those of an earlier SIESTA study that used the local spin density approximation, are qualitatively in keeping with the result obtained by VASP/GGA calculations.     

Communication: Effect of the orbital-overlap dependence in the meta generalized gradient approximation

We study for the first time the effect of the dependence of meta generalized gradient approximation (MGGA) for the exchange-correlation energy on its input, the kinetic energy density, through the dimensionless inhomogeneity parameter, α, that characterizes the extent of orbital overlap. This leads to a simple MGGA exchange functional, which interpolates between the single-orbital regime, where α = 0, and the slowly varying density regime, where α ≈ 1, and then extrapolates to α → ∞. When combined with a variant of the Perdew-Burke-Ernzerhof GGA correlation, the resulting MGGA performs equally well for atoms, molecules, surfaces, and solids.     

NO adsorption on MoS(x) clusters: a density functional theory study

The density functional theory (DFT) method has been used to investigate NO probe molecule adsorption on the stoichiometric (Mo(16)S(32)) and nonstoichiometric (Mo(16)S(34) and Mo(16)S(29)) clusters. The calculated adsorption energies indicate that the stoichiometric cluster has stronger NO affinity than the nonstoichiometric surfaces. It is also found that mononitrosyl adsorption is favored at low NO coverage, while dinitrosyl (germinal) and (NO)(2) dimer adsorption at high NO coverage are possible. Strong repulsive interaction has been found for the adsorbed dinitrosyl and (NO)(2) dimer species. In addition, the computed NO stretching frequencies for the mononitrosyl and dinitrosyl species agree well with the experimental data, while those of the dimer species are much lower than the suggested experimental data.     

Microhydration of Cs+ ion: a density functional theory study on Cs+-(H2O)n clusters (n=1-10)

Structure, energy enthalpy, and IR frequency of hydrated cesium ion clusters, Cs+-(H2O)n (n=1-10), are reported based on all electron calculations. Calculations have been carried out with a hybrid density functional, namely, Becke's three-parameter nonlocal hybrid exchange-correlation functional B3LYP applying cc-PVDZ correlated basis function for H and O atoms and a split valence 3-21G basis function for Cs atom. Geometry optimizations for all the cesium ion-water clusters have been carried out with several possible initial guess structures following Newton-Raphson procedure leading to many conformers close in energy. The calculated values of binding enthalpy obtained from present density functional based all electron calculations are in good agreement with the available measured data. Binding enthalpy profile of the hydrated clusters shows a saturation behavior indicating geometrical shell closing in hydrated structure. Significant shifts of O-H stretching bands with respect to free water molecule in IR spectra of hydrated clusters are observed in all the hydrated clusters.     

Sum Rules for the spin surface response of metal clusters within the time dependent local spin density approximation

We study surface spin collective modes of simple metallic clusters using a sum rule approach. We derive analytical expressions for the energy and cubic energy-weighted moments of the time-dependent local-spin-density approximation (TDLSDA) strength function. With these two moments we obtain a mean excitation energy of multipole spin modes and study its dispersion with cluster size up to very large sizes. The crucial roles played by the surface diffuseness and by correlation and kinetic energy effects on the resonance energies are stressed.     

Electronic structure of spheroidal metal containing carbon shells: study of the LaC60 and C60 clusters and their ions within the local density approximation

Ionization thresholds and electron affinities have been calculated within the Local Density Approximation (LDA) for the neutral, positively and negatively charged clusters of LaC₆₀ and C₆₀. The evaluated energies are found to be in good agreement with available experimental data. More accurate measurements are however necessary to verify the suggested spheroidal cage structure for these molecules.