A necessary first step in the development of technologies such as artificial photosynthesis is understanding the photoexcitation process within the basic building blocks of naturally occurring light harvesting complexes (LHCs). The most important of these building blocks in biological LHCs such as LHC II from green plants are the chlorophyll a (Chl a) and chlorophyll b (Chl b) chromophores dispersed throughout the protein matrix. However, efforts to describe such systems are still hampered by the lack of computationally efficient and accurate methods that are able to describe optical absorption in large biomolecules. In this work, we employ a highly efficient linear combination of atomic orbitals (LCAOs) to represent the Kohn–Sham (KS) wave functions at the density functional theory (DFT) level and perform time-dependent density functional theory (TDDFT) calculations in either the reciprocal space and frequency domain (LCAO-TDDFT-k-ω) or real space and time domain (LCAO-TDDFT-r-t) of the optical absorption spectra of Chl a and b monomers and dimers. We find that our LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t calculations reproduce results obtained with a plane-wave (PW) representation of the KS wave functions (PW-TDDFT-k-ω) but with a significant reduction in computational effort. Moreover, by applying the Gritsenko, van Leeuwen, van Lenthe, and Baerends solid and correlation derivative discontinuity correction Δx to the KS eigenenergies, with both LCAO-TDDFT-k-ω and LCAO-TDDFT-r-t methods, we are able to semiquantitatively reproduce the experimentally measured photoinduced dissociation results. This work opens the path to first principles calculations of optical excitations in macromolecular systems.

No abstract

No abstract

Understanding, optimizing, and controlling the optical absorption process, exciton gemination, and electron-hole separation and conduction in low dimensional systems is a fundamental problem in materials science. However, robust and efficient methods capable of modelling the optical absorbance of low dimensional macromolecular systems and providing physical insight into the processes involved have remained elusive. We employ a highly efficient linear combination of atomic orbitals (LCAOs) representation of the Kohn-Sham (KS) orbitals within time dependent density functional theory (TDDFT) in the reciprocal space (k) and frequency (ω) domains, as implemented within our LCAO-TDDFT-k-ω code, applying either a priori or a posteriori the derivative discontinuity correction of the exchange functional Δ x to the KS eigenenergies as a scissors operator. In so doing we are able to provide a semi-quantitative description of the photoabsorption cross section, conductivity, and dielectric function for prototypical 0D, 1D, 2D, and 3D systems within the optical limit ( q → 0 + ) as compared to both available measurements and from solving the Bethe-Salpeter equation with quasiparticle G 0 W 0 eigenvalues (G 0 W 0 -BSE). Specifically, we consider 0D fullerene (C 60 ), 1D metallic (10, 0) and semiconducting (10, 10) single-walled carbon nanotubes, 2D graphene (Gr) and phosphorene (Pn), and 3D rutile (R-TiO 2 ) and anatase (A-TiO 2 ). For each system, we also employ the spatially and energetically resolved electron-hole spectral density to provide direct physical insight into the nature of their optical excitations. These results demonstrate the reliability, applicability, efficiency, and robustness of our LCAO-TDDFT-k-ω code, and open the pathway to the computational design of macromolecular systems for optoelectronic, photovoltaic, and photocatalytic applications in silico.