Handbook of Aggregation-Induced Emission, Volume 1. Группа авторов

2.4, k_IC is mainly ruled by three factors: nonadiabatic electronic coupling, transition energy, and the vibronic relaxation energy. Based on the electronic structure information obtained from first‐principle calculations, the radiative and nonradiative rate constants were calculated by solving Equations 2.1 and 2.3 through the TVCF rate theory in the MOMAP program [26]. The difference between the PESs of two electronic states is considered by Q_i = SQ_f + D, where S is the Duschinsky rotation matrix and D is the displacement vector, that is, Duschinsky rotation effect (DRE) [13, 14].

      2.2.2 Computational Details

      To unravel the AIE mechanism, the photophysical properties of compounds in both the isolated (dilute solution) and the aggregated states (in amorphous aggregate and crystal) need to be considered. Different models were built according to the specific environment of the studied systems.

For the dispersed molecule in dilute solution, the polarizable continuum model (PCM) model is chosen to include bulk solvent effect implicitly (see Figure 2.1a). The polarities of different solvents are represented by dielectric constants.

      In crystal, the hybrid quantum mechanics/molecular mechanics (QM/MM) model can be set up based on the experimental crystal structure to consider the environment effect. First, a large cluster is extracted from the crystal. To balance the accuracy and computational cost, the cluster is spatially divided into a region where the QM description is essential, with the remainder represented by MM classically (Figure 2.1b). The general Amber force field (GAFF) [27] or the united force field (UFF) was used for the MM treatment. In the QM part, density functional theory (DFT)/time‐dependent DFT (TDDFT) is chosen to deal with the properties of the ground and excited states, respectively. During the calculation, the QM molecule is active, and all molecules in the MM region are frozen. The electrostatic embedding scheme was adopted to consider the interactions between QM and MM regions [28, 29].

For AIEgens in amorphous aggregate, the molecular conformation is not available in the experiment. Therefore, the large time‐scale molecular dynamics (MD) simulations need to be performed in solution to obtain the conformations of amorphous aggregates of AIEgens. The topology of the AIEgens could be set up by GAFF. The initial conformations were generated by randomly placing tens or hundreds of AIEgens in cubic simulation boxes with enough edge length. Then, the randomly dispersed AIEgen system was solvated by specific solvents reported in the experiment. The energy minimization was firstly performed using the steepest descent algorithm, followed by a 5‐ns pre‐equilibrated MD simulation under the NPT ensemble. Finally, several long time‐scale (more than 20 ns) production MD simulations were performed with different initial velocities coupled by the velocity rescaling thermostat [30] and Parrinello–Rahman barostat [31]. The time step was 2 fs. The configurations were stored at a specific time interval (e.g. 20 ps) for data analysis and the following calculation of photophysical properties. Then, we randomly extracted equilibrated aggregates from MD trajectories and set up QM/MM models accordingly to calculate the optical properties. Different from crystal, the molecular packing of amorphous HPS aggregates is irregular, resulting in distinct local environments for each molecule in amorphous aggregate. For example, the molecule embedded in aggregate is much different from the exposed one. In addition, conformations of different frames or trajectories are also different. Thus, several MD conformations need to be extracted from MD trajectories to set up QM/MM models to consider various molecular packings in amorphous aggregate. In the QM/MM model, the solvents could be removed from the extracted MD conformations to reduce the computational cost. In each QM/MM model for amorphous aggregate, one molecule is chosen as the QM region, and the others are treated as the MM region (see Figure 2.1c). The other details of QM/MM calculations are similar to that in crystal.

Figure 2.1 (a) Setup of the PCM model for an isolated molecule in solution; setup of the QM/MM models for a molecule in (b) crystal and (c) amorphous aggregate. The molecule in green is the QM region and the others in silver are the MM region (taking HPS as an example).

The QM calculation in the PCM model can perform in Gaussian 09 [32], Turbomole 6.5 [33], NW‐Chem [34], and so on. QM/MM calculations could be carried out by the ONIOM model in Gaussian 16 [32] or through the ChemShell interface [35]. MD simulation can be done by the GROMACS package. The DRE, reorganization energy (λ), k_r, and k_ic can be calculated by the MOMAP program [26, 36].

2.3 Revealed AIE Mechanism

AIEgens are non‐ or low‐emissive in dilute solution but highly emissive in the aggregate phase. The direct characterization of the luminescent property is the fluorescence quantum efficiency (Φ_F), which is related to k_r and k_ic by the following equation: Φ_F = k_r/(k_r + k_ic). The competition between k_r and k_ic determines the emission behavior of AIEgens. As shown in Figure 2.2a, both PESs of the ground state (S0) and first excited state (S1) in the solution are very flat; many vibrational modes with high quantum numbers are activated at high temperatures. These vibrational modes are strongly coupled with the electron, generating large reorganization energy, and different vibrational modes are mixed seriously, the so‐called dramatic DRE. Both the behaviors produce more nonradiative decay channels from the excited state to the ground state and highly speed up the nonradiative decay rate, while in the solid state, the compact molecular packing and electrostatic interactions restrict the intramolecular vibrational motions, which cause more steep PESs. Accordingly, the active vibrational modes and their corresponding mixtures become less favor, and the couplings between the electron and vibration become weak. Consequently, the nonradiative decay channels are blocked or even are removed and thus the rate is decreased significantly. A lot of calculated results demonstrate that aggregation can restrict the couplings between the electron and a great variety of molecular vibration modes, that is, various kinds of nonradiative decay channels can be hindered, such as rotating (see Figure 2.2b) [37–40], stretching [41], bending [42], flipping [43], and twisting vibrations [9, 44]. In the following, we overview the revealed AIE mechanism related to each kind of vibrational normal mode.

      2.3.1 Rotating Vibrations of Intramolecular Aromatic Ring

Siloles (silacyclopentadines) are the first kind of AIEgens [4, 45, 46]. The first silole molecule where the Φ_F in both gas phase and crystal are quantitatively investigated is HPS (Figure 2.3) [37]. HPS is poorly luminescent in cyclohexane at room temperature with a Φ_F as low as 0.30%, while in thin film, the Φ_F reaches 78% [45]. Combining QM and QM/MM calculations, the structural change analysis of optimized structures between S₀ and S₁ reflects that the torsional angles between the central silacycle and phenyl groups at the 2,5‐positions show larger modifications (14.56° and 14.63°) in the isolated state than those (5.85° and 1.05°) in solid phase, suggesting that the geometric relaxations of the phenyl rings at the 2,5‐positions are largely hindered in the solid phase. Moreover, the dihedral angle between the central silacycle