Abstract

We investigate the energy relaxation of hot carriers in a CVD-grown graphene device with a top h-BN layer by driving the devices into the nonequilibrium regime. By using the magnetic field dependent conductance fluctuations of our graphene device as a self-thermometer, we can determine the effective carrier temperature at various driving currents while keeping the lattice temperature fixed. Interestingly, it is found that is proportional to I, indicating little electron-phonon scattering in our device. Furthermore the average rate of energy loss per carrier is proportional to (), suggesting the heat diffusion rather than acoustic phonon processes in our system. The long energy relaxation times due to the weak electron-phonon coupling in CVD graphene capped with h-BN layer as well as in exfoliated multilayer graphene can find applications in hot carrier graphene-based devices.

1. Introduction

Recently, researchers in the graphene community are much interested in hot carriers in graphene-based systems since they determine the performance of high-power and high-frequency electronics, thermal management of electronic devices, optoelectronic devices, the quantum Hall metrology, and bolometric detectors [16]. Most of the hot carrier graphene devices in the high carrier density limit (the Bloch-Gruneisen temperature () > the lattice temperature ()) show the dominant cooling power from a weak coupling of carriers to acoustic phonon processes [710], which is represented by a heat flow power law equation P = [11], where δ = 4 is a characteristic exponent, is the carrier temperature, is the lattice temperature, and is the coupling constant. Furthermore, disordered-enhanced properties in hot carrier graphene devices revealed the supercollision cooling processes, where δ = 3 [1214]. Recently, hexagonal Boron Nitride (h-BN) bottom substrate can be a great heat drained material for disordered graphene so as to extremely reduce the carrier-phonon scattering via Wiedemann-Franz law heat diffusion, which is δ = 2 [1517]. However, the carriers of energy relaxation in such h-BN/disordered graphene systems with less substrate phonon interactions following the Wiedemann-Franz law heat diffusion are rarely studied. A key parameter for discussing the cooling process is the energy relaxation time (), the characteristic time when the thermal energy is lost by the carriers [1820].

To date, CVD-grown graphene appears to be a good candidate for large-scale graphene-based applications. However, as such a system may not be air stable, it is highly desirable to cap the CVD graphene with an inert layer so as to experimentally realize stable, large-scale, and scalable devices. To this end, we have prepared CVD graphene with a top h-BN layer and studied the hot carrier effects in such a device. In this report, we studied the conductance fluctuations and hot carrier effects using a large measurement current (I) so as to raise the effective above on h-BN/CVD graphene and multilayer exfoliated graphene due to their disordered properties. Based on the self-thermometer linear properties between and due to Wiedemann-Franz law heat diffusion, we can investigate energy relaxation cooling processes under high carrier density limit () between the h-BN/CVD graphene and multilayer exfoliated graphene for future graphene-based applications. By determining from and input power , we can find that of the CVD graphene with a top h-BN layer is almost two orders of magnitude longer than those in exfoliated pristine monolayer/bilayer graphene. Such a result may find applications in hot carrier graphene-based transistors as a result of the weak electron-phonon coupling.

2. Experimental Section

2.1. Preparation of the Samples

As shown in the schematic diagram in Figure 1(a), we used the scotch tape method to mechanically exfoliate high purity and homogenous h-BN crystals synthesized by high pressure techniques [21] and transferred them by Gel-Pak polymer using the viscoelastic effect [22] on commercial chemical vapor deposition (CVD) graphene/285 nm SiO2/Si substrate [23]. The CVD graphene region outside h-BN sheet was etched by oxygen plasma so as to confine to our CVD graphene region under the h-BN sheet. CF4 gas was used to etch h-BN sheet protected by photoresist for 8 terminals and Cr/Au metal depositions as shown in Figure 1(b).

2.2. Electrical Measurements

An ac driving current from lock-in amplifiers passed into the source and drain contacts through the graphene devices for Hall-bar measurements. The magnetoresistance was measured in a He3 cryostat equipped with a superconducting magnet.

3. Results and Discussion

Figure 2(a) shows the longitudinal resistivity as a function of magnetic field with fixed current I = 20 nA at different temperatures that are equivalent to . The conductance fluctuations are observed, and they decrease as increase from = 0.32 to 50 K, which are typical properties in disordered mesoscopic graphene [2427]. Figure 2(a) inset shows the corresponding Hall resistivity at = 0.32 K, which can calculate the carrier density = 3.5 from the Hall slope and Hall mobility μ = 2092 cm2/Vs from carrier density and zero field resistivity. Figure 2(b) shows the results of at various driving currents from I = 20 nA to 30000 nA at a fixed of 0.32 K, which reveals similar conductance fluctuation characteristics by current heating in the nonequilibrium regime due to the hot carrier effects in disordered two-dimensional systems [18, 2831].

By using a conductance fluctuations-based thermometer between (I) and (T), we are able to reveal the clear conductance variations for lattice temperature and current dependence between B = 0.5 and 3 T by subtracting a smooth background that avoids the zero field weak localization peak and high field Shubnikov-de Haas-like oscillations as shown in Figures 3(a) and 3(b) [13]. Hence, we determine the root mean square (RMS) conductance fluctuation (, in units of ) for every () and (I) data set as a common graph (see Supporting Information (available here)) by averaging conductance fluctuations over the range of so as to assign an effective to the driving current as shown in Figure 4(a) [13, 18]. Interestingly, shows a linear dependence on . According to the seminal work done by Baker and coworkers about the energy transfer between carriers and the lattice [31], the following relation can be found:where α = 2/(p + 2) and is the exponent for the inelastic scattering rate . Therefore, Figure 4(a) shows that (I) of our sample A follows (1) with and , which suggested little carrier-phonon scattering in two-dimensional material heterostructure systems [26]. Furthermore, one would be interested in whether the heat dissipation was transferred by another mechanism rather than carrier-phonon scattering. Under the low-temperature limit (), , where is the sound velocity and is the Fermi radius. For graphene systems (s = 2.1 × 104 ms−1), ~ 54n1/2 K 93.4 K with the carrier density in units of in sample A [13, 32]. The energy loss (), the average rate of energy loss per carrier, is usually expressed to the carrier and lattice temperatures as with a constant and a characteristic exponent β. In general graphene systems, suggests the supercollision cooling mechanism in disordered systems [13, 18, 33] and suggests the two-dimensional acoustic phonon cooling processes in clean systems [34, 35]. Particularly, , the heat diffusion described by the Wiedemann-Franz law, was found in graphene on bottom h-BN substrate systems, where bottom h-BN substrate acts as a thermal conduction layer, effectively reducing the electron-phonon coupling [1416]. Such interesting bottom h-BN substrate underneath graphene can change graphene heat transport mechanism from electron-phonon interactions () to heat diffusion (). One might be interested in what heat transport mechanism will be in encapsulated graphene between bottom insulating SiO2 and top h-BN sheet as our sample A rather than traditional graphene on h-BN substrate [1416]. Based on the work of Baker et al. [31], β = p + 2, which suggested β in sample A should be close to 2 since the linear relation between and is p 0 (1 = α = 2/(p + 2)). Apparently, our results in sample A are highly consistent with this speculation for as shown in Figure 4(b), where = ( and are the width and length of our sample A, resp.). Interestingly, the heat transfer mechanism in our sample A with top h-BN sheet on graphene/SiO2 was dominated by heat diffusion (), which is the same as traditional graphene/h-BN substrate devices [1416]. Such interesting results suggest that the top h-BN sheet coupled with graphene is strong heat transfer medium in comparison with other coupled materials, like SiO2 or air. Also, such structures in the top h-BN sheet on graphene not only protected graphene without air molecular adsorbing doping [236], but also screened graphene from electrostatic force [37, 38], a great advancement for graphene-based devices.

In order to further discuss the linear relation between and that indicates little carrier-phonon scattering in our sample A, we fabricate multilayer exfoliated graphene (sample B) that showed the same linear relation between and by zero field resistance as a thermometer and conductance fluctuations due to its disordered property as in our previous reports [2326]. Again, we are able to measure the conductance fluctuations for temperature and current dependence as shown in Figure 5(a) and inset [13]. Consequently, we found the linear relation between and by utilizing () and (I) as a thermometer for determining (see Supporting Information) [13] as shown in the inset of Figure 5(b), which suggested p 0 and β 2. Under the low-temperature limit (),  K 188.6 K with the carrier density in units of 1012 in sample B [13, 32]. Consistently, is proportional to () in sample B as shown in Figure 5(b), which belonged to Wiedemann-Franz law heat diffusion and corresponded to the transport situation for little carrier-phonon scattering for α 1 [7, 26]. Interestingly, the multilayer exfoliated graphene (sample B) has the same heat transfer mechanism (β 2) as the h-BN/graphene (sample A) systems [1416].

Although the heat transfer mechanism of multilayer exfoliated graphene is the same as h-BN/CVD graphene for heat diffusion (β 2), the energy relaxation time can describe this cooling process about the typical time on which energy is lost from carriers [39, 40]. Therefore, the energy relaxation time can be expressed by Figure 6 compares the energy relaxation time between the h-BN/CVD graphene (sample A) and the multilayer exfoliated graphene (sample B). Interestingly, we could obviously find that of h-BN/CVD graphene (sample A) is about two orders of magnitude faster than that of multilayer graphene (sample B). Both of these two samples are the same n-type of carrier density and far from the Dirac point, which cannot be ascribed to the fast energy relaxation of hot carriers near the Dirac point of graphene [18]. These results suggest that the carrier-phonon scattering is absent in the multilayer graphene and the h-BN/graphene. The extremely long energy relaxation times in both devices (at least two orders of magnitudes longer than those in pristine exfoliated graphene [27, 28] and graphene on SiC [30, 34]) can be advantageous for applications in graphene-based hot carrier transistors [41] since carriers can maintain their high kinetic energy (and hence the high effective temperature) with a relatively low driving current.

4. Conclusion

We have studied conductance fluctuations and hot carrier effects caused by current heating on h-BN/CVD graphene and multilayer graphene as a self-thermometer. It has been shown that () is linearly proportional to I () in both of the disordered graphene devices, suggesting the little electron-phonon scattering and heat diffusion due to Wiedemann-Franz law. The extremely long energy relaxation time may find applications in graphene-based hot carrier devices. Such a result may be useful for other material systems [42, 43].

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

As an International Research Fellow of the Japan Society for the Promotion of Science (JSPS), C. Chuang acknowledges a JSPS Postdoctoral Fellowship and the MOST and Institute of Atomic and Molecular Sciences, Academia Sinica, for a postdoctoral fellowship under the Contract no. MOST 106-2811-M-001-164. C.-T. Liang would like to acknowledge the Ministry of Science and Technology (MOST), Taiwan, for financial support (Grant no. MOST 105-2119-M-002-048-MY3). N. Aoki acknowledges funding from JSPS KAKENHI (Grant no. JP16H00899 on Innovative Areas “Science of Atomic Layers”).

Supplementary Materials

This includes the method of determining the effective electron temperature using amplitudes of the observed conductance fluctuations. (Supplementary Materials)