Charge carrier coherence and Hall effect in organic semiconductors

Hall effect measurements are important for elucidating the fundamental charge transport mechanisms and intrinsic mobility in organic semiconductors. However, Hall effect studies frequently reveal an unconventional behavior that cannot be readily explained with the simple band-semiconductor Hall effect model. Here, we develop an analytical model of Hall effect in organic field-effect transistors in a regime of coexisting band and hopping carriers. The model, which is supported by the experiments, is based on a partial Hall voltage compensation effect, occurring because hopping carriers respond to the transverse Hall electric field and drift in the direction opposite to the Lorentz force acting on band carriers. We show that this can lead in particular to an underdeveloped Hall effect observed in organic semiconductors with substantial off-diagonal thermal disorder. Our model captures the main features of Hall effect in a variety of organic semiconductors and provides an analytical description of Hall mobility, carrier density and carrier coherence factor.

Unique combination of processability, interesting optoelectronic properties and mechanical flexibility makes molecular semiconductors very promising for development of organic electronics. However many fundamental aspects of charge carrier transport in these materials are still not fully understood. Weak van der Waals intermolecular interactions in these materials may lead to formation of rather narrow (0.1-0.3 eV) electronic bands of extended states 1 , that can be relatively easily destroyed by thermal molecular fluctuations 2,3 . Such fluctuations in molecular crystals can create significant off-diagonal disorder, leading to localized states in the gap below the mobility edge (the tail states) [2][3][4] . Thus, the intrinsic factors that compete in determination of the dominant transport mechanism in organic semiconductors include intermolecular interactions, governed by the equilibrium molecular positions in the crystal (the transfer integrals) 1 , carrier self-localization due to polaron formation (local electron-phonon coupling) 5,6 , and off-diagonal thermal disorder (non-local electron-phonon coupling) 2,3 . This leads to a host of small-molecule organic semiconductors with the intrinsic (that is, not dominated by static disorder) charge transport mechanism varying from an incoherent hopping in localized states to a pure band transport in extended states (see, e.g. 7 and refs therein). For instance, several high-performance organic semiconductors, with a band transport and carrier mobilities in the range μ ~ 1-20 cm 2 V −1 s −1 , have been recently identified 7 . On the other hand, there are crystals that seem to be at the borderline between a coherent band transport and an incoherent hopping. Thermal disorder, undoubtedly detrimental for robust band transport, is affected by specific molecular structure and crystal packing, and it can be suppressed at low temperatures [2][3][4] . Besides these intrinsic factors, static disorder (chemical impurities and structural defects) also plays an important role in organic semiconductors by leading to in-gap trap states that immobilize charge carriers at various time scales 8,9 . It is important to emphasize that localized tail states in pristine crystalline organic semiconductors can originate just from off-diagonal thermal disorder (thermal fluctuation of molecules with respect to each other). Such tail states are therefore to be considered as an intrinsic phenomenon occurring even in defect-free highly ordered single crystals. The important distinction, however, between thermally-induced and defect-induced disorder is that the former has a dynamic and homogeneous character (that is, it occurs uniformly throughout the entire crystal between every pair of lattice sites), while the latter is "frozen" at specific diluted spatial locations. The concentration of physical defects in high-purity molecular crystals is usually much smaller than the density of molecules. For instance, molecular density of crystalline rubrene is 1.43 × 10 21 cm −3 , while the density of rubrene peroxide defects at the surface of severely photooxidized crystals is ~0.2 × 10 21 cm −3 and only ~0.7 × 10 19 cm −3 in pristine crystals (in the bulk, this concentration is even lower) 10 . Thus, the density of tail states due to thermal disorder in such crystals is expected to be much greater than that due to physical defects, and therefore thermal disorder should play a dominant role in determining the charge transport mechanism in these materials Scientific RepoRts | 6:23650 | DOI: 10.1038/srep23650 (provided, of course, that it wins over the π -π intermolecular interactions). Based on this consideration, one would anticipate that materials, whose charge transport is dominated by thermal disorder (or, conversely, by band forming π -π interactions), will likely continue to behave as hopping (or, band-like) systems, even if an additional static disorder is introduced. Such disorder can be intentionally added in the form of chemical impurities or structural defects, produced, for instance, by photooxidation 10,11 , absorption of atmospheric gasses 12,13 , addition of molecular impurities during the crystal growth 14 , or ionizing radiation 15 . On the contrary, variation of temperature should have a much stronger effect on the transport properties of materials dominated by off-diagonal thermal disorder.
A conventional band-semiconductor Hall effect is an important signature in the studies of intrinsic charge transport in organic semiconductors, because its observation signals the presence of delocalized band-like states. Indeed, the classic Lorenz force considerations in the Hall bar geometry are readily applicable to delocalized band carriers (that is, carriers that have a well-defined microscopic drift velocity), but not to substantially localized carriers 7,16 . In such a simplistic picture, Hall effect would mainly probe instantaneously mobile charge carriers in extended states, thus yielding their concentration, n Hall , and intrinsic trap-free mobility, μ Hall . It is then remarkable that Hall effect measurements in organic field-effect transistors (OFETs) frequently lead to a surprising observation of not well-understood improper Hall effect, in which Hall carrier density, n Hall ≡ (e·R H ) −1 , where e is the electron charge and R H is a Hall coefficient (see below), appears to be greater than the total carrier density electrostatically induced in OFETs (excluding those trapped in deep traps), n FET ≡ e −1 ·C i (V G − V th ), where C i is a gate-channel capacitance per unit area, V G is the gate voltage, and V th is a threshold voltage in the linear-regime measurements of OFETs 17 . Simultaneously, Hall mobility, μ Hall , appears to be smaller than the longitudinal drift mobility, μ FET (n Hall > n FET , μ Hall < μ FET ). This behavior, which is also sometimes called an underdeveloped Hall effect, contradicts the simple interpretation of Hall effect as probing only instantaneously mobile charge carriers and therefore giving the highest possible, trap-free mobility (μ Hall ≥ μ FET ) 7 . An improper Hall effect has recently been observed in several systems, including pentacene 17 , 6,13-bistriisopropyl-silylethynyl pentacene (TIPS-pentacene) 18 , and even conjugated polymers 19,20 . It was empirically assigned to a partial carrier coherence, caused in these systems by strong off-diagonal thermal disorder, and parameterized by a carrier coherence factor, α , defined as the ratio of the total carrier density to that determined from Hall effect measurements, α ≡ n FET /n Hall 17 .
The postulated empirical relation between the carrier coherence factor and improper Hall effect does not offer a mechanistic picture of the phenomenon. It should also be noted that recently reported high-resolution ac Hall measurements have led to an observation which cannot be easily rationalized within that empirical framework. The ac technique allows for reliable measurements of Hall mobility and carrier density in OFETs with very low drift mobilities μ < 1 cm 2 V −1 s −1 , in which conventional Hall effect measurements are hardly possible even in very high dc magnetic fields of up to 12 T 21 . The ac Hall effect studies show that certain crystals expected to have a strong thermal disorder at room temperature, such as for instance tetracene 22 , still exhibit a normal (fully developed) Hall effect, even though their carrier mobility is as low as μ ~ 0.3 cm 2 V −1 s −1 21 . This observation thus does not support the view that less coherent carriers always lead to an underdeveloped Hall effect with α < 1.
Here, we propose a simple picture and the corresponding model description that could readily accommodate many salient features observed in recent Hall effect studies in organic semiconductors. The picture is based on the coexistence of band-like (delocalized) and hopping (localized) charge carriers in the accumulation channel of OFETs that respond differently to the applied magnetic field. Specifically, while the delocalized carriers exhibit the conventional Lorentz-force effect, the localized carriers are assumed to have a negligible response. It should be mentioned that Hall effect in purely hopping transport regime, originating from quantum-mechanical interference effects on the closed-loop trajectories, was theoretically predicted 23 , but apparently thus far has eluded clear experimental observations 24 . One may speculate that local hopping trajectories in systems with strong off-diagonal disorder can have a more one-dimensional character, thereby substantially decreasing the loop interference effects in comparison with those predicted for hopping carriers on regular lattices. This assumption appears to be an interesting subject for a dedicated theoretical study, which is beyond the scope of this paper. Of course, the hopping carriers would still experience a transverse Hall electric field, produced by the band carriers, and drift in the transverse direction opposite to the Lorentz force acting on the band carriers ( Fig. 1). In equilibrium, the resulting Hall voltage is correspondingly modified by the hopping carriers. The resultant equations derived below, relating Hall effect and FET mobilities and carrier densities, appear to capture most of the recently observed experimental trends.
For clarity, let us start with a simple textbook discussion of an accumulation channel comprising only mobile band carriers with density n and mobility μ. In the geometry of Hall effect measurements in OFETs, these carriers experience a Lorentz force, |F L | = |e[v × B]| = evB, where v = μ·E is a longitudinal terminal drift velocity of the carrier motion, occurring with a drift mobility μ in the longitudinal source-drain electric field E = V/L, V is the source-drain voltage applied to the accumulation channel of length L, and B is the magnetic field perpendicular to the channel. This force causes the carrier trajectories to bend toward one side of the channel, leading to a Hall voltage, V Hall , developing across the channel. The transverse Hall electric field, V Hall /W (W is the channel width), balances the Lorentz force, thus leading to an equilibrium according to the equation: eμ(V/L)B = eV Hall /W. Here, V/L and V Hall /W are the two orthogonal electric fields corresponding to the source-drain and Hall voltages, respectively. Such a system is described by the following relationship between Hall voltage V Hall , source-drain (longitudinal) current I, magnetic field B, and mobile carrier density n:

Hall H
where σ ≡ σ x = enμ is the longitudinal channel conductivity per square. It can be rewritten as: Scientific RepoRts | 6:23650 | DOI: 10.1038/srep23650 Hall H These formulas are normally used to extract the Hall carrier mobility, μ Hall , and Hall carrier density, n Hall , from experiments as: Note that in these equations, the current I, conductivity σ (or σ x ) and carrier density n are meant to represent only the mobile band carriers, because only these carriers experience the Lorentz force. If the transport is purely band-like, μ Hall and n Hall determined from Hall measurements according to Eq. 3 would then represent the actual mobility and density of mobile carriers. In a mixed transport regime of our model, however, when both band and hopping carriers are present in the accumulation channel and contribute to the longitudinal conductivity, direct application of Eq. 3 might result in an error. Indeed, the longitudinal conductivity in this case is: where σ band and σ hopping are the portions of conductivity due to band and hopping carriers, with the corresponding densities and mobilities of these carriers, n band , n hopping and μ band , μ hopping . This shows that using the total conductivity σ or the total current I in Eqs 1-3 might lead to at least one problem: an overestimated mobile carrier density.
For illustration of such an error, consider applying Eq. 1 to a case of pure hopping conduction, when V Hall = 0 and I ≠ 0, which leads to an erroneous result n → ∞ . This would not occur, however, if one used the correct value for I, the current carried by mobile band charges, which is zero in this example. There are further important considerations to be taken into account as we describe below. Despite the absence of a transverse magnetic field effect on hopping carriers (within the model discussed in this work), these carriers should still affect the total Hall voltage presented across the channel, because they do "feel" the transverse electric field, E y ≡ V Hall /W, associated with the Hall voltage. This field pushes the hopping carriers in the direction opposite to the Lorentz force and results in the corresponding compensation of the electric charge at the opposite side of the channel (Fig. 1). This will cause Hall voltage to be smaller than it would be in the case of only band carriers present in the channel. Such a compensation effect can be significant, if sufficient number of hopping carriers coexist with band carriers (see below). It can be easily seen from Eq. 3 that this compensation of V Hall would lead to an underestimated mobility, that is, μ Hall will be smaller than the actual mobility of band carriers, μ Hall < μ band . It can also be seen from Eq. 3 that in such a system, the Hall carrier density will be overestimated and in principle can be greater than the total carrier density (see below). Indeed, in the case of a mixed transport regime, one cannot experimentally identify what fraction of longitudinal channel conductivity σ is due to band carriers and what is due to hopping carriers, and thus only the total longitudinal conductivity σ can be used for parameter extraction from Hall effect measurements. In addition, we have no choice but to use the experimental V Hall , which is already reduced due to the compensation effect described above. Thus, Eqs 1-3 yield: In Hall effect measurements, these carriers are skewed towards the opposite sides of the channel, because the band carriers experience both the Lorentz and electric forces, while the hopping carriers only respond to the electric forces. The motion of hopping carriers in the transverse direction leads to a (partial) compensation of the Hall voltage, which results in underestimated mobility and overestimated carrier density in Hall effect measurements, as compared to those obtained in longitudinal FET measurements (μ Hall < μ FET , n Hall > n FET ).
Scientific RepoRts | 6:23650 | DOI: 10.1038/srep23650 In Eq. 5, as we showed above, the first ratio is greater than unity, μ band /μ Hall > 1, while the second ratio is usually much smaller than unity, μ hopping /μ Hall ≪ 1. Thus, depending on the particular combination of the relative densities and mobilities of band and hopping carriers, n band /n hopping , and μ band /μ hopping , the resultant carrier density determined from Hall measurements within conventional model (Eqs 1-5) might be greater than the actual total carrier density. We remind the reader that in OFETs, the total gate-induced carrier density in the accumulation channel, excluding those trapped in deep trap states, is: 7 . This compensatory effect might lead to an improper (underdeveloped) Hall effect. For a self-consistent analysis of the problem, we will now suppose that accumulation channel contains two types of charge carriers that have different carrier densities and mobilities: n 1 and μ 1 for the 1 st type, and n 2 and μ 2 for the 2 nd type, respectively. Let's also assume for certainty that carriers of the 2 nd type have a smaller mobility, μ 2 < μ 1 . If both types of carriers experienced the conventional Lorentz and electric forces in the geometry of Hall effect measurements, the Lorentz force acting on the 2 nd type of carriers would be smaller, and these carriers would contribute to the transverse current in the direction opposite to the transverse current flow of the 1 st type of carriers. The transverse current densities corresponding to these two types of carriers can be written as: where E x = V/L and E y = V Hall /W are the longitudinal source-drain and the transverse Hall electric fields. Here we assumed that the carrier mobility is isotropic, which is not necessarily the case in organic single crystals. However, more accurate derivation that takes into account mobility anisotropy leads to similar final equations, especially when the anisotropy is not too strong, or to identical equations in the case when the mobility anisotropy (the ratios μ a /μ b ) for hopping and band-like subsystems of carriers are similar. The dynamic equilibrium in the channel is achieved, when the total transverse current density is zero: j 1y + j 2y = 0.
This would give a relationship between the Hall electric field E y and the longitudinal current density j x for the case of coexisting two types of carriers: In the limiting case of a single type of carriers with an equivalent total carrier density n ≡ n 1 + n 2 and band mobility μ 1 , this formula gives the familiar result for the Hall effect in a band-like conductor: In the picture we explore in this work, however, the 2 nd type of charge carriers does not respond to the magnetic field. In such a case, balancing the transverse currents j 1y and j 2y leads to a modified formula (6): which can be rewritten as: We note that the Hall electric field E y given by Eq. (8) is always smaller than µ = E BE y x 1 in a purely band transport regime with the same total carrier concentration n ≡ n 1 + n 2 , irrespectively of particular values of carrier densities and mobilities of the two components (as long as μ 2 < μ 1 ). The measured Hall mobility in a system of mixed band and hopping carriers will then become: which is evidently always smaller than the intrinsic mobility of the band-like carriers, μ 1 . Here we introduced two dimensionless parameters: γ ≡ n 1 /(n 1 + n 2 ) = n 1 /n (0 < γ < 1) that represents the fraction of band-like carriers, and β ≡ μ 2 /μ 1 (0 < β < 1) that gives the ratio of hopping and band mobilities. The corresponding carrier density extracted from Hall measurements is: To compare the results of Hall and FET measurements, consider now the longitudinal channel conductivity, σ = e n 1 ·μ 1 + en 2 ·μ 2 = en·μ 1 ·(γ − γ β + β ), where the total carrier density n induced by application of a gate voltage, V G , is given by en ≡ en FET = C i (V G − V th ). The fraction of band carriers γ ≡ n 1 /n can be considered roughly independent of (or weakly dependent on) V G . Of course, in real systems with high density of localized tail states, this may not be fulfilled precisely, as the relative fraction of mobile carriers may increase with gate voltage. However, in the cases when super-linearities in FET's transfer characteristics are not too strong, we can assume that γ is almost V G independent. In such a case, longitudinal FET mobility μ FET extracted from transistor measurements can be expressed through μ 1 , γ and β as: FET 1 It is easy to see from Eqs (9-11) that: Hall FET FET Hall 2 which shows that there is an inverse relationship between the ratio of Hall and FET carrier densities and that of the corresponding mobilities: an overestimation of the carrier density is linked to an underestimation of mobility, or vice versa. Next, we show that depending on the combination of parameters, μ 1 , μ 2 , n 1 and n 2 (or, γ and β ) in Eqs 9-11, the Hall carrier density, n Hall , can be greater, comparable to, or smaller than the total carrier density, n, while the corresponding Hall mobility, μ Hall , can be smaller, comparable to, or greater than the FET mobility, μ FET , respectively. Indeed, these three cases can be derived directly from the above equations: (a) "Strong hopping contribution", with a high concentration of hopping carriers n 2 and relatively high hopping mobility μ 2 (that is, small γ , and/or not too small β ): These conditions directly follow from Eqs 9-11. For an example of the case (a) (Eq. 13), if the band-like carriers constitute only 1% of the total carrier population (γ = 0.01), and the mobility of the hopping carriers is 0.13 of the band mobility (β = 0.13), the measured Hall carrier density according to Eq. 10 will be about twice the total carrier density, n Hall ≈ 2n, the longitudinal FET mobility according to Eq. 11 will be μ FET ≈ 0.14μ 1 , while the Hall mobility according to Eq. 9 will be μ Hall ≈ 0.071 μ 1 , that is, μ Hall will be almost twice smaller than the longitudinal FET mobility, μ Hall ≈ ½·μ FET .
For an example of the case (b) (Eq. 14), if the band carriers are 50 times more mobile than the hopping ones (β = 0.02, which can occur for instance when μ 1 = 5 and μ 2 = 0.1 cm 2 V −1 s −1 ), but their number is small, say γ = n 1 /n = 4.16 × 10 −4 , then the above equations yield μ FET ≈ μ Hall ≈ 0.2 cm 2 V −1 s −1 < μ 1 , and n Hall ≈ n FET . That is, from the Hall measurements it will look like there is a good match between FET and Hall mobilities and carrier densities, but in fact the measured mobilities are far lower than the intrinsic band mobility.
Finally, for an example of the case (c) (Eq. 15), if γ = 0.25 (a quarter of carriers are band-like) and β = 0.1 (the band-like carriers are 10 times more mobile than hopping ones), we would have μ Hall = 0.77μ 1 , μ FET = 0.325μ 1 , and n Hall = 0.42n FET .
Note that case (12b), n Hall ≈ n FET and μ Hall ≈ μ FET , can be easily confused with the truly intrinsic band transport, in which all the carriers are band-like, and there is a good correspondence between Hall effect and FET measurements (such as, for instance, in pristine rubrene OFETs 16,25 ). It is clear that case (12b) is far from the intrinsic fully band-like regime, and it is still considerably affected by hopping (μ Hall ≈ μ FET < μ 1 ).
The above consideration shows that Hall effect measurements in a mixed transport regime may result in an incorrect estimate of the charge carrier mobility and density. Hall mobility can be underestimated, compared to the true mobility of band-like carriers (as well as relatively to the longitudinal FET mobility), while the Hall carrier density can be overestimated. In particular, the Hall carrier density can appear to be even greater than the total carrier density. This situation, sometimes reported in the literature, corresponds to the case of a "strong hopping contribution" described by Eq. 13. It has to be noted that in samples with substantial intrinsic or extrinsic disorder, we are likely to have a continuous mobility distribution, instead of just two "discrete" types of carriers, fast and slow, as we considered above. Thus, the simplified formulas derived here may not necessarily describe the actual disordered systems precisely. However, this model correctly captures the tendencies observed in Hall effect measurements in such systems. This simple result emphasizes the caution that should be exercised while interpreting Hall effect data in disordered systems with coexisting band and hopping carriers.
It would be instructive to plot the function f(γ , β ) = n Hall /n FET = μ FET /μ Hall given by Eq. 12 as a 3D plot of two parameters, γ and β (the Origin file with this function can be downloaded from Supplementary Materials). The function can be rotated around three axes by clicking on the plot and dragging it with a mouse, which gives a good idea about the 3D shape of the function. Figure 2a shows the function (the blue surface) and the constant n Hall /n FET = μ FET /μ Hall = 1 (the pink horizontal plane) given for comparison. The curve, along which the blue f(γ , β ) surface and the pink f = 1 plane intercept, describes the conditions in Eq. 14, for which Hall effect gives exactly the same result as FET measurements (n Hall /n FET = μ FET /μ Hall = 1).
We note that due to the specific shape of this function, a rather good correspondence between Hall effect and FET measurements is observed not only for the parameters along this intercept, but in a much broader parameter space. Indeed, the blue surface grazes very closely along the horizontal pink plane at large γ = n 1 /n, while going almost perpendicular to it for small γ and β near the origin. This shows that there is a large parameter space, where the function f(γ , β ) deviates little from unity. To visualize this, we have added two more horizontal planes, f = 1 ± 0.2, above and below the pink plane f = 1, to represent the range where discrepancy between Hall effect and FET measurements is within ± 20%. The result is shown as a top view of the function in Fig. 2b. The yellow and orange areas correspond to the parameter space, where Hall effect and FET measurements approximately match each other (within ± 20%). This suggests that once a system is in the so-called "ideal" or "fully developed" Hall effect regime (that is, when n Hall ≈ n FET and μ FET ≈ μ Hall , as in the case of Eq. 14, or in the case of a pure band transport γ = 1), additional variations in the relative concentrations of hopping and band carriers or their mobilities will not easily lead to a significant change of the type of Hall effect (or a sizable discrepancy between Hall and FET measurements). For instance, when the fraction of band carriers in Fig. 2d is γ > 0.8, Hall effect will approximately match FET measurements irrespectively of particular mobility ratio β . This explains the empirical observation that if there is a good correspondence between Hall and FET results, it is rather difficult to disturb it by simply adding extrinsic disorder (physical defects).
To illustrate this experimentally, we have prepared a series of rubrene OFETs, based on single crystals subjected to various degrees of photooxidation (see Methods for details). It is well established that photooxidation creates charge traps in organic semiconductors, which is detrimental for carrier mobility in these materials 10,11 . We have performed ac Hall effect measurements in this series of devices, including a pristine (nearly trap-free) rubrene OFET, an OFET based on a lightly photooxidized rubrene crystal, and finally an OFET based on a severely photooxidized crystal (Fig. 3). All the devices had the same type of FET structure prepared using parylene-N as a gate dielectric (for details on this type of OFETs, see e.g., ref. 26). The Hall measurements reveal a seemingly surprising behavior: while, as expected, the drift carrier mobility in this series of OFETs decreases, following a sequence μ = 4, 2.7 and ~0.4 cm 2 V −1 s −1 , not only can we measure Hall effect, but we do observe a nearly proper (fully developed) Hall effect in all of these devices, including the one with the greatest amount of added disorder and the smallest mobility of 0.4 cm 2 V −1 s −1 . Indeed, a good match between μ Hall and μ FET , as well as n Hall and n FET , is roughly maintained from pristine to severely photooxidized devices, even though μ is intentionally reduced well below 1 cm 2 V −1 s −1 (Fig. 3). Such "resilience" of Hall effect to the added disorder is understandable, given the shape of the function f(γ , β ) depicted in Fig. 2.
It can be seen from the above equations that the Hall voltage compensation has some interesting consequences for the temperature dependence of carrier mobility, μ Hall (T). In the case of a "strong hopping contribution" (Eq. 13), when n Hall > n FET and μ Hall < μ FET < μ 1 (that is, when γ ≡ n 1 /n << 1 is small, and β < 1 but not necessarily small), we obtain the following approximate expression for the Hall mobility from Eq. 9: Here, the band mobility μ 1 is expected to have a typical for band-like transport inverse temperature dependence μ 1 ∝ T −η with a power exponent η ≈ 1.5-2, which has been experimentally observed and theoretically modeled in band-like organic semiconductors 27,28 . The hopping mobility μ 2 is thermally activated, μ 2 ∝ exp(− Δ h /k B T), where Δ h > 0 is a characteristic energy barrier for hopping in localized tail states, and k B is Boltzmann's constant. In the "intermediate hopping contribution" case (Eq. 14), when n Hall ~ n FET , μ Hall ~ μ FET < μ 1 , and γ = (γ − γ β + β ) 2 from Eq. 12, Eqs 9 and 11 lead to the following approximate expression for the mobility: First of all, these equations show that, contrary to the common expectation that temperature dependence of Hall mobility should be representative of the intrinsic band mobility, μ 1 (T), this is not necessarily the case in a mixed transport regime: μ Hall (T) above is the band mobility μ 1 modified by the factors γ /β or γ 1/2 in the cases described by Eqs 16 and 17, respectively. Correspondingly, the measured temperature dependence of Hall mobility will differ from μ 1 (T), because the factors are temperature dependent themselves. Second, the temperature dependences of the prefactors in front of n 1 /n and (n 1 /n) 1/2 in Eqs 16 and 17 are rather different. The prefactor in Eq. 16 has a strong "band-like" character that tends to make the temperature dependence of mobility band-like (dμ/dT < 0). In Eq. 17, the prefactor also has a "band-like character", but weaker.
The fraction of band carriers γ = n 1 /n on the other hand must have a thermally activated temperature dependence. Indeed, in the model of carrier localization by dynamic disorder, we can think about individual carriers as intermittently existing in the two states: delocalized (band) or localized (hopping) states. The fraction of time an average carrier spends in one state or the other (τ band and τ tail , respectively) depends on the strength of thermal fluctuations, thermal energy, position of the Fermi level with respect to the mobility edge, etc. Thus, a crystal should always have populations of band and hopping carriers, with their relative densities related to these time scales as: n band /n total ≡ n 1 /n = τ band /(τ band + τ tail ). Given the typical total carrier densities in OFETs (n = 10 11 -10 13 cm −2 ) and the density of tail states in molecular crystals (10 14 -10 18 cm −3 eV −1 ) 4,8,9 , the Fermi level in the channel of an operating OFET typically occurs in the tail, separated from the mobility edge by an appreciable energy of a few hundred meV: E ME − E F ~ 0.1-0.3 eV, which is much greater than k B T at room temperature. Thus, generation of band-like carriers will mainly rely on thermal excitation from the tail states to (and above) the mobility edge. Therefore, temperature dependence of γ = n 1 /n can be described by an activation exponent: that competes with the "band-like" prefactors in Eqs 16 and 17. A counterintuitive prediction arising from Eqs 16-18 is that such a competition shifts the character of μ Hall (T) dependence toward a band-like type more significantly in the case of a strong hopping contribution (Eq. 16), compared to the intermediate hopping case (Eq. 17). Indeed, the extrema of functions in Eqs 16 and 17 (with the parameters taken to be η = 2, E ME − E F = 0.1 eV and Δ h = 0.026 eV for certainty) occur at T max (16) = (E ME − E F − Δ h )/ 4k B ~ 210 K in the case of a strong hopping contribution, and T max (17) = (E ME − E F )/4k B ~ 290 K in the intermediate case. That is, T max (16) < T max (17) , suggesting that when hopping contribution is strong, a "band-like" μ(T) dependence persists to lower temperatures than that at moderate hopping. Given the Hall voltage compensation effect described above, such behavior makes sense, because the motion of hopping carriers quickly "freezes out" with cooling (and so does the effect of V Hall compensation), leading to a stronger apparent "band-like" tendency in μ Hall (T).
Such a behavior has indeed been observed by J.-F. Chang et al. (see Fig. 1c,d in 18 ), where an apparently fully developed Hall effect, with n Hall = n FET and μ Hall = μ FET , was found for 1,4,8,11-tetramethyl-6,13-trie thylsilylethynyl pentacene (TMTES-pentacene), suggesting that all the charge carriers in that system are band-like (and experience a classic Lorentz force), thus leading to a fully developed Hall effect. On the contrary, TIPS-pentacene exhibited an underdeveloped Hall effect, with n Hall > n FET and μ Hall < μ FET , suggesting that the charge carriers are only partially coherent. Yet, surprisingly, the mobility μ Hall (T) in TMTES-pentacene OFETs was purely thermally activated, which is inconsistent with the band-like character suggested by the ideal Hall effect, while TIPS-pentacene exhibited μ Hall (T) closer to a band-like at high temperatures, with a crossover to a thermally activated regime at about T max ~ 220 K. If μ Hall (T) of TMTES-pentacene exhibited a maximum, it would likely appear at a higher temperature (outside of the measured range). Our model shows that this apparent inconsistency is superficial, as it is resolved by considering the Hall voltage compensation effect described above.
We next study in more detail the effect of mixed transport regime on the temperature dependence of mobility μ(T) in the case of moderate hopping contribution (Eq. 14). In this important case, an approximate match between Hall and FET carrier densities and mobilities is observed. We have carried out μ(T) measurements in OFETs based on intentionally photooxidized rubrene and pristine tetracene single crystals (Fig. 4). Since pristine rubrene OFETs have already been extensively studied, and a band-like (non-activated) μ(T) has been confirmed in these devices (see e.g. 7 and refs therein, and 16,27,28 ), we will not be reproducing this case here. It's worth noting that recently developed high-resolution ac Hall measurements allowed resolving Hall effect in single-crystal tetracene OFETs 21 . These measurements demonstrated that Hall effect in this system is proper (fully developed), even though the carrier mobility is rather low, μ ~ 0.3-0.7 cm 2 V −1 s −1 , similarly to the case of the photooxidized rubrene described above (Fig. 3b,c). Thus, both of these systems (photooxidized rubrene and pristine tetracene) seem to belong to the same category of systems with an intermediate hopping contribution (Eq. 14), with n Hall ≈ n FET and μ Hall ≈ μ FET < μ 1 . In this case, while the simplistic interpretation of Hall effect suggests that all the carriers are band-like, the experiment shows that the temperature dependence of mobility is, nevertheless, thermally activated (Fig. 4). This apparent inconsistency can be explained by the model considered above and specifically by Eqs 17-18, according to which the mobility in this regime is a product of the two competing factors, a OFETs based on pristine, as-grown tetracene crystals. In these two cases, Hall effect is proper (n Hall ≈ n FET and μ Hall ≈ μ FET ), as illustrated in the insets by matching Hall (triangles) and FET (solid line) carrier densities (see also ref. 21), yet μ(T) is thermally activated (solid black symbols in the main panels). Such behavior can be rationalized within the model developed in this work: μ(T) is fitted with Eqs 17-18 (solid red lines in the main panels), yielding reasonable activation energy, E ME − E F ~ 0.3 eV, and inverse-T dependence power exponent, η ~ 2. All T-variable measurements have been carried our using a gated 4-probe technique 26 .
It's worth noting that the carrier coherence factor, α ≡ n FET /n Hall , introduced by T. Uemura et al. 17 can be expressed in our model in terms of parameters γ and β as: It's temperature dependence can be analytically or graphically investigated (Fig. 5). Plotting α as a 3D function of two variables, γ and β , shows that when the initial value of α < 1 (corresponding to a point on the turquoise surface below the pink plane α = 1 in Fig. 5a), and temperature is reduced, thus causing the parameters β = μ 2 /μ 1 ∝ T η · −∆ e k T / h B and γ = n 1 /n ∝ − − e E E k T ( ) / ME F B to decrease (that is, move toward the origin), the point will typically move toward the spike (occurring near the origin), which corresponds to α increasing with cooling. A carrier coherence factor growing from 0.5 to ~0.65, when temperature is decreased from 300 to 155 K, has indeed been recently reported by T. Fukami et al. in Hall effect measurements of pentacene transistors 22 .
In conclusion, we show that Hall effect measurements in organic semiconductors in a mixed transport regime (with coexisting band and hopping carriers) may result in an incorrect estimate of the charge carrier mobility and density. The Hall mobility can be underestimated, compared to the true mobility of band-like carriers (as well as relative to the longitudinal FET mobility), while the Hall carrier density can be overestimated. In particular, the Hall carrier density can appear to be even greater than the total carrier density. This occurs because of a (partial) Hall voltage compensation due to the transverse drift of hopping carriers in the direction opposite to the Lorentz force. We have developed an analytical model of the Hall effect in organic semiconductors that expresses the Hall mobility μ Hall and carrier density n Hall , as well as the carrier coherence factor α in terms of the two microscopic parameters, the relative density of band carriers and the ratio of mobilities of hopping and band carriers. It has to be noted that in samples with substantial intrinsic or extrinsic disorder, we are likely to have a continuous mobility distribution, instead of just two "discrete" types of carriers, fast and slow, as we considered above. Thus, the simplified formulas derived here may not necessarily describe the actual disordered systems precisely. However, this model correctly captures the main tendencies observed in Hall effect measurements in such systems, including an underdeveloped Hall effect, peculiarities in the temperature dependence of carrier mobility and modifications of the carrier coherence factor α with temperature. This simple result emphasizes the caution that should be exercised while interpreting Hall effect data in disordered systems with coexisting band and hopping carriers.

Methods
Device fabrication. Rubrene and tetracene single crystals were grown using a physical vapor transport method in a flow of ultra-high purity He gas, at sublimation temperatures of 320 and 280 °C, respectively. Electrical contacts of OFETs were prepared on (a,b) facets of the crystals with an aqueous solution of colloidal graphite (Ted Pella). Parylene-N was used as a gate insulator. In tetracene FETs, a very thin under-layer of Cytop (< 30 nm) was spin-coated on the crystals before depositing a much thicker layer of parylene-N (C i = 1.7 nF•cm −2 ) (for details see ref. 21). For preparation of photooxidized rubrene OFETs, bare single crystals were photooxidized in a chamber filled with an atmosphere of ultra-high purity O 2 gas for 12 min (lightly photooxidized) and 2 hours (severely photooxidized) under a white-light illumination with a 85 mW·cm −2 Xenon lamp placed 10 inches away from the samples. After the photooxidation and contact deposition, the crystals were coated with parylene-N (C i = 0.8 and 1 nF•cm −2 for lightly and severely photooxidized crystals, respectively). 35 nm-thick Ag gate was thermally evaporated through a shadow mask on top of the parylene dielectric defining a channel with typical aspect ratio of L/W = 1-3, D/W = 0.3-0.7, where L is the channel length, W is the channel width, and D is the longitudinal distance between the 4-probe voltage probes.