Modulating Efficiency and Color of Thermally Activated Delayed Fluorescence by Rationalizing the Substitution Effect

Thermally activated delayed fluorescence (TADF) constitutes the process by which third-generation organic light-emitting diodes (OLEDs) are being designed and produced. Despite several years of trial-and-error attempts, mainly driven by chemical intuition about how to improve a certain aspect of the process, few studies focused on the in-depth description of its two key properties: efficiency of the T1 → S1 intersystem crossing and further S1 → S0 emission. Here, by means of a newly developed theoretical formalism, we propose a systematic rationalization of the substituent effect in a paradigmatic class of OLED compounds, based on phenothiazine-dibenzothiophene-S,S-dioxide, known as PTZ-DBTO2. Our methodology allows to discern among geometrical and electronic effects induced by the substituent, deeply understanding the relationships existing between charge transfer, spin density, geometrical deformations, and energy modulations between electronic states. By our results, we can finally elucidate, depending on the substituent, the fate of the overall TADF process, quantitatively assessing its efficiency and predicting the color emission. Moreover, the general terms by which this methodology was developed allow its application to any chromophore of interest.

An Algorithm Predicting the Optimal Mechanical Response of Electronic Energy Difference

The use of mechanical forces at the molecular level has been shown to be an interesting tool for modulating different chemical and physical molecular properties. The so-called covalent mechanochemistry deals with the application of precise mechanical forces that induce specific changes in the structure, stability, reactivity, and other physical properties. The use of this kind of force to modulate photophysical properties and photochemical reactivity has also been studied. Nevertheless, the general problem of mechanical modulation of the energy gap between two electronic states has been addressed only with the development of simple theoretical models. Here, we develop and implement an algorithm providing the Largest energy Gap variation with Minimal mechanical Force (LGMF) that allows the determination of the optimal mechanical forces tuning the electronic energy gap, as well as to identify the maximum mechanical response of a molecular system to the application of any mechanical stimulus. The algorithm has been implemented for diverse molecular systems showing different degrees of flexibility. The phyton code of the algorithm is available in a public repository.

Mechanical Activation of Forbidden Photoreactivity in Oxa-di-π-methane Rearrangement

In this work, we demonstrate that the forbidden oxirane-type photoproduct (the cyclopropyl ketone photoproduct is the allowed one) of the oxa-di-π-methane photorearrangement can be obtained by mechanochemical control of the photoreactions. This control is achieved by the application of simple force pairs rationally chosen. By analyzing in detail the effect of the applied forces on this photoreaction, it comes to light that the mechanical action affects the diverse properties of the oxa-di-π-methane rearrangement, modifying all the steps of the reaction: (i) the initial ground-state conformers’ distribution becomes affected; (ii) the new conformational population makes the triplet excitation process to be changed, responding to the magnitude of the applied force; (iii) the stability of the different intermediates along the triplet pathway also becomes affected, changing the dynamical behavior of the system and the reaction kinetics; and (iv) the intersystem crossing also becomes strongly affected, making the forbidden oxirane-type photoproduct to decay more efficiently to the ground state. All these changes provide a complex scenario where a detailed study of the effect of applied forces is necessary in order to predict its overall effect on the photoreactivity.

Thermal and Mechanochemical Tuning of the Porphyrin Singlet-Triplet Gap for Selective Energy Transfer Processes: A Molecular Dynamics Approach

Molecular dynamics simulations provide fundamental knowledge on the reaction mechanism of a given simulated molecular process. Nevertheless, other methodologies based on the "static" exploration of potential energy surfaces are usually employed to firmly provide the reaction coordinate directly related to the reaction mechanism, as is the case in intrinsic reaction coordinates for thermally activated reactions. Photoinduced processes in molecular systems can also be studied with these two strategies, as is the case in the triplet energy transfer process. Triplet energy transfer is a fundamental photophysical process in photochemistry and photobiology, being for instance involved in photodynamic therapy, when generating the highly reactive singlet oxygen species. Here, we study the triplet energy transfer process between porphyrin, a prototypical energy transfer donor, and different biologically relevant acceptors, including molecular oxygen, carotenoids, and rhodopsin. The results obtained by means of nanosecond time-scale molecular dynamics simulations are compared to the "static" determination of the reaction coordinate for such a thermal process, leading to the distortions determining an effective energy transfer. This knowledge was finally applied to propose porphyrin derivatives for producing the required structural modifications in order to tune their singlet-triplet energy gap, thus introducing a mechanochemical description of the mechanism.

The concept of substituent-induced force in the rationale of substituent effect

Controlling the thermochemistry and kinetics of chemical reactions is a central problem in chemistry. Among factors permitting this control, the substituent effect constitutes a remarkable example. Here, we develop a model accounting for the effect of a substituent on the potential energy surface of the substrate (i.e., substituted molecule). We show that substituents affect the substrate by exerting forces on the nuclei. These substituent-induced forces are able to develop a work when the molecule follows a given reaction path. By applying a simple mechanical model, it becomes possible to quantify this work, which corresponds to the energy variation due to the effect of the substituent along a specific pathway. Our model accounts for the Hammett equation as a particular case, providing the first non-empirical scale for the σ and ρ constants, which, in the developed model, are related to the forces exerted by the substituents (σ) and the reaction path length (ρ), giving their product (σ · ρ) the well-known variation on the reaction energy due to the substituent.

How mechanical forces can modulate the metal affinity and selectivity of metal binding sites in proteins

The results obtained reveal that applying mechanical forces with a given strength and directionality can modulate the metal affinity and selectivity of metal binding sites in metalloproteins.

Photoreactivity Control Mediated by Molecular Force Probes in Stilbene

We report theoretical and experimental evidence showing that photochemical reactivity of a chromophore can be modified by applying mechanical forces via molecular force probes. This mechanical action permits us to modulate main photochemical properties, such as fluorescence yield, excited-state lifetime, or photoisomerization quantum yield. The effect of molecular force probes can be rationalized in terms of simple mechanochemical models, establishing a qualitative framework for understanding the mechanical control of photoreactivity in stilbenes.

C-H Functionalization of BN-Aromatics Promoted by Addition of Organolithium Compounds to the Boron Atom

Addition of an organolithium compound to a BN-phenanthrene with embedded B and N atoms is proposed to result in coordination of RLi to the boron atom. This coordination, supported by NMR spectroscopy and DFT calculations, increases the nucleophilicity of the system in the β position to the N atom and is therefore a useful tool for promoting regioselective C-H functionalization of BN aromatics.

How Do Methyl Groups Enhance the Triplet Chemiexcitation Yield of Dioxetane?

Chemiluminescence is the emission of light as a result of a nonadiabatic chemical reaction. The present work is concerned with understanding the yield of chemiluminescence, in particular how it dramatically increases upon methylation of 1,2-dioxetane. Both ground-state and nonadiabatic dynamics (including singlet excited states) of the decomposition reaction of various methyl-substituted dioxetanes have been simulated. Methyl-substitution leads to a significant increase in the dissociation time scale. The rotation around the O-C-C-O dihedral angle is slowed; thus, the molecular system stays longer in the "entropic trap" region. A simple kinetic model is proposed to explain how this leads to a higher chemiluminescence yield. These results have important implications for the design of efficient chemiluminescent systems in medical, environmental, and industrial applications.

Synthesis, Optical Properties, and Regioselective Functionalization of 4a-Aza-10a-boraphenanthrene

4a-Aza-10a-boraphenanthrene has been synthesized in only four steps from commercially available materials with a remarkable overall yield of 62%. In contrast to other BN-isosteres of phenathrene, this isomer is weakly fluorescent, which has been explained by means of computational studies that found a low energy conical intersection for the nonradiative deactivation of the excited state. Moreover, a completely regioselective functionalization of 4a-aza-10a-boraphenanthrene at C-1 by reaction with activated electrophiles has been achieved.

Study of Model Systems for Bilirubin and Bilin Chromophores: Determination and Modification of Thermal and Photochemical Properties

Bilin chromophores and bilirubin are involved in relevant biological functions such as light perception in plants and as protective agents against Alzheimer and other diseases. Despite their extensive use, a deep rationalization of the main factors controlling the thermal and photochemical properties has not been performed yet, which in turn hampers further applications of these versatile molecules. In an effort to understand those factors and allow control of the relevant properties, a combined experimental and computational study has been carried out for diverse model systems to understand the interconversion between Z and E isomers. In this study, we have demonstrated the crucial role of steric hindrance and hydrogen-bond interactions in thermal stability and the ability to control them by designing novel compounds. We also determined several photochemical properties and studied the photodynamics of two model systems in more detail, observing a fast relaxation of the excited state shorter than 2 ps in both cases. Finally, the computational study allowed us to rationalize the experimental evidence.

Probing the Photodynamics of Rhodopsins with Reduced Retinal Chromophores

While the light-induced population dynamics of different photoresponsive proteins has been investigated spectroscopically, systematic computational studies have not yet been possible due to the phenomenally high cost of suitable high level quantum chemical methods and the need of propagating hundreds, if not thousands, of nonadiabatic trajectories. Here we explore the possibility of studying the photodynamics of rhodopsins by constructing and investigating quantum mechanics/molecular mechanics (QM/MM) models featuring reduced retinal chromophores. In order to do so we use the sensory rhodopsin found in the cyanobacterium Anabaena PCC7120 (ASR) as a benchmark system. We find that the basic mechanistic features associated with the excited state dynamics of ASR QM/MM models are reproduced using models incorporating a minimal (i.e., three double-bond) chromophore. Furthermore, we show that ensembles of nonadiabatic ASR trajectories computed using the same abridged models replicate, at both the CASPT2 and CASSCF levels of theory, the trends in spectroscopy and lifetimes estimated using unabridged models and observed experimentally at room temperature. We conclude that a further expansion of these studies may lead to low-cost QM/MM rhodopsin models that may be used as effective tools in high-throughput in silico mutant screening.

A biomimetic molecular switch at work: coupling photoisomerization dynamics to peptide structural rearrangement

A combined experimental and computational study of a peptide-linked retinal-like molecular switch shows the effects on photoreactivity and the α-helix structure.

Molcas 8: New capabilities for multiconfigurational quantum chemical calculations across the periodic table

In this report, we summarize and describe the recent unique updates and additions to the Molcas quantum chemistry program suite as contained in release version 8. These updates include natural and spin orbitals for studies of magnetic properties, local and linear scaling methods for the Douglas-Kroll-Hess transformation, the generalized active space concept in MCSCF methods, a combination of multiconfigurational wave functions with density functional theory in the MC-PDFT method, additional methods for computation of magnetic properties, methods for diabatization, analytical gradients of state average complete active space SCF in association with density fitting, methods for constrained fragment optimization, large-scale parallel multireference configuration interaction including analytic gradients via the interface to the Columbus package, and approximations of the CASPT2 method to be used for computations of large systems. In addition, the report includes the description of a computational machinery for nonlinear optical spectroscopy through an interface to the QM/MM package Cobramm. Further, a module to run molecular dynamics simulations is added, two surface hopping algorithms are included to enable nonadiabatic calculations, and the DQ method for diabatization is added. Finally, we report on the subject of improvements with respects to alternative file options and parallelization.

Hydantoin-based molecular photoswitches

A new family of molecular photoswitches based on arylidenehydantoins is described together with their synthesis and photochemical and photophysical studies. A series of hydantoin derivatives have been prepared as single isomers using simple and versatile chemistry in good yields. Our studies show that the photostationary states of these compounds can be easily controlled by means of external factors, such as the light source or filters. Moreover, the detailed investigations proved that these switches are efficient (i.e., they make efficient use of the light energy, are high fatigue resistant, and are very photostable). In some cases, the switches can be completely turned on/off, a desirable feature for specific applications. A series of theoretical calculations have also been carried out to understand the photoisomerization mechanism at the molecular level.

Directionality of Double-Bond Photoisomerization Dynamics Induced by a Single Stereogenic Center

In light-driven single-molecule rotary motors, the photoisomerization of a double bond converts light energy into the rotation of a moiety (the rotor) with respect to another (the stator). However, at the level of a molecular population, an effective rotary motion can only be achieved if a large majority of the rotors rotate in the same, specific direction. Here we present a quantitative investigation of the directionality (clockwise vs counterclockwise) induced by a single stereogenic center placed in allylic position with respect to the reactive double bond of a model of the biomimetic indanylidene-pyrrolinium framework. By computing ensembles of nonadiabatic trajectories at 300 K, we predict that the photoisomerization is >70% unidirectional for the Z → E and E → Z conversions. Most importantly, we show that such directionality, resulting from the asymmetry of the excited state force field, can still be observed in the presence of a small (ca. 2°) pretwist or helicity of the reactive double bond. This questions the validity of the conjecture that a significant double-bond pretwist (e.g., >10°) in the ground state equilibrium structure of synthetic or natural rotary motors would be required for unidirectional motion.

Shape of Multireference, Equation-of-Motion Coupled-Cluster, and Density Functional Theory Potential Energy Surfaces at a Conical Intersection

We report and characterize ground-state and excited-state potential energy profiles using a variety of electronic structure methods along a loop lying on the branching plane associated with a conical intersection (CI) of a reduced retinal model, the penta-2,4-dieniminium cation (PSB3). Whereas the performance of the equation-of-motion coupled-cluster, density functional theory, and multireference methods had been tested along the excited- and ground-state paths of PSB3 in our earlier work, the ability of these methods to correctly describe the potential energy surface shape along a CI branching plane has not yet been investigated. This is the focus of the present contribution. We find, in agreement with earlier studies by others, that standard time-dependent DFT (TDDFT) does not yield the correct two-dimensional (i.e., conical) crossing along the branching plane but rather a one-dimensional (i.e., linear) crossing along the same plane. The same type of behavior is found for SS-CASPT2(IPEA=0), SS-CASPT2(IPEA=0.25), spin-projected SF-TDDFT, EOM-SF-CCSD, and, finally, for the reference MRCISD+Q method. In contrast, we found that MRCISD, CASSCF, MS-CASPT2(IPEA=0), MS-CASPT2(IPEA=0.25), XMCQDPT2, QD-NEVPT2, non-spin-projected SF-TDDFT, and SI-SA-REKS yield the expected conical crossing. To assess the effect of the different crossing topologies (i.e., linear or conical) on the PSB3 photoisomerization efficiency, we discuss the results of 100 semiclassical trajectories computed by CASSCF and SS-CASPT2(IPEA=0.25) for a PSB3 derivative. We show that for the same initial conditions, the two methods yield similar dynamics leading to isomerization quantum yields that differ by only a few percent.

Chiral Hydrogen Bond Environment Providing Unidirectional Rotation in Photoactive Molecular Motors

Generation of a chiral hydrogen bond environment in efficient molecular photoswitches is proposed as a novel strategy for the design of photoactive molecular motors. Here, the following strategy is used to design a retinal-based motor presenting singular properties: (i) a single excitation wavelength is needed to complete the unidirectional rotation process (360°); (ii) the absence of any thermal step permits the process to take place at low temperatures; and (iii) the ultrafast process permits high rotational frequencies.

First principles study of photostability within hydrogen-bonded amino acids

The photochemistry and photophysics of a two-glycine minimal model is studied at the CASPT2//CASSCF level of theory. Different photoinduced processes are discussed, on the basis of the calculated minimum energy paths and the characterization of the electronic state crossings. Two main processes could provide UV-photostability to the hydrogen-bonded peptide system: (i) forward-backward photoinduced electron/proton transfer involving the H in the hydrogen bond, (ii) singlet-singlet energy transfer between two amino acids, providing ultrafast population of the low-energy n,π* state.

Modeling, preparation, and characterization of a dipole moment switch driven by Z/E photoisomerization

We report the results of a multidisciplinary research effort where the methods of computational photochemistry and retrosynthetic analysis/synthesis have contributed to the preparation of a novel N-alkylated indanylidene-pyrroline Schiff base featuring an exocyclic double bond and a permanent zwitterionic head. We show that, due to its large dipole moment and efficient photoisomerization, such a system may constitute the prototype of a novel generation of electrostatic switches achieving a reversible light-induced dipole moment change on the order of 30 D. The modeling of a peptide fragment incorporating the zwitterionic head into a conformationally rigid side chain shows that the switch can effectively modulate the fluorescence of a tryptophan probe.

Photoinduced electron and proton transfer in the hydrogen-bonded pyridine-pyrrole system

We present here a detailed analysis of the mechanism of photoinduced electron and proton transfer in the planar pyrrole-pyridine hydrogen-bonded system, a model for the photochemistry of hydrogen bonds in DNA base pairs. Two different crossings, an avoided crossing and a conical intersection, are the key steps for forward and backward electron and proton transfer providing to the system photostability against UV radiation by restoring the system in its initial electronic and geometric structure.

Tracking the excited-state time evolution of the visual pigment with multiconfigurational quantum chemistry

The primary event that initiates vision is the photoinduced isomerization of retinal in the visual pigment rhodopsin (Rh). Here, we use a scaled quantum mechanics/molecular mechanics potential that reproduces the isomerization path determined with multiconfigurational perturbation theory to follow the excited-state evolution of bovine Rh. The analysis of a 140-fs trajectory provides a description of the electronic and geometrical changes that prepare the system for decay to the ground state. The data uncover a complex change of the retinal backbone that, at approximately 60-fs delay, initiates a space saving "asynchronous bicycle-pedal or crankshaft" motion, leading to a conical intersection on a 110-fs time scale. It is shown that the twisted structure achieved at decay features a momentum that provides a natural route toward the photoRh structure recently resolved by using femtosecond-stimulated Raman spectroscopy.

The role of the intersection space in the photochemistry of tricyclo[3.3.0.0(2,6)]octa-3,7-diene

CASSCF and CASPT2 methods were used to study the photochemistry of tricyclo[3.3.0.0(2,6)]octa-3,7-diene (TOD). The analysis of different S1 reaction paths as well as the topology of the S1/S0 intersection space allows us to establish two novel properties associated with the photochemical behavior of this compound: (i) simple low-lying intersection space domains can mediate different photoproducts, and (ii) TOD photochemistry is probably mediated by two disconnected intersection space domains, related to the formation of cyclooctatetraene and semibulvalene in different time-scales. It is shown that these domains are chemically distinct since the first, leading to COT, mediates barrierless pericyclic reactions while the second, leading to SBV, is accessed through the formation of an excited-state biradical intermediate. To the best of our knowledge, in the domain of single molecule photochemistry, TOD represents the first example where a different chemical role of distinct low-lying intersection spaces has been computationally documented. The observed photoproducts can be rationalized in terms of branching space diagrams, constructed by determining the branching space (derivative coupling and gradient difference vectors) for each conical intersection involved in the photochemical process.

A theory of nonvertical triplet energy transfer in terms of accurate potential energy surfaces: the transfer reaction from pi,pi* triplet donors to 1,3,5,7-cyclooctatetraene

Triplet energy transfer (TET) from aromatic donors to 1,3,5,7-cyclooctatetraene (COT) is an extreme case of "nonvertical" behavior, where the transfer rate for low-energy donors is considerably faster than that predicted for a thermally activated (Arrhenius) process. To explain the anomalous TET of COT and other molecules, a new theoretical model based on transition state theory for nonadiabatic processes is proposed here, which makes use of the adiabatic potential energy surfaces (PES) of reactants and products, as computed from high-level quantum mechanical methods, and a nonadiabatic transfer rate constant. It is shown that the rate of transfer depends on a geometrical distortion parameter gamma=(2g(2)/kappa(1))(1/2) in which g stands for the norm of the energy gradient in the PES of the acceptor triplet state and kappa(1) is a combination of vibrational force constants of the ground-state acceptor in the gradient direction. The application of the model to existing experimental data for the triplet energy transfer reaction to COT from a series of pi,pi(*) triplet donors, provides a detailed interpretation of the parameters that determine the transfer rate constant. In addition, the model shows that the observed decrease of the acceptor electronic excitation energy is due to thermal activation of C=C bond stretchings and C-C bond torsions, which collectively change the ground-state COT bent conformation (D(2d)) toward a planar triplet state (D(8h)).

Correlated MO Study of the Low-Barrier Intramolecular Motions in Donor−Acceptor Ethenes

Correlated MP2 and MCSCF MO calculations of several model push-pull ethenes in most cases indicate no great participation of excited singlet and triplet electronic configurations in either the minima or the transition structures for the suggested facilitated intramolecular rotation about the polarized C=C bond. This situation changes significantly only in molecules with sulfur atoms in the molecule as either donors or acceptors. The outstanding contribution of sulfur atoms as either donors or acceptors is a significant increase of push-pull ethene molecular polarizabilities. Thus, within the studied small series of mostly planar push-pull ethenes, polarizability appears a better indicator of rapid intramolecular motions about the C=C bond than straight polarity. Substituents with larger steric demands around the C=C bond are shown to likely preclude its complete turnaround, thus prompting a ramification of the interpretations of dynamic NMR phenomena in sterically constrained push-pull ethenes as large-amplitude librations resulting from strong rovibrational and relatively weak electronic coupling. These librations, as shown by complete vibrational mode analysis of corresponding rotational transition structures, cover in fact certain sectors of the intuitively suggested full rotations similar to those about C-C single bonds.

A new algorithm for predicting triplet-triplet energy-transfer activated complex coordinate in terms of accurate potential-energy surfaces

The new algorithm presented here allows, for the first time, the determination of the optimal geometrical distortions that an acceptor molecule in the triplet-triplet energy-transfer process undergoes, as well as the dependence of the activation energy of the process on the triplet energy difference of donor and acceptor molecules. This algorithm makes use of the complete potential-energy surfaces (singlet and triplet states), and contrasts with the first-order approximation already published [L. M. Frutos, O. Castano, J. L. Andres, M. Merchan, and A. U. Acuna, J. Chem. Phys. 120, 1208 (2004)] in which an expansion of the potential-energy surfaces was used. This algorithm is gradient based and finds the best trajectory for the acceptor molecule, starting from S(0) ground-state equilibrium geometry, to achieve the maximum variation of the singlet-triplet energy gap with the minimum energy of activation on S(0). Therefore, the algorithm allows the determination of a "reaction path" for the triplet-triplet energy-transfer processes. Also, the algorithm could also serve eventually to find minimum-energy crossing (singlet-triplet) points on the potential-energy surface, which can play an important role in the intersystem crossing process for the acceptor molecules to recover their initial capacity as acceptors. Also addressed is the misleading use of minimum-energy paths in T(1) to describe the energy-transfer process by comparing these results with those obtained using the new algorithm. The implementation of the algorithm is illustrated with different potential-energy surface models and it is discussed in the frame of nonvertical behavior.

Cyclooctatetraene computational photo- and thermal chemistry: a reactivity model for conjugated hydrocarbons

We use ab initio CASSCF and CASPT2 computations to construct the composite multistate relaxation path relevant to cycloocta-1,3,5,7-tetraene singlet photochemistry. The results show that an efficient population of the dark excited state (S(1)) takes place after ultrafast decay from the spectroscopic excited state (S(2)). A planar D(8)(h)-symmetric minimum represents the collecting point on S(1). Nonadiabatic transitions to S(0) appear to be controlled by two different tetraradical-type conical intersections, which are directly accessible from the S(1) minimum following specific excited-state reaction paths. The higher-energy conical intersection belongs to the same type of intersections previously documented in linear and cyclic conjugated hydrocarbons and features a triangular -(CH)(3)- kink. This point mediates both cis --> trans photoisomerization and cyclopropanation reactions. The lowest energy conical intersection has a boat-shaped structure. This intersection accounts for production of semibullvalene or for double-bond shifting. The mapping of both photochemical and thermal reaction paths (including also Cope rearrangements, valence isomerizations, ring inversions, and double-bond shifting) has allowed us to draw a comprehensive reactivity scheme for cyclooctatetraene, which rationalizes the experimental observations and documents the complex network of photochemical and thermal reaction path interconnections. The factors controlling the selection and accessibility of a number of conjugated hydrocarbon prototype conical intersections and ground-state relaxation channels are discussed.

Role of bifurcation in the bond shifting of cyclooctatetraene

The present study of the cyclooctatetraene potential energy surface shows the presence of a bifurcation (valley ridge inflection point) in the intrinsic reaction coordinate path between the two transition states of D(8h) and D(4h) symmetries. This result is of capital importance for the correct understanding of the bond shifting and ring inversion processes in this compound.