## Abstract

We report on the development of a terahertz time-domain technique for measuring the momentum relaxation time of charge carriers in ultrathin semiconductor layers. Making use of the Drude model, our phase sensitive modulation technique directly provides the relaxation time. Time-resolved THz experiments were performed on n-doped GaAs and show precise agreement with data obtained by electrical characterization. The technique is well suited for studying novel materials where parameters such as the charge carriers’ effective mass or the carrier density are not known a priori.

©2009 Optical Society of America

## 1. Introduction

Terahertz (THz) time-domain spectroscopy (TDS) has developed towards an alternative technique for characterizing electrical transport in condensed matter. Many studies on semiconductors demonstrate excellent agreement with the Drude model and characteristic properties such as the complex conductivity have been deduced [1–5]. Recent reports on organic semiconductors [6], low doped inorganic semiconductors [7], and nanoparticles [8] reported deviations from a pure Drude behavior. Experimental data were reproduced with excellent precision using modified Drude models such as the Drude-Smith model [9], the Cole-Davidson model [10], or by incorporating the Maxwell Garnett theory [11].

All Drude models have in common that scattering is described by the momentum relaxation time *τ* which makes this time constant one of the most fundamental quantities that control charge transport. Deducing *τ* from optical data, however, demands knowledge of further parameters, such as of the carriers' density *N*, the effective mass *m**, the background permittivity *ε*
_{∞}, and the thickness *d* of the sample [1,12,13]. The thickness *d* can be deduced in time-resolved experiments when the optical thickness exceeds the wavelength of the radiation [14–16]. However, the progress in material science provides many novel materials, which are extremely thin and have unknown charge carrier densities and effective masses. Besides THz TDS, far-infrared as well as microwave techniques are used for characterizing such modern materials [17,18]. But these approaches require precise knowledge of material parameters or use fitting procedures.

In this work, we present an optical method for probing the Drude relaxation time without using any further parameters. In our technique we modulate the optical properties of the sample and record the differential response signal at THz frequencies. We will show that the phase of this signal directly provides the relaxation time. It is worth pointing out, that our method is not restricted to the THz band. The method can be adapted to any spectroscopic technique, which provides the optical phase of transmitted or reflected radiation.

Small modulations of the sample's optical properties can be achieved most conveniently by electromodulation techniques similar to the method described by Allen et al. [19]. The fact that the modulation of the charge carrier density resembles the operation of a field-effect transistor illustrates the application potential of the technique. Currently, numerous novel materials such as nanostructures or organic semiconductors have been developed for similar devices. Many of these materials, however, are inhomogeneous, include a high density of traps and grain boundaries, or have a broadened density of electronic states. Dispersive transport and hopping dominate electrical properties and mask the fundamental carrier dynamics on microscopic scales [20,21]. In such cases THz modulation spectroscopy may provide further insight into the fundamental carrier relaxation dynamics.

## 2. Method

The concept of our method is illustrated in Fig. 1a
). We study a thin sheet of a semiconductor, which has the thickness *d* and contains mobile carriers. Within the part *Δd* the charge carrier density can be modulated, e.g. the density can be reduced to zero. A fraction of an incident THz electric field *E*_{inc} is transmitted through the doped layer, and generates a current density ** j** which is proportional to the complex conductivity

*σ*. In the following, it is convenient to discuss differential quantities such as the differential field Δ

*E*=*E*_{1}

**-**

*E*_{0}where

*E*_{1}and

*E*_{0}are the transmitted electric fields in the case of an interaction across the entire layer and in the case the layer thickness is reduced by Δ

*d*, respectively.

The differential current density Δ** j** can be described by a change of the thickness

**on the reflection and transmission of radiation, we consider that Δ**

*j***causes a differential magnetic field Δ**

*j***, which can be calculated by Ampere's law:**

*H*In layers which are thin compared with the wavelength ($d\ll \lambda $), the x-components of Δ** H** dominate the above equation whereas the

*z*-components can be neglected. We also consider that Δ

*d*is smaller than the skin depth. In consequence, the modulated current density Δ

**does not change significantly across the depleted layer. The magnetic field Δ**

*j***causes a differential electrical field given by Δ**

*H***=**

*E**Z*

_{0}Δ

*H*/*n*, where

*Z*

_{0}is the impedance of vacuum and

*n*is the refractive index of the surrounding medium. Figure 1b) displays the fields in the complex plane. Using

*E*

_{1y}as reference provides:

At this point we assume the conductivity of the original Drude model:

*N*is the carrier density and

*m** is the carriers’ effective mass [22,23]. Relating the imaginary part of

*S*to the real part yields the optical phase of the signal:

Thus, the momentum relaxation time can be directly deduced from *φ* without the use of further parameters. This includes the change of the layer thickness Δ*d*. Alternatively, the amplitude or intensity of the signal *S* can be analyzed in case that *ωτ* ≈1 [24]. In the Hagen-Rubens range (ωτ$\ll 1$), however, the spectra become structureless because the amplitude is proportional to 1/(1 + *ω*
^{2}
*τ*
^{2})^{1/2}. In contrast, the tangent of the phase angle is proportional to *τ*, which illustrates the strength of our technique when characterizing materials in which ultrafast scattering occurs.

## 3. Experimental

We tested the method on a GaAs structure grown by molecular beam epitaxy (MBE). It comprises an n-doped layer grown on a semi-insulating GaAs substrate. The layer is doped at 2·10^{16} cm^{−3} and has a thickness of 2 µm. For electronic control we fabricated a Schottky contact on top of the structure by depositing a 10 nm thick film of Cr. The second terminal is an ohmic AuGe contact. The doped layer in the GaAs can be partially depleted from electrons by the application of a negative bias *V* to the Schottky contact. The width of the depletion region follows a square root dependence on *V* [25], from which the modulated depth Δ*d* can be calculated. In our experiments Δ*d* is of the order of 100 nm. This is much smaller than the skin depth of the electron gas, which is about 20 μm at THz frequencies. Thus, the assumptions made in section 2 are justified. They can be applied to many other devices where conductive sheets are modulated, as for instance, in field-effect transistors.

In our experiments we apply time-resolved THz transmission spectroscopy. The sample is mounted in a cryostat, which allows for low-temperature measurements. A titanium-sapphire laser delivers pulses of 80 fs duration at 780 nm center wavelength and at a repetition rate of 80 MHz. About 700 mW of the laser power is used for the photoexcitation of an interdigitated THz emitter [26,27]. Parabolic mirror optics with NA = 0.3 focus the THz radiation onto the GaAs structure. The transmitted radiation is time-resolved by standard electro-optic sampling in a [110] ZnTe crystal of 1 mm thickness [28]. The setup has a bandwidth of 2.8 THz and the signal to noise ratio is about 10^{5} Hz^{1/2}. We measured the field strength of the THz pulses using the procedure described in Ref [29]. and deduced a peak field of 70 V/cm. The curve in Fig. 2a
) shows a THz pulse transmitted through the GaAs structure, whereas Fig. 2b) displays the differential signal, which we obtained by modulating the structure with a bias of −5 V. This corresponds to a change of the depletion width of about 400 nm.

## 4. Results and discussion

In order to prove the concept of our technique we measured the differential transmission through the GaAs structure at room temperature. The electron gas was modulated by switching the bias to the Schottky contact between 0 V and −1 V. This corresponds to a modulated layer thickness of 120 nm or to a modulated sheet density of n_{2}
* _{D}* = 2.4· 10

^{11}cm

^{−2}. We recorded the THz signals

*E*

_{1y}and

*E*

_{0y}[30]. After Fourier transform of the differential THz signal Δ

*E*=

_{y}*E*

_{1y}-

*E*

_{0y}, the optical phase was obtained as described above. The data depicted in Fig. 3 show the expected dependence of the phase on the frequency. Using a relaxation time of 198 ± 2 fs provides an excellent agreement with Eq. (5). The linear behavior of tan

*φ*confirms the validity of the original Drude model in our case. We do not observe deviations from the Drude model as reported for poly(3-hexylthiophene) [6], low doped silicon [7], and nanoparticles [8].

In section 2 we assumed that only the lateral components of Δ** H** contribute to Ampere's law whereas the perpendicular components can be neglected. We also assumed that the exciting THz field is virtually constant across the modulated layer. We validate these assumptions by measuring the THz signals in dependence on the modulation depth Δ

*d*by increasing the modulation bias

*V*. As displayed in Fig. 4a ) the measured phase shows no dependence on

*V*, which proves the validity of the assumptions. In contrast, the amplitude plotted in Fig. 4b) follows a square root dependence as expected for the increase of the depletion width. The calculation, however, shows deviations for higher voltages. We attribute these differences to imperfections of the Schottky contact and to inhomogeneities of the doping profile. This illustrates the difficulties that arise in deducing transport properties from amplitude data and highlights the reliability of phase measurements.

In the case of our well-defined MBE grown structure, parameters such as the doping density and the effective carrier mass are well known. Thus in this particular case, the momentum relaxation time can be extracted from four-point measurements [31]. Figure 5 shows temperature resolved data. Electrical and phase measurement differ by about 3%, which illustrates that our technique provides precise data. At low temperatures, however, freeze-out of the donors may lead to a reduced electrical conductivity, which explains the differences at about 100 K.

We emphasize that our approach is not advantageous in the case of homogeneous semiconductors with well known properties. In this case electrical characterization certainly is less elaborate. But many modern materials exhibit unknown charge carrier densities and effective masses. Additionally, many new materials are inhomogeneous. In such systems the mesoscopic morphology may limit electrical transport and in consequence extremely low mobilities are observed by electrical characterization. Examples are polycrystalline semiconductors, percolated networks of carbon nanotubes, and many organic substances where hopping transport is the most limiting factor. As mentioned in the introduction, such systems can be well described by modified Drude models. For the technique described in this work, we so far assumed the original Drude conductivity. It is worth pointing out that similar results as Eq. (5) can be achieved for the Drude-Smith approach [9] and the Cole-Davidson model [10]. This indicates that the presented method of phase analysis can be applied in a much broader framework and may open up new opportunities for characterizing charge transport in modern materials.

## 5. Conclusions

We developed a technique for directly measuring the Drude momentum relaxation time *τ* of charge carriers using the optical phase of a differential THz signal. The modulation of the structure under investigation provides a phase signal, which is proportional to *τ*. In contrast to other spectroscopic approaches, our technique is advantageous particularly in the Hagen-Rubens range. Thus, we are confident that our technique will provide new insights into many modern materials where scattering occurs on ultrafast time scales and electronic properties are not directly accessible by electronic means.

## Acknowledgments

This work is partially supported by the Deutsche Forschungsgemeinschaft (DFG) through the Nanosystems Initiative Munich (NIM), and by the Deutsche Forschungsgemeinschaft (DFG), contract Ke 516/1-1. The authors acknowledge technical support by A. Guggenmos and S. Niedermaier.

## References and links

**1. **M. van Exter and D. Grischkowsky, “Optical and electronic properties of doped silicon from 0.1 to 2 THz,” Appl. Phys. Lett. **56**(17), 1694–1696 (1990). [CrossRef]

**2. **M. van Exter and D. Grischkowsky, “Carrier dynamics of electrons and holes in moderately doped silicon,” Phys. Rev. B **41**(17), 12140–12149 (1990). [CrossRef]

**3. **T.-I. Jeon, D. Grischkowsky, A. K. Mukherjee, and R. Menon, “Electrical characterization of conducting polypyrrole by THz time-domain spectroscopy,” Appl. Phys. Lett. **77**(16), 2452–2454 (2000). [CrossRef]

**4. **O. Ostroverkhova, D. G. Cooke, S. Shcherbyna, R. F. Egerton, F. A. Hegmann, and J. E. Anthony, “Bandlike transport in pentacene and functionalized pentacene thin films revealed by subpicosecond transient photoconductivity measurements,” Phys. Rev. B **71**(3), 035204 (2005). [CrossRef]

**5. **P. Parkinson, J. Lloyd-Hughes, Q. Gao, H. H. Tan, C. Jagadish, M. B. Johnston, and L. M. Herz, “Transient THz conductivity of GaAs nanowires,” Nano Lett. **7**(7), 2162–2165 (2007). [CrossRef]

**6. **X. Ai, M. C. Beard, K. P. Knutsen, S. E. Shaheen, G. Rumbles, and R. J. Ellingson, “Photoinduced charge carrier generation in a poly(3-hexylthiophene) and methanofullerene bulk heterojunction investigated by time-resolved terahertz spectroscopy,” J. Phys. Chem. B **110**(50), 25462–25471 (2006). [CrossRef] [PubMed]

**7. **T.-I. Jeon and D. Grischkowsky, “Observation of a Cole-Davidson type complex conductivity in the limit of very low carrier densities in doped silicon,” Appl. Phys. Lett. **72**(18), 2259–2261 (1998). [CrossRef]

**8. **E. Hendry, M. Koeberg, B. O’Regan, and M. Bonn, “Local field effects on electron transport in nanostructured TiO_{2} revealed by terahertz spectroscopy,” Nano Lett. **6**(4), 755–759 (2006). [CrossRef] [PubMed]

**9. **N. V. Smith, “Classical generalization of the Drude formula for the optical conductivity,” Phys. Rev. B **64**(15), 155106 (2001). [CrossRef]

**10. **D. W. Davidson and R. H. Cole, “Dielectric relaxation in glycerol, propylene glycol, and n-propanol,” J. Chem. Phys. **19**(12), 1484–1490 (1951). [CrossRef]

**11. **J. C. M. Garnett,, “Colours in metal glasses and in metallic films,” Philos. Trans. R. Soc. London Sect. A **203**(1), 385–420 (1904). [CrossRef]

**12. **M. Schall and P. U. Jepsen, “Photoexcited GaAs surfaces studied by transient terahertz time-domain spectroscopy,” Opt. Lett. **25**(1), 13–15 (2000). [CrossRef]

**13. **In fact, the THz response scales with the sample's plasma frequency given by carrier density, effective mass, and background permittivity.

**14. **L. Duvillaret, F. Garet, and J.-L. Coutaz, “Highly precise determination of optical constants and sample thickness in terahertz time-domain spectroscopy,” Appl. Opt. **38**(2), 409–415 (1999). [CrossRef]

**15. **T. D. Dorney, R. G. Baraniuk, and D. M. Mittleman, “Material parameter estimation with terahertz time-domain spectroscopy,” J. Opt. Soc. Am. A **18**(7), 1562–1570 (2001). [CrossRef]

**16. **M. Scheller, C. Jansen, and M. Koch, “Analyzing sub-100-μm samples with transmission terahertz time domain spectroscopy,” Opt. Commun. **282**(7), 1304–1306 (2009). [CrossRef]

**17. **M. Fischer, M. Dressel, B. Gompf, A. K. Tripathi, and J. Pflaum, “Infrared spectroscopy on the charge accumulation layer in rubrene single crystals,” Appl. Phys. Lett. **89**(18), 182103 (2006). [CrossRef]

**18. **M. P. de Haas, R. J. O. M. Hoofman, L. D. A. Siebbeles, and J. M. Warman,“Highly mobile electrons and holes on isolated chains of the semiconducting polymer poly(phenylene vinylene),” Nature **392**(6671), 54–56 (1998). [CrossRef]

**19. **S. J. Allen Jr, D. C. Tsui, and F. DeRosa, “Frequency dependence of the electron conductivity in the silicon inversion layer in the metallic and localized regimes,” Phys. Rev. Lett. **35**(20), 1359–1362 (1975). [CrossRef]

**20. **H. Scher and E. W. Montroll, “Anomalous transit-time dispersion in amorphous solids,” Phys. Rev. B **12**(6), 2455–2477 (1975). [CrossRef]

**21. **M. Pollak, “On dispersive transport by hopping and by trapping,” Philos. Mag. **36**(5), 1157–1169 (1977). [CrossRef]

**22. **P. Drude, “Zur Elektronentheorie der Metalle,” Annalen der Physik **306**(3), 566–613 (1900). [CrossRef]

**23. **N. Ashcroft, and D. Mermin, *Solid State Physics* (Harcourt, 1976).

**24. **M. Dressel, and G. Grüner, *Electrodynamics of Solids - Optical Properties of Electrons in Matter* (Cambridge University press, 2002).

**25. **S. Sze, *Semiconductor Devices* (John Wiley & Sons, 1985).

**26. **A. Dreyhaupt, S. Winnerl, T. Dekorsy, and M. Helm, “High-intensity terahertz radiation from a microstructured large-area photoconductor,” Appl. Phys. Lett. **86**(12), 121114 (2005). [CrossRef]

**27. **G. Acuna, F. Buersgens, C. Lang, M. Handloser, A. Guggenmos, and R. Kersting, “Interdigitated Terahertz Emitters,” Electron. Lett. **44**(3), 229–231 (2008). [CrossRef]

**28. **Q. Wu and X.-C. Zhang, “Ultrafast electro-optic field sensors,” Appl. Phys. Lett. **68**(12), 1604–1606 (1996). [CrossRef]

**29. **P. C. M. Planken, H.-K. Nienhuys, H. J. Bakker, and T. Wenckebach, “Measurement and calculation of the orientation dependence of terahertz pulse detection in ZnTe,” J. Opt. Soc. Am. B **18**(3), 313–317 (2001). [CrossRef]

**30. **To be more specific, *E*_{0y} and *E*_{1y} were measured after passage through the GaAs-wafer and the optical detection system. This does not affect differential quantities.

**31. **F. M. Smits, “Measurement of sheet resistivities with the four-point probe,” Bell Syst. Tech. J. **37**, 711–718 (1958).