Exchange energy gradients with respect to atomic positions and cell parameters within the Hartree-Fock Gamma-point approximation

Recently, linear scaling construction of the periodic exact Hartree-Fock exchange matrix within the Gamma-point approximation has been introduced [J. Chem. Phys. 122, 124105 (2005)]. In this article, a formalism for evaluation of analytical Hartree-Fock exchange energy gradients with respect to atomic positions and cell parameters at the Gamma-point approximation is presented. While the evaluation of exchange gradients with respect to atomic positions is similar to those in the gas phase limit, the gradients with respect to cell parameters involve the accumulation of atomic gradients multiplied by appropriate factors and a modified electron repulsion integral (ERI). This latter integral arises from use of the minimum image convention in the definition of the Gamma-point Hartree-Fock approximation. We demonstrate how this new ERI can be computed with the help of a modified vertical recurrence relation in the frame of the Obara-Saika and Head-Gordon-Pople algorithm. As an illustration, the analytical gradients have been used in conjunction with the QUICCA algorithm [K. Nemeth and M. Challacombe, J. Chem. Phys. 121, 2877 (2004)] to optimize periodic systems at the Hartree-Fock level of theory.     

Linear scaling computation of the Fock matrix. VII. Periodic density functional theory at the Gamma point

Linear scaling quantum chemical methods for density functional theory are extended to the condensed phase at the Gamma point. For the two-electron Coulomb matrix, this is achieved with a tree-code algorithm for fast Coulomb summation [M. Challacombe and E. Schwegler, J. Chem. Phys. 106, 5526 (1997)], together with multipole representation of the crystal field [M. Challacombe, C. White, and M. Head-Gordon, J. Chem. Phys. 107, 10131 (1997)]. A periodic version of the hierarchical cubature algorithm [M. Challacombe, J. Chem. Phys. 113, 10037 (2000)], which builds a telescoping adaptive grid for numerical integration of the exchange-correlation matrix, is shown to be efficient when the problem is posed as integration over the unit cell. Commonalities between the Coulomb and exchange-correlation algorithms are discussed, with an emphasis on achieving linear scaling through the use of modern data structures. With these developments, convergence of the Gamma-point supercell approximation to the k-space integration limit is demonstrated for MgO and NaCl. Linear scaling construction of the Fockian and control of error is demonstrated for RBLYP6-21G* diamond up to 512 atoms.     

Electrical effects on the vibrational transitions of hydrogen fluoride due to hydrogen bonding and applied fields

Derivative Hartree-Fock (DHF) theory has been used to study the transition dipole moment of hydrogen fluoride in several hydrogen-bonded complexes and in the presence of applied fields. DHF is an open-ended, analytical means for finding energy derivatives with respect to any number of parameters, and large basis set calculations going through the seventh dipole hyperpolarizability (ninth derivative) are reported. Using multipole moments, multipole polarizabilities and hyperpolarizabilities, the intermolecular electrical influence on vibrational transitions is analyzed.     

Gradients for two-component quasirelativistic methods. Application to dihalogenides of element 116

The authors report the implementation of geometry gradients for quasirelativistic two-component Hartree-Fock and density functional methods using either the zero-order regular approximation Hamiltonian or spin-dependent effective core potentials. The computational effort of the resulting program is comparable to that of corresponding nonrelativistic calculations, as it is dominated by the evaluation of derivative two-electron integrals, which is the same for both types of calculations. Besides the implementation of derivatives of matrix elements of the one-particle Hamiltonian with respect to nuclear displacements, the calculation of the derivative exchange-correlation energy for the open shell case involves complicated expressions because of the noncollinear approach chosen to define the spin density. A pilot application to dihalogenides of element 116 shows how spin-orbit coupling strongly affects the chemistry of the superheavy p-block elements. While these molecules are bent at a scalar-relativistic level, spin-orbit coupling is so strong that only the 7p3/2 atomic orbitals of element 116 are involved in bonding, which favors linear molecular geometries for dihalogenides with heavy terminal halogen atoms.     

Efficient evaluation of analytic vibrational frequencies in Hartree-Fock and density functional theory for periodic nonconducting systems

We report a method for the efficient evaluation of analytic energy second derivatives with respect to in-phase nuclear coordinate displacements within Hartree-Fock and Kohn-Sham density functional theories using Gaussian orbitals and periodic boundary conditions. The use of an atomic orbital formulation for all computationally challenging steps allows us to adapt the direct space fast multipole method for the Coulomb-type infinite summations. Our implementation also exploits the local character of the exact Hartree-Fock exchange in nonconducting systems. Exchange-correlation contributions are computed using extensive screening and fast numerical quadratures. We benchmark our scheme for in-phase vibrational frequencies of a trans-polyacetylene chain, a two-dimensional boron nitride sheet, and bulk diamond with the 6-31G** basis set and various density functionals. A study of computational scaling with the size of the unit cell for trans-polyacetylene reveals subquadratic scaling for our scheme.     

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.     

Linear-scaling calculation of static and dynamic polarizabilities in Hartree-Fock and density functional theory for periodic systems

We present a linear-scaling method for analytically calculating static and dynamic polarizabilities with Hartree-Fock and density functional theory, using Gaussian orbitals and periodic boundary conditions. Our approach uses the direct space fast multipole method to evaluate the long-range Coulomb contributions. For exact exchange, we use efficient screening techniques developed for energy calculations. We then demonstrate the capabilities of our approach with benchmark calculations on one-, two-, and three-dimensional systems.     

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.     

The calculation of the vibrational frequencies of crystalline compounds and its implementation in the CRYSTAL code

The problem of numerical accuracy in the calculation of vibrational frequencies of crystalline compounds from the hessian matrix is discussed with reference to alpha-quartz (SiO(2)) as a case study and to the specific implementation in the CRYSTAL code. The Hessian matrix is obtained by numerical differentiation of the analytical gradient of the energy with respect to the atomic positions. The process of calculating vibrational frequencies involves two steps: the determination of the equilibrium geometry, and the calculation of the frequencies themselves. The parameters controlling the truncation of the Coulomb and exchange series in Hartree-Fock, the quality of the grid used for the numerical integration of the Exchange-correlation potential in Density Functional Theory, the SCF convergence criteria, the parameters controlling the convergence of the optimization process as well as those controlling the accuracy of the numerical calculation of the Hessian matrix can influence the obtained vibrational frequencies to some extent. The effect of all these parameters is discussed and documented. It is concluded that with relatively economical computational conditions the uncertainty related to these parameters is smaller than 2-4 cm(-1). In the case of the Local Density Approximation scheme, comparison is possible with recent calculations performed with a Density Functional Perturbation Theory method and a plane-wave basis set.     

Linear-scaling formation of Kohn-Sham Hamiltonian: application to the calculation of excitation energies and polarizabilities of large molecular systems

We present calculations of excitation energies and polarizabilities in large molecular systems at the local-density and generalized-gradient approximation levels of density-functional theory (DFT). Our results are obtained using a linear-scaling DFT implementation in the program system DALTON for the formation of the Kohn-Sham Hamiltonian. For the Coulomb contribution, we introduce a modification of the fast multipole method to calculations over Gaussian charge distributions. It affords a simpler implementation than the original continuous fast multipole method by partitioning the electrostatic Coulomb interactions into "classical" and "nonclassical" terms which are explicitly evaluated by linear-scaling multipole techniques and a modified two-electron integral code, respectively. As an illustration of the code, we have studied the singlet and triplet excitation energies as well as the static and dynamic polarizabilities of polyethylenes, polyenes, polyynes, and graphite sheets with an emphasis on the trends observed with system size.     

Analytical Hartree–Fock gradients with respect to the cell parameter for systems periodic in three dimensions

Analytical Hartree–Fock gradients with respect to the cell parameter have been implemented in the electronic structure code CRYSTAL, for the case of three-dimensional periodicity. The code is based on Gaussian-type orbitals, and the summation of the Coulomb energy is performed with the Ewald method. It is shown that a high accuracy of the cell gradient can be achieved.     

Efficient implementation of the fast multipole method

A number of computational techniques are described that reduce the effort related to the continuous fast multipole method, used for the evaluation of Coulomb matrix elements as needed in Hartree-Fock and density functional theories. A new extent definition for Gaussian charge distributions is proposed, as well as a new way of dividing distributions into branches. Also, a new approach for estimating the error caused by truncation of multipole expansions is presented. It is found that the use of dynamically truncated multipole expansions gives a speedup of a factor of 10 in the work required for multipole interactions, compared to the case when all interactions are computed using a fixed multipole expansion order. Results of benchmark calculations on three-dimensional systems are reported, demonstrating the usefulness of our present implementation of the fast multipole method.     

Analytical Hessian of electronic excited states in time-dependent density functional theory with Tamm-Dancoff approximation

We present the analytical expression and computer implementation for the second-order energy derivatives of the electronic excited state with respect to the nuclear coordinates in the time-dependent density functional theory (TDDFT) with Gaussian atomic orbital basis sets. Here, the Tamm-Dancoff approximation to the full TDDFT is adopted, and therefore the formulation process of TDDFT excited-state Hessian is similar to that of configuration interaction singles (CIS) Hessian. However, due to the replacement of the Hartree-Fock exchange integrals in CIS with the exchange-correlation kernels in TDDFT, many quantitative changes in the derived equations are arisen. The replacement also causes additional technical difficulties associated with the calculation of a large number of multiple-order functional derivatives with respect to the density variables and the nuclear coordinates. Numerical tests on a set of test molecules are performed. The simulated excited-state vibrational frequencies by the analytical Hessian approach are compared with those computed by CIS and the finite-difference method. It is found that the analytical Hessian method is superior to the finite-difference method in terms of the computational accuracy and efficiency. The numerical differentiation can be difficult due to root flipping for excited states that are close in energy. TDDFT yields more exact excited-state vibrational frequencies than CIS, which usually overestimates the values.     

On the calculation of dipole moment and polarizability derivatives by the analytical energy gradient method: Application to the formaldehyde molecule

Following the suggestion of Komornicki and McIver we have implemented an efficient computational scheme for the evaluation of dipole moment and polarizability derivatives at the Hartree-Fock SCF level. The derivatives are obtained by utilizing the analytical gradients of the molecular energy, calculated in the presence of an external electric field, with respect to the atomic cartesian coordinates, which are differentiated numerically with respect to the field. The implementation of the method within the framework of the MOLECULE program is discussed, concentrating on such aspects as numerical accuracy, utilization of molecular symmetry and computational efficiency. As an application, the dipole moment and polarizability derivatives of the formaldehyde molecule have been calculated, yielding infrared intensities and Raman scattering activities in the double harmonic approximation. The theoretical results are compared with the available experimental data; the agreement is satisfactory given the inherent restrictions of the SCF model.     

Step structure in the atomic Kohn-Sham potential

In this work we analyze the exchange-correlation potentialv _xc within the Kohn-Sham approach to density functional theory for the case of atomic systems. The exchange-correlation potential is written as the sum of two potentials. One of these potentialsv _xc,scr is the long-range. Coulombic potential of the coupling constant integrated exchange-correlation hole which represents the screening of the two-particle interactions due to exchange-correlation effects. The other potentialv _xc,scr ^resp contains the functional derivative with respect to the electron density of the coupling constant integrated pair-correlation function representing the sensitivity of this exchange-correlation screening to density variations. As explicit expression of the exchange-part of this functional derivative is derived using an approximation for the Greens function of the Kohn-Sham system and is shown to display a distinct atomic shell structure. The corresponding potentialv _xc,scr ^resp has a clear step structure and is constant within the atomic shells and changes rapidly at the atomic shell boundaries. Numerical examples are presented for the Be and Kr atoms using the Optimized Potential Model (OPM).     

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.     

On-the-fly localization of electronic orbitals in Car-Parrinello molecular dynamics

The ab initio molecular-dynamics formalism of Car and Parrinello is extended to preserve the locality of the orbitals. The supplementary term in the Lagrangian does not affect the nuclear dynamics, but ensures "on the fly" localization of the electronic orbitals within a periodic supercell in the Gamma-point approximation. The relationship between the resulting equations of motion and the formation of a gauge-invariant Lagrangian combined with a gauge-fixing procedure is briefly discussed. The equations of motion can be used to generate a very stable and easy to implement numerical integration algorithm. It is demonstrated that this algorithm can be used to compute the trajectory of the maximally localized orbitals, known as Wannier orbitals, in ab initio molecular dynamics with only a modest increase in the overall computer time. In the present paper, the new method is implemented within the generalized gradient approximation to Kohn-Sham density-functional theory employing plane wave basis sets and atomic pseudopotentials. In the course of the presentation, we briefly discuss how the present approach can be combined with localized basis sets to design fast linear scaling ab initio molecular-dynamics methods.     

Real-space post-Hartree-Fock correlation models

Real-space models of nondynamical correlation between electrons of opposite spin and of parallel spin in multicenter (molecular) systems are discussed. These models are designed to be partnered with Hartree-Fock or exact Kohn-Sham exchange. Thus the numerous and well-known problems of local density-functional exchange-correlation approximations, especially in stretched odd-electron systems, are circumvented.