School of Biological and Chemical Sciences , Queen Mary University of London , Mile End Road , London E1 4NS , U.K.
J Chem Theory Comput. 2019 Apr 9;15(4):2504-2516. doi: 10.1021/acs.jctc.8b01180. Epub 2019 Mar 21.
Understanding photoinduced processes in molecular crystals is central to the design of highly emissive materials such as organic lasers and organic light-emitting diodes. The modeling of such processes is, however, hindered by the lack of excited state methodologies tailored for these systems. Embedding approaches based on the Ewald sum can be used in conjunction with excited state electronic structure methods to model the localized excitations which characterize these materials. In this article, we describe the implementation of a two-level ONIOM(QM:QM') point charge embedding approach based on the Ewald method, the ONIOM Ewald embedded cluster (OEEC) model. An alternative self-consistent method is also considered to simulate the response of the environment to the excitation. Two molecular crystals with opposing photochemical behavior were used to benchmark the results with single reference and multireference methods. We observed that the inclusion of an explicit ground state cluster surrounding the QM region was imperative for the exploration of the excited state potential energy surfaces. Using OEEC, accurate absorption and emission energies as well as S-S conical intersections were obtained for both crystals. We discuss the implications of the use of these embedding schemes considering the degree of localization of the excitation. The methods discussed herein are implemented in an open source platform (fromage, https://github.com/Crespo-Otero-group/fromage ) which acts as an interface between popular electronic structure codes (Gaussian, Turbomole, and Molcas).
了解分子晶体中的光诱导过程对于设计高辐射材料(如有机激光器和有机发光二极管)至关重要。然而,这些过程的建模受到缺乏针对这些系统定制的激发态方法的阻碍。基于 Ewald 求和的嵌入方法可以与激发态电子结构方法结合使用,以模拟这些材料的局域激发。在本文中,我们描述了一种基于 Ewald 方法的两级 ONIOM(QM:QM')点电荷嵌入方法的实现,即 ONIOM Ewald 嵌入簇(OEEC)模型。还考虑了一种替代的自洽方法来模拟环境对激发的响应。我们使用两种具有相反光化学反应性的分子晶体来用单参考和多参考方法对结果进行基准测试。我们观察到,包含围绕 QM 区域的显式基态簇对于探索激发态势能表面至关重要。使用 OEEC,我们为两种晶体都获得了准确的吸收和发射能量以及 S-S 锥形交叉。我们讨论了这些嵌入方案的使用对激发的局域化程度的影响。本文讨论的方法在一个开源平台(fromage,https://github.com/Crespo-Otero-group/fromage)中实现,该平台充当流行电子结构代码(Gaussian、Turbomole 和 Molcas)之间的接口。