Search tips
Search criteria 


Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
Sci Rep. 2017; 7: 10625.
Published online 2017 September 6. doi:  10.1038/s41598-017-09809-7
PMCID: PMC5587735

Contactless graphene conductivity mapping on a wide range of substrates with terahertz time-domain reflection spectroscopy


We demonstrate how terahertz time-domain spectroscopy (THz-TDS) operating in reflection geometry can be used for quantitative conductivity mapping of large area chemical vapour deposited graphene films on sapphire, silicon dioxide/silicon and germanium. We validate the technique against measurements performed with previously established conventional transmission based THz-TDS and are able to resolve conductivity changes in response to induced back-gate voltages. Compared to the transmission geometry, measurement in reflection mode requires careful alignment and complex analysis, but circumvents the need of a terahertz transparent substrate, potentially enabling fast, contactless, in-line characterisation of graphene films on non-insulating substrates such as germanium.


The development of scalable integrated manufacturing pathways for graphene is crucial to all its emerging applications and industrial development1. Chemical vapour deposition (CVD) has become the dominating technique to synthesise large area “electronic-quality” graphene2, with the size of single mono-layer graphene crystals now on the cm-scale3 and continuous films now routinely produced roll-to-roll or at a size just limited by the reactor4. In fact, progress in growth reached a level where detailed, adequate characterisation over such large areas has become a key challenge. Prevailing electrical characterisation for instance is based on the fabrication of field-effect or Hall bar devices typically combined with Raman spectroscopy. This is time-consuming for large samples and for statistically relevant sample numbers in particular also as the as-grown graphene typically has to be transferred away from the growth substrate. Among a range of emerging contactless characterisation methods, terahertz time-domain spectroscopy (THz-TDS) operating in transmission geometry has been demonstrated to allow the direct, accurate mapping of graphene conductivity and mobility over large areas, producing data consistent with the Drude model to describe graphene intra-band transitions510. Graphene’s complex conductivity is determined by the terahertz pulse transmitted through the graphene film relative to the support, and analysed using Fresnel coefficients where graphene is modelled as an infinitely thin conducting film. Drift and field-effect mobilities can then be extracted by fitting the conductivity spectra to the Drude model9, and measuring graphene conductivity changes as a function of the applied back-gate voltages on a gate-stack support8, respectively. By repeating the measurement and analysis across the entire graphene area, a conductivity or mobility map can be reliably produced. While these demonstrations could potentially enable a rapid in-line graphene monitoring and large-area characterisation, to date THz-TDS has been carried out in transmission mode that necessitates a terahertz transparent support.

Here we overcome this restriction and demonstrate, as a proof-of-concept, the quantitative contactless measurement of the electrical conductivity of CVD graphene on a range of application relevant supports using THz-TDS operating in reflection geometry. Reflection based THz-TDS has previously been used to characterise optically dense materials where due to the high energy loss within the sample, transmission based measurement cannot be used1115. As schematically outlined in Fig. 1, we first validate our measurements in reflection geometry by performing the more well-established terahertz transmission conductivity mapping5, 7, 10 on the same graphene sample. For this we use CVD grown graphene synthesised on commercial Cu foils and transferred to a sapphire substrate, which has the required transparency in the THz frequencies. We then show that by sample back-gating, our proposed method can be used to resolve graphene conductivity changes and hence to directly determine the graphene mobility on a p-doped Si substrate with a 300 nm thick SiO2 layer, which is one of the most commonly used substrates for graphene device manufacturing. To illustrate that this technique has potential as a tool for in-line graphene quality monitoring, we show that the graphene conductivity can also be directly mapped on a substrate like Ge, which has a considerable number of intrinsic carriers at room temperature (42 Ω·cm) and exhibits a Drude like carrier absorption16. In particular we demonstrate THz-TDS mapping for graphene that has been directly grown on Ge.

Figure 1
Schematic of THz-TDS setup where pulsed terahertz radiation is directed at graphene on a range of substrates: (a) sapphire support to allow terahertz reflection and transmission measurements used for the initial validation of our method, (b) boron-doped ...


Graphene growth and transfer

As reference process and material, we use well-established graphene CVD on commercial Cu foil and subsequent PMMA transfer, which is widely used in literature17. We use standard Cu foils (25 μm thick, Alfa Aesar purity 99.8%) and CH4 as the carbon precursor18. For transfer, PMMA (poly methyl-methacrylate) was used as support, followed by FeCl3 chemical etching to remove the Cu. As target substrates we used sapphire (430 μm thick), Si/SiO2 wafer (300 nm/525 ± 25 μm thick), and intrinsic single crystalline Ge (110) wafer (50 Ω·cm, from MTI Corporation). Raman spectroscopy was performed using a 532 nm laser for characterising the transferred graphene. The Si wafer support was boron-doped (100 Ω·cm) to allow back-gating. The process conditions for graphene CVD directly on Ge were similar to recent reports in literature19 and we used Ge (110) wafer substrates and a 1:52 precursor gas mixture of CH4/H2 at a growth temperature of 920 °C in a Aixtron Black Magic cold wall CVD reactor.

Terahertz time-domain reflection spectroscopy

Reflection based THz-TDS experiments were conducted with a Terahertz Pulsed Imaging (TPI) Imaga 2000 system (TeraView, Cambridge, UK), as schematically shown in Fig. 1. The terahertz radiation used here is broadband, covering a spectral range of 0.15–3 THz in free-space. Terahertz radiation is generated by pumping a biased photoconductive antenna with an ultrashort laser pulse from a Ti:Sapphire laser. The emitted terahertz pulse is collected, collimated, and then focused onto the sample with a focal length of 7 mm at an incident angle of 30°. The reflected terahertz pulse is then collected and focused onto an unbiased photoconductive antenna for the laser-gated terahertz detection. The TPI achieves a spatial resolution of approximately 400 μm at 1 THz20 that in turn allows us to estimate the spot size of approximately 420 μm with Gaussian beam optics21. One of the main barriers for accurately extracting optical parameters in reflection geometry is the great sensitivity to any phase misalignment between the sample and reference measurements. In general, the phase misalignment can be mitigated by performing additional measurements with a slab of material of known optical constants14, maintaining same path length change for both terahertz and optical beam15 and numerical correct phase correction using the maximum entropy method22. Attention was given to precise sample positioning by having both the Al mirror and the sample mounted on a motorised stage and positioning the respective front surfaces in order to ensure that the particular reflecting plane reflects the incident wave into the detector at maximum level12. It should also be noted that the substrate refractive index measured is not adversely affected by phase misalignment, as in the case for the extinction coefficient13. As an experimental check, the measured substrate refractive index for s-polarisation23 was always compared against the literature values. The reflection coefficient depends on the polarisation of the incident terahertz wave and the angle of incidence. Given the substrate refractive index and using Tinkham’s formulae to describe the effect of a thin conducting film24, the equation for obtaining conductivity from s-polarisation reflection measurements can be derived as:


where r˜ is the Fourier transformed ratio of the reflected wave’s complex electric field from the sample to the incident wave from the mirror, z0 is the vacuum impedance (376.7 Ω), θ i and θ t are the incident and transmitted angles, respectively, and n˜1 and n˜2are the complex refractive indices of air and substrate, respectively25, 26. At normal incidence and assuming zero conductivity, we get the well-known Fresnel coefficient for the reflection at the interface:


The analysis becomes simpler at normal incidence but such a reflection system in turn would require a beamsplitter that is not commonly implemented due to reflection losses27.

Using Snell’s law cosθt=1sin2θin˜22, the complex conductivity can be expressed as


We performed analytical simulations to double check the derived expressions and optical constants, whereby the simulations involved (1) simulating a terahertz pulse generated from a photoconductive antenna switch with realistic parameter settings such as 120 fs, 300 fs and 180 fs being the laser pulse duration, emitter/detector carrier recombination and collision time, respectively28; (2) determining reflections from graphene and the substrate using the Fresnel coefficients, the real part of the materials refractive index in the literature16 and a constant graphene conductivity value for CVD graphene grown on Cu foil7, 10, 29; and (3) applying the derived expressions to obtain conductivity for comparison against the conductivity defined in step 2. The use of only the real part of the materials refractive index is justified by the use of substrate materials such as sapphire and intrinsic Ge that have a negligible extinction coefficient in the relevant frequency range between 0.5 to 1 THz16. The simulation highlighted that, where there is a phase shift deliberately introduced between the sample reflection and the reference reflection, the slope of the derived conductivity spectra is no longer zero (see Supporting Information). This fact is exploited in the automated analysis of our measurement for phase correction. It should be noted that phase misalignment can be due to several reasons, such as fibre drifts and mechanical jittering of the optical delay stage, and that these problems are most severe in fibre-coupled terahertz systems as used in this study. Our phase correction method therefore shifts the acquired reflection pulse with respect to the reference measurement in the time-domain in order for the real part of the calculated conductivity spectra to have a slope close to zero. A single step shift in time here corresponds to a sample being placed approximately 2.5 μm with respect to the reference mirror, and a positive value means that the sample position is shifted in the direction of the incident beam. We note that the phase compensation scheme can be alternatively implemented in the frequency domain by multiplying the phase shift term6.

Transferred graphene on sapphire was measured with TPI at a step size of 200 μm, where 15 waveform traces were averaged to represent a measurement for one single pixel. Here the number of waveforms acquired for averaging is relatively small as would be the requirement for potential in-line applications. Before the measurement, however, terahertz reflection from the Al mirror placed nominally at the same position as the sample was used as reference measurement. An Al mirror generally works as an almost perfect reflector in the terahertz regime30. From the raster scanned measurement, the region of the substrate covered by graphene was isolated by masking the data with an intensity threshold value given that the graphene covered area corresponds to regions of higher reflectivity relative to the plain sapphire substrate. As the primary reflection was well separated from the first reflection in the time-domain, a time windowing function was used to process the acquired waveforms in the regions of interest in order to remove Fabry-Perot or etalon effects that would otherwise corrupt the conductivity measurement. A conductivity map was subsequently generated. In order to validate our proposed method, the same sample was scanned with THz-TDS operating in transmission mode. In particular, measurements were acquired with the Tera K15 T-Light (Menlo Systems GmbH, Germany) where the 60 mW pump pulse was focused to a 40 μm spot onto the terahertz photoconductive antenna, generating terahertz radiation with a beam diameter of approximately 1 mm at 1 THz. The sample was placed between the terahertz emitter and the detector at normal incidence without nitrogen purge, and data were acquired at an integration time of 10 ms at 200 μm step intervals. The integration time constant corresponded to an average of 700 waveform traces. The region of the substrate covered by graphene was again obtained by intensity masking, where graphene covered areas correspond to regions with a reduction in transmitted intensity due to the higher carrier absorption in graphene10. By performing the analysis detailed in recent literature7, 10, a comparative conductivity map for approximately the same sample area was generated.

Data availability

Additional data sets related to this publication are available from the Cambridge University data repository at

Results and Discussion

Graphene conductivity on sapphire substrate

Figure 2a–c shows the Raman graphene D/G ratio map, frequency distribution, and 2D/G ratio map of the transferred graphene on sapphire. Most measured points show a D/G ratio of roughly 10%, highlighting the presence of defects that have been introduced during the transfer procedure. The 2D/G ratio map shows an average value of more than 1 highlighting that the film is predominately monolayer graphene. The substrate refractive index measured in transmission mode TDS is 3.1, in close agreement with literature16. This value was then used to obtain flat conductivity spectra7, 10, 29 (see Supporting Information) without any phase correction, which in turn was used to generate the transmission conductivity map shown in Fig. 2d. The imaginary conductivity measured here is no longer negligible due to the fact that phase was not accounted for. Graphene film conductivity measurements with THz-TDS in transmission geometry have previously been benchmarked against micro-four point probe, micro Raman spectroscopy and optical imaging7. Here we use transmission measurements to validate our measurements in reflection geometry. The substrate refractive index measured in reflection mode was approximately 3. The slight discrepancy may be due to sample surface related imperfections, leading to scattering losses and measurement with a focused beam as opposed to a collimated beam31. In a manner similar to the literature7, 10, 29 and our transmission measurements, the graphene conductivity spectra on sapphire for a well-aligned reflection measurement have a real part characterised by a flat spectral response near to its DC value well below the Drude roll-off frequency or the inverse scattering time, while the imaginary part is close to zero (see Supporting Information). The spectrally resolved conductivity in turn can be represented by a single real number. It has been shown that phase misalignment would only significantly affect the imaginary part of the conductivity6 and therefore emphasis is placed only on the real part of the conductivity. As the conductivity spectra are approximately constant over the spectral range, the representative conductivity is taken as the average between 0.6 to 0.9 THz because at higher frequencies, the transmission conductivity spectra (see Supporting Information) become affected by water vapour absorption under ambient conditions. This effect become visibly more pronounced with increasing optical path length. Similarly, for the conductivity measurement in reflection mode, the representative pixel conductivity is taken as the average over the same spectral range of 0.6 to 0.9 THz.

Figure 2
Raman map of graphene on sapphire substrate, (a) D/G ratio map, (b) D/G frequency distribution, (c) Raman 2D/G map. For the same region a conductivity map of graphene on sapphire substrate was measured with THz-TDS between 0.6–0.9 THz ...

Figure 2 compares the conductivity map and histogram obtained with terahertz transmission and reflection mode TDS, respectively. The conductivity map acquired in reflection mode contains intermittent interlacing artefacts between alternate rows on the image. This is due to the small signal fluctuations that propagate to the conductivity calculation. Sources of signal fluctuations include fibre drifts, laser instability, optical and electronic noise. Nevertheless, qualitatively, an agreement between the two measurement geometries can be seen, for instance by looking at the regions of low local conductivity. When comparing the histogram in Fig. 2g, the conductivity frequency distribution is generally in agreement despite differences in the spot sizes of the THz-TDS systems used. In order to allow for a comparison that accounts for the differences in terahertz spot size, spatial filtering was applied to the reflection conductivity map to emulate a spot size approximately 2.4 times greater in order to generate a spatially averaged conductivity map and histogram shown in Fig. 2f and h, respectively. This results in a much closer agreement against the transmission measurements. The low conductivity values in the histogram correspond to the pixels positioned close to the graphene boundary and the slightly higher conductivity in reflection compared to transmission may be due to the averaging of the aforementioned signal fluctuations. The demonstration of the reflection measurement here means that, in principle, conductivity can be reasonably estimated directly from the first transmitted pulse, recently demonstrated6, as opposed to the first echo in a transmission measurement, even though the measurement robustness is lower for direct transmission analysis7, 10. Overall, we were able to validate our proposed approach with good agreement against established transmission measurements.

Graphene mobility

In order to demonstrate that our proposed method can be used to resolve graphene conductivity changes, we implemented simultaneous back-gating19, 32 via a graphene film transferred to a p-doped Si substrate (Fig. 1, see methods). Unlike previous work where complex transistor and Hall bar devices were realised to measure the mobility33, 34, here we attached a piece of Cu foil to the Si side of the graphene/wafer-stack and another Cu foil to the edge of the graphene. This rather simplistic setup defined the electrodes for back-gating and grounding the graphene film, respectively. The Fermi level of graphene was electrically tuned via the back-gate voltage supplied by a DC variable voltage supply (Keithley, Model 2400) where various different DC voltages between −160 to 160 V were applied. The back-gate leakage current was negligible. At a fixed point on graphene and at every 20 V back-gate voltage increments, reflection measurements were acquired on the TPI at 15 waveforms per average. The acquired waveforms were then analysed and the sequence of steps were then repeated for 2 other randomly selected locations on the graphene. The measurements were performed in ambient conditions at room temperature. Figure 3a shows a selection of the reflected waveforms, and the expanded views on the signal near the peaks against the applied back-gate voltage. From Equation 1, it is expected that terahertz reflection increases with increasing graphene conductivity induced by electrical gating. The reflection changes were not due to the influence of the substrate as no changes to the waveforms were observed when the terahertz spot was on the substrate under similar gating voltages. At Vg between 40 and 60 V, reflection is the weakest indicating that the Fermi energy at this gate voltage is closest to the Dirac point. At all other voltages below 40 V terahertz reflection increases monotonically with Vg as shown in Fig. 3b. Based on the acquired terahertz reflection waveforms, the real part of the frequency dependent conductivity was determined, as shown in Fig. 3b. This shows a strong dependence on the applied gate voltage. It can also be observed that there is a slight increase in the conductivity with increasing frequency. This may be due to a small degree of preferential back-scattering of charge carriers in the graphene8. Due to the aforementioned reasons, the spectrally resolved conductivity is represented by a single real-value conductivity value taken as the average conductivity between 0.4–0.9 THz. Figure 3c shows the measured conductivity for 3 randomly selected points on the graphene as a function of the applied back-gate voltage. Conductivity shows a linear dependence on the applied gate voltage in the range of −160  Vg  60 V representing the field-effect mobility of the hole-carriers. For Vg > 60 V, a slow conductivity change is observed indicating that electron-mobility is impaired in this graphene film. Linear curve-fitting was used to extract the line slope to determine hole-mobility. A R2 value of 0.986 was achieved in the fit that corresponded to a conductivity change of 0.0083 mS Vg −1. This value was then used to determine a mobility of 721 cm2V−1s−1, where Cox=εoxtox=11.5×109Fcm2 for the 300 nm thick SiO2 gate oxide. The extracted back-gated mobility value is comparable to the mobility values extracted from graphene transistor measurements of graphene films prepared in a similar manner33, 34. Air-exposed graphene on SiO2 is typically found to be p-doped34, consistent with the charge neutrality point VCNP at a positive bias34. Even though not demonstrated here, mobility mapping with terahertz reflection measurement across the sample could in principle be performed by raster scanning the sample at different back-gate voltages.

Figure 3
The measured terahertz reflection from a randomly selected position on the graphene on boron-doped Si/SiO2 substrate as function of applied back-gate voltage in (a) and an expanded view on the peak of the terahertz reflection in (b). The corresponding ...

Graphene conductivity mapping on Ge

Utilising a reflection rather than a transmission geometry allows contactless mapping of the graphene conductivity on a wider range of substrates. In terms of integrated graphene CVD manufacturing and process optimisation, a clear need is to probe graphene at each process step. We here so far discussed two examples of graphene conductivity measurements after graphene transfer on substrates like sapphire and Si/SiO2, but other substrates such as flexible polymers are in principle possible too provided sufficient contrast can be observed, which is generally the case as the real part of the refractive index is less than 2 and the extinction coefficient is negligible at less than 2 THz35, 36. This has been highlighted already for transmission measurements7. The challenge of adequate characterisation of the graphene still on the growth catalyst/support however remains. For instance, even individual Raman measurements of graphene on transition metals can be challenging, not to speak of mapping large graphene areas. As a first step towards direct contactless graphene conductivity mapping on technologically relevant growth substrates, we here focus on Ge substrates. As a starting point we first transferred graphene on intrinsic Ge(110), and then applied our method to graphene directly grown on Ge(110). The Ge (110) orientation is chosen because it yields higher quality graphene under the CVD process used here19, 37. For our initial analysis, the imaginary part of the complex refractive index of Ge can still be assumed to be negligible. For the undoped Ge(110) reference substrate, the graphene conductivity can be expected to be close to the values measured on sapphire, except some minor deviations based on, for instance, charge transfer between graphene and Ge38. Hence we can also assume that the real conductivity spectra remain constant over the THz spectral range7, 10, 29. For reflection measurements on highly doped Ge, the assumption of a negligible extinction coefficient will no longer hold39 and therefore the frequency dependent extinction coefficient needs to be accounted for in the conductivity analysis. Accurate measurement of the extinction coefficient in turn would require very high precision alignment between the sample and reference13. Analogous to our measurements of graphene on sapphire, the graphene film transferred on Ge was measured with TPI at a step size of 200 μm where 15 waveform traces were averaged. A conductivity map was then generated from the raster scanned measurement by applying the aforementioned analysis where the region of graphene coverage was obtained by intensity masking the data, where the graphene covered area corresponds to regions with higher reflection relative to Ge. The measured substrate refractive index of Ge was approximately 4 in close agreement with literature16. It should be noted that the contrast i.e. reflection change from bare Ge substrate to graphene covered area is approximately 5% for graphene with conductivity of 1.2 mS. This can be compared against graphene on sapphire where for the same conductivity, the contrast is approximately 10%, allowing a clear discrimination between graphene covered areas and the bare substrate. With a lower contrast, the measurement becomes more susceptible to noise and therefore we defined the threshold for intensity masking by considering an average of the substrate reflection. Figure 4a shows the measured conductivity map of transferred graphene on Ge(110) (see Supporting Information for corresponding optical microscope image). Here the blacked out pixels in the region of interest correspond to pixels where the conductivity was not computed because the reflection was below the defined threshold. For all other pixels, the average conductivity over 0.6 to 0.9 THz was estimated (see Supporting Information). On the conductivity map, there are visible spots of low conductivity in part due to the small signal fluctuations as a result of the lower contrast or from the graphene film, though the former is more likely. The conductivity histogram in Fig. 4b shows that the distribution does not closely follow a Gaussian function, but is centred around 1 mS in agreement with the conductivity value measured on sapphire support (Fig. 2). The lower average conductivity may be due to the aforementioned reduced contrast and hence increased sensitivity to the noise inherent to commercial fibre-coupled THz-TDS systems, as well as a smaller graphene covered area leading to the increased influence of boundary areas where lower graphene conductivity is generally observed. There is also qualitative agreement in the shape of the histogram compared to Fig. 2g where a sharp roll-off is observed for high conductivities as opposed to the lower conductivity values. Overall, the Ge measurements show that graphene film conductivity mapping by terahertz reflection spectroscopy on lightly doped substrates is possible, with the obtained conductivity values in agreement to our sapphire measurements. Note that the conductivity value is expected to vary slightly when comparing graphene on sapphire and Ge due to substrate induced doping38.

Figure 4
Conductivity map (a) of graphene film transferred on Ge(110) substrate measured with TPI and (b) the corresponding histogram. (c) The peak terahertz electric field map of patterned CVD graphene synthesised directly on Ge(110).

As final validation, we apply our method to graphene directly grown on Ge(110). After CVD on Ge of a uniform graphene coverage (see Methods), we subsequently used drop casted PMMA photoresist and oxygen plasma etching (with oxygen partial pressure 50 mbar at 50 W for 8 s) to define graphene patterns to allow for reference contrast against the plain underlying Ge. After oxygen plasma etching, the remaining photoresist was removed with acetone. We then measured with TPI at a step size of 75 μm where 50 waveform traces were averaged for one pixel. Figure 4c shows the peak electric field map where a contrast of approximately 5% can be observed between the graphene and the exposed underlying Ge. After graphene CVD, the Ge substrate appears to be highly doped (0.25 Ω·cm) as verified by standard conductivity measurements. This therefore means that although qualitative analysis of the intensity map is possible, a quantitative analysis to obtain the conductivity value cannot be reliably performed on this sample because the underlying assumption of negligible extinction coefficient is no longer valid. For the accurate determination of the extinction coefficient, high precision sample alignment would be required for substrate Drude model fitting that is outside the scope of this paper. The proposed technique, however, would be possible for lightly doped substrates (1–10 Ω·cm) where the extinction coefficient still remains negligible for frequencies greater than 0.5 THz as in the case of Si40. In principle, these substrates could also work with transmission THz-TDS, though the measurement may slightly overestimate the conductivity due to increased carrier absorption as the terahertz pulse makes a return path inside the substrate. The fact that there is an observable contrast in the reflection measurement is promising for future conductivity mapping of graphene synthesised directly on Ge. It can be expected that in samples where contrast is lower, either because of a higher substrate refractive index at terahertz frequencies, a lower graphene conductivity, or a combination thereof, the conductivity measurement in reflection geometry becomes increasingly susceptible to the effect of signal fluctuations and hence become less reliable. For instance, on growth substrates such as Cu2, no contrast could be observed relative to Cu for a sample of similar graphene conductivity because of the high Cu conductivity (59 MS/m) leading to a refractive index (~730) at least two orders of magnitude larger than that observed for sapphire or Ge at 1 THz. At the same time, it should also be pointed out that with the possibility of making changes to the graphene manufacturing process, such as intercalation of a micron-thick oxide layer41, these changes potentially could be useful to increase the measurement contrast for future in-line characterisation of graphene on e.g. Cu.


In summary, we have demonstrated the feasibility and potential of measuring the electrical conductivity of CVD graphene with THz-TDS in reflection geometry. Using terahertz transparent sapphire support, we have validated the technique against current state-of-the-art THz-TDS transmission measurements, where after taking into account the differences in terahertz spot sizes, we find a close agreement between the conductivity histograms. Using back-gated Si/SiO2 support, we have further demonstrated the sensitivity of the technique to resolve conductivity changes during electrostatic gating and hence the ability to directly determine graphene mobilities, whereby our values were consistent with standard electrical measurements. To illustrate that this technique has the potential as a tool for in-line graphene quality monitoring, we showed that the graphene conductivity can also be directly mapped on a substrate like Ge, where only half the contrast is seen relative to the sapphire substrate due to an approximately 30% increase in refractive index in the relevant frequency range. For the graphene films synthesised directly on Ge(110), the Ge substrate became highly doped after graphene synthesis and hence high precision alignment is required to accurately determine graphene conductivity. It should be noted that the proposed technique also assumes flat conductivity spectrum, which is not the case for all graphene samples, such as graphene grown on a single Cu crystal. For these samples, provided the shape of the conductivity spectra is known, a similar slope fitting method could be used. A more robust approach, however, would involve high precision alignment that is the subject of ongoing research. Our data shows that while measurement and analysis are understandably more complex and can be less robust when measuring in reflection compared to transmission, especially considering that there is less terahertz interaction with the sample, terahertz time-domain reflection spectroscopy is a highly interesting contactless, quantitative characterisation technique with clear potential to complement existing characterisation techniques for graphene and other related 2D materials and to open new opportunities for the rapid screening of large-area 2D crystals and films, which is crucial towards emerging applications and industrial development of these materials.

Electronic supplementary material


H.L. and J.A.Z. acknowledge financial support from the EPSRC (Grant No. EP/L019922/1). P.B.W. acknowledges EPSRC Cambridge NanoDTC EP/G037221/1. R.D., H.E.B. and D. R. acknowledge financial support from the EPSRC (Grant No. EP/J017671/1, Coherent Terahertz Systems). S.H. acknowledges funding from the EPSRC (Grant No. EP/K016636/1, GRAPHTED). H.L. and P.B.W also thank Dr Patrick Whelan for useful discussion.

Author Contributions

Author Contributions

H.L., P.B.W., J.A.Z. and S.H. conceived and designed the experiments. P.B.W. prepared all the samples and performed Raman and optical measurements. H.L. and P.B.W. performed reflection THz-TDS measurements. V.S.K., D.S.J. and R.D. performed transmission THz-TDS measurements. H.L. analysed the experimental results. H.L., P.B.W., J.A.Z., R.D. and S.H. interpreted the results. H.L. and P.B.W. wrote the main manuscript prepared all the figures. All authors reviewed the manuscript.


Competing Interests

The authors declare that they have no competing interests.


Hungyen Lin and Philipp Braeuninger-Weimer contributed equally to this work.

Electronic supplementary material

Supplementary information accompanies this paper at doi:10.1038/s41598-017-09809-7

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Contributor Information

Hungyen Lin,

Philipp Braeuninger-Weimer,


1. Things you could do with graphene. Nat Nano9, 737 (2014).
2. Hofmann S, Braeuninger-Weimer P, Weatherup RS. CVD-Enabled Graphene Manufacture and Technology. J Phys Chem Lett. 2015;6:2714–2721. doi: 10.1021/acs.jpclett.5b01052. [PMC free article] [PubMed] [Cross Ref]
3. Li XS, et al. Large-Area Synthesis of High-Quality and Uniform Graphene Films on Copper Foils. Science. 2009;324:1312–1314. doi: 10.1126/science.1171245. [PubMed] [Cross Ref]
4. Bae S, et al. Roll-to-roll production of 30-inch graphene films for transparent electrodes. Nat Nanotechnol. 2010;5:574–578. doi: 10.1038/nnano.2010.132. [PubMed] [Cross Ref]
5. Sensale-Rodriguez B, et al. Broadband graphene terahertz modulators enabled by intraband transitions. Nat Commun. 2012;3 doi: 10.1038/ncomms1787. [PubMed] [Cross Ref]
6. Whelan PR, et al. Robust mapping of electrical properties of graphene from terahertz time-domain spectroscopy with timing jitter correction. Opt Express. 2017;25:2725–2732. doi: 10.1364/OE.25.002725. [Cross Ref]
7. Buron JD, et al. Graphene Conductance Uniformity Mapping. Nano Letters. 2012;12:5074–5081. doi: 10.1021/nl301551a. [PubMed] [Cross Ref]
8. Buron JD, et al. Graphene mobility mapping. Sci Rep. 2015;5 doi: 10.1038/srep12305. [PMC free article] [PubMed] [Cross Ref]
9. Buron JD, et al. Terahertz wafer-scale mobility mapping of graphene on insulating substrates without a gate. Opt Express. 2015;23:30721–30729. doi: 10.1364/OE.23.030721. [PubMed] [Cross Ref]
10. Tomaino JL, et al. Terahertz imaging and spectroscopy of large-area single-layer graphene. Opt Express. 2011;19:141–146. doi: 10.1364/OE.19.000141. [PubMed] [Cross Ref]
11. Jeon T-I, Grischkowsky D. Characterization of optically dense, doped semiconductors by reflection THz time domain spectroscopy. Applied Physics Letters. 1998;72:3032–3034. doi: 10.1063/1.121531. [Cross Ref]
12. Hashimshony D, et al. Characterization of the electrical properties and thickness of thin epitaxial semiconductor layers by THz reflection spectroscopy. Journal of Applied Physics. 2001;90:5778–5781. doi: 10.1063/1.1412574. [Cross Ref]
13. Nagashima T, Hangyo M. Measurement of complex optical constants of a highly doped Si wafer using terahertz ellipsometry. Applied Physics Letters. 2001;79:3917–3919. doi: 10.1063/1.1426258. [Cross Ref]
14. Nashima S, Morikawa O, Takata K, Hangyo M. Measurement of optical properties of highly doped silicon by terahertz time domain reflection spectroscopy. Applied Physics Letters. 2001;79:3923–3925. doi: 10.1063/1.1413498. [Cross Ref]
15. Pashkin A, Kempa M, Němec H, Kadlec F, Kužel P. Phase-sensitive time-domain terahertz reflection spectroscopy. Review of Scientific Instruments. 2003;74:4711–4717. doi: 10.1063/1.1614878. [Cross Ref]
16. Grischkowsky D, Keiding S, Vanexter M, Fattinger C. Far-Infrared Time-Domain Spectroscopy with Terahertz Beams of Dielectrics and Semiconductors. J Opt Soc Am B. 1990;7:2006–2015. doi: 10.1364/JOSAB.7.002006. [Cross Ref]
17. Reina A, et al. Large Area, Few-Layer Graphene Films on Arbitrary Substrates by Chemical Vapor Deposition. Nano Letters. 2009;9:30–35. doi: 10.1021/nl801827v. [PubMed] [Cross Ref]
18. Braeuninger-Weimer P, Brennan B, Pollard AJ, Hofmann S. Understanding and Controlling Cu-Catalyzed Graphene Nucleation: The Role of Impurities, Roughness, and Oxygen Scavenging. Chemistry of Materials. 2016;28:8905–8915. doi: 10.1021/acs.chemmater.6b03241. [PMC free article] [PubMed] [Cross Ref]
19. Lee JH, et al. Wafer-Scale Growth of Single-Crystal Monolayer Graphene on Reusable Hydrogen-Terminated Germanium. Science. 2014;344:286–289. doi: 10.1126/science.1252268. [PubMed] [Cross Ref]
20. Shen YC, Taday PF. Development and application of terahertz pulsed imaging for nondestructive inspection of pharmaceutical tablet. IEEE J Sel Top Quant. 2008;14:407–415. doi: 10.1109/JSTQE.2007.911309. [Cross Ref]
21. Saleh, B. E. A. & Teich, M. C. Fundamentals of photonics. (Wiley, 1991).
22. Vartiainen EM, et al. Numerical phase correction method for terahertz time-domain reflection spectroscopy. Journal of Applied Physics. 2004;96:4171–4175. doi: 10.1063/1.1786345. [Cross Ref]
23. Born, M. & Wolf, E. Principles of optics: electromagnetic theory of propagation, interference and diffraction of light. 7th edn, (Cambridge University Press, 1999).
24. Tinkham M. Energy Gap Interpretation of Experiments on Infrared Transmission through Superconducting Films. Phys Rev. 1956;104:845–846. doi: 10.1103/PhysRev.104.845. [Cross Ref]
25. Mcknight SW, Stewart KP, Drew HD, Moorjani K. Wavelength-Independent Anti-Interference Coating for the Far-Infrared. Infrared Phys. 1987;27:327–333. doi: 10.1016/0020-0891(87)90074-1. [Cross Ref]
26. Thoman A, Kern A, Helm H, Walther M. Nanostructured gold films as broadband terahertz antireflection coatings. Phys Rev B. 2008;77 doi: 10.1103/PhysRevB.77.195405. [Cross Ref]
27. Dhillon SS, et al. The 2017 terahertz science and technology roadmap. J Phys D Appl Phys. 2017;50 doi: 10.1088/1361-6463/50/4/043001. [Cross Ref]
28. Duvillaret L, Garet F, Roux JF, Coutaz JL. Analytical modeling and optimization of terahertz time-domain spectroscopy experiments using photoswitches as antennas. IEEE J Sel Top Quant. 2001;7:615–623. doi: 10.1109/2944.974233. [Cross Ref]
29. Buron JD, et al. Electrically Continuous Graphene from Single Crystal Copper Verified by Terahertz Conductance Spectroscopy and Micro Four-Point Probe. Nano Letters. 2014;14:6348–6355. doi: 10.1021/nl5028167. [PubMed] [Cross Ref]
30. Gatesman AJ, Giles RH, Waldman J. High-Precision Reflectometer for Submillimeter Wavelengths. J Opt Soc Am B. 1995;12:212–219. doi: 10.1364/JOSAB.12.000212. [Cross Ref]
31. Naftaly, M. An international intercomparison of THz time-domain spectrometers. Int Conf Infra Milli (2016).
32. Ren L, et al. Terahertz and Infrared Spectroscopy of Gated Large-Area Graphene. Nano Letters. 2012;12:3711–3715. doi: 10.1021/nl301496r. [PubMed] [Cross Ref]
33. Alexander-Webber JA, et al. Encapsulation of graphene transistors and vertical device integration by interface engineering with atomic layer deposited oxide. 2d Mater. 2017;4:1–9.
34. Van Veldhoven ZA, Alexander-Webber JA, Sagade AA, Braeuninger-Weimer P, Hofmann S. Electronic properties of CVD graphene: The role of grain boundaries, atmospheric doping, and encapsulation by ALD. Phys Status Solidi B. 2016;253:2321–2325. doi: 10.1002/pssb.201600255. [Cross Ref]
35. Cunningham PD, et al. Broadband terahertz characterization of the refractive index and absorption of some important polymeric and organic electro-optic materials. Journal of Applied Physics. 2011;109 doi: 10.1063/1.3549120. [Cross Ref]
36. Wietzke S, et al. Terahertz spectroscopy on polymers: A review of morphological studies. J Mol Struct. 2011;1006:41–51. doi: 10.1016/j.molstruc.2011.07.036. [Cross Ref]
37. Dai JY, et al. How Graphene Islands Are Unidirectionally Aligned on the Ge(110) Surface. Nano Letters. 2016;16:3160–3165. doi: 10.1021/acs.nanolett.6b00486. [PubMed] [Cross Ref]
38. Cavallo F, et al. Exceptional Charge Transport Properties of Graphene on Germanium. ACS Nano. 2014;8:10237–10245. doi: 10.1021/nn503381m. [PubMed] [Cross Ref]
39. Auston DH, Cheung KP. Coherent Time-Domain Far-Infrared Spectroscopy. J Opt Soc Am B. 1985;2:606–612. doi: 10.1364/JOSAB.2.000606. [Cross Ref]
40. Hangyo M, Nagashima T, Nashima S. Spectroscopy by pulsed terahertz radiation. Meas Sci Technol. 2002;13:1727–1738. doi: 10.1088/0957-0233/13/11/309. [Cross Ref]
41. Wang RZ, et al. Catalyst Interface Engineering for Improved 2D Film Lift-Off and Transfer. ACS Appl Mater Inter. 2016;8:33072–33082. doi: 10.1021/acsami.6b11685. [PMC free article] [PubMed] [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group