Density functional theory for molecular and periodic systems using density fitting and continuous fast multipole method: Analytical gradients

A full implementation of analytical energy gradients for molecular and periodic systems is reported in the TURBOMOLE program package within the framework of Kohn–Sham density functional theory using Gaussian‐type orbitals as basis functions. Its key component is a combination of density fitting (DF) approximation and continuous fast multipole method (CFMM) that allows for an efficient calculation of the Coulomb energy gradient. For exchange‐correlation part the hierarchical numerical integration scheme (Burow and Sierka, Journal of Chemical Theory and Computation 2011, 7, 3097) is extended to energy gradients. Computational efficiency and asymptotic O(N) scaling behavior of the implementation is demonstrated for various molecular and periodic model systems, with the largest unit cell of hematite containing 640 atoms and 19,072 basis functions. The overall computational effort of energy gradient is comparable to that of the Kohn–Sham matrix formation. © 2016 Wiley Periodicals, Inc.     

Analytical nuclear excited‐state gradients for the Tamm–Dancoff approximation using uncoupled frozen‐density embedding

We report the derivation and implementation of analytical nuclear gradients for excited states using time‐dependent density functional theory using the Tamm–Dancoff approximation combined with uncoupled frozen‐density embedding using density fitting. Explicit equations are presented and discussed. The implementation is able to treat singlet as well as triplet states and functionals using the local density approximation, the generalized gradient approximation, combinations with Hartree–Fock exchange (hybrids), and range‐separated functionals such as CAM‐B3LYP. The new method is benchmarked against supermolecule calculations in two case studies: The solvatochromic shift of the (vertical) fluorescence energy of 4‐aminophthalimide on solvation, and the first local excitation of the benzonitrile dimer. Whereas for the 4‐aminophthalimide–water complex deviations of about 0.2 eV are obtained to supermolecular calculations, for the benzonitrile dimer the maximum error for adiabatic excitation energies is below 0.01 eV due to a weak coupling of the subsystems. © 2017 Wiley Periodicals, Inc.     

对含重元素体系的接合二分量-标量相对论密度泛函计算方法

在相对论密度泛函ZORA方法的基础上,提出一种用于含重元素体系的接合二分量-标量相对论密度泛函计算方法.对于只含少数几个重元素的较大体系,仅对其中旋轨耦合作用强的重元素作二分量相对论计算,而对体系的其余部分则作标量相对论计算,通过对动能矩阵元的近似处理实现两种计算的接合.对一系列含6p区重元素分子进行计算的结果表明,当非重元素是第三周期以前的元素时,此方法与二分量ZO-RA方法的计算结果吻合得很好.当非重元素为第四周期元素时,计算结果有一定偏差,表明在后一种情况下旋轨耦合作用已比较显著,但误差仍在目前近似密度泛函计算的精度范围内.此方法可以有效地节省计算量,而且避免了Dyall方法的缺点.     

Atomistic Model of Uranium

The electronic state and potential data of U₂ molecules are performed by first principle calculations with B3LYP hybrid exchange-correlation functional, the valence electrons of U atom are treated with the (5s4p3d4f)/[3s3p2d2f] contraction basis sets, and the cores are approximated with the relativistic effective core potential. The results show that the ground electronic state is X⁹∑⁺_g . The pair potential data are fitted with a Murrell-Sorbie analytical potential function. The U-U embedded atom method (EAM) interatomic potential is deter-mined based on the generalized gradient approximation calculation within the framework of the density functional theory using Perdew-Burke-Ernzerhof exchange-correlation functional at the spin-polarized level. The physical properties, such as the cohesive energy, the lattice constant, the bulk modulus, the shear modulus, the sc/fcc relative energy, the hcp/fcc rela-tive energy, the shear modulus and the monovacancy formation energy are used to evaluate the EAM potential parameters. The U-U pair potential determined by the first principle calculations is in agreement with that defined by the EAM potential parameters. The EAM calculated formation energy of the monovacancy in the fcc structure is also found to be in close agreement with DFT calculation.     

Density functional theory for molecular and periodic systems using density fitting and continuous fast multipole method: Stress tensor

A full implementation of the analytical stress tensor for periodic systems is reported in the TURBOMOLE program package within the framework of Kohn–Sham density functional theory using Gaussian-type orbitals as basis functions. It is the extension of the implementation of analytical energy gradients (Lazarski et al., Journal of Computational Chemistry 2016, 37, 2518–2526) to the stress tensor for the purpose of optimization of lattice vectors. Its key component is the efficient calculation of the Coulomb contribution by combining density fitting approximation and continuous fast multipole method. For the exchange-correlation (XC) part the hierarchical numerical integration scheme (Burow and Sierka, Journal of Chemical Theory and Computation 2011, 7, 3097–3104) is extended to XC weight derivatives and stress tensor. The computational efficiency and favorable scaling behavior of the stress tensor implementation are demonstrated for various model systems. The overall computational effort for energy gradient and stress tensor for the largest systems investigated is shown to be at most two and a half times the computational effort for the Kohn–Sham matrix formation. © 2019 Wiley Periodicals, Inc.     

Scalar relativistic all-electron density functional calculations on periodic systems

Scalar relativistic effects are included in periodic boundary conditions calculations with Gaussian orbitals. This approach is based on the third-order Douglas-Kroll-Hess approximation, allowing the treatment of all electrons on an equal footing. With this methodology, we are able to perform relativistic all-electron density functional calculations using the traditional local spin-density and generalized gradient approximations (GGA), as well as meta-GGA and hybrid density functionals. We present benchmark results for the bulk metals Pd, Ag, Pt, and Au, and the large band gap semiconductors AgF and AgCl.     

几种密度泛函理论公式用于镧系硫属化合物计算的比较

以镧系元素La、Gd、Lu的硫属化合物为对象,系统考察几种密度泛函理论公式对镧系化合物计算的适用情况,考虑了相对论效应的影响. 计算结果显示,相对论效应引起的键长变化在十2Pm到-3Pm之间,引起的键能减小为0. 4～0. 6eV,与采用的密度泛函公式关系不明显. 不同的密度泛函公式对键长的计算结果影响也不太大,但对键能有显著影响,其中LDA(VWN)＋PW86X公式给出最好结果. 交换能梯度校正明显改善键能计算结果,而相关能梯度校正反而使之变差. 简单的X_α公式给出相当好的键能计算结果. 在考虑相对论效应和梯度校正以后,密度泛函理论方法给出比较可靠的键长数值,键能则仍然偏高,但不超过20％.     

Density functional theory of field theoretical systems

The field theoretical background of relativistic density functional theory is emphasized and its consequences for relativistic Kohn-Sham equations are shown. The local density approximation for the exchange energy functional is reviewed and the importance of relativistic corrections for an accurate representation of the exchange functional is demonstrated. © 1995 John Wiley & Sons, Inc.     

Calculation of electric-field gradients based on higher-order generalized Douglas-Kroll transformations

In this paper, the calculation of electric-field-like properties based on higher-order Douglas-Kroll-Hess (DKH) transformations is discussed. The electric-field gradient calculated within the Hartree-Fock self-consistent field framework is used as a representative property. The properties are expressed as an analytic first derivative of the four-component Dirac energy and the nth-order DKH energy, respectively. The differences between a "forward" transformation of the relativistic energy or the "back transformation" of the wave function is discussed in some detail. Detailed test calculations were carried out on the electric-field gradient at the halogen nucleus in the series HX (X=F,Cl,Br,I,At) for which extensive reference data are available. The DKH method is shown to reproduce (spin-free) four-component Dirac-Fock results to an accuracy of better than 99% which is significantly closer than previous DKH studies. The calculations of both the Hamiltonian and the property operator are shown to be essentially converged after the second-order transformation, even for elements as heavy as At. In addition, we have obtained results within the density-functional framework using the DKHZ and zeroth-order regular approximation (ZORA) methods. The latter results included picture-change effects at the scalar relativistic variant of the ZORA-4 level and were shown to be in quantitative agreement with earlier results obtained by van Lenthe and Baerends. The picture-change effects are somewhat smaller for the ZORA method compared to DKH. For heavier elements significant differences in the field gradients predicted by the two methods were found. Based on comparison with four-component Dirac-Kohn-Sham calculations, the DKH results are more accurate. Compared to the spin-free Dirac-Kohn-Sham reference values, the ZORA-4 formalism did not improve the results of the ZORA calculations.     

Adsorption of superheavy elements on metal surfaces

Fully relativistic four-component density functional theory with the general gradient approximation calculations have been performed to determine the adsorption energy and position of the superheavy element 112 on a Au surface. Extended cluster as well as embedded cluster calculations were used to simulate the surface which allow for the top, bridge, and hollow adsorption sites without losing the advantage of symmetry considerations. Comparison with analogous calculations of the adsorption of the homologue element Hg allows to predict the adsorption of element 112 at a bridge site with a binding energy of 0.67 eV.     

Analytic derivatives for perturbatively corrected "double hybrid" density functionals: theory, implementation, and applications

A recently proposed new family of density functionals [S. Grimme, J. Chem. Phys. 124, 34108 (2006)] adds a fraction of nonlocal correlation as a new ingredient to density functional theory (DFT). This fractional correlation energy is calculated at the level of second-order many-body perturbation theory (PT2) and replaces some of the semilocal DFT correlation of standard hybrid DFT methods. The new "double hybrid" functionals (termed, e.g., B2-PLYP) contain only two empirical parameters that have been adjusted in thermochemical calculations on parts of the G2/3 benchmark set. The methods have provided the lowest errors ever obtained by any DFT method for the full G3 set of molecules. In this work, the applicability of the new functionals is extended to the exploration of potential energy surfaces with analytic gradients. The theory of the analytic gradient largely follows the standard theory of PT2 gradients with some additional subtleties due to the presence of the exchange-correlation terms in the self-consistent field operator. An implementation is reported for closed-shell as well as spin-unrestricted reference determinants. Furthermore, the implementation includes external point charge fields and also accommodates continuum solvation models at the level of the conductor like screening model. The density fitting resolution of the identity (RI) approximation can be applied to the evaluation of the PT2 part with large gains in computational efficiency. For systems with approximately 500-600 basis functions the evaluation of the double hybrid gradient is approximately four times more expensive than the calculation of the standard hybrid DFT gradient. Extensive test calculations are provided for main group elements and transition metal containing species. The results reveal that the B2-PLYP functional provides excellent molecular geometries that are superior compared to those from standard DFT and MP2.     

First-order relativistic corrections to MP2 energy from standard gradient codes: Comparison with results from density functional theory

The evaluation of the first-order scalar relativistic corrections to MP2 energy based on either direct perturbation theory or the mass–velocity and Darwin terms is discussed. In a basis set of Lévy-Leblond spinors the one- and two-electron matrix elements of the relativistic Hamiltonian can be decomposed into a nonrelativistic part and a relativistic perturbation. Thus, a program capable of calculating nonrelativistic energy gradients can be used to calculate the cross-term between relativity and correlation. The method has been applied to selected closed-shell atoms (He, Be, Ne, and Ar) and molecules (CuH, AgH, and AuH). The calculated equilibrium distances and harmonic frequencies were compared with results from first-order relativistic density functional calculations. It was found that the cross-term is not the origin of the nonadditivity of relativistic and correlation effects. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 1596–1603, 1998     

Density functional investigations of the properties and thermochemistry of UF6 and UF5 using valence-electron and all-electron approaches

The structural properties and thermochemistry of UF6 and UF5 have been investigated using both Hartree-Fock and density functional theory (DFT) approximations. Within the latter approach, the local spin-density approximation, the generalized gradient approximation, and hybrid density functionals were considered. To describe the uranium atom we employed small-core (60 electrons) and large-core (78 electrons) relativistic effective core potentials (RECPs), as well as the all-electron approximation based on the two-component third-order Douglas-Kroll-Hess Hamiltonian. For structural properties, we obtained very good agreement with experiment with DFT and both large and small-core RECPs. The best match with experiment is given by the hybrid functionals with the small-core RECP. The bond dissociation energy (BDE) was obtained from the relative energies of the fragments [UF6 --> UF5 + F], corrected for zero-point energy and spin-orbit interaction. Very good agreement was found between the BDE obtained from all-electron calculations and those calculated with the small-core RECP, while those from the large-core RECP are off by more than 50%. In order to obtain good agreement with experiment in the BDE it is imperative to work with hybrid density functionals and a small-core RECP.     

Analytical evaluation of Fukui functions and real-space linear response function

Many useful concepts developed within density functional theory provide much insight for the understanding and prediction of chemical reactivity, one of the main aims in the field of conceptual density functional theory. While approximate evaluations of such concepts exist, the analytical and efficient evaluation is, however, challenging, because such concepts are usually expressed in terms of functional derivatives with respect to the electron density, or partial derivatives with respect to the number of electrons, complicating the connection to the computational variables of the Kohn-Sham one-electron orbitals. Only recently, the analytical expressions for the chemical potential, one of the key concepts, have been derived by Cohen, Mori-Sánchez, and Yang, based on the potential functional theory formalism. In the present work, we obtain the analytical expressions for the real-space linear response function using the coupled perturbed Kohn-Sham and generalized Kohn-Sham equations, and the Fukui functions using the previous analytical expressions for chemical potentials of Cohen, Mori-Sánchez, and Yang. The analytical expressions are exact within the given exchange-correlation functional. They are applicable to all commonly used approximate functionals, such as local density approximation (LDA), generalized gradient approximation (GGA), and hybrid functionals. The analytical expressions obtained here for Fukui function and linear response functions, along with that for the chemical potential by Cohen, Mori-Sánchez, and Yang, provide the rigorous and efficient evaluation of the key quantities in conceptual density functional theory within the computational framework of the Kohn-Sham and generalized Kohn-Sham approaches. Furthermore, the obtained analytical expressions for Fukui functions, in conjunction with the linearity condition of the ground state energy as a function of the fractional charges, also lead to new local conditions on the exact functionals, expressed in terms of the second-order functional derivatives. We implemented the expressions and demonstrate the efficacy with some atomic and molecular calculations, highlighting the importance of relaxation effects.     

Density functional calculations of large systems containing heavy elements by means of the regionalization algorithm

The compounds containing transition metals orother heavy elements play important roles in chemistry.They possess diverse properties which have been util-ized in numerous functional materials. Some of theseproperties are also involved in biological processes.However, it is very difficult to study these compoundsdeeply by means of quantum chemical methods, be-cause in these systems usually there are a large amountof electrons, the distribution of the frontier orbitals isdense, the correlation be…     

闪锌矿MTe (M=Zn/Mg)的几何结构、弹性性质、电子结构和光学性质

采用基于密度泛函理论(DFT)框架下广义梯度近似(GGA)平面波超软赝势(PP-PW)方法, 计算了闪锌矿型MTe (M=Zn/Mg)的几何结构、弹性性质、电子结构和光学性质. 同时采用杂化密度泛函调准了带隙. 结果表明, 立方相ZnTe和MgTe均为直接带隙半导体材料. 所得晶格参数、弹性常数及体模量与实验数据基本吻合. 由弹性常数推导出ZnTe、MgTe的德拜温度分别为758、585 K. 研究了MTe的复介电函数、折射率、反射率和能量损失系数等光学性质, 并基于电子能带结构和态密度对光学性质进行了解释.     

Computational strategies for a four-component Dirac-Kohn-Sham program: implementation and first applications

An implementation of the generalized gradient approximation within the four-component formulation of relativistic density-functional theory using G-spinor basis sets is presented. This approach is based on the direct evaluation of the relativistic density and its gradient from the G-spinor amplitudes and gradients without explicit reference to the total density matrix. This proves to be a particularly efficient scheme, with an intrinsic computational cost that scales linearly with the number of G-spinor basis functions. In order to validate this new implementation, incorporated in the parallel version of the program BERTHA, a detailed study of the diatomic system CsAu is also reported. The spectroscopic constants D(e),r(e),omega(e), and x(e)omega(e) and the dipole moment mu have been calculated and compared with the best available theoretical and experimental data. The sensitivity of our results to the details of the numerical schemes used to evaluate the matrix elements is analyzed in detail. Also presented is a comparative study of molecular properties in the alkali auride series which have been obtained using several standard non-relativistic density functionals.     

Homogeneous nucleation in vapor-liquid phase transition of Lennard-Jones fluids: a density functional theory approach

Density functional theory (DFT) with square gradient approximation for the free energy functional and a model density profile are used to obtain an analytical expression for the size-dependent free energy of formation of a liquid drop from the vapor through the process of homogeneous nucleation, without invoking the approximations used in classical nucleation theory (CNT). The density of the liquid drop in this work is not the same as the bulk liquid density but it corresponds to minimum free energy of formation of the liquid drop. The theory is applied to study the nucleation phenomena from supersaturated vapor of Lennard-Jones fluid. The barrier height predicted by this theory is significantly lower than the same in CNT which is rather high. The density at the center of the small liquid drop as obtained through optimization is less than the bulk density which is in agreement with other earlier works. Also proposed is a sharp interface limit of the proposed DFT of nucleation, which is as simple as CNT but with a modified barrier height and this modified classical nucleation theory, as we call it, is shown to lead to improved results.     

Analytical time-dependent density functional derivative methods within the RI-J approximation, an approach to excited states of large molecules

Time-dependent density functional theory (TDDFT) is now well established as an efficient method for molecular excited state treatments. In this work, we introduce the resolution of the identity approximation for the Coulomb energy (RI-J) to excited state gradient calculations. In combination with nonhybrid functionals, the RI-J approximation leads to speed ups in total timings of an order of magnitude compared to the conventional method; this is demonstrated for oligothiophenes with up to 40 monomeric units and adamantane clusters. We assess the accuracy of the computed adiabatic excitation energies, excited state structures, and vibrational frequencies on a set of 36 excited states. The error introduced by the RI-J approximation is found to be negligible compared to deficiencies of standard basis sets and functionals. Auxiliary basis sets optimized for ground states are suitable for excited state calculations with small modifications. In conclusion, the RI-J approximation significantly extends the scope of applications of analytical TDDFT derivative methods in photophysics and photochemistry.     

Calculation of indirect nuclear spin-spin coupling constants within the regular approximation for relativistic effects

A new method for calculating the indirect nuclear spin-spin coupling constant within the regular approximation to the exact relativistic Hamiltonian is presented. The method is completely analytic in the sense that it does not employ numeric integration for the evaluation of relativistic corrections to the molecular Hamiltonian. It can be applied at the level of conventional wave function theory or density functional theory. In the latter case, both pure and hybrid density functionals can be used for the calculation of the quasirelativistic spin-spin coupling constants. The new method is used in connection with the infinite-order regular approximation with modified metric (IORAmm) to calculate the spin-spin coupling constants for molecules containing heavy elements. The importance of including exact exchange into the density functional calculations is demonstrated.