Sabine Richert1, Bart Limburg2, Harry L Anderson2, Christiane R Timmel1. 1. Centre for Advanced Electron Spin Resonance (CAESR), University of Oxford , South Parks Road, Oxford OX1 3QR, United Kingdom. 2. Chemistry Research Laboratory, Department of Chemistry, University of Oxford , 12 Mansfield Road, Oxford OX1 3TA, United Kingdom.
Abstract
The extent of triplet state delocalization is investigated in rigid linear zinc porphyrin oligomers as a function of interporphyrin bonding characteristics, specifically in meso-meso singly linked and β,meso,β fused structures, using electron paramagnetic resonance techniques. The results are compared with those of earlier measurements on porphyrin oligomers with alkyne linkers exhibiting different preferred conformations. It is shown that dihedral angles near 90° between the porphyrin planes in directly meso-to-meso linked porphyrin oligomers lead to localization of the photoexcited triplet state on a single porphyrin unit, whereas previous work demonstrated even delocalization over two units in meso-to-meso ethyne or butadiyne-bridged oligomers, where the preferred dihedral angles amount to roughly 30° and 0°, respectively. The triplet states of fused porphyrin oligomers (i.e., porphyrin tapes) exhibit extended conjugation and even delocalization over more than two porphyrin macrocycles, in contrast to meso-to-meso ethyne or butadiyne-bridged oligomers, where the spin density distribution in molecules composed of more than two porphyrin units is not evenly spread across the oligomer chain.
The extent of triplet state delocalization is investigated in rigid linear zinc porphyrin oligomers as a function of interporphyrin bonding characteristics, specifically in meso-meso singly linked and β,meso,β fused structures, using electron paramagnetic resonance techniques. The results are compared with those of earlier measurements on porphyrin oligomers with alkyne linkers exhibiting different preferred conformations. It is shown that dihedral angles near 90° between the porphyrin planes in directly meso-to-meso linked porphyrin oligomers lead to localization of the photoexcited triplet state on a single porphyrin unit, whereas previous work demonstrated even delocalization over two units in meso-to-meso ethyne or butadiyne-bridged oligomers, where the preferred dihedral angles amount to roughly 30° and 0°, respectively. The triplet states of fused porphyrin oligomers (i.e., porphyrin tapes) exhibit extended conjugation and even delocalization over more than two porphyrin macrocycles, in contrast to meso-to-meso ethyne or butadiyne-bridged oligomers, where the spin density distribution in molecules composed of more than two porphyrin units is not evenly spread across the oligomer chain.
Triplet excitons are
important intermediates in many processes
of biological or technological relevance.[1−5] The kinetics of their formation and their properties
have a major influence on the performance of devices relying on singlet
or triplet exciton fission or fusion, such as organic semiconductors.[6,7] Triplet formation is a common loss mechanism in organic solar cells.[8−10] Furthermore, the extent of delocalization of the photoexcited triplet
state has recently been shown to impact the performance of organic
photovoltaic devices and organic light emitting diodes.[11,12]A thorough understanding of the factors influencing triplet
state
delocalization could advance the development of structure–property
relations guiding the design of new functional materials. A targeted
modification of the triplet state properties should allow control
over the loss mechanisms and therefore enable enhanced device performance.In this study, the influence of the linker connecting the individual
units and the structural rigidity of the oligomeric chain on triplet
state delocalization is explored in linear singly linked and fused
porphyrin assemblies. The type of the linker determines the preferred
dihedral angle between the porphyrin planes, which approaches 90°
as the linker is shortened from a butadiyne group to a direct meso-to-meso bond.[13,14] At the same time, the increased steric hindrance upon shortening
of the linker renders the oligomer chain more rigid and reduces the
conformational spread.Exploiting the paramagnetic nature of
triplet states, electron
paramagnetic resonance (EPR) spectroscopy has been shown to be particularly
suitable for the study of triplet delocalization in porphyrin[15−19] and other[20] oligomers. If the orientation
of the zero-field splitting (ZFS) tensor is known, information about
delocalization in a series of oligomers can be obtained from changes
in the ZFS D-value by comparison with a suitable
monomer spectrum. The cube of the average electron–electron
distance, r, is inversely proportional to the D-value, related to the width of the transient EPR triplet
spectrum[19]where θ refers to the angle between
the dipolar Z-axis and the spin–spin vector, ge is the electronic g-factor,
βe is the Bohr magneton, μ0 is the
vacuum permeability, and the angular brackets indicate the expectation
value. The angle θ characterizes the shape of the spin density
distribution and is equal to 90° in the case of an oblate spin
density distribution and equal to 0° in the case of a prolate
one. As long as the spin distribution stays invariant and there are
no other significant contributions to the ZFS tensor apart from the
spin–spin contribution, increased delocalization should manifest
itself in a reduction of the width of the triplet EPR spectrum, since
the width of the spectrum is equal to 2|D|.Further information on triplet state delocalization can be obtained
from triplet state electron nuclear double resonance (ENDOR) measurements,
especially in cases where the orientation of the ZFS tensor is likely
to be different for individual molecules within a series.[18,21,22] The magnitudes of the measured
hyperfine coupling constants are directly proportional to the triplet
spin density. Therefore, if the triplet spin density is equally distributed
in an oligomer with N identical units, the hyperfine
coupling constants of the oligomer should be reduced by a factor of
1/N with respect to those of the monomer.Previous
EPR studies on triplet state delocalization in linear
ethyne and butadiyne-bridged zinc porphyrin oligomers have shown that
even delocalization of the triplet state is limited to about two porphyrins
units.[18,22] For even delocalization to be observed,
these two porphyrin units additionally need to be equivalent.[22] These experimental results from EPR and ENDOR
are in excellent qualitative agreement with density functional theory
(DFT) calculations (B3LYP/EPRII). For the longer oligomers, the calculations
reveal that the majority of the spin density is located on the central
porphyrin units of the structures, as illustrated in Figure for the example of meso-to-meso ethyne and butadiyne bridged
linear zinc porphyrin trimers. The corresponding ZFS and hyperfine
coupling parameters are however in agreement with a description of
Frenkel-type exciton delocalization (i.e., delocalization of the position
of the center of mass of the electron and hole of a strongly bound
electron–hole pair).[20,23]
Figure 1
Structures of previously
investigated[18,22] linear ethyne and butadiyne-bridged
zinc porphyrin trimers (left)
and corresponding spin density distributions calculated by DFT (right).
The spin density is unevenly distributed with the majority of the
spin density on the central porphyrin unit. “Ar” represents
a solubilizing aryl side group, 3,5-bis(2-(2-(2-methoxyethoxy)ethoxy)ethoxy)phenyl
(top) or 3,5-bis(trihexylsilyl)phenyl (bottom).
Structures of previously
investigated[18,22] linear ethyne and butadiyne-bridged
zinc porphyrin trimers (left)
and corresponding spin density distributions calculated by DFT (right).
The spin density is unevenly distributed with the majority of the
spin density on the central porphyrin unit. “Ar” represents
a solubilizing aryl side group, 3,5-bis(2-(2-(2-methoxyethoxy)ethoxy)ethoxy)phenyl
(top) or 3,5-bis(trihexylsilyl)phenyl (bottom).DFT calculations furthermore predict that the preferred dihedral
angle between two adjacent porphyrin units amounts to approximately
0° in butadiyne-bridged and roughly 30° in ethyne-bridged
oligomers. While the barrier to rotation of the porphyrin units about
the linker is comparatively low in butadiyne-bridged oligomers leading
to a broad distribution of conformations,[24] this conformational distribution is expected to be narrower for
ethyne-bridged oligomers due to a higher rotational barrier resulting
from increased steric constraints and increased π-conjugation.[25,26]Results from excitation wavelength dependent EPR measurements
on
butadiyne-bridged porphyrins[27] suggest
that the excited triplet state localizes on a single porphyrin unit
if conformations with interporphyrin dihedral angles near 90°
are predominantly probed. However, since the fraction of excited molecules
with angles close to 90° was low in these measurements and the
corresponding changes in signal intensity fairly small, only tentative
conclusions could be drawn.Here we use transient EPR and ENDOR
techniques to investigate triplet
state delocalization in a series of linear zinc porphyrin oligomers,
directly linked in meso position as shown in Figure . Monomer, dimer,
and trimer will be referred to as P1, P2, and P3, respectively. In these oligomers, the
interporphyrin dihedral angles have previously been shown to amount
to roughly 90° with very little conformational flexibility.[28−30] The results obtained for the singly linked oligomers will be contrasted
with those obtained for their fused counterparts, where two additional
covalent bonds in β-position introduce additional π-conjugated
connections and force the individual porphyrin macrocycles into a
planar, tapelike structure.
Figure 2
Chemical structures of the investigated linear
zinc porphyrin oligomers P1, P2, and P3. The porphyrin
side group “Ar” represents an aryl substituent, 3,5-bis(trihexylsilyl)phenyl,
which provides high solubility and prevents aggregation. Further details
on the synthesis and characterization of the compounds are given in
the SI.
Chemical structures of the investigated linear
zinc porphyrin oligomers P1, P2, and P3. The porphyrin
side group “Ar” represents an aryl substituent, 3,5-bis(trihexylsilyl)phenyl,
which provides high solubility and prevents aggregation. Further details
on the synthesis and characterization of the compounds are given in
the SI.
Results and Discussion
The transient cw EPR spectra of P1, P2, and P3 are shown in Figure (left). All EPR
data were acquired at X-band
(9.75 GHz) in frozen 2-methyltetrahydrofuran (MTHF) solution at 20
K. Unless otherwise stated, the samples were excited at 532 nm (∼2
mJ), using the second harmonic of a Nd:YAG laser operated at 20 Hz.
Further details on the sample preparation and experimental setup for
the EPR measurements are given in the Supporting Information (SI).
Figure 3
Overlay of the transient cw EPR spectra of the
singly linked compounds
(left), and visualization of the calculated hyperfine coupling tensors
and definition of the ZFS tensor axes in top view (center) and side
view (right). Proton hyperfine tensors are shown in green, whereas
nitrogen hyperfine tensors are depicted in blue. The field positions
corresponding to the X–, Y– and Z– orientations are indicated in the spectra. A large hyperfine coupling
exceeding −10 MHz is calculated for the proton in meso-position of the porphyrin core. This large coupling is omitted in
the figure for clarity. The EPR data were acquired in frozen MTHF
at 20 K.
Overlay of the transient cw EPR spectra of the
singly linked compounds
(left), and visualization of the calculated hyperfine coupling tensors
and definition of the ZFS tensor axes in top view (center) and side
view (right). Proton hyperfine tensors are shown in green, whereas
nitrogen hyperfine tensors are depicted in blue. The field positions
corresponding to the X–, Y– and Z– orientations are indicated in the spectra. A large hyperfine coupling
exceeding −10 MHz is calculated for the proton in meso-position of the porphyrin core. This large coupling is omitted in
the figure for clarity. The EPR data were acquired in frozen MTHF
at 20 K.Since the EPR spectra were recorded
using direct detection, a positive
signal corresponds to an absorptive (a) transition,
whereas negative signals result from emissive (e)
transitions. The observed spin polarization in P1, aaaeee, is consistent with
results obtained for similar zinc porphyrin monomers and indicates
a selective population of the out-of-plane triplet sublevel (Z) by intersystem crossing. A positive D-value is assumed, as confirmed for zinc porphyrin monomers with
terminal ethyne groups in the meso-positions by magnetophotoselection
experiments and DFT calculations.[21]The comparison of the triplet state EPR spectra of the singly linked
oligomers (Figure , left) shows that no significant changes in the ZFS parameter D or the relative triplet sublevel populations occur within
the series of molecules since the spectral width and the overall spectral
shape remain almost invariant. The ZFS parameters and relative triplet
state sublevel populations obtained from a numerical simulation of
the spectra are shown in the SI (Figure
S4). The similarity of the triplet spectra of P1, P2, and P3 might indicate that the delocalization
of the triplet excited state is similar in all of the investigated
singly linked compounds and does therefore not exceed a single porphyrin
unit. However, since the shape of the spin-density distribution in
the longer oligomers is a priori not known, a confirmation of this
tentative interpretation by pulse ENDOR spectroscopy, if attainable
(fast relaxation and/or short triplet state lifetimes are limiting
factors, see below), is desirable.The canonical orientations
in the spectra are assigned to the three
ZFS tensor axes in Figure under the assumption of a positive D-value.
At the field positions corresponding to these three orientations,
proton Mims ENDOR spectra were taken and are shown in the SI for reference. To interpret the orientation
selective ENDOR spectra and assign the couplings to the individual
protons of the structure, DFT calculations of the hyperfine coupling
tensors were performed using ORCA (B3LYP/EPRII).[31] Technical details on these calculations are provided in
the SI. A visualization of the calculated
proton and nitrogen hyperfine tensors for P1 is shown in Figure (center and right) together with an assignment of
the ZFS tensor axes to the molecular axes. The proton hyperfine tensors
are shown in green and their relative size reflects the magnitude
of the couplings. The couplings measured using ENDOR should correspond
to a projection of the hyperfine tensors onto the indicated ZFS tensor
axes. In our assignment, the Y axis corresponds to
the axis along the aryl groups, whereas Z is assigned
to the out-of-plane axis in analogy to similar porphyrins.[21] A pronounced proton coupling with a magnitude
of roughly −4 MHz is observed in the ENDOR spectra recorded
at the X- and Z-orientations (cf. SI, Figure S7, right) and corresponds to the
proton labeled H1 in Figure .To confirm the extent of triplet state delocalization
in the longer
oligomers, the ENDOR spectra recorded at the Z– orientation were compared. The detection of the ENDOR
couplings at the outmost position of the triplet spectrum (Z) has the advantage that no spurious contributions from
other orientations need to be considered. Furthermore, the comparatively
large coupling of ∼ –4 MHz resulting from H1 is quite pronounced at this orientation and should therefore
be easily comparable between the oligomers. An overlay of the corresponding
spectra for P1, P2, and P3 is shown in Figure (left). No significant change in the magnitude of the hyperfine
coupling of ∼ –4 MHz is observed, confirming
that the excited triplet state remains localized on a single porphyrin
unit in these singly linked oligomers where the dihedral angle between
adjacent porphyrin planes is close to 90°. DFT calculations of
the spin density distribution further confirm that the triplet state
localizes only on one of the porphyrin units in the dimer and trimer,
as shown in Figure (right).[32]
Figure 4
Comparison of the proton
Mims ENDOR spectra of P1, P2, and P3 recorded at magnetic field positions
corresponding to the Z– transition
and visualization of the spin density distributions in P1, P2, and P3 calculated by DFT.
The ENDOR spectra were recorded at 20 K and pulse lengths of 16 ns
were applied. To compensate for blind spots, three ENDOR spectra with
τ values of 120, 180, and 240 ns were averaged in all cases.
Although calculations predict the spin density to localize on the
central unit in P3 as shown
in the figure, solvent effects or changes of the geometry in solution
might lead to localization of the triplet state on any one of the
three units.
Comparison of the proton
Mims ENDOR spectra of P1, P2, and P3 recorded at magnetic field positions
corresponding to the Z– transition
and visualization of the spin density distributions in P1, P2, and P3 calculated by DFT.
The ENDOR spectra were recorded at 20 K and pulse lengths of 16 ns
were applied. To compensate for blind spots, three ENDOR spectra with
τ values of 120, 180, and 240 ns were averaged in all cases.
Although calculations predict the spin density to localize on the
central unit in P3 as shown
in the figure, solvent effects or changes of the geometry in solution
might lead to localization of the triplet state on any one of the
three units.Upon oxidation, the meso-meso singly linked porphyrin oligomers
can be transformed into triply
linked porphyrin tapes.[28] The structures
are shown in Figure (left), and the synthesis and characterization are presented in
the SI. Compared to the singly linked oligomers,
two additional covalent β–β bonds are formed between
the porphyrin macrocycles, introducing additional π-conjugated
connections and forcing the porphyrin π-systems into a coplanar
conformation with maximum π-overlap.
Figure 5
Chemical structures of
the fused porphyrin oligomers P1, , and (left), corresponding
transient cw EPR spectra
and simulations including simulation parameters (center), and visualization
of the calculated spin density distributions (right). Ar = 3,5-bis(trihexylsilyl)phenyl.
Chemical structures of
the fused porphyrin oligomers P1, , and (left), corresponding
transient cw EPR spectra
and simulations including simulation parameters (center), and visualization
of the calculated spin density distributions (right). Ar = 3,5-bis(trihexylsilyl)phenyl.The triplet yield of the fused
porphyrin dimer with terminal protons
in meso-positions has previously been shown to be
negligible.[33] To increase the triplet yield
for the measurements in this study, the series of fused porphyrin
oligomers was therefore prepared with −Br substituents in the meso-positions. The structures from monomer to trimer are
referred to as P1, , and . The triplet yield
for (with
different aryl side chains) has previously been estimated to amount
to roughly 0.12 by analysis of transient absorption spectra.[33] Due to the presence of the terminal −Br
substituents, the porphyrin in-plane triplet sublevels (X, Y) are primarily populated by intersystem crossing,
which results in the eeeaaa polarization pattern
observed in the triplet state cw EPR spectra of these compounds (cf. Figure , center).While the triplet yield, spin relaxation properties, and triplet
lifetime of P1 are similar
to those of a typical porphyrin monomer, the transient cw EPR experiments
on the fused porphyrins resulted in very low signal intensities due
to a very short triplet state lifetime[33] and possibly also fast spin relaxation. For these reasons, no pulse
EPR or ENDOR experiments could be carried out on and . Proton ENDOR experiments have
however been performed on the singly linked oligomers with terminal
−Br groups and are shown in the SI for reference. Compared to their counterparts with terminal protons,
the hyperfine spectra of these compounds are very similar and show
the same trends.For the series of fused porphyrin oligomers,
information about
delocalization could thus only be obtained by evaluation of the changes
in the ZFS parameter D. Since previous studies have
shown that precise information about the orientation of the ZFS tensor
is required in order to correctly interpret any changes in the spectral
width of the EPR spectra,[21,22] DFT calculations were
performed to support the experimentally observed trends. The latter
predict positive D-values for all investigated fused
porphyrin arrays, suggesting no marked changes of the orientation
of the ZFS tensor and thus the factor (1 – 3cos2 θ) in eq . If the orientation of the ZFS tensor remains invariant, delocalization
of the excited triplet state beyond a single unit should manifest
itself in a reduction of the width of the triplet EPR spectrum.The transient cw EPR spectra of , and , obtained when exciting the molecules
at their respective absorption maxima of 565 and 670 nm using an OPO
(optical parametric oscillator) system (2 mJ, 20 Hz, 5 ns pulse length),
are shown in Figure (center). To quantify the ZFS parameters as well as the relative
initial populations of the triplet sublevels, numerical simulations
were carried out using the MATLAB software package EasySpin[34] in combination with a home-written MATLAB fitting
routine. The triplet state parameters, D and E, as well as the relative populations of the triplet sublevels,
obtained as a result from the fit, are indicated in the figure.[35]A significant reduction in D from 1200 MHz to
about 590 MHz is observed experimentally between P1 and , suggesting increased delocalization of the
triplet state wavefunction. When the π-system is further extended
in , the D-value is decreased even further to about 370 MHz.To evaluate whether the decreased D-value is consistent
with an even delocalization of the triplet wavefunction, the experimentally
observed changes in D were contrasted with the calculated
values obtained from DFT. Only relative changes shall be compared
here, since it is known that DFT calculations generally underestimate
the ZFS parameters.[36] Comparisons of parameter
trends, however, have been shown to be reliable in previous studies
on similar systems.[18,22,27]The relative D-values in the series P1::, as calculated by DFT, amount to 1:0.58:0.33
and
correspond to the triplet spin density distributions shown in Figure (right). The absolute
values obtained from DFT are listed in the SI for reference. These calculated changes are in excellent agreement
with the experimental ratio of 1:0.49:0.31, indicating even delocalization
of the triplet state in fused porphyrin arrays as illustrated in Figure .In contrast
to meso-to-meso ethyne
or butadiyne-bridged porphyrins where uneven delocalization of the
triplet state was observed in oligomers with more than two porphyrin
units, this study showed that even delocalization can be achieved
in a fused, fully conjugated, porphyrin trimer. The observation of
even delocalization of the triplet spin density in tapelike porphyrin
molecules relates well to the extraordinary properties of their singlet
excited states such as absorption spectra reaching far into the near-infrared.[28]
Conclusions
In this study, the influence
of different bridges on triplet state
delocalization in linear zinc porphyrin oligomers was explored by
EPR spectroscopy. The results obtained on porphyrins linked by a direct meso-to-meso bond are consistent with previous
studies on ethyne and butadiyne-bridged oligomers[18,21,22,27] and clearly
indicate that the triplet state localizes on a single porphyrin unit
as the dihedral angle between adjacent porphyrin planes approaches
90°. No significant differences in the hyperfine coupling constants
measured by triplet state ENDOR are observed between P1 and the longer oligomers P2 and P3, showing that the extent of the triplet state wavefunction is approximately
identical within this series. The experimental results are fully consistent
with DFT calculations, suggesting that the triplet state localizes
on only one of the porphyrin units in P3.In linear porphyrin tapes, where the porphyrin units
are forced
to adopt a coplanar conformation by two additional covalent bonds
between the porphyrin macrocycles, even delocalization of the triplet
state was experimentally demonstrated by analysis of the zero-field
splitting parameters. Positive D-values were obtained
from DFT calculations for all molecules of the series P1, , and . With respect to the values obtained for the
monomer, P1, DFT predicts a
reduction in D by about 67% for . This is in excellent agreement
with the results from numerical simulations of the experimental spectra,
where a reduction by 69% was observed, and indicates even delocalization
of the triplet state wavefunction over the entire molecule.The even delocalization of the excited triplet state in fused porphyrin
oligomers was contrasted with the observed uneven delocalization in
ethyne or butadiyne-bridged oligomers with more than two porphyrin
units. Based on their singlet and triplet state properties, the investigated
linear, fully conjugated, porphyrin tapes seem to be ideal candidates
for molecular wires.[37] However, the investigation
of triplet state delocalization in these molecules is limited by small
triplet yields and short triplet state lifetimes.
Authors: Gita Sedghi; Louisa J Esdaile; Harry L Anderson; Santiago Martin; Donald Bethell; Simon J Higgins; Richard J Nichols Journal: Adv Mater Date: 2011-11-15 Impact factor: 30.849
Authors: Akshay Rao; Philip C Y Chow; Simon Gélinas; Cody W Schlenker; Chang-Zhi Li; Hin-Lap Yip; Alex K-Y Jen; David S Ginger; Richard H Friend Journal: Nature Date: 2013-08-07 Impact factor: 49.962
Authors: Claudia E Tait; Patrik Neuhaus; Martin D Peeks; Harry L Anderson; Christiane R Timmel Journal: Phys Chem Chem Phys Date: 2016-01-27 Impact factor: 3.676
Authors: Hannes Kraus; Michael C Heiber; Stefan Väth; Julia Kern; Carsten Deibel; Andreas Sperlich; Vladimir Dyakonov Journal: Sci Rep Date: 2016-07-06 Impact factor: 4.379
Authors: Gabriel Moise; Lara Tejerina; Michel Rickhaus; Harry L Anderson; Christiane R Timmel Journal: J Phys Chem Lett Date: 2019-09-12 Impact factor: 6.475
Authors: James O Thomas; Jakub K Sowa; Bart Limburg; Xinya Bian; Charalambos Evangeli; Jacob L Swett; Sumit Tewari; Jonathan Baugh; George C Schatz; G Andrew D Briggs; Harry L Anderson; Jan A Mol Journal: Chem Sci Date: 2021-07-26 Impact factor: 9.825
Authors: Bluebell H Drummond; Naoya Aizawa; Yadong Zhang; William K Myers; Yao Xiong; Matthew W Cooper; Stephen Barlow; Qinying Gu; Leah R Weiss; Alexander J Gillett; Dan Credgington; Yong-Jin Pu; Seth R Marder; Emrys W Evans Journal: Nat Commun Date: 2021-07-26 Impact factor: 14.919