ONIOM(QM:QM′) Electrostatic Embedding Schemes for Photochemistry in Molecular Crystals

Performance of point charge embedding schemes for excited states in molecular organic crystals

Modeling excited state processes in molecular crystals is relevant for several applications. A popular approach for studying excited state molecular crystals is to use cluster models embedded in point charges. In this paper, we compare the performance of several embedding models in predicting excited states and S1-S0 optical gaps for a set of crystals from the X23 molecular crystal database. The performance of atomic charges based on ground or excited states was examined for cluster models, Ewald embedding, and self-consistent approaches. We investigated the impact of various factors, such as the level of theory, basis sets, embedding models, and the level of localization of the excitation. We consider different levels of theory, including time-dependent density functional theory and Tamm-Dancoff approximation (TDA) (DFT functionals: ωB97X-D and PBE0), CC2, complete active space self-consistent field, and CASPT2. We also explore the impact of selection of the QM region, charge leakage, and level of theory for the description of different kinds of excited states. We implemented three schemes based on distance thresholds to overcome overpolarization and charge leakage in molecular crystals. Our findings are compared against experimental data, G0W0-BSE, periodic TDA, and optimally tuned screened range-separated functionals.

Accurate and Efficient Spin-Phonon Coupling and Spin Dynamics Calculations for Molecular Solids

Molecular materials are poised to play a significant role in the development of future optoelectronic and quantum technologies. A crucial aspect of these areas is the role of spin-phonon coupling and how it facilitates energy transfer processes such as intersystem crossing, quantum decoherence, and magnetic relaxation. Thus, it is of significant interest to be able to accurately calculate the molecular spin-phonon coupling and spin dynamics in the condensed phase. Here, we demonstrate the maturity of ab initio methods for calculating spin-phonon coupling by performing a case study on a single-molecule magnet and showing quantitative agreement with the experiment, allowing us to explore the underlying origins of its spin dynamics. This feat is achieved by leveraging our recent developments in analytic spin-phonon coupling calculations in conjunction with a new method for including the infinite electrostatic potential in the calculations. Furthermore, we make the first ab initio determination of phonon lifetimes and line widths for a molecular magnet to prove that the commonplace Born-Markov assumption for the spin dynamics is valid, but such "exact" phonon line widths are not essential to obtain accurate magnetic relaxation rates. Calculations using this approach are facilitated by the open-source packages we have developed, enabling cost-effective and accurate spin-phonon coupling calculations on molecular solids.

Organic compounds for solid state luminescence enhancement/aggregation induced emission: a theoretical perspective

Organic luminophores displaying one or more forms of luminescence enhancement in solid state are extremely promising for the development and performance optimization of functional materials essential to many modern key technologies. Yet, the effort to harness their huge potential is riddled with hurdles that ultimately come down to a limited understanding of the interactions that result in the diverse molecular environments responsible for the macroscopic response. In this context, the benefits of a theoretical framework able to provide mechanistic explanations to observations, supported by quantitative predictions of the phenomenon, are rather apparent. In this perspective, we review some of the established facts and recent developments about the current theoretical understanding of solid-state luminescence enhancement (SLE) with an accent on aggregation-induced emission (AIE). A description of the macroscopic phenomenon and the questions it raises is accompanied by a discussion of the approaches and quantum chemistry methods that are more apt to model these molecular systems with the inclusion of an accurate yet efficient simulation of the local environment. A sketch of a general framework, building from the current available knowledge, is then attempted via the analysis of a few varied SLE/AIE molecular systems from literature. A number of fundamental elements are identified offering the basis for outlining design rules for molecular architectures exhibiting SLE that involve specific structural features with the double role of modulating the optical response of the luminophores and defining the environment they experience in solid state.

Unraveling the Mechanism of ACQ-to-AIE Transformation of Fluorescein Derivatives

Although fluorescein derivatives have excellent properties and strong practicability, they are typical aggregation-induced quenching (ACQ) molecules, which are not conducive to working in the solid state. Recently, the fluorescein derivative Fl-Me with aggregation-induced emission (AIE) property was synthesized, which brought a new dawn for the research and development of fluorescein-based materials. In this study, the AIE mechanism of Fl-Me was investigated based on time-dependent density functional theory and the ONIOM method. The results revealed that an effective dark-state deactivation pathway leads to the fluorescence quenching of Fl-Me in a solution environment. Accordingly, the AIE phenomenon originates from the closure of the dark-state quenching channel. It is worth emphasizing that we found that the carbonyl group of molecular Fl-Me has intermolecular hydrogen bonding interaction with the adjacent molecules, which caused the increase of the dark-state energy in the crystalline state. Moreover, the restriction of the rotational motion and the nonexistence of the π-π stacking interaction are beneficial to the enhancement of fluorescence upon aggregation. Finally, the ACQ-to-AIE transformation mechanisms of fluorescein derivatives have been discussed. This work provides deeper insight into the photophysical mechanism for the fluorescein derivatives Fl-Me with AIE feature and eventually is expected to help researchers to develop more fluorescein-based AIE materials with remarkable properties for various fields.

Competition between ultralong organic phosphorescence and thermally activated delayed fluorescence in dichloro derivatives of 9-benzoylcarbazole

Optoelectronic materials based on metal-free organic molecules represent a promising alternative to traditional inorganic devices. Significant attention has been devoted to the development of the third generation of OLEDs which are based on the temperature-activated delayed fluorescence (TADF) mechanism. In the last few years, several materials displaying ultra-long organic phosphorescence (UOP) have been designed using strategies such as crystal engineering and halogen functionalisation. Both TADF and UOP are controlled by the population of triplet states and the energy gaps between the singlet and triplet manifolds. In this paper, we explore the competition between TADF and UOP in the molecular crystals of three dichloro derivatives of 9H-carbazol-3-yl(phenyl)methanone. We investigate the excited state mechanisms in solution and the crystalline phase and address the effects of exciton transport and temperature on the rates of direct and reverse intersystem crossing under the Marcus-Levich-Jortner model. We also analyse how the presence of isomeric impurities and the stabilisation of charge transfer states affect these processes. Our simulations explain the different mechanisms observed for the three derivatives and highlight the role of intramolecular rotation and crystal packing in determining the energy gaps. This work contributes to a better understanding of the connection between chemical and crystalline structures that will enable the design of efficient materials.

Excited-State Intramolecular Proton Transfer Dyes with Dual-State Emission Properties: Concept, Examples and Applications

Dual-state emissive (DSE) fluorophores are organic dyes displaying fluorescence emission both in dilute and concentrated solution and in the solid-state, as amorphous, single crystal, polycrystalline samples or thin films. This comes in contrast to the vast majority of organic fluorescent dyes which typically show intense fluorescence in solution but are quenched in concentrated media and in the solid-state owing to π-stacking interactions; a well-known phenomenon called aggregation-caused quenching (ACQ). On the contrary, molecular rotors with a significant number of free rotations have been engineered to show quenched emission in solution but strong fluorescence in the aggregated-state thanks to restriction of the intramolecular motions. This is the concept of aggregation-induced emission (AIE). DSE fluorophores have been far less explored despite the fact that they are at the crossroad of ACQ and AIE phenomena and allow targeting applications both in solution (bio-conjugation, sensing, imaging) and solid-state (organic electronics, data encryption, lasing, luminescent displays). Excited-State Intramolecular Proton Transfer (ESIPT) fluorescence is particularly suitable to engineer DSE dyes. Indeed, ESIPT fluorescence, which relies on a phototautomerism between normal and tautomeric species, is characterized by a strong emission in the solid-state along with a large Stokes’ shift, an enhanced photostability and a strong sensitivity to the close environment, a feature prone to be used in bio-sensing. A drawback that needs to be overcome is their weak emission intensity in solution, owing to detrimental molecular motions in the excited-state. Several strategies have been proposed in that regard. In the past few years, a growing number of examples of DSE-ESIPT dyes have indeed emerged in the literature, enriching the database of such attractive dyes. This review aims at a brief but concise overview on the exploitation of ESIPT luminescence for the optimization of DSE dyes properties. In that perspective, a synergistic approach between organic synthesis, fluorescence spectroscopy and ab initio calculations has proven to be an efficient tool for the construction and optimization of DSE-ESIPT fluorophores.

Mechanistic Investigation on the Initial Thermal Decomposition of Energetic Materials FOX-7 and RDX in the Crystal and Gas Phase: An MM/DFT-Based ONIOM Calculation

Interpreting the initial decomposition mechanism is important for evaluating the thermal stability of explosives. In this study, we theoretically investigated the initial thermal decomposition reactions for two typical energetic materials, FOX-7 and RDX, in both the gas phase and crystal phase. Single molecular decomposition pathways in the gas phase are calculated using the density functional theory (DFT) method, and the crystal phase reactions are simulated through the MM/DFT-based ONIOM method. The calculation results indicate that the crystal environment has a significant influence on the initial thermal decomposition mechanism of FOX-7 and RDX. The cage effect induced by the crystal environment greatly confines molecular mobility and diffusion, rendering the generated small molecules to react with the remaining fragment and yield new decomposition channels compared with the gas phase condition. The crystal packing structures and intermolecular interactions (hydrogen bonds/π-π stacking) significantly increase the reaction barriers of FOX-7 and RDX, leading to the crystal phase reactions being more difficult to occur than in the gas phase. Since the practical application of explosives is mostly in the crystal state, it is important to consider the environmental effects on the initial decomposition reactions. The same insight can also be relevant for other energetic materials.

Emission Quenching in Tetraphenylfuran Crystal: Why This Propeller-Shaped Molecule Does Not Emit in the Condensed Phase

Due to their substantial fluorescence quantum yields in the crystalline phase, propeller-shaped molecules have recently gained significant attention as potential emissive materials for optoelectronic applications. For the family of cyclopentadiene derivatives, light-emission is highly dependent on the nature of heteroatomic substitutions. In this paper, we investigate excited state relaxation pathways in the tetraphenyl-furan molecule (TPF), which in contrast with other molecules in the family, shows emission quenching in the solid-state. For the singlet manifold, our calculations show nonradiative pathways associated with C-O elongation are blocked in both vacuum and the solid state. A fraction of the population can be transferred to the triplet manifold and, subsequently, to the ground state in both phases. This process is expected to be relatively slow due to the small spin-orbit couplings between the relevant singlet-triplet states. Emission quenching in crystalline TPF seems to be in line with more efficient exciton hopping rates. Our simulations help clarify the role of conical intersections, population of the triplet states and crystalline structure in the emissive response of propeller-shaped molecules.

A quantum crystallographic approach to short hydrogen bonds

In this work we use high-resolution synchrotron X-ray diffraction for electron density mapping, in conjunction with ab initio modelling, to study short O-H⋯O and O+-H⋯O- hydrogen bonds whose behaviour is known to be tuneable by temperature. The short hydrogen bonds have donor-acceptor distances in the region of 2.45 Å and are formed in substituted urea and organic acid molecular complexes of N,N'-dimethylurea oxalic acid 2 : 1 (1), N,N-dimethylurea 2,4-dinitrobenzoate 1 : 1 (2) and N,N-dimethylurea 3,5-dinitrobenzoic acid 2 : 2 (3). From the combined analyses, these complexes are found to fall within the salt-cocrystal continuum and exhibit short hydrogen bonds that can be characterised as both strong and electrostatic (1, 3) or very strong with a significant covalent contribution (2). An additional charge assisted component is found to be important in distinguishing the relatively uncommon O-H⋯O pseudo-covalent interaction from a typical strong hydrogen bond. The electron density is found to be sensitive to the extent of static proton transfer, presenting it as a useful parameter in the study of the salt-cocrystal continuum. From complementary calculated hydrogen atom potentials, we attribute changes in proton position to the molecular environment. Calculated potentials also show zero barrier to proton migration, forming an 'energy slide' between the donor and acceptor atoms. The better fundamental understanding of the short hydrogen bond in the 'zone of fluctuation' presented in a salt-cocrystal continuum, enabled by studies like this, provide greater insight into their related properties and can have implications in the regulation of pharmaceutical materials.

Role of Conical Intersections on the Efficiency of Fluorescent Organic Molecular Crystals

Organic molecular crystals are attractive materials for luminescent applications because of their promised tunability. However, the link between the chemical structure and emissive behavior is poorly understood because of the numerous interconnected factors which are at play in determining radiative and nonradiative behaviors at the solid-state level. In particular, the decay through conical intersection dominates the nonadiabatic regions of the potential energy surface, and thus, their accessibility is a telling indicator of the luminosity of the material. In this study, we investigate the radiative mechanism for five organic molecular crystals which display a solid-state emission, with a focus on the role of conical intersections in their photomechanisms. The objective is to situate the importance of the accessibility of conical intersections with regards to emissive behavior, taking into account other nonradiative decay channels, namely, vibrational decay, and exciton hopping. We begin by giving a brief overview of the structural patterns of the five systems within a larger pool of 13 crystals for a richer comparison. We observe that because of the prevalence of sheet like and herringbone packing in organic molecular crystals, the conformational diversity of crystal dimers is limited. Additionally, similarly spaced dimers have exciton coupling values of a similar order within a 50 meV interval. Next, we focus on three exemplary cases, where we disentangle the role of nonradiative decay mechanisms and show how rotational minimum energy conical intersections in vacuum lead to puckered ones in the crystal, increasing their instability upon crystallization in typical packing motifs. In contrast, molecules with puckered conical intersections in vacuum tend to conserve this trait upon crystallization, and therefore, their quantum yield of fluorescence is determined predominantly by other nonradiative decay mechanisms.

The Mechanics of the Bicycle Pedal Photoisomerization in Crystalline cis,cis-1,4-Diphenyl-1,3-butadiene

Direct irradiation of crystalline cis,cis-1,4-diphenyl-1,3-butadiene (cc-DPB) forms trans,trans-1,4-diphenyl-1,3,-butadiene via a concerted two-bond isomerization called the bicycle pedal (BP) mechanism. However, little is known about photoisomerization pathways in the solid state and there has been much debate surrounding the interpretation of volume-conserving isomerization mechanisms. The bicycle pedal photoisomerization is investigated using the quantum mechanics/molecular mechanics complete active space self-consistent field/Amber force-field method. Important details about how the steric environment influences isomerization mechanisms are revealed including how the one-bond flip and hula-twist mechanisms are suppressed by the crystal cavity, the nature of the seam space in steric environments, and the features of the bicycle pedal mechanism. Specifically, in the bicycle pedal, the phenyl rings of cc-DPB are locked in place and the intermolecular packing allows a passageway for rotation of the central diene in a volume-conserving manner. In contrast, the bicycle pedal rotation in the gas phase is not a stable pathway, so single-bond rotation mechanisms become operative instead. Furthermore, the crystal BP mechanism is an activated process that occurs completely on the excited state; the photoproduct can decay to the ground state through radiative and non-radiative pathways. The present models, however, do not capture the quantitative activation barriers, and more work is needed to better model reactions in crystals. Last, the reaction barriers of the different crystalline conformations within the unit cell of cc-DPB are compared to investigate the possibility for conformation-dependent isomerization. Although some difference in reaction barriers is observed, the difference is most likely not responsible for the experimentally observed periods of fast and slow conversion.

fromage: A library for the study of molecular crystal excited states at the aggregate scale

The study of photoexcitations in molecular aggregates faces the twofold problem of the increased computational cost associated with excited states and the complexity of the interactions among the constituent monomers. A mechanistic investigation of these processes requires the analysis of the intermolecular interactions, the effect of the environment, and 3D arrangements or crystal packing on the excited states. A considerable number of techniques have been tailored to navigate these obstacles; however, they are usually restricted to in‐house codes and thus require a disproportionate effort to adopt by researchers approaching the field. Herein, we present the FRamewOrk for Molecular AGgregate Excitations (fromage), which implements a collection of such techniques in a Python library complemented with ready‐to‐use scripts. The program structure is presented and the principal features available to the user are described: geometrical analysis, exciton characterization, and a variety of ONIOM schemes. Each is illustrated by examples of diverse organic molecules in condensed phase settings. The program is available at https://github.com/Crespo-Otero-group/fromage.