Nematic Order Imposes Molecular Weight Effect on Charge Transport in Conjugated Polymers.

Nematic order, in the bulk or at interfaces, is ubiquitous for semiflexible conjugated polymers. Nevertheless, the effect of liquid crystalline order on charge transport remains unclear. Using an analytical model, we demonstrate that nematic order leads to an enhancement in charge mobilities when compared to isotropic chains. Furthermore, we predict a quadratic dependence of the charge mobility on molecular weight of the chains. Analysis of the probability of forming hairpin defects also shows how the persistence length affects charge transport in conjugated polymers. We speculate that the prevalence of nematic order in conjugated polymers explains the reported increase in charge mobilities with molecular weight.

Introduction

Conjugated polymers are promising materials for many flexible electronic applications, including displays and stretchable transistors.[1−3] Despite recent progress in the development of these soft semiconductors, charge transfer in conjugated polymers is still not fully understood. Indeed, models have been proposed based on experimental trends to qualitatively capture the dependence of charge mobility on morphological parameters, such as crystallinity and intercrystallite connectivity for semicrystalline conjugated polymers.[4,5] Still, the lack of mechanistic models inhibits quantitative predictions of charge transport behaviors based on molecular properties for conjugated polymers, such as the drastic molecular weight dependence of the charge mobility.[6]

Because many high-performance conjugated polymers are semicrystalline, or even nearly amorphous,[7,8] charge transport on amorphous chains plays an important role in overall electronic properties of conjugated polymers. Amorphous chains may orient isotropically or exhibit some degree of liquid crystalline order, depending on processing conditions and their abilities to access a liquid crystalline phase. Even for strongly crystalline polymers, accessing a liquid crystalline phase prior to crystallization can be important to enhance both local order and charge transport.[9] For semicrystalline conjugated polymers, macroscopic transport is likely limited by intercrystallite transport through amorphous tie chains that connect adjacent crystals.

For amorphous conjugated polymers, chain conformations can significantly affect electronic properties of samples. Using density functional theory (DFT) calculations, previous work has shown that distortion of backbone dihedral angles can localize charge carriers, and, in turn, hinder intrachain charge transfer in amorphous chains.[10,11] Under the assumption that chains are isotropic, Pearson et al. have predicted that charge mobilities scale linearly with molecular weights before leveling out at a high molecular weight.[12] In the model based on Marcus theory[13] and kinetic Monte Carlo (KMC) simulations developed by Spakowitz and co-workers, field-dependent charge mobilities increase with molecular weight for short isotropic chains as a result of reduction of chain ends that bound the length scale of fast intrachain charge transport.[14−16] When the chain length is greater than the Kuhn length, polymer backbones can bend noticeably, creating conformational traps for charge carriers, and, in turn, leading to constant charge mobilities at high molecular weights.

The local alignment of polymer backbones, such as in nematic order, can also be crucial for charge transport in amorphous chains. Uniaxially aligned polymer backbones can act as straight pathways for charge transfer and lead to enhanced charge mobility in the aligned direction. The anisotropic shapes of the nematic chains may also affect statistics of charge trajectories, which then leads to a dependence of charge transport on material properties, such as molecular weight and persistence length.

In fact, conjugated polymers can spontaneously exhibit nematic order in the bulk and at interfaces. In the bulk, the semiflexible backbones of conjugated polymers may be regarded as sequences of rod-like segments with length of about a Kuhn length, lk, dispersed either in a real solvent or in “bound solvent” composed of alkyl side chains. When the polymer concentration is high enough, randomly placing rod-like segments with a random orientation in the system is no longer entropically viable.[17−19] The system then undergoes a first-order transition from an isotropic phase to a nematic phase, in which chain segments align uniaxially. Experimental evidence for liquid crystalline phases abounds for many polythiophenes[9,20−23] and donor–acceptor conjugated polymers.[24,25]

Near impenetrable surfaces, such as gate dielectric layers in field-effect transistors (FET), conjugated polymers can exhibit surface-induced nematic order. Instead of bending chain trajectories sharply to avoid collisions with the walls, semiflexible polymers tend to align parallel to impenetrable interfaces, creating alignment layers of thickness of about a persistence length, lp.[26−30] Nematic interactions between backbone segments can further enhance surface-induced alignment. For chains that are near the isotropic-to-nematic phase boundary, surface-induced alignment is uniaxial, with the nematic director parallel to the interface.[29,30]

Investigating the effect of uniaxial nematic alignment on charge mobility for amorphous conjugated polymers is, therefore, important for understanding overall charge transport in electronic devices. For example, regardless of whether it is nematic in the bulk or not, conjugated polymers near interfaces may be uniaxially ordered. In FETs, charge carriers are localized within 5 nm thick layers at the dielectric surface.[31] Because the value of lp of high-performance conjugated polymers usually exceeds 3 nm,[32] we expect that amorphous chains align uniaxially in the conductive channel. Thus, modeling of charge mobilities of conjugated polymers requires a description of the role of nematic alignment on charge transport.

In this work, we use an analytical model to investigate charge transport for uniaxially aligned conjugated polymers. We show that nematic order may lead to the observed molecular weight dependence of charge mobility in conjugated polymers. We predict that charge mobilities are enhanced by nematic order and increase quadratically with molecular weight. The bending stiffness of chain backbones, characterized by the persistence length, governs formations of hairpin defects, which in turn alter the shapes of intrachain trajectories and the macroscopic charge mobility of nematic chains.

Results

Charge Transfer in Isotropic Chains

In our model for charge transport, intrachain charge motion is regarded as a one-dimensional diffusion of charge carriers with reflecting boundaries at the chain ends. In the absence of applied electric fields, charge carries diffuse randomly along polymer backbones. The probability density for a charge carrier that is injected at time 0 into monomer s to be found at monomer s′ at time t iswhere N is the chain length and τi is the characteristic time for a charge to diffusively spread along a chain of length N, defined as τi = N2π–1ki–1. The parameter ki is the mean local intrachain charge transfer rate. The effect of local distortions of backbone dihedral angles is not explicitly included in the model; rather, we assume that disorder in dihedral angles leads to a smaller ki.

On average, the relative distance (measured in number of monomers) between the location of a charge carrier at t and its initial position is

After diffusing on the same chain for a while, the charge carrier can hop to another chain. We assume that the time required for interchain hopping follows an exponential distribution:so that the mean lifetime for intrachain charge transfer is τh = ∫0∞dtth(t). The value of τh is inversely proportional to the interchain charge transfer rate, kh.

Between two interchain hopping events, the trajectory of charge motion is governed by the conformation of a semiflexible chain segment with length ⟨ΔN(τh)⟩:where γ = π(τi/τh)1/2 = N(kh/ki)1/2. The value of γ2 represent the ratio of the characteristic time for intra- and interchain charge transport. If chains are long enough and hopping between chains is frequent enough, then γ is large, and ⟨ΔN⟩ is much shorter than (and independent of) the total chain length N.

Using eq , we can obtain the mean square displacement per hopping time τh for charge carriers in the conducting direction as the chain dimension ⟨R∥2(τh)⟩. By regarding the long-time scale charge transfer as a one-dimensional random walk of step length ⟨R∥2(τh)⟩1/2, we obtain the charge mobility from the diffusion coefficient of charge carriers using the Einstein relation:

For chains in the isotropic phase, Pearson et al. obtained the charge mobility by treating conjugated polymers as random walks of Kuhn segments:[12]in which a is the length of a monomer and lp is the persistence length. The charge mobility scales linearly with persistence length, which governs the straightness of the intrachain charge trajectories for isotropic chains.

When interchain charge hopping is slow or molecular weight is low (τi ≪ τh), the charge carrier can explore the whole chain between two hopping events so that ⟨ΔN(τh)⟩ ∼ N. In this case, the charge mobility scales linearly with the molecular weight. When interchain charge transfer is fast or molecular weight is high (τi ≫ τh), a charge carrier will hop to another chain before it realizes how long the polymer is. For this case, the overall charge mobility is independent of molecular weight.

As a consequence, the charge mobility for isotropic chains increases linearly with molecular weight and reaches a plateau at high molecular weight. By normalizing the charge mobility with respect to the plateau value, one can obtain the molecular weight dependence of charge mobility:Depending on the ratio of intra and interchain charge transfer rates ki/kh, the molecular weight dependence may be significant for isotropic polymers.

Effect of Nematic Order on Charge Mobility

We now consider charge transport in semiflexible chains with uniaxial alignment. Assuming the intrachain charge transfer rate, ki, is not affected by uniaxial alignment, we can combine the chain dimension in the alignment direction with eq to predict the anisotropic charge mobility for conjugated polymers.

For conjugated polymers with nematic order, chain backbones align along the nematic director . The degree of alignment can be described using a quadrupolar order parameter q:where θ is the deflection angle between a backbone tangent and the nematic director.

Thermal fluctuations can lead to deviations of chain paths from the aligning direction . The nematic field, however, can constrain transverse fluctuations of chain paths Δr⊥2. We can define a deflection length λ below which the transverse fluctuations of nematic chains are governed by the bending stiffness of semiflexible chains and are unaffected by the nematic field. The average transverse fluctuation of the ends of the deflection segments (of length λ) isEffectively, the uniaxial alignment field acts as a cylindrical confinement of radius Δr⊥ and axis in the direction.[33−36]

When lp is greater than λ, the deflection segments are rather straight, and the nematic order in the system can be approximated aswhere θλ is the angle between the deflection segment and the nematic director .

Combining eqs and 10, one can also write the deflection length, λ, as a function of the nematic order, q, and the chain persistence length, lp:

In the uniaxial alignment direction, charge carriers travel along some rod-like paths with sequences of length, λ⟨cos θλ⟩. As a consequence, the intrachain charge trajectory between two interchain hopping events follows the statistics of rod-like chains:

By treating long-time charge transfer as a one-dimensional random walk in the alignment direction, we predict the charge mobility in the aligning direction for uniaxially ordered chains:The charge mobility of uniaxially ordered (q > 0) chains can be about 3 orders of magnitude higher than the charge mobility of isotropic chains (Figure a). The rod-like intrachain charge trajectories lead to the drastic enhancement of charge transport. Indeed, the nematic order, q, may also depend on chain length, N, because longer chains have fewer chain ends and may interact more effectively with the aligning field. Once chains transition into the nematic phase and become uniaxially aligned, however, further increases in the degree of alignment, q, result only in a slight enhancement of charge mobility, about a factor of 2. As such, we expect that our model, in which q is independent of N, can capture the molecular weight dependence of charge mobility for nematic conjugated polymers.

Figure 1

(a) Charge mobility for isotropic chains and nematic chains with different order q. (b) Charge mobility for nematic chains and nematic chains with hairpin defects. Reference mobility is and a is the monomer size. Relative charge transfer rate is ki/kh = 105 and persistence length is Np = 5. (c) Normalized charge mobility for isotropic chains, nematic chains, and nematic chains with hairpin defects. Relative charge transfer rate is ki/kh = 105. Np = 5 (green) and 10 (purple).

We also obtained the molecular weight dependence of charge mobility for uniaxially aligned semiflexible chains. When τi ≪ τh, the charge mobility in the alignment direction scales with N2. This molecular weight dependence is the same as the prediction by Pearson and co-workers for charge transport on nematic rod-like chains.[12] When interchain charge transfer is very fast or the molecular weight is high (τi ≫ τh), the charge mobility of aligned polymers is independent of molecular weight. We write the normalized charge mobility in the alignment direction asThe molecular weight dependence is more significant for the uniaxially aligned polymers than for the isotropic chains (Figure c).

In the direction perpendicular to the aligning field, the intrachain charge trajectory is governed by the random transverse fluctuations of the deflection segments:The resulting charge mobility isThe molecular weight dependence of the charge mobility of nematic chains in the perpendicular direction is the same as that of isotropic chains. The value of μ⊥ also decreases as the nematic order, q, increases. We expect that the charge mobility in the aligning direction is much higher than the mobility in the perpendicular direction.

Role of Hairpin Defects on Charge Transport

In the above analysis, we assume that every segment of the aligned chains uniaxially points in the alignment direction. The assumption is true when the chain contour length, l, is comparable to the persistence length, lp. For short semiflexible chains, forming a hairpin defect involves sharply bending the stiff polymer backbone against the nematic field. As a consequence, the energy penalty for forming short hairpins is large and the probability for finding those defects is negligible. When chains are long enough, however, semiflexible polymers can form hairpin defects by bending the chain trajectory more smoothly over a length greater than lp. In this way, the elastic energy penalty for bending a chain is smaller and the probability for forming hairpins is higher.

The existence of these “U-turns” can affect the statistics of charge trajectories. If intrachain charge transport occurs on a length scale shorter than the mean distance between two hairpin defects, Na (Figure a), the charge trajectory is rod-like (eq ). When the contour length of the intrachain charge trajectory is much longer than Na, intrachain charge motions are essentially one-dimensional random walks with step length Na.

Figure 2

(a) Hairpin defects for nematic chains. Na is the mean chain length between two hairpin defects. Nb is the length of the elbow for a hairpin defect. (b) The chain path for a semiflexible hairpin with length N = 20. Ratio of aligning field and stiffness ω = A/Np.

To reveal the effect of hairpin defects on charge transport for polymer nematics, we construct a simple model to estimate the bending energy cost of a hairpin, Eh, and the mean chain length, Na, between two defects. In the model, we align semiflexible chains using an external aligning field, A, and write the single chain Hamiltonian asin which Np is the persistence length measured in a number of monomers (Np = lp/a), and (s) is the normalized chain path, in which each segment is a unit vector. The interaction of a chain segment, s, with the external aligning field iswhere θ is the angle between the sth chain segment and the aligning axis, and P2 is the second order Legendre polynomial.

For nematic polymers, we expect that the energy cost for forming a hairpin defect is expensive, of at least several kT. As a consequence, the conformation of the hairpin must be dominated by a classical chain path that possesses the lowest free energy. Assuming the fluctuations around the most favorable chain path are negligible, we obtain the hairpin conformation by minimizing the free energy:Because the alignment is uniaxial, we are only interested in the deflection angle θ of the chain path from the alignment axis:in which . Using eq , we can solve the most favorable conformation of a hairpin by constraining the orientation of the chain ends so that θ(0) = 0 and θ(N) = π.

In absence of the aligning field, A (ω = 0), the solution of the above equation is evidently θ(s) = πs/N. The hairpin is a semicircle so that the deflection angle θ linearly varies over the chain length (Figure b). The stiff polymer backbone bends uniformly to reduce the elastic energy cost.

For systems with uniaxial alignment (so that A and ω are greater than zero), we analytically solve eq by realizing the differential equation is essentially the equation of motion for a particle subjected to a force given by −ω sin 2θ(s):When the strength of the uniaxial aligning field, A, increases (and so ω = A/Np also increases), the radius of the hairpin decreases (Figure b). Although the bending energy is higher for smaller hairpins, the cost of alignment against the field, A, is lower.

Using the most favorable chain path, θ(s) (eq ), we can obtain the characteristic size of the hairpin elbow, Nb (see Figure a):We can also compute the bending energy cost per hairpin defect for long nematic chains:The value of Eh arises from the optimal size of the hairpin elbow (Nb) and increases with both the aligning field, A, and the chain stiffness, Np.

Because the bending energy of a hairpin can be as large as several kT, the probability for finding these defects is low. We may consider the hairpin defects as a one-dimensional ideal gas, distributed randomly on the uniaxially aligned chains. The “volume fraction” of the hairpin defects can be approximated asand we can write the free energy for the “hairpin gas molecules” asBy minimizing the free energy with respect to Na, we obtain the mean chain length between two defects:

To demonstrate the typical value of Na for uniaxially aligned conjugated polymers, we approximate the aligning field strength, A, for long semiflexible chains in the nematic phase. We estimate the aligning field that a chain segment experiences as the product of the average order in the system, q, and the nematic coupling parameter, α. For nematic chains near the isotropic-to-nematic phase boundary, we expect that A = αcqc, where αc is the critical nematic coupling constant and qc is the critical order parameter for the system to undergo a first order transition from the isotropic phase to the nematic phase. In our pervious work, we estimated the values of αc and qc for chains with differing stiffness, Np.[23] For long chains (N ≫ Np), the critical coupling parameter, αc ∼ Np–1, and qc are rather constant for chains with differing stiffness. As a consequence, we estimate the aligning field strength, A ≈ Np–1 (a coefficient of order unity is neglected here). The mean chain length between adjacent hairpin defects now depends only on the chain stiffness:The linear relation between Na and Np agrees with the previous prediction by Vroege and Odijk, which is obtained by applying a Gaussian trial function to approximate the narrow distribution of backbone tangent orientation near θ = 0 and π for nematic chains.[37] For mildly stiff chains, we expect Na to be as large as 102 (Figure ).

Figure 3

Mean distance between two hairpin defects, Na vs persistence length Np, for chains in the nematic phase.

When the chain length is much longer than the typical distance between two defects (N ≫ Na), the intrachain charge trajectory is no longer rod-like. We expect charge motion along the polymer backbone is now a one-dimensional random walk with step length equal to Naa(2 + q)/3:The resulting charge mobility depends linearly on the persistence length, lp, comparable to the charge mobility for isotropic chains. For nematic chains without hairpin defects, however, the charge mobility is independent of the persistence length (eq ).

Our simple model captures the chain dimensions of uniaxially aligned semiflexible polymers in both the ballistic (R∥2 ∼ N2) and random walk (R∥2 ∼ N) regimes. Our results are consistent with previous calculations by Spakowitz and Wang.[38] Using self-consistent field theory and the ground state dominance approximation, these authors analytically compute a correlation length, ξ∥, that separates the ballistic and random walk regimes for chain conformations in strong aligning fields. The correlation length, ξ∥, represents the critical length for hairpin formation. Similar to our Na, the value of ξ∥ linearly scales with Np at the isotropic-to-nematic phase boundary (assuming the critical value of ANp is a constant for the phase transition).

We write the molecular weight dependence of charge mobility using a piecewise function:in which the crossover point, Nc, can be solved by equating the two portions of eq . The value of Nc depends on the arm length of hairpin defects, Na, and the relative charge transfer rate, ki/kh.

To observe the transition of the intrachain charge trajectory from the rod-like scaling regime (μ ∼ N–2) to the coil-like regime (μ ∼ N–1), the relative charge transfer rate, ki/kh, needs to be large enough thatIf ki/kh is small, a charge carrier can hop to another chain before it realizes the existence of hairpin defects.

Nonetheless, the existence of hairpins does not significantly alter the overall molecular weight dependence and the value of charge mobility (Figure b). For uniaxially aligned chains with hairpin defects, we expect that the charge mobility increases with N2 when the molecular weight is low, comparable to the mobility of defect-free chains. When the chain length exceeds Nc, which is much greater than the hairpin arm length, Na, the dependence on molecular weight becomes small, similarly to our predictions for isotropic coils. Because Nc is rather large, the charge mobility for N > Nc can be close to the plateau value and the coil-like behaviors of the charge trajectories may not be obvious. In fact, hairpin defects only reduce the charge mobility of a high molecular weight sample by a factor of about two.

Comparison with Experiments

Because the ratio of intra and interchain charge transfer rates, ki/kh, can be large, the enhancement in charge mobility can be many orders of magnitude as the molecular weight increases for nematic chains (Figure c). For example, by fitting the experimental mobility data, Spakowitz et al. suggest that the ratio, ki/kh, can be as large as 105 for amorphous poly(spirobifluorene) and 104 for poly(3-hexylthiophene) (P3HT) in amorphous regions.[14,16] Using microwave conductivity measurements, previous authors have found that the intrachain charge transfer can be 103–104 times faster than overall charge transport for fluorene-thiophene copolymers.[39] We expect the ratio, ki/kh, to be smaller for high-performance materials, but still much greater than unity.

Our analytical charge transfer model for nematic chains may explain the strong molecular weight dependence of charge mobility for conjugated polymers. Current measurements of charge mobilities rely on thin-film transistors, where the conductive channel is localized to about 5 nm from the dielectric interface. In the conductive region, the dielectric interface can induce uniaxial order to semiflexible conjugated polymers.[29,30] Indeed, the nematic order of real conjugated polymers may not be truly long-ranged. Any quenched disorder can disturb long-range uniaxial order, and lead to nematic domains in which the liquid crystalline directors point in different directions.[41] The sizes of nematic domains may affect long-range charge transfer in the sample. For example, using resonant scattering of polarized soft X-rays, Collins et al. showed that charge mobility can be enhanced by an order of magnitude as the nematic domain size, ξ, increases for liquid crystalline poly(2,5-bis(3-alkylthiophen-2-yl)thieno[3,2-b]thiophene (PBTTT).[40] When the sizes of ordered regions, ξ, are large or charge transfer across domain boundaries is sufficient, we expect that the molecular weight dependence of charge mobility arises from the “short-range” uniaxial order of chains in the nematic domains.

To show the effects of domain sizes and interdomain charge transport on net charge mobilities in nematic conjugated polymer samples, we make a simple argument as follows. In the absence of applied fields, the net charge mobility, μnet, is proportional to the long time scale diffusion constant of charge carriers, Dnet. Assuming the orientations of nematic domains are random, we can write the diffusion constant, Dnet ∼ ξD2/τd, in which ξD is the domain size and τd is the time for a charge carrier move from one domain to another. Inside the local nematic domains, the charge motion is a one-dimensional random walk, with diffusion constant Dlocal(N) ∼ ⟨ΔN(τh)⟩2. For intradomain charge transport, we consider diffusion in the direction parallel to the nematic field because charge transport in the perpendicular direction is much slower. In the direction perpendicular to the nematic field, the intrachain charge trajectory is a random walk with a step size decreasing with the nematic order, q. The resulting mobility, μ⊥, is much smaller than μ∥ for long chains. We expect charge carriers move more effectively to the domain boundaries in the aligning direction, followed by interdomain charge transport. In this way, we can write τd by summing the time for charges to move across a domain, ξD2/Dlocal(N), and an additional time, τb, for charges to hop across the domain boundaries. Thus, the net charge mobility scales asThe value of μnet increases as the nematic domain size increases and the interdomain transfer time, τb, decreases. When ξD is large or τb is small, the net charge mobility is proportional to ⟨ΔN(τh)⟩2, the same as our nematic model prediction.

Charge transport across misaligned domain boundaries may involve charge motion in disordered chains so that (D(N) ∼ ⟨ΔN(τh)⟩) and the required time τb ∼ ⟨ΔN(τh)⟩–1. When τb is much greater than the intradomain charge diffusion time, τd, we expect that the net charge mobility, μnet(N), is proportional to ⟨ΔN(τh)⟩, comparable to the charge mobility for isotropic chains. When τb and τd are comparable, we expect an intermediate scaling for the net charge mobility, μnet(N) ∼ ⟨ΔN(τh)⟩, where the exponent ν is between unity and two.

We compare the nematic model (eq ) with the experimentally measured charge mobilities for semicrystalline P3HT, obtained by different groups using various processing conditions (Figure a).[4] Because the nematic phase has been reported for P3HT,[20−23] we expect chains in amorphous regions to align uniaxially. The different degrees of crystallinity and the various sizes of nematic domains in the amorphous regions may lead to the rather scattered experimental charge mobilities. Nonetheless, our model can qualitatively capture the molecular weight dependence of charge mobility for P3HT by assuming a relative transfer rate, ki/kh, of about 103–104.

Figure 4

Charge mobility vs molecular weight for P3HT. Experimental data (dots). Analytical model for nematic chains (solid curves). (a) Mobilities obtained from various groups.[4] The average mobility of the four high molecular weight samples used as plateau mobility μ∥(∞). (b) Mobilities reported by Chang et al.[6] Curve obtained from nematic model for charge mobilities with μ∥(∞) = 0.095 cm2/V–1 s–1 and ki/kh = 1240.

We also compare our model with the experimental mobilities obtained by Chang et al. for P3HT samples with similar degrees of crystallinity.[6] Using eq and the relative transfer rate, ki/kh, and the plateau mobility, μ∥(∞), as fitting parameters, the nematic transfer model nicely captures the molecular weight dependence of charge mobility for P3HT. Nematic-like order may be more important than semicrystallinity in P3HT for governing macroscopic charge transport. We propose that the straight backbones of nematic conjugated polymers may give rise to the significant molecular weight dependence of charge mobility.

Discussion

We propose an analytical model to describe the effect of uniaxial alignment on charge transport for conjugated polymers. We show that charge mobilities of uniaxially aligned chains depend more significantly on molecular weight than charge mobilities of isotropic chains. When the molecular weight is low, we predict that the charge mobility in the alignment direction, μ∥, for nematic chains quadratically scales with the chain length (∼ N2). Indeed, nematic conjugated polymer samples are composed of domains with different orientations. Still, when interdomain charge transport is sufficient or the domain size, ξ, is large, we expect that the local rod-like charge trajectories lead to a quadratic molecular weight dependence of the overall charge mobility. For isotropic chains, however, μ∥ in the conducting direction only scales linearly with N. This is because the intrachain charge trajectories are rod-like for nematic chains in the alignment direction and coil-like for isotropic chains. For both nematic and isotropic chains, their charge mobilities reach plateau values when the time for a charge carrier to explore the whole chain is comparable to the interchain hopping time.

The charge mobility, μ∥, for nematic chains without any defects is expected to be independent of the persistence length, lp, although μ∥ of isotropic chains depends linearly on lp. For isotropic chains, the straightness of local charge trajectories is governed by the bending stiffness of the polymer backbone so that the charge mobility depends on lp. For nematic chains, however, the rod-like intrachain charge trajectories are controlled by the nematic order, q, instead of lp. As a consequence, μ∥ scales with (2 + q)2 and is independent of lp for nematic chains.

Together with the alignment field, the persistence length, however, governs the probability for forming hairpin defects in polymer nematics, and, in turn, they affect the overall charge mobility. For chains with longer lp, the bending energy penalty for forming a hairpin is larger, and the concentration of these defects is lower. As a consequence, the distance between two intrachain hairpin defects, Na, increases with increasing persistence length, Np. When the chain length, N, is smaller than Na, the intrachain charge trajectory is not affected by the existence of hairpin defects and is still rod-like. When N is greater than Na, the nematic polymer backbones become one-dimensional random walks of rod-like pieces of length about Naa. With a large enough relative charge transfer rate, ki/kh > Na2, we expect intrachain charge trajectories to also behave like one-dimensional random walks. Because Na scales linearly with Np in the nematic phase, we expect the charge mobility of nematic polymers in the high molecular weight limit to also scale linearly with the persistence length.

For many conjugated polymers, however, hairpin defects may not affect the charge mobility, because their relative charge transfer rates, ki/kh, are not large enough, only about 103–104, smaller than the typical value of Na2. A charge carrier can hop to another chain before realizing the existence of hairpin “U-turns”.

Using a nematic charge transfer model, we demonstrate that nematic order and the rather large relative transfer rate, ki/kh, may give rise to the significant molecular weight dependence of the charge mobility for semicrystalline P3HT. Unfortunately, the molecular weight dependence of the mobility in other conjugated polymers has not been characterized in the same level of detail. Further efforts are warranted to compare our model with experimental results from both amorphous and stiffer conjugated polymers.