Brooke
Kuei
a and
Enrique D.
Gomez
*bc
aDepartment of Materials Science and Engineering, The Pennsylvania State University, University Park, PA 16802, USA
bDepartment of Chemical Engineering, The Pennsylvania State University, University Park, PA 16802, USA. E-mail: edg12@psu.edu
cMaterials Research Institute, The Pennsylvania State University, University Park, PA 16802, USA
First published on 1st August 2016
Conjugated polymers may play an important role in various emerging optoelectronic applications because they combine the chemical versatility of organic molecules and the flexibility, stretchability and toughness of polymers with semiconducting properties. Nevertheless, in order to achieve the full potential of conjugated polymers, a clear description of how their structure, morphology, and macroscopic properties are interrelated is needed. We propose that the starting point for understanding conjugated polymers includes understanding chain conformations and phase behavior. Efforts to predict and measure the persistence length have significantly refined our intuition of the chain stiffness, and have led to predictions of nematic-to-isotropic transitions. Exploring mixing between conjugated polymers and small molecules or other polymers has demonstrated tremendous advancements in attaining the needed properties for various optoelectronic devices. Current efforts continue to refine our knowledge of chain conformations and phase behavior and the factors that influence these properties, thereby providing opportunities for the development of novel optoelectronic materials based on conjugated polymers.
The electronic properties of conjugated polymers and hence their performance in electronic applications are intimately tied to the conformations of polymer chains and how chains are ordered within films.20–23 Theoretical descriptions of polymers in general begin with the chain conformation, or assumptions of the stiffness of the chain. For example, in solution, chain conformations likely impact assembly in the film and chain conformations in the melt (films) affect the formation of liquid crystalline phases, crystallization, and packing of chains. Conjugated polymers are stiffer than common flexible polymers such as polyethylene, but there are only a few measurements of the persistence length in the literature. We dedicate the first section of this article to chain conformations of conjugated polymers.
In the second half of this review, we discuss the phase behavior of conjugated polymers; specifically, we highlight the crystalline and liquid crystalline phases found in homopolymers, and discuss mixing with solvents, fullerene and other polymers. Solvent choice is a common way to tune the properties of polymeric active layers and therefore the performance of organic electronic devices, and mixing with small molecules such as fullerene lead to significantly enhanced optoelectronic properties, yet open questions remain on determining the solubilities of these materials.24–26 For example, the phase behavior in polymer–fullerene active layers strongly affects charge photogeneration efficacy in photovoltaic devices,27,28 and it is hypothesized that an optimum miscibility exists between the components of the active layer to maximize photovoltaic device performance.29,30 Blending of polymers with conjugated and non-conjugated polymers also has the potential to enhance optoelectronic properties;31,32 for example, blending of different kinds of polymers in the active layer has improved the photoluminescence and electroluminescence efficiencies in LED devices.33 Mixtures containing conjugated polymers, however, are complex because they potentially contain semicrystalline and liquid crystalline components, making the description of the phase behavior challenging.
In this review article, we summarize current knowledge of chain conformations and phase behavior of conjugated polymers. In the first section, we focus on how to measure and predict persistence lengths and what factors influence chain conformations. In the second section, we discuss the phase behavior of conjugated polymers, particularly for homopolymers and mixtures of polymers with small molecules and other polymers. We highlight both emerging concepts to describe such systems, as well as experimentally observed phase behavior.
The size of polymer chains is described by the average end-to-end distance , where h is the sum of the vectors representing monomers in a polymer chain. We can model monomer units making up a polymer backbone by vectors ri with fixed lengths and fixed deflection angles θi between them, as shown in Fig. 1a. The vectors rotate by a dihedral angle ϕi as allowed by the dihedral potential. Thus, if the dihedral potential is known, the persistence length can be extracted by integrating over the distribution of dihedral angles;36 this is known as the hindered rotation model.37 A simplified version is the freely rotating chain model, in which the dihedral potential is simply a constant. The Kratky–Porod chain model based on freely rotating chains gives the relationship for worm-like polymers between 〈h2〉, lp, and the contour length L, which is the length of the polymer chain:
〈h2〉 = 2lpL − 2lp2(1 − exp[−L/lp]) | (1) |
Fig. 1 (a) Depiction of the monomer vectors and angles used to construct the freely-rotating and hindered rotation chain models. Monomer units are represented by vectors ri and have fixed lengths and fixed deflection angles θi. They can rotate in a cone with dihedral angle ϕi. Reprinted with permission ref. 36. (b) Average of the end-to-end distance squared (normalized by the persistence length squared) as a function of number of persistence lengths per chain using the freely-rotating worm-like chain model. The rod limit has a slope of 2 and the coil limit a slope of 1. The cross indicates where the chain length is one persistence length. |
The persistence lengths of conjugated polymers are governed by the backbone deflection angles and torsional potentials.36 For example, alternating double and single bonds lead to small backbone deflection angles, potentially leading to stiff chain conformations. Also, five-member rings (such as thiophene and selenophene) and other moieties (i.e. pyrrolopyrrole) lead to small deflection angles, circa 15 degrees, and smaller lp. Polymers with six-member rings that have straight backbones such as polyphenylene vinylene (PPV) have zero deflection angles, such that persistence lengths are due to bond-bending along the backbone and are on the order of tens of nanometers.38 As an example, Fig. 2 shows a chain conformation of poly(3-hexylthiophene-2,5-diyl) (P3HT), with a persistence length of 2.8 nm (7 repeat units) and molecular weight of 16.6 kg mol−1 (100 repeat units). The chain is not rod-like, but cannot adopt a random coil conformation either; instead, P3HT behaves as a worm-like chain. The chain can loop and make significant backbone deflections, but over distances larger than seven repeat units.
Fig. 2 Depiction of a P3HT chain obtained from a dihedral distribution derived using density functional theory. The chain is worm-like, where it can bend significantly over large distances. Reprinted (adapted) with permission from ref. 36. |
The next sub-sections highlight the significant progress made on understanding the chain conformations of conjugated polymers. In addition to providing a table of the available persistence lengths, we discuss methods for prediction of the persistence length and present an overview of the different factors that influence chain conformations.
Polymer | Solvent | Method | Persistence length lp (nm) | Freely rotating chain model |
---|---|---|---|---|
a Freely rotating chain model returns lp > 1 μm due to rod-like backbones. Actual lp should be governed by bond bending. | ||||
P3BT | Nitrobenzene | SANS | 2.7541 | 2.8 |
rra-P3HT | DCB | SANS | 0.9 ± 0.139 | 2.8 |
Toluene | SANS | 1.1 ± 0.439 | ||
rr-P3HT | DCB | SANS | 2.9 ± 0.139 | 2.8 |
Toluene | SANS | 2.7545 | ||
Toluene | Viscosity | 2.5 ± 0.245 | ||
Tetrahydrofuran | Viscosity | 2.2 ± 0.246 | ||
Tetrahydrofuran | Light scattering | 2.6 ± 0.246 | ||
DFT dihedral potential | 4.036 | |||
MD dihedral potential | 3.236 | |||
MD tangent–tangent correlations | 2.847 | |||
P3OT | Tetrahydrofuran | Viscosity | 2.2 ± 0.248 | 2.8 |
Tetrahydrofuran | Light scattering | 2.6 ± 0.548 | ||
P3BT | Nitrobenzene | SANS | 2.7541 | 2.8 |
P3EHT | Toluene | SANS | 2.7545 | 2.8 |
Toluene | Viscosity | 2.5 ± 0.245 | ||
P3DDT | DCB | SANS | 1.6 ± 0.139 | 2.8 |
Toluene | SANS | 1.5 ± 0.139 | ||
PBTTT | Chlorobenzene | Viscosity | 9.049 | 4.1 |
PFTBT (APFO-3) | MD dihedral potential | 5.936 | 6.6 | |
PFO | Toluene | Light scattering | 7.052 | 10.9 nm |
Tetrahydrofuran | Light scattering | 8.653 | ||
cis-PBO | MSA | Light scattering | 20–3050 | >1 μma |
MSA | Viscosity | 20–2550 | ||
MD | 32.551 | |||
PPP | Chloroform | Light scattering | 2854 | >1 μma |
Toluene | Viscosity | 15.655 | ||
MD | 2251 | |||
trans-PBZT | MD | 21.551 | >1 μma | |
MEH-PPV | Tetrahydrofuran | Viscosity | 2.2 ± 0.246 | >1 μma |
Tetrahydrofuran | Light scattering | 2.6 ± 0.22 | ||
DFT dihedral potential | 4.036 | |||
BEH-PPV | p-Xylene | Light scattering | 11.02 | >1 μma |
DP10-PPV | MD tangent–tangent correlations | 2.847 | >1 μma | |
Chloroform | Light scattering | 6.556 | ||
Chloroform | MD (coarse-grained) | 25 ± 0.557 | ||
BCHA-PPV | p-Xylene | Light scattering | 402 | >1 μma |
PPE | Tetrahydrofuran | Light scattering | 15.058 | >1 μma |
Experimental techniques that can measure the persistence length include small angle neutron scattering (SANS), light scattering, and viscometry. For persistence length measurements of conjugated polymers using SANS, scattering contrast is obtained from either deuterated chains within a nondeuterated melt or between solvent and polymer chains where one component is deuterated.39 The resultant scattering curves can be fitted to several models of scattering from polymer chains,39–41 such as the Debye model,42 the Des Cloizeaux model,43 and the Sharp and Bloomfield model, to determine the radius of gyration.44 The persistence length can then be extracted if the molecular weight of the polymer is known.
Light scattering experiments can also measure the size of polymer chains. Because the wavelength is much larger than that of neutrons and X-rays used in scattering experiments, the accessible range is limited to small scattering vectors. Furthermore, light scattering experiments to determine the size of the chain require polymer chains dispersed in a solvent or other media to generate the needed contrast. Using the Zimm equation, the radius of gyration Rg can be extracted from light scattering data and the persistence length can be inferred.2 Similarly, viscometry of polymer solutions can be used to determine the viscosity-averaged radius of gyration from the Flory–Fox equation.46,59 In both viscometry and light scattering experiments, it is crucial that chains are fully dispersed in the solvent. An assumption is that chains are not aggregated or have crystallized, as that can lead to overestimation of the stiffness.51 A common artifact is then the appearance of different persistence lengths at different polymer concentrations in solution.56
In principle, SANS is the only technique that allows for measurements of chain dimensions in the melt, by deuterating only a small portion of the sample, but no reports of such measurements exist for conjugated polymers. In all of the above techniques, the molecular weight must be known accurately to extract the persistence length. Even though both light scattering and viscometry yield estimates of the molecular weight, these techniques are most useful in combination with gel permeation chromatography to extract a robust measurement of the molecular weight with a combination of refractive index, viscometry and light scattering detectors.37 Alternatively, the molecular weight can be determined using 1H NMR by counting chain ends or using matrix-assisted laser desorption/ionization with a time of flight mass spectrometer (MALDI-TOF) to measure the molecular weights of the individual chains, but these techniques are only accurate at low molecular weights, typically less than 10 kg mol−1.
Table 1 shows the persistence length for various conjugated polymers, in order of the persistence lengths predicted from the freely rotating chain model. Despite the limited set of values currently available, it is apparent that most conjugated polymers are semiflexible. If the smallest persistence length is about that of P3HT (3 nm), that is significantly larger than most repeat units even for push–pull copolymers, and is larger than the persistence length for polyethyelene (lp of 0.7 nm).60 Thus, conjugated polymers span a range in chain stiffness from nearly that of polyethylene to the persistence lengths of liquid-crystalline polyesters (10–20 nm)61 and DNA (approximately 50 nm).62,63
Table 1 also highlights some limitations in our knowledge of the persistence lengths of conjugated polymers. Although all of the persistence lengths of poly(3-alkylthiophene)s should be about 3 nm, experimental measurements vary from about 1 nm for regiorandom P3HT (rra-P3HT) to 3 nm for regioregular P3HT (rr-P3HT). Similarly, for poly(phenylene vinylene)s (PPVs) with various side chains, the persistence length should be on the order of tens of nanometers, yet experimental values can be as low as 2 nm. One possibility is that chemical defects cause significant deflections of the chain and therefore reduce the persistence length; for example, for PPV such defects are likely present and consequently kink the chain.20,64–67
If the dihedral potentials are flat enough, the freely rotating chain model can be used to predict the persistence length from the backbone deflection angles. The barriers for torsion of thiophene backbones are within a few kT,36,72 not too dissimilar than for polyethylene. Recipes for applying the freely-rotating worm-like chain model analytically for alternating copolymers with different monomers and multiple types of dihedral angles analytically are available.36 Unfortunately, comparing the values for the persistence length obtained from experiment and simulations to the values predicted from the freely rotating chain model shown in Table 1 is challenging, because variability among the experimental values suggests further measurements are warranted. Nevertheless, it is not unusual for the freely rotating chain model to provide a good estimate of the chain stiffness, as it appears for P3HT and PFTBT, for example.
Refinements beyond the freely-rotating chain model are needed if the dihedral potentials include two (or more) torsional angles that are important. For example, the rotational isomeric state model can be used to predict the persistence length from the Boltzmann-averaged populations of distinct dihedral angles.73 More generally, the hindered rotation model can take the entire dihedral potential to generate a distribution of dihedral angles and therefore chain conformations. This approach can be done numerically, where an average is taken over many backbone conformations, with each conformation corresponding to a set of dihedral angles that are randomly selected from a dihedral distribution.36 For each set of dihedral angles, a backbone conformation is generated by rotating it from a trans conformation to the selected dihedral angle. For each backbone conformation, the tangent–tangent correlation function can be calculated, from which the persistence length is extracted. This approach relies on dihedral potentials as an input, which can be derived from DFT or MD simulations, and can be generalized for polymers that contain multiple types of monomers and different dihedral potentials.
Crystallization can lead to significantly more planar backbone conformations. When conjugated polymers crystallize, clear signatures of strong coupling between chains is apparent in the absorption spectra.74 Furthermore, detailed analysis of the vibronic peaks in polymers such as P3HT suggest an enhancement in chain planarity;75–78 thiophene rings in the unit cell of P3HT are tilted by about 5 degrees,68 and the dihedral distribution is likely within 10% of this angle. Thus, chain conformations are likely extended when compared to the amorphous phase, even though chain folding, or the formation of loops, can occur.79–81 Nevertheless, the planar conformation can have a larger backbone deflection than freely rotating chains, potentially leading to a drop in the persistence length with aggregation or crystallization. In contrast, the radius of gyration of non-conjugated polymers such as polyethylene remain invariant during crystallization as apparent from neutron scattering studies.82–84 Such experiments have not been explored for conjugated polymers.
Chemical defects can kink the chain, decrease the persistence length and disrupt chain packing; an example is the aforementioned PPV derivatives with relatively short persistence lengths (Table 1). The role of various types of defects on chain conformations have been systematically explored using single molecule spectroscopy and MD simulations of MEH-PPV chains (Fig. 4).85 Linear para-terphenyl defects maintain the linearity of the PPV backbone, while ortho-terphenyl defects, even at low concentration, kink the chain and decrease the persistence length. The effect of saturated defects, such as biphenylethane, is concentration dependent; low inclusions do not affect conformations much other than by providing a “hinge” to allow the polymer to fold unto itself. High amounts of inclusions result in a wide distribution of conformations, where both rod-like but folded structures and more random worm-like conformations with shorter persistence lengths are apparent. Similarly, systematically introducing flexible linkers in the PPV backbone makes the chain more flexible.86 These results demonstrate the potential for synthetic strategies based on incorporating defects to tune the chain stiffness and conformations.
Fig. 4 Two representative chain conformations (right and left) for 15 repeat unit oligomers of (A) MEH-PPV, (B) ∼33% para-terphenyl MEH-PPV, (C) ∼33% ortho-terphenyl MEH-PPV, and (D) ∼33% saturated MEH-PPV. Different types of chemical defects perturb chain conformations to different degrees. Results from replica exchange molecular dynamics in the gas phase. Reprinted with permission from ref. 85. |
Side chains or groups can also affect the chain stiffness by affecting dihedral potentials or introducing steric constraints.23,39,87–89 For example, bulky side groups can prevent planar conformations due to sterics and therefore kink the chains. An example is regiorandom P3HT (rra-P3HT), where the head-to-head addition of the side chains could prevent planar conformations and therefore reduce lp. Indeed, SANS measurements place the persistence length of rra-P3HT at about 1 nm,39 and more compact conformations are observed in MD simulations,88 but spacing the side chains to every other ring (such as in regiorandom poly(3-hexyl-2,5-thienylene vinylene)) recovers the stiffer chain conformations.23 The fact that the freely rotating chain model works well for rr-P3HT (see Table 1) suggests that the dihedral potentials must be severely disrupted for rra-P3HT to have a significantly lower persistence length.
Overall, the role of side chains on the persistence length of conjugated polymers remains unclear. SANS experiments on saturated polymers poly-1-octene and poly-1-decene with different numbers of CH2 units in their side groups show that the largest side group leads to the largest persistence length,89 and the authors have argued that a similar dependence might be found in conjugated polymers.90 Nevertheless, atomistic simulations suggest that the longer side chains of poly(2,3-diphenylphenylene vinylene) with decyl side chains (DP10-PPV) do not affect chain packing and persistence lengths in a melt of oligomers as compared to PPV with hexyl side chains (DP6-PPV),57 even though the persistence lengths of PPVs appear to vary with side chains (Table 1). The discrepancy is perhaps due to chemical defects in the PPV chain, as discussed earlier. P3DDT also has a significantly lower persistence length than P3HT in solution, at 1 nm versus 3 nm.39 Perhaps the longer side chains can contribute significant steric effects on the dihedral potentials, although it is challenging to imagine such large contributions.36,72,91
Non-bonding interactions can also be exploited to control chain conformations.92–101 In principle, organic heterocylic groups may be able to use attractive interactions to “lock-in” specific chain conformations, such as highly planar backbones. Detailed calculations have assessed the strength of various non-bonding interactions, and have suggested promising targets for stabilization of chain conformations. For example, as shown in Table 2, CH⋯N and CH⋯O interactions are about three to four times that of kT at room temperature (0.593 kcal mol−1). Other interactions are weaker, such as CH⋯F, while O⋯S, N⋯S, and F⋯S have interaction energies of about kT at room temperature or less, and are not expected to perturb molecular conformations.100
Interaction energies (kcal mol−1) | |
---|---|
CH–N | 2.20 |
CH–O | 1.86 |
CH–F | 0.94 |
CH–S | 0.74 |
S–S | 0.72 |
O–S | 0.51 |
N–S | 0.46 |
F–S | 0.44 |
N–F | Nonbinding |
O–F | Nonbinding |
O–N | Nonbinding |
Interestingly, many polymers in the active layer of high performance thin-film transistors or photovoltaics are very stiff. Table 3 shows a few examples of high-performance conjugated polymers and their corresponding persistence lengths, as calculated from the freely rotating chain model.36 PTB7 and PffBT4T-20D are used as donors in organic photovoltaic devices with power conversion efficiencies beyond 10%,107,108 and DPP-DTT, C16IDT-BT, and PDVT-10 are used as the active layer in thin-film transistors with charge mobilities greater than 1 cm2 V−1 s−1.109–111 Even though it is unclear how the persistence length matters, the large persistence lengths of these materials, at 8 nm or greater, or more than twice that of P3HT, suggests it may be important for optoelectronic properties.
Inducing a liquid-crystalline phase prior to crystallization has been shown to improve order and enhance charge mobilities. For example, a semiconducting polythiophene was designed to assemble into large crystalline domains on crystallization from a liquid-crystal phase, resulting in a system with close intermolecular π–π distances that led to 0.2–0.6 cm2 V−1 s−1 field-effect mobilities.102 It is hypothesized that the liquid-crystalline phase facilitates the alignment of chains and that this order is maintained during crystallization. As a consequence, regions between crystalline domains after this method of processing are likely more organized than amorphous regions of a material crystallized from an isotropic melt.
As alluded to in Section 2.4, the alignment of rigid-rod chains in liquid crystals can improve charge transport. For instance, polyfluorene copolymers in a nematic glassy state with chains aligned in the direction parallel to the transport direction in thin-film transistors (TFTs) through the use of an alignment layer showed enhanced charge carrier mobilities and good operating stability.104 Similarly, improved mobility was observed in chains that were aligned using three molecular design components: concentration-induced chain planarization, a tetrahedral carbon linker with out-of-plane bonding, and bulk side chains that prevent side-chain interaction.127 Thus, predicting the propensity for liquid crystalline phases in conjugated polymers is important.
As a starting point, we can consider the formation of liquid crystalline phases by rigid rods. Onsager128 showed that the isotropic to nematic transition in rigid rods occurs when the volume fraction of rods is high enough that a random placement and orientation of rods would lead to approximately one collision per rod. As the isotropic phase becomes less entropically favorable at high concentrations, the rods will align with each other and form a nematic phase. In semiflexible chains, the critical volume fraction at which the isotropic to nematic transition is predicted to occur is influenced largely by d/lp, where d is the diameter of a rod-like chain and lp is the persistence length.129,130 Although this scaling has been verified for polymers with modest persistence length, recent computer simulations suggest a different scaling of (d/lp)2/3 when lp is large; the origin of such scaling is currently unknown.131
As temperature increases, entropy dominates and order gives way to amorphous phases. The nematic-to-isotropic transition temperature is often described in terms of Maier–Saupe theory, where the Maier–Saupe interaction, or nematic coupling parameter, describes the propensity for rods or stiff polymers to form nematic phases.132–137 A combination of self-consistent field theory (SCFT) and MD simulations was used to calculate the nematic coupling constant of P3HT and the isotropic-to-nematic transition temperature TIN as a function of chain length n.47 The combination of these predictions and experimentally determined crystal melting temperatures69,138 give a comprehensive picture of the phase behavior for P3HT.139
An example of a framework to generalize the behavior of semiconducting polymers that have both liquid crystalline and crystalline phases can be found in Fig. 5.114 The isotropic–nematic phase boundary comes from the computational work of Zhang et al. mentioned above, and the equilibrium melting temperature line is confirmed with values reported in the literature.69,116 The nematic phase is most prevalent at low N, where Tm is about 20–60 degrees below TIN. Fig. 5 also shows that regiodefects decrease Tm and potentially increases the window for the nematic phase.114 Altogether, considering the effects of regioregularity and N help explain seemingly contradictory results in the literature, where nematic phases are visible only at low molecular weights and for regiodefective materials.
Fig. 5 Phase diagram for P3HT. TIN denotes the isotropic–nematic phase boundary, and the dashed curves denote the melting temperatures as a function of N and the regioregularity p. Reproduced with permission from ref. 114. |
Despite tremendous efforts and progress, the crystallization of conjugated polymers is not yet fully understood. So far, the structures of only a few conjugated polymers have been determined, such as for P3HT,140 (poly(3-2,5-dioctylphenyl)thiophene) (PDOPT),141 and poly(2,5-bis(3-tetradecylthiophene-2-yl)thieno[3,2-b]thiophene) (PBTTT-C14).142 Early X-ray diffraction (XRD) studies showed that, in general, poly(3-alkylthiophene)s (P3ATs) exhibit a well-organized lamellar structure consisting of stacks of planar thiophene main chains that are uniformly spaced by alkyl side chains. Structural analysis revealed that this organization occurs even in the presence of local side-chain disorder, an observation that suggests that a secondary structural ordering occurs at the alkyl-chain interface between alkyl-thiophene layers; this process could be related to the liquid crystalline behavior of P3ATs.68,140 More recently, the crystal structure of P3HT has been refined through a combination of X-ray diffraction, solid state NMR spectroscopy, and density functional theory calculations.143 The unit cell of P3HT is monoclinic with a space group P21/c and two chains per cell (a = 1.62 nm, b = 0.76 nm, c = 0.77 nm, γ = 87 degrees). This approach is promising for the determination of the crystal structure of conjugated polymers.
In principle, the degree of crystallinity should affect charge mobilities within conjugated polymers. Strong coupling within crystalline regions enhances delocalization of charge carriers, such that transport is enhanced when compared to amorphous regions. As crystallization progresses, mobilities in thin-film transistors can increase by a factor of about 50; this is likely a result of the formation of a crystalline network within an amorphous phase in the active layer.144 Furthermore, accessing different polymorphs with tighter π–π stacking in the active layer can lead to enhancements of the field-effect mobility in thin-film transistors.145 But, introducing unsubstituted thiophene rings to P3HT can enhance coupling between chains and change mobilities by an order of magnitude, despite a drop in the crystallinity.146 Charge mobilities in thin-film transistors are not always correlated with the crystallinity of the active layer.147,148
The crystallization of P3HT is often limited by the diffusion of the chains, such that low molecular weight chains crystallize more readily than high molecular weight P3HT. Nevertheless, charge mobilities increase with molecular weight.78,147,148 The leading hypothesis to explain this dependence is that charge transport is limited by interconnectivity between crystalline domains, such as through tie chains. Thus, varying the crystallization kinetics can lead to differences in tie chain densities, such that faster crystallization leads to more tie chains in the active layer and higher charge mobilities in thin-film transistors.149 Overall, the lack of correlation between the crystallinity and crystal texture (crystal orientation) in the active layer with charge mobilities in thin-film transistors suggests a one-dimensional character to transport,110 where occasional strong coupling between chains is sufficient to achieve high charge mobilities, to date up to 10 cm2 V−1 s−1 in transistors with polymer semiconductors as the active layer.150
State-of-the-art solution processed organic solar cells rely on mixtures of conjugated polymers and fullerene as the active layer.113,152–154 If the polymer and fullerene are strongly immiscible, the two components phase separate, leading to structures with length scales larger than hundreds of nanometers107,155 that are too large for efficient exciton harvesting due to exciton diffusion being limited to about 10 nm.156,157 If phase separation is suppressed through either miscibility or a compatibilizing agent, then the system relies on crystallization of one or both of the components to create morphologies and interfaces amenable for charge photogeneration.
Typically, both the conjugated polymer donor and fullerene can crystallize, such that the equilibrium state is likely separate fully crystalline phases. Polymers, however, seldom fully crystallize, making amorphous phases ubiquitous. Furthermore, the crystallization of fullerene is also retarded in mixtures with conjugated polymers.162,163 As a consequence, if the crystallization of the polymer component occurs first, then the microstructure is set by the crystallization motif, as shown in Fig. 6. Polymers that crystallize into fibril-like structures with nanometer dimensions create nanoscale morphologies amenable for exciton dissociation and charge extraction,29,158,160,161,164–168 and polymers that crystallize into plate-like structures can allow the domain sizes to grow with processing.29
Fig. 6 Sulfur elemental maps generated from core-loss energy-filtered transmission electron microscopy demonstrating that polymer crystallization dominates the morphology in donor–acceptor mixtures used as the active layer of solar cells. (a) 1:1 rr-P3HT/PCBM films as a function of annealing temperature. Reproduced with permission from ref. 158. Copyright 2011 American Chemical Society. (b) 1:4 poly[2,5-bis(3-hexadecylthiophen-2-yl)thieno[3,2-b]thiophene]/phenyl-C71-butyric acid methyl ester (PC71BM) films, as a function of annealing time and temperature. Adapted from ref. 159 with permission from The Royal Society of Chemistry. (c) 1:1 poly[(4,4′-bis(2-ethylhexyl)dithieno[3,2-b:2′,3′-d]germole)-2,6-diyl-alt-(2,1,3-benzothiadiazole)-4,7-diyl]/PC71BM films as-cast and annealed. Reproduced with permission from ref. 160. (d) 1:2.5 rr-P3HT:N,N′-di((thiophen-2-yl)methyl)-2,6-di(N-cyclohexylamino)-1,4,5,8-naphthalenetetracarboxydiimide films with and without diiodooctane (DIO). Reprinted from ref. 161 with permission from Elsevier. |
Crystallization leads to pure polymer phases, unless the polymer co-crystallizes with fullerene.169–171 The remaining amorphous phase is then enriched in fullerene, such that phase separation driven by limited solubility can then occur. In this context, describing liquid–liquid phase separation becomes relevant, to predict what will occur in the mixed phases between polymer and fullerene. How to treat the miscibility within the non-crystalline domains remains an open challenge.
Perhaps the first question is how to treat fullerene in polymer/fullerene mixtures. Commonly used fullerene derivatives such as [6,6]-phenyl-C61-butyric acid methyl ester (PCBM) have been assumed to behave as nanoparticles, given that they have a diameter of approximately 1 nm.172–174 Others have assumed PCBM is small enough to be treated as a solvent molecule.158,175,176 If particles (or molecules) are large enough, phase separation in mixtures can be driven purely entropically, because the nanoparticle surfaces limit chain conformations and therefore promote separation of the components.
McGarrity and coworkers177 have examined the critical size needed for nanoparticle effects to be important using simulations. In a polymer–nanoparticle mixture with only neutral interactions, they examined the presence of surface segregation as a function of the size of the nanoparticle. They found that when the diameter of the nanoparticle is 1.6 times larger than the statistical segment length of the polymer, the loss of configurational entropy in the polymer becomes important and nanoparticles segregate to the surface entropically (Fig. 7). P3HT has one of the smallest monomer lengths of all conjugated polymers and PC71BM is one of the largest fullerene molecules commonly employed. The ratio of the cube root of the monomer volumes of P3HT and PC71BM (239 and 1080 Å3, respectively178) is about 1.7, right at the onset of entropically-driven phase separation. Thus, although fullerene likely behaves as a small molecule in most mixtures with conjugated polymers, more work is needed to determine the importance of nanoparticle effects in selected systems, such as in P3HT/PC71BM mixtures.
Fig. 7 Excess surface energy for a polymer/nanoparticle composite as a function of the ratio of nanoparticle size to the statistical segment length of the polymer (σn/σp). If nanoparticles are larger than σn/σp of 1.6, surface segregation occurs (regions B and C), suggesting entropic effects can drive phase separation. Reprinted with permission from ref. 177. |
The miscibility of liquid mixtures of small molecules and polymers can be described using Flory–Huggins theory,37,60,179,180 which describes the phase behavior in terms of the Flory–Huggins interactions parameter χ. Essentially, each lattice site is taken to be the volume of one solvent molecule and we assume that each polymer occupies N lattice sites (number of statistical segment lengths). As a consequence, the statistical segment size is calculated at a reference volume of the solvent. For volume fractions φ1 and φ2, we obtain an expression for the free energy of mixing per site of a polymer solution:
(2) |
Flory–Huggins can be useful to predict the behavior of mixtures of conjugated polymers and fullerene. If crystallization of fullerene is suppressed, as is often in PCBM and PCBM derivatives, and the polymer crystallizes, leaving behind some amorphous chains enriched in fullerene content, then one can predict whether phase separation will occur if χ is known. Fig. 8 shows the phase behavior predicted from Flory–Huggins for a polymer mixture of P3HT of 50 kg mol−1 and PC61BM. The phase diagram is clearly asymmetric with composition, where the critical volume fraction is fullerene-rich. The asymmetry is a direct result of the large differences in molecular sizes between polymers and small molecules, leading to a smaller contribution to the entropy of mixing by the polymer, as seen in eqn (2). It may be challenging to dissolve even a small amount of polymer in some solvent (or fullerene), while a small amount of small molecules can always infiltrate into amorphous polymer phases. Thus, when phase separation in polymer/fullerene mixtures occurs, fullerene phases will be nearly pure while amorphous polymer phases will contain a significant amount of fullerene unless χ is very high.
χ at the critical point goes to 0.5 as N gets large, such that for systems where χ is below 0.5 the components are miscible at all compositions. As we discuss in Section 3.4, Kozub et al. measured the Flory–Huggins χ parameter to be 0.86, and demonstrated that phase separation in mixtures of amorphous regiorandom P3HT and PC61BM roughly follow the predictions by Flory–Huggins.158,181
Because miscibility governs phase separation, it also affects charge transport.182 Electron charge mobilities through fullerene domains are suppressed in miscible systems, because the fullerenes are dispersed. Transport through these systems is governed by percolation, such that the mobility decreases with decreasing concentration of fullerene as a power law with an exponent of two.181 If the polymer and fullerene phase separate, then charge mobilities decrease more slowly with dilution, linearly with concentration or volume fraction. As a consequence, an optimum miscibility exists for the active layer in organic solar cells. Strongly immiscible systems lead to large scale phase separation and prevent the polymer donor from establishing a nanoscale morphology. Highly miscible systems suffer from poor electron mobilities in the mixed phase between amorphous polymer and fullerene.181 Treat and coworkers demonstrated this concept experimentally, where they optimized the performance of devices with poly(3-hexylselenophene-2,5-diyl) (P3HS)/fullerene mixtures in the active layer by tuning the miscibility with P3HS using various fullerene derivatives, and they found that indeed an optimum polymer–fullerene miscibility exists for application in the active layer of solar cell devices.30 As a consequence, polymer–fullerene phase behavior in the active layer can have a significant impact on device performance.183–187
Despite its utility, there are limitations to Flory–Huggins.37,188 When crystals are present, the chemical potential of the liquid solution has to be compared to that of the crystalline phases.151 Furthermore, the mean-field assumption within Flory–Huggins fails for dilute polymer solutions, because chain connectivity inevitably causes the concentration of solute to be non-uniform throughout the solution. Also, Flory–Huggins assumes no volume change occurs upon mixing. Other problems lie in lumping all interactions into a single χ parameter, which often requires a temperature-independent contribution to χ; χ often has limited applicability to systems with highly polar groups or specific interactions.37 The role of backbone stiffness on polymer miscibility has also been explored.136,189–191
In mixtures where the microstructure can evolve with time, time temperature transformation (TTT) diagrams are useful for describing the microstructure and microstructural evolution.192–196 Although such diagrams do not currently exist for conjugated polymers or mixtures containing conjugated polymers, we propose a rough sketch of what a time temperature transformation diagram might look like for P3HT/PCBM in Fig. 9. Crystallization is delayed as either the glass transition or melting temperature for each component is approached. Furthermore, if the two materials are immiscible, phase separation is expected to occur faster at higher temperatures due to enhanced diffusivities. Clearly, the interplay between phase separation and crystallization can significantly complicate TTT diagrams, but such effects are not included in our schematic shown in Fig. 9. Nevertheless, we propose that generating TTT diagrams would be useful for understanding the morphological evolution in polymer/fullerene mixtures.
Fig. 9 Rough sketch of a proposed time temperature transformation diagram for P3HT and PCBM. The melting and glass transition temperatures are taken from various studies.154,178,197–201 The red and blue curves outline the boundary of P3HT and PCBM crystallization, while the purple curve denotes the onset of liquid–liquid phase separation. |
Blends of conjugated polymers have demonstrated unique properties. A few decades ago, blends of poly(phenylene-vinylene) (PPV) with an insulating polymer, poly(ethylene oxide) were explored.213 The conductive properties of the blend depend on the compositions and processing conditions, but a composition of greater than 60 wt% PPV gives flexibility, toughness, and the highest conductivity. Others214 made blends of polythiophene (PTh) and polypyrrole (PPy) with insulators polystyrene (PS) and polycarbonate (PC). Interestingly, PTh/PS and PPy/PS blends have high conductivities, but PPy/PC conductivities are lower, likely due to hydrogen bonding between polypyrrole and polycarbonate. Blends of donor and acceptor conjugated polymers are used as the active layer of organic solar cells achieving nearly 8% power conversion efficiencies.215–222 As expected, most of these polymer blends phase separate, as exemplified by the scanning transmission X-ray microscopy (STXM) images shown in Fig. 10 for mixtures of poly(9,9′-dioctylfluorene-cobis(N,N′-(4,butylphenyl))bis(N,N′-phenyl-1,4-phenylene)-diamine) (PFB) and poly(9,9′-dioctylfluorene-co-benzothiadiazole) (F8BT).223
Fig. 10 Composition maps of PFB/F8BT blends generated from STXM experiments. Coarse phase separation is apparent. From ref. 223. |
Various approaches have been demonstrated to enhance mixing of semiflexible polymers. Hydrogen bonding has been employed as specific interactions to promote miscibility between mixtures of semi-flexible and flexible polymers.224,225 The addition of flexible side chains to stiff polymers has also been shown to help form a miscible polymer–polymer blend of two stiff polymers, polyester and polyimide.226,227 In some cases, the flexible side chains can induce cocrystallization of the polymers at a molecular level.203,204 A common way to compatibilize immiscible polymer blends is to include a surfactant.227–233 This approach has only been demonstrated for a few conjugated polymer blends. For example, by adding 3-hexylthiophene to mixtures of P3HT and poly[2,7-(9,9-didodecylfluorene)-alt-5,5-(4′,7′-bis(2-thienyl)-2′,1′,3-benzothiadiazole)] (PF12TBT), phase separation can be suppressed such that a decrease in the domain size is apparent.
Another approach to mix conjugated polymers is to rely on alkyl side chains to promote miscibility without decreasing conductivity. For example, a miscibility window was identified for a blend of poly(3-octylthiophene) (P3OT) with poly(ethylene-co-vinylacetate) (EVA).234 The interaction between P3OT and EVA that resulted in miscibility also induced a conformational change in the P3OT. The role of the side chains on mixing has been systematically examined by mapping out the phase diagrams of P3ATs with four different n-alkyl side chains lengths (n = 4, 6, 8, 12) blended with PS. A non-monotonic relationship between the alkyl side chain length and miscibility was found, such that from n = 4 (P3BT) to n = 6 (P3HT), the miscibility in the blends decreases, but from n = 8 to n = 12 the miscibility increases. The phase diagrams of these blends are shown in Fig. 11.31
Fig. 11 Predicted (solid lines) and experimentally determined (data points) binodal curves for (a) P3BT/PS, (b) P3HT/PS, (c) P3OT/PS, and (d) P3DDT/PS. Reprinted from ref. 31. |
Beyond polymer blends, conjugated block copolymers provide a strategy to control the microstructure using the chemical architecture.235–237 The mesostructures in which block copolymers self-assemble are on the length scale of tens of nanometers, providing an opportunity to tune the microstructure in the active layer of organic electronic devices. Furthermore, covalently bonded donor and acceptor blocks may offer a means to manipulate the molecular structure such that donor–acceptor interfaces can be tuned and used to study charge and energy transfer processes.238–243 Incorporating the wormlike chain model into self-consistent field theory calculation of block copolymers have demonstrated that the chain stiffness affects the stability of various diblock copolymer mesostructures, such that the phase space for lamellae grows with stiffness over other typical morphologies such as cylinders and spheres.191 We refer readers to recent reviews on conjugated block copolymers for more in-depth discussions.6,244–246
Another approach to estimate χ is through Hildebrand solubility parameters (δi) introduced by Hildebrand and Scott250,251 that are defined as the square root of the cohesive energy density. The Flory–Huggins interaction parameter χ can then be estimated from the Hildebrand solubility parameters of the two components:
(3) |
(4) |
A perturbation by Hansen to the Hildebrand solubility parameters is to split the cohesion energy in terms of three components: dispersion forces (atomic), permanent dipole–permanent dipole forces (molecular), and hydrogen bonding. This adds more terms to the previously stated definition of the solubility parameter such that the square of the total (Hildebrand) solubility parameter is the sum of the squares of the Hansen dispersive, polar, and hydrogen bonding solubility parameters.254,255 Using this approach, the solubility parameters of P3HT, poly[2,6-(4,4-bis-(2-ethylhexyl)-4H-cyclopenta[2,1-b;3,4-b′]dithiophene)-alt-4,7(2,1,3-benzothiadiazole)] (PCPDTBT), and PCBM in various solvents have been estimated.25,256 Furthermore, χ for P3HT/PCBM was extracted using mixtures of solvents and non-solvents to track the transition from solubility to non-solubility.26
Monte Carlo and MD simulations have also been used to estimate χ between P3ATs.31 These estimates do not include entropic contributions to χ, as only the pair-wise interaction energies and coordination numbers between 3-alkylthiophene monomers and styrene dimers are considered. Nevertheless, the values of χ are reasonable when applied in Flory–Huggins and compared to experimental phase diagrams.
Table 4 compares estimates from various techniques of the Flory–Huggins interaction parameter χ for mixtures containing conjugated polymers. Discrepancies exist between values reported for the same mixtures, even by a factor of two or more. The majority of χ values for polymer/small molecule mixtures are high, many above the critical value for phase separation in polymer solutions (χ = 0.5). All of the values reported for polymer/polymer blends are high as well, between 0.2 and 2; for comparison, χ between polystyrene and polyisoprene is about 0.05 at 140 °C. Thus, because phase separation occurs for χN > 2, blends should phase separate even at small N.
Polymer–solvent (and polymer–fullene) interaction parameters χi-solvent | |||||||
---|---|---|---|---|---|---|---|
i | χ i-chloroform | χ i-THF | χ i-cyclohexane | χ i-water | χ i-chlorobenzene | χ i-PC61BM | Temperature (°C) |
rr-P3HT | ∼0.2–0.4257 | 150 | |||||
rr-P3HT | 0.92 ± 0.1258 | 150 | |||||
rr-P3HT | 0.86 ± 0.09158 | ∼210 | |||||
rr-P3HT | 0.99 ± 0.0124 | 1.04 ± 0.0424 | 0.078256 | 21 | |||
rr-P3HT | 0.27253 | 0.51253 | 1.02253 | 3.86253 | 30 | ||
rr-P3DDT | 0.54253 | 0.75253 | 1.02253 | 3.86253 | 30 | ||
rra-P3BT | −0.37253 | −0.17253 | −0.17253 | 3.21253 | 30 | ||
rra-P3DDT | −0.69253 | −0.48253 | −0.32253 | 2.09253 | 30 | ||
rra-P3HT | 0.37253 | 0.41253 | 30 | ||||
PF12TBT | 0.063256 | 21 |
Polymer–polymer interaction parameters (reference volume of 0.1 nm3) | ||
---|---|---|
Blend | χ | Temperature (K) |
P3BT/PS (without annealing) | 0.4831 | 423 |
P3BT/PS (with annealing) | 0.3231 | 423 |
P3BT/PS | −0.0181 + (32.8 − 13.7 ϕP3BT)/T31 | — |
P3HT/PS | 0.48 ± 0.0624 | 296 |
P3HT/PS (without annealing) | 0.631 | 423 |
P3HT/PS (with annealing) | 0.3731 | 423 |
P3HT/PS | −0.0433 + (49.3 − 11.8 ϕP3HT)/T31 | — |
P3OT/PS (without annealing) | 0.7431 | 423 |
P3OT/PS (with annealing) | 0.3531 | 423 |
P3OT/PS | −0.0464 + (44.0 − 2.86 ϕP3OT)/T31 | — |
P3DDT/PS (without annealing) | 1.4931 | 423 |
P3DDT/PS (with annealing) | 0.231 | 423 |
P3DDT/PS | −0.0276 + (31.7 − 0.489 ϕP3DDT)/T31 | — |
We put a particular emphasis on how chain conformations are governed by the persistence length and described how this parameter can be both experimentally determined and computationally predicted. Besides the stiffness of the polymer backbone, other factors such as side chains, non-bonding interactions, and defects have also been shown to influence chain conformations. In discussing the phase behavior of conjugated polymers, we reviewed Flory–Huggins theory and described the phase behavior of homopolymers, polymer–polymer blends, and polymer–fullerene blends. In both cases, the field has seen significant progress in the last few years, although experiments, simulations and theory are warranted to resolve discrepancies between reported values for fundamental parameters such as the persistence length and the Flory–Huggins interaction parameter χ.
While conjugated polymers have many attractive qualities, their underlying physics is complex, with some areas still under debate. For example, despite recent progress, a direct connection between the chain conformation, microstructure and charge transport efficacy remains elusive. The phase behavior of conjugated polymers is not fully mapped out, yet it is important as creating mixtures is a route for achieving unique functionalities. As methods for probing the fundamentals of these interesting materials continue to improve, their structure–property relationships can be better understood, and such knowledge will inevitably lead to novel materials.
This journal is © The Royal Society of Chemistry 2017 |