Introduction

Recently, the rapid development of wearable electronics and sensors has led to the concept of the “Internet of Things” (IoT).1,2 In this idea, a network of billions of smart devices that connect people, systems and other applications collect and share data, providing key information for artificial intelligence (AI) as well as realizing a smart society. However, most wearable devices are still powered by batteries that are subject to frequent recharging and replacement. Therefore, this powering problem is one of the serious limiting factors for a future IoT society. One possible solution for wearable devices without a battery is energy harvesting technologies, and among them, thermoelectric energy conversion based on flexible materials has attracted much attention over the years.1,2 To date, flexible large-scale thermoelectric conversion devices have been investigated in van der Waals materials, such as organic conducting polymers, carbon nanotubes (CNTs) and their composites,3,4,5,6,7,8,9,10,11 and the best power factors among them are 2.71 and 2.46 mW m−1 K−2 for p-type and n-type devices, respectively.7,11

Another promising van der Waals material is graphene. In micron-scale single-crystal graphene films, extremely high maximum power factors of 36.0 and 19.0 mW m−1 K−2 have already been reported in p-type and n-type devices, respectively, although micron-scale single crystals are too small to apply in wearable devices.12,13,14 Moreover, centimeter-scale polycrystalline graphene films have been produced by chemical vapor deposition (CVD), opening a pathway for the development of flexible large-scale thermoelectric conversion devices.15,16 Therefore, we can conclude that a CVD-grown large-area graphene film is one of the potential candidates for the power generator for IoT devices. However, to the best of our knowledge, there is no report on the thermoelectric properties of large-area graphene films on a flexible substrate.

Here, we investigated the thermoelectric properties of CVD-grown polycrystalline graphene films on plastic substrates. To perform p- and n-type doping, we fabricated electric double layer transistors (EDLTs) of graphene films.17,18 By combining thermoelectric measurements and the EDLT technique, we continuously controlled the Seebeck coefficient and electrical conductivity. We observed maximum power factors of 6.93 and 3.29 mW m−1 K−2 for p- and n-type carrier doping, respectively, which exceeds the best previously reported results.

Results

Thermoelectric properties of CVD-grown large-area graphene films

We used CVD-grown large-area graphene films on a polyethylene terephthalate (PET) substrate.16 The details of the synthetic procedure are provided in Experimental section. The average layer number and grain size of the resultant graphene films are estimated to be two (see Supplementary Information S1) and ~10 nm, respectively.19,20,21 Figure 1a shows a schematic illustration of device fabrication. The original centimeter-scale graphene film on the PET substrate was cut into smaller pieces. Then, to pattern the graphene film, we used an adhesive tape as a shadow mask and removed excess graphene film by an O2 plasma treatment. Afterwards, Au/Ni electrodes were thermally deposited onto the graphene samples, forming the four-point-probe configuration. Importantly, we introduced a long channel length of 0.4 mm and a channel width of 2 mm to introduce a clear temperature difference between the source and drain electrodes for reliable thermoelectric measurements. For the gate dielectric layer, we prepared an ion gel film, which is a mixture of an ionic liquid, N,N-diethyl-N-(2-methoxyethyl)-N-methylammonium bis-(trifluoromethylsulfonyl)-imide ([DEME][TFSI]), and an organic polymer, poly(vinylidene fluoride-co-hexafluoropropylene) (P(VDF-HFP)).

Fig. 1
figure 1

a Schematic illustration of a procedure for fabricating an EDLT sample using CVD-grown large-area graphene films on a PET substrate. An optical image of the graphene-based EDLT sample with the four-point-probe configuration is also indicated. The ion gel film is a mixture of an ionic liquid, [DEME][TFSI], and an organic polymer, P(VDF-HFP). The scale bar corresponds to 1 cm. b The drain current (ID) – voltage (V4pp) curves at various gate voltages. A drain voltage ranging from −0.1 to 0.1 V was applied. c The four-point-probe electrical conductivity (σ) as a function of VGVCNP

Figure 1b shows the current-voltage curves using the four-point-probe configuration in the graphene-based EDLT on the PET substrate. We used a semiconductor parameter analyzer (Agilent Technologies, Inc. B1505A) and nanovoltmeter (Keysight, 34420 A) to collect the data, and all measurements were performed in an N2-filled glove box at room temperature. We applied drain voltages ranging from −0.1 to 0.1 V and simultaneously monitored both the drain current (ID) and the voltage drop (V4pp) at each gate voltage. During the measurement, the gate leakage current was less than 5 nA, indicating that the electrochemical reduction and/or oxidization of ions at the interface of the graphene film was negligible. The four-point-probe electrical conductivity was estimated by linear fitting of the IDV4pp curves and thickness of the graphene bilayer (0.7 nm), which was approximated from the Raman spectra results (see Supplementary Information S1). The average standard error of the four-point-probe electrical conductivity is 571 S cm−1, which is smaller than 6% of each fitting result. Additionally, Fig. 1c shows the gate voltage dependence of the four-point-probe electrical conductivity. The minimum electrical conductivity indicates the charge neutrality point (CNP), in which the hole carrier density and electron carrier density are equal. Therefore, in Fig. 1c, we represent the applied gate voltage, VG, in the relative gate voltage from VCNP (i.e., VGVCNP). As shown in Fig. 1c, ambipolar characteristics are clearly observed, suggesting the carrier doping of holes and electrons into the graphene films by the formation of electric double layers (EDLs).

Along with the measurements of the IDV4pp curves, we also investigated the VG dependence of the Seebeck coefficient (S) in the graphene-based EDLT using the method established for transition metal dichalcogenide (TMDC) monolayers.22 A schematic illustration of the thermoelectric measurement system is shown in Fig. 2a.22,23,24,25 The measurement details are provided in the Experimental section. As shown in Fig. 2a, two thermocouples were directly placed onto the surface of the graphene-based EDLT sample, and the temperature difference (ΔT = TA − TB) was induced by two Peltier elements. Both the ΔT and generated thermoelectromotive force (ΔV) were simultaneously measured. The S was calculated as the slope of the ΔV − ΔT profiles based on the definition of S = ΔVT. At a fixed gate voltage, we continuously measured the four-point-probe electrical conductivity and the Seebeck coefficient. Next, we repeated these measurements at various gate voltages.

Fig. 2
figure 2

a A schematic of the experimental setup for the thermopower measurements. The right-side temperature (TA) and the left-side temperature (TB) are measured by K-type (alumel–chromel) thermocouples. b The ΔV – ΔT profiles measured at various gate voltages. The red and blue colors indicate the hole and electron carriers, respectively. c The Seebeck coefficient in the EDLT of large-area graphene on the flexible substrate as a function of VGVCNP. d The power factor (S2σ) in the EDLT of large-area graphene on the flexible substrate as a function of VGVCNP

Figure 2b represents the ΔV − ΔT profiles, in which we obtained linear ΔV − ΔT responses for all measurements. In particular, the average standard error in all linear fitting results is 0.28 µV K–1, indicating that the derived S from our results is reliable. Because the sign of S depends on the type of majority carriers, both positive and negative S values should be observed in the ambipolar EDLTs for p- and n-type transistor operations, respectively.22,24,25 Figure 2c shows the VGVCNP dependences of S in the EDLT composed of large-area graphene on a flexible substrate. As shown in Fig. 2c, we achieved continuous control of S in the graphene-based EDLT on the PET substrate, and the shift in the polarity in S was reasonably observed. Interestingly, we observed the maximum S value in the hole-doping (57.1 μV K–1) and electron-doping (−36.6 μV K–1) regions. Notably, the difference in the maximum S between hole and electron doping has been previously reported, and this is most likely due to the effect of charged impurities on the substrates.12,13,14 Importantly, the S – σ relationship indicated no hysteresis behavior, meaning that we succeeded in the measurement of intrinsic thermoelectric properties in these CVD-grown large-area graphene films (see Supplementary Information S2). Because the gate voltage dependences of the electrical conductivity (σ) and the Seebeck coefficient (S) were collected, we evaluated the gate voltage dependence of the power factor (S2σ). As shown in Fig. 2d, similar to the results of S, we obtained the peak power factor in the hole-doping (6.93 ± 0.29 mW m–1 K–2) and electron-doping (3.29 ± 0.10 mW m–1 K–2) regions.

Discussion

To demonstrate the performance of our samples, in Tables 1 and 2, we summarize the previous reports on flexible large-scale van der Waals materials, organic conducting polymers, CVD-grown polycrystalline graphene films, and single-crystal graphene films.3,4,5,6,7,8,9,10,11,12,13,14,26,27,28 As shown in Fig. 3a, b, the peak power factors of our samples are higher than the best previously reported values among flexible large-scale van der Waals materials.7,11 Therefore, we can conclude that the CVD-grown large-area polycrystalline graphene film is a strong candidate for the power generator of IoT devices.

Table 1 Comparison of the maximum hole power factors (S2σ) of van der Waals materials. The results of the electrical conductivity (σ) and the Seebeck coefficient (S) at the maximum S2σ are also listed
Table 2 Comparison of the maximum electron power factors (S2σ) of van der Waals materials. The results of the electrical conductivity (σ) and the Seebeck coefficient (S) at the maximum S2σ are also listed
Fig. 3
figure 3

Summary of the power factors (S2σ) in van der Waals materials. Comparison of the a holes and b electrons in flexible large-scale van der Waals materials

We also compare our results with other CVD-grown polycrystalline graphene films.26,27,28 Although there has been no report on the thermoelectric properties of millimeter-scale polycrystalline graphene films prior to this study, the power factors of several hundred micron-scale graphene films have already been investigated. However, it should be noted that these previous studies did not perform thermoelectric measurements using the transistor configuration. Therefore, the carrier densities of previous studies are not optimized, and a fair comparison of maximum power factors among the other polycrystalline graphene samples is impossible. Interestingly, as summarized in Table 1 and Fig. 3a, our result is approximately one order of magnitude larger than that of previous reports. This difference might be understood by the optimized carrier doping in our samples. As shown in Fig. 3, we have succeeded in tuning the thermoelectric properties of CVD-grown polycrystalline graphene films via an EDLT.

Because there are several reports on optimized carrier doping in single-crystal graphene films, we also compare our results with reports on single crystals to understand the reason for the high performance of our samples.12,13,14 As shown in Tables 1 and 2, the Seebeck coefficient values at the maximum power factors in our polycrystalline films are slightly lower than those of single-crystal graphene samples. In contrast, very surprisingly, we obtained single-crystal-like high electrical conductivity at the maximum power factors in our samples. Consequently, the maximum power factors in our CVD-grown polycrystalline graphene were close to the values in single-crystal graphene. According to the above arguments, interestingly, we suggest that S and σ are not seriously affected by the grain-boundary effect in our samples. Therefore, the extremely high electrical conductivity in our films is expected to be an important point for understanding their high power factors.

To clarify the mechanism of the high electrical conductivity, we measured the temperature dependences of the electrical conductivity in the graphene-based EDLT sample. Figure 4a shows the temperature dependences of the four-point-probe electrical conductivity at gate voltages of −1.5, 0.3, and 1.5 V. The temperature range is set below 220 K, which is lower than the freezing temperature of the ionic gel film. As shown in Fig. 4a, when the gate voltage was close to VCNP (VG = 0.3 V), the electrical conductivity decreased with decreasing temperature. This insulating behavior suggests that the transport property near the CNP is possibly determined by the trap and/or mid-gap states induced by the polycrystalline nature of our films. On the other hand, as shown in Fig. 4a, metallic behavior was observed in both the p- (VG = –1.5 V) and n-doping (VG = 1.5 V) regions, indicating that the charge transport mechanism changed from insulating to metallic by carrier doping. To gain more details about the change in the transport property, Fig. 4b summarizes the electrical conductivities at various temperatures as a function of VGVCNP. Except for the result of VGVCNP ~0 V, the electrical conductivity increased with decreasing temperature at all gate voltages. Furthermore, as shown in Fig. 4c, the Sσ relationships for both the hole and electron carriers follow the trend of S σ−1 in the high conductivity region. This S σ −1 relationship is explained by the Mott equation, which describes the thermoelectric properties in metallic states (see Supplementary Information S3).29 Therefore, we can conclude that the excellent electrical conductivity and high thermoelectric performances in CVD-grown large-area graphene films on the PET substrate originate from the filling of trap states by the EDLT technique.

Fig. 4
figure 4

a The temperature dependences of electrical conductivity at gate voltages of –1.5, 0.3, and 1.5 V. The temperature range is set below 220 K, which is lower than the freezing temperature of the ionic gel film. b The electrical conductivities at various temperatures as a function of VGVCNP. c The |S| – σ relationships in hole and electron carriers. The trend of S σ –1, which is explained by the Mott equation, is also indicated

Finally, we should comment on future challenges to realize flexible thermoelectric batteries based on two-dimensional (2D) materials. In general, the conversion efficiency of thermoelectric materials is evaluated by the dimensionless figure of merit, zT = S2σT/κ, where T and κ are the temperature and the thermal conductivity, respectively. Based on our experimental results and previous work,30 the zT of our films is estimated to be 0.003 at 300 K. On the other hand, 2D TMDCs are expected to be promising materials because of their high Seebeck coefficient and low thermal conductivity.31,32,33 By combining the experimental and theoretical values, a zT value of 0.5 is predicted in TMDC materials.31 However, we would like to suggest that these trends are not simply applicable to 2D materials. Since the atomically thin thermoelectric materials are prepared on thicker thermally insulating flexible substrates, the induced temperature difference is mainly determined by the substrate. Although the thermal conductivity of graphene is extremely high, in this work, a large ∆T was induced in the graphene films by the PET substrate. Therefore, we can expect that the power factor of the active materials and the thermal transport properties including the substrate are more crucial for future 2D thermoelectric generators.

Furthermore, we calculated the output power of flexible thermoelectric generators based on the optimized graphene films (see Supplementary Information S4). Under the simplified conditions of no contact resistance, a temperature difference of 5 K, and an active area of 50–100 cm2, the maximum output power reaches ~1–10 µW, which is enough to operate wearable IoT devices, sensors, and processors. We should also compare our results with other energy harvesting technologies, such as photovoltaic, piezoelectric, and triboelectric techniques. Although these technologies can generate power greater than 10 µW, it is very difficult for them to operate continuously for a long time. For instance, photovoltaic batteries cannot generate electric power at night and on rainy days, and piezoelectric and triboelectric batteries require frequent mechanical motion for power generation. Therefore, thermoelectric energy harvesting is more suitable for wearable IoT devices.

In summary, we investigated the thermoelectric properties of CVD-grown large-area graphene films on a PET substrate by electric double layer gating. We achieved ambipolar transfer characteristics and realized the continuous modulation of positive and negative Seebeck coefficients. Moreover, the power factors in both the hole- and electron-doped regions indicated a dome-like structure by tuning the Fermi level of the polycrystalline large-area graphene films. Finally, their maximum power factors exceed the best recorded values among flexible large-scale van der Waals materials. These results strongly suggest that the CVD-grown large-area polycrystalline graphene film is a promising candidate for thermoelectric power generators.

Methods

Sample preparation

Polycrystalline Cu foil samples (23 cm × 20 cm × 33 μm thick) were treated with Ar/H2 plasma at 5 Pa for 20 min to clean their surfaces before the CVD process. After cleaning, CVD was performed (CVD conditions: 3 Pa, CH4/Ar/H2 = 30/20/10 standard cubic centimeters per minute, 4.5 kW for each microwave generator). The synthesized graphene films were transferred onto a PET substrate (188 μm thick). The device fabrication process was as follows. First, we cut the original CVD-grown large-area graphene film on the PET substrate into smaller pieces (1 cm × 1 cm) and patterned them by oxygen plasma treatment for 300 s using adhesive tape as a shadow mask. Second, Au/Ni (30 nm/5 nm) electrodes were thermally deposited onto the graphene samples. Finally, for the gate dielectric, we selected an ion gel, which is a mixture of an ionic liquid ([DEME][TFSI]) and organic polymer (P(VDF-HFP)). The gate dielectric layer was prepared by drop-casting onto the graphene samples.

Experimental details of the measurements of thermoelectric properties

The transfer characteristics of the graphene EDLT samples were evaluated using a semiconductor parameter analyzer (Agilent, B1505A) in an N2-filled glove box at room temperature. We applied a constant drain voltage of –0.1 V and gate voltages ranging from 1.5 V to −1.5 V. The four-point-probe electrical conductivity was estimated from the current-voltage curves by applying a drain voltage from –0.1 to 0.1 V. The voltage drop between the two voltage probes in the four-point-probe configuration was monitored using a nanovoltmeter (Keysight 34420 A). Along with the measurements of the transistor characteristics, thermoelectric measurements were performed. The EDLT sample was fixed on two Peltier elements (Laird Technologies CP1.0-31-05L), which were wired in series with a DC power supply (Keysight B2902A). To measure the temperature difference between the source and drain electrodes, we placed two K-type (alumel–chromel) thermocouples (Sakaguchi EH VOC Corp. TCKT0051) on the surface of the sample. The measured temperature was monitored by a temperature controller (Lakeshore model 335). In the measurement of the Seebeck coefficient, we turned on two Peltier elements to cause a temperature difference, and the generated thermoelectromotive force was simultaneously measured using a nanovoltmeter (Keysight 34420A).

Measurements of charge transport properties using the graphene EDLT on the PET substrate

Low-temperature experiments were performed in a physical property measurement system (PPMS, Quantum Design, model 6000). The gate voltage was applied at 300 K by applying the –10 mV drain voltage. The transistor characteristics were monitored by a semiconductor parameter analyzer (Agilent Technologies, Inc. B1500A). Simultaneously, the voltage drop between the two voltage probes in the four-point-probe configuration was measured using a nanovoltmeter (Hewlett-Packard 34420 A). We monitored the drain current in the temperature range from 300 to 100 K at each gate voltage.