Local ab initio methods for calculating optical band gaps in periodic systems. I. Periodic density fitted local configuration interaction singles method for polymers.

We present a density fitted local configuration interaction singles (CIS) method for calculating optical band gaps in 1D-periodic systems. The method is based on the Davidson diagonalization procedure, carried out in the reciprocal space. The one-electron part of the matrix-vector products is also evaluated in the reciprocal space, where the diagonality of the Fock matrix can be exploited. The contraction of the CIS vectors with the two electron integrals is performed in the direct space in the basis of localized occupied (Wannier) and virtual (projected atomic) orbitals. The direct space approach allows to utilize the sparsity of the integrals due to the local representation and locality of the exciton. The density fitting approximation employed for the two electron integrals reduces the nominal scaling with unit cell size to O(N(4)). Test calculations on a series of prototypical systems demonstrate that the method in its present stage can be used to calculate the excitonic band gaps of polymers with up to a few dozens of atoms in the cell. The computational cost depends on the locality of the exciton, but even relatively delocalized excitons occurring in the polybiphenyl in the parallel orientation, can be routinely treated with this method.

[1]  Elizabeth A. Costner,et al.  Fundamental optical properties of linear and cyclic alkanes: VUV absorbance and index of refraction. , 2009, The journal of physical chemistry. A.

[2]  Frederick R. Manby,et al.  Fast local-MP2 method with density-fitting for crystals. I. Theory and algorithms , 2007 .

[3]  Manabu Oumi,et al.  A doubles correction to electronic excited states from configuration interaction in the space of single substitutions , 1994 .

[4]  Richard A Friesner,et al.  Pseudospectral time-dependent density functional theory. , 2008, The Journal of chemical physics.

[5]  Philippe Y. Ayala,et al.  Atomic orbital Laplace-transformed second-order Møller–Plesset theory for periodic systems , 2001 .

[6]  Tatiana Korona,et al.  Local CC2 electronic excitation energies for large molecules with density fitting. , 2006, The Journal of chemical physics.

[7]  S. Suhai.  Green's-function study of optical properties of polymers: Charge-transfer exciton spectra of polydiacetylenes , 1984 .

[8]  Daniel Kats,et al.  Local CC2 response method for triplet states based on Laplace transform: excitation energies and first-order properties. , 2010, The Journal of chemical physics.

[9]  Asbjörn M. Burow,et al.  Resolution of identity approximation for the Coulomb term in molecular and periodic systems. , 2009, The Journal of chemical physics.

[10]  F. Manby.  Accurate condensed-phase quantum chemistry , 2010 .

[11]  Poul Jørgensen,et al.  The second-order approximate coupled cluster singles and doubles model CC2 , 1995 .

[12]  Cesare Pisani,et al.  Symmetry-adapted Localized Wannier Functions Suitable for Periodic Local Correlation Methods , 2006 .

[13]  S. B. Trickeya,et al.  Chapter 6 - Methods and Implementation of Robust, High-Precision Gaussian Basis DFT Calculations for Periodic Systems: the GTOFF Code , 2004 .

[14]  Study of excited states of polyethylene in the Hartree–Fock, Tamm–Dancoff, and random‐phase approximations , 1995 .

[15]  M. Head‐Gordon,et al.  Configuration interaction singles, time-dependent Hartree-Fock, and time-dependent density functional theory for the electronic excited states of extended systems , 1999 .

[16]  J. D. Cloizeaux,et al.  Energy Bands and Projection Operators in a Crystal: Analytic and Asymptotic Properties , 1964 .

[17]  Brett I. Dunlap,et al.  Robust and variational fitting , 2000 .

[18]  N. Marzari,et al.  Maximally localized generalized Wannier functions for composite energy bands , 1997, cond-mat/9707145.

[19]  H Germany,et al.  Correlation-induced corrections to the band structure of boron nitride: a wave-function-based approach. , 2009, The Journal of chemical physics.

[20]  P. Knowles,et al.  Poisson equation in the Kohn-Sham Coulomb problem. , 2001, Physical review letters.

[21]  M. Preuss,et al.  Optical absorption of water: coulomb effects versus hydrogen bonding. , 2005, Physical review letters.

[22]  M. Challacombe,et al.  Molecular-orbital-free algorithm for excited states in time-dependent perturbation theory. , 2008, The Journal of chemical physics.

[23]  Poul Jørgensen,et al.  Response functions from Fourier component variational perturbation theory applied to a time-averaged quasienergy , 1998 .

[24]  B. Paulus.  The method of increments—a wavefunction-based ab initio correlation method for solids , 2006 .

[25]  J. D. Cloizeaux.  Analytical Properties of n-Dimensional Energy Bands and Wannier Functions , 1964 .

[26]  E. G. Wilson,et al.  Intrinsic Photoconduction and Photoemission in Polyethylene , 1973 .

[27]  L. Reining,et al.  Electronic excitations: density-functional versus many-body Green's-function approaches , 2002 .

[28]  A Marini,et al.  The Bethe–Salpeter equation: a first-principles approach for calculating surface optical spectra , 2004 .

[29]  B. Roos,et al.  Theoretical Studies of the Electronic Spectra of Organic Molecules , 1995 .

[30]  Lu J. Sham,et al.  Local-field and excitonic effects in the optical spectrum of a covalent crystal , 1975 .

[31]  Frederick R. Manby,et al.  Fast linear scaling second-order Møller-Plesset perturbation theory (MP2) using local and density fitting approximations , 2003 .

[32]  Georg Kresse,et al.  Second-order Møller-Plesset perturbation theory applied to extended systems. II. Structural and energetic properties. , 2010, The Journal of chemical physics.

[33]  Christof Hättig,et al.  Structure Optimizations for Excited States with Correlated Second-Order Methods: CC2 and ADC(2) , 2005 .

[34]  Artur F Izmaylov,et al.  Resolution of the identity atomic orbital Laplace transformed second order Møller-Plesset theory for nonconducting periodic systems. , 2008, Physical chemistry chemical physics : PCCP.

[35]  G. Kresse,et al.  Accurate quasiparticle spectra from self-consistent GW calculations with vertex corrections. , 2007, Physical review letters.

[36]  R. Dovesi,et al.  A general method to obtain well localized Wannier functions for composite energy bands in linear combination of atomic orbital periodic calculations , 2001 .

[37]  E. Davidson.  The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices , 1975 .

[38]  G. Wannier.  The Structure of Electronic Excitation Levels in Insulating Crystals , 1937 .

[39]  M. Schütz,et al.  Second Order Local Møller-Plesset Perturbation Theory for Periodic Systems: the CRYSCOR Code , 2010 .

[40]  H. Katagiri.  Equation-of-motion coupled-cluster study on exciton states of polyethylene with periodic boundary condition. , 2005, The Journal of chemical physics.

[41]  R. Evarestov,et al.  The hartree-fock method and density-functional theory as applied to an infinite crystal and to a cyclic cluster , 2002 .

[42]  J. Mintmire,et al.  Fitting the Coulomb potential variationally in linear-combination-of-atomic-orbitals density-functional calculations , 1982 .

[43]  N. Doltsinis.  Computational nanoscience: do it yourself! , 2006 .

[44]  Lucia Reining,et al.  Efficient ab initio calculations of bound and continuum excitons in the absorption spectra of semico , 2007, 0705.3140.

[45]  Jozef Noga,et al.  Density fitting of two-electron integrals in extended systems with translational periodicity: the Coulomb problem. , 2006, The Journal of chemical physics.

[46]  L. Maschio,et al.  Fitting of local densities in periodic systems , 2008 .

[47]  S. Suhai.  First principles charge transfer exciton theory of the UV spectrum of DNA , 1984 .

[48]  J.-H. Chem,et al.  SoC verification: a comprehensive, flexible methodology for complex, multimillion gate desiguis , 2003, ASIC, 2003. Proceedings. 5th International Conference on.

[49]  Jochen Schirmer,et al.  An efficient polarization propagator approach to valence electron excitation spectra , 1995 .

[50]  Claudia Ambrosch-Draxl,et al.  Lowest optical excitations in molecular crystals: bound excitons versus free electron-hole pairs in anthracene. , 2004, Physical review letters.

[51]  R. Bartlett,et al.  Many-body Green's-function calculations on the electronic excited states of extended systems , 2000 .

[52]  M. Zaider,et al.  A calculation of exciton energies in periodic systems with helical symmetry: Application to a hydrogen fluoride chain , 1992 .

[53]  J. Paier,et al.  Second-order Møller-Plesset perturbation theory applied to extended systems. I. Within the projector-augmented-wave formalism using a plane wave basis set. , 2009, The Journal of chemical physics.

[54]  Stoll,et al.  Correlation energy of diamond. , 1992, Physical review. B, Condensed matter.

[55]  Peter Pulay,et al.  Local configuration interaction: An efficient approach for larger molecules , 1985 .

[56]  Steven G. Louie,et al.  Optical Excitations in Conjugated Polymers , 1999 .

[57]  Stephen R. Langhoff,et al.  Quantum mechanical electronic structure calculations with chemical accuracy , 1995 .

[58]  Cesare Pisani,et al.  Periodic local MP2 method for the study of electronic correlation in crystals: Theory and preliminary applications , 2008, J. Comput. Chem..

[59]  Walter Kohn,et al.  Analytic Properties of Bloch Waves and Wannier Functions , 1959 .

[60]  Steven G. Louie,et al.  Electron-Hole Excitations in Semiconductors and Insulators , 1998 .

[61]  Frederick R. Manby,et al.  The Poisson equation in density fitting for the Kohn-Sham Coulomb problem , 2001 .

[62]  S. Hirata,et al.  Fast second-order many-body perturbation method for extended systems , 2009 .

[63]  Georg Hetzer,et al.  Low-order scaling local electron correlation methods. I. Linear scaling local MP2 , 1999 .

[64]  M. Schütz,et al.  On the use of the Laplace transform in local correlation methods. , 2008, Physical chemistry chemical physics : PCCP.

[65]  Per E. M. Siegbahn,et al.  The Direct Configuration Interaction Method from Molecular Integrals , 1977 .

[66]  S. Suhai.  On the excitonic nature of the first UV absorption peak in polyene , 1986 .

[67]  M. Gillan,et al.  Calculation of properties of crystalline lithium hydride using correlated wave function theory , 2009 .

[68]  M. Zaider,et al.  A study of excited states in trans‐polyacetylene in the Hartree–Fock, Tamm–Dancoff, and random‐phase approximation , 1993 .

[69]  G. Scuseria,et al.  Why are time-dependent density functional theory excitations in solids equal to band structure energy gaps for semilocal functionals, and how does nonlocal Hartree-Fock-type exchange introduce excitonic effects? , 2008, The Journal of chemical physics.