Stacking Principles on π- and Lamellar Stacking for Organic

Dec 28, 2018 - ... lamellar stacking, are generally found for organic semiconductors, ... The stacking principles are associated with the closest pack...
1 downloads 0 Views 1MB Size
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.

Article Cite This: ACS Omega 2018, 3, 18656−18662

http://pubs.acs.org/journal/acsodf

Stacking Principles on π- and Lamellar Stacking for Organic Semiconductors Evaluated by Energy Decomposition Analysis Yu-Ying Lai,* Vi-Hsiang Huang, Hao-Ting Lee, and Hau-Ren Yang Institute of Polymer Science and Engineering, National Taiwan University, Taipei 10617, Taiwan

Downloaded via 46.161.63.119 on January 1, 2019 at 07:41:21 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

S Supporting Information *

ABSTRACT: Two stacking manners, that is, π- and lamellar stacking, are generally found for organic semiconductors, in which the π-stacking occurs between conjugated groups and the lamellar stacking refers to the separation of the conjugated and aliphatic moieties. The stacking principles are yet not well-defined. In this work, extended transition state−natural orbitals for chemical valence (ETS−NOCV), an energy decomposition analysis, is utilized to examine the π- and lamellar stacking for a series of naphthalenetetracarboxylic diimide (R-NDI) crystals. The crucial role of dispersion is validated. The perception that π-stacking is merely determined by the conjugated moiety is challenged. The stacking principles are associated with the closest packing model. Nanoscopic phase separation of conjugated and aliphatic moieties and the formation of lamellar and herringbone motifs in the R-NDIs can thus be clarified. Moreover, the interactions between NDI and the alkyl chain are investigated, revealing that the interactions can be significant, being contradictory to the conventional point of view. Along with R-NDIs, additional organic crystals consisting of various conjugated functionalities and substituents are also investigated by ETS−NOCV. The sampling scope is up to 108 conjugated molecules. The dominant role of dispersion force irrespective of the variation in the conjugated moieties and substituents is further confirmed. It is envisaged that the established principles are applicable to other organic semiconductors. The perspective toward the π- and lamellar stacking might be modified, paving the way for ultimate morphological control.



Scheme 1. Illustrations for π- and Lamellar Stacking: (a) Lamellar Motif and (b) Herringbone Motifa

INTRODUCTION Organic semiconductors have been used extensively in organic photovoltaics,1,2 organic field-effect transistors,3−5 organic light-emitting diodes,6−8 organic biomedical devices,9,10 and so on. The device performance is associated with numerous factors, such as chemical functionality,11 chemical structure,12,13 molecular weight,14−16 morphology,17 and device architecture.18,19 Among all, the effect of morphology is often prominent. The morphology can be quantified by the parameter of order. One extreme is amorphous, and the other is crystalline. It is common for organic semiconductors to have both amorphous and crystalline regions, especially in solution-processed organic thin films. Two stacking manners, that is, π- and lamellar stacking, are generally found for the order domains, in which the π-stacking occurs between conjugated groups and the lamellar stacking usually refers to the regular separation of the conjugated and aliphatic moieties (Scheme 1).20 Moreover, two packing motifs lamella (Scheme 1a) and herringbone (Scheme 1b) are typically utilized to describe the relationship between the individual π-channels along the direction of the lamellar stacking. To avoid confusion, it has to be noted that the lamellar motif is used to represent the pattern, being different from the lamellar stacking (Scheme 1). Electrostatic interaction and direct interaction have been proposed to account for the formation of π-stacking. In the © 2018 American Chemical Society

a

Other lamellar motifs exist.

former, the electron density of conjugated systems is modulated by the attaching substituent, thus giving rise to the variation in the electrostatic interactions between conjugated systems.21−25 The latter stresses that direct interactions exist between the substituent and the conjugated Received: October 9, 2018 Accepted: December 14, 2018 Published: December 28, 2018 18656

DOI: 10.1021/acsomega.8b02713 ACS Omega 2018, 3, 18656−18662

ACS Omega

Article

system along with the interactions induced from the change in bond polarity because of the presence of the substituent.26−30 It has been examined experimentally and theoretically that the direct interaction model is usually dominant for the πstacking.31,32 The direct interaction could be associated with the enhanced London dispersion force.32−41 On the other hand, the occurrence of lamellar stacking is often rationalized by the compatibility principle, that is, like dissolves like. On account of the polarity difference between the conjugated and aliphatic moieties, the separation takes place accordingly. It is instinctive to relate the polarity effect with the electrostatic interaction. Combination of the extended transition-state (ETS) method42−44 and the natural orbitals for chemical valence (NOCV) theory,45−47 coined as ETS−NOCV, has been employed to decompose the interaction energy (ΔEint) between two molecular entities into three chemically meaningful componentselectrostatic interactions (Vele), Pauli repulsion (EPauli), and orbital interactions (Eoi). Vele is the conventional electrostatic interaction between two molecules. EPauli is the repulsive Pauli interaction between the occupied orbitals of two different molecules. Eoi stands for the stabilizing interactions between the occupied molecular orbitals in one molecule and the unoccupied molecular orbitals in the other molecule along with the interaction of occupied and virtual orbitals within the same molecule in the final geometry. Because traditional density functionals often fail in describing the dispersion force accurately, an additional dispersion correction (Edis) can be included in the above energy decomposition scheme with the aim of taking long-range interactions into consideration.48−56 Herein, numerous N,N′dialkyl-1,4,5,8-naphthalenetetracarboxylic diimides (R-NDIs) are chosen for the study, such as R = C14H29,57 C12H25,58 C6H13,59 C5H11,60 C4H9,58 C2H5,61 and 1-CH3C6H1258 (Scheme 2). The π- and lamellar stacking extracted from

Figure 1. Central C14H29-NDI molecule surrounded by six counterparts in the π-stacking.

surroundings can be classified into three types and the corresponding intermolecular interactions between two molecules are labeled as Ta, Tp, and Tq (Figure 1a). Figure 1b demonstrates the distinct separation of the conjugated and aliphatic moieties in the π-stacking from another viewpoint. Originally, the BP86/TZP level of theory was used to estimate Ta, Tp, and Tq. However, in the absence of Edis, the interactions are commonly repulsive, being contrary to the formation of these single crystals. BP86-BJDAMP/TZP with the inclusion of Grimme’s BJDAMP dispersion correction was applied instead.50 ETS−NOCV results of π-stacking in C14H29-NDI are summarized in Table 1. In the parenthesis following Tx (x = Table 1. ETS−NOCV Results of the π-Stacking in C14H29NDI in kcal mol−1

Scheme 2. Chemical Structure of R-NDI

selected single-crystal structures are examined by ETS− NOCV. The interactions between NDI and the alkyl chain are also investigated computationally. The stacking principles on π- and lamellar stacking for R-NDIs are thus established. Moreover, in conjunction with R-NDIs, additional organic crystals containing a variety of conjugated functionalities and substituents are examined as well. The total sampling scope is up to 108 conjugated molecules. It is envisioned that the principles should apply to other organic semiconductors, paving the way for ultimate morphological control.

a

ΔEinta

EPauli

Vele

Eoi

Edis

1

Ta

−37.00

31.20

2

Ta(NDI)

−18.02

13.54

3

Ta(R)

−15.02

16.19

4

Tp

−28.97

23.46

5

Tp(NDI)

−6.96

5.41

6

Tp(R)

−16.79

12.33

7

Tq

−7.79

5.34

8

Tq(NDI)

−5.96

4.95

9

Tq(R)

−0.94

0.19

−10.71 (16%)b −6.36 (20%) −3.28 (11%) −8.09 (15%) −2.12 (17%) −3.44 (12%) −5.27 (40%) −5.40 (49%) 0.11

−8.04 (12%) −4.36 (14%) −3.46 (11%) −6.48 (12%) −1.67 (14%) −3.22 (11%) −2.06 (16%) −1.95 (18%) −0.10 (8%)

−49.45 (73%) −20.84 (66%) −24.47 (78%) −37.86 (73%) −8.58 (69%) −22.46 (77%) −5.80 (44%) −3.56 (33%) −1.14 (92%)

entry

type

ΔEint = EPauli + Vele + Eoi + Edis. bPercentage of total stabilization.

a, p, or q), NDI denotes that the aliphatic chains are replaced by four hydrogens and R denotes that the NDI groups are substituted by four hydrogens. Tx(NDI) and Tx(R) are utilized to give rough estimation about the individual contribution of conjugated and aliphatic moieties to the ΔEint. Several conclusions can be drawn from Table 1: (1) Edis > Vele > Eoi is found for most decompositions, (2) Edis is the principal stabilization for either Ta, Tp, or Tq, (3) Ta has the strongest stabilization (−37.00 kcal mol−1, entry 1, Table 1), followed by Tp (−28.97 kcal mol−1, entry 4, Table 1) and Tq (−7.79 kcal mol−1, entry 7, Table 1), (4) comparison of Tx(NDI) and Tx(R) reveals that the contribution of the tetradecanyl group to the ΔEint is comparable to that of the NDI functionality, and



RESULTS AND DISCUSSION π-Stacking. It can be concluded from all the examined crystal structures that in the π-stacking, each R-NDI molecule is basically enclosed by 6 equiv. C14H29-NDI is taken for illustration (Figure 1). Symmetry operation reveals that the basic alignments between the central C14H29-NDI and the 18657

DOI: 10.1021/acsomega.8b02713 ACS Omega 2018, 3, 18656−18662

ACS Omega

Article

Table 2. Summary of ETS−NOCV Results for π-Stacking in kcal mol−1 C14H29-NDI C12H25-NDI C6H13-NDI C5H11-NDI C4H9-NDI C2H5-NDI 1-CH3C6H12-NDI

Edis(Ta)a

Ta

Tp

Tq

−37.00 −35.46 −23.99 −22.59 −22.10 −20.83 −21.64

−28.97 −26.14 −13.92 −15.92 −9.31 −8.45 −21.26

−7.79 −7.76 −7.44 −8.62 −6.57 −5.96 −3.29

−49.45 −47.36 −28.65 −26.65 −26.39 −24.26 −26.40

Edis(Tp)a −37.86 −34.42 −16.74 −19.62 −10.92 −10.14 −25.87

b

(73%) (71%) (70%) (66%) (62%) (63%) (70%)

Edis(Tq)a −5.80 −5.74 −7.09 −7.75 −4.33 −3.77 −4.92

(73%) (71%) (72%) (68%) (66%) (65%) (70%)

(44%) (42%) (50%) (47%) (39%) (34%) (69%)

a

Edis(Ta) represents the Edis in Ta, Edis(Tp) represents the Edis in Tp, and Edis(Tq) represents the Edis in Tq. bPercentage of total stabilization.

Table 3. ETS−NOCV Results (kcal mol−1) of the Two C14H29-NDI Molecules Aligned in a Face-to-Face Manner

(5) sum of Tx(NDI) and Tx(R) is all less than Tx, suggesting that the mutual interactions between the tetradecanyl and NDI moieties are not trivial. It is envisaged that the packing mode should be decided by intermolecular interactions. In the case of C14H29-NDI, the dispersion force is the dominant contribution. In theory, dispersion is treated as atom-pairwise empirical function. In other words, the dispersion should be highly related to the molecular contact. The increase in the molecular contact should encourage the dispersion force. It is thus rational to see Ta > Tp > Tq because the molecular contact decreases from Ta to Tp and Tq. Moreover, given that Edis is the principal stabilization rather than Vele, the relatively apolar tetradecanyl group can then result in greater stabilization than the NDI unit (see entries 5−6, Table 1), challenging the common belief that the π-stacking is merely determined by the conjugated moiety. The ETS−NOCV results of C12H25-NDI, C6H13-NDI, C5H11-NDI, C4H9-NDI, C2H5-NDI, and 1-CH3C6H12-NDI are listed in the Supporting Information (Tables S1−S6) and summarized in Table 2. The principal character of dispersion to the ΔEint is generally observed. On the whole, Ta, Tp, and Tq decrease gradually in magnitude from the tetradecanyl group to the ethyl group (Table 2). Tables S1−S6 also reveal that the proportion of Tx(R) to the ΔEint is diminished as the chain gets shorter. Both observations can be justified by the reduction in the molecular contact as the shortening of side chain. Structurally speaking, R-NDI consists of two primary units, that is, the R and NDI groups (Scheme 2). The R group has a linear and flexible molecular shape. The NDI group has yet a rectangular and rigid configuration. For a given space with a finite number of R-NDI molecules, the molecules would arrange themselves to gain maximal intermolecular interactions, of which the dispersion is the most dominant. As mentioned previously, the dispersion is linked with the molecular contact. It is thus envisioned that the dispersion stabilization of the whole system can be maximized by the closest packing. For R-NDIs, as a result of their disparity in the molecular shape, the R and NDI groups would undergo nanoscopic phase separation to accomplish the closest packing (Figure 1b). Along this line, the formation of π-stacking might stem essentially from the similarity in the shape rather than the polarity. It can be regarded as the consequence deriving from the governing role of dispersion over electrostatic interaction. However, it is recognized that molecular steric hindrance is associated with EPauli, which would prevent two molecules from establishing too compact stacking. For instance, in Table 3, two C14H29-NDI molecules are forced to align face-to-face and the resultant ΔEints are all repulsive irrespective of the variation in the intermolecular distance from 3.4 to 4.2 Å, being attributed to the drastically large value of EPauli. Moreover, as

entry

distance (Å)a

ΔEint

EPauli

Vele

Eoi

Edis

1 2 3 4 5

3.4 3.6 3.8 4.0 4.2

1847 944 496 256 130

3800 2166 1227 716 437

−1161 −703 −415 −249 −153

−669 −409 −216 −123 −76

−123 −111 −99 −88 −78

a

Distance between the two C14H29-NDI molecules.

judged by the degree of stabilization, Vele and Eoi should still be able to regulate the molecular position in the π-stacking to a certain extent. Nonetheless, their dominance is just much inferior to Edis. Lamellar Stacking. Single-crystal structures confirm that C14H29-NDI, C12H25-NDI, C6H13-NDI, C4H9-NDI, and 1CH3C6H12-NDI fall into the category of the lamellar motif (Scheme 1a), whereas C5H11-NDI and C2H5-NDI exhibit the herringbone motif (Scheme 1b). With the aim of investigating the intermolecular interactions in the lamellar stacking, we intentionally set a nonhydrogen intermolecular distance of 4 Å as a criterion. For a R-NDI molecule in one π-stacking channel, only the R-NDI molecule in the adjacent π-stacking locating close to or within 4 Å is considered to likely give rise to noticeable intermolecular interactions. On the basis of this definition, in the lamellar stacking, there are two corresponding interacting modes for C14H29-NDI, C12H25-NDI, and C6H13NDI, three modes for C4H9-NDI and 1-CH3C6H12-NDI, and four modes for C5H11-NDI and C2H5-NDI (Figure S1). For illustration, the ETS−NOCV results of C14H29-NDI are listed in Table 4. Edis is still the major stabilization in the ΔEint. Nevertheless, modes a and b both exhibit much weaker stabilization than Ta, Tp, or Tq in the π-stacking (Table 1). Similar results are also observed for C12H25-NDI, C6H13-NDI, C5H11-NDI, C4H9-NDI, C2H5-NDI, and 1-CH3C6H12-NDI (Tables S7−S12). Overall, most of the interactions in the 18658

DOI: 10.1021/acsomega.8b02713 ACS Omega 2018, 3, 18656−18662

ACS Omega

Article

Table 4. ETS−NOCV Results of the Lamellar Stacking in C14H29-NDI in kcal mol−1

mode

ΔEint

EPauli

Vele

Eoi

Edis

a b

−3.17 −0.55

3.16 1.53

−0.76 (12%) −0.24 (12%)

−0.85 (13%) −0.40 (19%)

−4.72 (75%) −1.44 (69%)

lamellar stacking are below −3.00 kcal mol−1 in magnitude, revealing that the contributions of the lamellar stacking to the total stabilization are much smaller than those resulting from the π-stacking. In order to achieve maximal stabilization, one would expect that to reduce the length of lamellar stacking and leave more space for the π-stacking is a rational approach, which would lead to the herringbone motif. C5H11-NDI and C2H5-NDI belong to this group. Figure 2 clearly shows that in

stabilization, accounting for the construction of a lamellar motif for C4H9-NDI. Interactions between NDI and the Alkyl Chain. The πstacking and lamellar stacking are investigated individually hitherto. It is intriguing to further look into the interactions between NDI and the alkyl chain. Nevertheless, the corresponding crystal structures are not experimentally accessible. Theoretical approaches are employed instead. Geometry optimizations were performed between the NDI molecule and various alkanes at the BP86-BJDAMP/TZP level of theory. The minimum nature of each optimized geometry was confirmed by frequency calculation. The computational results reveal that the NDI molecule is able to form stable adducts with different alkanes and Edis still plays the most important role (Table 5). The results strongly suggest that the interactions between NDI and the alkyl chain can be substantial, especially for the extensive alkyl chain, which again challenges the conventional point of view that a conjugated moiety interacts barely with an aliphatic chain. As expected, the ΔEint increases from ethane (entry 1, Table 5) to all C6H14 isomers (entries 2−6, Table 5) and tetradecane (entry 7, Table 5) because of the increase in the molecular contact. Comparison of all C6H14 isomers (entries 2−6, Table 5) indicates that the ΔEint is sensitive to the molecular shape. Branched structures could lead to greater stabilization. Overall, these results provide justification for the use of long and branch side chains to increase the solubility for conjugated molecules and polymers by hindering the π-stacking.62−64 For comparison, the NDI−NDI adduct is optimized, and its ΔEint is calculated to be −24.52 kcal mol−1 (entry 8, Table 5), which is much higher in magnitude than the others (entries 1−7, Table 5). It is envisaged that in order to maximize the intermolecular interactions, the NDI molecule would prefer stacking with itself rather than the aliphatic group, being consistent with the examined single-crystal structures and our closest packing model. However, the disturbance in the πstacking would be enhanced gradually as the increase of linear or branch hydrocarbon in the alkyl chain. General Applicability. With the intention of examining the general applicability of the closest packing model derived from the governing character of dispersion, additional 101 organic crystals comprising various conjugated functionalities and substituents are investigated by ETS−NOCV as well. Molecular pairs extracted from the crystal structures were calculated, and the ΔEints are summarized in Table S13. Along with the seven molecules listed in Table 2, the sampling scope is up to 108. The results reveal that in nearly all the cases, the contribution of dispersion in the ΔEint is larger than 60%, validating that the dominant role of dispersion force is independent of the variation in the conjugated moieties and substituents. It can thus be envisaged that the derived principles in this work should have broad applicability.

Figure 2. Deviation of the pentanyl group from the stable zigzag configuration in the C5H11-NDI crystal.

the crystal structure of C5H11-NDI, the pentanyl group deviates from the stable zigzag configuration to reduce the dimension of lamellar stacking, leading to the formation of a herringbone motif. In contrast, the long or bulky alkyl chains, such as C14H29, C12H25, C6H13, and 1-CH3C6H12, would preclude R-NDI molecules from the herringbone motif and the lamellar motif would be realized instead. Nevertheless, C4H9NDI seems to be an exception. Given that the butyl group is shorter than the pentanyl one, C4H9-NDI is expected to adopt a herringbone motif, which is yet not the experimental observation. In Table S10, the interactions of C4H9-NDI in the lamellar stacking can be up to −7.70 kcal mol−1 in magnitude, which is significantly higher than the others. Figure 3 clearly shows the intercalation of the butyl chains in the lamellar stacking for C4H9-NDI, which not only contracts the lamellarstacking dimension but also enhances the dispersion

Figure 3. Intercalation of the butyl chains in the lamellar stacking in the C4H9-NDI crystal. 18659

DOI: 10.1021/acsomega.8b02713 ACS Omega 2018, 3, 18656−18662

ACS Omega

Article

Table 5. ETS−NOCV Results of the Optimized Adducts between the NDI and Various Molecules in kcal mol−1



entry

molecule

ΔEint

EPauli

1 2 3 4 5 6 7 8

ethane hexane 2-methyl-pentane 3-methyl-pentane 2,2-dimethyl-butane 2,3-dimethyl-butane tetradecane NDI

−4.79 −8.96 −8.45 −9.28 −8.21 −9.06 −13.36 −24.52

6.4 11.65 11.22 12.81 11.36 12.32 17.22 28.25

CONCLUSIONS ETS−NOCV has been performed to study the π- and lamellar stacking for a series of R-NDI crystals. The energy decomposition results indicate that the dispersion is the most important stabilization in the intermolecular interactions. The increase in the molecular contact should in principle strengthen the dispersion. The closest packing model for elucidating the stacking principles for R-NDIs is thus developed. It is rational to see that the R and NDI groups undergo nanoscopic phase separation in the π-stacking. Furthermore, the contribution of R to the formation of πstacking can be significant, challenging the common belief that the π-stacking is predominantly determined by the conjugated moiety. As for the lamellar stacking, the corresponding interactions are much smaller than those resulting from the π-stacking. In order to achieve maximal stabilization, to reduce the length of lamellar stacking and leave more space for the πstacking is a straightforward approach, which would lead to the herringbone motif. In contrast, the long or bulky alkyl chains would prohibit R-NDIs from the herringbone motif, and the lamellar motif would be realized instead. It has to be emphasized that Pauli repulsion would restrain two molecules from establishing too compact stacking and the closest packing model has its limitation. Moreover, the interactions between NDI and the alkyl chain are examined. The NDI molecule is able to form stable adducts with different alkanes, and the interactions between NDI and the alkyl chain can be substantial, especially for the extensive alkyl chain. The results provide rationalization for the use of long and branch side chains to increase the solubility for conjugated molecules and polymers by hindering the π-stacking. Along with R-NDIs, additional organic crystals comprising various conjugated functionalities and substituents are also investigated by ETS−NOCV. The sampling scope is up to 108 conjugated molecules. The dominant role of dispersion force irrespective of the variation in the conjugated moieties and substituents is further confirmed. It is envisaged that the established principles are applicable to other organic semiconductors. The traditional perspective toward the π- and lamellar stacking might be modified, paving the way for ultimate morphological control.

Vele



−2.15 −4.35 −3.93 −4.61 −4.03 −4.40 −7.51 −11.5

(19%) (21%) (20%) (21%) (21%) (21%) (25%) (22%)

Eoi

Edis

−1.9 (17%) −3.28 (16%) −3.19 (16%) −3.59 (16%) −3.31 (17%) −3.48 (16%) −4.88 (16%) −8.31 (16%)

−7.14 (64%) −12.98 (63%) −12.55 (64%) −13.89 (63%) −12.23 (62%) −13.50 (63%) −18.19 (59%) −32.96 (62%)

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsomega.8b02713.



ETS−NOCV results of π- and lamellar stacking in C12H25-NDI, C 6H13-NDI, C5H11-NDI, C4H9-NDI, C2H5-NDI, and 1-CH3C6H12-NDI in kcal mol−1; interacting modes in the lamellar stacking for C14H29NDI, C 12 H 25 -NDI, C 6 H 13 -NDI, C 4 H 9 -NDI, 1CH3C6H12-NDI, C5H11-NDI, and C2H5-NDI; ETS− NOCV results in kcal mol−1 for the molecular pairs in 101 organic crystals; and Cartesian coordinates of the optimized structures (PDF)

AUTHOR INFORMATION

Corresponding Author

*E-mail: yuyinglai@ntu.edu.tw (Y.-Y.L.). ORCID

Yu-Ying Lai: 0000-0002-1921-6923 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS We acknowledge the support of this work by the Ministry of Science and Technology, Taiwan, and National Taiwan University. We are grateful to the National Center for Highperformance Computing for computer time and facilities.



REFERENCES

(1) Wu, J.-S.; Cheng, S.-W.; Cheng, Y.-J.; Hsu, C.-S. Donor− acceptor conjugated polymers based on multifused ladder-type arenes for organic solar cellsDonor-acceptor conjugated polymers based on multifused ladder-type arenes for organic solar cells. Chem. Soc. Rev. 2015, 44, 1113−1154. (2) Hedley, G. J.; Ruseckas, A.; Samuel, I. D. W. Light Harvesting for Organic Photovoltaics. Chem. Rev. 2016, 117, 796−837. (3) Muccini, M. A bright future for organic field-effect transistors. Nat. Mater. 2006, 5, 605. (4) Mas-Torrent, M.; Rovira, C. Role of Molecular Order and SolidState Structure in Organic Field-Effect Transistors. Chem. Rev. 2011, 111, 4833−4856. (5) Torsi, L.; Magliulo, M.; Manoli, K.; Palazzo, G. Organic fieldeffect transistor sensors: a tutorial review. Chem. Soc. Rev. 2013, 42, 8612−8628. (6) Xu, R.-P.; Li, Y.-Q.; Tang, J.-X. Recent advances in flexible organic light-emitting diodes. J. Mater. Chem. C 2016, 4, 9116−9142. (7) Yang, Z.; Mao, Z.; Xie, Z.; Zhang, Y.; Liu, S.; Zhao, J.; Xu, J.; Chi, Z.; Aldred, M. P. Recent advances in organic thermally activated delayed fluorescence materials. Chem. Soc. Rev. 2017, 46, 915−1016.



THEORETICAL METHODS All calculations were performed with ADF 2016.65 BP86 density functional is selected as the computational method. Triple zeta with 1 polarization function (TZP) is the basis set. Grimme’s BJDAMP is used for dispersion correction. Geometry optimization was carried out at the same level of theory mentioned above. The minimal nature of all optimized structures is examined by frequency computation. Cartesian coordinates are listed in the Supporting Information. 18660

DOI: 10.1021/acsomega.8b02713 ACS Omega 2018, 3, 18656−18662

ACS Omega

Article

Substituent effects on aromatic stacking interactions. Org. Biomol. Chem. 2007, 5, 1062−1080. (26) Bloom, J. W. G.; Wheeler, S. E. Taking the Aromaticity out of Aromatic Interactions. Angew. Chem., Int. Ed. 2011, 50, 7847−7849. (27) Wheeler, S. E. Local Nature of Substituent Effects in Stacking Interactions. J. Am. Chem. Soc. 2011, 133, 10262−10274. (28) Bloom, J. W. G.; Raju, R. K.; Wheeler, S. E. Physical Nature of Substituent Effects in XH/π Interactions. J. Chem. Theory Comput. 2012, 8, 3167−3174. (29) Raju, R. K.; Bloom, J. W. G.; Wheeler, S. E. Broad Transferability of Substituent Effects in π-Stacking Interactions Provides New Insights into Their Origin. J. Chem. Theory Comput. 2013, 9, 3479−3490. (30) Wheeler, S. E. Understanding Substituent Effects in Noncovalent Interactions Involving Aromatic Rings. Acc. Chem. Res. 2012, 46, 1029−1038. (31) Hwang, J.; Li, P.; Carroll, W. R.; Smith, M. D.; Pellechia, P. J.; Shimizu, K. D. Additivity of Substituent Effects in Aromatic Stacking Interactions. J. Am. Chem. Soc. 2014, 136, 14060−14067. (32) Parrish, R. M.; Sherrill, C. D. Quantum-Mechanical Evaluation of π−π versus Substituent−π Interactions in π Stacking: Direct Evidence for the Wheeler−Houk Picture. J. Am. Chem. Soc. 2014, 136, 17386−17389. (33) Sinnokrot, M. O.; Sherrill, C. D. Substituent Effects in π−π Interactions: Sandwich and T-Shaped Configurations. J. Am. Chem. Soc. 2004, 126, 7690−7697. (34) Ringer, A. L.; Sinnokrot, M. O.; Lively, R. P.; Sherrill, C. D. The Effect of Multiple Substituents on Sandwich and T-Shaped π−π Interactions. Chem.Eur. J. 2006, 12, 3821−3828. (35) Ringer, A. L.; Sherrill, C. D. Substituent Effects in Sandwich Configurations of Multiply Substituted Benzene Dimers Are Not Solely Governed By Electrostatic Control. J. Am. Chem. Soc. 2009, 131, 4574−4575. (36) Hohenstein, E. G.; Sherrill, C. D. Effects of Heteroatoms on Aromatic π−π Interactions: Benzene−Pyridine and Pyridine Dimer. J. Phys. Chem. A 2009, 113, 878−886. (37) Hohenstein, E. G.; Duan, J.; Sherrill, C. D. Origin of the Surprising Enhancement of Electrostatic Energies by ElectronDonating Substituents in Substituted Sandwich Benzene Dimers. J. Am. Chem. Soc. 2011, 133, 13244−13247. (38) Watt, M.; Hardebeck, L. K. E.; Kirkpatrick, C. C.; Lewis, M. Face-to-Face Arene−Arene Binding Energies: Dominated by Dispersion but Predicted by Electrostatic and Dispersion/Polarizability Substituent Constants. J. Am. Chem. Soc. 2011, 133, 3854− 3862. (39) Sherrill, C. D. Energy Component Analysis of π Interactionsr. Acc. Chem. Res. 2012, 46, 1020−1028. (40) Ehrlich, S.; Moellmann, J.; Grimme, S. Dispersion-Corrected Density Functional Theory for Aromatic Interactions in Complex Systems. Acc. Chem. Res. 2012, 46, 916−926. (41) Stephan, E.; Bettinger, H. F.; Stefan, G. Dispersion-Driven Conformational Isomerism in σ-Bonded Dimers of Larger Acenes. Angew. Chem., Int. Ed. 2013, 52, 10892−10895. (42) Morokuma, K. Why do molecules interact? The origin of electron donor-acceptor complexes, hydrogen bonding and proton affinity. Acc. Chem. Res. 1977, 10, 294−300. (43) Ziegler, T.; Rauk, A. Carbon monoxide, carbon monosulfide, molecular nitrogen, phosphorus trifluoride, and methyl isocyanide as .sigma. donors and .pi. acceptors. A theoretical study by the HartreeFock-Slater transition-state method. Inorg. Chem. 1979, 18, 1755− 1759. (44) Ziegler, T.; Rauk, A. A theoretical study of the ethylene-metal bond in complexes between copper(1+), silver(1+), gold(1+), platinum(0) or platinum(2+) and ethylene, based on the HartreeFock-Slater transition-state method. Inorg. Chem. 1979, 18, 1558− 1565. (45) Mitoraj, M. P.; Michalak, A.; Ziegler, T. A Combined Charge and Energy Decomposition Scheme for Bond Analysis. J. Chem. Theory Comput. 2009, 5, 962−975.

(8) Liu, Y.; Li, C.; Ren, Z.; Yan, S.; Bryce, M. R. All-organic thermally activated delayed fluorescence materials for organic lightemitting diodes. Nat. Rev. Mater. 2018, 3, 18020. (9) Park, S.; Kang, Y. J.; Majd, S. A Review of Patterned Organic Bioelectronic Materials and their Biomedical Applications. Adv. Mater. 2015, 27, 7583−7619. (10) Simon, D. T.; Gabrielsson, E. O.; Tybrandt, K.; Berggren, M. Organic Bioelectronics: Bridging the Signaling Gap between Biology and Technology. Chem. Rev. 2016, 116, 13009−13041. (11) Lee, C.-H.; Lai, Y.-Y.; Cao, F.-Y.; Hsu, J.-Y.; Lin, Z.-L.; Jeng, U.S.; Su, C.-J.; Cheng, Y.-J. Synthesis, molecular and photovoltaic/ transistor properties of heptacyclic ladder-type di(thienobenzo)fluorene-based copolymers. J. Mater. Chem. C 2016, 4, 11427−11435. (12) Lee, C.-H.; Lai, Y.-Y.; Hsu, J.-Y.; Huang, P.-K.; Cheng, Y.-J. Side-chain modulation of dithienofluorene-based copolymers to achieve high field-effect mobilities. Chem. Sci. 2017, 8, 2942−2951. (13) Chiou, D.-Y.; Cao, F.-Y.; Hsu, J.-Y.; Tsai, C.-E.; Lai, Y.-Y.; Jeng, U.-S.; Zhang, J.; Yan, H.; Su, C.-J.; Cheng, Y.-J. Synthesis and sidechain isomeric effect of 4,9-/5,10-dialkylated-β-angular-shaped naphthodithiophenes-based donor−acceptor copolymers for polymer solar cells and field-effect transistors. Polym. Chem. 2017, 8, 2334− 2345. (14) Xiao, Z.; Sun, K.; Subbiah, J.; Qin, T.; Lu, S.; Purushothaman, B.; Jones, D. J.; Holmes, A. B.; Wong, W. W. H. Effect of molecular weight on the properties and organic solar cell device performance of a donor−acceptor conjugated polymer. Polym. Chem. 2015, 6, 2312− 2318. (15) Kang, H.; Uddin, M. A.; Lee, C.; Kim, K.-H.; Nguyen, T. L.; Lee, W.; Li, Y.; Wang, C.; Woo, H. Y.; Kim, B. J. Determining the Role of Polymer Molecular Weight for High-Performance All-Polymer Solar Cells: Its Effect on Polymer Aggregation and Phase Separation. J. Am. Chem. Soc. 2015, 137, 2359−2365. (16) Katsouras, A.; Gasparini, N.; Koulogiannis, C.; Spanos, M.; Ameri, T.; Brabec, C. J.; Chochos, C. L.; Avgeropoulos, A. Systematic Analysis of Polymer Molecular Weight Influence on the Organic Photovoltaic Performance. Macromol. Rapid Commun. 2015, 36, 1778−1797. (17) Hung, K.-E.; Tsai, C.-E.; Chang, S.-L.; Lai, Y.-Y.; Jeng, U.-S.; Cao, F.-Y.; Hsu, C.-S.; Su, C.-J.; Cheng, Y.-J. BispentafluorophenylContaining Additive: Enhancing Efficiency and Morphological Stability of Polymer Solar Cells via Hand-Grabbing-Like Supramolecular Pentafluorophenyl−Fullerene Interactions. ACS Appl. Mater. Interfaces 2017, 9, 43861−43870. (18) Lai, Y.-Y.; Shih, P.-I.; Li, Y.-P.; Tsai, C.-E.; Wu, J.-S.; Cheng, Y.J.; Hsu, C.-S. Interface Engineering to Enhance the Efficiency of Conventional Polymer Solar Cells by Alcohol-/Water-Soluble C60 Materials Doped with Alkali Carbonates. ACS Appl. Mater. Interfaces 2013, 5, 5122−5128. (19) Tsai, C.-E.; Liao, M.-H.; Chen, Y.-L.; Cheng, S.-W.; Lai, Y.-Y.; Cheng, Y.-J.; Hsu, C.-S. Triarylamine-based crosslinked hole-transporting material with an ionic dopant for high-performance PEDOT:PSS-free polymer solar cells. J. Mater. Chem. C 2015, 3, 6158−6165. (20) Yao, Z.-F.; Wang, J.-Y.; Pei, J. Control of π−π Stacking via Crystal Engineering in Organic Conjugated Small Molecule Crystals. Cryst. Growth Des. 2017, 18, 7−15. (21) Hunter, C. A.; Sanders, J. K. M. The nature of .pi.-.pi. interactions. J. Am. Chem. Soc. 1990, 112, 5525−5534. (22) Hunter, C. A.; Lawson, K. R.; Perkins, J.; Urch, C. J. Aromatic interactions. J. Chem. Soc., Perkin Trans. 2 2001, 651−669. (23) Cockroft, S. L.; Hunter, C. A.; Lawson, K. R.; Perkins, J.; Urch, C. J. Electrostatic Control of Aromatic Stacking Interactions. J. Am. Chem. Soc. 2005, 127, 8594−8595. (24) Cockroft, S. L.; Hunter, C. A. Chemical double-mutant cycles: dissecting non-covalent interactions. Chem. Soc. Rev. 2007, 36, 172− 188. (25) Cockroft, S. L.; Perkins, J.; Zonta, C.; Adams, H.; Spey, S. E.; Low, C. M. R.; Vinter, J. G.; Lawson, K. R.; Urch, C. J.; Hunter, C. A. 18661

DOI: 10.1021/acsomega.8b02713 ACS Omega 2018, 3, 18656−18662

ACS Omega

Article

(46) Mitoraj, M. P.; Michalak, A.; Ziegler, T. On the Nature of the Agostic Bond between Metal Centers and β-Hydrogen Atoms in Alkyl Complexes. An Analysis Based on the Extended Transition State Method and the Natural Orbitals for Chemical Valence Scheme (ETS-NOCV). Organometallics 2009, 28, 3727−3733. (47) Mitoraj, M. P.; Michalak, A. σ-Donor and π-Acceptor Properties of Phosphorus Ligands: An Insight from the Natural Orbitals for Chemical Valence. Inorg. Chem. 2010, 49, 578−582. (48) Stefan, G. Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787−1799. (49) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. (50) Stefan, G.; Stephan, E.; Lars, G. Effect of the damping function in dispersion corrected density functional theory. J. Comput. Chem. 2011, 32, 1456−1465. (51) Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB). J. Phys. Chem. Lett. 2014, 5, 1785−1789. (52) Liao, M.-H.; Tsai, C.-E.; Lai, Y.-Y.; Cao, F.-Y.; Wu, J.-S.; Wang, C.-L.; Hsu, C.-S.; Liau, I.; Cheng, Y.-J. Morphological Stabilization by Supramolecular Perfluorophenyl-C60 Interactions Leading to Efficient and Thermally Stable Organic Photovoltaics. Adv. Funct. Mater. 2014, 24, 1418−1429. (53) Grimme, S.; Hansen, A.; Brandenburg, J. G.; Bannwarth, C. Dispersion-Corrected Mean-Field Electronic Structure Methods. Chem. Rev. 2016, 116, 5105−5154. (54) Chen, Y.-L.; Hsu, J.-Y.; Lin, F.-Y.; Lai, Y.-Y.; Chou, H.-C.; Cheng, Y.-J. Synthesis and Isomeric Effects of Ladder-Type Alkylated Terbenzodithiophene Derivatives. J. Org. Chem. 2016, 81, 2534− 2542. (55) Cao, F.-Y.; Lai, Y.-Y.; Chen, Y.-L.; Cheng, Y.-J. Self-assembled tri-, tetra- and penta-ethylene glycols as easy, expedited and universal interfacial cathode-modifiers for inverted polymer solar cells. J. Mater. Chem. A 2016, 4, 8707−8715. (56) Tsai, C.-E.; Yu, R.-H.; Lin, F.-J.; Lai, Y.-Y.; Hsu, J.-Y.; Cheng, S.-W.; Hsu, C.-S.; Cheng, Y.-J. Synthesis of a 4,9-Didodecyl AngularShaped Naphthodiselenophene Building Block To Achieve HighMobility Transistors. Chem. Mater. 2016, 28, 5121−5130. (57) D. E. Lynch, D. G. H. CCDC 230429: Experimental Crystal Structure Determination 2014. (58) Alvey, P. M.; Reczek, J. J.; Lynch, V.; Iverson, B. L. A Systematic Study of Thermochromic Aromatic Donor−Acceptor Materials. J. Org. Chem. 2010, 75, 7682−7690. (59) Shukla, D.; Nelson, S. F.; Freeman, D. C.; Rajeswaran, M.; Ahearn, W. G.; Meyer, D. M.; Carey, J. T. Thin-Film Morphology Control in Naphthalene-Diimide-Based Semiconductors: High Mobility n-Type Semiconductor for Organic Thin-Film Transistors. Chem. Mater. 2008, 20, 7486−7491. (60) Andric, G.; Boas, J. F.; Bond, A. M.; Fallon, G. D.; Ghiggino, K. P.; Hogan, C. F.; Hutchison, J. A.; Lee, M. A.-P.; Langford, S. J.; Pilbrow, J. R.; Troup, G. J.; Woodward, C. P. Spectroscopy of Naphthalene Diimides and Their Anion Radicals. Aust. J. Chem. 2004, 57, 1011−1019. (61) Pandeeswar, M.; Khare, H.; Ramakumar, S.; Govindaraju, T. Biomimetic molecular organization of naphthalene diimide in the solid state: tunable (chiro-) optical, viscoelastic and nanoscale properties. RSC Adv. 2014, 4, 20154−20163. (62) Fang, L.; Zhou, Y.; Yao, Y.-X.; Diao, Y.; Lee, W.-Y.; Appleton, A. L.; Allen, R.; Reinspach, J.; Mannsfeld, S. C. B.; Bao, Z. Side-Chain Engineering of Isoindigo-Containing Conjugated Polymers Using Polystyrene for High-Performance Bulk Heterojunction Solar Cells. Chem. Mater. 2013, 25, 4874−4880. (63) Hwang, Y.-J.; Earmme, T.; Subramaniyan, S.; Jenekhe, S. A. Side chain engineering of n-type conjugated polymer enhances photocurrent and efficiency of all-polymer solar cells. ChemComm 2014, 50, 10801−10804.

(64) Qi, J.; Zhou, X.; Yang, D.; Qiao, W.; Ma, D.; Wang, Z. Y. Optimization of Solubility, Film Morphology and Photodetector Performance by Molecular Side-Chain Engineering of Low-Bandgap Thienothiadiazole-Based Polymers. Adv. Funct. Mater. 2014, 24, 7605−7612. (65) ADF2016, SCM. Theoretical Chemistry; Vrije Universiteit: Amsterdam, The Netherlands. http://www.scm.com.

18662

DOI: 10.1021/acsomega.8b02713 ACS Omega 2018, 3, 18656−18662