A R T I C L E S
Coropceanu et al.
Table 2. Electronic Coupling Hab (eV) for Systems 1+-6+ As
symmetric modes Q+.6,7,25,26 As a consequence, the effective
electronic coupling depends on the coordinate of the symmetric
vibration(s)7 (see Figure 2b); therefore, using the ∆IP value will
give a lower limit estimate for Hab.
Obtained from Optical and UPS Measurements
optical
∆IP/2
(∆IP/2) + La
1+
2+
3+
4+
5+
6+
1.01b,c
0.79
0.42
0.33
0.35
0.15
0.19
-
-
0.61b,c
The gas-phase photoelectron spectra of 1-6 are shown in
Figure 3; the details of the first ionizations are shown in Figure
4. From these figures and Table 1, it is seen that replacing the
methyl groups of 1 and 2 with phenyl groups in 3 and 5 has
little effect on the first ionizations IP1 but leads to very
significant reductions in ∆IP splitting; on the other hand,
methoxy substitution on the terminal rings when going to 4 and
6 causes a 0.3 eV decrease in IP1 but has little effect on ∆IP
with respect to 3 and 5. The replacement of the benzene bridge
with a biphenyl bridge of approximately twice the length reduces
the ∆IP splitting, and thus the electronic coupling, in the same
proportion, i.e., by about 50%.
0.70 (0.45)d,e
0.53 (0.40)d,f
0.36 (0.20)d,g
0.31 (0.19)d,h
0.66
0.55
0.30
0.34
a Using the following values of L (obtained either from electronic
structure calculations or the fit of the profile of IV-CT bands): 0.3 eV
(3+),7 0.2 eV (4+),6b 0.15 eV (5+), and 0.15 eV (6+).6b,7 b Calculated taking
Hab ) Eop/2. c From ref 3f; in CH3CN. d Calculated according to eq 1 using
R ) R(TD-DFT) and, between parentheses, R ) R(NN); both R(TD-DFT)
and R(NN) values are obtained using the DFT-optimized geometries of the
radical-cation state.6,7 e From ref 7: Eop ) 1.50 eV; µab ) 8.1 D in
propionitrile/0.2 M (TBA)H. f From ref 5b: Eop ) 1.18 eV; µab ) 9.2 D in
CH2Cl2. g Present work: Eop ) 0.89 eV; µab ) 11.0 D in CH2Cl2. h From
ref 5b: Eop ) 0.79 eV; µab ) 11.6 D in CH2Cl2.
We have recently evaluated ∆0 (and thus ∆IP) using several
computational approaches: (a) from the energy spectrum
obtained at the TD-DFT level, (b) from DFT calculations by
computing the energy of the lowest state of each symmetry
representation, and (c) by means of Koopmans’ theorem (KT).6
It was found that the best agreement with experimental data
was obtained from KT-AM1 calculations; the KT-AM1 values
are, therefore, given in Table 1. In the context of electronic
structure calculations, it is also particularly interesting to
examine the reliability of the TD-DFT results, since this method
is currently widely used for excited-state computations. In
contrast to the simple semiempirical AM1 calculations, B3LYP-
TD-DFT overestimates the ∆0 value for all systems; for instance
in 5+ and 6+, the TD-DFT estimates are about 100% larger
than the experimental value. This result is in line with current
findings that, due to inherent problems in the applied standard
exchange-correlation functionals (see for instance ref 27 for
detailed explanations), TD-DFT yields substantial errors for
systems exhibiting charge-transfer states.27 We note, however,
that for all systems studied here, the error of the TD-DFT
estimate is systematic, being about 0.4 eV; this finding suggests
that it should be possible to obtain more accurate TD-DFT
estimates via an empirical correction.28
significantly lower than the estimates based on the direct
geometric distances between redox centers. For instance, it was
shown very recently3g that in delocalized nitro-centered MV
systems, R is only about 25-40% of the direct distance between
the nitrogen atoms, R(NN). Since for 3+-6+ there are no
experimental electron-transfer distances available, Table 2 gives
values of Hab obtained using the following: (i) TD-DFT
estimates for the ET distances R ) R(TD-DFT)6,7 and (ii) the
commonly applied approach, where R is approximated to the
nitrogen-nitrogen distance, R ) R(NN). In 3+-6+, the ratio
R(TD-DFT)/R(NN) is about 0.6;6,7 we expect, therefore, the Hab
estimates based on the R(TD-DFT) values to be more realistic.
The data in Table 2 indicate that the Hab values estimated
optically show trends similar to the ∆IP/2 values and are
comparable in magnitude. However, as expected from the
neglect of the interaction with symmetric vibrations (see Figure
2b), the UPS data yield values lower than those derived from
optical data. The strength of the electron-vibration couplings
with symmetric modes can actually be estimated from the data
for 1+ and 2+: since 1+ and 2+ are class-III systems, the
geometries of the ground states of both cation and neutral species
possess a symmetric configuration; the difference32 of 0.38-
0.44 eV between Eop and ∆IP can then be related to the
relaxation processes along symmetric vibrations.
We now turn to a comparison of the electronic couplings as
derived from UPS measurements and obtained from optical
data;29 see Table 2. Hab for 1+ and 2+ was determined as Hab
)
To obtain more accurate estimates of Hab from UPS data,
the ∆IP/2 values should be corrected to account for the
interaction with symmetric vibrations. The most accurate way
to accomplish such a correction would be to fit the shape of
the ionization spectrum to a vibronic coupling model. However,
a simple but effective estimate of the vibronic contribution can
Eop/2,30 while for systems 3+-6+ Hab was derived using eq 1.
In the latter cases, the question of the appropriate value of R
arises; several experimental and theoretical studies3g,6,7,31 indicate
that the actual electron-transfer distances in organic systems are
(25) Similar results have been recently found for other delocalized organic MV
systems (see refs 3a,g and 8).
be provided by correcting the UPS values according to Hab
)
(∆IP/2) + L, where L is the vibrational relaxation energy of
the symmetric modes associated with the excitation of the cation
from its ground state to its first excited state (see Figure 2b).
This approximation assumes that the minimum of the cation
(26) (a) Interestingly, the fact that the interaction with symmetric vibrations
could play a major role in the charger-transfer process in strongly
delocalized MV systems was also predicted by Hush.26b,c (b) Hush, N. S.
In Mixed-Valence Compounds: Theory and Applications in Chemistry,
Physics, Geology and Biology; Brown, D. B., Ed.; D. Reidel Academic
Publishers: Dordrecht, The Netherlands, 1980. (c) Reimers, R.; Hush, N.
S. Chem. Phys. 1996, 208, 177.
(27) (a) Drew, A.; Weisman, J. L.; Head-Gordon, M. J. Chem. Phys. 2003, 119,
2943. (b) Grimme, S.; Parac, M. ChemPhysChem 2003, 3, 292.
(28) (a) Recent28b investigations of (neutral) heterocyclic conjugated polymers
show that after an empirical correction via a linear regression, TD-DFT
yields accurate results for band gap. (b) Hutchison, G. R.; Ratner, M. A.;
Marks, T. J., J. Phys. Chem. A 2002, 106, 1059.
(29) It should be noted that the optical values apply to solution, whereas the
UPS values are determined in the gas phase. However, in strongly coupled
systems, the IV-CT bands show typically little dependence on solvent
polarity; this suggests that the medium has little effect on Hab and, therefore,
that solution and gas-phase values should be similar.
(30) (a) In the case of strongly coupled MV systems (class III), the energy of
the intervalence transition becomes a direct measure of the electronic
coupling: Eop ) 2Hab.10c The observed vibrational structure in the IV-CT
band of 1+ 3f and 2+ 30b and the absence of solvent dependence3f are features
characteristic of a class-III system. (b) Shida, T. Electronic Absorption
Spectra of Radical Ions; Elsevier: Amsterdam, 1988.
(31) (a) Nelsen, S. F.; Newton, M. D. J. Phys. Chem. A 2000, 104, 10023. (b)
Johnson, R. C.; Hupp, J. T. J. Am. Chem. Soc. 2001, 123, 2053.
(32) A similar difference between Eop and ∆IP values has been recently reported
for a delocalized phthalocyanine MV dimer: Binstead, R. A.; Reimers, J.
R.; Hush, N. S. Chem. Phys. Lett. 2003, 378, 654.
9
2730 J. AM. CHEM. SOC. VOL. 126, NO. 9, 2004