Theoretical Study of Internal Conversion Decay Rates Associated with the Three Lowest Singlet Electronic States in Pyrazine

Intramolecular rate-constant calculations based on the correlation function using temperature dependent quantum Green's functions

A theoretical method for calculating rate constants for internal conversion (IC), intersystem crossing (ISC) and radiative (R) electronic transitions is presented. The employed method uses temperature-dependent quantum Green's functions, which give the opportunity to consider almost any nth-order polynomial perturbation operator and the influence of external electromagnetic fields on the rate constants. The rate constants of the IC, ISC and R processes are calculated for two important indocyanine molecules namely indocyanine green (ICG) and heptamethine cyanine (IR808) at the Franck-Condon level using the temperature-dependent quantum Green's function approach. Calculations at the time-dependent density functional theory level with the MN15 functional show that ICG and IR808 have only one triplet state below the S1 state. The main deactivation channel of the S1 state is the IC process with a large (kIC(S1 → S0)) rate constant of ∼109-1011 s-1. The estimated quantum yield of fluorescence (φfl) is ∼0.001-0.24 for the two studied molecules, which agrees rather well with experimental values. Thus, the present approach enables calculations of the three kinds of rate constants and the quantum yield of fluorescence using the same computational methodology.

Mechanism of Enhanced Triplet-Triplet Upconversion in Organic Molecules

We use time-dependent density functional theory (TDDFT) to investigate the mechanism of efficient triplet-triplet upconversion (TTU) in certain organic materials. In particular, we focus on materials where some singlets are generated in a two-step spin-nonconserving process (T1 + T1 → T2 → S1). For this mechanism to contribute significantly, the intersystem crossing (ISC) from the high-lying triplet to the singlet (T2 → S1) must outcompete the internal conversion (IC) to the low-lying triplet (T2 → T1). By considering multiple families of materials, we show that the T2 → S1 ISC can be enhanced in a number of ways: the substitution of electron-donating (ED) and electron-withdrawing (EW) groups at appropriate positions; the substitution of bulky groups that distort the molecular geometry; and the substitution of heavy atoms that enhance the spin-orbit coupling (SOC). In the first two cases, the enhancements are consistent with El-Sayed's rule in that rapid T2 → S1 ISC requires significant differences in the characters of the S1 and the T2 wavefunctions. Together, these effects enable a wide tunability of T2 → S1 ISC rates over at least 5 orders of magnitude. Meanwhile, the T2 → T1 IC is inhibited in these systems due to the large T2 - T1 energy gap >0.5 eV, which entails a high energy barrier to the T2 → T1 IC and the prediction of a slow rate regardless of the substituents or the presence of heavy atoms. In this way, tuning the T2 → S1 ISC appears to provide an effective strategy to achieve systematic improvement of TTU materials.

First-principles calculations of anharmonic and deuteration effects on the photophysical properties of polyacenes and porphyrinoids

A new method for calculating internal conversion rate constants (k[combining low line]IC), including anharmonic effects and using the Lagrangian multiplier technique, is proposed. The deuteration effect on k[combining low line]IC is investigated for naphthalene, anthracene, free-base porphyrin (H2P) and tetraphenylporphyrin (H2TPP). The results show that anharmonic effects are important when calculating k[combining low line]IC for transitions between electronic states that are energetically separated (ΔE) by more than 20 000-25 000 cm-1. Anharmonic effects are also important when ΔE < 20 000-25 000 cm-1 and when the accepting modes are X-H stretching vibrations with a frequency larger than 2000 cm-1. The calculations show that there is mixing between the S1 and S2 states of naphthalene induced by non-adiabatic interactions. The non-adiabatic interaction matrix element between the S1 and S2 states is 250 cm-1 and 50 cm-1 for the normal and fully deuterated naphthalene structure and this difference significantly affects the estimated fluorescence quantum yield. Besides aromatic hydrocarbons H2P and H2TPP, the k[combining low line]IC rate constant is also calculated for pyrometene (PM567) and tetraoxa[8]circulene (4B) with a detailed analysis of the effect of the vibrational anharmonicity.

Dynamics and spectroscopy of molecular ensembles in a lossy microcavity

The radiative and nonradiative relaxation dynamics of an ensemble of molecules in a microcavity are investigated with emphasis on the impact of the cavity lifetime on reactive and spectroscopic properties. Extending a previous study [I. S. Ulusoy et al., J. Phys. Chem. A 123, 8832-8844 (2019)], it is shown that the dynamics of the ensemble and of single molecules are influenced by the presence of a cavity resonance as long as the polariton splitting can be resolved spectroscopically, which critically depends on the lifetime of the system. Our simulations illustrate how the branching between nonradiative intersystem crossing and radiative decay through the cavity can be tuned by selecting specific cavity photon energies resonant at specific molecular geometries. In the case of cavity-photon energies that are not resonant at the Franck-Condon geometry of the molecules, it is demonstrated numerically and analytically that collective effects are limited to a handful of molecules in the ensemble.

Franck–Condon factors—Computational approaches and recent developments

Algorithms and methodologies for the calculation of Franck–Condon integrals are reviewed. Starting from the standard approach based on the Cartesian representation of the normal modes and the use of Duschinsky's transformation and recursion formulas, methods for treating large displacements of the equilibrium positions and anharmonic and non-Condon effects are considered. Finally, focusing attention on problems arising from the use of recurrence relations, some of the proposed solutions are critically reviewed along with new alternative approaches.

Quantum chemical insights into the aggregation induced emission phenomena: a QM/MM study for pyrazine derivatives

There have been intensive studies on the newly discovered phenomena called aggregation induced emission (AIE), in contrast to the conventional aggregation quenching. Through combined quantum mechanics and molecular mechanics computations, we have investigated the aggregation effects on the excited state decays, both via radiative and nonradiative routes, for pyrazine derivatives 2,3-dicyano-5,6-diphenylpyrazine (DCDPP) and 2,3-dicyanopyrazino phenanthrene (DCPP) in condensed phase. We show that for DCDPP there appear AIE for all the temperature, because the phenyl ring torsional motions in gas phase can efficiently dissipate the electronic excited state energy, and get hindered in aggregate; while for its "locked"-phenyl counterpart, DCPP, theoretical calculation can only give the normal aggregation quenching. These first-principles based findings are consistent with recent experiment. The primary origin of the exotic AIE phenomena is due to the nonradiative decay effects. This is the first time that AIE is understood based on theoretical chemistry calculations for aggregates, which helps to resolve the present disputes over the mechanism.

Theory of excited state decays and optical spectra: application to polyatomic molecules

General formalism of absorption and emission spectra, and of radiative and nonradiative decay rates are derived using a thermal vibration correlation function formalism for the transition between two adiabatic electronic states in polyatomic molecules. Displacements, distortions, and Duschinsky rotation of potential energy surfaces are included within the framework of a multidimensional harmonic oscillator model. The Herzberg-Teller (HT) effect is also taken into account. This formalism gives a reliable description of the Q(x) spectral band of free-base porphyrin with weakly electric dipole-allowed transitions. For the strongly dipole-allowed transitions, e.g., S(1) --> S(0) and S(0) --> S(1) of linear polyacenes, anthracene, tetracene, and pentacene, the HT effect is found to enhance the radiative decay rates by approximately 10% compared to those without the HT effect. For nonradiative transition processes, a general formalism is presented to extend the application scope of the internal conversion theory by going beyond the promoting-mode approximation. Numerical calculations for the nonradiative S(1) --> S(0) decay rate of azulene well explain the origin of the violation of Kasha's rule. When coupled with first-principles density functional theory (DFT) calculations, the present approach appears to be an effective tool to obtain a quantitative description and detailed understanding of spectra and photophysical processes in polyatomic molecules.