Improving the Stability of Organic Semiconductors: Distortion Energy

Nov 15, 2016 - Kingdom of Saudi Arabia. ‡. Department of Polymer Science & Engineering, Conte Polymer Research Center, University of Massachusetts, ...
0 downloads 0 Views 3MB Size
Article pubs.acs.org/cm

Improving the Stability of Organic Semiconductors: Distortion Energy versus Aromaticity in Substituted Bistetracene Simil Thomas,† Jack Ly,‡ Lei Zhang,‡ Alejandro L. Briseno,‡ and Jean-Luc Bredas*,† †

Laboratory for Computational and Theoretical Chemistry of Advanced Materials, Solar and Photovoltaics Engineering Research Center, Division of Physical Science and Engineering, King Abdullah University of Science and Technology, Thuwal 23955-6900, Kingdom of Saudi Arabia ‡ Department of Polymer Science & Engineering, Conte Polymer Research Center, University of Massachusetts, Amherst, Massachusetts 01003, United States S Supporting Information *

ABSTRACT: Polycyclic aromatic hydrocarbons (PAHs) have been widely explored as molecular semiconductors in organic electronic devices such as field-effect transistors or solar cells. However, their tendency to undergo photooxidation is a primary limitation to their practical applications. Bistetracene derivatives have recently been demonstrated to possess much larger photooxidation stability than the widely investigated pentacene and rubrene, while maintaining high charge-carrier mobilities. Here, using several levels of density functional theory, we identify the origin of the increased stability of bistetracene with respect to molecular oxygen by systematically investigating the [4 + 2] cycloaddition (Diels−Alder) photooxidation reaction mechanism. Importantly, our computational results indicate that endoperoxide formation in bis(2-(trimethylsilyl)ethynyl) bistetracene (BT) occurs not on the ring with least aromaticity, but rather on the ring with smallest distortion energy. This feature was subsequently confirmed by experimental NMR analyses. The oxidation activation barriers of bistetracene, pentacene, and rubrene are found to be 17.7, 13.6, and 14.4 kcal/mol, respectively, in agreement with the observed order of stability of these molecules with respect to oxidation reactions in solution. In the cases of BT and pentacene, the rates of electron transfer to create charged species (PAH+ and O2−) are at least two orders of magnitude lower than that of the charge recombination process (back to PAH and O2); for rubrene, both of these processes are calculated to be of the same order of magnitude, in agreement with experimental electron paramagnetic resonance spectroscopy observations.

1. INTRODUCTION Molecular semiconductors such as pentacene and rubrene have garnered considerable interest as model systems since they serve as active materials in organic photovoltaic (OPV) devices,1 organic field-effect transistors (OFETs),2−5 Hall effect devices, etc.6,7 Moreover, hole mobilities greater than 5 cm2/(V s)2,3 in pentacene and 15 cm2/(V s)4,5,8 in rubrene have been reported in transistor-type configurations; these values are significantly larger than that in hydrogenated amorphous silicon transistors (1 cm2/(V s)), a reference material for thin-film transistors.9,10 However, when exposed to light and oxygen in solution, pentacene and rubrene quickly degrade via a photooxidation reaction to produce transannular endoperoxides within a few minutes.11,12 The increased reactivity of the central ring of polyacenes toward Diels−Alder reactions with molecular oxygen in the presence of light has been attributed to a reduction in local aromaticity of the individual rings toward the central ring.13,14 Since the highest occupied molecular orbital (HOMO) is stabilized by the oxidation of pentacene (due to reduced conjugation length), the generated endoperoxide © 2016 American Chemical Society

molecule cannot easily exchange a hole with an adjacent pentacene, which represents a barrier for hole transport. Also, the endoperoxide structure is nonplanar and creates a local lattice deformation that acts as a scattering center. Hence, charge mobility in the device decreases with an increase in the number of oxidized species.15,16 Similar features have been observed in rubrene.12 Understanding the origin of the factors preventing oxidative degradation can thus represent a major step toward the rational design of chemically stable organic semiconductors.12,17−20 Polycyclic aromatic hydrocarbons (PAHs) react photochemically with O2 via two mechanistic pathways: (1) In the electron transfer mechanism (referred to as type I), a photoexcited PAH molecule transfers an electron to a neighboring O2 to create a PAH cation (PAH+) and an oxygen anion (O 2− ); these two charged species subsequently react to form an endoperoxide.21 Received: June 23, 2016 Revised: October 16, 2016 Published: November 15, 2016 8504

DOI: 10.1021/acs.chemmater.6b02552 Chem. Mater. 2016, 28, 8504−8512

Article

Chemistry of Materials (2) In the energy transfer mechanism (referred to as type II), excited PAHs undergo intersystem crossing (ISC) from an excited singlet state (S1) to the triplet state (T1), which has a longer lifetime. If the triplet energy of the PAH is greater than that of the singlet−triplet gap (ΔST) of O2, the PAH molecule can transfer its energy to the 3 O2 molecule, which produces an excited 1O2 molecule.21 Singlet oxygen (1O2) is highly reactive with a lifetime of 45 min in vacuum, 4.4 μs in water,22 and 240 μs in chloroform.23 Generated 1O2 reacts with PAH to form an endoperoxide through a concerted or a biradical stepwise mechanism. The type-I and type-II mechanisms are illustrated in Figure 1. Earlier theoretical studies, based on the G3(MP2) computa-

Figure 2. Chemical structures of: (A) TIPS-BT, (B) pentacene, and (C) rubrene, along with the ring numbering considered throughout this work.

2. MATERIALS AND METHODS 2.1. Computational Methods. The DFT calculations were performed using the Gaussian 09 Revision D.01 suite of programs29 and generally the 6-31G(d) basis.30,31 We have initially considered the ωB97X-D,32 B3LYP,33 CAM-B3LYP,34 and M06-2X35 functionals. Restricted and unrestricted DFT methodologies were used for closedshell and open-shell species, respectively. The B3LYP and CAMB3LYP functionals were found to be inadequate to locate the concerted transition state during endoperoxide formation of pentacene, rubrene, and TIPS-BT, whereas those transition states are well located with the ωB97X-D and M06-2X functionals. Frequency calculations were carried out: (i) to confirm that the geometries correspond to local minima (no imaginary frequency) or transition states (one imaginary frequency pertaining to the reaction coordinates, which was confirmed in addition by intrinsic-reactioncoordinate (IRC) calculations); and (ii) to provide the thermochemical data at standard temperature (298.15 K) and pressure (1 atm), which include entropy contributions. The energies were further refined by considering the 6-311+G(d,p) basis. The influence of solvent was modeled by using the integral equation formalism of the polarizable continuum model (IEF-PCM), a solvation model within the selfconsistent reaction field (SCRF) framework (taking into account a dielectric constant, ε = 4.71, equivalent to chloroform, to maintain consistency with experiment).36 We have found that optimizing the geometry within the IEF-PCM has only a marginal effect on geometrical parameters; thus, gas-phase optimized geometries are reported throughout this work. The energetics of the reactions are presented on the basis of relative Gibbs energies (ΔGrel) with respect to infinitely separated reactants. The ΔST energy difference between the triplet ground state and singlet excited state of O2, as calculated via the ωB97X-D functional with a spin projection technique,37 is 23.5 kcal/mol (1.02 eV), in good agreement with the experimental value of 22.5 kcal/mol (0.98 eV).38 When calculated using the M06-2X functional, the ΔST of the O2 molecule is significantly overestimated (28.5 kcal/mol), which may result in an underestimation of the barrier heights of the reaction. Therefore, the properties calculated using the ωB97X-D functional with spin contamination removal are considered throughout this study unless otherwise specified. The distortion/interaction model,39−41 known as the activation strain model,42,43 has been investigated to understand the origin of regioselectivity of oxidation on various rings of TIPS-BT. In this model, the activation energy (ΔEact) is decomposed into the distortion energy (ΔEdist) of each fragment and the interaction energy between the fragments (ΔEint):

Figure 1. Electron transfer (type I) and energy transfer (type II) pathways for the photooxidation reaction of PAHs with molecular oxygen along with the general chemical structure of the formed endoperoxide of triisopropylsilylethynyl-substituted (TIPS) bistetracene (TIPS-BT).

tional method, reported a concerted mechanism for 1,4addition reactions between 1O2 and PAHs such as benzene, naphthalene, anthracene, tetracene, and pentacene.24 Subsequent studies on polyacenes, using density functional theory (DFT) calculations at the B3LYP/6-31G* level, showed that the smaller analogs, such as benzene and naphthalene, undergo a concerted reaction pathway, whereas the larger analogs undergo a biradical stepwise mechanism.25 However, since B3LYP fails to account for dispersion interaction effects, the concerted transition state cannot be properly located when using such a functional. Therefore, here, we have chosen to systematically investigate the PAH degradation mechanism in more detail at several levels of DFT including dispersioncorrected functionals. Triisopropylsilylethynyl (TIPS)-substituted pentacene has been reported to be more stable than pentacene itself.26,27 This property has been attributed to a reduced singlet−triplet energy gap ΔST, which thereby makes sensitization difficult and the energy-transfer mechanism inaccessible.11 We recently reported the synthesis of an ethynylsilyl-substituted bistetracene (BT)28 with a hole mobility as large as 6.1 cm2/(V s) and more stable than pentacene in chlorinated solvents. A recent study on the Diels−Alder reactivity of pentacene and bistetracene with C60 by Cao et al. showed that these reactivities are related to distortion energies.14 Here, our goal is to investigate the origin of the increased stability of this substituted bistetracene compared to pentacene and rubrene with respect to the oxidation reaction, which is one of the major limitations for practical applications of organic semiconductors. The molecular structures are given in Figure 2.

ΔEact = ΔE int + ΔEdist

(1)

The distortion energy is the energy required to distort the reactants from their ground-state structures to those in the transition state without allowing for interaction between the reactants. To calculate ΔEdist, each transition state structure is separated into two fragments, 8505

DOI: 10.1021/acs.chemmater.6b02552 Chem. Mater. 2016, 28, 8504−8512

Article

Chemistry of Materials and single-point calculations are carried out on each of the fragments. The energy difference between the distorted structure in the transition state (or at any point in the IRC path) and the optimized ground-state structures is defined as the distortion energy. The interaction energy (ΔEint) is then the energy difference between the activation energy and distortion energy. The local aromaticity indices of TIPS-BT, pentacene, and rubrene were calculated at the B3LYP/6-31G* level of theory. We have calculated aromaticity indices such as nuclear independent chemical shift (NICS), 44 harmonic oscillator model of aromaticity (HOMA),45,46 and multicenter index (MCI).47,48 MCI values were calculated within the Becke-rho atomic partition49 and computed using the APOST-3D and ESI-3D programs.50,51 The rate of electron transfer was estimated within the framework of semiclassical Marcus theory52 as k=

2 4π 2 2 1 VRP e[−(ΔG + λ) /4λkBT ] h 4πkBT

ring in the TIPS-BT, pentacene, and rubrene molecules, see Table 2.57 The trends in aromaticity predicted by the NICS, Table 2. NICS (in ppm), HOMA, and MCI Atomicity Indices of Various Rings in Pentacene, Rubrene, and TIPSBTa ring

NICS

HOMA

MCI

1 2 3 1 2 3 1 2 3 4

−7.04 −12.04 −13.49 −7.07 −9.61 −9.11 −10.60 −14.16 −10.81 −6.67

0.48 0.54 0.56 0.52 0.44 0.97 0.63 0.71 0.43 0.53

0.020 0.016 0.016 0.022 0.016 0.057 0.024 0.022 0.011 0.016

rubrene

TIPS-BT

(2)

where ΔG is the Gibbs free energy of the reaction; VRP is the electronic coupling between the final and initial states, and λ is the reorganization energy. The ΔG values for the forward electron-transfer reaction (ΔGf) (PAH* + O2 → PAH+ + O2−) and the backward electron-transfer reaction (ΔGb) (PAH+ + O2− → PAH+ + O2) were calculated as the energy difference between the final and initial adiabatic states.53,54 The electronic couplings between the molecules were evaluated for the lowest energy conformer using the fragment molecular orbital (FMO) approach.55 Approximately 2000 conformations of the systems of interest with O2 were generated through a scan of O2 over the molecule to determine the lowest energy conformer. 2.2. Experimental Methods. For the kinetic NMR study, a solution of n-octyldiisopropylsilyl-ethynylene substituted bistetracene (NODISPA-BT) (4 mg) in CDCl3 (1 mL) was stirred in ambient laboratory light/air at room temperature. 1H NMR measurements were performed with a Bruker Avance 400 at 3, 12, and 24 days of exposure.

a

Table 1. Calculated ΔST Values of O2, TIPS-BT, Pentacene, and Rubrene at the ωB97X-D/6-31G* Levela molecule

ΔST(vertical)

ΔST(adiabatic)

O2 TIPS-BT pentacene rubrene

1.02 1.10 1.22 1.51

1.02 0.70 0.85 1.08

Ring numbering is given in Figure 2.

HOMA, and MCI values, when compared to available experimental observations for pentacene and rubrene, suggest that the MCI approach is a more reliable aromaticity index (see Supporting Information for a more detailed discussion). We note that, in addition to the aforementioned aromaticity indices, in which both MCI and HOMA predict that ring 1 has the largest aromatic character and ring 3 the lesser, we have also calculated additional aromaticity indices proposed in the literature including the para-delocalization index (PDI)58 as well as the Iring,59 INB,60 and ING indexes.60 In the case of TIPSBT, we find all these indices to be in agreement with the MCI results (Table 1 in the Supporting Information). Hence, on the basis of aromaticity criteria grounded on MCI results, ring 3 is expected to be more prone to an oxidation reaction than rings 1, 2, or 4. To assess the role of aromaticity in controlling the stability of the molecules under investigation, we have evaluated the activation energy for endoperoxide formation on each ring of the aromatic molecule using transition state theory. We have investigated both the concerted and stepwise reaction mechanisms of TIPS-BT, pentacene, and rubrene for each of their individual rings. We find that the concerted Diels−Alder oxidation reaction is favorable compared to the stepwise mechanism (see Figure 7 and Figure S1 in the Supporting Information), with the exception of rubrene (Figure S2 in the Supporting Information), where the concerted and stepwise mechanisms are found to have similar activation barriers. For TIPS-BT, the oxidation of ring 2 is found to be the most favorable as it corresponds to the lowest activation energy barrier (ΔG‡ = 17.7 kcal/mol via a concerted path, BTTSr2, as shown in Figure 3); it occurs in a concerted fashion, with the reaction being exergonic with a free energy gain of 19.4 kcal/ mol to form a stable endoperoxide product (BT4r2; Figure 3). The ΔG‡ values for the concerted pathway of the oxidation reaction at ring 1, ring 3, and ring 4 of TIPS-BT are 27.0, 19.3, and 31.0 kcal/mol, respectively (see Figure S1 in the Supporting Information). The stepwise and concerted mechanisms of endoperoxide formation of ring 2 of TIPS-BT, the most energetically favorable, are illustrated in Figure 3; the stepwise mechanism proceeds via a radical coupling reaction with a biradical nature for BT2r2 and BT3r2, which has been verified via a spin density analysis, see Figure 4.

3. RESULTS AND DISCUSSIONS To investigate the possibility of energy transfer from TIPS-BT, pentacene, and rubrene toward O2, both the vertical and adiabatic ΔST values of O2, TIPS-BT, pentacene, and rubrene were calculated at the ωB97X-D/6-31G* level; the results are collected in Table 1. As mentioned above, the ωB97X-D/6-

a

system pentacene

All energies are given in eV.

31G* ΔST value of O2 (1.02 eV) is in excellent agreement with the experimental excitation energy (0.98 eV). The long lifetime of the excited triplet state and the large range of ΔST values (0.4 eV) of TIPS-BT, pentacene, and rubrene, see Table 1, enable them to undergo resonance energy transfer with O2 followed by photo-oxidation reaction. From a purely qualitative standpoint, the increased stability of the TIPS-BT compared to that of pentacene has been previously attributed to the fact that the former has two Clar aromatic sextets, whereas the latter has only one (we recall that a larger number of Clar aromatic sextets is considered to increase the stability of aromatic molecules).56 In this context, we have chosen to calculate the local aromaticity indices of each 8506

DOI: 10.1021/acs.chemmater.6b02552 Chem. Mater. 2016, 28, 8504−8512

Article

Chemistry of Materials

Figure 3. Concerted (solid lines) and stepwise (dashed lines) reaction pathways of 1O2 with TIPS-BT leading to endoperoxide formation on ring 2. Gibbs energies (in kcal/mol) are given for each state.

Figure 4. Illustration of the ωB97X-D spin densities in the BT2r2 and BT 3r2 states of Figure 3 (the red and blue colors represent spin-up and spin-down electron densities, respectively).

Interestingly, according to the aromaticity analysis described above, ring 3 is the least aromatic (as predicted by MCI and HOMA values) and would be expected to undergo oxidation more readily than ring 2, while the activation free energy barriers suggest that the endoperoxide formation on ring 2 is more favorable by 1.6 kcal/mol. (We note that results for activation free energies obtained using the M06-2X functional are close to ωB97X-D values, see Table S2 in the Supporting Information.) To understand the origin of this regioselectivity of endoperoxide formation, we have carried out distortion/ interaction analyses for the concerted transition states of endoperoxide formation on rings 1, 2, 3, and 4 of TIPS-BT, see Figure 5.40,41 From Figure 5, it appears that the distortion energy is minimum on ring 2 (11.9 kcal/mol) and is 2.2 kcal/ mol smaller than that on ring 3 (14.1 kcal/mol). However, the interaction energy is similar on ring 2 and ring 3. This is in agreement with the recent study on Diels−Alder reaction of bis(2-(trimethylsilyl)ethynyl) bistetracene with C60 where Cao et al. found that the lowest distortion energy is on ring 2.14

Figure 5. Distortion/interaction (activation strain) analysis for the concerted transition state of 1O2 with ring 1, ring 2, ring 3, and ring 4 of TIPS-BT. ΔEdist, ΔEint, and ΔEact are given in kcal/mol.

These findings again highlight that aromaticity is not the main factor controlling the stability; the distortion energy required to access the transition state plays a critical role. The distortion energies of the transition states for the TIPS-BT oxidation on various rings correlate well with the transition state energies, which points to the most favorable sites for TIPS-BT oxidation. Moreover, the distortion energy difference (2.2 kcal/mol) calculated between ring 2 and ring 3 is consistent with the activation free-energy difference (1.6 kcal/mol). These results underline that a balance between distortion and interaction energies, instead of aromaticity, actually determines the regioselectivity of endoperoxide formation in TIPS-BT. (We note that we also carried out similar calculations with water and 8507

DOI: 10.1021/acs.chemmater.6b02552 Chem. Mater. 2016, 28, 8504−8512

Article

Chemistry of Materials

Figure 6. (A) Proposed endoperoxide formation upon photooxidation of n-octyldiisopropylsilyl-ethynylene substituted bistetracene (NODISPABT). (B) 1H NMR spectra monitoring the photooxidation of NODISPSA-BT with various time intervals.

Figure 7. Concerted (solid lines) and stepwise (dashed lines) reaction pathways of 1O2 with pentacene. Free energies are given in kcal/mol.

8508

DOI: 10.1021/acs.chemmater.6b02552 Chem. Mater. 2016, 28, 8504−8512

Article

Chemistry of Materials

Figure 8. Concerted (solid lines) and stepwise (dashed lines) reaction pathways of 1O2 with rubrene. Free energies are given in kcal/mol.

with oxygen and focused on the central ring; both the concerted and stepwise reaction pathways for the Diels−Alder reaction of pentacene with 1O2 are illustrated in Figure 7. The activation free energy barrier for the concerted reaction is 13.6 kcal/mol (via PTSr3), which is 16.2 kcal/mol lower than that of the stepwise mechanism (via PTS1r3, ΔG = 29.8 kcal/mol), and the reaction is exergonic (ΔGrel= −31.6 kcal/mol, Figure 7). This result suggests that a concerted mechanism is preferred over the stepwise mechanism that was previously reported.25 The activation energy barrier (ΔG‡) for the endoperoxide formation on pentacene (13.6 kcal/mol) is smaller than that on TIPS-BT (17.7 kcal/mol), which is in agreement with the experimentally observed higher stability (half-life) of pentacene than that of TIPS-BT toward photooxidation reaction. In the case of rubrene, our calculations indicate that the Diels−Alder reaction with 1O2 occurs on ring 2, in agreement with experimental observations.22 The concerted path has a ΔG‡ of 14.4 kcal/mol (via RTSr3), which is lower than that of TIPS-BT and higher than that of pentacene (Figure 8). It is interesting to note that ΔG‡ of the stepwise mechanism (via R TS1r3) for the endoperoxide formation at ring 2 of rubrene (15.3 kcal/mol) is similar to that of the concerted mechanism; hence, both of these mechanisms appear to be energetically possible in rubrene. The resulting rubrene endoperoxide is stable by about 20.7 kcal/mol. The biradical nature of the intermediates and transition states of the stepwise mechanism has also been verified via a spin density analysis, see Figure S3 in the Supporting Information. For the concerted pathways of rubrene oxidation at ring 1 and ring 3, see Figure S2 in the Supporting Information. The experimentally observed half-lives of TIPS-BT, pentacene, and rubrene in solution (5 days,28 3 min,62 and 40 min,12 respectively) correlate well with our calculated ΔG‡ values of 17.7, 13.6, and 14.4 kcal/mol, respectively. As stated above, the activation barrier for the oxidation of TIPS-BT is 4.1 kcal/mol higher (at ring 2) than that of pentacene (at ring 3); this corresponds to an approximately 1000-times higher stability of TIPS-BT than pentacene and is well in line with

found that it does not affect the photostability of the neutral TIPS-BT molecule.) The proposed regioselectivity of endoperoxide formation based on the activation energy analysis was subsequently confirmed by the results of 1H NMR experiments on the noctyldiisopropylsilyl-ethynylene substituted bistetracene (NODIPSA-BT), whose structural backbone is similar to that of TIPS-BT, see Figure 6 (we note that NODIPSA-BT was chosen for this experiment as it is more soluble in chloroform than TIPS-BT itself). Starting with the pristine PAH in the bottom spectrum,16 this compound slowly converts into the endoperoxide after several days, which indicates that NODIPSA-BT has a large activation barrier against oxidation. The new peaks arising at 6.75 and 7.0 ppm with a relative intensity ratio of 1:1 are indicative of the protons on ring 2 being shifted further upfield, a consequence of having the endoperoxide forming on that ring. As the oxidized product becomes more concentrated in solution with time, peaks in the further downfield positions broaden as a result of the lower molecular symmetry; discrete individual peaks are no longer observed as they become wider and merge. The spectra at 12 and 24 days are a zoom-in view to clearly observe this peak broadening, which is the reason why the solvent peak (at ∼7.3 ppm) as well as all other peaks appear to grow in intensity. Isolating the endoperoxide product did not prove to be possible, but the observation of two new peaks arising and the simultaneous broadening of all other aromatic peaks further downfield gave proof that our species only decomposed after an extended amount of time in solution. We note that the recent study by Cao et al. on the Diels−Alder reaction of C60 with TIPS-BT also showed that the lowest barrier for the reaction is on ring 2,14 which is in line with our experimental observation of Diels−Alder reaction of 1O2 with NODIPSA-BT. To provide a baseline for comparison, we have also considered pentacene and rubrene. From earlier experimental61 and theoretical studies25 of pentacene, it has been shown that endoperoxide formation occurs at the central ring. Therefore, we have limited our investigations of the reaction of pentacene 8509

DOI: 10.1021/acs.chemmater.6b02552 Chem. Mater. 2016, 28, 8504−8512

Article

Chemistry of Materials

Table 3. Gibbs Free Energies of Forward Electron Transfer Reaction (ΔGf), Backward Electron Transfer Reaction (ΔGb), Electronic Couplings between Donor and Acceptor Molecule (VRP), Reorganization Energy (λ), and Rate of Electron Transfer (k) Reactionsa molecule

ΔGf [ΔGb] (eV)

|VRP|f [|VRP|b] (meV)

λf [λb] eV

kf [kb] s−1

TIPS-BT pentacene rubrene

−0.46 [−1.67] −0.84 [−1.50] −1.74 [−0.60]

7 [17] 0.004 [0.37] 104 [442]

1.24 [1.26] 1.26 [1.31] 1.26 [1.31]

6.5 × 109 [1.1 × 1012] 6.1 × 104 [1.5 × 109] 2.8 × 1013 [6.6 × 1013]

a

Parameters for the forward electron transfer reaction are denoted by f subscript, and those for the backward electron recombination reaction are denoted by b subscript.

NMR chemical shift analyses validate the theoretical prediction by confirming that endoperoxide formation occurs on the ring with smallest amount of distortion energy. Our work underlines that to enhance the stability of polyaromatic hydrocarbons toward oxidation: (i) the singlet− triplet energy gap (ΔST) of the organic molecules should be substantially smaller than that of oxygen; and (ii) the distortion energy of the aromatic rings should be increased, for example, via incorporation of a two-dimensional PAH framework, which is consistent with the increased stability of the 2D-like bistetracene.

the experimental observation. Moreover, the distortion energy of the concerted oxidation on ring 2 of TIPS-BT (11.9 kcal/ mol) is higher than that on ring 3 of pentacene (6.1 kcal/mol) and on ring 2 of rubrene (9.2 kcal/mol). Even though a similar photodegradation mechanism is expected in the solid-state environment, the crystal packing in the solid state has to be taken into account for a proper understanding of device degradation.63 Finally, we have also investigated the electron transfer mechanisms shown in Figure 1. The rates of electron transfer between TIPS-BT, pentacene, or rubrene with 3O2 for the forward and backward reactions are calculated using eq 2 and presented in Table 3, along with the parameters entering the above equation. For pentacene and TIPS-BT, the rate of forward reaction to create the charged species (PAH+ + O−2 ) is at least two orders of magnitude smaller than that of their recombination (backward reaction). Thus, we expect that degradation in pentacene and TIPS-BT does not involve an electron transfer mechanism (type I). However, in rubrene, both the forward and backward reactions are of the same order, again in agreement with previous experimental detection of the rubrene radical cation and oxygen radical anion using EPR spectroscopy, indicating that an electron transfer process is possible.19



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.chemmater.6b02552. Local aromaticity indices; oxidation reaction pathways of TIPS-BT and rubrene at various rings; spin density plots (PDF)



AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected].

4. CONCLUSION We have described a detailed investigation, at the long-range corrected density functional theory level, of the [4 + 2] cycloaddition (Diels−Alder) reaction mechanism for a recently synthesized bistetracene derivative. Our aim was to elucidate the origin of the increased stability toward oxidation displayed by the substituted bistetracene compared to pentacene and rubrene. We find that the singlet−triplet gaps of pentacene, rubrene, and TIPS-BT are similar to the excitation energy of molecular oxygen and can lead to singlet oxygen generation, which opens up a pathway for degradation. The calculated activation free energy barrier of the oxidation reaction in TIPS-BT is 4.1 kcal/ mol higher than that of pentacene, and 3.2 kcal/mol higher than that of rubrene. Local aromaticity indices calculated with electron sharing indices (MCIs) are found to be more reliable than those based on NICS. Interestingly, if aromaticity were the main factor controlling the reaction, degradation of TIPS-BT with singlet oxygen would be expected to take place on ring 3 (which is one of the rings shared by the two tetracene subunits). However, our results show that the free energy barrier for oxidation on ring 2 (pertaining to a single tetracene subunit) is lower by 1.6 kcal/mol than on ring 3. We were able to rationalize these results on the basis of a distortion/interaction energy analysis, which points out that the difference in activation energy primarily originates in distortion energy. Experimental 1H

ORCID

Simil Thomas: 0000-0002-8069-4940 Jean-Luc Bredas: 0000-0001-7278-4471 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This work was supported by King Abdullah University of Science and Technology (KAUST) and by ONR-Global (Award No. N62909-15-1-2003). We are grateful to the KAUST Supercomputing Laboratory and the KAUST IT Research Computing Team for providing continuous assistance as well as computational and storage resources. J.L. acknowledges funding from the National Science Foundation (DMR1508627); L.Z. thanks the Office of Naval Research (N00014714-1-0053); S.T. thanks Dr. Manjaly J. Ajitha for useful discussions.



REFERENCES

(1) Yoo, S.; Domercq, B.; Kippelen, B. Efficient thin-film organic solar cells based on pentacene/C60 heterojunctions. Appl. Phys. Lett. 2004, 85, 5427−5429. (2) Kelley, T. W.; Muyres, D. V.; Baude, P. F.; Smith, T. P.; Jones, T. D. High Performance Organic Thin Film Transistors. MRS Online Proc. Libr. 2003, 771, 169.

8510

DOI: 10.1021/acs.chemmater.6b02552 Chem. Mater. 2016, 28, 8504−8512

Article

Chemistry of Materials (3) 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. (4) Podzorov, V.; Menard, E.; Borissov, A.; Kiryukhin, V.; Rogers, J. A.; Gershenson, M. E. Intrinsic Charge Transport on the Surface of Organic Semiconductors. Phys. Rev. Lett. 2004, 93, 086602. (5) Sundar, V. C.; Zaumseil, J.; Podzorov, V.; Menard, E.; Willett, R. L.; Someya, T.; Gershenson, M. E.; Rogers, J. A. Elastomeric Transistor Stamps: Reversible Probing of Charge Transport in Organic Crystals. Science 2004, 303, 1644−1646. (6) Podzorov, V.; Menard, E.; Rogers, J. A.; Gershenson, M. E. Hall Effect in the Accumulation Layers on the Surface of Organic Semiconductors. Phys. Rev. Lett. 2005, 95, 226601. (7) Schön, J. H.; Kloc, C.; Batlogg, B. Fractional Quantum Hall Effect in Organic Molecular Semiconductors. Science 2000, 288, 2338−2340. (8) Reyes-Martinez, M. A.; Crosby, A. J.; Briseno, A. L. Rubrene crystal field-effect mobility modulation via conducting channel wrinkling. Nat. Commun. 2015, 6, 6948. (9) Dimitrakopoulos, C. D.; Malenfant, P. R. L. Organic Thin Film Transistors for Large Area Electronics. Adv. Mater. 2002, 14, 99−117. (10) Nomura, K.; Ohta, H.; Takagi, A.; Kamiya, T.; Hirano, M.; Hosono, H. Room-temperature fabrication of transparent flexible thinfilm transistors using amorphous oxide semiconductors. Nature 2004, 432, 488−492. (11) Maliakal, A.; Raghavachari, K.; Katz, H.; Chandross, E.; Siegrist, T. Photochemical Stability of Pentacene and a Substituted Pentacene in Solution and in Thin Films. Chem. Mater. 2004, 16, 4980−4986. (12) Uttiya, S.; Raimondo, L.; Campione, M.; Miozzo, L.; Yassar, A.; Moret, M.; Fumagalli, E.; Borghesi, A.; Sassella, A. Stability to photooxidation of rubrene and fluorine-substituted rubrene. Synth. Met. 2012, 161, 2603−2606. (13) Zade, S. S.; Bendikov, M. Reactivity of acenes: mechanisms and dependence on acene length. J. Phys. Org. Chem. 2012, 25, 452−461. (14) Cao, Y.; Liang, Y.; Zhang, L.; Osuna, S.; Hoyt, A.-L. M.; Briseno, A. L.; Houk, K. N. Why Bistetracenes Are Much Less Reactive Than Pentacenes in Diels−Alder Reactions with Fullerenes. J. Am. Chem. Soc. 2014, 136, 10743−10751. (15) 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. (16) Zschieschang, U.; Ante, F.; Yamamoto, T.; Takimiya, K.; Kuwabara, H.; Ikeda, M.; Sekitani, T.; Someya, T.; Kern, K.; Klauk, H. Flexible Low-Voltage Organic Transistors and Circuits Based on a High-Mobility Organic Semiconductor with Good Air Stability. Adv. Mater. 2010, 22, 982−985. (17) Yousif, E.; Haddad, R. Photodegradation and photostabilization of polymers, especially polystyrene: review. SpringerPlus 2013, 2, 1− 32. (18) Hintz, H.; Egelhaaf, H. J.; Lüer, L.; Hauch, J.; Peisert, H.; Chassé, T. Photodegradation of P3HT−A Systematic Study of Environmental Factors. Chem. Mater. 2011, 23, 145−154. (19) Maliakal, A. J.; Chen, J. Y. C.; So, W.-Y.; Jockusch, S.; Kim, B.; Ottaviani, M. F.; Modelli, A.; Turro, N. J.; Nuckolls, C.; Ramirez, A. P. Mechanism for Oxygen-Enhanced Photoconductivity in Rubrene: Electron Transfer Doping. Chem. Mater. 2009, 21, 5519−5526. (20) Abdou, M. S. A.; Holdcroft, S. Mechanisms of photodegradation of poly(3-alkylthiophenes) in solution. Macromolecules 1993, 26, 2954−2962. (21) Adam, W.; Saha-Möller, C. R.; Schönberger, A. Type I and Type II Photosensitized Oxidative Modification of 2‘-Deoxyguanosine (dGuo) by Triplet-Excited Ketones Generated Thermally from the 1,2-Dioxetane HTMD. J. Am. Chem. Soc. 1997, 119, 719−723. (22) Nardello, V.; Marti, M.-J.; Pierlot, C.; Aubry, J.-M. Photochemistry without Light: Oxidation of Rubrene in a Microemulsion with a Chemical Source of Singlet Molecular Oxygen (1O2, 1Dg). J. Chem. Educ. 1999, 76, 1285. (23) Hurst, J. R.; McDonald, J. D.; Schuster, G. B. Lifetime of singlet oxygen in solution directly determined by laser spectroscopy. J. Am. Chem. Soc. 1982, 104, 2065−2067.

(24) Chien, S.-H.; Cheng, M.-F.; Lau, K.-C.; Li, W.-K. Theoretical Study of the Diels−Alder Reactions between Singlet (1Δg) Oxygen and Acenes. J. Phys. Chem. A 2005, 109, 7509−7518. (25) Reddy, A. R.; Bendikov, M. Diels-Alder reaction of acenes with singlet and triplet oxygen - theoretical study of two-state reactivity. Chem. Commun. 2006, 1179−1181. (26) Anthony, J. E.; Brooks, J. S.; Eaton, D. L.; Parkin, S. R. Functionalized Pentacene: Improved Electronic Properties from Control of Solid-State Order. J. Am. Chem. Soc. 2001, 123, 9482−9483. (27) Okamoto, T.; Senatore, M. L.; Ling, M. M.; Mallik, A. B.; Tang, M. L.; Bao, Z. Synthesis, Characterization, and Field-Effect Transistor Performance of Pentacene Derivatives. Adv. Mater. 2007, 19, 3381− 3384. (28) Zhang, L.; Fonari, A.; Liu, Y.; Hoyt, A.-L. M.; Lee, H.; Granger, D.; Parkin, S.; Russell, T. P.; Anthony, J. E.; Brédas, J.-L.; Coropceanu, V.; Briseno, A. L. Bistetracene: An Air-Stable, High-Mobility Organic Semiconductor with Extended Conjugation. J. Am. Chem. Soc. 2014, 136, 9248−9251. (29) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A.; Nakatsuji, H.; Caricato, M.; Li, X.; Hratchian, H. P.; Izmaylov, A. F.; Bloino, J.; Zheng, G.; Sonnenberg, J. L.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M.; Heyd, J. J.; Brothers, E.; Kudin, K. N.; Staroverov, V. N.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Rega, N.; Millam, J. M.; Klene, M.; Knox, J. E.; Cross, J. B.; Bakken, V.; Adamo, C.; Jaramillo, J.; Gomperts, R.; Stratmann, R. E.; Yazyev, O.; Austin, A. J.; Cammi, R.; Pomelli, C.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Zakrzewski, V. G.; Voth, G. A.; Salvador, P.; Dannenberg, J. J.; Dapprich, S.; Daniels, A. D.; Farkas, O.; Foresman, J. B.; Ortiz, J. V.; Cioslowski, J.; Fox, D. J. Gaussian 09, revision D.01; Gaussian, Inc.: Wallingford, CT, 2009. (30) Hariharan, P. C.; Pople, J. A. The influence of polarization functions on molecular orbital hydrogenation energies. Theoret. Chim. Acta 1973, 28, 213−222. (31) Francl, M. M.; Pietro, W. J.; Hehre, W. J.; Binkley, J. S.; Gordon, M. S.; DeFrees, D. J.; Pople, J. A. Self-consistent molecular orbital methods. XXIII. A polarization-type basis set for second-row elements. J. Chem. Phys. 1982, 77, 3654−3665. (32) Chai, J.-D.; Head-Gordon, M. Long-range corrected hybrid density functionals with damped atom-atom dispersion corrections. Phys. Chem. Chem. Phys. 2008, 10, 6615−6620. (33) Becke, A. D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 1993, 98, 5648−5652. (34) Yanai, T.; Tew, D. P.; Handy, N. C. A new hybrid exchange− correlation functional using the Coulomb-attenuating method (CAMB3LYP). Chem. Phys. Lett. 2004, 393, 51−57. (35) Zhao, Y.; Truhlar, D. G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215−241. (36) Scalmani, G.; Frisch, M. J. Continuous surface charge polarizable continuum models of solvation. I. General formalism. J. Chem. Phys. 2010, 132, 114110. (37) Goldstein, E.; Beno, B.; Houk, K. N. Density Functional Theory Prediction of the Relative Energies and Isotope Effects for the Concerted and Stepwise Mechanisms of the Diels−Alder Reaction of Butadiene and Ethylene. J. Am. Chem. Soc. 1996, 118, 6036−6043. (38) Xu, X.; Goddard, W. A. Peroxone chemistry: Formation of H2O3 and ring-(HO2) (HO3) from O3/H2O2. Proc. Natl. Acad. Sci. U. S. A. 2002, 99, 15308−15312. (39) Legault, C. Y.; Garcia, Y.; Merlic, C. A.; Houk, K. N. Origin of Regioselectivity in Palladium-Catalyzed Cross-Coupling Reactions of Polyhalogenated Heterocycles. J. Am. Chem. Soc. 2007, 129, 12664− 12665. 8511

DOI: 10.1021/acs.chemmater.6b02552 Chem. Mater. 2016, 28, 8504−8512

Article

Chemistry of Materials (40) Ess, D. H.; Houk, K. N. Theory of 1,3-Dipolar Cycloadditions: Distortion/Interaction and Frontier Molecular Orbital Models. J. Am. Chem. Soc. 2008, 130, 10187−10198. (41) Green, A. G.; Liu, P.; Merlic, C. A.; Houk, K. N. Distortion/ Interaction Analysis Reveals the Origins of Selectivities in IridiumCatalyzed C−H Borylation of Substituted Arenes and 5-Membered Heterocycles. J. Am. Chem. Soc. 2014, 136, 4575−4583. (42) Fernandez, I.; Bickelhaupt, F. M. The activation strain model and molecular orbital theory: understanding and designing chemical reactions. Chem. Soc. Rev. 2014, 43, 4953−4967. (43) Fernández, I.; Bickelhaupt, F. M. Alder-ene reaction: Aromaticity and activation-strain analysis. J. Comput. Chem. 2012, 33, 509−516. (44) Schleyer, P. v. R.; Maerker, C.; Dransfeld, A.; Jiao, H.; Hommes, N. J. R. v. E. Nucleus-Independent Chemical Shifts: A Simple and Efficient Aromaticity Probe. J. Am. Chem. Soc. 1996, 118, 6317−6318. (45) Kruszewski, J.; Krygowski, T. M. Definition of aromaticity basing on the harmonic oscillator model. Tetrahedron Lett. 1972, 13, 3839−3842. (46) Krygowski, T. M. Crystallographic studies of inter- and intramolecular interactions reflected in aromatic character of.pi.electron systems. J. Chem. Inf. Model. 1993, 33, 70−78. (47) Bultinck, P.; Ponec, R.; Van Damme, S. Multicenter bond indices as a new measure of aromaticity in polycyclic aromatic hydrocarbons. J. Phys. Org. Chem. 2005, 18, 706−718. (48) Feixas, F.; Matito, E.; Poater, J.; Sola, M. Quantifying aromaticity with electron delocalisation measures. Chem. Soc. Rev. 2015, 44, 6434−6451. (49) Becke, A. D. A multicenter numerical integration scheme for polyatomic molecules. J. Chem. Phys. 1988, 88, 2547−2553. (50) Salvador, P. R.-C. E. APOST-3D; Universitat de Girona, 2013. (51) Matito, E. ESI-3D: Electron Sharing Indices Program for 3D Molecular Space Partitioning; Institute of Computational chemistry and Catalysis (IQCC), University of Girona: Catalonia, Spain, 2006. (52) Marcus, R. A. Electron transfer reactions in chemistry. Theory and experiment. Rev. Mod. Phys. 1993, 65, 599−610. (53) Lemaur, V.; Steel, M.; Beljonne, D.; Brédas, J.-L.; Cornil, J. Photoinduced Charge Generation and Recombination Dynamics in Model Donor/Acceptor Pairs for Organic Solar Cell Applications: A Full Quantum-Chemical Treatment. J. Am. Chem. Soc. 2005, 127, 6077−6086. (54) Northrop, B. H.; Houk, K. N.; Maliakal, A. Photostability of pentacene and 6,13-disubstituted pentacene derivatives: a theoretical and experimental mechanistic study. Photochem. Photobiol. Sci. 2008, 7, 1463−1468. (55) Valeev, E. F.; Coropceanu, V.; da Silva Filho, D. A.; Salman, S.; Brédas, J.-L. Effect of Electronic Polarization on Charge-Transport Parameters in Molecular Organic Semiconductors. J. Am. Chem. Soc. 2006, 128, 9882−9886. (56) Clar, E. The Aromatic Sextet; John Wiley & Sons: London, 1972. (57) Garcia-Borras, M.; Osuna, S.; Luis, J. M.; Swart, M.; Sola, M. The role of aromaticity in determining the molecular structure and reactivity of (endohedral metallo)fullerenes. Chem. Soc. Rev. 2014, 43, 5089−5105. (58) Poater, J.; Fradera, X.; Duran, M.; Solà, M. The Delocalization Index as an Electronic Aromaticity Criterion: Application to a Series of Planar Polycyclic Aromatic Hydrocarbons. Chem. - Eur. J. 2003, 9, 400−406. (59) Giambiagi, M.; Segre de Giambiagi, M.; dos Santos Silva, C. D.; Paiva de Figueiredo, A. Multicenter bond indices as a measure of aromaticity. Phys. Chem. Chem. Phys. 2000, 2, 3381−3392. (60) Cioslowski, J.; Matito, E.; Solà, M. Properties of Aromaticity Indices Based on the One-Electron Density Matrix†. J. Phys. Chem. A 2007, 111, 6521−6525. (61) Kaur, I.; Jia, W.; Kopreski, R. P.; Selvarasah, S.; Dokmeci, M. R.; Pramanik, C.; McGruer, N. E.; Miller, G. P. Substituent Effects in Pentacenes: Gaining Control over HOMO−LUMO Gaps and Photooxidative Resistances. J. Am. Chem. Soc. 2008, 130, 16274− 16286.

(62) Kim, H. G.; Choi, H. H.; Song, E.; Cho, K.; Choi, E. J. Synthesis, stability and electrical properties of new soluble pentacenes with unsaturated side groups. RSC Adv. 2015, 5, 8070−8076. (63) Yamada, M.; Ikemoto, I.; Kuroda, H. Photooxidation of the Evaporated Films of Polycyclic Aromatic Hydrocarbons Studied by XRay Photoelectron Spectroscopy. Bull. Chem. Soc. Jpn. 1988, 61, 1057−1062.

8512

DOI: 10.1021/acs.chemmater.6b02552 Chem. Mater. 2016, 28, 8504−8512