## Observation of vacancy-induced suppression of electronic cooling in defected graphene

Qi Han, ∗ Yi Chen, ∗ Gerui Liu, Dapeng Yu, and Xiaosong Wu † State Key Laboratory for Artificial Microstructure and Mesoscopic Physics, Peking University, Beijing 100871, China

Collaborative Innovation Center of Quantum Matter, Beijing 100871, China

## Abstract

Previous studies of electron-phonon interaction in impure graphene have found that static disorder can give rise to an enhancement of electronic cooling. We investigate the effect of dynamic disorder and observe over an order of magnitude suppression of electronic cooling compared with clean graphene. The effect is stronger in graphene with more vacancies, confirming its vacancyinduced nature. The dependence of the coupling constant on the phonon temperature implies its link to the dynamics of disorder. Our study highlights the effect of disorder on electron-phonon interaction in graphene. In addition, the suppression of electronic cooling holds great promise for improving the performance of graphene-based bolometer and photo-detector devices.

## INTRODUCTION

In recent years, there has been considerable interest in utilizing graphene as photodetectors[1-8]. Most of these detectors are based on a hot electron effect, i.e. the electronic temperature being substantially higher than the lattice temperature. Two properties of graphene strongly enhance the effect. First, low carrier density gives rise to a very small electron specific heat. Second, weak electron-phonon (e-p) interaction reduces the heat transfer from the electron gas to the lattice. Thus, it is of practical interest to understand the e-p interaction in graphene. Both theoretical and experimental efforts have been devoted to this topic. Earlier work was mainly focused on clean graphene and considered the Dirac spectrum of electrons[9-15]. As the important role of impurities in electronic transport has been revealed, its effects on the e-p interaction began to draw attention[16-18]. For instance, due to the chiral nature of electrons, long range and short range potentials scatter electrons differently in graphene[19-21]. Recently, a strong enhancement of electronic cooling via e-p interaction in presence of short range disorder has been predicted[18]. This is achieved via a so-called supercollision process. When the carrier density is low, the Bloch-Gr¨neisen u temperature T BG can be quite small. Since T BG sets the maximum wave vector of phonons that can exchange energy with electrons, when T BG &lt; T , only a portion of phonons can contribute to the energy relaxation. Interestingly, in presence of short range potentials, the theory has found that a disorder-assisted scattering process can occur, in which all available phonons are able to participate. As a result, the energy relaxation is strongly enhanced. Shortly, two experiments confirmed the supercollision[22, 23], although long range potential scattering usually dominates in such samples[24, 25]. In the case of long range potentials, Chen and Clerk have also predicted an increase of electronic cooling at low temperature for weak screening[17]. Note that besides the different potential profiles, e.g. long range or short range, disorder can be static or dynamic. Despite these studies, in which only static disorder was considered, the dynamics of disorder has not been addressed.

Here, we present an experimental investigation of the effect of vacancy on electronic cooling in both monolayer and bilayer defected graphene. In contrast to typical scattering potentials previously treated in theories or encountered in experiments, which are static, vacancies in our defected graphene are dragged by phonons, hence highly dynamic. By studying the nonlinear electric transport of defected graphene, a strong suppression of e-

p energy relaxation, instead of an enhancement in the case of static potentials, has been observed. The more disordered the graphene film, the stronger the suppression is. Our work provides new experimental insight on the effect of scattering potential on e-p interaction. Moreover, the suppression suggests that the performance of graphene hot electron photodetectors can be further improved by introducing vacancies.

## EXPERIMENT

In this work, we have investigated four exfoliated graphene samples on Si/SiO 2 substrates. Thickness of all the monolayer (SM1 and SM2) and bilayer (SB1, SB2) samples were estimated by optical contrast and confirmed by Raman spectroscopy[26]. Graphene flakes were patterned into ribbons, using e-beam lithography. 5 nm Ti/80 nm Au were e-beam deposited, followed by lift-off to form electrodes. Typical sample geometry can be seen in the inset of Fig. 1a. In order to introduce vacancies, samples were then loaded into a Femto plasma system and subject to Argon plasma treatment for various periods (from 1 to 5 s)[27]. Four-probe electrical measurements were carried out in a cryostat using a standard lock-in technique. Room temperature π -filters were used to avoid heating of electrons by radio frequency noise. Information for four samples are summarized in Table I.

## RESULTS AND DISCUSSION

Previously, we have already demonstrated a hot electron bolometer based on disordered graphene[28]. It has been shown that the divergence of the resistance at low temperature can be utilized as a sensitive thermometer for electrons. By applying Joule heating, the energy transfer rate between the electron gas and the phonon gas can be obtained. The same method has been employed in this work. As showing in Fig. 1a, the resistance of defected graphene exhibits a sharp increase as the temperature decreases. The divergence becomes stronger as one approaches the CNP. The R -T behavior can be well fitted to variable range hopping transport, described as R ∝ exp[( T /T 0 ) 1 3 / ][29]. Here, the characteristic temperature T 0 = 12 [ / πk B ν E ( F ) ξ 2 ], with k B the Boltzmann constant, ν E ( F ) the density of states at the Fermi level E F , and ξ the localization length. By fitting to this formula, the localization length ξ is determined. It is employed as a measure of the degree of disorder. ξ

near the CNP for all samples are listed in Table I.

In the steady state of Joule heating, the electron cooling power equals to the heating power. The corresponding thermal model is sketched in Fig. 1c. Two thermal energy transfer pathways are indicated, i.e. via electron diffusion into electrodes or e-p interaction into the lattice. In our strongly disordered graphene, the former is significantly suppressed due to a very low carrier diffusivity. It has been found that e-p interaction dominates the energy dissipation in such devices[28]. Then, the electronic temperature can be directly inferred from the resistance. Furthermore, it is estimated that the thermal conductance between the graphene lattice and the substrate is much higher than that due to e-p interaction. Thus, the phonon temperature T ph is approximately equal to the substrate temperature T [5, 22, 30]. Under these conditions, the energy balance at the steady state of Joule heating can be written as

$$P = A ( T _ { \text{e} } ^ { \delta } - T _ { \text{ph} } ^ { \delta } )$$

where P is the Joule Heating power, A is the coupling constant and T e is the electronic temperature. δ ranges from 2 to 6, depending on the detail of the e-p scattering process[12].

Upon Joule heating, the electronic temperature is raised, leading to decrease of the resistance, depicted in Fig. 1b. Based on the resistance as a function of temperature, we obtain the P -T e relation at different carrier densities, plotted in the insets of Fig. 2. P is also plotted against T 3 e -T 3 ph . The linear behavior agrees well with Eq. (1) with δ = 3 for both monolayer and bilayer graphene at all carrier densities. It has been theoretically shown that both clean monolayer and bilayer graphene can be described by Eq. (1) with δ = 4 at low temperature [9, 12]. In presence of disorder, e-p interaction is enhanced and δ is reduced to 3[17, 18]. δ obtained in our result is consistent with these theories, indicating the effect of defects. T 3 dependence has also been reported in some other experiments. In the following, we will compare our results in detail with previous theoretical and experimental results.

The e-p interaction is usually considered in two distinct regimes, high temperature and low temperature. In normal metals, Debye temperature θ D demarcates two regimes. Below θ D , the phase space of available phonons increases with temperature, while it becomes constant above it(all modes are excited). In graphene, because of its low carrier density, the BlochGr¨ uneisen temperature T BG becomes the relevant characteristic temperature. It is defined as 2 k B T BG = 2 hck F . Here k B is the Boltzmann constant, h the Plank constant, c the sound velocity of graphene and k F the Fermi wave vector. T BG stems from the momentum

conservation in e-p scattering. Because of it, when T ph &gt; T BG , only a portion of phonons can participate in the process[31]. Considering the band structure of graphene, we have T BG = 2( c/v F ) E /k F B in monolayer graphene and T BG = 2( c/v F ) √ γ E 1 F /k B in bilayer graphene[12]. Here v F ≈ 10 6 m/s is the Fermi velocity, c ≈ 2 × 10 4 m/s and γ 1 ≈ 0 4 eV is the interlayer . coupling coefficient. Taking into account a residual carrier density n 0 ≈ 4 × 10 11 cm 2 due to charge puddles[32, 33], it can be readily estimated that even at the CNP, T BG &gt; 34 K. It is much higher than T e = 1 5 K in our experiment. . Consequently, we are well in the low temperature regime.

In the low temperature regime, the whole population of phonons can interact with electrons. Thus, the disorder-assisted supercollision is negligible[18], which rules out it as the origin of the observed T 3 dependence. It has been theoretically shown that in the case of weak screening, static charge impurities leads to enhanced e-p cooling power over clean graphene and δ = 3[17]. For comparison, we plot our data, the theoretical cooling power of clean graphene in Fig. 3. The theoretical prediction of the cooling power per unit area in clean monolayer graphene is [12]

$$P _ { \text{clean} } = \frac { \pi ^ { 2 } D ^ { 2 } E _ { \text{F} } k _ { B } ^ { 4 } } { 1 5 \rho \hbar { ^ { 5 } } v _ { \text{F} } ^ { 3 } c ^ { 3 } } ( T _ { e } ^ { 4 } - T _ { \text{ph} } ^ { 4 } )$$

where ρ ≈ 0 76 . × 10 -6 kg/m 2 is the mass density of graphene and D is the deformation potential chosen as a common value 18 eV [23, 34, 35] (this choice will be discussed later). The theoretical cooling power P clean as a function of the carrier density and electron temperature is depicted as a transparent surface (with T ph =1.5 K) in Fig. 3a. It can be clearly seen that the cooling power of our disordered samples SM1 and SM2 (green and blue lines) are well below the surface at all carrier densities. For comparison, we also plot the data from two other experiments in which T 3 -dependence were observed at low temperatures[35, 36]. These results (with similar T ph ) are either on or above the surface. The suppression of the cooling in Fig. 3a is considerable. For instance, at n = 4 × 10 11 cm 2 and T e = 20 K, the theory predicts P clean =4.7 nW/ µ m . In Ref. 36 the cooling power was found to be 27 nW/ 2 µ m . 2 In sharp contrast, our experiment gives a cooling power of 0.33 nW/ µ m for SM1, over an 2 order of magnitude lower than that in clean graphene. For the more disordered sample, SM2, it is even smaller.

Similar suppression occurs in bilayer graphene samples, too. The cooling power per unit

area in clean bilayer graphene is given by [12]

$$P _ { \text{clean} } = \frac { \pi ^ { 2 } D ^ { 2 } \gamma _ { 1 } k _ { B } ^ { 4 } } { 6 0 \rho \hbar { ^ { 5 } } v _ { \text{F} } ^ { 3 } c ^ { 3 } } \sqrt { \frac { \gamma _ { 1 } } { E _ { \text{F} } } } ( T _ { \text{e} } ^ { 4 } - T _ { \text{ph} } ^ { 4 } )$$

Fig. 3a shows the plot of Eq. (3), the cooling power of the bilayer samples SB1, SB2 and the data from Ref. 36. Although not as pronounced as monolayer graphene, our data still below the theoretical surface. The weaker suppression may result from the fact that the bottom layer of bilayer graphene has experienced less damage by our low energy plasma than the top one[27]. Therefore, this less disordered layer provides a channel of substantial cooling.

The e-p coupling strength depends on the deformation potential D , which characterizes the band shift upon lattice deformation[37-39]. For the theoretical cooling power surface in Fig. 3, we use D = 18 eV. Note that D for graphene ranges from 10 to 70 eV in various experiments, but 18 eV is the most common value for graphene[35]. If the suppression is due to an over-estimated D , to account for the small cooling power, one would require D to be only about 5 eV, one-half of the lowest value reported. Therefore, we believe that the suppression cannot be explained by a small D .

By linear fits of P versus T 3 e -T 3 ph , the coupling constant A can be obtained. In Fig. 4, A is plotted as a function of carrier density n . A for all samples decreases when approaching the CNP. This is because fewer carriers at Fermi level could contribute to total cooling power of the sample.

We now take a look at the dependence of the coupling constant on the degree of disorder. As listed in Table I, the samples have been subject to various periods of plasma treatment. Consequently, the degree of disorder is different, indicated by the localization length ξ . For instance, ξ for SM1 and SM2 is 156 nm and 21 nm, respectively. As plotted in Fig. 4a, the coupling constant A of the less disordered SM1 is only about one-third of the value for the more disordered SM2. The dependence of A on ξ is consistent with the suppression of the e-p scattering by disorder. For the two bilayer samples, SB1 and SB2, the localization lengths are close. The n dependence of A for both samples aligns reasonably well and is consistent with the monolayer samples, see Fig. 4b.

The Joule heating experiment has also been carried out at different phonon temperatures T ph . In Fig. 4c, the coupling constant A is plotted as a function of T ph . Usually, A is independent of T ph , which is actually seen at low temperature for SB1. However, as the temperature goes above 7 K, A is enhanced. Later, we will show that the unexpected

T -dependence is likely related to the dynamic nature of vacancies.

At first glance, the suppression of electronic cooling by vacancies seems surprising, in that previous theories have predicted that disorder would enhance the cooling[17, 18]. Most of earlier experimental results have confirmed the enhancement[22, 23, 35, 36]. However, there is a key difference between those earlier studies and ours. In the former, disorder is theoretically considered to be static. This is indeed true in other experimental work, in which the dominant disorder is due to charge impurities[24, 25]. However, in our samples, the dominant disorder is vacancies, which are completely dragged by phonons. The effect of disorder on the e-p interaction has been studied in disordered metals and found to depend on the character of disorder[40-43]. In the case of static disorder, diffusive motion of electrons increases the effective interacting time between an electron and a phonon, leading to an enhancement of interaction. However, dynamic disorder modifies the quantum interference of scattering processes[41]. As a result, the interaction is suppressed, in accordance with the famous Pippard's inefficient condition[44]. It is reasonable to believe that the observed suppression results from dynamic disorder, vacancies. Furthermore, since the dynamics of disorder apparently depends on T ph , the dependence of the coupling constant A on the phonon temperature T ph is then conceivable. As described in Schmid's theory[40, 41], the e-p scattering is suppressed due to strong disorder. The resultant energy relaxation rate τ -1 e-p is of the order of ( q T l ) τ -1 0 where τ -1 0 ∝ T 3 is the relaxation rate in pure material, q T is the wave vector of a thermal phonon and l is the mean free path. As q T ∝ T ph , the relaxation rate increases with T ph , in agreement with our result.

It is also worthy to note that charge impurities are long range potentials that preserve the sublattice symmetry. This is in contrast to vacancies, which are short range potentials and break the sublattice symmetry. The theory for supercollision models disorder as short range potential[18], while in Ref. 17, disorder potential is long-ranged. This character of disorder strongly affects scattering of chiral electrons in graphene. Our samples represent a graphene system that is quite different from what was commonly seen, in that dynamic and short-ranged potentials dominate. Therefore, the quantitative understanding of our experimental results, including the power index δ , relies on future theory that takes both the dynamics and the symmetry of disorder into account.

## CONCLUSION

In conclusion, we have observed significant suppression of electronic cooling in defected graphene. The cooling power of both monolayer and bilayer graphene samples show T 3 e dependence, consistent with disorder-modified electron-phonon coupling in graphene [17, 18]. However, the magnitude of the cooling power is over an order of magnitude smaller than that of clean graphene predicted by theory[9, 12] and also less than other experiments [35, 36]. The more disordered a graphene film is, the lower cooling power is observed, confirming the effect of disorder. The suppression of electronic cooling is attributed to the dynamic nature of vacancies, which has not been studied in graphene. This effect can be utilized to further improve the performance of graphene-based bolometer and photo-detector devices.

This work was supported by National Key Basic Research Program of China (No. 2012CB933404, 2013CBA01603) and NSFC (project No. 11074007, 11222436, 11234001).

∗ These two authors contributed equally.

† xswu@pku.edu.cn

- [1] X. Xu, N. M. Gabor, J. S. Alden, A. M. van der Zande, and P. L. McEuen, Nano Lett. 10 , 562 (2009).
- [2] N. M. Gabor, J. C. W. Song, Q. Ma, N. L. Nair, T. Taychatanapat, K. Watanabe, T. Taniguchi, L. S. Levitov, and P. Jarillo-Herrero, Science 334 , 648 (2011).
- [3] N. G. Kalugin, L. Jing, W. Bao, L. Wickey, C. D. Barga, M. Ovezmyradov, E. A. Shaner, and C. N. Lau, Appl. Phys. Lett. 99 , 013504 (2011).
- [4] D. Sun, G. Aivazian, A. M. Jones, J. S. Ross, W. Yao, D. Cobden, and X. Xu, Nat Nano 7 , 114 (2012).
- [5] J. Yan, M.-H. Kim, J. A. Elle, A. B. Sushkov, G. S. Jenkins, M. H. M., M. S. Fuhrer, and H. D. Drew, Nat Nano 7 , 472 (2012).
- [6] H. Vora, P. Kumaravadivel, B. Nielsen, and X. Du, Appl. Phys. Lett. 100 , 153507 (2012).
- [7] K. Yan, D. Wu, H. Peng, L. Jin, Q. Fu, X. Bao, and Z. Liu, Nat Commun 3 , 1280 (2012).
- [8] X. Cai, A. B. Sushkov, R. J. Suess, G. S. Jenkins, J. Yan, T. E. Murphy, H. D. Drew, and M. S. Fuhrer (2013), arXiv: 1305.3297.

- [9] S. S. Kubakaddi, Phys. Rev. B 79 , 075417 (2009).
- [10] R. Bistritzer and A. H. MacDonald, Phys. Rev. Lett. 102 , 206410 (2009).
- [11] W.-K. Tse and S. Das Sarma, Phys. Rev. B 79 , 235406 (2009).
- [12] J. K. Viljas and T. T. Heikkil¨, Phys. Rev. B a 81 , 245404 (2010).
- [13] A. C. Betz, F. Vialla, D. Brunel, C. Voisin, M. Picher, A. Cavanna, A. Madouri, G. F`ve, e J.-M. Berroir, B. Pla¸ais, et al., Phys. Rev. Lett. c 109 , 056805 (2012).
- [14] A. M. R. Baker, J. A. Alexander-Webber, T. Altebaeumer, and R. J. Nicholas, Phys. Rev. B 85 , 115403 (2012).
- [15] A. M. R. Baker, J. A. Alexander-Webber, T. Altebaeumer, S. D. McMullan, T. J. B. M. Janssen, A. Tzalenchuk, S. Lara-Avila, S. Kubatkin, R. Yakimova, C.-T. Lin, et al., Phys. Rev. B 87 , 045414 (2013).
- [16] F. T. Vasko and V. V. Mitin, Phys. Rev. B 84 , 155445 (2011).
- [17] W. Chen and A. A. Clerk, Phys. Rev. B 86 , 125443 (2012).
- [18] J. C. W. Song, M. Y. Reizer, and L. S. Levitov, Phys. Rev. Lett. 109 , 106602 (2012).
- [19] E. McCann, K. Kechedzhi, V. I. Fal'ko, H. Suzuura, T. Ando, and B. L. Altshuler, Phys. Rev. Lett. 97 , 146805 (2006).
- [20] K. Nomura and A. H. MacDonald, Phys. Rev. Lett. 98 , 076602 (2007).
- [21] X. S. Wu, X. B. Li, Z. M. Song, C. Berger, and W. A. de Heer, Phys. Rev. Lett. 98 , 136801 (2007).
- [22] A. C. Betz, S. H. Jhang, E. Pallecchi, R. Ferreira, G. Feve, J.-M. Berroir, and B. Placais, Nat Phys 9 , 109 (2013).
- [23] M. W. Graham, S.-F. Shi, D. C. Ralph, J. Park, and P. L. McEuen, Nat Phys 9 , 103 (2013).
- [24] S. Adam, E. H. Hwang, V. M. Galitski, and S. Das Sarma, Proc. Natl. Acad. Sci. USA 104 , 18392 (2007).
- [25] J. H. Chen, C. Jang, S. Adam, M. S. Fuhrer, E. D. Williams, and M. Ishigami, Nat Phys 4 , 377 (2008).
- [26] A. C. Ferrari, J. C. Meyer, V. Scardaci, C. Casiraghi, M. Lazzeri, F. Mauri, S. Piscanec, D. Jiang, K. S. Novoselov, S. Roth, et al., Phys. Rev. Lett. 97 , 187401 (2006).
- [27] J. Chen, T. Shi, T. Cai, T. Xu, L. Sun, X. Wu, and D. Yu, Appl. Phys. Lett. 102 , 103107 (2013).

- [28] Q. Han, T. Gao, R. Zhang, Y. Chen, J. Chen, G. Liu, Y. Zhang, Z. Liu, X. Wu, and D. Yu, Sci. Rep. 3 , 3533 (2013).
- [29] N. Mott, J. Non-Cryst. Solids 1 , 1 (1968).
- [30] I. V. Borzenets, U. C. Coskun, H. T. Mebrahtu, Y. V. Bomze, A. I. Smirnov, and G. Finkelstein, Phys. Rev. Lett. 111 , 027001 (2013).
- [31] M. S. Fuhrer, Physics 3 , 106 (2010).
- [32] Q. Li, E. H. Hwang, and S. Das Sarma, Phys. Rev. B 84 , 115442 (2011).
- [33] Y. Zhang, V. W. Brar, C. Girit, A. Zettl, and M. F. Crommie, Nat Phys 5 , 722 (2009).
- [34] J.-H. Chen, C. Jang, S. Xiao, M. Ishigami, and M. S. Fuhrer, Nat Nano 3 , 206 (2008).
- [35] K. C. Fong, E. E. Wollman, H. Ravi, W. Chen, A. A. Clerk, M. D. Shaw, H. G. Leduc, and K. C. Schwab, Phys. Rev. X 3 , 041008 (2013).
- [36] R. Somphonsane, H. Ramamoorthy, G. Bohra, G. He, D. K. Ferry, Y. Ochiai, N. Aoki, and J. P. Bird, Nano Lett. 13 , 4305 (2013).
- [37] J. Bardeen and W. Shockley, Phys. Rev. 80 , 72 (1950).
- [38] C. Herring and E. Vogt, Phys. Rev. 101 , 944 (1956).
- [39] H. Suzuura and T. Ando, Phys. Rev. B 65 , 235412 (2002).
- [40] A. Schmid, Zeitschrift f¨r Physik u 259 , 421 (1973).
- [41] A. Sergeev and V. Mitin, Phys. Rev. B 61 , 6041 (2000).
- [42] J. J. Lin and J. P. Bird, J. Phys.: Condens. Matter 14 , R501 (2002).
- [43] Y. L. Zhong, J. J. Lin, and L. Y. Kao, Phys. Rev. B 66 , 132202 (2002).
- [44] A. Pippard, Philosophical Magazine Series 7 46 , 1104 (1955).

## FIGURES

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0167v1-figure-955f588e-acd1-4c70-a027-93c7224c4564.jpg)

FIG. 1. Resistance of defected graphene. (a) Temperature dependence of resistance in sample SM1 at different gate voltages, showing divergence at low temperature. Inset: Optical micrograph of a typical device configuration. (b) Resistance of SM1 as a function of Joule heating current at different gate voltages at T = 1 5 K. The CNP is at 14.5 V. (c)Thermal model for the structure. . The pathways of heat dissipation are indicated by red arrows.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0167v1-figure-990d788c-e6cd-4332-b259-acf36c2a514c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0167v1-figure-33afdcda-d414-4402-8b79-00bdd7c59aab.jpg)

FIG. 2. Cooling power of monolayer and bilayer defected graphene. (a)(b) Cooling power P against T 3 e -T 3 ph shows a linear dependence for both monolayer and bilayer samples. Inset: P versus T e .

FIG. 3. Suppression of electronic cooling in defected graphene. (a) Cooling power of clean monolayer graphene is depicted as a transparent surface, with a logarithmic z-axis scale and as a function of T e and n . P for SM1 and SM2 are over an order of magnitude smaller than clean graphene at all carrier densities, while the data from others' work is either on or above the surface. n 0 is chosen as 4 × 10 11 cm 2 to account for charge puddles near CNP. T ph is 1.5 K in SM1 and the theoretical surface, 7K in SM2, 0.8 K in the data from Ref. 35 and 1.8 K in the data from Ref. 36. (b) Similar suppression is observed in bilayer defected graphene samples. n 0 is chosen as 4 × 10 11 cm . 2 T ph is 1.5K in SB1, SB2 and the theoretical surface, and 1.8 K in the data from Ref. 36.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0167v1-figure-affe7e01-09c1-445a-af1a-66d42a2bf641.jpg)

FIG. 4. Coupling constant A . (a) Extracted coupling constant A as a function of carrier density n in monolayer samples. The more defective sample, SM1, exhibits a smaller coupling constant. (b) Dependence of A on carrier density n in bilayer graphene samples. The curves for two samples with similar degree of disorder align reasonably well. (c) Dependence of A on phonon temperature T ph . The data of SM2 are plotted with respect to the right y-axis.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0167v1-figure-ed775671-7b82-4e41-b516-74993d81caec.jpg)

## TABLES

TABLE I. Sample information of four investigated devices. Different Ar gas flow rates and plasma treatment times have been applied to produce different amount of vacancies. V CNP is the charge neutrality point (CNP) of samples and ξ is the localization length near the CNP.

| Devices   |   Length( µ m) |   Width( µ m) |   Ar flow rate(sccm) |   Plasma treatment period(s) |   V CNP (V) |   ξ (nm) |
|-----------|----------------|---------------|----------------------|------------------------------|-------------|----------|
| SM1       |            2   |           3   |                    3 |                          1   |        14.5 |      156 |
| SM2       |            6.7 |           2.7 |                    4 |                          3   |        30   |       21 |
| SB1       |            3   |           2.7 |                    4 |                          3.5 |        70   |       50 |
| SB2       |            6   |           2.7 |                    4 |                          5   |        57   |       54 |
Qi Han Yi Chen Gerui Liu Dapeng Yu Xiaosong Wu
Observation of vacancy-induced suppression of electronic cooling in defected graphene
Previous studies of electron-phonon interaction in impure graphene have found
that static disorder can give rise to an enhancement of electronic cooling. We
investigate the effect of dynamic disorder and observe over an order of
magnitude suppression of electronic cooling compared with clean graphene. The
effect is stronger in graphene with more vacancies, confirming its
vacancy-induced nature. The dependence of the coupling constant on the phonon
temperature implies its link to the dynamics of disorder. Our study highlights
the effect of disorder on electron-phonon interaction in graphene. In addition,
the suppression of electronic cooling holds great promise for improving the
performance of graphene-based bolometer and photo-detector devices.