The combined density functional theory and multireference configuration interaction (DFT/MRCI) method of Grimme and Waletzke [J. Chem. Phys. 111, 5645 (1999)] is a well-established semi-empirical quantum chemical method for efficiently computing excited-state properties of organic molecules. As it turns out, the method fails to treat bi-chromophores owing to the strong dependence of the parameters on the excitation class. In this work, we present an alternative form of correcting the matrix elements of a MRCI Hamiltonian which is built from a Kohn-Sham set of orbitals. It is based on the idea of constructing individual energy shifts for each of the state functions of a configuration. The new parameterization is spin-invariant and incorporates less empirism compared to the original formulation. By utilizing damping techniques together with an algorithm of selecting important configurations for treating static electron correlation, the high computational efficiency has been preserved. The robustness of the original and redesigned Hamiltonians has been tested on experimentally known vertical excitation energies of organic molecules yielding similar statistics for the two parameterizations. Besides that, our new formulation is free from artificially low-lying doubly excited states, producing qualitatively correct and consistent results for excimers. The way of modifying matrix elements of the MRCI Hamiltonian presented here shall be considered as default choice when investigating photophysical processes of bi-chromophoric systems such as singlet fission or triplet-triplet upconversion.

In the past two decades, the combined density functional theory and multireference configuration interaction (DFT/MRCI) method has developed from a powerful approach for computing spectral properties of singlet and triplet excited states of large molecules into a more general multireference method applicable to states of all spin multiplicities. In its original formulation, it shows great efficiency in the evaluation of singlet and triplet excited states which mainly originate from local one-electron transitions. Moreover, DFT/MRCI is one of the few methods applicable to large systems that yields the correct ordering of states in extended π-systems where double excitations play a significant role. A recently redesigned DFT/MRCI Hamiltonian extends the application range of the method to bichromophores such as hydrogen-bonded or π-stacked dimers and loosely coupled donor-acceptor systems. In conjunction with a restricted-open shell Kohn-Sham optimization of the molecular orbitals, even electronically excited doublet and quartet states can be addressed. After a short outline of the general ideas behind this semi-empirical method and a brief review of alternative approaches combining density functional and multireference wavefunction theory, formulae for the DFT/MRCI Hamiltonian matrix elements are presented and the adjustments of the two-electron contributions are discussed. The performance of the DFT/MRCI variants on excitation energies of organic molecules and transition metal compounds against experimental or ab initio reference data is analyzed and case studies are presented which show the strengths and limitations of the method. Finally, an overview over the properties available from DFT/MRCI wavefunctions and further developments is given.ABBREVIATIONS: ACRXTN, 3-(9,9-dimethylacridin-10[9H]-yl)-9H-xanthen-9-one; AO, atomic orbital; CASSCF, Complete active space self-consistent field; CASPT2, Complete active space second-order perturbation theory; CCSD(T), Coupled-cluster with singles and doubles and perturbative treatment of triples; CC2, Coupled-cluster with approximate treatment of doubles; CD, Circular dichroism; CI, Configuration interaction; CIPSI, configuration interaction by perturbation with multiconfigurational zeroth-order wave functions selected by iterative

During the past decade the one-center mean-field approximation has proven to be a very appropriate framework for the accurate description of spin-orbit effects at the correlated all-electron level. Here, a new efficient code, SPOCK, is introduced that calculates spin-orbit matrix elements in the one-center mean-field approximation for multireference CI wave functions. For the first time, the computation of spin-dependent interactions within a Kohn-Sham orbital based CI (DFT/MRCI) scheme1 is made possible. The latter approach is suitable for large scale systems with up to 100-200 valence electrons. Test calculations are performed on well-known diatomic molecules and the thiocarbonyl pyranthione. Spin-orbit matrix elements show good agreement with their Hartree-Fock orbital based counterparts but are obtained at considerably lower expense, thus demonstrating the power of the method. As an application singlet-triplet couplings in thiophene are investigated that are important for the photophysics and photochemistry. Spin-orbit matrix elements between all pi --> pi* excited states are found to be small. Considerably larger spin-orbit matrix elements are observed only for cases in which pi --> sigma* excited configurations are involved.

(Time-dependent) Kohn-Sham density functional theory and a combined density functional/multi-reference configuration interaction method (DFT/MRCI) were employed to explore the ground and low-lying electronically excited states of thiophene. Spin-orbit coupling was taken into account using an efficient, nonempirical mean-field Hamiltonian. Phosphorescence lifetimes were calculated by means of spock.ci, a selecting direct multi-reference spin-orbit configuration interaction program. Throughout this paper, we use the following nomenclature: S1, S2,..., T1, T2,..., denominate electronic structures in their energetic order at the ground state minimum geometry, whereas S1, S2,..., T1, T2,..., refers to the actual order of electronic states at a given nuclear geometry. Multiple minima were found on the first excited singlet (S1) potential energy hypersurface with electronic structures S1 (piHOMO-1-->pi+piHOMO-->pi), S2 (piHOMO-->pi), and S3 (piHOMO-->sigma*) corresponding to the 2 1A1 (S1), 1 1B2 (S2), and 1 1B1 (S3) states in the vertical absorption spectrum, respectively. The S1 and S2 minimum geometries show out-of-plane deformations of the ring. The S3 electronic structure yields the global minimum on the S1 surface with an adiabatic excitation energy of merely 3.81 eV. It exhibits an asymmetric planar nuclear arrangement with one significantly elongated C-S bond. A constrained minimum energy path calculation connecting the S1 and S3 minima suggests that even low-lying vibrational levels of the S1 potential well can access the global minimum of the S1 surface. Nonradiative decay of the electronically excited singlet population to the electronic ground state via a close-by conical intersection will be fast. According to our work, this ring opening mechanism is most likely responsible for the lack of fluorescence in thiophene and the ultrafast decay of the S1 vibrational levels, as observed in time-resolved pump-probe femtosecond multiphoton ionization experiments. An alternative relaxation pathway leads from the S1 minimum via vibronic coupling to the S2 potential well followed by fast inter-system crossing to the T2 state. For an estimate of individual rate constants a quantum dynamical treatment will be required. The global minimum of the T1 surface has a chair-like nuclear conformation and corresponds to the T1 (1 3B2, piHOMO-->pi) electronic structure. Phosphorescence is weak here with a calculated radiative lifetime of 0.59 s. For the second potential well on the T1 surface with T3 (1 3B1, piHOMO-->sigma*) electronic structure, nonradiative processes are predicted to dominate the triplet decay.

We present a parallelized version of a direct selecting multireference configuration interaction (MRCI) code [S. Grimme and M. Waletzke, J. Chem. Phys. 111, 5645 (1999)]. The program can be run either in ab initio mode or as semiempirical procedure combined with density functional theory (DFT/MRCI). We have investigated the efficiency of the parallelization in case studies on carotenoids and porphyrins. The performance is found to depend heavily on the cluster architecture. While the speed-up on the older Intel Netburst technology is close to linear for up to 12-16 processes, our results indicate that it is not favorable to use all cores of modern Intel Dual Core or Quad Core processors simultaneously for memory intensive tasks. Due to saturation of the memory bandwidth, we recommend to run less demanding tasks on the latter architectures in parallel to two (Dual Core) or four (Quad Core) MRCI processes per node. The DFT/MRCI branch has been employed to study the low-lying singlet and triplet states of mini-n-beta-carotenes (n=3, 5, 7, 9) and beta-carotene (n=11) at the geometries of the ground state, the first excited triplet state, and the optically bright singlet state. The order of states depends heavily on the conjugation length and the nuclear geometry. The (1)B(u) (+) state constitutes the S(1) state in the vertical absorption spectrum of mini-3-beta-carotene but switches order with the 2 (1)A(g) (-) state upon excited state relaxation. In the longer carotenes, near degeneracy or even root flipping between the (1)B(u) (+) and (1)B(u) (-) states is observed whereas the 3 (1)A(g) (-) state is found to remain energetically above the optically bright (1)B(u) (+) state at all nuclear geometries investigated here. The DFT/MRCI method is seen to underestimate the absolute excitation energies of the longer mini-beta-carotenes but the energy gaps between the excited states are reproduced well. In addition to singlet data, triplet-triplet absorption energies are presented. For beta-carotene, where these transition energies are known from experiment, excellent agreement with our calculations is observed.