We present the explicit derivation of an approach to the multiscale description of molecules in complex environments that combines frozen-density embedding (FDE) with continuum solvation models, in particular the conductor-like screening model (COSMO). FDE provides an explicit atomistic description of molecule-environment interactions at reduced computational cost, while the outer continuum layer accounts for the effect of long-range isotropic electrostatic interactions. Our treatment is based on a variational Lagrangian framework, enabling rigorous derivations of ground- and excited-state response properties. As an example of the flexibility of the theoretical framework, we derive and discuss FDE + COSMO analytical molecular gradients for excited states within the Tamm-Dancoff approximation (TDA) and for ground states within second-order Møller-Plesset perturbation theory (MP2) and a second-order approximate coupled cluster with singles and doubles (CC2). It is shown how this method can be used to describe vertical electronic excitation (VEE) energies and Stokes shifts for uracil in water and carbostyril in dimethyl sulfoxide (DMSO), respectively. In addition, VEEs for some simplified protein models are computed, illustrating the performance of this method when applied to larger systems. The interaction terms between the FDE subsystem densities and the continuum can influence excitation energies up to 0.3 eV and, thus, cannot be neglected for general applications. We find that the net influence of the continuum in presence of the first FDE shell on the excitation energy amounts to about 0.05 eV for the cases investigated. The present work is an important step toward rigorously derived ab initio multilayer and multiscale modeling approaches. © 2017 Wiley Periodicals, Inc.

Excitonic coupling plays a key role for the understanding of excitonic energy transport (EET) in, for example, organic photovoltaics. However, the calculation of realistic systems is often beyond the applicability range of accurate wavefunction methods so that lower‐scaling semi‐empirical methods are used to model EET events. In the present work, the distance and angle dependence of excitonic couplings of dimers of selected organic molecules are evaluated for the semi‐empirical long‐range corrected density functional based tight binding (LC‐DFTB) method and spin opposite scaled second order approximate coupled cluster singles and doubles (SOS‐CC2). While semi‐empirically scaled methods can lead to slightly increased deviations for excitation energies, the excitonic couplings and their dependence on the dimer geometry are reproduced. LC‐DFTB yields a similar accuracy range as density‐functional theory (DFT) employing the ωB97X functional while the computation time is reduced by several orders of magnitude. The dependence of the exchange contributions to the excitonic couplings on the dimer geometry is analyzed assessing the calculation of Coulombic excitonic couplings from monomer local excited states only, which reduces the computational effort significantly. The present work is a necessary first step toward the simulation of excitonic energy transport using semi‐empirical methods.

No abstract

We report the extended explicitly correlated approximate coupled-cluster singles and doubles CC2(F12*)-XSP method suitable for response properties. Equations are derived using an automated approach and have subsequently been hand-coded into the computer program KOALA, in which for all two-electron integrals, density fitting is employed. Numerical results are presented for the lowest two vertical singlet excitation energies of a set of selected molecules. The results show that the CC2(F12*)-XSP method provides the correct basis-set limit with no bias to the ground state, and an excellent agreement with reference CC2 values using large basis sets is found. Using Dunning’s aug-cc-pVTZ basis, the CC2(F12*)-XSP method yields excitation energies which are converged within 1 mEh to the basis-set limit for valence excitations.

We report the implementation of a Fock-operator complete-active space self-consistent field (CAS-SCF) method combined with frozen-density embedding (FDE) into the KOALA quantum-chemistry program. The implementation is based on configuration interaction from an unrestricted reference determinant and is able to treat electronic configurations such as singlet, triplet, or quintet states embedded in a molecular environment. In order to account for possible spin polarization effects, the FDE contribution is extended to the unrestricted case. We assess the convergence obtained with the implementation at the example of a stretched lithium dimer with significant multi-reference character. The efficiency of the implementation enables the orbital optimization for 25 states in a state-average SA[S0–S10,T1–T12,Q1–Q2]-CAS(10,10)-SCF calculation for the retinal molecule using a def2-TZVP basis. The FDE ansatz leads to orbitals localized by definition on the target system, thus facilitating the orbital selection required for CAS methods in complex environments.