A Large Anisotropic Enhancement of the Charge Carrier Mobility of Flexible Organic Transistors with Strain: A Hall Effect and Raman Study

Abstract Utilizing the intrinsic mobility–strain relationship in semiconductors is critical for enabling strain engineering applications in high‐performance flexible electronics. Here, measurements of Hall effect and Raman spectra of an organic semiconductor as a function of uniaxial mechanical strain are reported. This study reveals a very strong, anisotropic, and reversible modulation of the intrinsic (trap‐free) charge carrier mobility of single‐crystal rubrene transistors with strain, showing that the effective mobility of organic circuits can be enhanced by up to 100% with only 1% of compressive strain. Consistently, Raman spectroscopy reveals a systematic shift of the low‐frequency Raman modes of rubrene to higher (lower) frequencies with compressive (tensile) strain, which is indicative of a reduction (enhancement) of thermal molecular disorder in the crystal with strain. This study lays the foundation of the strain engineering in organic electronics and advances the knowledge of the relationship between the carrier mobility, low‐frequency vibrational modes, strain, and molecular disorder in organic semiconductors.


Introduction
Organic semiconductors are considered for numerous applications in flexible, wearable and disposable electronics due to their www.advancedscience.com Charge transport in "soft" van der Waals molecular solids is governed by the three key factors. The first factor is the equilibrium position of molecules in a crystal lattice that defines the average transfer integrals between the adjacent molecules. Large transfer integrals promote charge carrier delocalization and formation of electronic bands, thus leading to higher carrier mobilities. [7] The second factor is the thermal molecular motion responsible for the off-diagonal thermal disorder in transfer integrals. [8,9] In molecular crystals, this intrinsic disorder may lead to a breakdown or weakening of the extended band and occurrence of localized states near the band edges (the tail states), thus causing transient carrier localization. [10] It has been recently suggested that molecular vibrations most significantly limiting the charge carrier transport and reducing μ are among the low-frequency (1 -300 cm −1 ) modes in Raman spectra of organic semiconductors (see, e.g., ref. [11] and references therein). However, a direct link between the charge carrier mobility and these Raman modes has not been experimentally established. In rubrene, these modes were theoretically shown to be the result of a combination of inter-and intramolecular vibrations. [12] Finally, formation of small polarons (the carriers "dressed" in a polarization cloud) can also lead to a breakdown of the band and a reduced charge mobility. [13] In small-molecule organic semiconductors, these competing factors may depend on molecular structure, molecular packing, as well as external conditions (temperature, pressure, strain). As a consequence, we have systems with different transport properties: some materials exhibit a nearly perfect band-like transport with high mobility, such as, for example, rubrene (5,6,11, 12-tetraphenyltetracene, C 42 H 28 ), [14] as well as alkyl end-capped benzothieno-benzothiophene and dinaphto-thienothiophene (C 8 -BTBT and C 10 -DNTT), [15,16] some are characterized with a purely hopping conduction and low μ, and others show a "mixed" transport, in which band-like and hopping states coexist, as has been recently shown via high-resolution ac-Hall effect measurements of OFETs. [17] Both π-π interactions and thermal disorder can be delicately affected by intermolecular distances in crystals. For instance, in rubrene, the equilibrium transfer integral is a steep function of the relative displacement of the molecules, [7] and it was experimentally shown that application of a high hydrostatic pressure results in a significant increase of the carrier mobility of singlecrystal rubrene OFETs due to a reduction of the intermolecular distances. [18,19] It was also hypothesized that hydrostatic pressure or strain can reduce off-diagonal disorder by suppressing thermal fluctuations. [20] The most convenient experimental approach to study these intrinsic effects in organic semiconductors is to apply uniaxial mechanical strain to flexible single-crystal transistors. The advantages of strain-variable measurements over the experiments under a hydrostatic pressure are as follows: a) they are technically much easier to perform (a high-pressure anvil cell is not required), b) one can apply either compressive or tensile strains, c) strain can be applied along different crystallographic directions within the basal plane of the device, and d) there is practically no limitation on device size. Reyes-Martinez et al. recently demonstrated that wrinkling instability could be used as a tool to modulate the conventional two-probe mobility (the so-called longitudinal mobility), μ FET , in rubrene single-crystal OFETs with strain. [21,22] Strain induced modulations of μ FET were also recently demonstrated in solutionprocessed crystalline 3,11-didecyldinaphto[2,3-d:2′,3′-d′] benzo[1,2-b:4,5-b′]dithiophene (C 10 -DNBDT-NW) OFETs. [20] While such strain variable measurements in single-crystal OFETs have started to emerge, [20][21][22][23][24] to our knowledge there are no studies of the Hall effect versus strain in organic semiconductors. Such studies would reveal the response of the intrinsic (trap free) carrier mobility to strain, because Hall effect in highly ordered crystalline organic semiconductors tends to be sensitive primarily to the mobile carriers moving in extended states. [15,17,25,26] We emphasize that while applied studies show that polycrystalline or amorphous OFETs can withstand a considerable bending without mobility degradation, [3,6] there are only a handful of studies of the impact of strain on the charge transport in single-crystal organic semiconductors. [20][21][22][23][24] Given the relevance of high-mobility crystalline OFETs for organic electronics, [2,4,20,27,28] addressing the intrinsic effect of strain on μ due to the elastic lattice deformation is very important. Rubrene, in particular, stands out as a relevant test system because of its well-defined surface morphology, [29] crystal structure, [30] interesting optoelectronic properties, [4] ideal field-effect transistor behavior, [14,31] as well as a band-like charge transport in extended states with a high room-temperature mobility. [25,26] In addition, it has been recently shown that extremely thin (0.1-1 μm) and flexible rubrene crystals of macroscopic lateral dimensions (up to 1 cm) can be grown by physical vapor transport and integrated into OFETs. [21,22] Here, we quantify the strain effect on mobility by a strain factor, g, which is a relative change of the carrier mobility, Δμ/μ 0 , per unit of the applied strain ε, occurring in a FET on bending: g  |δ(Δμ/μ 0 )/δε|. Here, Δμ(ε) = μ(ε) − μ 0 is the difference between the mobility μ(ε) of a device under strain ε and the mobility μ 0 of the unstrained device. Both Δμ/μ 0 and ε in this definition are expressed either in percent or as decimal fractions. For instance, a strain factor of 100 would mean that the device under a compressive strain of 1% has a carrier mobility 100% greater than that of the unstrained device. Recent measurements of strain factors in crystalline OFETs have led to significant ambiguity in results. For instance, Wang et al. and Matta et al. observed a relatively weak modulation of μ with strain in twoprobe transistor measurements, yielding g = 20 and 6, respectively. [23,24] In contrast, Reyes-Martinez et al. reported a much stronger enhancement of μ with compressive strain in measurements of "wrinkled" two-probe FETs, with g up to 250. [22] Besides rubrene, the mobility-strain relationship was also evaluated in solution-processed single-crystalline C 8 -DNBDT-NW OFETs, also via longitudinal FET transport measurements, in which g ≈ 20 was reported. [20] These widely varying results on the apparent strain-induced mobility variations could be attributed to potential drawbacks of two-probe FET measurements, related for instance to the modulations of the metal contact work function with strain, recently observed by Wu et al. [32] In wrinkled devices, extraction of the net strain is not straightforward and has to rely on complex modeling. Finally, the conventional (longitudinal) FET transport measurements cannot easily distinguish between the intrinsic (not dominated by trapping) and extrinsic (dominated by trapping) conduction regimes. [25] To address these issues, here we employ the Hall effect studies coupled with four-probe FET measurements of ultrathin singlecrystal organic transistors subjected to a calibrated uniaxial strain. These flexible four-probe FETs allow us to avoid the artifacts caused by disorder and reveal the intrinsic mobility-strain relationship (for details on four-probe OFETs, see, e.g., ref. [14] and references therein). The Hall effect measurements were carried out via the high-resolution ac-Hall methodology capable of reliably detecting Hall signal and assessing the intrinsic (Hall) carrier mobilities in OFETs with μ well below 1 cm 2 V −1 s −1 . [33] Overall, we found that both the longitudinal FET mobility, μ FET , and the Hall mobility, μ Hall , of rubrene single-crystal OFETs exhibit a strong anisotropic dependence on strain, ε, with the anisotropy in μ(ε) response consistent with the molecular packing of the crystal. For charge transport along the high-mobility axis of rubrene (here, the b-axis), we measured the intrinsic strain factor of g = 70-110 for the strain applied along the same transport direction, and g ≈ 0 for the strain applied along the orthogonal low-mobility direction (here, the a-axis) of the crystal. These values are among the highest straininduced modulations of μ (and the greatest strain effect anisotropy) reported for organic semiconductors.

Flexible Single-Crystal OFETs and Measurement Approaches
The experimental arrangement is shown in Figure 1. The devices consist of ultrathin flexible rubrene single crystals laminated on semirigid 40 μm thick polyethylene terephthalate (PET) substrates precoated with 1 μm thick polydimethylsiloxane (PDMS) layer for better crystal adhesion. Before laminating a crystal, Ti/Au contacts in four-probe geometry are sputtered on PDMS surface via a contact shadow mask. A parylene-N gate insulator and Ag gate are deposited on top of the structures (Figure 1a,b). The (orthorhombic) rubrene crystals are aligned to the contacts, so that their high-mobility b-axis is the transport direction (Figure 1b,c). The crystal axes are identified via a polarized optical transmission microscopy (Section S1, Supporting Information). For Hall measurements, an ac magnetic field of an r.m.s. magnitude B = 0.23 T is applied perpendicular to the crystal's (a,b)-facet, when the FET is on, and ac-Hall voltage is measured by a lock-in amplifier ( Figure 1d). [33] Uniaxial strain parallel (or perpendicular) to the charge transport direction can be applied by bending the semirigid PET substrate supporting the device in a custom-designed miniature strain stage ( Figure 1e). The substrate can be bent inward (ε < 0, compressive strain) or outward (ε > 0, tensile strain) from its initial planar geometry (ε = 0, unstrained), as shown in Figure 1f. The strain magnitude is defined as where t is the net thickness of the substrate (t = 41 μm), and R is the bending radius (Section S2, Supporting Information). [34] The Hall mobility, μ Hall , is obtained from the Hall voltage, V Hall , and the longitudinal four-probe voltage, V 4p , simultaneously measured in a device carrying a source-drain current, I SD , in the accumulation channel, while the device is in the magnetic field B, according to where D is the distance between the voltage probes (in the longitudinal direction) in the four-probe contact geometry and W is the channel width (roughly the same as the distance Adv. Sci. 2020, 7, 1901824 between the Hall voltage probes). The strain factor g is determined from the experimental μ(ε) dependence as where subscript "0" refers to the value for the unstrained (ε = 0) device, and the strain ε is in %.

Measurements of FET and Hall Mobility under Strain
FET and Hall effect measurements in these flexible rubrene OFETs as a function of uniaxial strain are shown in It is worth noting that the device is functioning as a classic highperformance OFET with a linear transconductance and a nearly zero onset voltage. The mobilities, μ FET and μ Hall , are extracted from the extended linear portion of the curve at sufficiently high negative V G , where the mobility is gate voltage and carrier density independent. For more details on reliable extraction of carrier mobility in FETs, see ref. [14]. It can be clearly seen from this inset that the transconductance slope changes with strain. The inset in the lower panel shows examples of as-measured Hall voltage in this OFET for two levels of I SD excitation and the same three values of strain, also clearly showing that V Hall is sensitive to strain. In addition to the strain along the transport direction (along the molecular stacks), we have also evaluated the response of the FET and Hall mobilities to the strain applied perpendicular to the transport direction (perpendicular to the stacks). Several important observations can be made from these μ FET (ε) and μ Hall (ε) measurements. Both Hall and FET measurements yield consistent results, with similar absolute values of μ, which is increasing (decreasing) with compressive (tensile) strain applied along the transport direction. The mobilitystrain relationship μ(ε) shows a weak non-linearity around ε = 0 that may look like a threshold. Currently, the origin of this plateau-like feature at |ε| ≤ 0.15% is unclear, but it might be due to a partial stress relief in the PDMS elastomer used here for better crystal adhesion, or due to a very small amount of the initial strain being accommodated by residual structural defects of the crystal. Nevertheless, this should not affect our measurements of the strain factor that rely on the derivative, δμ/δε, of the quasi-linear portion of μ(ε) curves. The maximum strain in our measurements (±0.40%) is dictated by higher noise in V Hall occurring for ε outside of this range, which is perhaps due to the generation of additional defects under greater strains (for as-measured V Hall within this ε range see Section S3 in the Supporting Information). We have verified the elastic nature of the strain effect on mobility by first carrying out the measurements up to the maximum value of |ε| and then reducing |ε| step-by-step to several lower values, which yielded similar μ(ε) data (Figure 2), thus showing that within the errorbar of these measurements the data are reproducible on strain reversal. Anisotropy of the strain effect, consistent with the molecular packing in rubrene, is observed. The strain factors for FET and Hall mobilities are g FET ≈ 70 and g Hall ≈ 110 for the strain applied along the transport direction (along the molecular stacks), that is, μ increases by nearly a factor of two at a compressive strain of only 1% applied along the high-mobility b-axis. On the contrary, no changes in μ are observed for the strain applied perpendicular to the transport direction (perpendicular to the stacks, which is the low-mobility a-axis). Overall, μ Hall and μ FET are close to each other, which is indicative of a nearly trap-free (intrinsic) charge transport operation of these OFETs.

Eliminating Possible Artifacts
Spurious effects are rigorously minimized in our experiment. For example, inhomogeneous strain distribution due to strain accommodation by defects and/or metal contacts is unlikely, because a) very thin (100-300 nm thick) and uniform single crystals are used, b) strain is governed by a well-defined macroscopic bending radius of a much thicker (41 μm thick) substrate. The use of monolithic single crystals eliminates possible effects of grain boundaries. A good match between the fourprobe FET and Hall measurements, as well as the reversibility of the observed mobility-strain dependence, indicate that our devices are neither contact-limited nor trap-dominated, and we mostly probe the intrinsic (trap-free) charge transport properties of the crystals. In addition, within the strain range used here (up to ±0.4%), no observable delamination, wrinkling or cracking were detected in our devices. Since Hall effect directly measures the mobile carrier density, without relying on certain gate-channel capacitance, possible modulation of C i with strain does not affect our measurements. Flexo-or piezoelectric effects can be also ruled out, because we did not observe any measurable voltage generated between the Hall and four-probe contacts of bent devices, which is not surprising given the fact that these effects require the crystal to be non-centrosymmetric, or the applied strain to be inhomogeneous, neither of which is the case here.

Low-Energy Raman Spectra as a Function of Strain
To better understand the relationship between the intrinsic mobility and strain in crystalline rubrene, we have also performed measurements of low-frequency Raman spectra of thin flexible rubrene single crystals as a function of uniaxial strain (Figure 3). Although selection rules make Raman scattering "blind" to molecular vibrations important for the charge transport, such as acoustic modes far from the center of the Brillouin zone (Γ-point at q = 0), strain-dependent Raman measurements can still serve as a proxy, assuming the modes far from the Γ-point exhibit a similar behavior with strain. [35] Specifically, the low-energy region of Raman spectra (0-300 cm −1 ) is believed to be important for charge transport properties of organic semiconductors, [11] as it contains low-frequency intra-and intermolecular vibrations reflecting the off-diagonal thermal disorder. [12] It was recently proposed theoretically that the higher the frequency, f, of the lowest-energy Raman peaks of an organic semiconductor, the better its charge carrier mobility. [36] Alternatively, relatively low intensities (or a complete absence) of Raman peaks in this low-frequency region would suggest  In both panels, the charge transport direction is along the high-mobility b-axis of the crystal, while strain is applied parallel (red and blue symbols) or perpendicular (green and brown symbols) to this transport direction. Solid and open symbols correspond to the measurements taken on increasing or decreasing strain magnitude, respectively. Negative or positive ε corresponds to a compressive or tensile strain, respectively. μ FET is obtained in the linear regime of FET operation from the extended linear portions of transconductance characteristics (at |V G | >> V SD = 2 V). μ Hall is measured at V G = -40 V. Device parameters: parylene-N thickness is 1.9 μm (C i = 1.235 nF cm −2 ), channel length L = 0.3 mm, channel width W = 0.1 mm, the center-to-center distance between the voltage probes is D = 0.1 mm, and the probe width is t = 48 μm. For details on the correction factors due to the longitudinal channel shunting effect, see ref. [46]. These data reveal a large, anisotropic effect of strain on the b-axis mobility of rubrene, with the strain factor of g = 70-110 for the strain along the b-axis, and g ≈ 0 for the strain along the a-axis. that the material might be a high-mobility organic semiconductor. If this proposal is correct, an increase (decrease) of μ of the crystal should be accompanied by an increase (decrease) of the frequency f of the low-energy Raman modes that have the greatest electron-phonon coupling. Figure 3a shows the low-frequency region of the Raman spectrum of a thin rubrene crystal laminated on a semiflexible substrate (for technical details see the Experimental Section), with the main resolved peaks labeled "peaks 2-5": 75, 104.5, 118, and 139.6 cm −1 . The low-frequency cutoff of our Raman spectrometer does not allow resolving features below ≈50 cm −1 (marked by the dotted gray line). Nevertheless, peaks 2-5 represent important molecular vibrations (for DFT simulations see Section S4 in the Supporting Information). For comparison, Figure 3b shows a low-energy Raman spectrum of rubrene adopted from Girlando et al., [37] where measurements reveal the lowest-energy peak at 35 cm −1 . The dependence of the Raman peak positions on compressive and tensile strain applied along the b-axis of the crystal is shown in Figure 3c-e. All the resolved peaks show qualitatively similar behavior (for simplicity, only the data for peaks 2-4 are shown in the main text, with other peaks shown in Section S5 in the Supporting Information). It is clear from these measurements that low-energy Raman peaks shift to higher (lower) frequencies with compressive (tensile) strain, which is indicative of the hardening (softening) of the corresponding phonon modes of the crystal.
To the best of our knowledge, this is the first measurement of Raman modes of an organic semiconductor as a function of uniaxial strain. In the past, Raman spectra of organic crystals were investigated under hydrostatic pressure in naphthalene, [38,39] anthracene, [38,40] tetracene, [41] and pentacene [42] (for a review, see ref. [11]). It was found that mode frequencies increase with pressure, which is consistent with the straindependent measurements performed here. These results suggest that interaction between the molecules strengthens, as they get closer with pressure. In the past, however, these Adv. Sci. 2020, 7, 1901824   Figure 3. Low-frequency region of Raman spectra of rubrene as a function of uniaxial strain. a) Raman spectrum of a thin rubrene crystal with the resolved peaks 2-5, marked by the thin dotted lines. Our spectrometer's resolution cutoff is at ≈50 cm −1 , below which peaks cannot be resolved (indicated by the dotted gray line). b) Raman spectrum of a rubrene crystal reproduced for comparison from Girlando et al. [37] The lowest-energy mode at 35 cm −1 (peak 1) is indicated by the red arrow. c-e) The dependence of the position of Raman peaks 2 (75 cm −1 ), 3 (104.5 cm −1 ), and 4 (118 cm −1 ) on strain applied along the high-mobility b-axis of rubrene. Compressive and tensile strains correspond to the negative and positive values, respectively. The solid red line is a linear fit, and the Grüneisen parameter, G, calculated from these fits via Equation (3) is indicated.
pressure-dependent measurements were not coupled with a concurrent characterization of the charge carrier mobility. In our work, strain-dependent Raman measurements, f(ε), are correlated with the measurements of the charge carrier mobility as a function of strain, μ(ε), obtained in the same single crystals.

Discussion
The effect of strain on the frequency of Raman modes can be quantified by the so-called Grüneisen parameter, G, which was originally defined as a fractional change of the mode's frequency, δf/f, per unit of the fractional change of the sample's volume, δν/ν, occurring under a hydrostatic pressure: f . [43] In measurements under uniaxial strain, the fractional change of the sample's volume is equivalent to the fractional change of its linear dimension, provided that the Poisson effect is not too strong. This allows us to redefine the Grüneisen parameter in terms of the applied strain ε (index i denotes a particular Raman mode) The Grüneisen parameter for peaks 2-4 is calculated from our data: G 2 = 3.9 ± 0.3, G 3 = 1.1 ± 0.1, and G 4 = 1.7 ± 0.1.
Besides the consistent qualitative correlation between the effects of strain on mobility, μ(ε), and the low-energy Raman modes, f(ε), our additional goal is to quantitatively analyze this correlation. In order to perform such an analysis, one needs to know the relationship between the charge carrier mobility and the frequency of the relevant Raman modes, μ( f ). In the absence of precise theoretical relationship, below we resort to an approximation valid at high temperatures that treats molecular vibrations classically and uses the experimental temperature dependence of the charge carrier mobility in rubrene, μ(T), previously obtained via four-probe FET and Hall effect measurements, [25] and theoretically analyzed based on the model of off-diagonal thermal disorder and transient carrier localization [9,10] Here, μ 0 is a constant, and exponent γ = -1.5 fits well the high-temperature region (175-300 K) of μ(T) data relevant for the current study (Section S6, Supporting Information).
At high temperatures, energy of molecular vibrations can be treated classically, and thus the mechanical energy of a given vibrational mode can be expressed via thermal energy  ω ( ) where M is the effective mass of the molecule in this vibrational mode, ω = 2πf(ε) is the strain-dependent angular frequency of the mode, Q is the effective coordinate of the mode (with the mode's amplitude be done by extrapolating the Grüneisen parameter measured for the resolved peaks. Such an extrapolation is possible via the well-known scaling behavior of Grüneisen parameter, originally observed in pressure-variable Raman measurements of molecular solids: [43] G was found to roughly scale as 1/f i 2 with the Raman mode's frequency f i , suggesting that lower energy peaks should typically exhibit a greater relative shift with pressure. This scaling is a consequence of the coexistence of strong intramolecular and weak intermolecular forces characteristic of all molecular crystals. Based on this scaling law and G i we measured for the resolved peaks (i = 2, 3, and 4), we can approximately predict the Grüneisen parameter for the 35 cm −1 peak: These estimates are consistent with the 23% mode shift for the lowest energy peak expected from the above preliminary analysis. Of course, more rigorous work is needed on modeling the influence of strain on mobility through the vibrational frequencies of relevant Raman modes.
It must be noted that the mobility-strain relationship, μ(ε), observed in our experiment does not agree well with either of the two recent theoretical papers modeling the mobilitystrain dependence in rubrene (Gali et al., [44] and Ruggiero et al. [45] ). Furthermore, these theoretical works are not in a good agreement with each other. While Gali et al. [44] do predict an increasing mobility for the charge transport along the high-mobility axis of rubrene (here, the b-axis) with compressive strain applied either along or perpendicular to this axis, the predicted effect is smaller in magnitude compared to our experiment and almost completely lacks the anisotropy observed here. Indeed, Figures 5 and 7 of ref. [44] suggest that for the orthorhombic rubrene the b-axis mobility should increase by only about 8% per 1% of compressive strain, irrespective of the direction of strain application. On the contrary, for the same transport direction (along the highmobility b-axis), Ruggiero et al. [45] predict that the effect of strain applied along the orthogonal direction (along the lowmobility a-axis) should be much stronger than that of strain applied along the transport direction (Figure 7 of ref. [45]). For the strain applied along the transport direction (the highmobility b-axis), Ruggiero et al. [45] predict a very small effect: in pristine crystals, it has a negative sign (that is, μ slightly decreases with compressive strain), but becomes positive (while remaining very small in magnitude) with extrinsic disorder incorporated into the model. For any amount of extrinsic disorder, however, Ruggiero et al. [45] predict a very small effect of strain applied along the high-mobility b-axis, negligible compared to the effect of strain applied along the low-mobility a-axis. Given the molecular packing of rubrene, such an anisotropy is counterintuitive, it is not agreeing with the predictions of Gali et al., [44] and it is also at odds with the anisotropy of the strain effect observed in our experiment. These disagreements, both on qualitative and quantitative levels, clearly indicate that adequate theoretical description of the transport properties of organic semiconductors is still lacking, and further theoretical and experimental work in this direction is needed.

Conclusions
In conclusion, we have carried out Hall effect and low-frequency Raman measurements in a crystalline organic semiconductor rubrene as a function of uniaxial strain. Our methodology allowed accessing for the first time the intrinsic dependences of the charge carrier mobility and low-frequency molecular vibrations as a function of strain. Measurements were performed with strain applied along or perpendicular to the high-mobility transport direction of rubrene (the direction of molecular stacks), revealing an anisotropic response of the FET and Hall mobilities to strain, consistent with the molecular packing of rubrene crystals. Both the (longitudinal) FET and Hall mobilities are found to be very sensitive to strain, with a significant strain factor g ≈ 70-100, for the strain along the high-μ direction, indicating that μ can be enhanced by nearly a factor of two with a small compressive strain of 1%. Raman measurements revealed hardening (softening) of the low-frequency molecular vibrations under compressive (tensile) strain, with the corresponding Raman modes at 75, 104.5, and 118 cm −1 shifting by about 1-4% per 1% of the applied strain, which is consistent with the observed changes in the charge carrier mobility with strain. To the best of our knowledge, this work is the first direct experimental demonstration of a correlation between the lowfrequency Raman modes and the charge carrier mobility of an organic semiconductor, obtained within the same system. Accurate magneto-transport and Raman measurements, targeting the intrinsic response of charge transport and molecular vibrations to strain, are of crucial importance for future development of organic and flexible electronics. Such studies lay the foundation of strain engineering in high-performance flexible organic electronics and can provide the necessary benchmark for theoretical studies, as well as help to introduce novel highperformance devices in applications.

Experimental Section
Device Fabrication: Single crystals of rubrene with thicknesses in the range 100-300 nm were grown by physical vapor transport following published procudure. [22] Selected crystals with well-defined facets were laminated onto a PDMS-coated PET substrates of dimensions 9 × 9 mm 2 . The thicknesses of the PDMS layer and the PET substrate are 1 and 40 μm, respectively. 5 nm thick Ti adhesion underlayer and 15 nm thick Au contacts were sputtered through a contact shadow mask on PDMS prior to the crystal lamination. The devices are loaded in a custom-designed miniature strain stage that can be used to apply a calibrated uniaxial mechanical strain of up to ±5% (Figure 1e). Device wiring is done using a gold wire of 25 μm in diameter and a conducting silver paint.
Electric Measurement: All electric measurements were carried out at room temperature in a coarse dry vacuum (10 −2 Torr). Keithley K6221 current source, Hewlett-Packard HP34401A multimeter, and Keithley K6514 electrometer were used to set the excitation source-drain current, to measure the source-drain voltage, and to measure the four-probe voltage, respectively. In ac-Hall measurements, in-phase and out-of-phase components of Hall voltage, V Hall ip and V Hall op , were registered by Stanford Research SR830 lock-in amplifier tuned into the frequency of an ac magnetic field applied perpendicular to the FET's channel carrying a dc source-drain current. The net Hall voltage, V Hall , is determined as the absolute value of the vector with these components: , following the previously developed procedure. [33] The ac-B field of frequency 0.5-1.5 Hz and an r.m.s. magnitude of www.advancedscience.com 0.23 T was generated by rotating an assembly of permanent Nd magnets. Finite width of the voltage probes in the four-probe contact geometry may lead to a longitudinal channel shunting and related inaccuracies in the four-probe conductivity and mobility calculations when the conventional equations are used, as investigated in detail by Choi et al. [46] For the structures used in this work, these corrections are minor. Nevertheless, they were taken into account by the procedure developed in ref. [46]. Application of callibrated uniaxial strain was achieved by using a custom-designed strain stage that holds a semirigid substrate of unstrained length, l 0 , between its stationary and sliding jaws and bends it with a vise-like motion, producing a displacement, Δl, of the sliding jaw controlled by a micrometric screw. The length l 0 of a flat (unbent) substrate and the edge-to-edge distance (along the chord) of a bent substrate, l 0 − Δl, are used to calculate the bending radius, R, near the center of the substrate (Section S2, Supporting Information). [34] The strain in % is then calculated as ε = 2 t R × 100%, where t = 41 μm is the net thickness of the composite substrate (PET coated with PDMS).
Raman Measurements: Flexible 125 μm thick PEN substrates were used for crystal support. To block off the Raman signal from PEN, the substrates were coated with a 50 nm thick Au film by thermal evaporation. To improve the adhesion of rubrene crystals, an ≈10 nm thick Cytop was spin-coated over gold. Raman measurements were performed with HORIBA LabRAM HR Evolution spectrometer using an incident laser of wavelength 632.8 nm, which provides slightly sub bandgap excitation of rubrene, necessary to achieve a good trade-off between the fluorescence of rubrene and Raman signal. The incident laser power was kept below 8 mW to avoid sample's heating. During measurements, the laser spot was continuously rastered over the area of 10 × 10 μm 2 to distribute the incident power and collect an averaged signal. The error in the Raman peak position measurements was ≈0.1 cm −1 , and the error in the determination of Grüneisen parameter was mostly due to the standard deviation of the slope of the linear fits in Figure 3c-e.

Supporting Information
Supporting Information is available from the Wiley Online Library or from the author.