We present an implementation of the damped coupled-cluster linear response function based on an asymmetric Lanczos chain algorithm for the hierarchy of coupled-cluster approximations CCS (coupled-cluster singles), CC2 (coupled. cluster singles and approximate doubles), and CCSD (coupled-cluster singles and doubles). Triple corrections to the excitation energies can be included via the CCSDR(3) (coupled-cluster singles and doubles with noniterative-triples-corrected excitation energies) approximation. The performance and some of the potentialities of the approach are investigated in calculations of the visible/ultraviolet absorption spectrum and the dispersion of the real polarizability in near-resonant regions of pyrimidine, the near-edge absorption fine structure (NEXAFS) of ammonia, and the direct determination of the C-6 dipole-dipole dispersion coefficient of the benzene dimer.
The influences of group 12 (Zn, Cd, Hg) metal-substitution on the valence spectra and phosphorescence parameters of porphyrins (P) have been investigated in a relativistic setting. In order to obtain valence spectra, this study reports the first application of the damped linear response function, or complex polarization propagator, in the four-component density functional theory framework [as formulated in J. Chem. Phys. 133, 064105 (2010)]. It is shown that the steep increase in the density of states as due to the inclusion of spin-orbit coupling yields only minor changes in overall computational costs involved with the solution of the set of linear response equations. Comparing single-frequency to multi-frequency spectral calculations, it is noted that the number of iterations in the iterative linear equation solver per frequency grid-point decreases monotonously from 30 to 0.74 as the number of frequency points goes from one to 19. The main heavy-atom effect on the UV/vis-absorption spectra is indirect and attributed to the change of point group symmetry due to metal-substitution, and it is noted that substitutions using heavier atoms yield small red-shifts of the intense Soret-band. Concerning phosphorescence parameters, the adoption of a four-component relativistic setting enables the calculation of such properties at a linear order of response theory, and any higher-order response functions does not need to be considered. For the substituted porphyrins, electronic coupling between the lowest triplet states is strong and results in theoretical estimates of lifetimes that are sensitive to the wave function and electron density parametrization. With this in mind, we report our best estimates of the phosphorescence lifetimes to be 460, 13.8, 11.2, and 0.00155 s for H_{2}P, ZnP, CdP, and HgP, respectively, with the corresponding transition energies being equal to 1.46, 1.50, 1.38, and 0.89 eV.
Predicting the polarizabilities of extended conjugated molecules with semilocal functionals has been a long-standing problem in density functional theory. These difficulties are due to the absence of a term in the typical semilocal Kohn-Sham exchange potentials that has been named "ultranonlocal". Such a term should develop in extended systems when an external electric field is applied, and it should counteract the field. We calculate the polarizabilities of polyacetylene molecules using the recently developed extended Becke-Johnson functional. Our results show that this functional predicts the polarizabilities with much better accuracy than typical semilocal functionals. Thus, the field-counteracting term in this functional, which is semilocal in the Kohn-Sham orbitals, can realistically describe real molecules. We discuss approaches of constructing an energy functional that corresponds to this potential functional, for example, via the Levy-Perdew virial relation.
The response equations as occurring in the Hartree-Fock, multiconfigurational self-consistent field, and Kohn-Sham density functional theory have identical matrix structures. The algorithms that are used for solving these equations are discussed, and new algorithms are proposed where trial vectors are split into symmetric and antisymmetric components. Numerical examples are given to compare the performance of the algorithms. The calculations show that the standard response equation for frequencies smaller than the highest occupied molecular orbital-lowest unoccupied molecular orbital gap is best solved using the preconditioned conjugate gradient or conjugate residual algorithms where trial vectors are split into symmetric and antisymmetric components. For larger frequencies in the standard response equation as well as in the damped response equation in general, the preconditioned iterative subspace approach with symmetrized trial vectors should be used. For the response eigenvalue equation, the Davidson algorithm with either paired or symmetrized trial vectors constitutes equally good options.
Molecular spectra can be determined from molecular response functions, by solving the so-called damped response equations using the complex polarization propagator approach. The overall structure of response equations is identical for variational wave functions such as the Hartree-Fock, multi-configuration self-consistent field, and Kohn-Sham density functional theory, and the key program module is the linear response equation solver. We present an implementation of the solver using the algorithm with symmetrized vectors, optimized for addressing spectral regions of a width of some 5-10 eV and a resolution below 0.1 eV. The work is illustrated by the consideration of UV-vis as well as near carbon K -edge absorption spectra of the C-60 fullerene. We demonstrate that it is possible to converge tightly response equations for hundreds of optical frequencies in resonance regions of the spectrum at a cost not much exceeding the solution of a single response equation in the nonresonant region. Our work is implemented in the molecular orbital based module of the Dalton program and serves as a documentation of the code distributed in the Dalton2013 release version.
Third-order nonlinear optical (NLO) properties of polymethine dyes have been widely studied for applications such as all-optical switching. However, the limited accuracy of the current computational methodologies has prevented a comprehensive understanding of the nature of the lowest excited states and their influence on the molecular optical and NLO properties. Here, attention is paid to the lowest excited-state energies and their energetic ratio, as these characteristics impact the figure-of-merit for all-optical switching. For a series of model polymethines, we compare several algebraic diagrammatic construction (ADC) schemes for the polarization propagator with approximate second-order coupled cluster (CC2) theory, the widely used INDO/MRDCI approach and the symmetry adapted cluster configuration interaction (SAC-CI) algorithm incorporating singles and doubles linked excitation operators (SAC-CI SD-R). We focus in particular on the ground-to-excited state transition dipole moments and the corresponding state dipole moments, since these quantities are found to be of utmost importance for an effective description of the third-order polarizability gamma and two-photon absorption spectra. A sum-overstates expression has been used, which is, found to quickly converge. While ADC(3/2) has been found to be the most appropriate method to calculate these properties, CC2 performs poorly.
We present a combination of the polarizable embedding (PE) scheme with the complex polarization propagator (CPP) method with the aim of calculating response properties including relaxation for large and complex systems. This new approach, termed PE-CPP, will benefit from the highly advanced description of the environmental electrostatic potential and polarization in the PE method as well as the treatment of near-resonant effects in the CPP approach. The PE-CPP model has been implemented in a Kohn-Sham density functional theory approach, and we present pilot calculations exemplifying the implementation for the UV/vis and carbon K-edge X-ray absorption spectra of the protein plastocyanin. Furthermore, technical details associated with a PE-CPP calculation are discussed.
We introduce a hybrid complex polarization propagator/molecular mechanics method for the calculation of near-resonant and resonant response properties of molecules in heterogeneous environments, which consist of a metallic surface, or nanoparticle, and a solvent. The applicability and performance of the method is demonstrated by computations of linear absorption spectra of p-nitroaniline physisorbed at a gold/dimethyl sulfoxide interface in the UV/vis and near carbon-K-edge regions of the spectrum. It is shown that the shift of absorption cross-section induced by the heterogeneous environment varies significantly depending on the nature,of the excited states encountered in the targeted frequency region as well as on the actual size of the resonant frequencies, and that the solvent component of the heterogeneous environment is responsible for the major part of the environmental shift, especially in the higher frequency range of the carbon K-edge region.
We present a study of mobility field and temperature dependence for C60 with Kinetic Monte-Carlo simulations. We propose a new scheme to take into account polarization effects in organic materials through atomic induced dipoles on nearby molecules. This leads to an energy correction for the single site energies and to an external reorganization happening after each hopping. The inclusion of polarization allows us to obtain a good agreement with experiments for both mobility field and temperature dependence.