The Discrete Reaction Field approach for calculating solvent effects. In: Canuto S, editor. Solvation Effects on Molecules and Biomolecules. Dordrecht: Springer Netherlands; 2008. pp. 39-102. [DOI: 10.1007/978-1-4020-8270-2_3]

Multistate, Polarizable QM/MM Embedding Scheme Based on the Direct Reaction Field Method: Solvatochromic Shifts, Analytical Gradients and Optimizations of Conical Intersections in Solution

We recently introduced a polarizable embedding scheme based on an integral-exact reformulation of the direct reaction field method (IEDRF) that accounts for the differential solvation of ground and excited states in QM/MM simulations. The polarization and dispersion interactions between the quantum-mechanical (QM) and molecular-mechanical (MM) regions are described by the DRF Hamiltonian, while the Pauli repulsion between explicitly treated QM electrons and the implicit electron density around MM atoms is modeled with effective core potentials. A single Hamiltonian is used for all electronic states so that Born-Oppenheimer states belonging to the same geometry are orthogonal and state crossings are well-defined. In this work, we describe the implementation of the method using graphical processing unit acceleration in TeraChem, where it is combined with multiple electronic structure methods, including Hartree-Fock, time-dependent density functional theory, and complete active space self-consistent field. In contrast with older implementations of the DRF method, integrals of the polarization operators are evaluated exactly. Expressions for ingredients needed to construct analytical gradients and nonadiabatic coupling vectors are derived and tested by optimizing a conical intersection between two excited states in the presence of a polarizable solvent shell. The method is applied to estimate the solvent shifts of absorption energies of a series of donor-acceptor dyes having low-lying charge-transfer states. Even for a nonpolar solvent such as n-hexane, the inclusion of its static polarizability leads to non-negligible shifts that improve the agreement to essentially quantitative levels (0.03 eV) with full-system calculations. Good agreement with the positions of the experimental absorption maxima measured in solution is also observed.

Conical Intersections in Solution with Polarizable Embedding: Integral-Exact Direct Reaction Field

A common strategy to exploring the properties and reactivity of complex systems is to use quantum mechanics/molecular mechanics (QM/MM) embedding, wherein a QM region is defined and treated with electronic structure theory, and the remainder of the system is treated with a force field. Important to the description of electronic excited states, especially those of charge-transfer character, is the treatment of the coupling between the QM and MM subsystems. The state of the art is to use a polarizable force field for the MM region and mutually couple the QM wavefunction and MM induced dipoles, in addition to the usual electrostatic embedding, yielding a polarizable embedding (QM/MM-Pol) approach. However, we showed previously that current popular QM/MM-Pol approaches exhibit issues of root flipping and/or incorrect descriptions of electronic crossings in multistate calculations [J. Chem. Theory Comput. 14, 2137 (2018)]. Here, we demonstrate a solution to these problems with an integral-exact reformulation of the direct reaction field approach of Thole and Van Duijnen (QM/MM-IEDRF). The resulting embedding potential includes one- and two-electron operators and many-body dipole-induced dipole interactions and thus includes a natural description of the screening of electron-electron interactions by the MM induced dipoles. Pauli repulsion from the environment is mimicked by effective core potentials on the MM atoms. Inherent to the DRF approach is the assumption that MM dipoles respond instantaneously to the positions of the QM electrons; therefore, dispersion interactions are captured approximately. All electronic states are eigenfunctions of the same Hamiltonian, while the polarization induced in the environment and the associated energetic stabilization are unique to each state. This allows for a consistent definition of transition properties and state crossings. We demonstrate QM/MM-IEDRF by exploring the influence of a (polarizable) inert xenon matrix environment on the conical intersection underlying the photoisomerization of ethylene.

Theoretical Study of the Charge Transfer Exciton Binding Energy in Semiconductor Materials for Polymer:Fullerene-Based Bulk Heterojunction Solar Cells

Recent efforts and progress in polymer solar cell research have boosted the photovoltaic efficiency of the technology. This efficiency depends not only on the device architecture but also on the material properties. Thus, insight into the design of novel semiconductor materials is vital for the advancement of the field. This paper looks from a theoretical viewpoint into two of the factors for the design of semiconductor materials with applications to bulk heterojunction solar cells: the charge transfer exciton binding energy and the nanoscale arrangement of donor and acceptor molecules in blend systems. Being aware that the exciton dissociation of local excitons in charge transfer states initiates the charge generation process, the excited state properties of four oligomers (one donor-type: PEO-PPV; and three donor-acceptor-types: PTFB, PTB7, and PTB7-Th) and two fullerene derivatives ([60]PCBM and [70]PCBM), previously reported in the literature as having high electrical conductance, are studied. With such a study, the donor molecules, either of donor-type or donor-acceptor type, are screened as candidates for [60]PCBM- and/or [70]PCBM-based bulk heterojunctions. The charge transfer energy and charge transfer exciton binding energy of suitable donor:acceptor bulk heterojunctions, some of them not yet fabricated, are studied. Further, the charge transfer exciton binding energies of [60]PCBM- and [70]PCBM-based blends are compared. A combination of molecular dynamics simulations with calculations based on Kohn-Sham density functional theory (KS-DFT) and its time-dependent extension (KS-TDDFT) is used. An important feature of this work is that it incorporates the effect of the environment of the quantum chemical system in KS-DFT or KS-TDDFT calculations through a polarizable discrete reaction field (DRF). Our predictions in terms of the influence of the nanoscale arrangement of donor and acceptor molecules on the performance of organic solar cells indicate that bulk heterojunction morphologies for donor-acceptor-type oligomers lead to their lowest excited states having charge transfer character. Further, we find that in terms of favorable charge transfer exciton binding energy, the PTB7-Th:[70]PCBM blends outperform the other blends.

Off-planar geometry and structural instability of EDO-TTF explained by using the extended debye polarizability model for bond angles

The geometry of ethylenedioxy-tetrathiafulvalene, EDO-TTF, plays an important role in the metal-insulator transition in the charge transfer salt (EDO-TTF)(2)PF(6). The planar and off-planar geometrical conformations of the EDO-TTF molecules are explained using an extended Debye polarizability model for the bond angle. The geometrical structure of EDO-TTF is dictated by its four sulfur bond angles and these are, in turn, determined by the polarizability of the sulfur atoms. With Hartree-Fock and second-order Møller-Plesset perturbation theory calculations on EDO-TTF, TTF, H(2)S, and their oxygen and selenium substituted counterparts we confirm this hypothesis. The Debye polarizability model for bond angles relates directly the optimum bond angle with the polarizability of the center atom. Considering the (EDO-TTF)(2)PF(6) material in this light proves to be very fruitful.