1. Introduction
Since graphene was first made available by mechanical exfoliation of graphite [
1], this two-dimensional material has continued to attract a significant amount of attention due to its remarkable mechanical, thermal, optical, and electronic properties [
2,
3]. It is predicted to have applications in a vast number of areas, including electronics, photonics, optoelectronics, and transistor- and composite-material manufacturing [
4,
5,
6,
7], and graphene material and device manufacturing continues to expand and thus creates a stronger need for reliable, fast, and non-destructive quality inspection techniques [
8].
Terahertz time-domain spectroscopy (THz-TDS) has been demonstrated to be a versatile technique for the non-destructive retrieval of the electrical properties of graphene [
9,
10,
11,
12]. While commercial THz-TDS devices are capable of providing both accurate and precise measurements, they are not easily implemented in production facilities, while THz systems aimed at high-speed commercial applications do not provide the deep level of information needed to establish key performance characteristics of graphene films. Other techniques for the retrieval of the electrical properties of graphene include microwave impedance microscopy and micro four-point probe measurements, where the former, such as THz-TDS, is a non-contact technique, and the latter requires direct physical contact with the sample [
13].
THz-TDS was first introduced in 1989 through the utilization of photoconductive antennas (PCAs) and femtosecond (fs) lasers [
14]. In a typical setup, a PCA emitter is optically excited by the laser pulse, which, through charge carrier acceleration, emits THz radiation that subsequently is coherently detected via photoconductive sampling in a PCA detector [
15]. Other methods for THz generation and detection include optical rectification in nonlinear crystals and third-order nonlinear effects in air plasmas, both methods utilizing fs lasers as well [
16,
17].
Instead of generating THz radiation using pulsed lasers, it was proposed in 1999 to excite PCAs with a temporally incoherent light source and measure the cross-correlation between the emitted incoherent radiation and the incoherent pump light [
18,
19]. Similar systems have since been fabricated employing different temporally incoherent light sources such as superluminescent diodes or multi-mode laser diodes [
20,
21,
22,
23]. In this work, we present a novel THz cross-correlation spectroscopy (CCS) device capable of obtaining results similar to commercial systems, while being smaller, easier to operate, and less complex, which allows for significantly cheaper production. We have previously demonstrated CCS applications of thickness measurements based on time-of-flight algorithms [
24,
25], and here we focus on the spectroscopic capabilities of the CCS system in sensing applications.
2. THz Cross-Correlation Spectroscopy
The main components of a typical CCS setup using PCAs is shown schematically in
Figure 1. A light source source with low temporal coherence generates an intensity
, which is split into an emitter and a detector arm, respectively. In the emitter arm, a biased PCA is illuminated, which generates THz radiation
with a similar low temporal coherence as the light source. Note that at any two points of equal optical distance from the beamsplitter, the incoherent light intensities in the two arms are identical.
Both emitter and detector operate by the standard principles of THz-frequency photomixing [
26]. Assuming a spatially uniform intensity profile and ignoring antenna effects, the emitted THz field is proportional to the temporal derivative of the photocurrent induced by free carriers generated by the incoherent light. The photocurrent is the product of the photoconductivity
and the constant emitter bias voltage
. Taking the impulse response
of the emitter into account, we find that
In the detector, as in the emitter, the incoherent light generates free carriers inducing photoconductance
,
where
is the impulse response of the detector antenna, and the incoming THz field drives a photocurrent
. Effectively, when integrating using, e.g., lock-in detection over a time period much greater than any characteristic time in the system, a current
proportional to the cross-correlation
between
and
is measured. Inserting Equations (
1) and (
2) in Equation (
3) and applying the Fourier transform, we find that
which is a coherent spectrum independent of the phase of the light source. Despite the use of an incoherent light source, the fact that the cross-correlation depends on the absolute square of the intensity
means that the measured cross-correlation has a pulse-like shape, when transformed back to the time domain. By properly choosing the bandwidth of the light source, the intensity
of it will contain frequency-mixed terms in the THz range.
3. Extraction of Sheet Conductivity of Thin Films in the THz Range
Upon insertion of a sample between the emitter and the detector, the field is modified. The interaction is described in the frequency domain by a transfer function
, which modifies the measured photocurrent to
The transfer function is set up by analyzing how plane waves of light propagate through the sample in a transmission experiment as displayed in
Figure 2. Especially in situations where the sample is placed on a thin, low-index substrate such as a polymer film, the inclusion of internal reflections becomes important for a quantitative analysis of the conductive properties of the thin film as the directly transmitted signals and the internal reflections may overlap in time [
27].
The total transmitted field
is the sum of the directly transmitted field and the fields having undergone
n internal reflections, i.e.,
where
and
are the standard Fresnel reflection and transmission coefficients for reflection/transmission between media
i and
j and
is the acquired phase through medium 2 of thickness
d and refractive index
.
When a thin film of a conducting material such as graphene is present on the interface between medium 1 and 2, the relevant transmission and reflection coefficients are described by the Tinkham equations [
28] given by
where
is the free-space impedance and
is the complex Drude conductivity given by
where
and
are the material-specific DC conductivity and momentum relaxation time, respectively.
To extract the conductivity experimentally, a sample measurement
of graphene on a substrate and a reference measurement
of the substrate are taken. When divided,
and
can be extracted by fitting the resulting expression
where subscripts
s, a, and
f refer to substrate, air, and film, respectively, using either deterministic or stochastic fitting algorithms such as basin hopping or differential evolution [
29,
30].
4. Experimental Setup
Two transmission mode experiments on the same sample of graphene at normal angle of incidence have been executed. The experimental setup is shown in
Figure 3a. The first experiment is performed with a commercial THz-TDS system (TOPTICA Teraflash Smart) and serves as a high-precision benchmark. The second experiment is conducted with our THz CCS device. A fiber-coupled C-band amplified spontaneous emission (ASE) light source is split into two arms. As displayed in
Figure 3b, the ASE output covers the range from 1525 to 1565 nm, which allows for the generation of photomixed frequencies of up to 5 THz. Due to the low-pass characteristics of the emitter and detector PCAs, bandwidths higher than 2 THz are not currently obtainable with the used component. Each beam is sent through one of the two fiber stretchers, which together are capable of creating a delay of 100 picoseconds (ps), before being coupled to an InGaAs-based emitter antenna and an InGaAs-based detector antenna, respectively (TOPTICA InGaAs Photomixers). The cross-correlation is measured via lock-in detection.
Key metrics for the used devices as operated are summarized in
Table 1. The TDS device has a bandwidth of 2.3 THz compared to 1.4 THz and a dynamic range of 50 dB compared to 58 dB of the CCS device at acquisition times of 3 s for the TDS device and 30 s for the CCS device. Due to the difference in acquisition time, the metrics cannot be directly compared—a relatively low acquisition time for the TDS device was used since it resulted in satisfactory results. In
Figure 4, example time traces (left) and spectra (right) of both devices are shown. A Tukey window function with a taper length of 5 ps has been applied to both traces. The CCS trace does not converge to zero in the end of the time window, which influences the spectral analysis.
Additionally, the spot sizes were determined through knife-edge scans in ambient air, shown in
Figure 5. The profiles were extracted from the series of full time-domain traces
as a function of the knife edge position
x by calculating the energy contained in each trace proportional to
. The beam intensity profile was then found by the derivative,
.
The spot sizes are found to be 1.3 mm for the TDS device and 5.6 mm for the CCS device, where the spot size is defined as twice the width extracted from fitting the beam intensity profile with a Gaussian function of the form
where
A is a normalization constant,
x is the spatial position, and
w is the beam width.
5. Sample Fabrication
The graphene sample is grown by chemical vapor deposition (CVD) on copper foil and then transferred to a transparent PET substrate with standard wet transfer techniques [
31]. During the transfer process, a polymer film model from a 1.3 wt.% solution of 2200 K Poly(methyl methacrylate) (PMMA) is used as a support. The sample is then placed in acetone for 1 h to remove the PMMA and finally cleaned by isopropanol.
Figure 6 shows an optical microscopy image (left) and the Raman spectra of the PET substrate and the graphene on top of the PET substrate (right). The peaks of the G and 2D bands located at 1588 cm
and 2720 cm
in the spectra are the most prominent features in the Raman spectra of graphene, which confirms the presence of graphene in the sample [
32]. The 2D band has a Lorentzian full width at half maximum of 32 cm
, which is consistent with a sample of 2–3 layers of graphene [
33,
34].
For better extraction of the graphene conductivity parameters, the thickness and dielectric permittivity of the substrate have been characterized experimentally with the TDS device and determined to be 208
m thick and well described (
) by a Debye model given by
with
,
and
fs. Many polymers display similar behavior in the low THz range, including a slightly decreasing real part of the permittivity and an increasing absorption towards higher frequencies [
35]. The real and imaginary parts of the extracted dielectric permittivity are shown in
Figure 7. The oscillatory feature at 0.8–1.3 THz was not observed in previous measurements [
36], and we attribute it to an incomplete cancellation of multiple reflections during the data extraction.
The Debye model describes the relaxation processes and displays the same qualitative behavior of the dielectric function. Hence, it is frequently used to model the dielectric properties of polymers in this range [
37,
38].
6. Experimental Results
A 9 cm × 9 cm sheet of graphene was raster-scanned in ambient air in a serpentine pattern with a stepsize of 3 mm, corresponding to 961 single-pixel measurements. For each pixel, the Drude conductivity was extracted by fitting with Equation (
10) using the stochastic optimization algorithm differential evolution [
39]. An example of the extracted conductivity using each device is shown in
Figure 8, where errors denote 68% confidence bounds. The conductivity spectra shown in
Figure 8 are not recorded at the same position on the sample, and hence the absolute values cannot be directly compared.
Since the TDS system has a small spot size compared to the pixel dimensions, >99.99% of the THz radiation is contained within each pixel, assuming a Gaussian beam profile. For the CCS device with a much larger spot size, only 35.6% of the radiation is contained, leading the extracted conductivities to be a weighted mean of the pixel in question and its neighbors. To provide a reasonable basis for comparison, the pixel values measured using the TDS device are therefore convolved with a Gaussian filter with a width equal to the beam width of the CCS device.
The spatial maps of the measured conductivities using both devices are shown in
Figure 9. The CCS measurements display the same conductivity pattern of spotted areas of higher and lower conductivity as the TDS measurements, and the measured DC conductivities have a linear correlation of
. The left-most vertical line corresponds to an area of the substrate without graphene, as indicated on the TDS conductivity plot. The CCS device is more biased towards measuring lower conductivities than the TDS device, and the bias displays spatial correlations as seen in the right panel of
Figure 9. Specifically, the mean error is
and the standard deviation on the relative errors given by
is
. The reason for the bias is discussed in the next section.
Since the transmission decreases when the DC conductivity increases, the peak-to-peak value is correlated with the DC conductivity. In
Figure 10, the DC conductivity determined with TDS is shown to display a strong linear correlation with the peak-to-peak values of the CCS measurements with a linear correlation of 89%. If the peak-to-peak values are used to predict
with a linear model (red) instead of retrieving it with a fitting procedure, the mean error is reduced to
and the standard deviation of the relative errors given by Equation (
13) is reduced to
. Similarly to the relative errors on the DC conductivities extracted with the fitting procedure shown in
Figure 9, the errors on the DC conductivities extracted with peak-to-peak procedure display spatial correlations as shown by coherent areas of positive or negative relative error.
7. Discussion of Results
The CCS device is capable of measuring the DC conductivity of graphene with a standard deviation on the error distribution of compared to the TDS benchmark. As shown in the previous section, patterns of high and low conductivity areas are captured and match the areas measured with the TDS device. However, as previously mentioned, the DC conductivities extracted with the fitting procedure using the CCS device are biased towards lower values. This bias is believed to primarily stem from drift in the optical path length in the fiber stretchers, where the length of the true time window increases over longer periods of operation, which is not reflected on the x-axis.
We believe that correctly compensating for the drift can remove the bias. This is currently being explored.
In the previous section, we demonstrated how the DC conductivity can be determined with higher precision from the peak-to-peak value of the trace. In
Figure 11, the linear correlation between the peak-to-peak value and the DC conductivity has been simulated for an identical experimental setup. For each value of the substrate thickness and refractive index, the corresponding pairs of
and peak-to-peak values are calculated for a simulated pulse with a bandwidth of 1 THz,
, and
, linearly varying between 0 and 12 mS. Each pixel value is then computed as the linear correlation between the peak-to-peak and
values, showing that it is possible to achieve correlations higher than 99%. The red cross denotes the point
, which is close to the settings at which the experiments presented here have been conducted.
Extracting from the peak-to-peak values removes the need to acquire the entire trace since the peak-to-peak value can be determined using only a fraction of the trace points. Since each CCS-trace consists of 200 data points, the acquisition time can ideally be reduced by up to a factor of 100 with a very moderate increase in measurement error merely by only scanning the maximum and minimum of the trace. Furthermore, it removes the error related to the signal not settling to zero within the limited time window and other issues with fit-based extraction procedures. At such rates, inline conductivity measurements of graphene become not only possible but viable.
8. Conclusions
We have demonstrated that our novel THz-CCS device is capable of measuring the electrical properties of graphene. The measurements have been compared to what is achievable with a state-of-the-art commercial THz-TDS device, resulting in a standard deviation on the relative errors of . Combined with a modern, industry-oriented system design, THz-CCS devices are becoming an attractive option for inline, non-contact, high throughput, high-precision quality assurance for 2D semiconductor materials.
Author Contributions
Conceptualization, P.B., S.J.L. and P.U.J.; data curation, B.H.M.; formal analysis, B.H.M.; funding acquisition, S.J.L. and P.U.J.; investigation, B.H.M. and T.B.; methodology, B.H.M. and T.B.; project administration, T.B.; resources, T.B., J.J. and S.J.L.; software, B.H.M.; supervision, P.U.J.; validation, B.H.M.; visualization, B.H.M. and T.B.; writing—original draft, B.H.M.; and writing—review & editing, B.H.M., T.B., S.J.L. and P.U.J. All authors have read and agreed to the published version of the manuscript.
Funding
We acknowledge partial support from the Independent Research Fund Denmark (project 9040-00360B, ULTRA-TED), the Innovation Fund Denmark (project TRIM), and DTU Elektro.
Data Availability Statement
Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.
Acknowledgments
The authors thank Mads Ehrhorn and Oscar Garcia Garcia for fruitful discussions related to the work described here.
Conflicts of Interest
The authors declare no conflict of interest.
Abbreviations
The following abbreviations are used in this manuscript:
THz-CCS | Terahertz cross-correlation spectroscopy |
THz-TDS | Terahertz time domain spectroscopy |
ASE | Amplified spontaneous emission |
PCA | Photoconductive antenna |
CVD | Chemical vapor deposition |
PMMA | Poly(methyl methacrylate) |
References
- Novoselov, K.S.; Geim, A.K.; Morozov, S.V.; Jiang, D.; Zhang, Y.; Dubonos, S.V.; Grigorieva, I.V.; Firsov, A.A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666–669. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Castro Neto, A.H.; Guinea, F.; Peres, N.M.R.; Novoselov, K.S.; Geim, A.K. The electronic properties of graphene. Rev. Mod. Phys. 2009, 81, 109–162. [Google Scholar] [CrossRef] [Green Version]
- Das Sarma, S.; Adam, S.; Hwang, E.H.; Rossi, E. Electronic transport in two-dimensional graphene. Rev. Mod. Phys. 2011, 83, 407–470. [Google Scholar] [CrossRef] [Green Version]
- Avouris, P.; Dimitrakopoulos, C. Graphene: Synthesis and applications. Mater. Today 2012, 15, 86–97. [Google Scholar] [CrossRef]
- Moreno, C.; Vilas-Varela, M.; Kretz, B.; Garcia-Lekue, A.; Costache, M.V.; Paradinas, M.; Panighel, M.; Ceballos, G.; Valenzuela, S.O.; Peña, D.; et al. Bottom-up synthesis of multifunctional nanoporous graphene. Science 2018, 360, 199–203. [Google Scholar] [CrossRef] [Green Version]
- Bonaccorso, F.; Sun, Z.; Hasan, T.; Ferrari, A.C. Graphene photonics and optoelectronics. Nat. Photonics 2010, 4, 611–622. [Google Scholar] [CrossRef] [Green Version]
- Stankovich, S.; Dikin, D.A.; Dommett, G.H.B.; Kohlhaas, K.M.; Zimney, E.J.; Stach, E.A.; Piner, R.D.; Nguyen, S.T.; Ruoff, R.S. Graphene-based composite materials. Nature 2006, 442, 282–286. [Google Scholar] [CrossRef]
- Backes, C.; Abdelkader, A.M.; Alonso, C.; Andrieux-Ledier, A.; Arenal, R.; Azpeitia, J.; Balakrishnan, N.; Banszerus, L.; Barjon, J.; Bartali, R.; et al. Production and processing of graphene and related materials. 2D Mater. 2020, 7, 022001. [Google Scholar] [CrossRef]
- Buron, J.D.; Petersen, D.H.; Bøggild, P.; Cooke, D.G.; Hilke, M.; Sun, J.; Whiteway, E.; Nielsen, P.F.; Hansen, O.; Yurgens, A.; et al. Graphene Conductance Uniformity Mapping. Nano Lett. 2012, 12, 5074–5081. [Google Scholar] [CrossRef]
- Buron, J.D.; Pizzocchero, F.; Jessen, B.S.; Booth, T.J.; Nielsen, P.F.; Hansen, O.; Hilke, M.; Whiteway, E.; Jepsen, P.U.; Bøggild, P.; et al. Electrically Continuous Graphene from Single Crystal Copper Verified by Terahertz Conductance Spectroscopy and Micro Four-Point Probe. Nano Lett. 2014, 14, 6348–6355. [Google Scholar] [CrossRef] [Green Version]
- Bøggild, P.; Mackenzie, D.M.A.; Whelan, P.R.; Petersen, D.H.; Buron, J.D.; Zurutuza, A.; Gallop, J.; Hao, L.; Jepsen, P.U. Mapping the electrical properties of large-area graphene. 2D Mater. 2017, 4, 042003. [Google Scholar] [CrossRef] [Green Version]
- Whelan, P.R.; Zhou, B.; Bezencenet, O.; Shivayogimath, A.; Mishra, N.; Shen, Q.; Jessen, B.S.; Pasternak, I.; Mackenzie, D.M.A.; Ji, J.; et al. Case studies of electrical characterisation of graphene by terahertz time-domain spectroscopy. 2D Mater. 2021, 8, 22003. [Google Scholar] [CrossRef]
- Whelan, P.R.; Panchal, V.; Petersen, D.H.; Mackenzie, D.M.A.; Melios, C.; Pasternak, I.; Gallop, J.; Østerberg, F.W.; Jepsen, P.U.; Strupinski, W.; et al. Electrical Homogeneity Mapping of Epitaxial Graphene on Silicon Carbide. ACS Appl. Mater. Interfaces 2018, 10, 31641–31647. [Google Scholar] [CrossRef] [Green Version]
- Exter, M.v.; Fattinger, C.; Grischkowsky, D. Terahertz time-domain spectroscopy of water vapor. Opt. Lett. 1989, 14, 1128–1130. [Google Scholar] [CrossRef] [PubMed]
- Jepsen, P.; Jacobsen, R.; Keiding, S. Generation and detection of terahertz pulses from biased semiconductor antennas. J. Opt. Soc. Am. B-Opt. Phys. 1996, 13, 2424–2436. [Google Scholar] [CrossRef]
- Xu, L.; Zhang, X.; Auston, D.H. Terahertz beam generation by femtosecond optical pulses in electro-optic materials. Appl. Phys. Lett. 1992, 61, 1784–1786. [Google Scholar] [CrossRef]
- Dai, J.; Xie, X.; Zhang, X.C. Detection of Broadband Terahertz Waves with a Laser-Induced Plasma in Gases. Phys. Rev. Lett. 2006, 97, 103903. [Google Scholar] [CrossRef]
- Morikawa, O.; Tonouchi, M.; Hangyo, M. A novel cross-correlation spectroscopy in sub-THz region using incoherent light source. In Proceedings of the Technical Digest. CLEO/Pacific Rim ’99. Pacific Rim Conference on Lasers and Electro-Optics (Cat. No.99TH8464), Seoul, Republic of Korea, 30 August–3 September 1999; Volume 4, pp. 1259–1260. [Google Scholar] [CrossRef]
- Morikawa, O.; Tonouchi, M.; Hangyo, M. A cross-correlation spectroscopy in subterahertz region using an incoherent light source. Appl. Phys. Lett. 2000, 76, 1519–1521. [Google Scholar] [CrossRef]
- Scheller, M.; Koch, M. Terahertz quasi time domain spectroscopy. Opt. Express 2009, 17, 17723–17733. [Google Scholar] [CrossRef]
- Brenner, C.; Hofmann, M.; Scheller, M.; Shakfa, M.K.; Koch, M.; Mayorga, I.C.; Klehr, A.; Erbert, G.; Trankle, G. Compact diode-laser-based system for continuous-wave and quasi-time-domain terahertz spectroscopy. Opt. Lett. 2010, 35, 3859–3861. [Google Scholar] [CrossRef]
- Probst, T.; Rehn, A.; Koch, M. Compact and low-cost THz QTDS system. Opt. Express 2015, 23, 21972–21982. [Google Scholar] [CrossRef] [PubMed]
- Molter, D.; Kolano, M.; Freymann, G.v. Terahertz cross-correlation spectroscopy driven by incoherent light from a superluminescent diode. Opt. Express 2019, 27, 12659–12665. [Google Scholar] [CrossRef] [PubMed]
- Hansen, R.E.; Bæk, T.; Lange, S.L.; Israelsen, N.M.; Mäntylä, M.; Bang, O.; Petersen, C.R. Non-Contact Paper Thickness and Quality Monitoring Based on Mid-Infrared Optical Coherence Tomography and THz Time Domain Spectroscopy. Sensors 2022, 22, 1549. [Google Scholar] [CrossRef] [PubMed]
- Knarreborg, J.Ø.; Hjortshøj-Nielsen, J.; Mølvig, B.H.; Bæk, T.; Jepsen, P.U.; Lange, S.J. Measuring Maritime Paint Thickness under Water Using THz Cross-Correlation Spectroscopy. Appl. Sci. 2022, 12, 12397. [Google Scholar] [CrossRef]
- Preu, S.; Döhler, G.H.; Malzer, S.; Wang, L.J.; Gossard, A.C. Tunable, continuous-wave Terahertz photomixer sources and applications. J. Appl. Phys. 2011, 109, 061301. [Google Scholar] [CrossRef] [Green Version]
- Whelan, P.R.; Shen, Q.; Luo, D.; Wang, M.; Ruoff, R.S.; Jepsen, P.U.; Bøggild, P.; Zhou, B. Reference-free THz-TDS conductivity analysis of thin conducting films. Opt. Express 2020, 28, 28819–28830. [Google Scholar] [CrossRef]
- Glover, R.E.; Tinkham, M. Conductivity of Superconducting Films for Photon Energies between 0.3 and 40 k T c. Phys. Rev. 1957, 108, 243–256. [Google Scholar] [CrossRef]
- Storn, R.; Price, K. Differential Evolution – A Simple and Efficient Heuristic for global Optimization over Continuous Spaces. J. Glob. Optim. 1997, 11, 341–359. [Google Scholar] [CrossRef]
- Wales, D.J.; Doye, J.P.K. Global Optimization by Basin-Hopping and the Lowest Energy Structures of Lennard-Jones Clusters Containing up to 110 Atoms. J. Phys. Chem. A 1997, 101, 5111–5116. [Google Scholar] [CrossRef] [Green Version]
- Reina, A.; Jia, X.; Ho, J.; Nezich, D.; Son, H.; Bulovic, V.; Dresselhaus, M.S.; Kong, J. Large Area, Few-Layer Graphene Films on Arbitrary Substrates by Chemical Vapor Deposition. Nano Lett. 2009, 9, 30–35. [Google Scholar] [CrossRef]
- Malard, L.M.; Pimenta, M.A.; Dresselhaus, G.; Dresselhaus, M.S. Raman spectroscopy in graphene. Phys. Rep. 2009, 473, 51–87. [Google Scholar] [CrossRef]
- Li, Z.; Kinloch, I.A.; Young, R.J.; Novoselov, K.S.; Anagnostopoulos, G.; Parthenios, J.; Galiotis, C.; Papagelis, K.; Lu, C.Y.; Britnell, L. Deformation of Wrinkled Graphene. ACS Nano 2015, 9, 3917–3925. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Li, X.; Cai, W.; An, J.; Kim, S.; Nah, J.; Yang, D.; Piner, R.; Velamakanni, A.; Jung, I.; Tutuc, E.; et al. Large-Area Synthesis of High-Quality and Uniform Graphene Films on Copper Foils. Science 2009, 324, 1312–1314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
- Wietzke, S.; Jansen, C.; Reuter, M.; Jung, T.; Kraft, D.; Chatterjee, S.; Fischer, B.M.; Koch, M. Terahertz spectroscopy on polymers: A review of morphological studies. J. Mol. Struct. 2011, 1006, 41–51. [Google Scholar] [CrossRef]
- Jin, Y.S.; Kim, G.J.; Jeon, S.G. Terahertz dielectric properties of polymers. J. Korean Phys. Soc. 2006, 49, 513–517. [Google Scholar]
- Nagai, N.; Fukasawa, R. Abnormal dispersion of polymer films in the THz frequency region. Chem. Phys. Lett. 2004, 388, 479–482. [Google Scholar] [CrossRef]
- Chen, X.L.; Chen, X.; Zhu, Z.Y. THz-TDS spectra study of polymer materials with different polarity. J. Infrared Millim. Waves 2013, 32, 150–154. [Google Scholar] [CrossRef]
- Virtanen, P.; Gommers, R.; Oliphant, T.E.; Haberland, M.; Reddy, T.; Cournapeau, D.; Burovski, E.; Peterson, P.; Weckesser, W.; Bright, J.; et al. SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python. Nat. Methods 2020, 17, 261–272. [Google Scholar] [CrossRef] [Green Version]
Figure 1.
A THz-CCS experiment utilizing PCAs and a temporally incoherent light source for detection and emission. The green trace on the lower black line indicates a typical electrical field profile as a function of time, extracted from the THz-CCS measurements.
Figure 1.
A THz-CCS experiment utilizing PCAs and a temporally incoherent light source for detection and emission. The green trace on the lower black line indicates a typical electrical field profile as a function of time, extracted from the THz-CCS measurements.
Figure 2.
Sketch of different contributions to the detected THz signal from multiple internal reflections. The gray box indicates the substrate and the green sheet the sample, such as graphene.
Figure 2.
Sketch of different contributions to the detected THz signal from multiple internal reflections. The gray box indicates the substrate and the green sheet the sample, such as graphene.
Figure 3.
(a) Illustration of the transmission geometry. The THz light is focused and transmitted through a sample of graphene on a dielectric substrate. (b) C-band ASE output power density.
Figure 3.
(a) Illustration of the transmission geometry. The THz light is focused and transmitted through a sample of graphene on a dielectric substrate. (b) C-band ASE output power density.
Figure 4.
Traces (
left) and spectra (
right) of the commercial THz-TDS system (black) and the novel THz-CCS device (red). The TDS trace has been offset. Key metrics are summarized in
Table 1.
Figure 4.
Traces (
left) and spectra (
right) of the commercial THz-TDS system (black) and the novel THz-CCS device (red). The TDS trace has been offset. Key metrics are summarized in
Table 1.
Figure 5.
Knife edge scans of the TDS device (blue) and the CCS device (yellow) with Gaussian fits to determine beam spot sizes. Black dotted lines indicate the spot sizes, which are found to be 1.3 mm for the TDS device and 5.6 mm for the CCS device.
Figure 5.
Knife edge scans of the TDS device (blue) and the CCS device (yellow) with Gaussian fits to determine beam spot sizes. Black dotted lines indicate the spot sizes, which are found to be 1.3 mm for the TDS device and 5.6 mm for the CCS device.
Figure 6.
Optical microscopy image (left) and Raman spectra (right) of the PET substrate and the graphene sample on top. The annotated peaks of the G and 2D bands of graphene are consistent with the presence of bi- or trilayer graphene.
Figure 6.
Optical microscopy image (left) and Raman spectra (right) of the PET substrate and the graphene sample on top. The annotated peaks of the G and 2D bands of graphene are consistent with the presence of bi- or trilayer graphene.
Figure 7.
Extracted real (top) and imaginary (bottom) parts of the dielectric permittivity of the PET substrate supporting the graphene films. The fitted values are given by a Debye model with parameters , , and fs.
Figure 7.
Extracted real (top) and imaginary (bottom) parts of the dielectric permittivity of the PET substrate supporting the graphene films. The fitted values are given by a Debye model with parameters , , and fs.
Figure 8.
Examples of real (blue) and imaginary (yellow) parts of the conductivity determined with the CCS system (
left) and the TDS system (
right) through fitting Equation (
10) with a stochastic optimization algorithm.
Figure 8.
Examples of real (blue) and imaginary (yellow) parts of the conductivity determined with the CCS system (
left) and the TDS system (
right) through fitting Equation (
10) with a stochastic optimization algorithm.
Figure 9.
Measured conductivity of a 9 cm × 9 cm sheet of graphene using THz-TDS (left), THz-CCS (center) and the relative error between them (right). The mean of the errors is , and the standard deviation is .
Figure 9.
Measured conductivity of a 9 cm × 9 cm sheet of graphene using THz-TDS (left), THz-CCS (center) and the relative error between them (right). The mean of the errors is , and the standard deviation is .
Figure 10.
measured using TDS against peak-to-peak measurements using CCS along with a linear fit (left) and the spatial relative error, when the peak-to-peak is used to predict (right).
Figure 10.
measured using TDS against peak-to-peak measurements using CCS along with a linear fit (left) and the spatial relative error, when the peak-to-peak is used to predict (right).
Figure 11.
Linear correlation between the peak-to-peak value and the DC conductivity in the range mS for a simulated 1 THz bandwidth pulse for different substrate thicknesses and refractive indices. The red cross denotes a parameter set close to the used experimental settings.
Figure 11.
Linear correlation between the peak-to-peak value and the DC conductivity in the range mS for a simulated 1 THz bandwidth pulse for different substrate thicknesses and refractive indices. The red cross denotes a parameter set close to the used experimental settings.
Table 1.
Summarization of key in situ properties used in the experiments for the THz-TDS and THz-CCS devices. Since different acquisition times are used, the performance of the devices cannot be compared directly.
Table 1.
Summarization of key in situ properties used in the experiments for the THz-TDS and THz-CCS devices. Since different acquisition times are used, the performance of the devices cannot be compared directly.
| Bandwidth | Dynamic Range | Delay | Spot Size | Acq. Time |
---|
TDS | 2.3 THz | 50 dB @ 500 GHz | 130 ps | 1.3 mm | 3 s |
CCS | 1.4 THz | 58 dB @ 100 GHz | 100 ps | 5.6 mm | 30 s |
| Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).