Density fitting in periodic systems: application to TDHF in diamond and oxides
arXiv:2005.09291 · doi:10.1063/5.0014106
Abstract
A robust density fitting method for calculating Coulomb matrix elements over Bloch functions based on calculation of two- and three-center matrix elements of the Ewald potential is described and implemented in a Gaussian orbital basis in the Exciton code. The method is tested by comparing Coulomb and exchange energies from density fitting to corresponding energies from SCF HF calculations for diamond, magnesium oxide and bulk Ne. Density fitting coefficients from the robust method are compared to coefficients from a variational method applied to wave function orbital products in bulk Ne. Four center Coulomb matrix elements from density fitting are applied to time dependent Hartree-Fock (TDHF) calculations in diamond, magnesium oxide and anatase and rutile polytypes of titanium dioxide. Shifting virtual states downwards uniformly relative to occupied states and scaling the electron-hole attraction term in the TDHF Hamiltonian by 0.4 yields good agreement with either experiment and/or Bethe-Salpeter equation calculations. This approach mirrors similar 'scissors' adjustments of occupied and virtual states and introduction of a scaled electron-hole attraction term in some time dependent DFT calculations.
References in corpus (5)
- Quasiparticle and Optical Properties of Rutile and Anatase TiO
- Gaussian and plane-wave mixed density fitting for periodic systems
- Optical to UV spectra and birefringence of SiO and TiO: First-principles calculations with excitonic effects
- A simple, efficient, and general treatment of the singularities in Hartree-Fock and exact-exchange Kohn-Sham methods for solids
- Efficient Evaluation of Exact Exchange for Periodic Systems via Concentric Atomic Density Fitting