Full Paper
of the intensity of the two lowest energy bands[13,14] as ob-
served for o2-CP and r2-CP homodimers in Figure 1a,b. Ac-
cordingly, the simple model described above has to be extend-
ed by considering rotational displacements as well as vibra-
tional coupling between the chromophores. Following our ear-
lier work on homoaggregates of PBIs[13] we can be confident
to obtain a conclusive understanding of the absorption bands
of o2-CP and r2-CP homodimers and novel insights into the
hitherto unexplored impact of exciton-vibrational couplings on
the spectra of hetero-PBI cyclophane or-CP.
Towards this goal we have performed time-dependent quan-
tum mechanical simulations. We solve the time-dependent
Schrçdinger equation for a Hamiltonian of Holstein type,[29,30]
which in many cases gives an excellent description of vibronic
spectra[4a,31] and has been applied to self-assembled PBI aggre-
gates before.[13] The absorption spectra are then obtained by
a Fourier transform of the calculated time-dependent autocor-
relation function.[32] A more detailed description of the proce-
dure and the model Hamiltonian can be found in the Support-
ing Information.
lar values for the vibrational frequency and the Huang–Rhys
factor indicate that the substituents attached to the red PBI
monomer do not influence the vibrational motion significantly.
Turning to the homodimers, we fix the angle a between the
transition dipole moments to 28.68 for o2-CP and 24.58 for r2-
CP, respectively (see Figure 2). The absolute values of the tran-
sition dipole moments result, at least in the case of homodim-
ers, only in a global prefactor in the spectrum, and is therefore
set equal to one. We find best agreement for the simulated
spectra with experiment for an exciton coupling energy of J=
613 cmꢀ1 for o2-CP and 500 cmꢀ1 for r2-CP. These values are
comparable to the coupling strength of 524 cmꢀ1 found for
self-assembled homodimers of core-unsubstituted PBI mole-
cules.[13]
The two exciton states illustrated in Figure 3a overlap, be-
cause the coupling strength is small compared to the vibra-
tional frequency. Since the S0!S1‘ transition is only weakly al-
lowed for angles of roughly 258, we attribute the line-shape
mostly to vibrational levels in the S1’’ exciton state. There, the
intensities of the peaks marked with an asterisk in the figure
are strongly dependent on the exciton coupling energy.[31,33]
We note that the Huang–Rhys factors are increased to x=0.88
(orange) and 0.83 (red) to achieve a better agreement with the
experiment. This may indicate that in the dimers a strain acts
on the PBI scaffold, and thus changes the excited-state equilib-
rium geometry of the monomer units. Such possible geometry
deformations are not the focus of this work, and more extend-
ed electronic structure calculations involving excited-state ge-
ometry optimizations would have to be conducted in the
future.
In Figure 4, we compare calculated and measured spectra
for the monochromophoric reference PBIs 4 and 6 (panels a,c)
and the o2-CP and r2-CP homodimers (panels b,d). We use the
monomer spectra to fix the parameters that enter the model
Hamiltonian. These are the frequencies w, Huang–Rhys-factors
x, and excitation energies DE. Best agreement with experiment
is achieved for the sets: w=1411 cmꢀ1, x=0.79, DE=
18954 cmꢀ1 for the orange monomer and w=1331 cmꢀ1, x=
0.78 and DE=17164 cmꢀ1 for the red monomer. The very simi-
We now apply the same model to the or-CP heterodimer.
The calculated spectra for the sum of the monomers and the
heterodimer are shown in Figure 5. In contrast to the homo-
dimers, the relative strength of the transition dipole moments
is important here. Their absolute values will, however, still
result in a global factor. We will therefore only give the
ꢀ
ratio m=m(r)/m(o) of the individual absolute values. The sum of
the monomer spectra is best matched by a model spectrum
for an uncoupled dimer (J=0) when we set the dipole
ꢀ
moment ratio m=0.69. This, however, most probably underes-
timates the value of m(r) because the oscillator strength corre-
sponds to the total area of the spectral band and the spectrum
for the red monomer is broader for reasons discussed above.
All other parameters have been kept at the values determined
from the individual monomer spectra. For the spectrum of the
or-CP heterodimer, we fix the angle a between the transition
dipole moments to a value of 10.68 as indicated in Figure 2c
and the Huang–Rhys factors at the values determined from
the homodimer spectra. From the comparison of calculated
and measured spectra, we arrive at a dipole moment ratio
ꢀ1
ꢀ
of m=0.81 and an exciton coupling energy of J=427 cm . It
Figure 4. Comparison of experimental (d, measured in dichloromethane)
and calculated spectra (c): a) reference PBI 4 (oPBI), b) o2-CP homodimer,
c) reference PBI 6 (rPBI), and d) r2-CP homodimer. The calculated spectra are
normalized to match the measured peak of highest intensity in each case.
The value of the exciton coupling energy J is given for the two dimers. It es-
sentially influences the intensity of the low energy peaks (*), which are
marked in panels b and d. The model does not account for the S0!S2 transi-
tion seen around 22000 cmꢀ1 in panels c and d.
is somewhat smaller than that found for the homodimers, but
still comparable in its numerical value. Again, the coupling
strength has a large effect on the intensity of the 0,0 transition
peak of the oPBI chromophore (marked with an asterisk in Fig-
ure 5b). The ratio of the dipole moments further determines
the intensity of the low-energy band originating mainly from
&
&
Chem. Eur. J. 2016, 22, 1 – 9
6
ꢁ 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim
ÝÝ These are not the final page numbers!