Towards sensitive terahertz detection via thermoelectric manipulation using graphene transistors

Graphene has been highly sought after as a potential candidate for hot-electron terahertz (THz) detection benefiting from its strong photon absorption, fast carrier relaxation, and weak electron-phonon coupling. Nevertheless, to date, graphene-based thermoelectric THz photodetection is hindered by low responsivity owing to relatively low photoelectric efficiency. In this work, we provide a straightforward strategy for enhanced THz detection based on antenna-coupled CVD graphene transistors with the introduction of symmetric paired fingers. This design enables switchable photodetection modes by controlling the interaction between the THz field and free hot carriers in the graphene-channel through different contacting configurations. Hence a novel “bias-field effect” can be activated, which leads to a drastic enhancement in THz detection ability with maximum responsivity of up to 280 V/W at 0.12 THz relative to the antenna area and a Johnson-noise limited minimum noise-equivalent power (NEP) of 100 pW/Hz0.5 at room temperature. The mechanism responsible for the enhancement in the photoelectric gain is attributed to thermophotovoltaic instead of plasma self-mixing effects. Our results offer a promising alternative route toward scalable, wafer-level production of high-performance graphene detectors. Controlling how finger-like electrodes are positioned on graphene transistors makes it simpler to detect terahertz radiation for applications including medical diagnostics. When photons strike graphene sheets, multiple electrons with high kinetic energy are released. Changlong Liu at the Chinese Academy of Sciences, Shanghai, and coworkers have exploited these electrons to improve detection of terahertz light, which is useful for seeing features underneath human skin. The team’s device places single-layer graphene into a channel flanked by two microscale antennas, and then adds a pair of thin metal fingers on top of the 2D carbon film. By manipulating the separation of the finger contacts, the team unlocked an enhanced response to terahertz light at room temperature. Experiments and simulations revealed that temperature gradients between contact points generated a terahertz-sensitive photocurrent through conversion of heat into electricity. Graphene is very promising for THz application, especially in the fields desiring fast THz imaging, like security screening, communication, biomedical and pharmacy control. However, current graphene-based THz detectors are severely hindered by its poor switching behavior, and lack of internal gain to boost up the responsivity in the framework of plasma-wave self-mixing or thermoelectric detections. A proper means to selectively trigger the photoelectric conversion is highly desirable for practical applications. This work offers alternative way for room temperature THz detection via manipulating the hot carriers either electrically and electromagnetically. Internal photoelectric gain is achieved via the bias-field effect, which leads to the giant enhancement of graphene-based detector’s responsivity over 200 V/W. Furthermore, switching behavior between photoconductive and photovoltaic modes can be given rise in a single device, being promising for scalable THz imaging.


Introduction
Recently, two-dimensional (2D) materials (transition metal dichalcogenides (TMDs), graphene, black phosphorus, etc.) have ignited intensive interest due to their novel electronic and photonic properties, which are distinct from those of their bulk counterparts, such as defectfree topological transport 1-3 , linear dichorism 4,5 , superconductivity 6,7 , and plasmonics 8,9 . It has been proved experimentally that few-layer 2D materials can strongly interact with light 10,11 , which indicates they are promising materials with great potential in optoelectronic devices for communication, sensing, and imaging, particularly in the visible and near-infrared range. However, it is rather difficult to extend the spectral response of 2D materials to the mid-/far-infrared and terahertz (THz) regime owing to the relatively large bandgap of these materials. Many efforts, such as bandgap engineering and heterogeneous integration with quantum dots or other 2D materials (van der Waals heterostructures), have been dedicated to resolving this bottleneck [12][13][14] . Nevertheless, light absorption in narrow-gap 2D materials such as TMDs decreases significantly when the photon energy is far below the interband threshold, similar to the behavior of other bulk semiconductors.
On the other hand, gapless graphene, in principal, can interact strongly with electromagnetic fields over a much wider spectral range. The photoresponsivity of graphene is usually quite low due to the weak light absorption of 2.3% through the interband transition as well as a large dark current 15 . Nevertheless, by utilizing the Drude response of free electrons in combination with heavy chemical doping, one can achieve 40% photon absorption in a single graphene layer in the mid-/far-infrared band of the electromagnetic spectrum [16][17][18][19] . Therefore, graphene is one of the most preferred materials for THz optoelectronics. Graphene-based THz detectors have been reported to exhibit suitable performance in the framework of plasma self-mixing and thermoelectric effects. The plasma self-mixing effect can be achieved via the asymmetrical antenna coupling of a THz field between a source/drain and gate, during which the nonlinear rectification of carrier transport leads to direct detection with a responsivity of approximately 0.04~20 V/W [19][20][21][22] . Moreover, thermoelectric THz detection in graphene can be achieved by asymmetric contacting, which produces a Seebeck coefficient difference and direct photocurrent along the channel due to metal-induced asymmetrical doping 18 , and the responsivity can reach as high as 15 V/ W. While high-performance graphene-based devices have become increasingly critical to satisfying the requirements of everyday applications, the current implementation of direct THz detection is hindered by the lack of sufficient photoelectric gain 23 . The superior performance of graphene can be achieved by engineering the bias-field and electromagnetic size effects to manipulate the hot-carrier distribution. The results of the current study demonstrate the important role of carrier-THz field interactions in detection, which enable switching between the thermophotovoltaic (PV) and photoconductive (PC) modes, thereby producing previously unattainable levels of photoelectric gain.

Device fabrication
The detector was fabricated using chemical vapor deposition (CVD)-grown single-layer graphene, which was transferred onto a highly resistive Si/SiO 2 (300 nm) substrate for antenna patterning. The purpose of using a highly resistive Si substrate instead of a highly doped one was to eliminate the reflective loss of the incident wave, as the THz wavelength is at least two orders of magnitude longer than the dielectric layer thickness. The graphenechannel was contacted symmetrically by a log-periodic antenna to maximize the field confinement along the channel rather than perpendicular to it (see Fig. 1 and 2 (a), (b)). Thus, particular asymmetrical boundary conditions related to previously reported photothermoelectric (PTE) 18,19 and plasma-wave effects 21,22 were not considered. Oxygen-based plasma etching was performed after the antenna patterning and lift-off processes to form a graphene-channel (for more details, see Fig. S1, Supplementary Information). Additional paired contact fingers with different slit widths s were then fabricated in the middle of channel by electron beam lithography, resulting in a four-terminal configuration 24,25 . The terminals are labeled A, B, C, and D in Fig. 1b. Thus, the device response depended specifically on the different electrical connections among these terminals. The homogeneity of the formed graphene-channel was confirmed by Raman scattering spectroscopy, and the current-voltage output characteristics with terminals A and D connected were determined for all the fabricated devices. It could also be inferred from the blueshift of the optical phonon-related G peak in the Raman spectrum that the Fermi level relative to the Dirac point of the transferred sample is approximately 0.3 eV (see Fig. S7, Supplementary  Information).

Simulations
Because the initial motivation of our work was to employ the electromagnetic size effect to enhance THz absorption, we performed a full-vector Finite Difference Time Domain simulation to better understand the electromagnetic size effect of the fabricated devices before further performance characterization. The geometrical parameters of the model strictly follow the structure of the fabricated devices. The metal layer was considered a perfect conductor, with perfectly matched layer (PML) absorbing boundary conditions around the device to remain consistent with the experimental conditions. show the corresponding amplitudes of the THz field at 0.12 THz for the same device without and with the paired contacts. The comparison of these two configurations shows that the amplitude of the THz field for the latter configuration is much larger than that for the former one, which indicates that the additional paired contacts can significantly improve the electromagnetic coupling. The spectral characteristics of the devices without and with paired fingers are depicted in Fig. 2(c), (d), respectively, from which it can be determined that the localized near field in the center of the antenna is more than doubled with the incident frequencies being tuned Simulated near-field distribution at 0.12 THz in the center of devices without and with paired fingers, respectively. c, d Respective spectral profiles of the designed devices without and with paired fingers. e Comparison of photocurrent response between the fabricated devices without and with contact fingers measured at 0.12 THz and a power density of 0.6 mW/cm 2 from 0.04 to 0.2 THz. The largest field amplitude is observed at approximately 0.04 THz (λ~8 mm) when the diameter of the antenna (D~2 mm) is consistent with the wavelength λ/n sub of the electromagnetic wave in the substrate (the refractive index n sub is approximately 3.5 for highly resistive silicon). On the other hand, the incorporation of the paired fingers does not clearly alter the spectral profile of the antenna, and a broad photoresponse from 0.08 to 0.15 THz is well preserved. The above-mentioned simulated results were confirmed experimentally with the fabricated device, as shown in Fig. 2e, which shows that the photocurrent responses at 0.12 THz for the devices with and without paired contacts differ by a factor of four. Thus, our strategy of using paired contacts to enhance performance is valid.

Results
We characterized two types of devices with different resistances and structures. The I-V characteristics of the fabricated devices were measured beforehand using an Agilent B2912A (see more details in Figs. S2 and S4, Supplementary Information). The electromagnetic response properties of the fabricated devices were characterized using a lock-in technique with a reference signal generated from an Agilent E8257D microwave source, which was subsequently connected to a VDI tripler to cover photon frequencies from 0.04 to 0.12 THz. The THz wave was outer-coupled with a horn antenna at a power density of 0.6 mW/cm 2 . To study the detection mechanism and response speed of the devices, 0.12 THz radiation was used under fast on/off modulations. Figure 3a, d show 3D schematics of the wide-(W~60 μm, with R AD~1 00 Ω) and narrow-(W~6 μm, with R AD 1 kΩ) channel devices, respectively. Figure 3b, e show the measured photocurrent I ph-AD vs. the electrical bias U AD for the corresponding devices. In addition, for the narrow-channel device shown in Fig. 3d, the slit width s between the two fingers was varied from 0.5 to 1.5 μm to further verify the electromagnetic size effect. The resistance R AD remained nearly constant even as the slit was varied (see Fig. S8a, Supplementary Information). The photocurrent response in Fig. 3b is only approximately two times higher than that in Fig. 3e, even though the resistance of the wide-channel device was an order of magnitude smaller than that of the narrow-channel one. Therefore, the voltage responsivity of the narrow-channel device should have been much higher, due to improved THz photon funneling in the narrow-channel device (more circular teeth can be achieved in a narrow-channel device, which is beneficial for THz field focusing; see the simulated results in Fig. S11a, b, Supplementary Information). Figure 3b, e show that non-zero crossings consistently occurred in the studied devices, particularly for those with a narrower channel (Fig. 3e). Such behavior is attributed to the non-equilibrium hot-carrier diffusion triggered by the asymmetrical THz near-field distribution along the channel due to the dislocation of paired fingers during fabrication (see Figs. 3g and Fig. S11c, Supplementary Information).
Formally, the thermoelectric effect can be given by 26 where S ¼ À ficient retrieved from the Mott relation, which depends on the doping level of graphene 27,28 , and dT h /dx is the temperature gradient along the channel depending on local photon absorption. Equation 1 indicates that the thermoelectric current can be generated either by a temperature gradient within a non-zero Seebeck coefficient or by the local Seebeck coefficient difference. The doping distribution of the graphene-channel in our devices was homogeneous, with a Fermi-level difference less than 0.05 eV (see Figs. S2 and S8 in Supplementary Information) 27 ; therefore, the difference was negligible for producing a sizable Seebeck coefficient difference. In this respect, thermoelectric photocurrent arises from the temperature gradient induced by the asymmetrical THz field coupling along the channel, rather than from differences in the Seebeck coefficient 28 . Because of carrier (hole) diffusion from the hot to the cool side, as shown in Fig. 3h, a net current flows along the channel in the absence of the electrical bias U AD . The temperature rise of carriers after local photon absorption can be estimated following the zero-bias photocurrent I ph-AD based on the relationship ΔT h = ΔU ph-AD /S = I ph-AD R AD /S, where ΔU ph-AD is the open-circuit photovoltage. Therefore, the temperature rise ΔT h for hot carriers should be <10 and 60 K for the wide-and narrow-channel devices, respectively.
Furthermore, the photocurrent I ph-AD grows linearly by more than an order of magnitude under a finite electrical bias U AD . Because the PTE effect induces direct photocurrent generation and does not change with the static electric field, the increase in the photocurrent with the bias voltage can be attributed to the onset of additional resistance changes in the channel under THz radiation. Two competing potential photoresponse mechanisms require an electrical bias: (i) the bolometric effect, which is caused by changes in carrier mobility due to lattice heating, and (ii) the photoconductive effect, which is caused by changes in carrier number. Thus, the total dc current can be formulated as I AD (+) = I AD-off (+) -|I B | + I PE + I pc (+) . The last term I pc~e μn*U AD is the photoconductive contribution due to the change in carrier density n* and |I B | is the bolometric current summed with the opposite sign due to the metallic behavior of graphene under lattice heating 29 .
To determine the dominant mechanism of linear photocurrent growth, the absorbed power must be accurately determined first 30 . The absorption rate of graphene in this work is mainly attributed to the Drude response and can be expressed analytically as where n sub~3 is the Si index and Z 0 = 377 Ω is the impedance of free space), following the methodology of Ref. 31 . Because the incident power at 0.12 THz onto the active channel is approximately 0.4 μW, based on the electromagnetic simulation (see Fig. S10, Supplementary Information), the absorbed power P abs can be determined to be approximately 0.05 μW. The thermal resistance of graphene is proved to be on the order of (T/5) −3 K/nW4 .6 × 10 3 K/W (T~300 K is room temperature) 32 , and thus, the temperature rise ΔT L is on the order of 10 −4 and 10 −5 K, corresponding to the incident and absorbed power, respectively. The temperature coefficient of resistance in our device was approximately 4 × 10 −3 %/K; thus, the resistance change should be <0.004 and 0.04 Ω for the wide-and narrow-channel devices, even if we suppose the temperature rise ΔT L is approximately 1 K. On the other hand, based on the experimental results shown in Fig. 3b, e, the resistance change can be determined to be approximately 0.1 and 5 Ω for the wide-and narrowchannel devices respectively, based on the relationship ΔR AD~Iph-AD R AD 2 /U AD . Such a large discrepancy between the experimental and calculated results implies that the bolometric effect can be ruled out, and the resistance change under THz radiation is attributed to the change in carrier number. Although the photocarriers cannot be excited directly via interband transition, they can be generated by the charge transfer driven by the thermoelectric potential, as shown Fig. 3i. Under electrical bias U AD , the thermoelectric photocurrent can be neglected, and photoconduction becomes the dominant effect when U AD > ΔU ph-AD (ΔU ph-AD is the open-circuit photovoltage), as shown in Fig. 3i.
To characterize the response speed, the devices were subjected to THz radiation with fast on/off modulation, as shown in Fig. 3c, f. The pulsed shape was well preserved with a good signal-to-noise ratio for both the wide-and narrow-channel devices. The response time, defined as the timescale between 10 and 90% on the rising edge or, analogously, on the falling edge, was approximately 4 μs (see Fig. 3f and Fig. S6 in Supplementary Information). This response speed is already adequate for video-rate THz imaging applications and is actually limited by the transmission loss of current readout electronic systems.
The results in Fig. 4 are based on the electrical contacting configuration shown in Fig. 4a. The following analysis and discussion focus on the devices with a narrow-channel due to their higher performance. The output photocurrent I ph-AD was measured with application of a finite bias voltage U AB between contacts A and B. Figure 4b shows that I AD can be pinched off and its sign reversed when U AB is sufficiently large. During this process, I AD shows a nonlinear relationship with the bias voltage U AB . Owing to an electrical connection similar to that of a triode, the current will change significantly depending on the sign of the bias voltage. Therefore, we referred to the process as a "bias-field effect", as no strong PN junction effects such as those in a triode can arise in the device due to the absence of a bandgap. However, such nonlinear behavior is favorable to the rectification of THz signals (see, for example, the self-mixing in Fig. S3, Supplementary Information) depending on the incident frequency and metal-graphene interface. To clarify the THz detection mechanism with this configuration, here we consider a more general case with the detuned frequency far from the self-mixing regime. Figure 4d displays the measured signal vs. the U AB bias voltage, where the photocurrent I ph-AD experiences its minimum at U AB0 .5 V and then increases almost linearly as U AB continues to increase. During this process, no additional bias voltage U AD , such as that in Fig. 3, is needed, as shown in Fig. 4d, which indicates the existence of a direct photocurrent generation mechanism.
Generally, there are also two competing mechanisms for direct photocurrent generation: (i) the above-mentioned thermoelectric effect and (ii) the plasma-wave self-mixing effect, which is related to the rectification of the nonlinear IV characteristics. The photocurrent associated with the self-mixing effect is proportional to the first derivative of conductance 22 : To generate photocurrent, both competing mechanisms require asymmetrical boundary conditions. In previously reported graphene FETs 22 , the self-mixing effect can arise due to asymmetrical antenna coupling between the source and gate, leading to a longitudinal direct electric field and photocurrent/photovoltage generation. The photocurrent derived from the self-mixing effect reverses its sign following the derivative dσ AD /dU AB near the Dirac point. Figure 4c depicts the derivative of the conductance σ AD / dU AB as a function of U AB . The derivative changes sign found at U AB~0 .5 V in Fig. 4c. However, the linear shape of the photocurrent curve in Fig. 4d departs significantly from that in Fig. 4c, and there are no hints of sign change in the photocurrent at any bias points. To this case, the plasma self-mixing effect should be ruled out. Further proof of the role played by the self-mixing effect in photocurrent generation can be found in the Supplementary Information, with paired contact fingers replaced by other types of metallic contacts (e.g., Ti/Au, see Fig. S2, Supplementary Information), where the critical point of U AB for photocurrent polarity inversion coincides with the sign change of the derivative dσ AD /dU AB . These results clearly demonstrate that the tuning conditions required by the self-mixing effect depend on the properties of metal-graphene contacts, such as parasitic capacitance and resistance. Furthermore, the results indicate an alternative route to achieving nonlinear rectification of THz waves via bias-field configuration. In the absence of self-mixing, the hot carriers still prevail due to the Drude response, and the thermoelectric response can generally be activated. Figure 3e illustrates the physics underlying the photocurrent response: (i) the asymmetrical boundary is satisfied when the bias voltage U AB is applied because of charge transfer at the metal-graphene interface (Fig. 3e) 33,34 . (ii) Diffusion of non-equilibrium hot carriers under THz excitation between contacts A and D cannot be neglected due to the non-symmetrical carrier distribution (Fig. 3f, i.e., ∫ AD S(x)dT h ≠0), which would lead to direct photocurrent generation. As the bias voltage U AB increases, the Fermi level of graphene continues to shift upward, and the Seebeck coefficient difference among the local channels also increases, facilitating thermoelectric photocurrent generation, as shown in Fig. 4d. Following this process, the photocurrent is improved by more than an order of magnitude when U AB is only 1 V. Figure 4d shows that the bias-field-induced photocurrent is at least eight times higher than that generated through the photoconductive effect (Fig. 3e). A similar effect can be observed under THz radiation at other frequencies, e.g., 0.3 THz (1 mW/cm 2 ), as shown in the inset of Fig. 4d, and at a lower frequency of 0.04 THz (see Figs. S2 and S5 in Supplementary Information), all of which exhibit similarly drastic enhancements in direct photocurrent under higher bias fields. In this study, the photocurrent of 30 nA observed at 0.3 THz was an order of magnitude weaker than that recorded at 0.12 THz, likely due to the decreased coupling efficiency of electromagnetic waves according to the simulated results shown in Fig. 2c, d. To evaluate device sensitivity, the voltage/current responsivity and noise-equivalent power (NEP) of the detector can be obtained via the following equation: where R AD is the resistance and P in~Pout S A /S b is the incident power on the device; here, S A (~3 mm 2 ) is the total area of the device, including the antenna area that overlaps with the plane-wave output in the beam spot area S b (S b~1 0 cm 2 ) and P out is the output power (6 mW at 0.12 THz). Figure 5a shows the U AD bias dependence of Fig. 4 The direct photocurrent/photovoltaic response I ph-AD under 0.12 THz radiation (power density of 0.6 mW/cm 2 ), tuned by the bias voltage U AB between contacts A and B. a Schematic of the bias-field configuration. The photocurrent I ph-AD is readout as a function of the bias voltage U AB . b U AB dependence of the dark current I AD ; the bias voltage U AD is approximately 100 mV. c Derivation of conductance σ AD vs. U AB : σ AD −1 dσ AD /dU AB . d Direct photocurrent I ph-AD vs. U AB under 0.12 THz radiation. The inset shows a similar trend of photocurrent response at 0.3 THz. e Band diagram of the "bias-field effect" induced by direct photocurrent generation. f The symmetry of carrier distribution is broken when the electrical bias U AB is applied due to charge transfer between the metal and graphene 33,34 , which ultimately facilitates non-equilibrium hot-carrier diffusion in our devices.
T(x) denotes the hot-carrier temperature distribution; S(x) is the Seebeck coefficient the photoconductive responsivity for the narrow-channel devices. Although much lower photocurrent occurred in the narrow-channel devices than in the wider-channel devices, the voltage responsivity was improved by at least a factor of five, reaching a value of 15 V/W, and higher responsivity could be anticipated at a higher bias voltage U AD. Because the antenna size was the same for both the wide-and narrow-channel devices, the higher responsivity can be attributed to the tailorable near field in close proximity to the channel. In this regard, it is clear that the effective area for funneling the electromagnetic wave into the channel can be much smaller than the actual antenna area used in this work (see Figs. S9 and S10 in Supplementary Information). Therefore, the responsivity can be improved greatly by exploring the electromagnetic size effect of a subwavelength antenna with a smaller gap between the antenna arms to improve the dipole oscillation across the channel. Based on the above-mentioned results, the maximum voltage responsivity of our device reached as high as 280 V/W, by employing the "bias-field effect" electrical configuration shown in Fig. 5d. Such remarkable enhancement in device performance is similar to that achieved by the plasma self-mixing effect, which produces a responsivity maximum near the threshold (see Supporting Information, Fig. S3) and can be understood as an interior loading effect (ΔU~R ch SΔT/R 1 , where R ch is the total resistance and R 1 is the hot spot channel resistance for thermoelectric generation) 35 . A similar voltage dependence could be observed at 0.3 THz, which showed the lowest efficiency due to the mismatched THz coupling with the current antenna structure, as shown in Figs. 2 and 4. Nevertheless, the responsivity observed at 0.3 THz was modest and could exceed 10 V/W, indicating performance comparable to that of other optimized detectors operating at a similar frequency 36 .
Another important parameter for representing the sensitivity of detectors is the NEP, which is extracted from the ratio V N /R v by assuming that the main contribution to the noise is the thermal Johnson-Nyquist noise (N th ) associated with the non-zero resistance of the FET channel (V N = ffiffiffiffiffiffiffiffiffiffiffiffi 4k B Tr p ) 37 , where k B is Boltzmann's constant and T is room temperature. Although this hypothesis neglects shot noise, it provides a lower limit for the NEP. Figure 5b, e show the NEP of our devices operating in the photoconductive and photovoltaic modes, respectively. In addition, the minimal NEP of our devices was 0.1nW/Hz 0.5 when operating in photovoltaic mode. For the photoconductive mode, the NEP was approximately five times larger (~0.5 nW/ Hz 0.5 ), which is also moderate compared with recently reported results 18,20 . Therefore, switching between different working modes can be achieved within a single device, depending on the configuration of the electrical contacts and the subsequent readout set-up. Figure 5c, f show the response spectra of the devices operating in different modes. The average responsivities of the devices operating in the photovoltaic (U AB~0 .5 V) and photoconductive mode (U AD~0 .6 V) were 100 and 10 V/W, respectively 38 .

Discussion
Based on the foregoing discussion, it can be concluded that the considerable improvement in the NEP is ultimately the result of the enhanced responsivity originating from careful antenna design, interface control, and THz field focusing. The generation of PTE current can be formally boosted by either asymmetrical electromagnetic heating or dissimilar electric contacts with significant Seebeck coefficient differences. The main limitation of the PTE effect is the lack of sufficient internal gain, as can be obtained in a normal interband photon detector in which multiple charge carriers can be generated with only one absorbed photon 35 . Such difficulties can now be circumvented in a photoconductive mode with weak electrical or structural asymmetry. The hot carriers diffused from the metal-graphene interface are efficiently extracted under DC electric field, leading to a large photoconductive gain due to both the high carrier mobility and Fermi velocity of graphene. Furthermore, the bias-field-dependent effect clearly modifies the current-voltage characteristics or hot-carrier distribution, enabling a drastic enhancement in thermoelectric efficiency for direct photocurrent generation.
As previously mentioned, the responsivity observed for the devices fabricated in this work are still rather conservative estimations owing to the over-sized antenna and rather long channel length employed in the current design. On one hand, the device response mainly relies on the power density focused onto the channel rather than the antenna area. In our case, the antenna was an order of magnitude larger than previously reported ones 16,[18][19][20] . Indeed, the antenna can retain a superior radiation funneling capability with a much smaller effective antenna area, even though it is smaller than the diffraction limit area defined by S~(λ/2n sub ) 2 (λis the wavelength and n sub 3.5 is the refractive index of the substrate; also see the Supporting Information for a detailed calculation of the effective antenna area). On the other hand, the response also depends on the near-field tailoring of the received radiation from the antenna to the active channel. Thus, the performance can be improved by further exploring the electromagnetic size effect by shortening the channel length or enhancing the dipole interaction between the two opposite antenna arms, which is exactly what we proposed. Therefore, it is reasonable to expect that a re-designed antenna with a scaled area in conjunction with a distinct size effect can substantially boost the interaction between graphene and THz photons, which is essential to obtaining both higher responsivity and higher chip-level pixel intensity for potential THz imaging devices.

Conclusions
In summary, we demonstrated enhanced roomtemperature THz detection using antenna-coupled graphene transistors via the introduction of paired finger contacts, which enable different detecting modes driven by different mechanisms, including photoconductive, thermophotovoltaic and plasma self-mixing effects. In particular, through the manipulation of the electrical connection, a novel "bias-field effect" could be activated, leading to a drastic enhancement in THz detection ability with a responsivity of up to 280 V/W and a Johnson-noise limited minimum noise-equivalent power (NEP) of 100 pW/Hz 0.5 at room temperature, indicating a performance level competitive with that of the best mechanically exfoliated graphene thermoelectric detectors. In addition to well-established strategies, our work shows that the metal-graphene interface, workfunction engineering, and localized electromagnetic engineering can play prominent roles in improving THz detection. The room-temperature operation, fast response, high sensitivity, and low NEP of the devices reported in this work offer a promising alternative route to scalable, wafer-level production of high-performance graphene detectors.