Lithium adsorption on graphite from density functional theory calculations

The structural, energetic, and electronic properties of the Li/graphite system are studied through density functional theory (DFT) calculations using both the local spin density approximation (LSDA), and the gradient-corrected Perdew-Burke-Ernzerhof (PBE) approximation to the exchange-correlation energy. The calculations were performed using plane waves basis, and the electron-core interactions are described using pseudopotentials. We consider a disperse phase of the adsorbate comprising one Li atom for each 16 graphite surface cells, in a slab geometry. The close contact between the Li nucleus and the graphene plane results in a relatively large binding energy (larger than 1.1 eV). A detailed analysis of the electronic charge distribution, density difference distribution, and band structures indicates that one valence electron is entirely transferred from the atom to the surface, which gives rise to a strong interaction between the resulting lithium ion and the cloud of pi electrons in the substrate. We show that it is possible to explain the differences in the binding of Li, Na, and K adatoms on graphite considering the properties of the corresponding cation/aromatic complexes.     

Diminished gradient dependence of density functionals: constraint satisfaction and self-interaction correction

The Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation for the exchange-correlation energy functional has two nonempirical constructions, based on satisfaction of universal exact constraints on the hole density or on the energy. We show here that, by identifying one possible free parameter in exchange and a second in correlation, we can continue to satisfy these constraints while diminishing the gradient dependence almost to zero (i.e., almost recovering the local spin density approximation or LSDA). This points out the important role played by the Perdew-Wang 1991 nonempirical hole construction in shaping PBE and later constructions. Only the undiminished PBE is good for atoms and molecules, for reasons we present, but a somewhat diminished PBE could be useful for solids; in particular, the surface energies of solids could be improved. Even for atoms and molecules, a strongly diminished PBE works well when combined with a scaled-down self-interaction correction (although perhaps not significantly better than LSDA). This shows that the undiminished gradient dependence of PBE and related functionals works somewhat like a scaled-down self-interaction correction to LSDA.     

Density functional energy decomposition into one- and two-atom contributions

The present work provides a generalization of Mayer's energy decomposition for the density-functional theory (DFT) case. It is shown that one- and two-atom Hartree-Fock energy components in Mayer's approach can be represented as an action of a one-atom potential V(A) on a one-atom density rho(A) or rho(B). To treat the exchange-correlation term in the DFT energy expression in a similar way, the exchange-correlation energy density per electron is expanded into a linear combination of basis functions. Calculations carried out for a number of density functionals demonstrate that the DFT and Hartree-Fock two-atom energies agree to a reasonable extent with each other. The two-atom energies for strong covalent bonds are within the range of typical bond dissociation energies and are therefore a convenient computational tool for assessment of individual bond strength in polyatomic molecules. For nonspecific nonbonding interactions, the two-atom energies are low. They can be either repulsive or slightly attractive, but the DFT results more frequently yield small attractive values compared to the Hartree-Fock case. The hydrogen bond in the water dimer is calculated to be between the strong covalent and nonbonding interactions on the energy scale.     

Density functional study of structural and electronic properties of AlnN (1 ≤ n ≤ 12) clusters

Low‐lying equilibrium geometric structures of Al_nN (n = 1–12) clusters obtained by an all‐electron linear combination of atomic orbital approach, within spin‐polarized density functional theory, are reported. The binding energy, dissociation energy, and stability of these clusters are studied within the local spin density approximation (LSDA) and the three‐parameter hybrid generalized gradient approximation (GGA) due to Becke–Lee–Yang–Parr (B3LYP). Ionization potentials, electron affinities, hardness, and static dipole polarizabilities are calculated for the ground‐state structures within the GGA. It is observed that symmetric structures with the nitrogen atom occupying the internal position are lowest‐energy geometries. Generalized gradient approximation extends bond lengths as compared with the LSDA lengths. The odd–even oscillations in the dissociation energy, the second differences in energy, the highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO–LUMO) gaps, the ionization potential, the electron affinity, and the hardness are more pronounced within the GGA. The stability analysis based on the energies clearly shows the Al₇N cluster to be endowed with special stability. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Potentials in density functional theory and the importance of sum rules

In an ab initio approach to density functional theory one needs to know the electronic pair-density averaged over the coupling strength of the pair-interaction. As this pair-function is not available without having solved the associated N-electron problem, one has to resort to universal properties of the pair-density that are independent of specific features of the ground-state wavefunction. By exploiting these universal properties and so-called sum rules for the pair-correlation factors we derive very simple approximate spin-dependent expressions for the exchange-correlation energy per particle and for the associated potential in the Kohn-Sham equations. There is some similarity of the resulting density functionals with those obtained from the widely applied local spin density approximation (LSDA) based on electron gas theory. As the application of the latter to exceedingly inhomogeneous gases in realistic systems is very debatable, the manifest similarity seems to suggest that LSDA can consistently be justified only via the above pair-density analysis, but the justification of certain electron gas refinements may remain questionable. We shortly review similar attempts made by other authors and particularly focus on the issue of self-interaction and the “overbinding problem”. We demonstrate for the 3d- and 4d-metals that our approximation yields density of states (DOS), magnetic moments and Stoner parameters that are practically identical with respective data obtained from up-to-date LSDA- or gradient corrected (GGA-)potentials. There is also an excellent agreement of the DOS for the insulators (semi-conductors) C, Si, Ge, and GaAs. We show that our approach yields cohesive energies for these materials that are very close to the GGA-values indicating a distinct improvement over the standard LSDA-values. The calculations have been performed with the aid of the WIEN 97 computer code based on the Full Potential Linear Augmented Plane Wave (FLAPW-) method.     

Average excitation energies from time-dependent density functional response theory

The authors present an occupation number averaging scheme for time-dependent density functional response theory (TD-DFRT) in frequency domain. The known problem that TD-DFRT within the local (spin) density approximation (LDA/LSDA) inaccurately predicts Rydberg and charge-transfer excitation energies has been reexamined from the methodology of linear response, without explicit correction of the exchange-correlation potential. The working equations of TD-DFRT are adapted to treat arbitrary difference of orbital occupation numbers, using the nonsymmetric matrix form of Casida's formulation of TD-DFRT [M. E. Casida, in Recent Advances in Density Functional Methods, edited by D. P. Chong (World Scientific, Singapore, 1995), Pt. I, p. 155]. The authors' scheme is applied to typical closed-shell and open-shell molecular systems by examining the dependence of excitation energies on the fraction of excited electron. Good performance of this modified linear response scheme is shown, and is consistent with the authors' previous examination by the real-time propagation approach, suggesting that the calculation of average excitation energies might be one of the ways to better decode excitation energies from LDA/LSDA. Different techniques for treating singlet, triplet, and doublet states are discussed.     

有机-无机杂化化合物[Cu(μ-cbdca)(H2O)]n的电子结构及铁磁性

利用基于密度泛函理论(DFT)的第一性原理的FP_LAPW方法, 对以铜离子为磁性中心的化合物[Cu(μ-cbdca)(H2O)]n(cbdca=cyclobutanedicarboxylate)的电子结构及磁性质进行了计算. 对该材料的铁磁性、反铁磁性和非磁性三种状态下的总能量进行了计算. 计算结果表明, [Cu(μ-cbdca)(H2O)]n的铁磁态能量最低, 该化合物为稳定的铁磁性物质, 该结果与实验吻合较好. 对原子磁矩的计算结果发现, 铜原子对化合物磁性的贡献较大, 双齿配体上的氧原子和碳原子的贡献相对较小.     

Design of low band gap polymers employing density functional theory—hybrid functionals ameliorate band gap problem

Band gaps in solids and excitation energies in finite systems are underestimated significantly if estimated from differences between eigenvalues obtained within the local spin density approximation (LSDA). In this article we present results on 20 small- and medium-sized π-systems which show that HOMO–LUMO energy differences obtained with the B3LYP, B3P86, and B3PW91 functionals are in good agreement with vertical excitation energies from UV-absorption spectra. The improvement is a result of the use of the exact Hartree–Fock exchange with hybrid methods. Negative HOMO energies and negative LUMO energies do not provide good estimates for IPs and EAs. In contrast to Hartree–Fock theory, where IPs are approximated well and EAs are given poorly, DFT hybrid methods underestimate IPs and EAs by about the same amount. LSDA yields reasonable EAs but poor IPs. © 1997 John Wiley & Sons, Inc. J Comput Chem 18 : 1943–1953, 1997     

有机-无机杂化化合物[Cu（μ-cbdca）（H2O）]n的电子结构及铁磁性

利用基于密度泛函理论(DVF)的第一性原理的FP_LAPW方法,对以铜离子为磁性中心的化合物[Cu(μ-cbdca)(H2O)]n(cbdca=cyclobutanedicarboxylate)的电子结构及磁性质进行了计算.对该材料的铁磁性、反铁磁性和非磁性三种状态下的总能量进行了计算.计算结果表明,[Cu(μ-cbdca)(H2O)]n的铁磁态能量最低,该化合物为稳定的铁磁性物质,该结果与实验吻合较好.对原子磁矩的计算结果发现,铜原子对化合物磁性的贡献较大,双齿配体上的氧原子和碳原子的贡献相对较小.     

用GGA密度泛函及其长程、色散校正方法计算各类氢键的结合能

应用广义梯度近似(GGA) (PW91和PBE)、含动能密度的广义梯度近似(meta-GGA) (M06-L)、杂化泛函(hyper-GGA)(M06-2X、X3LYP和B3LYP)及其长程校正泛函LC-DFT(CAM-B3LYP、LC-ωPBE和ωB97X)和色散校正密度泛函(DFT-D)(ωB97X-D和B97-D),用多种基函数对15种不同强度的传统氢键和非传统氢键体系的结合能进行了系统的计算与分析.并与高精度的CCSD(T)/aug-cc-pVQZ结果比较发现:在上述各类泛函中,对于氢键结合能的计算M06-2X和ωB97X-D泛函较为精确与可靠,且没有必要使用过大的基函数,6-311++G(2d,2p)或aug-cc-pVDZ水平的基组就已足够,各类泛函所计算结合能的基组重叠误差(BSSE)均较小,除ωB97X和ωB97X-D外,其它9种泛函不经BSSE校正也能得到同样甚至更准确的结果.     

Legendre-transform functionals for spin-density-functional theory

We provide a rigorous proof that the Hohenberg-Kohn theorem holds for spin densities by extending Lieb's Legendre-transform formulation to spin densities. The resulting spin-density-functional theory resolves several troublesome issues. Most importantly, the present paper provides an explicit construction for the spin potentials at any point along the adiabatic connection curve, thus providing a formal basis for the use of exchange-correlation functionals of the spin density in the Kohn-Sham density-functional theory (DFT). The practical implications of this result for unrestricted Kohn-Sham DFT calculations is considered, and the existence of holes below the Fermi level is discussed. We argue that an orbital's energy tends to increase as its occupation number increases, which provides the basis for a computational algorithm for determining the occupation numbers in Kohn-Sham DFT and helps explain the origin of Hund's rules and holes below the Fermi level.     

Comparison of DFT methods for molecular orbital eigenvalue calculations

We report how closely the Kohn-Sham highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) eigenvalues of 11 density functional theory (DFT) functionals, respectively, correspond to the negative ionization potentials (-IPs) and electron affinities (EAs) of a test set of molecules. We also report how accurately the HOMO-LUMO gaps of these methods predict the lowest excitation energies using both time-independent and time-dependent DFT (TD-DFT). The 11 DFT functionals include the local spin density approximation (LSDA), five generalized gradient approximation (GGA) functionals, three hybrid GGA functionals, one hybrid functional, and one hybrid meta GGA functional. We find that the HOMO eigenvalues predicted by KMLYP, BH&HLYP, B3LYP, PW91, PBE, and BLYP predict the -IPs with average absolute errors of 0.73, 1.48, 3.10, 4.27, 4.33, and 4.41 eV, respectively. The LUMOs of all functionals fail to accurately predict the EAs. Although the GGA functionals inaccurately predict both the HOMO and LUMO eigenvalues, they predict the HOMO-LUMO gap relatively accurately (approximately 0.73 eV). On the other hand, the LUMO eigenvalues of the hybrid functionals fail to predict the EA to the extent that they include HF exchange, although increasing HF exchange improves the correspondence between the HOMO eigenvalue and -IP so that the HOMO-LUMO gaps are inaccurately predicted by hybrid DFT functionals. We find that TD-DFT with all functionals accurately predicts the HOMO-LUMO gaps. A linear correlation between the calculated HOMO eigenvalue and the experimental -IP and calculated HOMO-LUMO gap and experimental lowest excitation energy enables us to derive a simple correction formula.     

Intermolecular dispersion energies from time-dependent density functional theory

Non-expanded dispersion energies are calculated from time-dependent coupled-perturbed density functional theory (DFT) employing various non-hybrid and hybrid exchange-correlation potentials and suitable adiabatic local density approximations for the exchange-correlation kernel. Considering the dimer systems He₂, Ne₂, Ar₂, NeAr, NeHF, ArHF, (H₂)₂, (HF)₂, and (H₂O)₂ it is shown that the effects of intramonomer electron correlation on the dispersion energy are accurately reproduced with the PBE0AC exchange-correlation potential. In contrast, the uncoupled sum-over-states approximation yields inacceptable errors. These are mainly due to neglect of the Coulomb and exchange-correlation kernels and therefore, not substantially improved through an asymptotic correction of the exchange-correlation potential.     

A dual-level approach to density-functional theory

An efficient approximate scheme for density-functional theory (DFT) calculations, which eliminates the time-consuming self-consistent-field (SCF) procedure, is proposed using a dual-level DFT approach. In this approach, dual levels of basis sets and exchange-correlation functionals are adopted. The dual-level DFT approach is based on the idea that the total electron density in the ground state can be represented in terms of the density evaluated using the low-quality basis set and the low-cost exchange-correlation functional. Since the SCF procedure is avoided in the total energy evaluation, the dual-level DFT approach drastically reduces the computational cost. The applications of several dual-level DFT calculations to molecular systems show that our approach is more efficient than the self-consistent DFT approach with a moderate accuracy.     

Quantum-chemical study of intramolecular spin-forbidden rearrangements of the transition metal chelate complexes

An overview of quantum-chemical methods of studying the spin-forbidden reactions is presented. It was shown that HF and MP2 methods overestimate stabilization of the high-spin states whereas DFT methods with the use of LSDA and CGA functionals underestimate these states. Hybrid exchange-correlation density functionals B3P86 and B3LYP give the best agreement with the experimental data. The intramolecular spin-forbidden rearrangement of bis-chelate Ni(II) and Co(II) complexes was studied using the B3LYP/6-311++G (d,p) method.     

Spin densities in two-component relativistic density functional calculations: noncollinear versus collinear approach

With present day exchange-correlation functionals, accurate results in nonrelativistic open shell density functional calculations can only be obtained if one uses functionals that do not only depend on the electron density but also on the spin density. We consider the common case where such functionals are applied in relativistic density functional calculations. In scalar-relativistic calculations, the spin density can be defined conventionally, but if spin-orbit coupling is taken into account, spin is no longer a good quantum number and it is not clear what the "spin density" is. In many applications, a fixed quantization axis is used to define the spin density ("collinear approach"), but one can also use the length of the local spin magnetization vector without any reference to an external axis ("noncollinear approach"). These two possibilities are compared in this work both by formal analysis and numerical experiments. It is shown that the (nonrelativistic) exchange-correlation functional should be invariant with respect to rotations in spin space, and this only holds for the noncollinear approach. Total energies of open shell species are higher in the collinear approach because less exchange energy is assigned to a given Kohn-Sham reference function. More importantly, the collinear approach breaks rotational symmetry, that is, in molecular calculations one may find different energies for different orientations of the molecule. Data for the first ionization potentials of Tl, Pb, element 113, and element 114, and for the orientation dependence of the total energy of I+2 and PbF indicate that the error introduced by the collinear approximation is approximately 0.1 eV for valence ionization potentials, but can be much larger if highly ionized open shell states are considered. Rotational invariance is broken by the same amount. This clearly indicates that the collinear approach should not be used, as the full treatment is easily implemented and does not introduce much more computational effort.     

Force fields for metallic clusters and nanoparticles

Atomic force fields for simulating copper, silver, and gold clusters and nanoparticles are developed. Potential energy functions are obtained for both monatomic and binary metallic systems using an embedded atom method. Many cluster configurations of varying size and shape are used to constrain the parametrization for each system. Binding energies for these training clusters were computed using density functional theory (DFT) with the Perdew-Wang exchange-correlation functional in the generalized gradients approximation. Extensive testing shows that the many-body potentials are able to reproduce the DFT energies for most of the structures that were included in the training set. The force fields were used to calculate surface energies, bulk structures, and thermodynamic properties. The results are in good agreement with the DFT values and consistent with the available experimental data.     

Assessing the performance of density functional theory for the dynamic polarizabilities of amino acids: Treatment of correlation and role of exact exchange

The theoretical determination of electric response properties of the biological systems is a field where the application of density functional theory (DFT) appears to be quite promising. In this work, the performance of 41 density functional methods is evaluated in predicting dynamic polarizabilities of an experimental benchmark set of 20 proteinogenic amino acids. The behavior of a large number of density functionals, including various types of the local spin density approximation (LSDA), generalized gradient approximation (GGA), meta‐GGA (m‐GGA), hybrid‐GGA (h‐GGA), hybrid meta‐GGA (hm‐GGA), and range‐separated hybrid‐GGA (rsh‐GGA), has been assessed for the purpose. Analyzing the results of our DFT benchmarking, we found that these computationally economical methods show very diverse predictive capability and a careful selection of DFT functionals is very important in the polarizability calculations. Considering the role of exchange, correlation, dispersion and long‐range corrections, it turned out that in the LSDA class, SVWN3 gives better results than SPL and SVWN5 toward the reference values. Of the GGA methods, OPBE outperforms all other functionals. The M06‐L is the best method of m‐GGA class. The B3LYP and TPSSh are the best functionals of h‐GGA and hm‐GGA lineages, respectively. Finally, CAM‐B3LYP is the best method of rsh‐GGA functionals that predicts the most accurate polarizability for amino acids by a large margin with respect to others. Overall, the best performing functionals turn out to be hm‐GGAs TPSSh, TPSS1KCIS, M05, tau‐HCTHhyb, and h‐GGA B3LYP. Hopefully, the results of this investigation might provide the useful guidance to propose a new exchange‐correlation functional for calculating the optical properties of biomolecular materials. © 2013 Wiley Periodicals, Inc.