# Interplay of creation, propagation, and relaxation of an excitation in a dimer

## A

Definitions of matrices in Eq. (20)

Appendix contains the definitions of the matrices $`𝒥_1`$, $`𝒥_2`$, $`𝒢_1(t)`$, $`𝒢_2(t)`$, $`_1(t)`$, $`_2(t)`$, $`_3(t)`$, and $`_4(t)`$ in Eq. (20). The influence of their elements on the excitation dynamics is also discussed.

The matrices $`𝒥_1`$ and $`𝒥_2`$ describe the dynamics of the free exciton system:

$`𝒥_1={\displaystyle \frac{1}{\mathrm{}}}\left[\begin{array}{ccccc}0& 0& 0& 2J& 0\\ 0& 0& 0& 2J& 0\\ 0& 0& 0& 2\epsilon & 0\\ J& J& 2\epsilon & 0& 0\\ 0& 0& 0& 0& 0\end{array}\right],𝒥_2={\displaystyle \frac{1}{\mathrm{}}}\left[\begin{array}{cccc}0& \epsilon & 0& J\\ \epsilon & 0& J& 0\\ 0& J& 0& \epsilon \\ J& 0& \epsilon & 0\end{array}\right].`$ (A10)   
The matrices $`𝒢_1(t)`$ and $`𝒢_2(t)`$ stemm from the exciton-phonon interaction:

$`𝒢_1(t)=\left[\begin{array}{ccccc}0& 0& 0& 0& 0\\ 0& 0& 0& 0& 0\\ A& C& E& F& 0\\ B& D& F& E& 0\\ 0& 0& 0& 0& 0\end{array}\right],𝒢_2(t)=\left[\begin{array}{cccc}A_1& B_1& C_1& D_1\\ B_1& A_1& D_1& C_1\\ C_2& D_2& A_2& B_2\\ D_2& C_2& B_2& A_2\end{array}\right].`$ (A20)   
The time-dependent coefficients $`A(t),\mathrm{},D_2(t)`$ are given as follows

$`A(t)`$ $`=`$ $`{\displaystyle \frac{J\epsilon }{\mathrm{\Delta }^2}}G^2\overline{g}_{2,1}+{\displaystyle \frac{J}{2\mathrm{\Delta }}}G^2\overline{g}_{3,2},`$ (A21)   
$`B(t)`$ $`=`$ $`{\displaystyle \frac{J\epsilon }{\mathrm{\Delta }^2}}G^2\overline{g}_{2,2}{\displaystyle \frac{J}{2\mathrm{\Delta }}}G^2\overline{g}_{3,1},`$ (A22)   
$`C(t)`$ $`=`$ $`{\displaystyle \frac{J\epsilon }{\mathrm{\Delta }^2}}G^2\overline{g}_{2,1}+{\displaystyle \frac{J}{2\mathrm{\Delta }}}G^2\overline{g}_{3,2},`$ (A23)   
$`D(t)`$ $`=`$ $`{\displaystyle \frac{J\epsilon }{\mathrm{\Delta }^2}}G^2\overline{g}_{2,2}{\displaystyle \frac{J}{2\mathrm{\Delta }}}G^2\overline{g}_{3,1},`$ (A24)   
$`E(t)`$ $`=`$ $`G^2\overline{g}_{1,1}2{\displaystyle \frac{J^2}{\mathrm{\Delta }^2}}G^2\overline{g}_{2,1},`$ (A25)   
$`F(t)`$ $`=`$ $`G^2\overline{g}_{1,3},`$ (A26)   
$`A_1(t)`$ $`=`$ $`G^2\overline{g}_{1,4}+{\displaystyle \frac{J^2}{\mathrm{\Delta }^2}}G^2\left[\overline{g}_{2,4}+\overline{g}_{2,5}+\overline{g}_{2,10}\right],`$ (A27)   
$`B_1(t)`$ $`=`$ $`G^2\overline{g}_{1,8}+{\displaystyle \frac{J^2}{\mathrm{\Delta }^2}}G^2\left[\overline{g}_{2,6}\overline{g}_{2,8}+\overline{g}_{2,9}\right],`$ (A28)   
$`A_2(t)`$ $`=`$ $`G^2\overline{g}_{1,7}+{\displaystyle \frac{J^2}{\mathrm{\Delta }^2}}G^2\left[\overline{g}_{2,5}\overline{g}_{2,7}\overline{g}_{2,10}\right],`$ (A29)   
$`B_2(t)`$ $`=`$ $`G^2\overline{g}_{1,11}+{\displaystyle \frac{J^2}{\mathrm{\Delta }^2}}G^2\left[\overline{g}_{2,6}+\overline{g}_{2,9}\overline{g}_{2,11}\right],`$ (A30)   
$`C_1(t)`$ $`=`$ $`{\displaystyle \frac{J\epsilon }{\mathrm{\Delta }^2}}G^2\left[\overline{g}_{2,4}\overline{g}_{2,5}\overline{g}_{2,10}\right]+{\displaystyle \frac{J}{2\mathrm{\Delta }}}G^2\left[\overline{g}_{3,6}+\overline{g}_{3,8}\overline{g}_{3,9}\right],`$ (A31)   
$`D_1(t)`$ $`=`$ $`{\displaystyle \frac{J\epsilon }{\mathrm{\Delta }^2}}G^2\left[\overline{g}_{2,6}+\overline{g}_{2,8}\overline{g}_{2,9}\right]+{\displaystyle \frac{J}{2\mathrm{\Delta }}}G^2\left[\overline{g}_{3,4}+\overline{g}_{3,5}+\overline{g}_{3,10}\right],`$ (A32)   
$`C_2(t)`$ $`=`$ $`{\displaystyle \frac{J\epsilon }{\mathrm{\Delta }^2}}G^2\left[\overline{g}_{2,5}\overline{g}_{2,7}\overline{g}_{2,10}\right]+{\displaystyle \frac{J}{2\mathrm{\Delta }}}G^2\left[\overline{g}_{3,6}\overline{g}_{3,9}+\overline{g}_{3,11}\right],`$ (A33)   
$`D_2(t)`$ $`=`$ $`{\displaystyle \frac{J\epsilon }{\mathrm{\Delta }^2}}G^2\left[\overline{g}_{2,6}+\overline{g}_{2,9}\overline{g}_{2,11}\right]+{\displaystyle \frac{J}{2\mathrm{\Delta }}}G^2\left[\overline{g}_{3,5}\overline{g}_{3,7}\overline{g}_{3,10}\right].`$ (A34)   
The functions

$`\overline{g}_{1,j}(t)`$ $`=`$ $`{\displaystyle _0^{tt_0}}𝑑\tau g_j(\tau ),`$ (A35)   
$`\overline{g}_{2,j}(t,\mathrm{\Delta }^{})`$ $`=`$ $`{\displaystyle _0^{tt_0}}𝑑\tau g_j(\tau )\mathrm{sin}^2(\mathrm{\Delta }^{}\tau ),`$ (A36)   
$`\overline{g}_{3,j}(t,\mathrm{\Delta }^{})`$ $`=`$ $`{\displaystyle _0^{tt_0}}𝑑\tau g_j(\tau )\mathrm{sin}(2\mathrm{\Delta }^{}\tau ),j=1,\mathrm{},11`$ (A37)   
describe the response of the exciton subsystem to the phonon one. The phonon subsystem is characterized by the functions In numerical calculations, we assume that $`\mathrm{}\mathrm{\Omega }_kG_k^i=G_i`$ in one half of the $`k`$–space and $`\mathrm{}\mathrm{\Omega }_kG_k^i=G_i^{}`$ in the remaining half of the $`k`$–space ($`i=1,2`$) (for details, see, Ref. ). This assumption is in agreement with hermiticity of $`\widehat{H}_{\mathrm{e}\mathrm{ph}}`$. Further, the mean numbers of phonons $`n_\mathrm{B}(\mathrm{}\mathrm{\Omega }_k)`$ are assumed to be $`k`$–independent ($`n_\mathrm{B}(\mathrm{}\mathrm{\Omega }_{k_0})=n_\mathrm{B}`$). The remaining summations $`\frac{1}{N}_k\mathrm{sin}(\mathrm{\Omega }_k\tau )`$ and $`\frac{1}{N}_k\mathrm{cos}(\mathrm{\Omega }_k\tau )`$ in Eqs. (A6) are replaced by the expressions $`\mathrm{sin}(\mathrm{\Omega }_{ph}\tau )\mathrm{exp}(\gamma _{ph}\tau )`$ and $`\mathrm{cos}(\mathrm{\Omega }_{ph}\tau )\mathrm{exp}(\gamma _{ph}\tau )`$, respectively. The frequency $`\mathrm{\Omega }_{ph}`$ then characterizes a mean phonon oscillation frequency and $`\gamma _{ph}`$ describes damping originating in dephasing.

$`g_1(\tau )`$ $`=`$ $`{\displaystyle \frac{1}{G^2N}}{\displaystyle \underset{k}{}}\mathrm{\Omega }_k^2|G_k^1G_k^2|^2\left[2n_\mathrm{B}(\mathrm{}\mathrm{\Omega }_k)+1\right]\mathrm{cos}(\mathrm{\Omega }_k\tau ),`$ (A38)   
$`g_2(\tau )`$ $`=`$ $`{\displaystyle \frac{1}{G^2N}}{\displaystyle \underset{k}{}}\mathrm{\Omega }_k^2|G_k^1G_k^2|^2\mathrm{sin}(\mathrm{\Omega }_k\tau ),`$ (A39)   
$`g_3(\tau )`$ $`=`$ $`{\displaystyle \frac{1}{G^2N}}{\displaystyle \underset{k}{}}\mathrm{\Omega }_k^2(G_k^1G_k^2)(G_k^1+G_k^2)\mathrm{sin}(\mathrm{\Omega }_k\tau ),`$ (A40)   
$`g_4(\tau )`$ $`=`$ $`{\displaystyle \frac{1}{G^2N}}{\displaystyle \underset{k}{}}\mathrm{\Omega }_k^2|G_k^1|^2\left[2n_\mathrm{B}(\mathrm{}\mathrm{\Omega }_k)+1\right]\mathrm{cos}(\mathrm{\Omega }_k\tau ),`$ (A41)   
$`g_5(\tau )`$ $`=`$ $`{\displaystyle \frac{1}{G^2N}}{\displaystyle \underset{k}{}}\mathrm{\Omega }_k^2\mathrm{Re}[G_k^1G_k^2]\left[2n_\mathrm{B}(\mathrm{}\mathrm{\Omega }_k)+1\right]\mathrm{cos}(\mathrm{\Omega }_k\tau ),`$ (A42)   
$`g_6(\tau )`$ $`=`$ $`{\displaystyle \frac{1}{G^2N}}{\displaystyle \underset{k}{}}\mathrm{\Omega }_k^2\mathrm{Im}[G_k^1G_k^2]\left[2n_\mathrm{B}(\mathrm{}\mathrm{\Omega }_k)+1\right]\mathrm{cos}(\mathrm{\Omega }_k\tau ),`$ (A43)   
$`g_7(\tau )`$ $`=`$ $`{\displaystyle \frac{1}{G^2N}}{\displaystyle \underset{k}{}}\mathrm{\Omega }_k^2|G_k^2|^2\left[2n_\mathrm{B}(\mathrm{}\mathrm{\Omega }_k)+1\right]\mathrm{cos}(\mathrm{\Omega }_k\tau ),`$ (A44)   
$`g_8(\tau )`$ $`=`$ $`{\displaystyle \frac{1}{G^2N}}{\displaystyle \underset{k}{}}\mathrm{\Omega }_k^2|G_k^1|^2\mathrm{sin}(\mathrm{\Omega }_k\tau ),`$ (A45)   
$`g_9(\tau )`$ $`=`$ $`{\displaystyle \frac{1}{G^2N}}{\displaystyle \underset{k}{}}\mathrm{\Omega }_k^2\mathrm{Re}[G_k^1G_k^2]\mathrm{sin}(\mathrm{\Omega }_k\tau ),`$ (A46)   
$`g_{10}(\tau )`$ $`=`$ $`{\displaystyle \frac{1}{G^2N}}{\displaystyle \underset{k}{}}\mathrm{\Omega }_k^2\mathrm{Im}[G_k^1G_k^2]\mathrm{sin}(\mathrm{\Omega }_k\tau ),`$ (A47)   
$`g_{11}(\tau )`$ $`=`$ $`{\displaystyle \frac{1}{G^2N}}{\displaystyle \underset{k}{}}\mathrm{\Omega }_k^2|G_k^2|^2\mathrm{sin}(\mathrm{\Omega }_k\tau ).`$ (A48)   
The constant $`G`$ having the meaning of the mean exciton–phonon interaction constant has been introduced into Eqs. (A3) and (A4) as well as into the definitions given in Eq. (A6) in order to get $`g_1(\tau ),\mathrm{},g_{11}(\tau )`$ dependent only on the dispersion of coupling constants. The symbols $`\mathrm{Re}`$ and $`\mathrm{Im}`$ denote real and imaginary parts. The new symbol $`\mathrm{\Delta }^{}=\mathrm{\Delta }/\mathrm{}`$ has been introduced here.

The time-dependent coefficients $`B(t)`$ and $`D(t)`$ given in Eq. (A3) renormalize the transfer integral $`J`$. The coefficients $`A(t)`$ and $`C(t)`$ are important for relaxation to equilibrium state (see, e.g., Refs. ).

The matrices $`_1(t)`$, $`_2(t)`$, $`_3(t)`$, and $`_4(t)`$ originate in the interaction with an optical field:

$`_1(t)`$ $`=`$ $`\left[\begin{array}{ccccc}2\overline{M}_1& 0& 2\overline{O}_1& 2\overline{O}_2& 2M_1\\ 0& 2\overline{N}_1& 2\overline{P}_1& 2\overline{P}_2& 2N_1\\ \overline{P}_1& \overline{O}_1& \overline{M}_1+\overline{N}_1& \overline{M}_2+\overline{N}_2& O_1P_1\\ \overline{P}_2& \overline{O}_2& \overline{M}_2\overline{N}_2& \overline{M}_1+\overline{N}_1& P_2O_2\\ 2\overline{M}_1& 2\overline{N}_1& 2\overline{P}_12\overline{O}_1& 2\overline{P}_22\overline{O}_2& 2M_1+2N_1\end{array}\right],`$ (A54)   
$`_2(t)`$ $`=`$ $`\left[\begin{array}{cccc}2K_1& 2K_2& 0& 0\\ 0& 0& 2L_1& 2L_2\\ L_1& L_2& K_1& K_2\\ L_2& L_1& K_2& K_1\\ 2K_1& 2K_2& 2L_1& 2L_2\end{array}\right],`$ (A60)   
$`_3(t)`$ $`=`$ $`\left[\begin{array}{ccccc}K_1& 0& L_1& L_2& K_1\\ K_2& 0& L_2& L_1& K_2\\ 0& L_1& K_1& K_2& L_1\\ 0& L_2& K_2& K_1& L_2\end{array}\right],`$ (A65)   
$`_4(t)`$ $`=`$ $`\left[\begin{array}{cccc}2M_1+N_12\stackrel{~}{M}_1& 2M_2+N_22\stackrel{~}{M}_2& O_1\stackrel{~}{O}_1\stackrel{~}{P}_1& O_2\stackrel{~}{O}_2\stackrel{~}{P}_2\\ 2M_2N_22\stackrel{~}{M}_2& 2M_1+N_1+2\stackrel{~}{M}_1& O_2\stackrel{~}{O}_2\stackrel{~}{P}_2& O_1+\stackrel{~}{O}_1+\stackrel{~}{P}_1\\ P_1\stackrel{~}{O}_1\stackrel{~}{P}_1& P_2\stackrel{~}{O}_2\stackrel{~}{P}_2& 2N_1+M_12\stackrel{~}{N}_1& 2N_2+M_22\stackrel{~}{N}_2\\ P_2\stackrel{~}{O}_2\stackrel{~}{P}_2& P_1+\stackrel{~}{O}_1+\stackrel{~}{P}_1& 2N_2M_22\stackrel{~}{N}_2& 2N_1+M_1+2\stackrel{~}{N}_1\end{array}\right].`$ (A70)   
The coefficients $`K_1(t)`$, $`K_2(t)`$, $`L_1(t)`$, and $`L_2(t)`$ describing the influence of the coherent part of an optical field have the form (the constants $`\stackrel{~}{F}_{K_0}^1`$ and $`\stackrel{~}{F}_{K_0}^2`$ are assumed to be real):

$`K_1(t)`$ $`=`$ $`\omega _{K_0}\stackrel{~}{F}_{K_0}^1\mathrm{Im}\left[\stackrel{~}{𝒜}(t)\mathrm{exp}(i\delta ^{}t)\right],`$ (A72)   
$`K_2(t)`$ $`=`$ $`\omega _{K_0}\stackrel{~}{F}_{K_0}^1\mathrm{Re}\left[\stackrel{~}{𝒜}(t)\mathrm{exp}(i\delta ^{}t)\right],`$ (A73)   
$`L_1(t)`$ $`=`$ $`\omega _{K_0}\stackrel{~}{F}_{K_0}^2\mathrm{Im}\left[\stackrel{~}{𝒜}(t)\mathrm{exp}(i\delta ^{}t)\right],`$ (A74)   
$`L_2(t)`$ $`=`$ $`\omega _{K_0}\stackrel{~}{F}_{K_0}^2\mathrm{Re}\left[\stackrel{~}{𝒜}(t)\mathrm{exp}(i\delta ^{}t)\right].`$ (A75)   
The symbol $`\delta ^{}`$ denotes the frequency mismatch ($`\delta ^{}=(E+\epsilon )/\mathrm{}\omega _{K_0}`$) and the envelope $`\stackrel{~}{𝒜}(t)`$ of the field is defined as follows,

$$\stackrel{~}{𝒜}(t)=𝒜(t)\mathrm{exp}(i\omega _{K_0}t).$$
(A76)   
The coefficients $`M_1(t),\mathrm{},P_2(t)`$ reflect statistical properties of the optical field (noise) and can be expressed in the form

$`M_1(t)`$ $`=`$ $`(\stackrel{~}{F}_{K_0}^1)^2i_1(\stackrel{~}{F}_{K_0}^1)^2{\displaystyle \frac{\epsilon }{\mathrm{\Delta }}}i_4\stackrel{~}{F}_{K_0}^1\stackrel{~}{F}_{K_0}^2{\displaystyle \frac{J}{\mathrm{\Delta }}}i_4,`$ (A77)   
$`M_2(t)`$ $`=`$ $`(\stackrel{~}{F}_{K_0}^1)^2i_3+(\stackrel{~}{F}_{K_0}^1)^2{\displaystyle \frac{\epsilon }{\mathrm{\Delta }}}i_2+\stackrel{~}{F}_{K_0}^1\stackrel{~}{F}_{K_0}^2{\displaystyle \frac{J}{\mathrm{\Delta }}}i_2,`$ (A78)   
$`N_1(t)`$ $`=`$ $`(\stackrel{~}{F}_{K_0}^2)^2i_1+(\stackrel{~}{F}_{K_0}^2)^2{\displaystyle \frac{\epsilon }{\mathrm{\Delta }}}i_4\stackrel{~}{F}_{K_0}^1\stackrel{~}{F}_{K_0}^2{\displaystyle \frac{J}{\mathrm{\Delta }}}i_4,`$ (A79)   
$`N_2(t)`$ $`=`$ $`(\stackrel{~}{F}_{K_0}^2)^2i_3(\stackrel{~}{F}_{K_0}^2)^2{\displaystyle \frac{\epsilon }{\mathrm{\Delta }}}i_2+\stackrel{~}{F}_{K_0}^1\stackrel{~}{F}_{K_0}^2{\displaystyle \frac{J}{\mathrm{\Delta }}}i_2,`$ (A80)   
$`O_1(t)`$ $`=`$ $`(\stackrel{~}{F}_{K_0}^1)^2{\displaystyle \frac{J}{\mathrm{\Delta }}}i_4+\stackrel{~}{F}_{K_0}^1\stackrel{~}{F}_{K_0}^2i_1+\stackrel{~}{F}_{K_0}^1\stackrel{~}{F}_{K_0}^2{\displaystyle \frac{\epsilon }{\mathrm{\Delta }}}i_4,`$ (A81)   
$`O_2(t)`$ $`=`$ $`(\stackrel{~}{F}_{K_0}^1)^2{\displaystyle \frac{J}{\mathrm{\Delta }}}i_2+\stackrel{~}{F}_{K_0}^1\stackrel{~}{F}_{K_0}^2i_3\stackrel{~}{F}_{K_0}^1\stackrel{~}{F}_{K_0}^2{\displaystyle \frac{\epsilon }{\mathrm{\Delta }}}i_2,`$ (A82)   
$`P_1(t)`$ $`=`$ $`(\stackrel{~}{F}_{K_0}^2)^2{\displaystyle \frac{J}{\mathrm{\Delta }}}i_4+\stackrel{~}{F}_{K_0}^1\stackrel{~}{F}_{K_0}^2i_1\stackrel{~}{F}_{K_0}^1\stackrel{~}{F}_{K_0}^2{\displaystyle \frac{\epsilon }{\mathrm{\Delta }}}i_4,`$ (A83)   
$`P_2(t)`$ $`=`$ $`(\stackrel{~}{F}_{K_0}^2)^2{\displaystyle \frac{J}{\mathrm{\Delta }}}i_2+\stackrel{~}{F}_{K_0}^1\stackrel{~}{F}_{K_0}^2i_3+\stackrel{~}{F}_{K_0}^1\stackrel{~}{F}_{K_0}^2{\displaystyle \frac{\epsilon }{\mathrm{\Delta }}}i_2.`$ (A84)   
The functions $`i_1(t,\mathrm{\Delta }^{},\delta ^{}),\mathrm{},i_4(t,\mathrm{\Delta }^{},\delta ^{})`$ characterize the response of the exciton subsystem to the photon field:

$`i_1(t,\mathrm{\Delta }^{},\delta ^{})`$ $`=`$ $`\omega _{K_0}^2{\displaystyle _{t_0}^t}𝑑\tau \mathrm{cos}\left[\mathrm{\Delta }^{}(t\tau )\right]\mathrm{Re}\left[\delta \stackrel{~}{N}(t,\tau )\mathrm{exp}\left[i\delta ^{}(t\tau )\right]\right],`$ (A85)   
$`i_2(t,\mathrm{\Delta }^{},\delta ^{})`$ $`=`$ $`\omega _{K_0}^2{\displaystyle _{t_0}^t}𝑑\tau \mathrm{sin}\left[\mathrm{\Delta }^{}(t\tau )\right]\mathrm{Re}\left[\delta \stackrel{~}{N}(t,\tau )\mathrm{exp}\left[i\delta ^{}(t\tau )\right]\right],`$ (A86)   
$`i_3(t,\mathrm{\Delta }^{},\delta ^{})`$ $`=`$ $`\omega _{K_0}^2{\displaystyle _{t_0}^t}𝑑\tau \mathrm{cos}\left[\mathrm{\Delta }^{}(t\tau )\right]\mathrm{Im}\left[\delta \stackrel{~}{N}(t,\tau )\mathrm{exp}\left[i\delta ^{}(t\tau )\right]\right],`$ (A87)   
$`i_4(t,\mathrm{\Delta }^{},\delta ^{})`$ $`=`$ $`\omega _{K_0}^2{\displaystyle _{t_0}^t}𝑑\tau \mathrm{sin}\left[\mathrm{\Delta }^{}(t\tau )\right]\mathrm{Im}\left[\delta \stackrel{~}{N}(t,\tau )\mathrm{exp}\left[i\delta ^{}(t\tau )\right]\right].`$ (A88)   
The photon field correlation function $`\delta \stackrel{~}{N}(t,\tau )`$ is of the form:

$$\delta \stackrel{~}{N}(t,\tau )=\delta N(t,\tau )\mathrm{exp}\left[i\omega _{K_0}(t\tau )\right].$$
(A89)   
The coefficients $`\overline{M}_1(t),\mathrm{},\overline{P}_2(t)`$ are defined similarly as the coefficients $`M_1(t),\mathrm{},P_2(t)`$ in Eqs. (A10) and (A11); only the photon field correlation function

$$\delta \stackrel{~}{N}_v(t,\tau )=\delta N_v(t,\tau )\mathrm{exp}\left[i\omega _{K_0}(t\tau )\right]$$
(A90)   
occurs in Eq. (A11) instead of $`\delta \stackrel{~}{N}(t,\tau )`$. Thus, the coefficients with bars include in addition effects of vacuum fluctuations.

Also the coefficients $`\stackrel{~}{M}_1(t),\mathrm{},\stackrel{~}{P}_2(t)`$ are defined similarly as the coefficients $`M_1(t),\mathrm{},P_2(t)`$ in Eqs. (A10) and (A11): only the expression $`\delta \stackrel{~}{N}(t,\tau )\mathrm{exp}\left[i\delta ^{}(t\tau )\right]`$ in Eq. (A11) must be replaced by the expression $`\delta \stackrel{~}{N}_a(t,\tau )\mathrm{exp}\left[i\delta ^{}(t+\tau )\right]`$, where

$$\delta \stackrel{~}{N}_a(t,\tau )=\delta N_a(t,\tau )\mathrm{exp}\left[i\omega _{K_0}(t+\tau )\right].$$
(A91)   
The coefficients $`\overline{O}_2(t)`$ and $`\overline{P}_2(t)`$ renormalize the transfer integral $`J`$ both at the positions (4,1), (1,4) and (4,2), (2,4) in the matrix $`𝒥_1`$, in contrast to the coefficients originating in the exciton–phonon interaction.

The influence of the coefficient $`\stackrel{~}{M}_1(t)`$ ($`\stackrel{~}{N}_1(t)`$) in the equations for $`\rho _{1r}(t)`$ and $`\rho _{1i}(t)`$ ($`\rho _{2r}(t)`$ and $`\rho _{2i}(t)`$) is remarkable. When, e.g., $`\stackrel{~}{M}_1(t)`$ ($`\stackrel{~}{N}_1(t)`$) is negative, it represents damping of $`\rho _{1r}(t)`$ ($`\rho _{2r}(t)`$), but at the same time amplification of $`\rho _{1i}(t)`$ ($`\rho _{2i}(t)`$). This property is connected with phase relations in the photon field reflected by $`\delta \stackrel{~}{N}_a(t,\tau )`$ given in Eq. (A14).