Principles of Aggregation‐Induced Emission: Design of Deactivation Pathways for Advanced AIEgens and Applications

Twenty years ago, the concept of aggregation‐induced emission (AIE) was proposed, and this unique luminescent property has attracted scientific interest ever since. However, AIE denominates only the phenomenon, while the details of its underlying guiding principles remain to be elucidated. This minireview discusses the basic principles of AIE based on our previous mechanistic study of the photophysical behavior of 9,10‐bis(N,N‐dialkylamino)anthracene ( BDAA ) and the corresponding mechanistic analysis by quantum chemical calculations. BDAA comprises an anthracene core and small electron donors, which allows the quantum chemical aspects of AIE to be discussed. The key factor for AIE is the control over the non‐radiative decay (deactivation) pathway, which can be visualized by considering the conical intersection (CI) on a potential energy surface. Controlling the conical intersection (CI) on the potential energy surface enables the separate formation of fluorescent (CI:high) and non‐fluorescent (CI:low) molecules [control of conical intersection accessibility ( CCIA )]. The novelty and originality of AIE in the field of photochemistry lies in the creation of functionality by design and in the active control over deactivation pathways. Moreover, we provide a new design strategy for AIE luminogens (AIEgens) and discuss selected examples.     

Ultrafast nonradiative decay of electronically excited States of malachite green: ab initio calculations

We have investigated the nonradiative deactivation process of malachite green in the singlet excited states, S(1) and S(2), by high-level ab initio quantum chemical calculations using the CASPT2//CASCF approach. The deactivation pathways connecting the Franck-Condon region and conical intersection regions are identified. The initial population in the S(1) state is on a flat surface and the relaxation involves a rotation of phenyl rings, which leads the molecule to reach the conical intersection between the S(1) and S(0) states, where it efficiently decays back to the ground state. There exists a small barrier connecting the Franck-Condon and conical intersection regions on the S(1) potential energy surface. The decay mechanism from the S(2) state also involves the twisting motion of phenyl rings. In contrast to the excitation to the S(1) state, the initial population is on a downhill ramp potential and the barrierless relaxation through the rotation of substituted phenyl rings is expected. During the course of relaxation, the molecule switches to the S(1) state at the conical intersection between S(2) and S(1), and then it decays back to the ground state through the intersection between S(1) and S(0). In relaxation from both S(1) and S(2), large distortion of phenyl rings is required for the ultrafast nonradiative decay to the ground state.     

Potential energy surface crossings and the mechanistic spectrum for intramolecular electron transfer in organic radical cations.

The structure of the potential energy surface for the intramolecular electron transfer (IET) of four different model radical cations has been determined by using reaction path mapping and conical intersection optimization at the ab initio CASSCF level of theory. We show that, remarkably, the calculated paths reside in regions of the ground-state energy surface whose structure can be understood in terms of the position and properties of a surface crossing between the ground and the first excited state of the reactant. Thus, in the norbornadiene radical cation and in an analogue compound formed by two cyclopentene units linked by a norbornyl bridge, IET proceeds along direct-overlap and super-exchange concerted paths, respectively, that are located far from a sloped conical intersection point and in a region where the excited-state and ground-state surfaces are well separated. A second potential energy surface structure has been documented for 1,2-diamino ethane radical cation and features two parallel concerted (direct) and stepwise (chemical) paths. In this case a peaked conical intersection is located between the two paths. Finally, a third type of energy surface is documented for the bismethyleneadamantane radical cation and occurs when there is, effectively, a seam of intersection points (not a conical intersection) which separates the reactant and product regions. Since the reaction path cannot avoid the intersection, IET can only occur nonadiabatically. These IET paths indicate that quite different IET mechanisms may operate in radical cations, revealing an unexpectedly enriched and flexible mechanistic spectrum. We show that the origin of each path can be analyzed and understood in terms of the one-dimensional Marcus-Hush model.     

Energetics of cytosine singlet excited-state decay paths--a difficult case for CASSCF and CASPT2

Three deactivation paths for singlet excited cytosine are calculated at the CASPT2//CASSCF (complete active space second-order perturbation//complete active space self-consistent field) level of theory, using extended active spaces that allow for a reliable characterization of the paths and their energies. The lowest energy path, with a barrier of approximately 0.1 eV, corresponds to torsion of the C5-C6 bond, and the decay takes place at a conical intersection analogous to the one found for ethylene and its derivatives. There is a further path with a low energy barrier of approximately 0.2 eV associated with the (n(N),pi*) state which could also be populated with a low energy excitation. The path associated with a conical intersection between the ground and (n(O),pi*) states is significantly higher in energy (> 1 eV). The presence of minima on the potential energy surface for the (n,pi*) states that could contribute to the biexponential decay found in the gas phase was investigated, but could not be established unequivocally.     

Algorithms for sampling a quantum microcanonical ensemble of harmonic oscillators at potential minima and conical intersections

Algorithms are presented for sampling quantum microcanonical ensembles for a potential energy minimum and for the conical intersection at the minimum energy crossing point of two coupled electronic states. These ensembles may be used to initialize trajectories for chemical dynamics simulations. The unimolecular dynamics of a microcanonical ensemble about a potential energy minimum may be compared with the dynamics predicted by quantum Rice-Ramsperger-Kassel-Marcus (RRKM) theory. If the dynamics is non-RRKM, it will be of particular interest to determine which states have particularly long lifetimes. Initializing a microcanonical ensemble for the electronically excited state at a conical intersection is a model for electronic nonadiabatic dynamics. The trajectory surface-hopping approach may be used to study the ensuing chemical dynamics. A strength of the model is that zero-point energy conditions are included for the initial nonadiabatic dynamics at the conical intersection.     

Efficient photochemical merocyanine-to-spiropyran ring closure mechanism through an extended conical intersection seam. A model CASSCF/CASPT2 study

A mechanism of the thermal and photochemical bleaching of merocyanine to spiropyran is proposed on the basis of CASSCF/CASPT2 calculations on the 6-(2-propenyliden)cyclohexadienone model system. Our results suggest that this photochemical transformation takes place in two steps. First, the initially pumped 1(pi-pi) S2 undergoes radiationless decay to 1(n-pi) S1 via an extended S2/S1 conical intersection seam that runs approximately parallel to the trans-to-cis isomerization coordinate, a few kilocalories per mole higher in energy. Thus, S2 --> S1 internal conversion is possible at all values of the S2 trans-to-cis reaction coordinate. Second, on the S1 potential energy surface, there is a barrierless ring closure reaction path from the S1 cis minimum that leads to a peaked S1/S0 conical intersection where the deactivation to the ground state takes place. The inertia of the moving nuclei then drives the system toward the ground-state minimum of the 2H-chromene product. Thus, the extended seam topology of the S2/S1 conical intersection and the coordinate of the branching space of the S1/S0 conical intersection are essential to explain the efficiency and high speed of this reaction.     

Ab initio study on deactivation pathways of excited 9H-guanine

The complete active space with second-order perturbation theory/complete active space self-consistent-field method was used to explore the nonradiative decay mechanism for excited 9H-guanine. On the 1pipi* (1L(a)) surface we determined a conical intersection (CI), labeled (S0pipi*)(CI), between the 1pipi* (1L(a)) excited state and the ground state, and a minimum, labeled (pipi*)min. For the 1pipi* (1L(a)) state, its probable deactivation path is to undergo a spontaneous relaxation to (pipi*)min first and then decay to the ground state through (S0pipi*)(CI), during which a small activation energy is required. On the 1n(N)pi* surface a CI between the 1n(N)pi* and 1pipi* (1L(a)) states was located, which suggests that the 1n(N)pi* excited state could transform to the 1pipi* (1L(a)) excited state first and then follow the deactivation path of the 1pipi* (1L(a)) state. This CI was also possibly involved in the nonradiative decay path of the second lowest 1pipi* (1L(b)) state. On the 1n(O)pi* surface a minimum was determined. The deactivation of the 1n(O)pi* state to the ground state was estimated to be energetically unfavorable. On the 1pisigma* surface, the dissociation of the N-H bond of the six-membered ring is difficult to occur due to a significant barrier.     

Optimizing conical intersections of solvated molecules: the combined spin-flip density functional theory/effective fragment potential method

Solvent effects on a potential energy surface crossing are investigated by optimizing a conical intersection (CI) in solution. To this end, the analytic energy gradient has been derived and implemented for the collinear spin-flip density functional theory (SFDFT) combined with the effective fragment potential (EFP) solvent model. The new method is applied to the azomethane-water cluster and the chromophore of green fluorescent protein in aqueous solution. These applications illustrate not only dramatic changes in the CI geometries but also strong stabilization of the CI in a polar solvent. Furthermore, the CI geometries obtained by the hybrid SFDFT/EFP scheme reproduce those by the full SFDFT, indicating that the SFDFT/EFP method is an efficient and promising approach for understanding nonadiabatic processes in solution.     

Carbene formation in its lower singlet state from photoexcited 3H-diazirine or diazomethane. A combined CASPT2 and ab initio direct dynamics trajectory study

The potential energy surfaces of the ground and valence excited states of both 3H-diazirine and diazomethane have been studied computationally by mean of the CASSCF method in conjunction with the cc-pVTZ basis set. The energies of the critical points found on such surfaces have been recomputed at the CASPT2/cc-pVTZ level. Additionally, ab initio direct dynamic trajectory calculations have been carried out on the S(1) and S(2) surfaces, starting each trajectory run at the region dominated by the conformational molecular rearrangement of diazomethane. It is found that both isomers are interconnected along a C(s)() reaction coordinate on each potential surface. Radiationless deactivation of the corresponding S(1) state of each isomer occurs through the same point on the surface, an S(1)/S(0) conical intersection. Thereafter, the system has enough energy to surmount the barrier which leads to dissociation products (CH(2) + N(2)) on S(0) state. Therefore, photoexcitation to S(1) state of either diazirine of diazomethane produces methylene in its lower singlet state on a very short time scale (ca. 100 fs). Furthermore, both isomers can generate excited singlet carbene when they are excited onto the S(2) surface; in this case, they lose the activation energy passing through another common S(2)/S(1) conical intersection and then proceed to dissociation into carbene and N(2) on the S(1) surface. For the special case of methylene, it rapidly experiences deexcitation to S(0) state.     

Restriction of Intramolecular Motion(RIM): Investigating AIE Mechanism from Experimental and Theoretical Studies

Mechanistic studies promote scientific development from phenomena to theories.Aggregation-induced emission(AIE),as an unusual photophysical phenomenon,builds the bridge between molecular science and aggregate mesoscience.With the twenty-year development of AIE,restriction of intramolecular motion(RIM)has been verified as the working mechanism of AIE effect.In this review,these mechanistic works about RIM are summarized from experimental and theoretical perspectives.Thereinto,the experimental studies are introduced from three parts:external rigidification,structural modification and structural characterization.In the theoretical part,calculations on the low-frequency motion of AIEgens have been performed to prove the RIM mechanism.By virtue of the theoretical calculations,some new mechanisms are proposed to supplement the RIM,such as restriction of access to conical intersection,suppression of Kasha transition,restriction of access to dark state,etc.It is foreseeable that the RIM mechanism will unify the photophysical theories for both molecules and aggregates,and inspire more progress in aggregate science.     

A three-state model for the photophysics of adenine

An ab initio theoretical study at the CASPT2 level is reported on minimum energy reaction paths, state minima, transition states, reaction barriers, and conical intersections on the potential energy hypersurfaces of two tautomers of adenine: 9H- and 7H-adenine. The obtained results led to a complete interpretation of the photophysics of adenine and derivatives, both under jet-cooled conditions and in solution, within a three-state model. The ultrafast subpicosecond fluorescence decay measured in adenine is attributed to the low-lying conical intersection (gs/pipi* La)(CI), reached from the initially populated 1(pipi* La) state along a path which is found to be barrierless only in 9H-adenine, while for the 7H tautomer the presence of an intermediate plateau corresponding to an NH2-twisted conformation may explain the absence of ultrafast decay in 7-substituted compounds. A secondary picosecond decay is assigned to a path involving switches towards two other states, 1(pipi* Lb) and 1(npi*), ultimately leading to another conical intersection with the ground state, (gs/npi*), with a perpendicular disposition of the amino group. The topology of the hypersurfaces and the state properties explain the absence of secondary decay in 9-substituted adenines in water in terms of the higher position of the 1(npi*) state and also that the 1(pipi* Lb) state of 7H-adenine is responsible for the observed fluorescence in water. A detailed discussion comparing recent experimental and theoretical findings is given. As for other nucleobases, the predominant role of a pipi*-type state in the ultrafast deactivation of adenine is confirmed.     

Ab initio molecular dynamics of excited-state intramolecular proton transfer around a three-state conical intersection in malonaldehyde

Excited-state potential energy surface (PES) characterization is carried out at the CASSCF and MRSDCI levels, followed by ab initio dynamics simulation of excited-state intramolecular proton transfer (ESIPT) on the S2(pipi*) state in malonaldehyde. The proton-transfer transition state lies close to an S2/S1 conical intersection, leading to substantial coupling of proton transfer with electronic relaxation. Proton exchange proceeds freely on S2, but its duration is limited by competition with twisting out of the molecular plane. This rotamerization pathway leads to an intersection of the three lowest singlet states, providing the first detailed report of ab initio dynamics around a three-state intersection (3SI). There is a significant energy barrier to ESIPT on S1, and further pyramidalization of the twisted structure leads to the minimal energy S1/S0 intersection and energetic terminal point of excited-state dynamics. Kinetics and additional mechanistic details of these pathways are discussed. Significant depletion of the spectroscopic state and recovery of the ground state is seen within the first 250 fs after photoexcitation.     

On the characterization of three state conical intersections using a group homomorphism approach: mapping the full N-5 dimensional seam space

A method for characterizing the degeneracy preserving seam space in the vicinity of a three state conical intersection is introduced. Second order degenerate perturbation theory is used to construct an approximately diabatic Hamiltonian whose eigenenergies and eigenstates accurately describe the vicinity of the three state conical intersection in its full dimensionality. The perturbative analysis enables the large number, 6(N(int)(N(int)+1)2), of unique second order parameters needed to construct this accurate Hamiltonian to be determined from ab initio data at a limited number of nuclear configurations, with (N(int)+10) being minimal. Using the minimum energy three state conical intersection of the pyrazolyl radical (N(int) = 18), the potential of this approach is illustrated. A Hamiltonian comprised of the ten characteristic (linear) parameters and over 1440 second order parameters is constructed and used to determine the locus of the conical intersection seam as well as to describe the 18 dimensional space in the vicinity of that point of intersection. Our results demonstrate the ability of this methodology to quantitatively reproduce the ab initio potential energy surfaces near a three state conical intersection.     

Water effect on the excited-state decay paths of singlet excited cytosine

Two low-energy deactivation paths for singlet excited cytosine, one through a S₁/S₀ conical intersection of the ethylene type, and one through a conical intersection that involves the (n_N, π^*) state, are calculated in the presence of water. Water is included explicitly for several cytosine monohydrates, and as a bulk solvent, and the calculations are carried out at the complete active space self-consistent field (CASSCF) and complete active space second order perturbation (CASPT2) levels of theory. The effect of water on the lowest-energy path through the ethylenic conical intersection is a lowering of the energy barrier. This is explained by stabilization of the excited state, which has zwitterionic character in the vicinity of the conical intersection due to its similarity with the conical intersection of ethylene. In contrast to this, the path that involves the (n_N, π^*) state is destabilized by hydrogen bonding, although the bulk solvent effect reduces the destabilization. Overall, this path should remain energetically accessible.     

A unified perspective on the hydrogen atom transfer and proton-coupled electron transfer mechanisms in terms of topographic features of the ground and excited potential energy surfaces as exemplified by the reaction between phenol and radicals

The relation between the hydrogen atom transfer (HAT) and proton-coupled electron transfer (PCET) mechanisms is discussed and is illustrated by multiconfigurational electronic structure calculations on the ArOH + R(*) --> ArO(*) + RH reactions. The key topographic features of the Born-Oppenheimer potential energy surfaces that determine the predominant reaction mechanism are the conical intersection seam of the two lowest states and reaction saddle points located on the shoulders of this seam. The saddle point corresponds to a crossing of two interacting valence bond states corresponding to the reactant and product bonding patterns, and the conical intersection corresponds to the noninteracting intersection of the same two diabatic states. The locations of mechanistically relevant conical intersection structures and relevant saddle point structures are presented for the reactions between phenol and the N- and O-centered radicals, (*)NH2 and (*)OOCH3. Points on the conical intersection of the ground doublet D0 and first excited doublet D1 states are found to be in close geometric and energetic proximity to the reaction saddle points. In such systems, either the HAT mechanism or both the HAT mechanism and the proton-coupled electron transfer (PCET) mechanism can take place, depending on the relative energetic accessibility of the reaction saddle points and the D0/D1 conical intersection seams. The discussion shows how the two mechanisms are related and how they blend into each other along intermediate reaction paths. The recognition that the saddle point governing the HAT mechanism is on the shoulder of the conical intersection governing the PCET mechanism is used to provide a unified view of the competition between the two mechanisms (and the blending of the two mechanisms) in terms of the prominent and connected features of the potential energy surface, namely the saddle point and the conical intersection. The character of the dual mechanism may be understood in terms of the dominant valence bond configurations of the intersecting states, which are zero-order approximations to the diabatic states.     

Exploring Potential Energy Surfaces for Aggregation‐Induced Emission—From Solution to Crystal

Aggregation‐induced emission (AIE) is a phenomenon where non‐luminescent compounds in solution become strongly luminescent in aggregate and solid phase. It provides a fertile ground for luminescent applications that has rapidly developed in the last 15 years. In this review, we focus on the contributions of theory and computations to understanding the molecular mechanism behind it. Starting from initial models, such as restriction of intramolecular rotations (RIR), and the calculation of non‐radiative rates with Fermi's golden rule (FGR), we center on studies of the global excited‐state potential energy surfaces that have provided the basis for the restricted access to a conical intersection (RACI) model. In this model, which has been shown to apply for a diverse group of AIEgens, the lack of fluorescence in solution comes from radiationless decay at a CI in solution that is hindered in the aggregate state. We also highlight how intermolecular interactions modulate the photophysics in the aggregate phase, in terms of fluorescence quantum yield and emission color.     

On the characterization of three state conical intersections: a quasianalytic theory using a group homomorphism approach

In this work, degenerate perturbation theory through second order is used to characterize the vicinity of a three state conical intersection. This report extends our recent demonstration that it is possible to describe the branching space (in which the degeneracy is lifted linearly) and seam space (in which the degeneracy is preserved) in the vicinity of a two state conical intersection using second order perturbation theory. The general analysis developed here is based on a group homomorphism approach. Second order perturbation theory, in conjunction with high quality ab initio electronic structure data, produces an approximately diabatic Hamiltonian whose eigenenergies and eigenstates can accurately describe the three adiabatic potential energy surfaces, the interstate derivative couplings, and the branching and seam spaces in their full dimensionality. The application of this approach to the minimum energy three state conical intersection of the pyrazolyl radical demonstrates the potential of this method. A Hamiltonian comprised of the ten characteristic (linear) parameters and over 300 second order parameters is constructed to describe the branching space associated with a point of conical intersection. The second order parameters are determined using data at only 30 points. In the vicinity of the conical intersection the energy and derivative couplings are well reproduced and the singularity in the derivative coupling is analyzed.     

Ultrafast deactivation of an excited cytosine-guanine base pair in DNA

Multiconfigurational ab initio calculations and QM/MM molecular dynamics simulations of a photoexcited cytosine-guanine base pair in both gas phase and embedded in the DNA provide detailed structural and dynamical insights into the ultrafast radiationless deactivation mechanism. Photon absorption promotes transfer of a proton from the guanine to the cytosine. This proton transfer is followed by an efficient radiationless decay of the excited state via an extended conical intersection seam. The optimization of the conical intersection revealed that it has an unusual topology, in that there is only one degeneracy-lifting coordinate. This is the central mechanistic feature for the decay both in vacuo and in the DNA. Radiationless decay occurs along an extended hyperline nearly parallel to the proton-transfer coordinate, indicating the proton transfer itself is not directly responsible for the deactivation. The seam is displaced from the minimum energy proton-transfer path along a skeletal deformation of the bases. Decay can thus occur anywhere along the single proton-transfer coordinate, accounting for the remarkably short excited-state lifetime of the Watson-Crick base pair. In vacuo, decay occurs after a complete proton transfer, whereas in DNA, decay can also occur much earlier. The origin of this effect lies in the temporal electrostatic stabilization of dipole in the charge-transfer state in DNA.     

Photophysics of intramolecularly hydrogen-bonded aromatic systems: ab initio exploration of the excited-state deactivation mechanisms of salicylic acid

Excited state reaction paths and the corresponding energy profiles of salicylic acid have been determined with the CC2 method, which is a simplified version of singles-and-doubles coupled cluster theory. At crucial points of the potential energy hypersurfaces, single-point energy calculations have been performed with the CASPT2 method (second-order perturbation theory based on the complete active space self-consistent field reference). Hydrogen transfer along the intramolecular hydrogen bond as well as torsion and pyramidization of the carboxy group have been identified as the most relevant photochemical reaction coordinates. The keto-type planar S(1) state reached by barrierless intramolecular hydrogen transfer represents a local minimum of the S(1) energy surface, which is separated by a very low barrier from a reaction path leading to a low-lying S(1)-S(0) conical intersection via torsion and pyramidization of the carboxy group. The S(1)-S(0) conical intersection, which occurs for perpendicular geometry of the carboxy group, is a pure biradical. From the conical intersection, a barrierless reaction path steers the system back to the two known minima of the S(0) potential energy surface (rotamer I, rotamer II). A novel structure, 7-oxa-bicyclo[4.2.0]octa-1(6),2,4-triene-8,8-diol, has been identified as a possible transient intermediate in the photophysics of salicylic acid.     

Ab initio studies on the photophysics of guanine tautomers: out-of-plane deformation and NH dissociation pathways to conical intersections

The radiationless decay mechanisms of the S1 excited states of the 7H-keto-amino, 7H-enol-amino, and 7H-keto-imino tautomers of guanine have been investigated with the CASPT2//CASSCF method. Out-of-plane deformation of the six-membered ring or the imino group as well as dissociation of NH bonds have been considered as photochemical pathways leading to conical intersections with the electronic ground state. It has been found that all three tautomers can reach S0-S1 conical intersections by out-of-plane deformation. However, only in the 7H-keto-amino tautomer the reaction path leading to the conical intersection is barrierless. This tautomer also has the lowest energy barrier for hydrogen detachment via the (1)pi sigma* state, whose potential energy surface intersects that of the (1)pi pi* state as well as that of the ground state. The other tautomers of guanine exhibit substantial energy barriers on their S1 potential energy surfaces with respect to both reaction mechanisms. These findings suggest that the 7H-keto-amino tautomer exhibits the shortest excited-state lifetime of the three tautomers due to particularly fast nonradiative deactivation processes through S0-S1 conical intersections. The computational results explain the remarkable observation that the energetically most stable 7H-keto-amino tautomer is missing in the resonant two-photon ionization spectrum of guanine in a supersonic jet. The results also explain that the energetically less stable 7H-enol-amino and 7H-keto-imino tautomers have longer excited-state lifetimes and are thus detectable by resonant two-photon ionization.