Grid-free density functional calculations on periodic systems

Density fitting scheme is applied to the exchange part of the Kohn-Sham potential matrix in a grid-free local density approximation for infinite systems with translational periodicity. It is shown that within this approach the computational demands for the exchange part scale in the same way as for the Coulomb part. The efficiency of the scheme is demonstrated on a model infinite polymer chain. For simplicity, the implementation with Dirac-Slater Xalpha exchange functional is presented only. Several choices of auxiliary basis set expansion coefficients were tested with both Coulomb and overlap metric. Their effectiveness is discussed also in terms of robustness and norm preservation.     

Relativistic all-electron two-component self-consistent density functional calculations including one-electron scalar and spin-orbit effects

We have implemented a Gaussian basis-set two-component self-consistent field method based on the fourth-order nuclear-only Douglas-Kroll-Hess approximation. Two-electron spin-orbit effects are included using Boettger's screened-nuclear spin-orbit approximation. In our two-component approach, the spin-orbit interaction is taken into account in a variational fashion employing a generalized Kohm-Sham scheme which allows one to work with hybrid density functionals. For open-shell systems we adopt the noncollinear spin-density approximation. Results are presented for equilibrium bond lengths, harmonic vibrational frequencies, and bond dissociation energies with local spin-density, generalized gradient approximation, and hybrid functionals in a set of benchmark molecules.     

Accurate density functional calculations on large systems

The theoretical underpinnings of the linear combination of Gaussian-type orbitals (LCGTO ) calculations of the density functional (DF ) energy of molecules and clusters are described. The generating function for three-center integrals of arbitrary angular momenta is given in the solid-spherical-harmonic basis. Variational fitting is described and its accuracy tested. The LCGTO-DF method is used to address questions related to the problem of how it is that the methods of cluster science, i.e., high-energy beams or currents, can be used to make C₆₀ in bulk quantities. In particular, it is shown that C₆₀ is neither especially stable nor is it the only large, stable, perfectly round, approximately sp² carbon molecule. © 1996 John Wiley & Sons, Inc.     

Analytical energy gradient evaluation in relativistic and nonrelativistic density functional calculations

The expressions of analytical energy gradients in density functional theory and their implementation in programs are reported. The evaluation of analytical energy gradients can be carried out in the fully 4-component relativistic, approximate relativistic, and nonrelativistic density functional calculations under local density approximation or general gradient approximation with or without frozen core approximation using different basis sets in our programs. The translational invariance condition and the fact that the one-center terms do not contribute to the energy gradients are utilized to improve the calculation accuracy and to reduce the computational effort. The calculated results of energy gradients and optimized geometry as well as atomization energies of some molecules by the analytical gradient method are in very good agreement with results obtained by the numerical derivative 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.     

Development of a new density functional program for all-electron calculation of proteins

In this article, we propose a new molecular orbital program for all-electron calculation of proteins which is based on density functional theory. To carry it out in a fully analytical way, we adopted the (pure-) analytical Xα method and modified it for saving a lot of memories for large-scale calculations. The recent software technology sophisticated in information science is inevitably applied to achieve calculations of large molecular systems. The program is coded by the object-oriented language C + +, its output is shown graphically, and the most of the procedures in this program are controlled through an efficient graphical user interface developed by ourselves. Such technology supports the safe construction of the huge software, the tidy representation of enormous data, and the ready control of complex calculations. Test calculations with various sizes of glycine polypeptides indicate that the computation time is proportional to the 1.7 powers of the number of residues. This result suggests that the all-electron calculations of proteins consisting of over 1000 atoms could be performed with distributed and/or massively parallel computers. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 245–256, 1997     

The heat of formation of the uranyl dication: theoretical evaluation based on relativistic density functional calculations

By using a set of model reactions, we estimated the heat of formation of gaseous UO2(2+) from quantum-chemical reaction enthalpies and experimental heats of formation of reference species. For this purpose, we performed relativistic density functional calculations for the molecules UO2(2+), UO2, UF6, and UF5. We used two gradient-corrected exchange-correlation functionals (revised Perdew-Burke-Ernzerhof (PBEN) and Becke-Perdew (BP)) and we accounted for spin-orbit interaction in a self-consistent fashion. Indeed, spin-orbit interaction notably affects the energies of the model reactions, especially if compounds of U(IV) are involved. Our resulting theoretical estimates for delta fH(o)0 (UO2(2+)), 365+/-10 kcal mol(-1) (PBEN) and 370+/-12 kcal mol(-1) (BP), are in quantitative agreement with a recent experimental result, 364+/-15 kcal mol(-1). Agreement between the results of the two different exchange-correlation functionals PBEN and BP supports the reliability of our approach. The procedure applied offers a general means to derive unknown enthalpies of formation of actinide species based on the available well-established data for other compounds of the element in question.     

Fractional occupation numbers and spin density functional calculations of degenerate systems

We implemented ab initio self‐consistent field (SCF) fractional occupation numbers (FON) calculation with Dunlap's interpolation scheme for the twisted ethylene, which is a prototype molecule of a σ–π biradical system. The calculational results are compared with those of complete‐active‐space (CAS) SCF and spin‐unrestricted Kohn–Sham (UKS) calculations on potential surfaces, occupation numbers of natural orbitals, and correlation entropies. It was found that the UKS methods gave similar results to CASSCF, while the FON solutions appeared in only the nearly complete degenerate region. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 93: 317–323, 2003     

Doping of polyaniline by acid-base chemistry: density functional calculations with periodic boundary conditions

Calculations with Gaussian orbitals and periodic boundary conditions using several density functionals are carried out to study the proton-doping of polyaniline. We explore previously proposed mechanisms to explain the spectacular increase of the electrical conductivity of polyaniline upon protonation. The structural and spectroscopic theoretical predictions for both polaron and bipolaron lattices agree quite well with the experimental data, and we find that the bipolaron structure is lower in energy.     

Periodic trends in indirect nuclear spin-spin coupling tensors: relativistic density functional calculations for interhalogen diatomics

There have been significant advances in the calculation and interpretation of indirect nuclear spin-spin coupling (J) tensors during the past few years; however, much work remains to be done, especially for molecules containing heavy atoms where relativistic effects may play an important role. Many J tensors cannot be explained based solely on a nonrelativistic Fermi-contact mechanism. In the present work, the relativistic zeroth-order regular approximation density-functional (ZORA-DFT) implementation for the calculation of J has been applied to the complete series of homonuclear and heteronuclear diatomic halogen molecules: F(2), Cl(2), Br(2), I(2), At(2), ClF, BrF, IF, ClBr, ClI, and BrI. For all of these compounds, the reduced isotropic coupling constant (K(iso)) is positive and the reduced anisotropic coupling constant (DeltaK) is negative. With the exception of molecular fluorine, the magnitudes of K(iso) and DeltaK are shown to increase linearly with the product of the atomic numbers of the coupled nuclei. ZORA-DFT calculations of J for F(2) and ClF are in excellent agreement with the results obtained from multiconfigurational self-consistent-field calculations. The relative importance of the various coupling mechanisms is approximately constant for all of the compounds, with the paramagnetic spin-orbit term being the dominant contributor to K(iso), at approximately 70-80%. Available experimental stimulated resonant Raman spectroscopy data are exploited to extract the complete J((127)I,(127)I) tensor for iodine in two rotational states. The dependence of K(iso) and DeltaK on bond length and rovibrational state is investigated by using calculated results in combination with available experimental data. In addition to providing new insights into periodic trends for J coupling tensors, this work further demonstrates the utility of the ZORA-DFT method and emphasizes the necessity of spin-orbit relativistic corrections for J calculations involving heavy nuclei.     

Relativistic all-electron configuration interaction calculations on the gold atom

We present the results of relativistic and non-relativistic self-consistent field and configuration interaction calculations for the gold atom, using the spin-free no-pair Hamiltonian in a basis set expansion. A new basis set for the gold atom is discussed and its results in relativistic and non-relativistic self-consistent field calculations are compared to those of numerical Dirac-Hartree-Focic and Hartree-Fock calculations, respectively. Excitation energies, electron affinities and ionization potentials were calculated using a multi-reference configuration interaction technique and are in reasonable agreement with experiment in the relativistic case.     

Implementation of empirical dispersion corrections to density functional theory for periodic systems

A recently developed empirical dispersion correction (Grimme et al., J. Chem. Phys. 2010, 132, 154104) to standard density functional theory (DFT‐D3) is implemented in the plane‐wave program package VASP. The DFT‐D3 implementation is compared with an implementation of the earlier DFT‐D2 version (Grimme, J. Comput. Chem. 2004, 25, 1463; Grimme, J. Comput. Chem. 2006, 27, 1787). Summation of empirical pair potential terms is performed over all atom pairs in the reference cell and over atoms in shells of neighboring cells until convergence of the dispersion energy is obtained. For DFT‐D3, the definition of coordination numbers has to be modified with respect to the molecular version to ensure convergence. The effect of three‐center terms as implemented in the original molecular DFT‐D3 version is investigated. The empirical parameters are taken from the original DFT‐D3 version where they had been optimized for a reference set of small molecules. As the coordination numbers of atoms in bulk and surfaces are much larger than in the reference compounds, this effect has to be discussed. The results of test calculations for bulk properties of metals, metal oxides, benzene, and graphite indicate that the original parameters are also suitable for solid‐state systems. In particular, the interlayer distance in bulk graphite and lattice constants of molecular crystals is considerably improved over standard functionals. With the molecular standard parameters (Grimme et al., J. Chem. Phys. 2010, 132, 154104; Grimme, J. Comput. Chem. 2006, 27, 1787) a slight overbinding is observed for ionic oxides where dispersion should not contribute to the bond. For simple adsorbate systems, such as Xe atoms and benzene on Ag(111), the DFT‐D implementations reproduce experimental results with a similar accuracy as more sophisticated approaches based on perturbation theory (Rohlfing and Bredow, Phys. Rev. Lett. 2008, 101, 266106). © 2012 Wiley Periodicals, Inc.     

The heat of formation of gaseous PuO(2)2+ from relativistic density functional calculations

Using a set of model reactions, we estimated the heat of formation of gaseous PuO2(2+) from quantum-chemical reaction enthalpies and experimental heats of formation of reference species. To this end, we carried out relativistic density functional calculations on the molecules PuO(2)2+, PuO2, PuF6, and PuF4. We used a revised variant (PBEN) of the Perdew-Burke-Ernzerhof gradient-corrected exchange-correlation functional, and we accounted for spin-orbit interaction in a self-consistent fashion. As open-shell Pu species with two or more unpaired 5f electrons are involved, spin-orbit interaction significantly affects the energies of the model reactions. Our theoretical estimate for the heat of formation DeltafH degree 0(PuO2(2+),g), 418+/-15 kcal mol-1, evaluated using plutonium fluorides as references, is in good agreement with a recent experimental result, 413+/-16 kcal mol-1. The theoretical value connected to the experimental heat of formation of PuO2(g) has a notably higher uncertainty and therefore was not included in the final result.     

An all-electron calculation of an antihypertensive protein with the Gaussian-based density functional method

An all-electron calculation of the antihypertensive and antiviral protein BDS-I from the sea anemone Anemonia sulcata was carried out using the density functional method based on Gaussian-type orbitals with workstations. Our computer program ProteinDF coded by the object-oriented language C++, was originally designed for all-electron calculations of proteins and controlled in parallel with high efficiency. The numbers of residues, atoms, orbitals and auxiliary functions of the calculated protein were 43, 629, 3615 and 6545, respectively. This Letter reports timing data of the computation and some calculated electronic properties. Our program will open the door for advanced studies of larger proteins.     

Application of semiempirical long-range dispersion corrections to periodic systems in density functional theory

Ewald summation is used to apply semiempirical long-range dispersion corrections (Grimme, J Comput Chem 2006, 27, 1787; 2004, 25, 1463) to periodic systems in density functional theory. Using the parameters determined before for molecules and the Perdew-Burke-Ernzerhof functional, structure parameters and binding energies for solid methane, graphite, and vanadium pentoxide are determined in close agreement with observed values. For methane, a lattice constant a of 580 pm and a sublimation energy of 11 kJ mol(-1) are calculated. For the layered solids graphite and vanadia, the interlayer distances are 320 pm and 450 pm, respectively, whereas the graphite interlayer energy is -5.5 kJ mol(-1) per carbon atom and layer. Only when adding the semiempirical dispersion corrections, realistic values are obtained for the energies of adsorption of C(4) alkenes in microporous silica (-66 to -73 kJ mol(-1)) and the adsorption and chemisorption (alkoxide formation) of isobutene on acidic sites in the micropores of zeolite ferrierite (-78 to -94 kJ mol(-1)). As expected, errors due to missing self-interaction correction as in the energy for the proton transfer from the acidic site to the alkene forming a carbenium ion are not affected by the dispersion term. The adsorption and reaction energies are compared with the results from M?ller-Plesset second-order perturbation theory with basis set extrapolation.     

Grid‐based density functional calculations of many‐electron systems

Exploratory variational pseudopotential density functional calculations are performed for the electronic properties of many‐electron systems in the 3D cartesian coordinate grid (CCG). The atom‐centered localized gaussian basis set, electronic density, and the two‐body potentials are set up in the 3D cubic box. The classical Hartree potential is calculated accurately and efficiently through a Fourier convolution technique. As a first step, simple local density functionals of homogeneous electron gas are used for the exchange‐correlation potential, while Hay‐Wadt‐type effective core potentials are employed to eliminate the core electrons. No auxiliary basis set is invoked. Preliminary illustrative calculations on total energies, individual energy components, eigenvalues, potential energy curves, ionization energies, and atomization energies of a set of 12 molecules show excellent agreement with the corresponding reference values of atom‐centered grid as well as the grid‐free calculation. Results for three atoms are also given. Combination of CCG and the convolution procedure used for classical Coulomb potential can provide reasonably accurate and reliable results for many‐electron systems. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Transition metal atoms on oxide supports density functional calculations

The adsorption properties of Cu, Ag, Ni, and Pd atoms on O^2?, F, and F⁺ sites of MgO, CaO, SrO, and BaO (001) surfaces have been studied by means of density functional calculations. The examined clusters were embedded in the simulated Coulomb fields that closely approximate the Madelung fields of the host surfaces. The adsorption properties have been analyzed with reference to the basicity and energy gap of the oxide support in addition to orbital interactions. While the free Ni d⁹s¹ triplet ground state is preserved on adsorption on the O^2? sites of MgO, CaO, and SrO surfaces, it is no longer preserved on the O^2? site of BaO. For all adsorbates considered, adsorption is found to be stronger on F⁺ sites compared with regular O^2? sites. While on the O^2? site, Pd and Ni form the most stable complexes, on the F site, Pd and Cu form the most stable counterparts. On the F⁺ site, the single valence electron of Cu and Ag atoms couples with the unpaired electron of the vacancy forming a covalent bond. As a result, the adsorption energies of these atoms on the F⁺ site are stronger than those on the F and O^2? sites. The adsorption properties correlate linearly with the basicity and energy gap of the oxide support in addition to orbital interactions. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Atomic multiplet energies from density functional calculations

Atomic multiplet term energies for dⁿ configurations have been estimated within density functional theory (DFT) exploiting symmetry to the largest possible extent. The electrostatic two-electron integrals, as well as term energies, are expressed in function of only three non-redundant single determinants (NRSDs), each of them being obtained from density functional calculations. The influence of correlation effects described with a gradient-corrected functional (GGA) is examined and discussed. Comparison with experimental data shows the reliability of this symmetry-based density functional approach.