SUPPLEMENTARY INFORMATION
|
|
- Evelyn Johns
- 5 years ago
- Views:
Transcription
1 In the format provided by the authors and unedited. DOI: /NMAT4970 A map of high-mobility molecular semiconductors S. Fratini 1, S. Ciuchi 2, D. Mayou 1, G. Trambly de Laissardière 3 & A. Troisi 4 1 Institut Néel CNRS and Université Grenoble Alpes, F Grenoble, France 2 Department of Physical and Chemical Sciences University of L Aquila, Via Vetoio, I L Aquila, Italy & CNR-ISC Via dei Taurini, I Rome, Italy 3 Laboratoire de Physique Théorique et Modélisation, CNRS, Université de Cergy-Pontoise, F Cergy-Pontoise, France and 4 Dept. Chemistry, University of Liverpool, Liverpool L69 7ZD, United Kingdom Contents I. Theoretical models for charge transport 2 A. Transient localisation theory 2 Factors controlling the mobility 2 Methods 4 Relation with IPR 4 B. Non-adiabatic hopping 5 C. Semiclassical band theory 6 II. Comparison with experiments 7 A. Mobility data from the literature 7 B. Experimental validation of transient localisation theory 8 III. ab initio determination of microscopic parameters 8 A. Evaluation of J and J/J 8 B. Correlations between Js 10 C. Fluctuation time τ 11 References 11 NATURE MATERIALS 1
2 I. THEORETICAL MODELS FOR CHARGE TRANSPORT A. Transient localisation theory The ensemble of tight-binding models defined in Fig. 1a is studied via the transient localisation theory as described in [1, 2]. For each choice of the pristine parameters J a,j b,j c, the first step consists in setting up a disordered lattice by sampling the transfer integrals from a gaussian distribution of width J around these average values, and computing the time dependent quantum spread X 2 (E,t)= [ ˆX(t) ˆX(0)] 2 E for states of energy E. From this quantity on can then directly calculate the mobility of the actual (dynamically disordered) material via the Kubo formalism, using the relaxation time approximation (RTA) of Ref. [1]. The spread attained by states of energy E after a typical inter-molecular oscillation time τ is defined as l 2 τ(e) = (1/τ) X 2 (E,t)e t/τ dt. This is illustrated in Fig. 1b of the main manuscript together with the density of states (DOS) of the random lattice for a given choice of microscopic parameters. The transient localisation length is defined by the statistical average L 2 τ = l 2 τ(e), assuming Boltzmann statistics for the carriers (since we are considering holes, the dominant contributions to the average come from the top edge of the band). The mobility is then expressed as in Eq. (1) of the main text, that we report here for convenience µ = e L 2 τ k B T 2τ. (1) The dependence of the transient localisation length on the fluctuation time is illustrated in Fig. S1(a). In the explored time range, this quantity increases as a power law L 2 τ τ α with α < 1, following the well-known subdiffusive behaviour of localised systems at intermediate times. Convergence to the actual localisation length of the statically disordered lattice is very slow for 2D isotropic band structures (estimated to L in the example shown, via a logarithmic extrapolation of the diffusivity at long times [3]) but becomes much faster for anisotropic lattices where localisation effects are stronger. Factors controlling the mobility We now analyse all possible factors that control the mobility within transient localisation theory. The different parameters of the theory are: the absolute values of the transfer integrals J, the ratios between the Js along different bond directions, their relative fluctuations J/J, the fluctuation time τ, the temperature T and the lattice parameters. Effect of the Js. This is the main focus of the main manuscript and is summarised in Figs. 1c and 1d. Our results show that the mobility is very much affected by the relative values of the Js, but only weakly by their absolute values. This is in sharp contrast with both hopping-like (µ J 2 ) and semiclassical band-like transport (µ 1/m J). The weak dependence of the mobility on the magnitude of J can be traced back to one of the most salient features of high-mobility organic semiconductors: the coexistence of highly conducting extended band states together with poorly conducting tail states that are generated by disorder near the band edges [6]. As illustrated in Fig. 1b of the main manuscript, the squared localisation length varies from values around one lattice NATURE MATERIALS 2
3 a L 2 L 2 / b L 2 / Table I FIG. S1: (a) Dependence of the statistically averaged transient localisation length on the fluctuation time τ, calculated at T =0.25J for an isotropic electronic structure with J a = J b = J c = J/ 3, J =0.1eV and disorder amplitude J/J =0.5 (length in units of the lattice parameter a, time in units of h/j). The gray arrow indicates the localisation length attained in the static disorder limit and the red bar is the interval of τj/ h ratios derived from the data in Table I. The dashed line is the ratio L 2 τ /τ that enters in Eq. (1). (b) The quantity L 2 τ /τ across the ensemble, evaluated at three different values of the fluctuation time τ = 10, 20, 40 h/j, corresponding to h/τ =2.5, 5, 10 mev. unit in the tails, to values 100 times larger within the band. As J increases, because the population of electronic states is controlled via the Boltzmann distribution by the ratio k B T/J, the proportion of extended states that can be thermally populated and participate to transport is progressively reduced [1]. This causes a reduction of the statistically averaged L 2 τ, cancelling the benefits of increasing J that would be predicted in the case of semiclassical transport. Effect of the transfer integral fluctuations J/J. It is well known by now that improved materials can be designed if the dynamical off-diagonal disorder J/J is overall reduced. This was proposed theoretically [1, 2] on the basis of Eq. (1), and verified experimentally in Refs. [4, 5]. Fig. 4 of the main manuscript shows that a reduction of dynamical disorder by 20% causes an increase of the mobility by a factor of about for all compounds in the studied ensemble. Values of the transfer integral fluctuations calculated from first principles for a number of compounds are reported in Table I. For all molecular compounds studied, the J/J averaged over all bond directions lies within the range (TIPS-pentacene being an exception, with anomalously large molecular disorder). Because our focus is on how different electronic structures respond to a given level of disorder, in the manuscript we have taken a constant value J/J =0.5 throughout the ensemble. Effect of τ. The dependence of the mobility on the molecular fluctuation time is weak, as illustrated in Fig. S1(b). This happens because L 2 τ increases with fluctuation time in the subdiffusive regime of relevance here, as shown in Fig. S1(a), which partly NATURE MATERIALS 3
4 cancels the explicit factor τ in the denominator. With L 2 τ τ α and α =0.7 in the example shown, the ratio L 2 τ/τ depends on the fluctuation time as a weak power law τ 0.3. Analogous exponents arise throughout the ensemble of models except for strictly one dimensional systems, where α drops to zero. Moreover, τ does not change very much across compounds as reported in section III.B and Table I of this SI. Effect of T. The temperature dependence is illustrated in Fig. 3a of the main manuscript. Effects of the lattice parameters. Finally, Eq. (2) shows that the mobility has a trivial explicit dependence on the inter-molecular distances. Because the dependence is quadratic, doubling all bond lengths will result in a fourfold increase in mobility. In the main manuscript we have considered for simplicity a structure where all bond lengths are the same. The effects of different lattice structures can be trivially included in the calculations by restoring the appropriate lengths as provided by Xray diffraction experiments. Interestingly, all compounds studied here can be separated in two clearly distinct categories. In systems where the molecules are standing perpendicular to the high mobility plane, the denser allowed packing leads to generally shorter bond distances. This is exemplified in Table I where we report the corresponding area per molecule: pentacene (24.28 Å 2 ), C10-DNTT (22.83), C10-DNBDT (23.62). In all other compounds, the molecules are organised with the long axis lying close to the high mobility plane, leading to larger inter-molecular distances. The corresponding areas per molecule in this second class are typically larger by a factor of 2 or more: rubrene (51.84 Å 2 ), TESADT (47.37), dif-tesadt (46.34), TMTES-Pn (54.80), TIPS-Pn (58.24). To compare all systems in a single universal plot and include these trivial (albeit quantitatively important) factors, in Fig. 2 we have reported the experimental values of the mobility divided by the corresponding area per molecule, which appropriately scales as the average bond distance squared. Methods The numerical simulations of the dynamical spread X 2 (E,t) performed in this work use a real-space method based on orthogonal polynomials of the Hamiltonian. The method takes into account exactly all quantum processes, including Anderson localisation effects caused by lattice randomness, and allows to treat efficiently system sizes of 10 7 sites and above, which were not attainable with the method used in Refs. [1, 2]. Historically it has been developed in several stages and applied to various models such as amorphous systems [7, 8], quasicrystals [9, 10] and graphene-based systems (see Ref. [11] and references therein). It is applied here for the first time to organic semiconductors. The present simulations are performed on clusters consisting of up to molecules, which is sufficient for convergence at the considered levels of disorder and temperature. Relation with IPR Because the calculation of L 2 τ in large samples requires non-trivial numerical methods as described above, we also propose here a simpler approach for a faster qualitative screening of NATURE MATERIALS 4
5 FIG. S2: Comparison of the squared transient localisation length L 2 τ at τ =0.05 h/j and the static IPR of disorder-induced tail states, calculated for the ensemble of models considered in the main manuscript, indexed by the angle θ. devices, which consists in calculating the inverse participation ratio (IPR) of the parent statically disordered lattice. This is defined as IPR= ( i i ψ i 4 ) 1, with ψ the wavefunction and i the state localized on each lattice point, and can be accessed straightforwardly by performing exact diagonalisations of the ensemble of tight-binding models defined in Fig. 1a assuming a static gaussian distribution of transfer integrals of width J. The diagonalisations shown in Fig. S2 are performed on clusters of 3200 molecules, sufficiently large not to affect the IPR at the edge region of the DOS (note instead that full convergence of the IPR within the band requires much larger sample sizes). The reported IPR is the value calculated at the Fermi energy for a fixed density of 1 carrier per 125 molecules, i.e. well within the disorder-induced tails. While this quantity does not allow to make quantitative predictions for the mobility (because, at variance with the L 2 τ introduced above, it does not contain information on the dynamics of molecules), the IPR of the tail states near the band edge is conceptually related to the (dynamical) transient localisation length introduced above (gray arrow in Fig. S1 (a)) and the main trends are correctly reproduced across the ensemble, as shown in Fig. S2. B. Non-adiabatic hopping We provide here further details on the calculations reported in Fig. 2a of the main manuscript, where we computed the relative variation of the charge mobility in the case of non-adiabatic hopping between nearest neighbours, which include for example Marcus or Jortner formulations of the charge hopping rate [12]. The hopping time should be much slower than the vibrational relaxation time for any hopping theory to be consistent. As the vibrational relaxation time is typically 1psin liquids and 10 ps in organic solids [13, 14] one concludes that the fluctuations of the transfer integrals with characteristic timescale of 1 ps are faster than the hopping time. In NATURE MATERIALS 5
6 this limit the hopping rate at a constant temperature is proportional to the squared average of the hopping integrals (in all non-adiabatic theories). In systems like liquid crystals or amorphous phases where the fluctuation time can become slower than the hopping rate one needs to proceed differently [15]. For the model considered in this work with a constant standard deviation σ = J/J =0.5 we have J 2 a = J a 2 + σ 2 =1.25J 2 a (and similarly for the couplings J b and J c ). Taking the hopping rates to be proportional to the squared coupling in all directions, we now evaluate the diffusivity tensor by performing a convolution of the diffusion in the three directions, as described in the SI of Ref. [16]. In brief, for a 2D hexagonal lattice, where the three hopping vectors form angles θ 1 = 0, θ 2 = π/3,θ 3 =2π/3 with respect to the x axis, the diffusivity tensor D can be expressed in terms of the diffusion coefficient along the three hopping directions D 1, D 2, D 3 (each proportional to the squared average transfer integral) as: where D 1 = V 1 (V 1 + V 2 ) 1 V 2 (V 1 (V 1 + V 2 ) 1 V 2 + V 3 ) 1 V 3 [ ] 1/D1 0 V 1 = 0 1/ε [ ][ ][ cos θ2 sin θ V 2 = 2 1/D2 0 cos θ2 sin θ 2 cos θ 2 0 1/ε sin θ 2 ] sin θ 2 cos θ 2 [ ][ ][ cos θ3 sin θ V 3 = 3 1/D3 0 cos θ3 sin θ 3 cos θ 3 0 1/ε sin θ 3 ] sin θ 3 cos θ 3 and ɛ is an arbitrary small number. The main manuscript reports the trace of the diffusivity tensor (proportional to the mobility tensor) multiplied by an arbitrary constant. C. Semiclassical band theory For the band theory calculations presented in the text we consider the same ensemble of tight-binding models and evaluate the mobility to lowest-order in the fluctuations of the transfer integral as described in the SI of Ref. [6]. We first express the carrier-lattice interaction Hamiltonian in second quantisation as, H I = (1/N ) k,q i=a,b,c g(i) k,q c+ k+q c kx (i) q, with N the number of molecules, c + k,c k the creation and annihilation operators for carriers, x (i) q the deformation mode corresponding to a given bond direction. Straightforward algebra allows to write the interaction matrix elements for uncorrelated bond disorder as [g (i) k,q ]2 = 4(dJ/dx) 2 cos[(k q/2) δ i ], with δ i the vectors connecting nearest-neighbours as shown in Fig. 1a and dj/dx the sensitivity of the transfer integrals to inter-molecular deformations. The semiclassical mobility is expressed as µ(t )=(e/k B T ) v 2 kτ k = e nk B T vkτ 2 k e +(ɛ k µ)/k B T where τ k and v k are respectively the scattering time and the band velocity for states of momentum k, µ is the chemical potential and n = k e+(ɛk µ)/t the thermally populated carrier density (the + sign in the exponent is for holes). In the quasi-elastic limit where the k (2) NATURE MATERIALS 6
7 intermolecular vibration frequency sets the smallest energy scale in the problem, hω 0 T,J the scattering time is obtained to second order in H I as 1/τ k = 2k BT dq gk,k+q 2 δ(ɛ k ɛ k+q ). (3) hω 0 Note that in the above expression we have omitted for simplicity the correction factor from the scattering angle between incoming and outgoing states that defines the transport scattering time, (1 cos θ k,k+q ), which however only amounts to a numerical correction factor of order 1 [6]. Since the considered band structures are anisotropic we consider vk 2 and consequently µ as tensors. To appropriately compare the semiclassical results with the RTA calculations we write the classical (thermal) fluctuation of the transfer integrals as ( J) 2 =(dj/dx) 2 k B T/K using the equipartition principle, with K the spring constant, which univocally fixes dj/dx and g (i) k,q for a given value of J/J. II. COMPARISON WITH EXPERIMENTS A. Mobility data from the literature In this section we justify the selected range of experimental mobility reported in the main manuscript for each material. Given the large spread of values that can be found in the literature, we have considered only materials for which a mobility decreasing with temperature has been demonstrated and showing good reproducibility between different groups, as these are as close to the ideal intrinsic behaviour as can be achieved today. A summary is given in Table I. Rubrene. Charge mobility data from different groups are similar despite different fabrication techniques. Considering few recent measurements, the Frisbie group [17] reported mobility in the cm 2 /Vs range (for a number of samples) along the high mobility direction with an air gap connection with the gate. The Batlogg group [18] reported an identical range with cytop as dielectric material. To compare across materials we have to consider the average mobility in the high mobility plane. The anisotropy ratio (µ xx /µ yy ) for rubrene was consistently estimated around 2.7 (see Ref. [19] and references therein). Thus a reasonable experimental mobility averaged over the plane is (12.5 ± 2.5) (1+1/2.7)/2 =8.6 ± 1.7 cm 2 /Vs. This value is very similar to the first high mobility of rubrene reported by Podzorov et al.[20]. Pentacene. Although this material is less broadly used now there are also fairly consistent reports of its FET mobility. On SiO 2 gate the values reported are 1 cm 2 /Vs (Ref. [21]), in the range of cm 2 /Vs depending on the orientation [22] and 2.2 cm 2 /Vs in Ref. [23]. Solution processed pentacene devices are reported to have mobilities in the cm 2 /Vs range [24]. Similar values coupled to Hall effect measurements (proving the close-tointrinsic nature of charge transport in pentacene) are reported in [25]. A fairly representative experimental range of mobility used to compare with the theory is 1.45 ± 0.85 cm 2 /Vs. TIPS-Pentacene (6,13-bis(triisopropylsilylethynyl)pentacene). The reported values are consistent in range. Ref. [4] reports µ<1 cm 2 /Vs; Ref. [26] reports 0.7 cm 2 /Vs; Ref. [27] reports 0.202±0.012 cm 2 /Vs. A range that captures these observations is therefore 0.6±0.4 cm 2 /Vs. NATURE MATERIALS 7
8 TMTES-Pentacene (1,4,8,11-Tetramethyl-6,13-triethylsilylethynylpentacene). We use 2.5± 0.5 cm 2 /Vs, considering the range 2 3 cm 2 /Vs given in Ref. [4] and the value of 2.5 cm 2 /Vs from a separate measurement of the same authors [28]. C 10 -DNTT (2,9-di-decyl-dinaphtho-[2,3-b:20,30-f]-thieno-[3,2-b]thiophene). We have used the value (reported with no error) of 8.5 cm 2 /Vs from Ref. [29]. C 10 -DNBDT (3,11-didecyldinaphto[2,3-d:2,3 -d ]benzo[1,2-b:4,5-b ]dithiophene). The Takeya group reports 12.1 cm 2 /Vs with standard deviation 1.4 cm 2 /Vs in Ref. [30] and 9.7 cm 2 /Vs in Ref. [5]. We use the range with two standard deviations 12.1 ± 2.8 cm 2 /Vs. TESADT (5,11-bis(triethylsilylethynyl) anthradithiophene). We take the value µ =1.5±0.5 cm 2 /Vs from Ref. [4]. dif-tesadt (2,8-Difluoro-5,11-bis(triethylsilylethynyl) anthradithiophene), which is the fluorinated analogue of TESADT. We take the value µ =3.5 ± 0.5 cm 2 /Vs from Ref. [4], consistent with the one reported in Ref. [31] B. Experimental validation of transient localisation theory In Figs. 1c and 1d of the main manuscript we have illustrated which patterns of transfer integrals are beneficial to the mobility of organic semiconductors. In doing so, and to address the whole class of materials on the same footing, we have kept a number of microscopic parameters fixed across the ensemble: the absolute scale J of the transfer integrals, their fluctuations J/J with respect to the mean value, the molecular fluctuation time τ and the lattice geometry. The effects of varying the disorder strength J/J and the lattice geometry have been investigated in the main text. Concerning the remaining parameters, their limited effect on the mobility was rationalised by the arguments presented in Sec. I A. The predictive power of transient localisation theory can however be demonstrated to a higher level of accuracy by supplementing it by material-specific ab initio calculations, releasing the assumptions of a constant parameter set. To this aim we have calculated all the microscopic parameters of the theory for a number of compounds (see Table I and Section III) and correspondingly determined the quantity L 2 τ/τ that governs the mobility. On the experimental side, in order to minimize the well known variations in mobility observed in devices realised with different methods, on different substrates [32] and measured by different groups [33], we have selected a homogeneous set of materials of the same chemical class (second block in Table I), that have been purified and measured in virtually the same way by Illig and collaborators [4]. The mobilities are thought by the authors to be close to intrinsic also on the basis of Hall mobility measurements performed on two of the materials (TIPS-Pn and TMTES-Pn) in a similar device configuration [28]. As shown in Fig. S3, by combining ab initio methods and transient localisation theory one can reach an excellent quantitative agreement with the transport experiments in field effect transistors when the uncertainties inherent to the experimental methods are appropriately controlled. III. AB INITIO DETERMINATION OF MICROSCOPIC PARAMETERS A. Evaluation of J and J/J Data on the average transfer integrals used to validate the theory in the main manuscript (Figure 1d) are summarized in Table I. All transfer integrals reported in Table I have been NATURE MATERIALS 8
9 dif-tesadt TMTES-Pn TESADT TIPS-Pn FIG. S3: Test of the proportionality between measured mobility and computed L 2 τ /τ using a homogeneous set of measurements on a group of molecules of the same chemical class appeared in Ref. [4]. L 2 τ was computed using the 2D lattice geometry of the actual crystal. τ was evaluated for each given material as 0.15, 0.22, 0.18 and 0.15 ps for TIPS-Pn, TMTES-Pn, TESADT and dif-tesadt respectively. The values of J and J were computed at the same level of theory (see Table I). This calculation, unlike the spherical map in the main manuscript, takes into account the different parameters J and J/J in different bond directions. computed with the method described in [34] using the B3LYP density functional [35] and the 6-31G* basis set. The comparison with experiments presented in Fig. S3 for a homogeneous set of experimental data with additional material-specific parameters required the evaluation of the transfer integral fluctuations J in all bond directions. This was performed as described several times in the past [36, 37] by carrying out a classical MD simulation using the MM3 force field, at 300 K, and at constant volume with lattice parameters taken from the crystallographic structure. For the molecules considered in Figure S3 the following supercells were considered for TMTES-Pn and dif-tesadt, and for TESADT or TIPS-Pn. To evaluate the characteristic fluctuations time τ, the transfer integrals have been computed for MD snapshots separated by 0.2 ps (0.12 for TIPS-P) along a dynamics of 12 ps (these are reported in Figure S5). For a more accurate evaluation of J, the transfer integrals have been computed for a longer dynamics of 40 ps or more with snapshots taken every 1 ps (0.6 ps for TIPS-Pn). The transfer integral for pentacene (at the same level B3LYP/6-31G*) derives from ref. [38], while J/J was evaluated with the semiempirical method ZINDO in [36] [42]. The data of rubrene have been computed as described in ref.[37] but at the B3LYP/6-31G* level for consistency with the other data presented here. J and J were computed with the same method for C10-DNBDT and C10-DNTT using and supercells respectively and sampling the MD trajectory every 1 ps. To verify that the transfer integrals (specifically their relative magnitude and J/J) are not too dependent on the computational details we have performed test calculations on a subset of structures for all molecules computed in this work with a much larger basis set (6- NATURE MATERIALS 9
10 TABLE I: Summary of the materials considered, the experimental range of mobility (in cm 2 /Vs), the hopping integrals (in mev), the relative disorder, the value of θ used for Fig. 1d of the main text, the fluctuation period 2πτ (in ps, when available) obtained from molecular dynamics simulations and the area per molecule in the plane of high mobility (in Å 2 ). Material µ exp J a ( J a /J a ) J b ( J b /J b ) J c ( J c /J c ) θ 2πτ area Rubrene 8.6 ± (0.282) 21.2 (0.497) 21.2 (0.497) Pentacene 1.45 ± (0.327) (0.330) (0.234) C 10 DNTT (0.210) (0.647) (0.647) C 10 DNBDT 12.1 ± (0.380) 59.1 (0.592) 59.1 (0.592) TIPS-pentacene 0.6 ± (0.569) 8.6 (3.25) 10.3 (5.07) TMTES-pentacene 2.5 ± (0.264) 26.9 (0.458) 27.3 (0.355) TESADT 1.5 ± (0.516) -165 (0.351) (0.862) dif-tesadt 3.5 ± (0.60) -185 (0.261) -55.9(0.650) G*) obtaining extremely consistent results (correlation r 2 between the transfer integrals was in the range for all materials). A similar comparison with calculations performed with the 6-31G* basis set and the PBE [39] functional also yields an excellent correlation (r 2 =0.9964) with the PBE transfer integrals typically smaller by 12% for TIPS- Pn. This level of agreement between electronic structure methods is common (see e.g. ref. [40]) and guarantees that the application of this methodology is not much influenced by the inaccuracies of electronic structure calculations. B. Correlations between Js The variations of the transfer integrals involving a common molecule and due to a localized (e.g. optical) mode are in principle correlated. However, in the main manuscript we have considered a model where the transfer integrals in the different bond directions are essentially uncorrelated, which we explain in the following. Different modes may have a positive or negative correlation, as illustrated with visually simple examples in Ref. [41]. If the number of modes is very large it seems reasonable that, due to the presence of both positive and negative correlation, the overall correlation is small. To carefully check this assumption we computed the r 2 correlation coefficient for transfer integrals between pairs of molecules sharing a common molecule for the four molecules considered in Figure S3. The largest values we found were 0.092, 0.091, 0.055, 0.022, for TESADT, TIPS-Pn, dif-tesadt, TMTES- Pn, suggesting that correlations can be neglected (Figure S4 shows the most correlated transfer integral pair). A similar lack of correlation was found in the past for pentacene [36]. A sample of the instantaneous correlations for TESADT along a MD trajectory is shown in Fig. S4. NATURE MATERIALS 10
11 J 13 / mev J 12 / mev FIG. S4: Instantaneous correlations between the transfer integrals in different bond directions and involving a common molecule, calculated for TESADT from the quantity J 12 (t)j 13 (t) along a given MD trajectory. C. Fluctuation time τ On the basis of previous studies showing a distribution of low frequency modes typically peaked at hω 0 5 mev in several materials, the characteristic time of the transfer integral fluctuation was set to a constant τ =1/ω 0 =0.13 ps in the main manuscript (this value corresponds to a period of molecular oscillation 2πτ =0.82 ps). To validate this assumption we have computed a time series of transfer integrals along a short portion of an MD simulation for 4 different compounds (transfer integrals computed every 0.2 ps, sufficient to address the fluctuation spectrum in the relevant frequency range). The time evolution of the transfer integrals as well as the fluctuation power spectrum are illustrated in Fig. S5. The period 2πτ for the different materials is reported in Table I and was derived from the frequency of the highest peak in the power spectrum for each molecule (for TESADT the average frequency between the two maxima of similar height was considered). [1] Ciuchi, S. & Fratini, S., Electronic transport and quantum localization effects in organic semiconductors. Phys. Rev. B 86, (2012). [2] Fratini, S., Mayou, D. & Ciuchi, S. The transient localization scenario for charge transport in crystalline organic materials. Adv. Funct. Mater. 26, , (2016). [3] Ciuchi, S., Fratini, S. & Mayou, D. Transient localization in crystalline organic semiconductors. Phys. Rev. B 83, (R) (2011). [4] Illig, S. et al. Reducing dynamic disorder in small-molecule organic semiconductors by suppressing large-amplitude thermal motions. Nature Comms 7, (2016). [5] Kubo, T. et al. Suppressing molecular vibrations in organic semiconductors by inducing strain Nature Comms 7, (2016). [6] Fratini, S. & Ciuchi, S. Bandlike motion and mobility saturation in organic molecular semi- NATURE MATERIALS 11
12 80 60 TIPS-Pn TMTES-Pn J / mev 20 0 J / mev t / ps t / ps 0-50 TESADT dif-tesadt J / mev J / mev t / ps t / ps dif-tesadt TESADT TIPS-Pn TMTES-Pn Intensity / a.u (2πτ) -1 / ps -1 FIG. S5: (top) Time evolution of the transfer integrals along MD trajectories and (bottom) the corresponding power spectrum calculated for TIPS-Pn, TMTES-Pn, TESADT and dif-tesadt. conductors. Phys. Rev. Lett. 103, (2009). [7] Mayou, D. Calculation of the conductivity in the short-mean-free-path regime. Europhys. Lett. 6, (1988). [8] Mayou, D. & Khanna, S.N. A real-space approach to electronic transport. Journal de Physique I 5, (1995). [9] Roche, S. & Mayou, D. Formalism for the computation of the RKKY interaction in aperiodic systems. Phys. Rev. B 60, (1999). [10] Triozon, F., Vidal, J., Mosseri, R. & Mayou, D. Quantum dynamics in two- and threedimensional quasiperiodic tilings. Phys. Rev. B 65, (R) (2002). [11] Trambly de Laissardière, G. & Mayou, D. Conductivity of graphene with resonant and non- NATURE MATERIALS 12
13 resonant adsorbates Phys. Rev. Lett. 111, (2013). [12] Barbara, P. F., Meyer, T. J. & Ratner, M. A. Contemporary issues in electron transfer research. J. Phys. Chem. 100, (1996). [13] Califano, S. & Schettino, V. Vibrational Relaxation in Molecular Crystals. Int. Rev. Phys. Chem. 7, (1988). [14] Bellows, J. C. & Prasad, P. N. Dephasing times and lineswidths of optical-transitions in molecular crystals - temperature dependence of line-shapes, lineswidths, and frequencies of raman active phonons in naphthalene. J. Chem. Phys 70, (1979). [15] Kirkpatrick, J., Marcon, V., Nelson, J., Kremer, K. & Andrienko, D. Charge mobility of discotic mesophases: A multiscale quantum and classical study. Phys. Rev. Lett 98, (2007). [16] Aragó, J. & Troisi, A. Regimes of Exciton Transport in Molecular Crystals in the Presence of Dynamic Disorder. Adv. Funct. Mater. 26, (2016). [17] Xie, W. et al. High-Mobility Transistors Based on Single Crystals of Isotopically Substituted Rubrene. J. Phys. Chem. C 117, (2013). [18] Blülle, B., Häusermann, R., & Batlogg, B. Approaching the Trap-Free Limit in Organic Single-Crystal Field-Effect Transistors. Phys. Rev. Appl. 1, (2014). [19] Reese, C. & Bao, Z. High-Resolution Measurement of the Anisotropy of Charge Transport in Single Crystals. Adv. Mater. 19, (2007). [20] Podzorov, V. Sysoev, S. E. Loginova, E. Pudalov, V. M. & Gershenson, M. E. Single-crystal organic field effect transistors with the hole mobility 8 cm 2 /Vs. Appl. Phys. Lett. 83, (2003). [21] Takeyama, Y., Ono, S. & Matsumoto, Y. Organic single crystal transistor characteristics of single-crystal phase pentacene grown by ionic liquid-assisted vacuum deposition. Appl. Phys. Lett. 101, (2012). [22] Lee, J. Y., Roth, S., & Park, Y. W. Anisotropic field effect mobility in single crystal pentacene. Appl. Phys. Lett. 88, (2006). [23] Roberson, L. B. et al. Pentacene disproportionation during sublimation for field-effect transistors. J. Am. Chem. Soc. 127, (2005). [24] Kimura, Y., Niwano, M., Ikuma, N., Goushi,K. & Itaya, K. Organic Field Effect Transistor Using Pentacene Single Crystals Grown by a Liquid-Phase Crystallization Process, Langmuir 25, (2009). [25] Uemura, T., Yamagishi, M., Soeda, J., Takatsuki, Y., Okada, Y., Nakazawa, Y. and Takeya, J. Temperature dependence of the Hall effect in pentacene field-effect transistors: Possibility of charge decoherence induced by molecular fluctuations. Phys. Rev. B 85, (2012) [26] Wade, J. et al. Charge mobility anisotropy of functionalized pentacenes in organic field effect transistors fabricated by solution processing. J. Mater. Chem. C 2, (2014). [27] Keum, C.-M. et al. Topography-guided spreading and drying of 6,13- bis(triisopropylsilylethynyl)-pentacene solution on a polymer insulator for the field-effect mobility enhancement. Appl. Phys. Lett. 102, (2013). [28] Chang, J.-F. et al. Hall-Effect Measurements Probing the Degree of Charge-Carrier Delocalization in Solution-Processed Crystalline Molecular Semiconductors. [29] Hofmockel, R. et al. High-mobility organic thin-film transistors based on a small-molecule semiconductor deposited in vacuum and by solution shearing. Organ. Elect. 12, (2013). [30] Mitsui, C. et al.. High-performance solution-processable N-shaped organic semiconducting NATURE MATERIALS 13
14 materials with stabilized crystal phase. Adv. Mater. 26, 4546?4551 (2014). [31] Niazi, N. R., et al. Solution-printed organic semiconductor blends exhibiting transport properties on par with single crystals. Nature Communications 6, 8598 (2015). [32] Hulea, I. N. et al. Tunable Fröhlich Polarons in Organic Single-Crystal Transistors. Nature Mater. 5, 982 (2006). [33] McCulloch, I., Salleo, A., & Chabinyc, M. Avoid the kinks when measuring mobility. Science 352, (2106). [34] Troisi, A. & Orlandi, G. The hole transfer in DNA: calculation of electron coupling between close bases. Chem. Phys. Lett. 344, (2001). [35] Becke, A. D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys (1993). [36] Troisi A. & Orlandi G. Dynamics of the intermolecular transfer integral in crystalline organic semiconductors. J. Phys. Chem. A 110, (2006). [37] Troisi, A. Prediction of the absolute charge mobility of molecular semiconductors: the case of rubrene. Adv. Mater. 19, 2000 (2007). [38] Troisi A. & Orlandi G. Band Structure of the four pentacene polymorphs and effect on the hole mobility at low temperature. J. Phys. Chem. B 109, (2005). [39] Perdew, J. P., Burke, K. and Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett (1996). [40] Huang, J. S., Kertesz, M. Validation of intermolecular transfer integral and bandwidth calculations for organic molecular materials. J. Chem. Phys (2005). [41] Li, Y., Yi, Y., Coropceanu, V. & Brédas, J.-L. Symmetry effects on nonlocal electron-phonon coupling in organic semiconductors. Phys. Rev. B 85, (2012). [42] In the main manuscript the experimental mobilities are compared with the prediction of our theory for a generic material with hopping integrals [J a, J a, J b ] whose coordinate on the semiconductor map is expressed by the coordinate of the angle θ defined in Fig. 1. Any two hopping integrals can be interchanged and/or multiplied (both) by 1 and an angle θ can be defined for all materials for which two hopping integrals are identical or close. For pentacene we have reported the value of θ =2.538 corresponding to the point on the cut that is closest to its actual position on the spherical map. Alternatively we could have used the value of θ =2.494 corresponding to a point on the map with the corresponding transient localisation length. NATURE MATERIALS 14
6.5 mm. ε = 1%, r = 9.4 mm. ε = 3%, r = 3.1 mm
Supplementary Information Supplementary Figures Gold wires Substrate Compression holder 6.5 mm Supplementary Figure 1 Picture of the compression holder. 6.5 mm ε = 0% ε = 1%, r = 9.4 mm ε = 2%, r = 4.7
More informationarxiv: v1 [cond-mat.mtrl-sci] 29 Mar 2019
From transient localization to transport: on quantum corrections and lack thereof S. Fratini, S. Ciuchi Université Grenoble Alpes, CNRS, Grenoble INP, Institut Néel, 38 Grenoble, France Dipartimento di
More informationDynamics of Electrons at Organic/Dielectric Interfaces
Dynamics of Electrons at Organic/Dielectric Interfaces S. Fratini 1 (simone.fratini@grenoble.cnrs.fr) G. Rastelli 1,2, S. Ciuchi 2, A. F. Morpurgo 3 1 LEPES, CNRS Grenoble, France 2 INFM-CNR SMC and Dipartimento
More informationarxiv: v1 [cond-mat.mtrl-sci] 11 Dec 2013
Electronic structure and transport in approximants of the Penrose tiling Guy Trambly de Laissardière, 1 Attila Szállás, 2 and Didier Mayou 3 arxiv:1312.325v1 [cond-mat.mtrl-sci] 11 Dec 213 1 Laboratoire
More informationSupplementary Figures
Supplementary Figures 8 6 Energy (ev 4 2 2 4 Γ M K Γ Supplementary Figure : Energy bands of antimonene along a high-symmetry path in the Brillouin zone, including spin-orbit coupling effects. Empty circles
More informationarxiv: v2 [cond-mat.mtrl-sci] 25 Aug 2015
The Transient Localization Scenario for Charge Transport in Crystalline Organic Materials. S. Fratini, D. Mayou, and S. Ciuchi August 26, 2015 arxiv:1505.02686v2 [cond-mat.mtrl-sci] 25 Aug 2015 Abstract
More informationSupporting Information
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2015 Supporting Information Single Layer Lead Iodide: Computational Exploration of Structural, Electronic
More information(a) (b) Supplementary Figure 1. (a) (b) (a) Supplementary Figure 2. (a) (b) (c) (d) (e)
(a) (b) Supplementary Figure 1. (a) An AFM image of the device after the formation of the contact electrodes and the top gate dielectric Al 2 O 3. (b) A line scan performed along the white dashed line
More informationBottom-up modelling of charge transport in organic semiconductors
Bottom-up modelling of charge transport in organic semiconductors David L. Cheung Department of Chemistry & Centre for Scientific Computing University of Warwick LCOPV 2010 Boulder 7th-10th August 2010
More informationMinimal Update of Solid State Physics
Minimal Update of Solid State Physics It is expected that participants are acquainted with basics of solid state physics. Therefore here we will refresh only those aspects, which are absolutely necessary
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION doi:10.1038/nature10683 Supplementary Methods and Discussion: Possibility of Toluene Inclusion in Unit Cell: The unit cell calculations show modest, non-systematic cell volume
More informationInitial Stages of Growth of Organic Semiconductors on Graphene
Initial Stages of Growth of Organic Semiconductors on Graphene Presented by: Manisha Chhikara Supervisor: Prof. Dr. Gvido Bratina University of Nova Gorica Outline Introduction to Graphene Fabrication
More informationSupplementary Figure 1 Experimental setup for crystal growth. Schematic drawing of the experimental setup for C 8 -BTBT crystal growth.
Supplementary Figure 1 Experimental setup for crystal growth. Schematic drawing of the experimental setup for C 8 -BTBT crystal growth. Supplementary Figure 2 AFM study of the C 8 -BTBT crystal growth
More informationTriplet state diffusion in organometallic and organic semiconductors
Triplet state diffusion in organometallic and organic semiconductors Prof. Anna Köhler Experimental Physik II University of Bayreuth Germany From materials properties To device applications Organic semiconductors
More information(002)(110) (004)(220) (222) (112) (211) (202) (200) * * 2θ (degree)
Supplementary Figures. (002)(110) Tetragonal I4/mcm Intensity (a.u) (004)(220) 10 (112) (211) (202) 20 Supplementary Figure 1. X-ray diffraction (XRD) pattern of the sample. The XRD characterization indicates
More informationLecture 4: Basic elements of band theory
Phys 769 Selected Topics in Condensed Matter Physics Summer 010 Lecture 4: Basic elements of band theory Lecturer: Anthony J. Leggett TA: Bill Coish 1 Introduction Most matter, in particular most insulating
More informationSECOND PUBLIC EXAMINATION. Honour School of Physics Part C: 4 Year Course. Honour School of Physics and Philosophy Part C C3: CONDENSED MATTER PHYSICS
A11046W1 SECOND PUBLIC EXAMINATION Honour School of Physics Part C: 4 Year Course Honour School of Physics and Philosophy Part C C3: CONDENSED MATTER PHYSICS TRINITY TERM 2015 Wednesday, 17 June, 2.30
More informationFrom Order to Disorder
ORGANIC ELECTRONICS Principles, devices and applications Charge Transport D. Natali Milano, 15-18 Novembre 011 From Order to Disorder From delocalized to localized states 1 The Two-Site approximation a,v
More informationElectronic structure and properties of a few-layer black phosphorus Mikhail Katsnelson
Electronic structure and properties of a few-layer black phosphorus Mikhail Katsnelson Main collaborators: Sasha Rudenko Shengjun Yuan Rafa Roldan Milton Pereira Sergey Brener Motivation Plenty of 2D materials
More informationSupplementary Materials
Supplementary Materials Sample characterization The presence of Si-QDs is established by Transmission Electron Microscopy (TEM), by which the average QD diameter of d QD 2.2 ± 0.5 nm has been determined
More informationSemiclassical formulation
The story so far: Transport coefficients relate current densities and electric fields (currents and voltages). Can define differential transport coefficients + mobility. Drude picture: treat electrons
More informationSupplementary Figure 1. Crystal packing of pentacene.
t 3 t 4 t 1 t 2 Supplementary Figure 1. Crystal packing of pentacene. The largestholecharge transfer integrals are shown in red:t 1 = 75 mev, t 2 = 32 mev, t 3 = 20 mev, t 4 = 6 mev. Note that IRactive
More informationSupplementary Figure 1: Spin noise spectra of 55 Mn in bulk sample at BL =10.5 mt, before subtraction of the zero-frequency line. a, Contour plot of
1 Supplementary Figure 1: Spin noise spectra of 55 Mn in bulk sample at BL =10.5 mt, before subtraction of the zero-frequency line. a, Contour plot of the spin noise spectra calculated with Eq. (2) for
More informationQuantum Confinement in Graphene
Quantum Confinement in Graphene from quasi-localization to chaotic billards MMM dominikus kölbl 13.10.08 1 / 27 Outline some facts about graphene quasibound states in graphene numerical calculation of
More informationSupporting Information
Supporting Information Fluorination of Metal Phthalocyanines: Single-Crystal Growth, Efficient N-Channel Organic Field-Effect Transistors, and Structure- Property Relationships Hui Jiang 1*, Jun Ye 2,
More informationCharge Transport in Organic Crystals
Charge Transport in Organic Crystals Karsten Hannewald European Theoretical Spectroscopy Facility (ETSF) & Friedrich-Schiller-Universität Jena Karsten Hannewald (Uni Jena) www.ifto.uni-jena.de/~hannewald/
More informationorganic semiconductors Henning Sirringhaus
Charge transport physics of highmobility organic semiconductors Henning irringhaus Organic electronics tatus and opportunities OLED Existing markets Emerging applications Advanced prototypes Next generation
More informationBreakdown of the semi-classical conduction theory in approximants of the octagonal tiling
Breakdown of the semi-classical conduction theory in approximants of the octagonal tiling Guy Trambly de Laissardière, Christophe Oguey, Didier Mayou To cite this version: Guy Trambly de Laissardière,
More informationBand-like Electron Transport with Record-High Mobility in the TCNQ family
Band-like Electron Transport with Record-High Mobility in the TCNQ family Yulia Krupskaya 1 *, Marco Gibertini 2 *, Nicola Marzari 2 *, and Alberto F. Morpurgo 1 * 1 Department of Quantum Matter Physics
More informationNiS - An unusual self-doped, nearly compensated antiferromagnetic metal [Supplemental Material]
NiS - An unusual self-doped, nearly compensated antiferromagnetic metal [Supplemental Material] S. K. Panda, I. dasgupta, E. Şaşıoğlu, S. Blügel, and D. D. Sarma Partial DOS, Orbital projected band structure
More informationarxiv: v3 [cond-mat.quant-gas] 5 May 2015
Quantum walk and Anderson localization of rotational excitations in disordered ensembles of polar molecules T. Xu and R. V. Krems 1 arxiv:151.563v3 [cond-mat.quant-gas] 5 May 215 1 Department of Chemistry,
More informationCarrier dynamics of rubrene single-crystals revealed by transient broadband terahertz
Supplemental Material Carrier dynamics of rubrene single-crystals revealed by transient broadband terahertz spectroscopy H. Yada 1, R. Uchida 1, H. Sekine 1, T. Terashige 1, S. Tao 1, Y. Matsui 1, N. Kida
More informationTunable Fröhlich Polarons in Organic Single-Crystal Transistors
Tunable Fröhlich Polarons in Organic Single-Crystal Transistors I. N. Hulea, S. Fratini 2, H. Xie, C.L. Mulder, N.N. Iossad, G. Rastelli 2,3, S. Ciuchi 3, and A. F. Morpurgo Kavli Institute of Nanoscience,
More informationCHARGE TRANSPORT (Katharina Broch, )
CHARGE TRANSPORT (Katharina Broch, 27.04.2017) References The following is based on these references: K. Seeger Semiconductor Physics, Springer Verlag, 9 th edition 2004 D. Jena Charge Transport in Semiconductors
More informationCondensed matter theory Lecture notes and problem sets 2012/2013
Condensed matter theory Lecture notes and problem sets 2012/2013 Dmitri Ivanov Recommended books and lecture notes: [AM] N. W. Ashcroft and N. D. Mermin, Solid State Physics. [Mar] M. P. Marder, Condensed
More informationSUPPLEMENTARY INFORMATION
DOI: 1.138/NNANO.211.214 Control over topological insulator photocurrents with light polarization J.W. McIver*, D. Hsieh*, H. Steinberg, P. Jarillo-Herrero and N. Gedik SI I. Materials and device fabrication
More informationPresented on American Physical Society March Meeting, March 16-20, 2009, Pittsburgh, PA,USA.
Presented on American Physical Society March Meeting, March 16-20, 2009, Pittsburgh, PA,USA. A highly-ordered, high mobility organic semiconductor grown from a mesophase: A test of polaron band theory
More informationSummary lecture IX. The electron-light Hamilton operator reads in second quantization
Summary lecture IX The electron-light Hamilton operator reads in second quantization Absorption coefficient α(ω) is given by the optical susceptibility Χ(ω) that is determined by microscopic polarization
More informationChapter 3 Properties of Nanostructures
Chapter 3 Properties of Nanostructures In Chapter 2, the reduction of the extent of a solid in one or more dimensions was shown to lead to a dramatic alteration of the overall behavior of the solids. Generally,
More informationStructure and Dynamics : An Atomic View of Materials
Structure and Dynamics : An Atomic View of Materials MARTIN T. DOVE Department ofearth Sciences University of Cambridge OXFORD UNIVERSITY PRESS Contents 1 Introduction 1 1.1 Observations 1 1.1.1 Microscopic
More informationTinselenidene: a Two-dimensional Auxetic Material with Ultralow Lattice Thermal Conductivity and Ultrahigh Hole Mobility
Tinselenidene: a Two-dimensional Auxetic Material with Ultralow Lattice Thermal Conductivity and Ultrahigh Hole Mobility Li-Chuan Zhang, Guangzhao Qin, Wu-Zhang Fang, Hui-Juan Cui, Qing-Rong Zheng, Qing-Bo
More informationSUPPLEMENTARY INFORMATION
DOI: 1.138/NMAT3449 Topological crystalline insulator states in Pb 1 x Sn x Se Content S1 Crystal growth, structural and chemical characterization. S2 Angle-resolved photoemission measurements at various
More informationOptical and transport properties of small polarons from Dynamical Mean-Field Theory
Optical and transport properties of small polarons from Dynamical Mean-Field Theory S. Fratini, S. Ciuchi Outline: Historical overview DMFT for Holstein polaron Optical conductivity Transport Polarons:
More informationSTRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY
STRUCTURAL AND MECHANICAL PROPERTIES OF AMORPHOUS SILICON: AB-INITIO AND CLASSICAL MOLECULAR DYNAMICS STUDY S. Hara, T. Kumagai, S. Izumi and S. Sakai Department of mechanical engineering, University of
More informationPolaron Transport in Organic Crystals: Theory and Modelling
Polaron Transport in Organic Crystals: Theory and Modelling Karsten Hannewald Institut für Physik & IRIS Adlershof Humboldt-Universität zu Berlin (Germany) Karsten Hannewald (HU Berlin/Germany) www.ifto.uni-jena.de/~hannewald/
More informationLinear temperature dependence of electron spin resonance linewidths in La 0.7 Ca 0.3 MnO 3 and YBaMn 2 O 6
Linear temperature dependence of electron spin resonance linewidths in La 0.7 Ca 0.3 MnO 3 and YBaMn 2 O 6 Abstract D. L. Huber Department of Physics, University of Wisconsin-Madison, Madison, WI 53706
More informationMicroscopic Ohm s Law
Microscopic Ohm s Law Outline Semiconductor Review Electron Scattering and Effective Mass Microscopic Derivation of Ohm s Law 1 TRUE / FALSE 1. Judging from the filled bands, material A is an insulator.
More informationBlack phosphorus: A new bandgap tuning knob
Black phosphorus: A new bandgap tuning knob Rafael Roldán and Andres Castellanos-Gomez Modern electronics rely on devices whose functionality can be adjusted by the end-user with an external knob. A new
More informationPBS: FROM SOLIDS TO CLUSTERS
PBS: FROM SOLIDS TO CLUSTERS E. HOFFMANN AND P. ENTEL Theoretische Tieftemperaturphysik Gerhard-Mercator-Universität Duisburg, Lotharstraße 1 47048 Duisburg, Germany Semiconducting nanocrystallites like
More informationCharge Extraction from Complex Morphologies in Bulk Heterojunctions. Michael L. Chabinyc Materials Department University of California, Santa Barbara
Charge Extraction from Complex Morphologies in Bulk Heterojunctions Michael L. Chabinyc Materials Department University of California, Santa Barbara OPVs Vs. Inorganic Thin Film Solar Cells Alta Devices
More informationBefore we go to the topic of hole, we discuss this important topic. The effective mass m is defined as. 2 dk 2
Notes for Lecture 7 Holes, Electrons In the previous lecture, we learned how electrons move in response to an electric field to generate current. In this lecture, we will see why the hole is a natural
More informationSummary lecture VI. with the reduced mass and the dielectric background constant
Summary lecture VI Excitonic binding energy reads with the reduced mass and the dielectric background constant Δ Statistical operator (density matrix) characterizes quantum systems in a mixed state and
More information2) Atom manipulation. Xe / Ni(110) Model: Experiment:
2) Atom manipulation D. Eigler & E. Schweizer, Nature 344, 524 (1990) Xe / Ni(110) Model: Experiment: G.Meyer, et al. Applied Physics A 68, 125 (1999) First the tip is approached close to the adsorbate
More informationLecture 12. Electron Transport in Molecular Wires Possible Mechanisms
Lecture 12. Electron Transport in Molecular Wires Possible Mechanisms In Lecture 11, we have discussed energy diagrams of one-dimensional molecular wires. Here we will focus on electron transport mechanisms
More informationPhysics 541: Condensed Matter Physics
Physics 541: Condensed Matter Physics Final Exam Monday, December 17, 2012 / 14:00 17:00 / CCIS 4-285 Student s Name: Instructions There are 24 questions. You should attempt all of them. Mark your response
More informationElectrical Transport. Ref. Ihn Ch. 10 YC, Ch 5; BW, Chs 4 & 8
Electrical Transport Ref. Ihn Ch. 10 YC, Ch 5; BW, Chs 4 & 8 Electrical Transport The study of the transport of electrons & holes (in semiconductors) under various conditions. A broad & somewhat specialized
More informationSUPPLEMENTARY INFORMATION
Dirac electron states formed at the heterointerface between a topological insulator and a conventional semiconductor 1. Surface morphology of InP substrate and the device Figure S1(a) shows a 10-μm-square
More informationSUPPLEMENTARY INFORMATION
Supplementary Methods Materials Synthesis The In 4 Se 3-δ crystal ingots were grown by the Bridgeman method. The In and Se elements were placed in an evacuated quartz ampoule with an excess of In (5-10
More informationSUPPLEMENTARY INFORMATION
Atomic structure and dynamic behaviour of truly one-dimensional ionic chains inside carbon nanotubes Ryosuke Senga 1, Hannu-Pekka Komsa 2, Zheng Liu 1, Kaori Hirose-Takai 1, Arkady V. Krasheninnikov 2
More informationGraphene for THz technology
Graphene for THz technology J. Mangeney1, J. Maysonnave1, S. Huppert1, F. Wang1, S. Maero1, C. Berger2,3, W. de Heer2, T.B. Norris4, L.A. De Vaulchier1, S. Dhillon1, J. Tignon1 and R. Ferreira1 1 Laboratoire
More informationLecture 10 Thin Film Growth
Lecture 10 Thin Film Growth 1/76 Announcements Homework: Homework Number 2 is returned today, please pick it up from me at the end of the class. Solutions are online. Homework 3 will be set Thursday (2
More informationSoft Carrier Multiplication by Hot Electrons in Graphene
Soft Carrier Multiplication by Hot Electrons in Graphene Anuj Girdhar 1,3 and J.P. Leburton 1,2,3 1) Department of Physics 2) Department of Electrical and Computer Engineering, and 3) Beckman Institute
More informationSemiclassical Electron Transport
Semiclassical Electron Transport Branislav K. Niolić Department of Physics and Astronomy, University of Delaware, U.S.A. PHYS 64: Introduction to Solid State Physics http://www.physics.udel.edu/~bniolic/teaching/phys64/phys64.html
More informationSUPPLEMENTARY INFORMATION
Supramolecular Spin Valves M. Urdampilleta, 1 J.-P. Cleuziou, 1 S. Klyatskaya, 2 M. Ruben, 2,3* W. Wernsdorfer 1,* 1 Institut Néel, associé á l Université Joseph Fourier, CNRS, BP 166, 38042 Grenoble Cedex
More informationSUPPLEMENTARY FIGURES
1 SUPPLEMENTARY FIGURES Supplementary Figure 1: Initial stage showing monolayer MoS 2 islands formation on Au (111) surface. a, Scanning tunneling microscopy (STM) image of molybdenum (Mo) clusters deposited
More informationEE301 Electronics I , Fall
EE301 Electronics I 2018-2019, Fall 1. Introduction to Microelectronics (1 Week/3 Hrs.) Introduction, Historical Background, Basic Consepts 2. Rewiev of Semiconductors (1 Week/3 Hrs.) Semiconductor materials
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Method: Epitaxial graphene was prepared by heating an Ir(111) crystal to 550 K for 100 s under 2 x 10-5 Pa partial pressure of ethylene, followed by a flash anneal to 1420 K 1.
More informationEE143 Fall 2016 Microfabrication Technologies. Evolution of Devices
EE143 Fall 2016 Microfabrication Technologies Prof. Ming C. Wu wu@eecs.berkeley.edu 511 Sutardja Dai Hall (SDH) 1-1 Evolution of Devices Yesterday s Transistor (1947) Today s Transistor (2006) 1-2 1 Why
More informationObservation of topological surface state quantum Hall effect in an intrinsic three-dimensional topological insulator
Observation of topological surface state quantum Hall effect in an intrinsic three-dimensional topological insulator Authors: Yang Xu 1,2, Ireneusz Miotkowski 1, Chang Liu 3,4, Jifa Tian 1,2, Hyoungdo
More informationCan we find metal-insulator transitions in 2-dimensional systems?
Can we find metal-insulator transitions in 2-dimensional systems? Marcelo Kuroda Term Essay for PHYS498ESM, Spring 2004 It has been almost a quarter of a century since the belief of the non existence metallic
More informationElectronic Quantum Transport in Mesoscopic Semiconductor Structures
Thomas Ihn Electronic Quantum Transport in Mesoscopic Semiconductor Structures With 90 Illustrations, S in Full Color Springer Contents Part I Introduction to Electron Transport l Electrical conductance
More informationQuantum Condensed Matter Physics Lecture 5
Quantum Condensed Matter Physics Lecture 5 detector sample X-ray source monochromator David Ritchie http://www.sp.phy.cam.ac.uk/drp2/home QCMP Lent/Easter 2019 5.1 Quantum Condensed Matter Physics 1. Classical
More informationEECS143 Microfabrication Technology
EECS143 Microfabrication Technology Professor Ali Javey Introduction to Materials Lecture 1 Evolution of Devices Yesterday s Transistor (1947) Today s Transistor (2006) Why Semiconductors? Conductors e.g
More informationElectronic structure of correlated electron systems. Lecture 2
Electronic structure of correlated electron systems Lecture 2 Band Structure approach vs atomic Band structure Delocalized Bloch states Fill up states with electrons starting from the lowest energy No
More informationA. Optimizing the growth conditions of large-scale graphene films
1 A. Optimizing the growth conditions of large-scale graphene films Figure S1. Optical microscope images of graphene films transferred on 300 nm SiO 2 /Si substrates. a, Images of the graphene films grown
More informationElectric Field-Dependent Charge-Carrier Velocity in Semiconducting Carbon. Nanotubes. Yung-Fu Chen and M. S. Fuhrer
Electric Field-Dependent Charge-Carrier Velocity in Semiconducting Carbon Nanotubes Yung-Fu Chen and M. S. Fuhrer Department of Physics and Center for Superconductivity Research, University of Maryland,
More informationCharge mobility of discotic mesophases: a multiscale quantum/
Charge mobility of discotic mesophases: a multiscale quantum/classical study Max Planck Institute for Polymer esearch LEA meeting Fundamental and Applied Macromolecular Science: Toward ext Generation Materials
More informationOutline. Introduction: graphene. Adsorption on graphene: - Chemisorption - Physisorption. Summary
Outline Introduction: graphene Adsorption on graphene: - Chemisorption - Physisorption Summary 1 Electronic band structure: Electronic properties K Γ M v F = 10 6 ms -1 = c/300 massless Dirac particles!
More informationLattice Vibrations. Chris J. Pickard. ω (cm -1 ) 200 W L Γ X W K K W
Lattice Vibrations Chris J. Pickard 500 400 300 ω (cm -1 ) 200 100 L K W X 0 W L Γ X W K The Breakdown of the Static Lattice Model The free electron model was refined by introducing a crystalline external
More informationA Combined Optical and EPR Spectroscopy Study: Azobenzene-Based Biradicals as Reversible Molecular Photoswitches
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2017 A Combined Optical and EPR Spectroscopy Study: Azobenzene-Based Biradicals as Reversible
More informationFrom Graphene to Nanotubes
From Graphene to Nanotubes Zone Folding and Quantum Confinement at the Example of the Electronic Band Structure Christian Krumnow christian.krumnow@fu-berlin.de Freie Universität Berlin June 6, Zone folding
More informationStrong Facet-Induced and Light-Controlled Room-Temperature. Ferromagnetism in Semiconducting β-fesi 2 Nanocubes
Supporting Information for Manuscript Strong Facet-Induced and Light-Controlled Room-Temperature Ferromagnetism in Semiconducting β-fesi 2 Nanocubes Zhiqiang He, Shijie Xiong, Shuyi Wu, Xiaobin Zhu, Ming
More informationLocalization and electron-phonon interactions in disordered systems
EUROPHYSICS LETTERS 20 February 1996 Europhys. Lett., 33 (6), pp. 459-464 (1996) Localization and electron-phonon interactions in disordered systems G. Kopidakis 1, C. M. Soukoulis 1 and E. N. Economou
More informationAnalysis of the ultrafast dynamics of the silver trimer upon photodetachment
J. Phys. B: At. Mol. Opt. Phys. 29 (1996) L545 L549. Printed in the UK LETTER TO THE EDITOR Analysis of the ultrafast dynamics of the silver trimer upon photodetachment H O Jeschke, M E Garcia and K H
More informationSUPPLEMENTARY INFORMATION
Supplementary Information Anisotropic conductance at improper ferroelectric domain walls D. Meier 1,, *, J. Seidel 1,3, *, A. Cano 4, K. Delaney 5, Y. Kumagai 6, M. Mostovoy 7, N. A. Spaldin 6, R. Ramesh
More information1) Institut d Electronique Fondamentale, CNRS, Univ. Paris- Sud, Université Paris- Saclay, Bâtiment 220, Rue André Ampère, F Orsay, France
Supporting information Direct band gap germanium microdisks obtained with silicon nitride stressor layers Moustafa El Kurdi, 1 Mathias Prost, 1 Abdelhamid Ghrib, 1 Sébastien Sauvage, 1 Xavier Checoury,
More informationPhysics 211B : Problem Set #0
Physics 211B : Problem Set #0 These problems provide a cross section of the sort of exercises I would have assigned had I taught 211A. Please take a look at all the problems, and turn in problems 1, 4,
More informationPhases of Na x CoO 2
Phases of Na x CoO 2 by Aakash Pushp (pushp@uiuc.edu) Abstract This paper deals with the various phases of Na x CoO 2 ranging from charge ordered insulator to Curie-Weiss metal to superconductor as the
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Aharonov-Bohm interference in topological insulator nanoribbons Hailin Peng 1,2, Keji Lai 3,4, Desheng Kong 1, Stefan Meister 1, Yulin Chen 3,4,5, Xiao-Liang Qi 4,5, Shou- Cheng
More informationA test of the method of images at the surface of molecular materials
JOURNAL OF CHEMICAL PHYSICS VOLUME 120, NUMBER 8 22 FEBRUARY 2004 A test of the method of images at the surface of molecular materials A. Eilmes Faculty of Chemistry, Jagiellonian University, Ingardena
More informationSupplementary Figure 1: A potential scheme to electrically gate the graphene-based metamaterial. Here density. The voltage equals, where is the DC
Supplementary Figure 1: A potential scheme to electrically gate the graphene-based metamaterial. Here density. The voltage equals, where is the DC permittivity of the dielectric. is the surface charge
More informationSUPPLEMENTARY INFORMATION
A Stable Three-dimensional Topological Dirac Semimetal Cd 3 As 2 Z. K. Liu, J. Jiang, B. Zhou, Z. J. Wang, Y. Zhang, H. M. Weng, D. Prabhakaran, S. -K. Mo, H. Peng, P. Dudin, T. Kim, M. Hoesch, Z. Fang,
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationLECTURE 11: Monte Carlo Methods III
1 LECTURE 11: Monte Carlo Methods III December 3, 2012 In this last chapter, we discuss non-equilibrium Monte Carlo methods. We concentrate on lattice systems and discuss ways of simulating phenomena such
More informationSupplementary Figures
Supplementary Figures Supplementary Figure 1 Measured versus calculated optical transitions in the CPX. The UV/Vis/NIR spectrum obtained experimentally for the 1:1 blend of 4T and F4TCNQ (red curve) is
More informationChapter 12: Semiconductors
Chapter 12: Semiconductors Bardeen & Shottky January 30, 2017 Contents 1 Band Structure 4 2 Charge Carrier Density in Intrinsic Semiconductors. 6 3 Doping of Semiconductors 12 4 Carrier Densities in Doped
More informationLimit of the electrostatic doping in two-dimensional electron gases of LaXO 3 (X = Al, Ti)/SrTiO 3
Supplementary Material Limit of the electrostatic doping in two-dimensional electron gases of LaXO 3 (X = Al, Ti)/SrTiO 3 J. Biscaras, S. Hurand, C. Feuillet-Palma, A. Rastogi 2, R. C. Budhani 2,3, N.
More informationSpins and spin-orbit coupling in semiconductors, metals, and nanostructures
B. Halperin Spin lecture 1 Spins and spin-orbit coupling in semiconductors, metals, and nanostructures Behavior of non-equilibrium spin populations. Spin relaxation and spin transport. How does one produce
More informationSupplementary Figure S1. AFM images of GraNRs grown with standard growth process. Each of these pictures show GraNRs prepared independently,
Supplementary Figure S1. AFM images of GraNRs grown with standard growth process. Each of these pictures show GraNRs prepared independently, suggesting that the results is reproducible. Supplementary Figure
More informationSupporting information. Realizing Two-Dimensional Magnetic Semiconductors with. Enhanced Curie Temperature by Antiaromatic Ring Based
Supporting information Realizing Two-Dimensional Magnetic Semiconductors with Enhanced Curie Temperature by Antiaromatic Ring Based Organometallic Frameworks Xingxing Li and Jinlong Yang* Department of
More information