# 12.11.1 Theory

(May 16, 2021)

So far we have only covered EDA methods for intermolecular interactions between molecules in their ground states. Since electronic excited states are associated with less strongly bound electrons, modified electrostatic multipole moments (due to electron transition), and often larger polarizabilities, effects imposed by other molecules can be even larger as well as less chemically intuitive than those on ground states. Furthermore, there exist systems that are weakly bound in the ground state but much more strongly bound in the electronic excited state (e.g. He${}_{2}$ vs. He${}_{2}^{\ast}$). Therefore, it is very desirable to develop an interpretation tool that can be utilized to study these important phenomena that are related to intermolecular interactions involving excited-state molecules.

Ge et al. recently extended the ALMO-EDA to treat exciplexes (where the excitation can be assigned to a single molecule within a complex) and excimers (where multiple fragments contribute to the excitation) computed at the CIS or TDDFT/TDA level of theory. Here we briefly overview the decomposition schemes. In the EDA for exciplexes, one first defines the interaction energy in the excited state ($\Delta E_{\mathrm{INT}}^{\ast}$) as

 $\Delta E_{\mathrm{INT}}^{\ast}=E^{\ast}-E_{\mathrm{frag}}^{\ast}$ (12.29)

where $E^{\ast}=E+\omega$ is the energy of the excited supersystem, and $E_{\mathrm{frag}}^{\ast}$ can be expressed as the sum of ground-state fragment energies and the excitation energy of one of the fragments (without losing generality, this excited fragment is denoted as fragment “1”):

 $E_{\mathrm{frag}}^{\ast}=\sum_{F}E_{F}+\omega_{1}$ (12.30)

Therefore, we can rewrite the excited-state interaction as

 $\Delta E_{\mathrm{INT}}^{\ast}=\Delta E_{\mathrm{INT}}+\Delta\omega_{\mathrm{% INT}}$ (12.31)

which contains contributions from the ground-state interaction energy ($\Delta E=E-\sum_{F}E_{F}$) and the excitation energy ($\Delta\omega_{\mathrm{INT}}=\omega-\omega_{1}$). Then, as in the first-generation ALMO-EDA for ground states , the excited-state interaction energy can be separated into contributions from frozen interaction (FRZ), polarization (POL), and charge transfer (CT):

 $\Delta E_{\mathrm{INT}}^{\ast}=\Delta E_{\mathrm{FRZ}}^{\ast}+\Delta E_{% \mathrm{POL}}^{\ast}+\Delta E_{\mathrm{CT}}^{\ast}$ (12.32)

Each term on the RHS of eq. 12.32 can be written in a similar form as eq. 12.31:

 \displaystyle\begin{aligned} \displaystyle\Delta E_{\mathrm{FRZ}}^{\ast}&% \displaystyle=\Delta E_{\mathrm{FRZ}}+\omega_{\mathrm{FRZ}}-\omega_{1}\\ &\displaystyle=\Delta E_{\mathrm{FRZ}}+\Delta\omega_{\mathrm{FRZ}}\\ \displaystyle\Delta E_{\mathrm{POL}}^{\ast}&\displaystyle=\Delta E_{\mathrm{% POL}}+\omega_{\mathrm{POL}}-\omega_{\mathrm{FRZ}}\\ &\displaystyle=\Delta E_{\mathrm{POL}}+\Delta\omega_{\mathrm{POL}}\\ \displaystyle\Delta E_{\mathrm{CT}}^{\ast}&\displaystyle=\Delta E_{\mathrm{CT}% }+\omega-\omega_{\mathrm{POL}}\\ &\displaystyle=\Delta E_{\mathrm{CT}}+\Delta\omega_{\mathrm{CT}}\end{aligned} (12.33)

$\Delta E_{\mathrm{FRZ}}$, $\Delta E_{\mathrm{POL}}$, and $\Delta E_{\mathrm{CT}}$ can be obtained by performing a ground-state ALMO-EDA for the supersystem. To compute $\Delta\omega_{\mathrm{FRZ}}$, $\Delta\omega_{\mathrm{POL}}$, and $\Delta\omega_{\mathrm{CT}}$, one needs to define $\omega_{\mathrm{FRZ}}$ and $\omega_{\mathrm{POL}}$, i.e., excitation energies associated with the frozen and polarized supersystem, respectively. The frozen intermediate state can be viewed as one excited fragment embedded in the environment formed by other ground-state fragments, whose effects on the excited fragment are only through the supersystem Fock matrix. The definition of the polarized intermediate state utilizes the ALMO-CIS model (see Sec. 12.19), where both MOs and excitation amplitudes are fragment-localized. We also note that the frozen contribution to the excited-state interaction energy, $\Delta E_{\mathrm{FRZ}}^{\ast}$, can be further partitioned into a classical electrostatics term (Coulomb interactions between isolated fragment charge distributions) and a non-electrostatic term [mostly Pauli repulsion if a non-dispersion-corrected model (e.g. CIS) is used]:

 $\Delta E_{\mathrm{FRZ}}^{\ast}=\Delta E_{\textrm{CLS-ELEC}}^{\ast}+\Delta E_{% \textrm{NON-ELEC}}^{\ast}$ (12.34)

Modifications are needed in order to extend this method to excimers, where different fragments are of degenerate or near-degenerate excited states. In such cases, we choose $M$ reference fragment states as the initial basis. Denote the $s^{\mathrm{th}}$ excited state on fragment $I$ as the $\kappa^{\mathrm{th}}$ reference state ($\kappa=1,2,\dots,M$). Similar to eqn. 12.30, we have

 $E_{\mathrm{frag}}^{\kappa}=\sum_{F}E_{F}+\omega_{I}^{s}$ (12.35)

The corresponding frozen excited-state wavefunction is then constructed by embedding this excited fragment into the environment formed by other fragments in their ground states:

 $\left|\Phi_{\mathrm{FRZ}}^{\kappa}\right\rangle=\left|\Psi_{1}\Psi_{2}\cdots% \Psi_{I}^{s}\cdots\Psi_{N}\right\rangle$ (12.36)

and the excited-state frozen interaction energy

 $\Delta E_{\mathrm{FRZ}}^{\kappa}=\Delta E_{\mathrm{FRZ}}+\Delta\omega_{\mathrm% {FRZ}}^{\kappa}=\Delta E_{\mathrm{FRZ}}+(\omega_{\mathrm{FRZ}}^{\kappa}-\omega% _{I}^{s})$ (12.37)

With $M$ degenerate or near-degenerate frozen excited states, a new intermediate state is then introduced to capture the pure excitonic-splitting (EXSP) effect in the formation of excimers, which can be expressed as a linear combination of the frozen states:

 $\left|\Phi^{\kappa}_{\mathrm{EXSP}}\right\rangle=\sum_{\kappa^{\prime}}^{M}c^{% \kappa\kappa^{\prime}}\left|\Psi_{1}\Psi_{2}\cdots\Psi_{I}^{s}\cdots\Psi_{N}\right\rangle$ (12.38)

The associated excitation energy $\omega^{\kappa}_{\mathrm{EXSP}}$ and the corresponding linear combination coefficients can be obtained by solving a secular equation in the basis of frozen states. As excitonic splitting is purely an excited-state phenomenon, we have

 $\Delta E^{\kappa}_{\mathrm{EXSP}}=\Delta\omega^{\kappa}_{\mathrm{EXSP}}=\omega% ^{\kappa}_{\mathrm{EXSP}}-\omega^{\kappa}_{\mathrm{FRZ}}$ (12.39)

Subsequently, polarization and charge transfer are handled in a similar way as in the excimer case:

 \displaystyle\begin{aligned} \displaystyle\Delta E^{\kappa}_{\mathrm{POL}}&% \displaystyle=\Delta E_{\mathrm{POL}}+\Delta\omega^{\kappa}_{\mathrm{POL}}=% \Delta E_{\mathrm{POL}}+(\omega^{\kappa}_{\mathrm{POL}}-\omega^{\kappa}_{% \mathrm{EXSP}})\\ \displaystyle\Delta E^{\kappa}_{\mathrm{CT}}&\displaystyle=\Delta E_{\mathrm{% CT}}+\Delta\omega^{\kappa}_{\mathrm{CT}}=\Delta E_{\mathrm{CT}}+(\omega^{% \kappa}-\omega^{\kappa}_{\mathrm{POL}})\end{aligned} (12.40)

One more complication compared to the EDA scheme for exciplexes is that since multiple ($M$) states are considered, extra caution needs be paid to the state-ordering at different stages (EXSP, POL and CT). In order to locate the states of interest (which can be most unambiguously identified at the EXSP stage) correctly during the entire EDA procedure, a state-tracking algorithm based on a maximum-overlap criterion is employed. The reader is referred to Ref. 345 for more details.