Polarized Absorption in Crystalline Pentacene: Theory vs Experiment

Sep 16, 2015 - As shown in earlier works(-41, 42, 61) strong Frenkel/CT mixing can mimic Coulomb coupling in creating H- and J-like exciton bands defi...
1 downloads 11 Views 3MB Size
Article pubs.acs.org/JPCC

Polarized Absorption in Crystalline Pentacene: Theory vs Experiment N. J. Hestand,† H. Yamagata,† Bolei Xu,† Dezheng Sun,‡ Yu Zhong,§ Avetik R. Harutyunyan,∥ Gugang Chen,∥ Hai-Lung Dai,† Yi Rao,*,† and F. C. Spano*,† †

Department of Chemistry, Temple University, Philadelphia, Pennsylvania 19122, United States Department of Physics and §Department of Chemistry, Columbia University, New York, New York 10027, United States ∥ Honda Research Institute USA, Inc., Columbus, Ohio 43212, United States ‡

S Supporting Information *

ABSTRACT: The polarized absorption spectra of crystalline pentacene are obtained for excitation normal to the ab herringbone plane by measuring transmitted light in ultrathin crystals. The spectral line shapes for excitation polarized along b and orthogonal to b are analyzed theoretically using a Holstein-like Hamiltonian which includes both Frenkel and charge transfer (CT) excitons represented in a multiparticle basis set. The model agrees with prior estimates regarding the strong CT contribution (≈45%) of the exciton responsible for the b-polarized lower Davydov component. The polarization resolution allows one to also establish the nature of the upper Davydov component, which is found to contain far less CT content (≈15%), as well as the natures of the higher-energy vibronic excitons, which are found to consist of a complex mixture of Frenkel one- and two-particle states and CT excitons. Generally, the spectrum polarized along b displays J-aggregatelike vibronic signatures while the spectrum polarized orthogonal to b displays H-aggregate-like vibronic signatures. The assignment is entirely consistent with the calculated exciton band dispersions which agree well with the measured ones.

I. INTRODUCTION The oligoacene series continues to attract considerable attention for electronic applications such as field-effect transistors, due to their exceptionally high crystal-phase charge mobilities.1−5 More recently, tetracene and pentacene have come under intense scrutiny for their ability to promote rapid singlet exciton fission,6−17 potentially increasing the efficiency with which light can be converted to electricity in solar cell applications.12,13 Many groups have attempted to obtain a more fundamental understanding of the fission process and the potential role of intermediate charge transfer (CT) states, which necessarily involves a deep appreciation for the nature of the low-energy photoexcited states. In this regard, controversy has erupted over the CT exciton composition of the excited state responsible for the absorption origin (the lower Davydov component) in pentacene.18 Several theories predict a predominant CT admixture19−22 while others claim a substantial, but minority, component of 10−30%,18,23,24 an essentially even mixture of CT and Frenkel character,17,25 and even a negligible CT contribution.26 Experimental evidence for substantial Frenkel/CT exciton mixing near the absorption origin in pentacene has come mainly from measurements of the exciton band dispersion using electron energy loss spectroscopy by Knupfer and co-workers.27,28 Despite the heightened attention granted to pentacene photophysics, very few studies have reported its polarized absorption spectrum in the crystal phase.29−31 Although it is not difficult to grow pentacene single crystals with dimensions © 2015 American Chemical Society

greater than 1 mm by means of modern physical methods, a direct measurement of the absorption spectrum remains a challenge since at such thicknesses essentially all incident light is absorbed. In fact, the optimal thickness for obtaining the transmission UV−vis absorption spectrum is generally on the order of 200−400 nm, making it a formidable challenge to prepare high-quality, optically transparent single crystals of pentacene. In order to circumvent the strong absorbance of crystalline pentacene, light reflection or ellipsometry measurements have been employed.32−36 However, such methods allow only an indirect acquisition of the spectrum, as they require certain models for the extraction of the linear optical absorbance. In this paper, we present polarized absorption spectra of crystalline pentacene using direct transmission measurements. Thin single crystals of pentacene have been synthesized and irradiated normal to the ab herringbone plane which contains the most strongly coupled molecules. Because of the small area (19 000 cm−1), where extra peaks occur arising mostly from CT states. The dashed red and blue spectra in Figure 5 were obtained by multiplying the homogeneous line width by a factor of 5 for all transitions above 18 500 cm−1, where our calculations indicate a rapid increase in the density of states. The increased linewidth 22141

DOI: 10.1021/acs.jpcc.5b07163 J. Phys. Chem. C 2015, 119, 22137−22147

Article

The Journal of Physical Chemistry C

importance of separated vibronic−vibrational excitations in the photophysics of pentacene. In marked contrast to the b-polarized spectrum, the first transition (A1) in the spectrum polarized mainly along the aaxisthe upper Davydov componentarises from an excited state which is dominated by one-particle Frenkel excitons, with a small CT admixture of less than 15%. The one-particle Frenkel excitations also play a majority role in the next higherenergy transition responsible for A2, although the CT admixture rises to almost 40%. Interestingly, the peak-to-peak separation between A1 and A2 is close to a vibrational quantum, showing that the excited states polarized ⊥b are more monomer-like than those polarized along b. The cluster of low-intensity transitions at 18 000 cm−1 (labeled A2′ in Figure 6) is responsible for the pronounced shoulder in the observed spectrum (see Figure 1b). Such states are dominated by CT excitons and two-particle Frenkel excitons. One of the most striking differences between the spectra polarized along b and orthogonal to b is the ratio of oscillator strengths involving the first two vibronic peakssee Figure 6. The measured B1/B2 ratio is ≈2.3, much larger than the 0−0/ 0−1 ratio of ≈1.2 in the monomer (solution) spectrum in Figure 1a. By contrast, the A1/A2 ratio is ≈0.75, significantly smaller than in the monomer spectrum. In past works involving Coulombically coupled Frenkel excitons we showed that a ratio increase (relative to the monomer) signifies J-aggregation, while a ratio decrease signifies H-aggregation.48 Hence, for pentacene, the b-polarized spectrum is J-like, while the ⊥b-polarized spectrum is H-like. Further clues into the J/H-aggregate nature of the polarized components can be obtained from the energy dispersions of the bands which contain the two Davydov components. We therefore evaluated the two lowest energy bands in crystalline pentacene within the reciprocal space stemming from the herringbone plane. In triclinic pentacene, the reciprocal lattice vectors Kb ≡ 2πa × c/(a × c·b) and Ka ≡ 2πb × c/(b × c·a) approximately point in the directions of b and a, respectively. For a herringbone layer containing N by N unit cells, the exciton wave vectors take on the values k = ka + kb with Kb = 0, ±Kb/N, ±2Kb/N, ..., Kb/2 and Ka = 0, ±Ka/N, ±2Ka/N, ..., Ka/ 2. Figure 7 shows the calculated energy bands; the point ka = kb = 0 on the red surface locates the energy of the optically allowed exciton polarized along b (and responsible for peak B1) while the corresponding point on the blue surface locates the optically allowed exciton polarized along the direction orthogonal to b (and responsible for peak A1). The energetic separation between the two surfaces at ka = kb = 0 is therefore the Davydov splitting. In Figure 7a only the Coulombic coupling is retained, by setting all electron and hole transfer integrals in the pentaceneparametrized Hamiltonian to zero, while in Figure 7b charge transfer is activated, corresponding to the complete Hamiltonian for pentacene used to evaluate the absorption spectra in Figure 5. In the presence of Coulomb coupling alone the DS is almost negligibleless than 100 cm−1and the polarizations of the two components are opposite to that observed experimentally. This was also found to be the case in tetracene.25 Activating Frenkel/CT mixing reverses the polarization assignments and dramatically increases the bandwidths. A large positive curvature develops in the lower band, mirrored by a large negative curvature induced in the upper band. The DS substantially increases to approximately 1000 cm−1, in excellent agreement with the measured value.29−31 The DS in

mimics enhanced population relaxation from the higher energy, primarily CT excitons. A further discrepancy between the experimental and calculated spectra is the ratio of the bpolarized to the ⊥b-polarized peak intensities (A1/B1), presumably due to our utilization of an isotropic dielectric constant. In order to further access the nature of the transitions responsible for the polarized spectra, Figure 5 also shows the polarized spectra calculated purely electronically (i.e., without vibronic coupling) by setting all HR factors to zero. The importance of vibronic coupling becomes immediately apparent. In particular, the DS more than doubles when vibronic coupling is neglected.

IV. NATURE OF THE TRANSITIONS The nature of the various transitions observed in the polarized absorption spectra of crystalline pentacene can be appreciated from Figure 6, which shows the dominant transitions

Figure 6. Calculated stick spectrum for excitation polarized along b (top) and normal to b (bottom) in crystalline pentacene.

represented in a bar spectrum, where the bar height is proportional to the oscillator strength of each transition. The relative CT contribution is indicated by the red shaded region of each bar, while the Frenkel contribution is divided between one-particle (black) and two-particle (blue) components. As is immediately apparent, the excited state responsible for the dominant b-polarized B1 transitionthe lower Davydov componenthas a large CT admixture, approximately 45%, in agreement with our earlier estimates.17,25 The Frenkel component is largely due to one-particle states. The higher energy transitions responsible for peaks B2 and B3 have a more complex composition. In particular, the main transition within the B2 cluster arises from an excited state with dominant CT admixtureroughly 65%with the smaller Frenkel contribution due almost entirely to two-particle states. Such states also contribute strongly to the main B3 transition, highlighting the 22142

DOI: 10.1021/acs.jpcc.5b07163 J. Phys. Chem. C 2015, 119, 22137−22147

Article

The Journal of Physical Chemistry C

Figure 8. Calculated energy dispersion parallel to the Ka (blue line) and Kb (red line) reciprocal lattice directions compared with the experimentally measured dispersions from ref 28 using EELS. The experimental data is shown as the error bars of the measurement. Note that in ref 28 the a and b crystallographic axis labels are interchanged relative to our definitions. Hence, a* in ref 28. corresponds to our Kb, while b* corresponds to our Ka.

V. NATURE OF THE EXCITON DISPERSION In this section we examine the origin of the band curvatures of the two lowest energy surfaces in crystalline pentacene. The band curvatures ultimately derive from the mixing between Frenkel and CT excitons as demonstrated in Figure 7, with practically no role played by Coulomb coupling. The essential physics responsible for the band shapes of Figures 7 and 8 is present in a much simpler model consisting of a onedimensional (1D) crystal with two molecules per unit cell, in which only the nearest-neighbor electron and hole transfer is included while Coulomb coupling and vibronic coupling are entirely neglected; see Figure 9. The 1D crystal has translational symmetry as well as a 2-fold screw axis along the aggregate axis (“b”) which interchanges the two sublattices. In the following analysis, the basis set is limited to Frenkel and nearest-neighbor CT excitons. The energy of a diabatic CT excitation is taken to be ω0−0 + ECT, where ω0−0 is the molecular S0 → S1 transition energy (recall ℏ is set to unity) . In the simple model outlined above, the diabatic Frenkel and CT exciton bands are flat and dispersionless. There are two Frenkel bands (red and blue dashed lines in Figure 9); one symmetric (red) and the other antisymmetric (blue) with respect to the screw rotation. The four CT bands, which are higher in energy by ECT, are not shown in Figure 9. The transition to the k = 0 state in the symmetric FE band is allowed and b-polarized, while the transition to the k = 0 state in the antisymmetric FE band is allowed and a-polarized. Activating te and th allows mixing to occur between the Frenkel and CT excitons of like symmetry, resulting in the bands shown by the solid curves in Figure 9. The excition with wave vector k in the symmetric (antisymmetric) band has energy given by

Figure 7. Calculated energy dispersion for the two lowest energy excitons in pentacene for directions within the Herringbone plane. In (a) Frenkel/CT exciton mixing is excluded by setting all CT integrals to zero. In (b) Frenkel/CT mixing is included. All other parameters are the same as in Figure 5.

pentacene is considerably greater than the value measured for tetracene (≈600 cm−1)60 due to the enhanced Frenkel/CT mixing. The calculated band dispersions in Figure 7b agree well with the measured momentum-dependent exciton energies obtained using EELS by Roth et al.,28 as demonstrated in Figure 8. The positive curvature of the lower band is characteristic of Jaggregates, since the point ka = kb = 0 defines the band minimum, and is the source of the J-like behavior observed in the b-polarized spectrum (enhanced B1/B2 ratio vs the monomer). Conversely, the negative curvature of the upper band, which is characteristic of H-aggregates, is responsible for the H-like behavior (attenuated A1/A2 ratio) for the spectrum polarized orthogonal to b. Note that although B1 is strongly red-shifted relative to the monomer (see Figure 1), consistent with its J-like behavior, the peak A1 is also red-shifted, in contrast to the blue shift expected for H-like behavior. This unusual feature is due to the manner in which the negative curvature of the upper surface develops in pentacene. In conventional Coulomb-coupled H-aggregates, the negative curvature arises because the energy of the k = 0 exciton is shifted positively relative to all other excitons. However, for pentacene as well as CT-mediated H-aggregates61 the negative curvature is instead the result of a strong negative shif t in the energies of the high wave vector excitons. This will be explored in greater detail in section V.

ES(AS)(k) = ω0 − 0 + ΔS(AS)(k)

(4a)

where the k-dependent energy shifts are given by 22143

DOI: 10.1021/acs.jpcc.5b07163 J. Phys. Chem. C 2015, 119, 22137−22147

Article

The Journal of Physical Chemistry C

th have the same sign, as in pentacene, the lower band is symmetric and the lower Davydov component is b-polarized. By contrast, opposite signs for te and th result in the lower band being antisymmetric, and the lower Davydov component being a-polarized. Hence, the polarization of the two Davydov components depends critically on the relative sign between te and th. We have confirmed the sign sensitivity in the polarization ordering of the Davydov components for the pentacene-parametrized Hamiltonian in section III: imposing opposite signs for te(l) and th(l) for l = ±(1/2,1/2) and for l = ±(1/2,−1/2) interchanges the polarizations of the two Davydov components. The resulting spectra are shown in Figure S.5 in the Supporting Information.

VI. DISCUSSION/CONCLUSION We have resolved the polarization dependence of the absorption spectrum in crystalline pentacene with excitation normal to the 001 herringbone plane, and subsequently analyzed the spectral line shapes using a theory based on vibronically coupled Frenkel and CT excitons. Overall, the calculated line shapes normal to and along the crystallographic b-axis are in very good agreement with experiment. Our combined experimental/theoretical approach yielded detailed information on the nature of several of the low-lying optical excitations. The high CT content (about 45%) for the exciton responsible for the b-polarized lower Davydov component (or B1) is consistent with previous analyses based on the unpolarized absorption spectrum.17,25 The polarization resolution of the current work allowed us to also probe the nature of the exciton responsible for the mainly a-polarized upper Davydov component (A1). In marked contrast to the B1 exciton, the exciton responsible for A1 has only a small CT admixture, ≈15% . The measured DS of ≈1000 cm−1 is accurately captured by the theory and is almost entirely due to Frenkel/CT exciton mixing, consistent with our earlier analysis.25 Coulomb coupling plays a very minor role (see Figure 7a), as is also the case for tetracene,25 since the oscillator strength of the short-axis polarized S0 → S1 transition is very weak.43 The theory also quantitatively accounts for the spectral positions and oscillator strengths of several of the higher-energy a- and b-polarized excitations arising, in part, from vibronic coupling involving the ubiquitous totally symmetric ringstretching mode at approximately 1400 cm−1. Such states are comprised of a complex mixture of one- and two-particle Frenkel excitons as well as CT excitons. Unlike in anthracene and tetracene, the b-polarized progression in pentacene (B1, B2, B3) is highly irregular and grossly distorted from the regular progression observed for the unaggregated monomer in solution (see Figure 1). The enhanced distortion and greater DS exhibited by pentacene is due to greater Frenkel/CT mixing. The latter is due to the relative ease for dissociating local Frenkel excitons into ion pair states in pentacene vs tetracene and anthracene. Generally, this is due to the diminishing value of Ip − EA (Ip = ionization potential, EA = electron affinity) expected with increasing oligoacene length. As estimated early on using electroabsorption spectroscopy, the lowest CT exciton in pentacene is approximately 0.29 eV higher in energy than the lowest (mainly) Frenkel exciton, considerably smaller than the corresponding values for anthracene (0.43 eV) and tetracene (0.35 eV).58,63 Indeed, Figure 6a shows a cluster of predominantly CT states approximately 0.25 eV above the lower Davydov component, consistent with the assignments

Figure 9. (top) One-dimensional crystal lattice with two molecules per unit cell. Arrows indicate the short-axis polarized transition dipole moments. The directions shown correspond to the b-polarized transition. (bottom) Dispersion curves for the diabatic Frenkel bands and for the lowest two excitons after Frenkel/CT exciton mixing. The higher-energy diabatic CT exciton bands with energy ECT = 800 cm−1 is not shown. The case shown corresponds to te and th having the same sign, with te = th = 200 cm−1.

ECT 2

ΔS(k) =



2(te 2 + th 2) + 4teth cos(k /2) + ECT 2 /4 (4b)

ΔAS(k) =

ECT 2 −

2(te 2 + th 2) − 4teth cos(k /2) + ECT 2 /4 (4c)

Here, the signs of te and th are consistent with the phase convention in which the HOMO (LUMO) on one sublattice is obtained from the HOMO (LUMO) on the other sublattice by the 2-fold screw rotation. Figure 9 shows that Frenkel/CT mixing results in a lower band with positive curvature and an upper band with negative curvature, similar to the case of pentacene in Figures 7 and 8. Interestingly, the upper band develops negative curvature, as defines H-aggregates, because the downward energy shift increases with exciton wave vector. In this simple model, where te and th were taken to be equal, the k = 0 exciton in the upper (antisymmetric) band undergoes no energy shift at all, as can be shown directly from eq 4c. However, if we were to make te and th unequal in magnitude, the k = 0 exciton in the upper band would undergo a red shift since ΔAS(k=0) < 0. This mechanism of creating the H-like band61 is therefore entirely different than the conventional mechanism of Kasha which is predicated entirely on Coulomb coupling.62 In the simple 1D aggregate of Figure 9, which band is symmetric or antisymmetric depends on the relative sign of te and th. One can readily appreciate from eqs 4a−4c that if te and 22144

DOI: 10.1021/acs.jpcc.5b07163 J. Phys. Chem. C 2015, 119, 22137−22147

Article

The Journal of Physical Chemistry C

the danger in assigning H- and J-like behavior based solely on spectral shifts. In the present case, the red shift of A1 is due to the peculiar nature of CT-mediated H-aggregates,61 as discussed in section V. We have found that in pentacene, and oligoacenes in general, the polarizations of the lower and upper Davydov components depend sensitively on the relative signs of the electron and hole transfer integrals connecting molecules in the two (inequivalent) sublattices. This is most easily appreciated for the simple 1D crystal in Figure 9 using eqs 4a−4c. When teth > 0, the lower Davydov component is b-polarized and the upper Davydov component is polarized perpendicular to b. Reversing the sign of teth interchanges the polarizations of the two components. The relationship between the band polarizations and the sign of teth assumes a phase convention for the HOMOs and LUMOs in the two sublattices based on a symmetry operation which exchanges the two sublattices but leaves the b-polarized transition dipole moment unchanged, such as the 2-fold screw rotation in Figure 9. In triclinic pentacene the two molecules in the unit cell are not related by a symmetry operation. However, one can still perform a 2-fold screw rotation along b to determine the MO phases and therefore the relative signs of te and th as described in ref 25 (see also the Supporting Information). In triclinic pentacene, Table 1 shows that te and th in the (1/2,1/2) direction have the same sign, as is also the case in the (1/2,−1/2) direction. Hence, the lower Davydov component must be b-polarized, as is found in experiment. When Coulomb coupling is dominant (which is not the case for pentacene), the sign of the coupling dictates the polarization assignments of the Davydov components. This has been discussed previously for chiral PDI bichromophores which are arranged in a way which resembles the molecular orientations present in a single unit cell of an oligoacene herringbone aggregate.69 The anisotropic properties of the low-lying excitations in crystalline pentacene should have immediate impact on singlet fission. In ref 17 it was shown that a similarly parametrized Frenkel/CT exciton Hamiltonian as the one employed here supports significant coupling between CT states and the multiexciton (“TT”) states, without disturbing the unpolarized absorption spectral line shape. With the more accurate parametrization obtained by reproducing both polarization components in hand, further investigation of singlet fission in pentacene is warranted, and will form the basis of future investigations.

made in ref 58. The correspondingly larger Frenkel/CT exciton mixing in pentacene is accommodated in the current model by allowing for extended charge-separated states with electron− hole separations out to 20 Å.17 Notably, in ref 25 only nearestneighbor CT states were incorporated, necessitating an unphysically small energy for ionizing a local Frenkel exciton into the (1/2,1/2) charge-separated pair. Interestingly, several groups have accounted for the satellite bands (B2,B3) in the b-polarized spectrum by considering only electronic excitations (Frenkel and CT) via the Bethe−Salpeter equation, and have attained satisfactory agreement with experiment.19−22 This has also been shown for TIPSpentacene.64 Hence, one may question the importance of vibronic coupling in accounting for B2 and B3. Figure 5a demonstrates that after removing vibronic coupling the bands are significantly attenuated but retain roughly the same spectral positions as observed in experiment, in rough agreement with refs 19−22. This is not the case, however, for the component polarized perpendicular to b. Without vibronic coupling the spectrum in Figure 5b no longer resembles the experimental spectrum; the majority of the ⊥b oscillator strength is retained by the lowest energy excitation (A1), which is so strongly blueshifted that the DS is roughly double the measured value. The greater impact of vibronic coupling on the ⊥b spectrum is consistent with the greater Frenkel character of the A1 and A2 peaks. Hence, vibronic coupling is necessary for achieving quantitative agreement with experiment. The higher energy B2, B3, ... and A2, A2′, ... transitions are a complex mixture of vibronic satellites and excitons of mixed Frenkel and CT character, in agreement with the assessment of Roth et al.28 As shown in earlier works41,42,61 strong Frenkel/CT mixing can mimic Coulomb coupling in creating H- and J-like exciton bands defined by the ordering of the excitons: in J-like exciton bands the lowest state has wave vector k = 0 consistent with a positive band curvature, while in H-like bands the highest exciton has k = 0 and the band curvature is negative. The greatly enhanced ratio of the B1 to B2 oscillator strengths observed in the b-polarized spectrum of pentacene, more than a factor of 2 larger than the 0−0/0−1 ratio for the monomer in solution, as well as the large red shift of B1, are spectral signatures of J-aggregation48 and arise from the positive curvature of the lowest-energy exciton band. The calculated curvature is indeed positive and the overall band dispersion is in excellent agreement with that measured by Roth et al. using EELS28 (see Figure 8). Furthermore, a large positive curvature of the lowest energy band also predicts J-like characteristics for the photoluminescence (PL) spectrum,48,52 namely a large 0− 0/0−1 ratio which decreases with temperature, as is the case for tetracene.65−67 In this regard, Anger et al. has shown that PL arising from the lower Davydov exciton in pentacene is characterized by a large 0−0/0−1 ratio (>10) which decreases precipitously with temperature.68 In stark contrast, the (mainly) a-polarized spectrum in pentacene maintains a progression in which the A1/A2 ratio is suppressed relative to the 0−0/0−1 ratio for a monomer in solution, as is characteristic of an H-aggregate.48 This arises from the negative curvature of the band containing the mainly apolarized upper Davydov component (see Figure 7b), which agrees well with the EELS measurements of Roth et al.28 We point out that although the lower Davydov component (B1) is red-shifted from the 0−0 peak in solution, as expected for Jaggregates, the upper Davydov component (A1) also remains slightly red-shifted, despite its H-like nature. This underscores



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.jpcc.5b07163. More complete descriptions of the experiment and pentacene samples; additional information about the parametrization of the pentacene Hamiltonian; details regarding the absorption spectrum calculations and phase convention used in choosing the signs of the electron and hole transfer integrals (PDF)



AUTHOR INFORMATION

Notes

The authors declare no competing financial interest. 22145

DOI: 10.1021/acs.jpcc.5b07163 J. Phys. Chem. C 2015, 119, 22137−22147

Article

The Journal of Physical Chemistry C



(19) Tiago, M. L.; Northrup, J. E.; Louie, S. G. Ab Initio Calculation of the Electronic and Optical Properties of Solid Pentacene. Phys. Rev. B: Condens. Matter Mater. Phys. 2003, 67, 115212. (20) Ambrosch-Draxl, C.; Nabok, D.; Puschnig, P.; Meisenbichler, C. The Role of Polymorphism in Organic Thin Films: Oligoacenes Investigated from First Principles. New J. Phys. 2009, 11, 125010. (21) Cudazzo, P.; Gatti, M.; Rubio, A.; Sottile, F. Frenkel Versus Charge-Transfer Exciton Dispersion in Molecular Crystals. Phys. Rev. B: Condens. Matter Mater. Phys. 2013, 88, 195152. (22) Sharifzadeh, S.; Darancet, P.; Kronik, L.; Neaton, J. B. LowEnergy Charge-Transfer Excitons in Organic Solids from FirstPrinciples: The Case of Pentacene. J. Phys. Chem. Lett. 2013, 4, 2197−2201. (23) Petelenz, B.; Petelenz, P.; Shurvell, H. F.; Smith, V. H. Reconsideration of the Electroabsorption Spectra of the Tetracene and Pentacene Crystals. Chem. Phys. 1988, 119, 25−39. (24) Petelenz, P.; Slawik, M.; Yokoi, K.; Zgierski, M. Z. Theoretical Calculation of the Electroabsorption Spectra of Polyacene Crystals. J. Chem. Phys. 1996, 105, 4427−4440. (25) Yamagata, H.; Norton, J.; Hontz, E.; Olivier, Y.; Beljonne, D.; Bredas, J. L.; Silbey, R. J.; Spano, F. C. The Nature of Singlet Excitons in Oligoacene Molecular Crystals. J. Chem. Phys. 2011, 134, 204703. (26) Zimmerman, P. M.; Bell, F.; Casanova, D.; Head-Gordon, M. Mechanism for Singlet Fission in Pentacene and Tetracene: From Single Exciton to Two Triplets. J. Am. Chem. Soc. 2011, 133, 19944− 19952. (27) Schuster, R.; Knupfer, M.; Berger, H. Exciton Band Structure of Pentacene Molecular Solids: Breakdown of the Frenkel Exciton Model. Phys. Rev. Lett. 2007, 98, 037402. (28) Roth, F.; Schuster, R.; Konig, A.; Knupfer, M.; Berger, H. Momentum Dependence of the Excitons in Pentacene. J. Chem. Phys. 2012, 136, 204708. (29) Zanker, V.; Preuss, J. Microspectrometrical Measurements on Single Crystals of Tetracene and Pentacene in Depedence of Thickness and Temperature. Z. Angew. Phys. 1969, 27, 363−365. (30) Aoki-Matsumoto, T.; Furuta, K.; Yamada, T.; Moriya, H.; Mizuno, K.; Matsui, A. H. Excitonic Photoluminescence in Pentacene Single Crystal. Int. J. Mod. Phys. B 2001, 15, 3753−3756. (31) Prikhotko, A. F.; Tsikora, L. I. Spectral Investigations of Pentacene. Opt. Spectrosc. USSR 1968, 25, 242. (32) Grobosch, M.; Schuster, R.; Pichler, T.; Knupfer, M.; Berger, H. Analysis of the Anisotropy of Excitons in Pentacene Single Crystals Using Reflectivity Measurements and Electron Energy-Loss Spectroscopy. Phys. Rev. B: Condens. Matter Mater. Phys. 2006, 74, 155202. (33) Dressel, M.; Gompf, B.; Faltermeier, D.; Tripathi, A. K.; Pflaum, J.; Schubert, M. Kramers-Kronig-Consistent Optical Functions of Anisotropic Crystals: Generalized Spectroscopic Ellipsometry on Pentacene. Opt. Express 2008, 16, 19770−19778. (34) Hinderhofer, A.; Heinemeyer, U.; Gerlach, A.; Kowarik, S.; Jacobs, R. M. J.; Sakamoto, Y.; Suzuki, T.; Schreiber, F. Optical Properties of Pentacene and Perfluoropentacene Thin Films. J. Chem. Phys. 2007, 127, 194705. (35) Faltermeier, D.; Gompf, B.; Dressel, M.; Tripathi, A. K.; Pflaum, J. Optical Properties of Pentacene Thin Films and Single Crystals. Phys. Rev. B: Condens. Matter Mater. Phys. 2006, 74, 125416. (36) Qi, D. C.; Su, H. B.; Bastjan, M.; Jurchescu, O. D.; Palstra, T. M.; Wee, A. T. S.; Rubhausen, M.; Rusydi, A. Observation of Frenkel and Charge Transfer Excitons in Pentacene Single Crystals Using Spectroscopic Generalized Ellipsometry. Appl. Phys. Lett. 2013, 103, 113303. (37) Hoffmann, M.; Soos, Z. G. Optical Absorption Spectra of the Holstein Molecular Crystal for Weak and Intermediate Electronic Coupling. Phys. Rev. B: Condens. Matter Mater. Phys. 2002, 66, 024305. (38) Hoffmann, M.; Schmidt, K.; Fritz, T.; Hasche, T.; Agranovich, V. M.; Leo, K. The Lowest Energy Frenkel and Charge-Transfer Excitons in Quasi-One-Dimensional Structures: Application to Meptcdi and Ptcda Crystals. Chem. Phys. 2000, 258, 73−96.

ACKNOWLEDGMENTS F.C.S. is supported by the National Science Foundation, Grant DMR-1505437. Y.R. thanks Dr. Colin Nuckolls for his support in the preparation of pentacene samples and Yajing Wu for her assistance in the experimental setup. The experimental work was supported by the Honda Research Institute, USA.



REFERENCES

(1) Zaumseil, J.; Sirringhaus, H. Electron and Ambipolar Transport in Organic Field-Effect Transistors. Chem. Rev. 2007, 107, 1296−1323. (2) Jurchescu, O. D.; Baas, J.; Palstra, T. T. M. Effect of Impurities on the Mobility of Single Crystal Pentacene. Appl. Phys. Lett. 2004, 84, 3061−3063. (3) Dimitrakopoulos, C. D.; Malenfant, P. R. L. Organic Thin Film Transistors for Large Area Electronics. Adv. Mater. 2002, 14, 99−117. (4) Lin, Y. Y.; Gundlach, D. J.; Nelson, S. F.; Jackson, T. N. Stacked Pentacene Layer Organic Thin-Film Transistors with Improved Characteristics. IEEE Electron Device Lett. 1997, 18, 606−608. (5) Klauk, H.; Halik, M.; Zschieschang, U.; Schmid, G.; Radlik, W.; Weber, W. High-Mobility Polymer Gate Dielectric Pentacene Thin Film Transistors. J. Appl. Phys. 2002, 92, 5259−5263. (6) Jadhav, P. J.; Brown, P. R.; Thompson, N.; Wunsch, B.; Mohanty, A.; Yost, S. R.; Hontz, E.; Van Voorhis, T.; Bawendi, M. G.; Bulovic, V.; et al. Triplet Exciton Dissociation in Singlet Exciton Fission Photovoltaics. Adv. Mater. 2012, 24, 6169−6174. (7) Smith, M. B.; Michl, J. Recent Advances in Singlet Fission. Annu. Rev. Phys. Chem. 2013, 64, 361−386. (8) Smith, M. B.; Michl, J. Singlet Fission. Chem. Rev. 2010, 110, 6891−6936. (9) Burdett, J. J.; Bardeen, C. J. The Dynamics of Singlet Fission in Crystalline Tetracene and Covalent Analogs. Acc. Chem. Res. 2013, 46, 1312−1320. (10) Burdett, J. J.; Bardeen, C. J. Quantum Beats in Crystalline Tetracene Delayed Fluorescence Due to Triplet Pair Coherences Produced by Direct Singlet Fission. J. Am. Chem. Soc. 2012, 134, 8597−8607. (11) Chan, W. L.; Ligges, M.; Zhu, X. Y. The Energy Barrier in Singlet Fission Can Be Overcome through Coherent Coupling and Entropic Gain. Nat. Chem. 2012, 4, 840−845. (12) Rao, A.; Wilson, M. W. B.; Hodgkiss, J. M.; Albert-Seifried, S.; Bässler, H.; Friend, R. H. Exciton Fission and Charge Generation Via Triplet Excitons in Pentacene/C-60 Bilayers. J. Am. Chem. Soc. 2010, 132, 12698−12703. (13) Congreve, D. N.; Lee, J. Y.; Thompson, N. J.; Hontz, E.; Yost, S. R.; Reusswig, P. D.; Bahlke, M. E.; Reineke, S.; Van Voorhis, T.; Baldo, M. A. External Quantum Efficiency above 100% in a Singlet-ExcitonFission-Based Organic Photovoltaic Cell. Science 2013, 340, 334−337. (14) Thorsmølle, V. K.; Averitt, R. D.; Demsar, J.; Smith, D. L.; Tretiak, S.; Martin, R. L.; Chi, X.; Crone, B. K.; Ramirez, A. P.; Taylor, A. J. Morphology Effectively Controls Singlet-Triplet Exciton Relaxation and Charge Transport in Organic Semiconductors. Phys. Rev. Lett. 2009, 102, 017401. (15) Chan, W. L.; Berkelbach, T. C.; Provorse, M. R.; Monahan, N. R.; Tritsch, J. R.; Hybertsen, M. S.; Reichman, D. R.; Gao, J. L.; Zhu, X. Y. The Quantum Coherent Mechanism for Singlet Fission: Experiment and Theory. Acc. Chem. Res. 2013, 46, 1321−1329. (16) Yost, S. R.; Lee, J.; Wilson, M. W. B.; Wu, T.; McMahon, D. P.; Parkhurst, R. R.; Thompson, N. J.; Congreve, D. N.; Rao, A.; Johnson, K.; et al. Transferable Model for Singlet-Fission Kinetics. Nat. Chem. 2014, 6, 492−497. (17) Beljonne, D.; Yamagata, H.; Bredas, J. L.; Spano, F. C.; Olivier, Y. Charge-Transfer Excitations Steer the Davydov Splitting and Mediate Singlet Exciton Fission in Pentacene. Phys. Rev. Lett. 2013, 110, 226402. (18) Pac, B.; Petelenz, P. Lowest Singlet Exciton in Pentacene: Modern Calculations Versus Classic Experiments. ChemPhysChem 2014, 15, 2801−2809. 22146

DOI: 10.1021/acs.jpcc.5b07163 J. Phys. Chem. C 2015, 119, 22137−22147

Article

The Journal of Physical Chemistry C

(61) Yamagata, H.; Pochas, C. M.; Spano, F. C. Designing J- and HAggregates through Wave Function Overlap Engineering: Applications to Poly(3-Hexylthiophene). J. Phys. Chem. B 2012, 116, 14494−14503. (62) Kasha, M. Energy Transfer Mechanisms and the Molecular Exciton Model for Molecular Aggregates. Radiat. Res. 1963, 20, 55−70. (63) Sebastian, L.; Weiser, G.; Peter, G.; Bässler, H. Charge-Transfer Transitions in Crystalline Anthracene and Their Role in Photoconductivity. Chem. Phys. 1983, 75, 103−114. (64) Sharifzadeh, S.; Wong, C. Y.; Wu, H.; Cotts, B. L.; Kronik, L.; Ginsberg, N. S.; Neaton, J. B. Relating the Physical Structure and Optoelectronic Function of Crystalline Tips-Pentacene. Adv. Funct. Mater. 2015, 25, 2038−2046. (65) Lim, S.-H.; Bjorklund, T. G.; Spano, F. C.; Bardeen, C. J. Exciton Delocalization and Superradiance in Tetracene Thin Films and Nanoaggregates. Phys. Rev. Lett. 2004, 92, 107402. (66) Voigt, M.; Langner, A.; Schouwink, P.; Lupton, J. M.; Mahrt, R. F.; Sokolowski, M. Picosecond Time Resolved Photoluminescence Spectroscopy of a Tetracene Film on Highly Oriented Pyrolytic Graphite: Dynamical Relaxation, Trap Emission, and Superradiance. J. Chem. Phys. 2007, 127, 114705. (67) Camposeo, A.; Polo, M.; Tavazzi, S.; Silvestri, L.; Spearman, P.; Cingolani, R.; Pisignano, D. Polarized Superradiance from Delocalized Exciton Transitions in Tetracene Single Crystals. Phys. Rev. B: Condens. Matter Mater. Phys. 2010, 81, 033306. (68) Anger, F.; Osso, J. O.; Heinemeyer, U.; Broch, K.; Scholz, R.; Gerlach, A.; Schreiber, F. Photoluminescence Spectroscopy of Pure Pentacene, Perfluoropentacene, and Mixed Thin Films. J. Chem. Phys. 2012, 136, 054701. (69) Kistler, K. A.; Pochas, C. M.; Yamagata, H.; Matsika, S.; Spano, F. C. Absorption, Circular Dichroism, and Photoluminescence in Perylene Diimide Bichromophores: Polarization-Dependent H- and JAggregate Behavior. J. Phys. Chem. B 2012, 116, 77−86.

(39) Gisslen, L.; Scholz, R. Crystallochromy of Perylene Pigments: Influence of an Enlarged Polyaromatic Core Region. Phys. Rev. B: Condens. Matter Mater. Phys. 2011, 83, 155311. (40) Gisslen, L.; Scholz, R. Crystallochromy of Perylene Pigments: Interference between Frenkel Excitons and Charge-Transfer States. Phys. Rev. B: Condens. Matter Mater. Phys. 2009, 80, 115309. (41) Yamagata, H.; Maxwell, D. S.; Fan, J.; Kittilstved, K. R.; Briseno, A. L.; Barnes, M. D.; Spano, F. C. Hj-Aggregate Behavior of Crystalline 7,8,15,16-Tetraazaterrylene: Introducing a New Design Paradigm for Organic Materials. J. Phys. Chem. C 2014, 118, 28842−28854. (42) Hestand, N. J.; Tempelaar, R.; Knoester, J.; Jansen, T. L. C.; Spano, F. C. Exciton Mobility Control through Sub-Angstrom Packing Modifications in Molecular Crystals. Phys. Rev. B: Condens. Matter Mater. Phys. 2015, 91, 195315. (43) Schlosser, D. W.; Philpott, M. R. ″Singlet Excitons in Crystalline Naphthalene, Anthracene, Tetracene and Pentacene. Chem. Phys. 1980, 49, 181−199. (44) Philpott, M. R. Optical Reflection Spectroscopy of OrganicSolids. Annu. Rev. Phys. Chem. 1980, 31, 97−129. (45) Laudise, R. A.; Kloc, C.; Simpkins, P. G.; Siegrist, T. Physical Vapor Growth of Organic Semiconductors. J. Cryst. Growth 1998, 187, 449−454. (46) Holmes, D.; Kumaraswamy, S.; Matzger, A. J.; Vollhardt, K. P. C. On the Nature of Nonplanarity in the [N]Phenylenes. Chem. - Eur. J. 1999, 5, 3399−3412. (47) Philpott, M. R. Theory of Coupling of Electronic and Vibrational Excitations in Molecular Crystals and Helical Polymers. J. Chem. Phys. 1971, 55, 2039−2054. (48) Spano, F. C. The Spectral Signatures of Frenkel Polarons in Hand J-Aggregates. Acc. Chem. Res. 2010, 43, 429−439. (49) Spano, F. C. Excitons in Conjugated Oligomer Aggregates, Films, and Crystals. Annu. Rev. Phys. Chem. 2006, 57, 217−243. (50) Spano, F. C. Absorption and Emission in Oligo-Phenylene Vinylene Nanoaggregates: The Role of Disorder and Structural Defects. J. Chem. Phys. 2002, 116, 5877−5891. (51) Stradomska, A.; Petelenz, P. Intermediate Vibronic Coupling in Sexithiophene Single Crystals. Ii. Three-Particle Contributions. J. Chem. Phys. 2009, 131, 044507. (52) Spano, F. C.; Yamagata, H. Vibronic Coupling in J-Aggregates and Beyond: A Direct Means of Determining the Exciton Coherence Length from the Photoluminescence Spectrum. J. Phys. Chem. B 2011, 115, 5133−5143. (53) Spano, F. C. The Fundamental Photophysics of Conjugated Oligomer Herringbone Aggregates. J. Chem. Phys. 2003, 118, 981− 994. (54) Spano, F. C.; Meskers, S. C. J.; Hennebicq, E.; Beljonne, D. Probing Excitation Delocalization in Supramolecular Chiral Stacks by Means of Circularly Polarized Light: Experiment and Modeling. J. Am. Chem. Soc. 2007, 129, 7044−7054. (55) Spano, F. C.; Meskers, S. C. J.; Hennebicq, E.; Beljonne, D. Using Circularly Polarized Luminescence to Probe Exciton Coherence in Disordered Helical Aggregates. J. Chem. Phys. 2008, 129, 024704. (56) Tempelaar, R.; Stradomska, A.; Knoester, J.; Spano, F. C. Anatomy of an Exciton: Vibrational Distortion and Exciton Coherence in H- and J-Aggregates. J. Phys. Chem. B 2013, 117, 457−466. (57) Tsiper, E. V.; Soos, Z. G. Electronic Polarization in Pentacene Crystals and Thin Films. Phys. Rev. B: Condens. Matter Mater. Phys. 2003, 68, 085301. (58) Sebastian, L.; Weiser, G.; Bässler, H. Charge-Transfer Transitions in Solid Tetracene and Pentacene Studied by ElectroAbsorption. Chem. Phys. 1981, 61, 125−135. (59) Coropceanu, V.; Cornil, J.; da Silva Filho, D. A.; Olivier, Y.; Silbey, R.; Bredas, J. L. Charge Transport in Organic Semiconductors. Chem. Rev. 2007, 107, 926−952. (60) Bree, A.; Lyons, L. E. The Intensity of Ultraviolet-Light Absorption by Monocrystals 0.4. Absorption by Naphthacene of Plane-Polarized Light. J. Chem. Soc. 1960, 5206−5212. 22147

DOI: 10.1021/acs.jpcc.5b07163 J. Phys. Chem. C 2015, 119, 22137−22147