Molecular properties of large systems based on the damped response theory approach
Title: Molecular properties of large systems based on the damped response theory approach
DNr: SNIC 2015/1-157
Project Type: SNIC Medium Compute
Principal Investigator: Joanna Kauczor <joaka@ifm.liu.se>
Affiliation: Linköpings universitet
Duration: 2015-05-01 – 2016-05-01
Classification: 10407
Keywords:

Abstract

Molecular properties for ground and excited states, can be determined from so-called molecular response functions and their poles and residues. In standard response theory, absorption spectra are obtained by solving a generalized eigenvalue problem. This computationally expensive iterative procedure starts from the energies of the lowest excited states, and only these lowest excitations are addressed. In consequence, it may not be possible to access all excitation energies of interest using this approach. The straightforward comparison between standard theory and experiment is therefore impossible in many interesting regions of the spectrum, e.g. the X-ray absorption region. This problem can be solved by using the complex polarization propagator approach [1], also known as damped (complex) response theory. In the past years, there has been interest in the development of an efficient solver for solving damped response equations, namely a complex polarization propagator solver (also in a linearly scaled fashion) at the Hartree-Fock and Kohn-Sham density functional level of theory. New solver, namely: the algorithm with symmetrized vectors [2] was developed and is implemented in DALTON2013 [3], a molecular electronic structure program, as well as LSDALTON [4], its linearly scaled counterpart. It may be used when dispersion coefficients [5], one-photon absorption, electronic circular dichroism, magnetic circular dichroism [6] and X-ray absorption spectra [7] are calculated using the complex polarization propagator approach. It is also possible to perform calculations of absorption spectra for large organic and biological systems in the presence of a solvent. When calculations of spectrum are performed in a chosen frequency interval, all response equations are solved simultaneously, greatly reducing the calculation time, however also increasing the memory requirements. The code has been modularized and can be exploited as an equation solver in other programs. 1. P. Norman, D. M. Bishop, H. J. Aa. Jensen and J. Oddershede, J. Chem. Phys. , 115, 10323 (2001). 2. J. Kauczor, P. Jørgensen and P. Norman, J. Chem. Theory Comput., 7, 1610 (2011). 3. DALTON, a molecular electronic structure program, Release Dalton2013 (2013), see http://daltonprogram.org/ 4. LSDALTON, a linear scaling molecular electronic structure program, Release LDalton2011 (2011), see http://daltonprogram.org/ 5. J. Kauczor, P. Norman, W. A Saidi, J. Chem. Phys.}, 138, 114107 (2013). 6. T. Fahleson, J. Kauczor, P. Norman, S. Coriani, Mol. Phys., (2013). 7. M. Ahrén, L. Selegård, F. Söderlind, M. Linares, J. Kauczor, P. Norman, P.-O. Käll and K. Uvdal, J. Nanopart. Res., 14, 1006 (2012).