Coherent Dynamics of Mixed Frenkel and Charge-Transfer Excitons in

Mar 24, 2016 - ... of Economics and Technology, Sogutozu, Ankara 06560, Turkey. ¶ Department of Chemistry and Chemical Biology, Harvard University, C...
1 downloads 11 Views 552KB Size
Subscriber access provided by Chicago State University

Letter

Coherent Dynamics of Mixed Frenkel and Charge Transfer Excitons in Dinaphtho[2,3-b:2#3#-f]thieno[3,2-b]-thiophene Thin Films: The Importance of Hole Delocalization Takatoshi Fujita, Sule Atahan-Evrenk, Nicolas Per Dane Sawaya, and Alán Aspuru-Guzik J. Phys. Chem. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.jpclett.6b00364 • Publication Date (Web): 24 Mar 2016 Downloaded from http://pubs.acs.org on March 25, 2016

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry Letters is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Coherent Dynamics of Mixed Frenkel and Charge Transfer Excitons in Dinaphtho[2,3-b:2′3′-f]thieno[3,2-b]-thiophene Thin Films: The Importance of Hole Delocalization Takatoshi Fujita,∗,† Sule Atahan-Evrenk,‡ Nicolas P. D. Sawaya,¶ and Alán Aspuru-Guzik∗,¶ Department of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan, TOBB University of Economics and Technology, Sogutozu, Ankara 06560, Turkey, and Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA E-mail: [email protected]; [email protected]

∗ To

whom correspondence should be addressed of Chemistry, Graduate School of Science, Kyoto University, Kyoto 606-8502, Japan ‡ TOBB University of Economics and Technology, Sogutozu, Ankara 06560, Turkey ¶ Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA † Department

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

Abstract Charge transfer states in organic semiconductors play crucial roles in processes such as singlet fission and exciton dissociation at donor/acceptor interfaces. Recently, a time-resolved spectroscopy study of dinaphtho[2,3-b:2′ 3′ -f]thieno[3,2-b]-thiophene (DNTT) thin films provided evidence for the formation of mixed Frenkel and charge-transfer excitons after the photoexcitation. Here we investigate optical properties and excitation dynamics of the DNTT thin films by combining ab initio calculations and a stochastic Schrödinger equation. Our theory predicts that the low-energy Frenkel exciton band consists of 8 to 47% CT character. The quantum dynamics simulations show coherent dynamics of Frenkel and CT states in 50 fs after the optical excitation. We demonstrate the role of charge delocalization and localization in

electron/hole delocalization

the mixing of CT states with Frenkel excitons as well as the role of their decoherence.

Time (fs)

e-h separation (Å)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

keywords electron-hole dynamics, organic semiconductor, charge-transfer exciton, exciton delocalization, organic crystals, stochastic Schrödinger equation, fragment molecular orbital

2 ACS Paragon Plus Environment

Page 2 of 22

Page 3 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Organic semiconductors (OSCs) are widely investigated as candidates for inexpensive and flexible materials for photovoltaics and other optoelectronic applications. 1–3 In recent years, new OSCs have been designed and improved through computational modeling 4–6 and virtual high-throughput screening. 7–9 Modeling energy and charge transport properties is also essential for understanding structure-property relationships and thus for rationally designing novel OSCs. 10–16 The low-energy optical excitations in OSCs lead to the formation of a bound electron–hole (e–h) pair, a Frenkel exciton. Excitonic properties of organic crystals are substantially different from those of isolated molecules, owing to excitonic couplings, near-field interactions between electronic transitions. 17–20 Interactions with a charge-transfer (CT) state—a state in which the electron and hole are located on spatially separated regions—can also affect the optical properties of an exciton, as has been demonstrated by several theoretical studies. 21–27 Experimentally, the degree of CT character has been studied by momentum-dependent electron-loss spectroscopy 28–30 or electroabsorption spectroscopy. 31–33 CT states can act as precursors for interfacial CT states; 34–36 mixing of CT states with Frenkel excitons would facilitate charge separation in photovoltaic materials. They have also gained recent attention due to their relevance to singlet fission, 37–39 in which singlet to triplet conversion can proceed through CT states. 40–43 Here we focus on dinaphtho[2,3-b:2′ 3′ -f]thieno[3,2-b]-thiophene (DNTT), a p-type OSC originally introduced by Takimiya and co-workers. 44 DNTT and its derivatives 5,45–51 have gained attention due to their high hole mobility values and air stability. A recent time-resolved spectroscopy study by Ishino et al. 50 concludes that mixed Frenkel and CT excitons are formed after the optical excitation. Although the degree of CT character in excited states has been reported as described in the earlier paragraph, its role in excitation dynamics has remained unclear. Furthermore, there has been growing interest in the effects of excited-state delocalization on charge photogeneration processes. 36,52–57 Therefore, we studied the excited-state dynamics in DNTT as a model system to provide insight into the role of the CT states in OSCs. In this Letter, we present a theoretical study of optical properties and dissipative excitonic dynamics of DNTT. The optical absorption spectra and associated CT character are computed by

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 22

introducing a tight-binding Hamiltonian that treats Frenkel and CT states on the same footing. The excitation dynamics of the DNTT thin films was simulated by a stochastic Schrödinger equation with spectral densities derived from molecular dynamics simulations and excited state quantum chemistry calculations. Our simulations show coherent dynamics of the Frenkel and CT excitons in 50 fs after optical excitation. The importance of charge delocalization in the mixing of CT states will be discussed. We adapt a tight-binding model of a system comprising of one electron and one hole in a molecular aggregate. 19,20,58,59 Each molecule has two frontier orbitals: the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). A diabatic state of one electron and one hole, with localized wavefunctions, is obtained from creation operators as fol( ) ⟩ ⟩ ⟩ lows: ei h j = √1 c†iα d†jβ + c†iβ d†jα 0 , where 0 is an electronic ground state of the molecular 2

aggregate, the operators c†iα and d†jβ create an electron of α spin on ith molecule (site) and a hole ⟩ of β spin on jth molecule, respectively. ei hi represents a localized Frenkel state on ith molecule, ⟩ while ei h j denotes a CT state of an electron and a hole being localized on ith and jth molecules, respectively. The electronic Hamiltonian in this one-electron one-hole basis can be written as follows: ⟨

⟩ (f) (e) ei h j He ek hl = δ jl tik + δik t(h) + δkl δi j (1 − δik )Wik − δik δ jl Wi(c) j . jl

(1)

ti(e/h) are composed of electron/hole transfer integrals for off-diagonal elements and electron affinij (f)

ties/ionization potentials for diagonals. Wi j are excitonic and Wi(c) j are e–h Coulomb interactions. To calculate the electronic Hamiltonian, we follow Fujimoto 24,60 and adapt that methodology to the fragment molecular orbital methods. 61–63 Representative state energies and electronic couplings are shown in Table 1 and 2. The transfer integrals and excitonic interactions were calculated from fragment molecular orbital calculations at the Hartree-Fock/6-31G* level. The localized excitation energy, ionization potential, and electron affinity were obtained from density functional theory at the B3LYP/aug-cc-pvdz level. Further details of the electronic structure calculations are shown in the supporting information. The present calculations of the electron-hole Hamiltonian 4 ACS Paragon Plus Environment

Page 5 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

is based on the mean-field approximation. A more accurate parameterization would be possible by solving the Bethe-Salpeter equation which explicitly takes account of electron-hole correlation effects. 25–27 We define the energy threshold of 4.50 eV; The CT states higher than the threshold were omitted because they did not mix with the lower Frenkel states. By diagonalizing the electronic Hamilto ⟩ ⟩ ∑ ⟩ ⟩ nian, we get an adiabatic wave function, He ψI = E I ψI , where ψI = i, j Ci,I j ei h j . We obtained ∑ I µi , where the absorption spectrum from the corresponding transition dipole moments, µI = i Ci,i µi denotes a transition dipole moment of the ith molecule. We use a 7 × 7 × 7 supercell containing 686 molecules for the model structure for the DNTT crystal. Each energy level is broadened by a convolution of a Lorentz function of one meV. We also calculated CT populations of the Ith adia∑ ⟨ ⟩ 2 , which quantifies the probability that the electron and batic wave function, pCT I = i, j ei h j |ψI hole are found on different molecules. Figure 1(d) shows the computed optical absorption spectrum of the DNTT crystal in comparison with the measurement from ref. 50 The absorption spectrum is composed of the main excitation with the peak position of 2.83 eV, and relatively weak excitation at around 3.36 and 3.59 eV. The lower energy bands consist of Frenkel states partially mixed with the CT states. The latter two weak bands have mostly CT character, the intensities of which are borrowed from the lower Frenkel energy bands. The CT character of the brightest state is 22%, and that of largest Frenkel band is 47%. The red edge of the main excitation has lower CT character due to the energy difference with the CT states. Our model predicts that the low energy exciton band has 8 to 47% of CT character. The computed absorption spectrum is in reasonable agreement with the experiment, while the peak at around 2.94 eV cannot be reproduced. This may be ascribed to vibrational progressions. 50 However, it lies outside the scope of this letter to incorporate the vibrational degrees of freedom into the absorption spectra. To investigate the quantum dynamics of the Frenkel–CT excitons in the DNTT, we apply the stochastic Schrödinger equation formulated by Zhao and coworkers. 64–67 In this formalism, time evolution of an excitonic wave function coupled to a phonon bath is obtained as follows:

5 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

] ∑ ∑[ ∫ t iHe τ iHe τ ∂ (0) − Ln iℏ ψ(t)⟩ = He ψ(t)⟩ + Ln un (t) ψ(t)⟩ − i dτCn (τ)e ℏ Ln e ℏ ψ(t)⟩, ∂t 0 n n

Page 6 of 22

(2)

where Ln = n⟩⟨n . Here effects of the phonon bath are incorporated through a stochastic force un (t) and a zero-temperature bath correlation function Cn(0) (t). Both of them are calculated from a spectral density jn (ω) for the nth diabatic state. 64,67 The spectral densities 68,69 have been obtained by combing molecular dynamics simulations for the DNTT crystal with excited-state calculations using time-dependent density functional theory at B3LYP/6-31G* level. The details for calculating spectral densities are presented in the supporting information. An ensemble average of stochastic realization of eq (2) yields the two-body electron–hole density matrix,

⟨ ⟩ ρeh (t) = ψ(t)⟩⟨ψ(t)

ens

.

(3)

A system of 18 DNTT molecules in the ab plane was taken as a model of thin films. The absorption spectrum and excited states of this system are shown in the supporting information. As illustrative ⟩ ⟩ examples, localized Frenkel ( e1 h1 ) and CT ( e1 h2 ) states are used for initial conditions for the simulations of dissipative quantum dynamics. See Figure 2(a) and (b) for the structure and site labeling. The state energies of the Frenkel and CT states are 2.85 and 3.49 eV, respectively. The Runge-Kutta method was used for numerical propagation with a time of 0.1 fs. The excitonic density matrix and physical properties are averaged over 5,000 trajectories at the temperature of 300 K. To see dissipative dynamics of the 18 DNTT molecular system, Frenkel and CT populations, and electronic energies defined by ⟨ψ(t) He ψ(t)⟩ are shown in Figure 2(c) and (d). In the equilibrium state, the electronic energy, Frenkel and CT populations are respectively, 2.74 eV, 0.90, and 0.10, which are obtained as ensemble averages with respect to the Boltzmann distribution at 300K. Thermal equilibration is achieved within 1.5 and 2.7 ps from the initial Frenkel and CT states, respectively, as they are converged to the equilibrium electronic energy. From the initial CT state, 6 ACS Paragon Plus Environment

Page 7 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

the CT populations follow almost the same dynamics as the electronic energies; these plots reflect the relaxation dynamics from the higher CT state toward lower Frenkel states. On the other hand, from the initial Frenkel state, the populations show the fast oscillation and then converge toward the equilibrium Frenkel and CT populations. As shown in Figure 2(d), the CT states are strongly mixed with the Frenkel states in the initial 50 fs—a mixed Frenkel–CT exciton is formed. In what follows, we investigate the ultrafast time scale mixing of the CT states. In addition to the Frenkel and CT populations, e–h separation, inverse participation ratios (IPRs) of the electron and hole density matrices, and coherence with the initial Frenkel state are shown in Figure 3 (a)– ⟩ (d). The e–h separation of ei h j is defined as the center of mass distance between ith and jth molecules. To calculate the IPRs, one-body electron or hole density matrices were first calculated by tracing out the hole or electron basis: ρe (t) = T rh ρeh (t) or ρh (t) = T rh ρeh (t). The IPRs of electron or hole density matrix is given by )2 e/h ρ mn mn Lρe/h = . ∑ e/h 2 N mn ρmn (∑

(4)

These measures describe the delocalization of the electron or hole density matrix. 70,71 In Figure 3(a), the CT populations are increased and show oscillatory behavior which implies coherence among the Frenkel and CT states. The similar mixing effect between Frenkel and CT states was computationally observed for the pentacene crystal. 43 The oscillation period between the Frenkel and CT populations is about 7 fs, which corresponds to the 4761 cm−1 . The energy difference between the highest Frenkel-dominant state and the smallest CT-dominant state is 3746 cm−1 . The oscillation frequency would be determined by the energy differences among Frenkeldominant states and CT-dominant states. The corresponding picture are also observed from the e–h separation: The e–h separation is first increased and subsequently decreased to the equilibrated value of around 0.5 Å. Because the hole transfer integrals are higher than those for the electrons, there is an increase in the e–h separation which is subsequently moderated by e–h Coulomb attractions. This competition between the higher hole transfer and e–h Coulomb attraction results in the 7 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

oscillation. The effect of higher hole transfer is also seen by comparing electron and hole IPRs as represented in Figure 3(c), indicating that the hole IPR is greater than that of electron. The delocal ⟩ ization of the hole wave function is also confirmed by coherence among diabatic states, e1 hi , as presented in Figure 3(d). This delocalization can effectively weaken the e–h Coulomb attraction and thus become an initial driving force for increasing e–h separation and the mixing of CT states. Localization of wave functions due to the nuclear vibrations works toward decreasing the e–h separations. Our finding is in line with the observations by Tamura and Burghardt 54,55 that charge separation at the donor/acceptor interface can be enhanced by charge delocalization. In the time-resolved spectroscopy experiment by Ishino et al., 50 the pump pulses excite to exciton manifolds at 3.1 eV, which leads to the optical excitation of the blue edge of the absorption maximum. The subsequent relaxation to the lowest exciton state occurs within about 2.1 ps. Our simulations shows that the relaxation time from the localized Frenkel state, which is resonant to the absorption maximum, to the lowest exciton state is about 1.5 ps. The relaxation time of 1.5 ps would qualitatively correspond to that observed by the experiments. However, simulating more extended systems would be necessary for quantitative comparison between the simulation and experiments. The mixing of CT states in optical excitations of the DNTT has been proposed based on a derivative-like feature of transient absorption spectra: 50 Upon the excitation pulse of 3.1 eV the transition absorption spectra exhibit derivative-like features which are similar to Stark spectra. This result suggests the existence of transient charged species exerting electric fields on surrounding molecules. It follows that the CT excitons are formed after photoexcitation, the local electric fields of which induce Stark shifts of surrounding molecules. Relating our simulations with the experiments, we argue that the derivative-like features would appear after Frenkel–CT decoherence. If the optically-excited state is a superposition of Frenkel and CT states, the absorption of them cannot be disentangled. In other words, Stark effects from CT states occur after the Frenkel-CT decoherence time. Our simulations predict that the derivative-like feature may appear about 50

8 ACS Paragon Plus Environment

Page 8 of 22

Page 9 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

fs after the optical excitation. A more rigorous comparison to the experimental observations has to be done by directly calculating the time-resolved spectroscopy signals, 72 which will be studied elsewhere. We have presented a theoretical study on the optical properties and excitation dynamics of the DNTT thin films. The tight-binding Hamiltonian combined with ab initio calculations reasonably reproduces the experimental absorption spectrum and predicts that the low-energy Frenkel exciton band consist of 8 to 47% CT character. The quantum dynamics simulations show the coherent dynamics of Frenkel and CT states about 50 fs after photoexcitation. The oscillation and lifetime of the e–h separations are also analyzed in terms of charge delocalization/localization and the e–h Coulomb attractions. Combined with large-scale quantum dynamics simulations, 73,74 the present approach can be applied to more extended systems such as organic bulk heterojunction solar cells. Such a device-level simulation will be useful to design and improve novel organic electronic materials.

Acknowledgement The authors thank Dr Joonsuk Huh for helpful discussion and implementation of the stochastic Schrödinger equation and Thomas Markovich for the calculation of spectral densities. T.F thanks Prof. Yuji Mochizuki for discussion on the electronic-structure calculations and Prof. Kazuya Watanabe and Shunsuke Tanaka for useful discussions on their experiments. The majority of our computation was carried out on Harvard University’s Odyssey cluster, supported by the Research Computing Group of the FAS Division of Science. T.F acknowledge support from JSPS Grant-inAid for Scientific Research on Innovative Areas "Dynamical ordering of biomolecular systems for creation of integrated functions" (Grant No. 25102002). T. F., N. S. and A. A.-G. acknowledge support from the Center for Excitonics, an Energy Frontier Research Center funded by the US Department of Energy, Office of Science and Office of Basic Energy Sciences under award DESC0001088.

9 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Supporting Information Available The computational details for calculations of electronic Hamiltonian, spectral densities, absorption spectrum and excited states of the 18 DNTT molecular system. This material is available free of charge via the Internet at http://pubs.acs.org/.

References (1) Bendikov, M.; Wudl, F.; Perepichka, D. F. Tetrathiafulvalenes, Oligoacenenes, and Their Buckminsterfullerene Derivatives: The Brick and Mortar of Organic Electronics. Chem. Rev. 2004, 104, 4891–4946. (2) Anthony, J. E. Functionalized Acenes and Heteroacenes for Organic Electronics. Chem. Rev. 2006, 106, 5028–5048. (3) Brédas, J.-L.; Norton, J. E.; Cornil, J.; Coropceanu, V. Molecular Understanding of Organic Solar Cells: The Challenges. Acc. Chem. Res. 2009, 42, 1691–1699. (4) Blouin, N.; Michaud, A.; Gendron, D.; Wakim, S.; Blair, E.; Neagu-Plesu, R.; Belletête, M.; Durocher, G.; Tao, Y.; Leclerc, M. Toward a Rational Design of Poly (2,7-Carbazole) Derivatives for Solar Cells. J. Am. Chem. Soc. 2008, 130, 732–742. (5) Sokolov, A. N.; Atahan-Evrenk, S.; Mondal, R.; Akkerman, H. B.; Sánchez-Carrera, R. S.; Granados-Focil, S.; Schrier, J.; Mannsfeld, S. C. B.; Zoombelt, A. P.; Bao, Z.; et al. From Computational Discovery to Experimental Characterization of a High Hole Mobility Organic Crystal. Nat. Commun. 2011, 2, 437. (6) Li, J.; Zhao, Y.; Tan, H. S.; Guo, Y.; Di, C.-A.; Yu, G.; Liu, Y.; Lin, M.; Lim, S. H.; Zhou, Y.; Su, H.; Ong, B. S. A Stable Solution-Processed Polymer Semiconductor with Record HighMobility for Printed Transistors. Sci. Rep. 2012, 2, 1–9. (7) Hachmann, J.; Olivares-Amaya, R.; Jinich, A.; Appleton, A. L.; Blood-Forsythe, M. A.; Seress, L. R.; Roman-Salgado, C.; Trepte, K.; Atahan-Evrenk, S.; Er, S.; et al. Lead Candidates 10 ACS Paragon Plus Environment

Page 10 of 22

Page 11 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

for High-Performance Organic Photovoltaics from High-Throughput Quantum Chemistry the Harvard Clean Energy Project. Energy Environ. Sci. 2014, 7, 698–704. (8) O’Boyle, N. M.; Campbell, C. M.; Hutchison, G. R. Computational Design and Selection of Optimal Organic Photovoltaic Materials. J. Phys. Chem. C 2011, 115, 16200–16210. (9) Pyzer-Knapp, E. O.; Suh, C.; Gómez-Bombarelli, R.; Aguilera-Iparraguirre, J.; AspuruGuzik, A. What is High-Throughput Virtual Screening? A Perspective from Organic Materials Discovery. Annu. Rev. Mater. Sci. 2015, 45, 195–216. (10) Brédas, J.-L.; Beljonne, D.; Coropceanu, V.; Cornil, J. Charge-Transfer and Energy-Transfer Processes in π-Conjugated Oligomers and Polymers: A Molecular Picture. Chem. Rev. 2004, 104, 4971–5004. (11) Wang, L.; Nan, G.; Yang, X.; Peng, Q.; Li, Q.; Shuai, Z. Computational Methods for Design of Organic Materials with High Charge Mobility. Chem. Soc. Rev. 2010, 39, 423–434. (12) Troisi, A. Charge Transport in High Mobility Molecular Semiconductors: Classical Models and New Theories. Chem. Soc. Rev. 2011, 40, 2347–2358. (13) Baumeier, B.; May, F.; Lennartz, C.; Andrienko, D. Challenges for In Silico Design of Organic Semiconductors. J. Mater. Chem. 2012, 22, 10971–10976. (14) Mei, J.; Diao, Y.; Appleton, A. L.; Fang, L.; Bao, Z. Integrated Materials Design of Organic Semiconductors for Field-Effect Transistors. J. Am. Chem. Soc. 2013, 135, 6724–6746. (15) Atahan-Evrenk, S.; Aspuru-Guzik, A. Prediction and Theoretical Characterization of p-Type Organic Semiconductor Crystals for Field-Effect Transistor Applications. In Prediction and Calculation of Crystal Structures; Atahan-Evrenk, S., Aspuru-Guzik, A., Eds.; Topics in Current Chemistry; Springer International Publishing, 2014; Vol. 345; pp 95–138.

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(16) Jiang, Y.; Zhong, X.; Shi, W.; Peng, Q.; Geng, H.; Zhao, Y.; Shuai, Z. Nuclear Quantum Tunnelling and Carrier Delocalization Effects to Bridge the Gap between Hopping and Bandlike Behaviors in Organic Semiconductors. Nanoscale Horiz. 2016, 1, 53–59. (17) Scholes, G. D.; Rumbles, G. Excitons in Nanoscale Systems. Nature Mater. 2006, 5, 683– 696. (18) Spano, F. C. The Spectral Signatures of Frenkel Polarons in H- And J-aggregates. Acc. Chem. Res. 2010, 43, 429–439. (19) Kühn, O.; Lochbrunner, S. Quantum Dynamics and Spectroscopy of Excitons in Molecular Aggregates. In Quantum Efficiency in Complex Systems, Part II From Molecular Aggregates to Organic Solar Cells; Uli Wüerfel, M. T., Weber, E. R., Eds.; Semiconductors and Semimetals; Elsevier, 2011; Vol. 85; p 47. (20) Saikin, S. K.; Eisfeld, A.; Valleau, S.; Aspuru-Guzik, A. Photonics Meets Excitonics: Natural and Artificial Molecular Aggregates. Nanophotonics 2013, 2, 21–38. (21) Hoffmann, M.; Soos, Z. G. Optical Absorption Spectra of the Holstein Molecular Crystal for Weak and Intermediate Electronic Coupling. Phys. Rev. B 2002, 66, 024305. (22) Gisslén, L.; Scholz, R. Crystallochromy of Perylene Pigments: Interference between Frenkel Excitons and Charge-Transfer States. Phys. Rev. B 2009, 80, 115309. (23) Yamagata, H.; Norton, J.; Hontz, E.; Olivier, Y.; Beljonne, D.; Brédas, J. L.; Silbey, R. J.; Spano, F. C. The Nature of Singlet Excitons in Oligoacene Molecular Crystals. J. Chem. Phys. 2011, 134, 204703. (24) Fujimoto, K. J.; Kitamura, C. A Theoretical Study of Crystallochromy: Spectral Tuning of Solid-State Tetracenes. J. Chem. Phys. 2013, 139, 084511.

12 ACS Paragon Plus Environment

Page 12 of 22

Page 13 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(25) Sharifzadeh, S.; Darancet, P.; Kronik, L.; Neaton, J. B. Low-Energy Charge-Transfer Excitons in Organic Solids from First-Principles: The Case of Pentacene. J. Phys. Chem. Lett. 2013, 4, 2197–2201. (26) Cudazzo, P.; Gatti, M.; Rubio, A.; Sottile, F. Frenkel versus Charge-Transfer Exciton Dispersion in Molecular Crystals. Phys. Rev. B 2013, 88, 195152. (27) 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. (28) 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. (29) Roth, F.; Schuster, R.; König, A.; Knupfer, M.; Berger, H. Momentum Dependence of the Excitons in Pentacene. J. Chem. Phys. 2012, 136, 204708. (30) Roth, F.; Mahns, B.; Hampel, S.; Nohr, M.; Berger, H.; Büchner, B.; Knupfer, M. Exciton Properties of Selected Aromatic Hydrocarbon Systems. Eur. Phys. J. B 2013, 86, 66. (31) Sebastian, L.; Weiser, G.; Bässler, H. Charge Transfer Transitions in Solid Tetracene and Pentacene Studied by Electroabsorption. Chem. Phys. 1981, 61, 125 – 135. (32) Möller, S.; Weiser, G.; Garnier, F. Electric Field Effects on the Davydov Components of a Strong Intramolecular Transition: α-Sexithiophene Single Crystals. Phys. Rev. B 2000, 61, 15749–15755. (33) Haas, S.; Matsui, H.; Hasegawa, T. Field-Modulation Spectroscopy of Pentacene Thin Films Using Field-Effect Devices: Reconsideration of the Excitonic Structure. Phys. Rev. B 2010, 82, 161301. (34) Zhu, X.-Y.; Yang, Q.; Muntwiler, M. Charge-Transfer Excitons at Organic Semiconductor Surfaces and Interfaces. Acc. Chem. Res. 2009, 42, 1779–1787. 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(35) Clarke, T. M.; Durrant, J. R. Charge Photogeneration in Organic Solar Cells. Chem. Rev. 2010, 110, 6736–6767. (36) Bakulin, A.; Rao, A.; Pavelyev, V. G.; van Loosdrecht, P. H. M.; Pshenichnikov, M. S.; Niedzialek, D.; Cornil, J.; Beljonne, D.; Friend, R. H. The Role of Driving Energy and Delocalized States for Charge Separation in Organic Semiconductors. Science 2012, 335, 1340–1344. (37) Smith, M. B.; Michl, J. Singlet Fission. Chem. Rev. 2010, 110, 6891–6936. (38) Chan, W.-L.; Berkelbach, T. C.; Provorse, M. R.; Monahan, N. R.; Tritsch, J. R.; Hybertsen, M. S.; Reichman, D. R.; Gao, J.; Zhu, X.-Y. The Quantum Coherent Mechanism for Singlet Fission: Experiment and Theory. Acc. Chem. Res. 2013, 46, 1321–1329. (39) Monahan, N.; Zhu, X.-Y. Charge Transfer-Mediated Singlet Fission. Annu. Rev. Phys. Chem. 2015, 66, 601–618. (40) Beljonne, D.; Yamagata, H.; Brédas, 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. (41) Berkelbach, T. C.; Hybertsen, M. S.; Reichman, D. R. Microscopic Theory of Singlet Exciton Fission. I. General Formulation. J. Chem. Phys. 2013, 138, 114102. (42) Berkelbach, T. C.; Hybertsen, M. S.; Reichman, D. R. Microscopic Theory of Singlet Exciton Fission. II. Application to Pentacene Dimers and the Role of Superexchange. J. Chem. Phys. 2013, 138, 114103. (43) Berkelbach, T. C.; Hybertsen, M. S.; Reichman, D. R. Microscopic Theory of Singlet Exciton Fission. III. Crystalline Pentacene. J. Chem. Phys. 2014, 141, 074705. (44) Yamamoto, T.; Takimiya, K. Facile Synthesis of Highly π-Extended Heteroarenes,

14 ACS Paragon Plus Environment

Page 14 of 22

Page 15 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Dinaphtho[2,3-b:2′ 3′ -f]chalcogenopheno[3,2-b]chalcogenophenes, and Their Application to Field-Effect Transistors. J. Am. Chem. Soc. 2007, 129, 2224–2225. (45) Sánchez-Carrera, R. S.; Atahan, S.; Schrier, J.; Aspuru-Guzik, A. Theoretical Characterization of the Air-Stable, High-Mobility Dinaphtho[2,3-b:2′ 3′ -f]thieno[3,2-b]-thiophene Organic Semiconductor. J. Phys. Chem. C 2010, 114, 2334–2340. (46) Kang, M. J.; Doi, I.; Miyazaki, E.; Takimiya, K.; Ikeda, M.; Kuwabara, H. Alkylated Dinaphtho[2,3-b:2′ 3′ -f ]thieno[3,2- b]thiophenes (Cn -DNTTs): Organic Semiconductors for High-Performance Thin Film Transistors. Adv. Mater. 2011, 23, 1222–1225. (47) Kang, M. J.; Miyazaki, E.; Osaka, I.; Takimiya, K. Two Isomeric Didecyl-dinaphtho[2,3b:2′ 3′ -f]thieno[3,2-b]thiophenes: Impact of Alkylation Positions on Packing Structures and Organic Field Effect Transistor Characteristics. Jpn. J. Appl. Phys. 2012, 51, 11PD04. (48) Xie, W.; Willa, K.; Wu, Y.; Häusermann, R.; Takimiya, K.; Batlogg, B.; Frisbie, C. D. Temperature-Independent Transport in High-Mobility Dinaphtho-thieno-thiophene (DNTT) Single Crystal Transistors. Adv. Mater. 2013, 25, 3478–3484. (49) Kinoshita, Y.; Tanaka, H.; Shimoi, Y.; Takimiya, K.; Kuroda, S.-i. Low-Temperature Carrier Dynamics in High-Mobility Organic Transistors of Alkylated Dinaphtho-thienothiophene as Investigated by Electron Spin Resonance. Appl. Phys. Lett. 2014, 105, 033301. (50) Ishino, Y.; Miyata, K.; Sugimoto, T.; Watanabe, K.; Matsumoto, Y.; Uemura, T.; Takeya, J. Ultrafast Exciton Dynamics in Dinaphtho[2,3-b:2′ 3′ -f]thieno[3,2-b]-thiophene Thin Films. Phys. Chem. Chem. Phys. 2014, 16, 7501–7512. (51) Miyata, K.; Tanaka, S.; Ishino, Y.; Watanabe, K.; Uemura, T.; Takeya, J.; Sugimoto, T.; Matsumoto, Y. Microscopic Hole-Transfer Efficiency in Organic Thin-Film Transistors Studied with Charge-Modulation Spectroscopy. Phys. Rev. B 2015, 91, 195306.

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(52) Kaake, L. G.; Moses, D.; Heeger, A. J. Coherence and Uncertainty in Nanostructured Organic Photovoltaics. J. Phys. Chem. Lett. 2013, 4, 2264–2268. (53) Chen, K.; Barker, A. J.; Reish, M. E.; Gordon, K. C.; Hodgkiss, J. M. Broadband Ultrafast Photoluminescence Spectroscopy Resolves Charge Photogeneration via Delocalized Hot Excitons in Polymer:Fullerene Photovoltaic Blends. J. Am. Chem. Soc. 2013, 135, 18502– 18512. (54) Tamura, H.; Burghardt, I. Potential Barrier and Excess Energy for Electron-Hole Separation from the Charge-Transfer Exciton at Donor-Acceptor Heterojunctions of Organic Solar Cells. J. Phys. Chem. C 2013, 117, 15020–15025. (55) Tamura, H.; Burghardt, I. Ultrafast Charge Separation in Organic Photovoltaics Enhanced by Charge Delocalization and Vibronically Hot Exciton Dissociation. J. Am. Chem. Soc. 2013, 135, 16364–16367. (56) Kaake, L. G.; Zhong, C.; Love, J. A.; Nagao, I.; Bazan, G. C.; Nguyen, T.-Q.; Huang, F.; Cao, Y.; Moses, D.; Heeger, A. J. Ultrafast Charge Generation in an Organic Bilayer Film. J. Phys. Chem. Lett. 2014, 5, 2000–2006. (57) Kaake, L. G.; Moses, D.; Heeger, A. J. Charge Transfer from Delocalized Excited States in a Bulk Heterojunction Material. Phys. Rev. B 2015, 91, 075436. (58) Meier, T.; Thomas, P.; Koch, S. W. Coherent Semiconductor Optics: From Basic Concepts to Nanostructure Applications; Springer-Verlag Berlin Heidelberg, 2007. (59) Abramavicius, D.; Mukamel, S. Energy-Transfer and Charge-Separation Pathways in the Reaction Center of Photosystem II Revealed by Coherent Two-Dimensional Optical Spectroscopy. J. Chem. Phys. 2010, 133, 184501. (60) Fujimoto, K. J. Transition-Density-Fragment Interaction Combined with Transfer Integral

16 ACS Paragon Plus Environment

Page 16 of 22

Page 17 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Approach for Excitation-Energy Transfer via Charge-Transfer States. J. Chem. Phys. 2012, 137, 034101. (61) Mochizuki, Y.; Koikegami, S.; Amari, S.; Segawa, K.; Kitaura, K.; Nakano, T. Configuration Interaction Singles Method with Multilayer Fragment Molecular Orbital Scheme. Chem. Phys. Lett. 2005, 406, 283–288. (62) Fedorov, D. G., Kitaura, K., Eds. The Fragment Molecular Orbital Method: Practical Applications to Large Molecular Systems; CRC Press: Boca Raton, FL, 2009. (63) Tanaka, S.; Mochizuki, Y.; Komeiji, Y.; Okiyama, Y.; Fukuzawa, K. Electron-Correlated Fragment-Molecular-Orbital Calculations for Biomolecular and Nano Systems. Phys. Chem. Chem. Phys. 2014, 16, 10310–10344. (64) Zhong, X.; Zhao, Y. Non-Markovian Stochastic Schrödinger Equation at Finite Temperatures for Charge Carrier Dynamics in Organic Crystals. J. Chem. Phys. 2013, 138, 014111. (65) Han, L.; Zhong, X.; Liang, W.; Zhao, Y. Energy Relaxation and Separation of a Hot ElectronHole Pair in Organic Aggregates from a Time-Dependent Wavepacket Diffusion Method. J. Chem. Phys. 2014, 140, 214107. (66) Han, L.; Ke, Y.; Zhong, X.; Zhao, Y. Time-Dependent Wavepacket Diffusion Method and Its Applications in Organic Semiconductors. Int. J. Quant. Chem. 2015, 115, 578–588. (67) Abramavicius, V.; Abramavicius, D. Excitation Transfer Pathways in Excitonic Aggregates Revealed by the Stochastic Schrödinger Equation. J. Chem. Phys. 2014, 140, 065103. (68) Damjanovíc, A.; Kosztin, I.; Kleinekathöfer, U.; Schulten, K. Excitons in a Photosynthetic Light-Harvesting System: A Combined Molecular Dynamics, Quantum Chemistry and Polaron Model Study. Phys. Rev. E 2002, 65, 031919. (69) Valleau, S.; Eisfeld, A.; Aspuru-Guzik, A. On the Alternatives for Bath Correlators and

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Spectral Densities from Mixed Quantum-Classical Simulations. J. Chem. Phys. 2012, 137, 224103. (70) Mukamel, S. Multidimensional Femtosecond Correlation Spectroscopies of Electronic and Vibrational Excitations. Annu. Rev. Phys. Chem. 2000, 51, 691–729. (71) Dahlbom, M.; Pullerits, T.; Mukamel, S.; Sundström, V. Exciton Delocalization in the B850 Light-Harvesting Complex: Comparison of Different Measures. J. Phys. Chem. B 2001, 105, 5515–5524. (72) Mukamel, S. Principles of Nonlinear Optical Spectroscopy; Oxford University Press: Oxford, U.K., 1999. (73) Huh, J.; Saikin, S. K.; Brookes, J. C.; Valleau, S.; Fujita, T.; Aspuru-Guzik, A. Atomistic Study of Energy Funneling in the Light-Harvesting Complex of Green Sulfur Bacteria. J. Am. Chem. Soc. 2014, 136, 2048–2057. (74) Sawaya, N. P. D.; Huh, J.; Fujita, T.; Saikin, S. K.; Aspuru-Guzik, A. Fast Delocalization Leads To Robust Long-Range Excitonic Transfer in a Large Quantum Chlorosome Model. Nano Lett. 2015, 15, 1722–1729.

18 ACS Paragon Plus Environment

Page 18 of 22

Page 19 of 22

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

Table 1: Electron and hole transfer integrals, and excitonic interactions in meV for molecular pairs in the DNTT crystal structure in the ab plane (see Figure 1(b) for pair labeling).

electron

hole

exciton

18 65 49

111 141 48

2 41 22

Pair 1 Pair 2 Pair 3

Table 2: Frenkel and CT energies in eV (see Figure 1(b) for site labeling).

State

e1 h1

e1 h2

e1 h3

e2 h1

e2 h3

e3 h1

e3 h2

Energy

2.85

3.49

3.80

3.53

3.48

3.77

3.51

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

(a)

(d)

(b) b

Pair 1 Site 1 Site 3 Pair 2 Pair 3

(e)

a

Site 2

(c) b

c

Figure 1: (a) Chemical structure of DNTT. Crystal structure of DNTT in (b) ab and (c) bc planes. The crystal parameters are as follows: a = 6.187, b = 7.662, c=16.62 Å; β=92.49◦ . 44 (d) Computed absorption spectrum of DNTT crystal in comparison with the experiment. 50 (e) The CT populations of excited states.

20 ACS Paragon Plus Environment

Page 20 of 22

Page 21 of 22

(c)

(b)

11 3 13 12 4 14 5 1 7 6 2 8 15 9 17 16 10 18

F/CT populations

(a)

Electronic Energy (eV)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry Letters

(d)

Time (fs) Figure 2: (a) A structure and (b) site indices of the 18 model. Relaxation dynamics from DNTT ⟩ ⟩ (c) CT (ψ(t = 0) = e1 h2 ) and (d) Frenkel (ψ(t = 0) = e1 h1 ) states. In (c) and (d), Frenkel and CT populations and electronic energies are shown.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry Letters

1.0 Electron Hole

(c)

4.0

0.8 3.0 0.6 F CT

IPR

Frenkel/CT characters

(a)

0.4

2.0 1.0

0.2 0.0

0.0 0

20

40

60

80

100

0

20

40

2.5 (b) 2.0 1.5 1.0 0.5 0.0

60

80

100

40

50

e1h2 e1h4 e1h5 e1h7 e1h8 e1h14

(d) Coherence with e1h1

e−h separation (Å)

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 22

0.2

0.1

0.0 0

20

40

60

80

100

0

10

Time(fs)

20

30

⟩ Figure 3: Ultrafast dynamics from Frenkel state ( e1 h1 ): (a) Frenkel and CT populations. (b) Electron–hole separation (c) inverse participation ratio of electron and hole density matrix, and (d) absolute values of coherence with the initial Frenkel state. The time range of (d) is different from others.

22 ACS Paragon Plus Environment