As the foundation for a new computational implementation, we survey the calculation of the complex electrical conductivity tensor based on the Kubo-Greenwood (KG) formalism (Kubo, 1957; Greenwood, 1958), with emphasis on derivations and technical aspects pertinent to use of projector augmented wave datasets with plane wave basis sets (BIlichl, 1994). New analytical results and a full implementation of the KG approach in an open-source Fortran 90 post-processing code for use with Quantum Espresso (Giannozzi et al., 2009) are presented. Named KGEC ([K]ubo [G]reenwood [E]lectronic [C]onductivity), the code calculates the full complex conductivity tensor (not just the average trace). It supports use of either the original KG formula or the popular one approximated in terms of a Dirac delta function. It provides both Gaussian and Lorentzian representations of the Dirac delta function (though the Lorentzian is preferable on basic grounds). KGEC provides decomposition of the conductivity into intra- and inter band contributions as well as degenerate state contributions. It calculates the dc conductivity tensor directly. It is MPI parallelized over k-points, bands, and plane waves, with an option to recover the plane wave processes for their use in band parallelization as well. It is designed to provide rapid convergence with respect to k-point density. Examples of its use are given.

This thesis concerns quadratic response properties and their application to properties in Jablonski diagrams such as resonant two-photon absorption and excited state absorption. Our main interest lies in optical power limiting applications, and in this context, molecules containing heavy metal atoms prove superior. Therefore, we are interested in how relativity affects these properties, and in order to assess this, a four-component relativistic framework is adopted. To properly address the molecular properties of interest, both relativistic effects and electron correlation need to be accounted for. These two properties are not additive, and, therefore, correlation needs to be incorporated into the four-component framework. We present the implementation of quadratic response properties at the four-component density functional level of theory. For second-harmonic generation, we have, with numerical examples, demonstrated that correlation and relativity are indeed not additive and that the inclusion of noncollinear magnetization is of little importance. We report that both electron correlation as well as relativity strongly affect results for second-harmonic generation. For example, relativity alone reduces the µβ-response signal by 62% and 75% for meta- and ortho-bromobenzene, respectively, and enhances the same response by 17% and 21% for meta- and ortho-iodobenzene, respectively. In the four-component framework, we present the implementations of single and double residues of the quadratic response function, which allows for the evaluation of resonant two-photon absorption cross sections and excited state properties. Using these tools, we discuss different levels of approximation to the relativistic Hamiltonian and we demonstrate that for two-photon absorption, a proper treatment of relativistic effects qualitatively alters the spectrum. For example, already for an element as light as neon, significant differences are seen between the relativistic and nonrelativistic spectra as triplet transitions acquire substantial absorption cross sections in the former case. Finally, quantum mechanics in conjunction with electrodynamics is applied to determine clamping levels in macroscopic samples. The microscopic properties of the optically active chromophores are determined by response theory, and then, electrodynamics is used to describe the interactions between the chromophores and incident laser pulses. Using this approach a series of molecules have been investigated and their performances have been compared and ranked in order to find novel materials for optical power limiting applications.