Journal of Applied Physics, 15 November 2008
J. Appl. Phys. 104, 103111 (2008) (7 pages)
©2008 American Institute of Physics. All rights reserved. Rightslink - Permissions for ReusePermissions for ReuseAbout Rightslink

Up: Issue Table of Contents
Go to: Previous Article | Next Article
Other formats: HTML (smaller files) | PDF ( kB)

Dependence on chirp rate and spectral resolution of the terahertz field pulse waveform measured by electro-optic detection using a chirped optical pulse and a spectrometer and its effect on terahertz spectroscopy

H. Murakami,1,2 *K. Shimizu,1,3 M. Katsurada,1,3 and S. Nashima3

1Kansai Photon Science Institute, Japan Atomic Energy Agency, Kizugawa, Kyoto 619-0215, Japan
2Research Unit for Quantum Beam Life Science Initiative, Japan Atomic Energy Agency, Kizugawa, Kyoto 619-0215, Japan
3Department of Applied Physics, Osaka City University, Sugimoto, Osaka 558-8585, Japan

(Received: 21 January 2008; accepted: 4 October 2008; published online: 20 November 2008)

Single-shot measurement of a terahertz field pulse waveform by electro-optic sampling using a chirped optical pulse and a spectrometer was demonstrated by and Jiang and Zhang [Appl. Phys. Lett. 72, 1945 (1998)]. We have performed an experimental and theoretical investigation into the dependence of the waveform thus measured on the chirp rate and spectral resolution. It was found that the waveform exhibits multicyclic behavior at a chirp rate of −0.24  THz2, which corresponds to a chirped-pulse width of over 10 ps, for the monocyclic original terahertz field, while it approaches the monocyclic behavior with decreasing pulse width. Further, broadening of the spectral resolution of the spectrometer gives rise to a monocyclic waveform in the chirp rate range where the waveform is expected to be multicyclic. In addition, we have derived an analytical expression for the terahertz field pulse waveform thus measured without using the method of stationary phase. The theoretical results were found to be consistent with measured ones. Finally, we examined the spectral bandwidth and resolution of terahertz spectroscopy using this method. ©2008 American Institute of Physics


Contents

INTRODUCTION

Optical constants such as complex permittivity have been measured for various materials since the development of terahertz time-domain spectroscopy (THz-TDS) using a femtosecond (fs) laser.1,2,3,4,5,6,7,8 This method allows simultaneous determination of both the real and imaginary parts of the optical constant across a spectral range from a few tens of gigahertz to several terahertz because it allows measurement of the change due to the optical constant in the amplitude and phase of a terahertz electric field.

THz-TDS is based on a pump-probe method. A fs optical pulse is first divided into two beams. One (pump) is used to generate a terahertz field pulse, which, for example, is emitted from a semiconductor such as InAs by irradiating its surface with the pump pulse. In measurements using the electro-optic (EO) effect, the other (probe) is overlapped with the terahertz field pulse within an EO crystal such as ZnTe, and the magnitude of the Pockels effect proportional to the amplitude of the terahertz field is obtained from measurement of the intensity of the probe pulse subject to the birefringence in the crossed-polarizer configuration. The temporal waveform of the terahertz field is measured by changing the time interval between the pump and probe pulses with an optical delay line. Thus, THz-TDS is not applicable to single-shot measurement of the temporal waveform of a terahertz field pulse.

A method that allows single-shot measurement of a terahertz field waveform was proposed by Jiang and Zhang.9 In the method, a fs probe pulse is linearly chirped and overlapped with a terahertz field pulse within the EO crystal. The chirped probe pulse is modulated by the terahertz field and dispersed onto a multichannel detector combined with a spectrometer. Since the wavelength axis can be converted to the time axis using the value of the chirp rate, the temporal waveform of the terahertz field is derived from the two spectra of the probe pulses with and without terahertz field modulation. They used a chirped pulse with a temporal width of approximately 30 ps and obtained a terahertz field pulse three times broader than the original pulse width measured by THz-TDS. Moreover, severely distorted terahertz field waveforms were observed even for chirp probe pulses with a time width of several picoseconds.10,11,12 Sun et al.13 analyzed the dependence of the temporal resolution of electro-optic detection using a chirped probe pulse (EODCP) on the chirp rate by assuming that the temporal widths of the chirped pulse and the terahertz field are much longer than the oscillation period of the optical beam and that the stationary phase method is applicable. The EODCP-derived terahertz field pulse waveform in their analysis was monocyclic for a monocyclic original terahertz field. It has recently been demonstrated that the temporal resolution of EODCP can be significantly enhanced by combination with an interferometric retrieval algorithm; the original terahertz field waveform is reproduced well by the algorithm.11,12

Single-shot measurement of the terahertz field pulse has been achieved by overlapping the terahertz radiation and the fs optical probe beam noncollinearly within an EO crystal, where the time window is determined by the spot size of the probe beam and the crossing angle between the two beams.14 Further, it has been demonstrated that for a time window of about 5 ps, the terahertz field waveform obtained from the single-shot measurement almost coincides with that from THz-TDS.

Complex systems such as liquids and proteins show low-frequency vibrations in the terahertz spectral range, and such vibrations are believed to play an important role in the dynamics of chemical reactions in solutions and proteins. Further, the origin of the boson peak observed in the low-frequency vibrational spectrum of the complex system is under debate.15,16 Thus, the low-frequency dynamics of the complex system has been extensively investigated by means of a variety of experimental methods such as due to inelastic neutron17,18,19 and light scatterings,20 as well as by THz-TDS.2,3,4,5,6,7,8 Most of the measurements, however, have been made for systems in thermal equilibrium. It is considered necessary to study how such low-frequency vibrations change during a reaction in order to explore the mechanism of the reaction in solutions and proteins.

This can be achieved by optical pump-terahertz probe (OPTP) spectroscopy. In dye solutions, the solvation dynamics after optical excitation of the dye molecule have been probed by a terahertz field pulse, and it has been suggested that librational modes localized close to the solute occur owing to the abrupt change caused by the excitation in the dipole moment of the dye molecule.21 In liquid hexane, the transient photoconductivity and recombination dynamics of quasifree electrons generated by two-photon ionization induced by the irradiation with a fs UV pulse have been probed by the terahertz pulse.22

THz-TDS is mostly adopted for the terahertz probe in OPTP spectroscopy. In such a case, a fs laser pulse is first divided into two beams, and one is used for exciting the sample. The other is used for THz-TDS; that is, it is further divided into two beams: one for the generation of the probe terahertz field pulse and the other for sampling the terahertz field. Hence, two optical delay lines are employed among the three fs pulses thus divided. OPTP spectroscopy with THz-TDS has several drawbacks. An intense fs laser such as that generated from a chirped-pulse-amplification system is often used in pump-probe spectroscopy to excite as many molecules as possible in the sample. The shot-to-shot fluctuation of such an intense laser deforms the terahertz field waveform obtained from THz-TDS. Even if the shot-to-shot fluctuation is small, the terahertz field pulse shape is deformed if the system under study does not relax within the period of pulse repetition or undergoes an irreversible process. Further, the use of two optical delay lines makes measurement overly time consuming. In fact, in the above mentioned OPTP studies, the amplitude of the terahertz field pulse was measured only at the specified detection delay time at which the waveform had its peak while the pump-probe delay time was scanned, but its temporal waveform was not measured.21,22 These measurements do not provide the so-called time-resolved terahertz spectrum. If EODCP is employed for the terahertz probe in OPTP spectroscopy, these drawbacks can be overcome.23,24 However, no experimental investigation has been performed to determine how the terahertz field waveform measured by EODCP depends on the chirp rate and spectral resolution of the spectrometer, although the chirp rate dependence of the terahertz field waveform has been analyzed under limited conditions using the stationary phase method.13 Such dependence determines the parameters such as the spectral resolution and bandwidth of terahertz spectroscopy using EODCP.

Thus, in this study, we measured the terahertz field waveform by EODCP as a function of the chirp rate and the spectral resolution of the spectrometer, and derived an expression for the temporal shape of the terahertz field pulse obtained from EODCP without using the stationary phase method. The experimental results were compared with the theoretical predictions. Further, in this paper, we examine the experimental conditions such as the chirp rate and the spectral resolution of the spectrometer in the application of EODCP to terahertz spectroscopy.

THEORETICAL BACKGROUND

In this section, we derive an expression of the temporal shape of the terahertz field pulse obtained from EODCP. Denoting the electric field waveform of a chirped probe pulse by C(t), we give the temporal shape of the chirped pulse modulated by a terahertz field ETHz(t) as follows:

<i>M</i>(<i>t</i>) = <i>C</i>(<i>t</i>)[<i>b</i> + <i>k</i><i>E</i><sub>THz</sub>(<i>t</i> − <i>tau</i>)],1

where k is a modulation coefficient and tau is the delay time between the terahertz field pulse and the chirped pulse. The first term on the right-hand side of Eq. (1) appears because the probe pulse passes partly through the two crossed polarizers, e.g., owing to the inherent residual birefringence of the EO crystal or the finite extinction ratio of the polarizers, where b is the transmission coefficient. By comparing the terahertz field waveform obtained from EODCP with numerical results, Yellampalle et al.10 showed that Eq. (1) is valid for an EO crystal with sufficient residual birefringence, as in the present study using a 1-mm-thick ZnTe crystal. However, the equation needs to be modified if analyzer detuning or an additional waveplate is used in the case where there is no residual birefringence.10

The waveform of the chirped pulse is defined as C(t)=exp(−t2T<sub><i>c</i></sub><sup>−2</sup>iat2iomega0t), where Tc, 2a, and omega0 are the width, chirp rate, and central frequency of the pulse, respectively. We represent a monocyclic original terahertz field waveform by ETHz(t)=(t/DeltaT)exp(−t2/DeltaT2), which is a symmetrically bipolar function and represents a global feature of the terahertz field waveform generated by a fs laser. If the modulated probe pulse is measured with a spectrometer with a multichannel detector, the spectral intensity I(omega)on is expressed by

<i>I</i>(<i>omega</i>)<sub>on</sub> [proportional] [integral]<sub>−[infinity]</sub><sup>[infinity]</sup><i>g</i>(<i>omega</i> − <i>omega</i><sup>[prime]</sup>)|[integral]<sub>−[infinity]</sub><sup>[infinity]</sup><i>M</i>(<i>t</i>)exp(<i>i</i> <i>omega</i><sup>[prime]</sup><i>t</i>)<i>d</i><i>t</i>|<sup>2</sup><i>d</i> <i>omega</i><sup>[prime]</sup>,2

where g(omegaomega[prime]) is the spectral response function of the spectrometer; we set g(omegaomega[prime])=delta(omegaomega[prime]) in this derivation. It is possible to solve Eq. (2) by using Gaussian integrals and dealing with complex numbers in the polar form. Further, we follow the procedure of Jiang and Zhang9 to derive the terahertz field waveform from EODCP, namely,

<i>E</i><sub>mes</sub>(<i>omega</i>) = ((<i>I</i>(<i>omega</i>)<sub>on</sub> − <i>I</i>(<i>omega</i>)<sub>off</sub>)/(<i>I</i>(<i>omega</i>)<sub>off</sub>)),3

where I(omega)off is obtained by setting M(t)=bC(t) in Eq. (2). Finally we obtain an expression of the terahertz field waveform as

<i>E</i><sub>mes</sub>(<i>omega</i>) [proportional] −{((2 <i>tau</i>)/<i>T</i><sub><i>c</i></sub><sup>2</sup>)cos[<i>theta</i>(<i>Delta</i> <i>omega</i>)] − (<i>Delta</i> <i>omega</i> + 2<i>a</i> <i>tau</i>)sin[<i>theta</i>(<i>Delta</i> <i>omega</i>)]}exp[−((<i>Delta</i> <i>omega</i><sup>2</sup>)/4)((1/(<i>beta</i> + <i>a</i><sup>2</sup><i>beta</i><sup>−1</sup>))−(1/(<i>T</i><sub><i>c</i></sub><sup>2</sup><i>gamma</i>)))−((<i>a</i> <i>tau</i> <i>Delta</i> <i>omega</i>)/(<i>Delta</i> <i>T</i><sup>2</sup>(<i>beta</i><sup>2</sup> + <i>a</i><sup>2</sup>)))],4

with

<i>theta</i>(<i>Delta</i> <i>omega</i>) = ((<i>a</i> <i>Delta</i> <i>omega</i><sup>2</sup>)/4)((1/(<i>beta</i><sup>2</sup> + <i>a</i><sup>2</sup>))−(1/(<i>gamma</i>)))−((<i>Delta</i> <i>omega</i> <i>tau</i> <i>beta</i>)/(<i>Delta</i> <i>T</i><sup>2</sup>(<i>beta</i><sup>2</sup> + <i>a</i><sup>2</sup>)))−((<i>a</i> <i>tau</i><sup>2</sup>)/(<i>Delta</i> <i>T</i><sup>4</sup>(<i>beta</i><sup>2</sup> + <i>a</i><sup>2</sup>)))+(1/2)arctan(<i>a</i><i>T</i><sub><i>c</i></sub><sup>2</sup>)−(3/2)arctan(<i>a</i>/<i>beta</i>),

where Deltaomega=omega0omega, beta=T<sub><i>c</i></sub><sup>−2</sup>+DeltaT−2, and gamma=T<sub><i>c</i></sub><sup>−4</sup>+a2. Here the quadratic term of k is neglected on the assumption that the modulation coefficient is very small and kETHz<<b. We confirmed that the waveform obtained from Eq. (4) agrees with that from Eq. (3) through the numerical integration of I(omega)on and I(omega)off with the same parameter values.

EXPERIMENTAL PROCEDURE

A light source was a regeneratively amplified Ti:sapphire laser (Clark-MXR CPA-1000). It generated an optical pulse with a temporal width of less than 150 fs [full width at half maximum (FWHM)], a wavelength of 800 nm, a repetition rate of about 1 kHz, and a pulse energy of 600  µJ. The fs laser beam was divided into two beams: one was irradiated onto an InAs(100) wafer to generate the terahertz field pulse, and the other was chirped by passing it through a pair of gratings. We used two types of grating pairs, i.e., 1200 and 600 grooves/mm, to change the chirp rate.

A schematic diagram of the experimental setup for EODCP is illustrated in Fig. 1. We employed a double-beam configuration to obtain I(omega)on and I(omega)off simultaneously. This configuration reduces the effect of the shot-to-shot fluctuation of the laser on the derived terahertz waveform using EODCP. The chirped beam was first passed through a polarizer and then divided into two beams. One beam, the reference R(omega), was led into a fiber that transferred the light to a spectrometer (Acton Research 300i) with three automatically exchangeable gratings (150, 1200, and 1800 grooves/mm), and was detected with a charge-coupled device (CCD) image sensor. The other, the signal S(omega), was aligned to travel collinearly with the terahertz field pulse, modulated by the terahertz field within a 1-mm-thick ZnTe crystal, and was passed through an analyzer into the other fiber. The reference and signal beams were detected simultaneously on two different areas of the image sensor. S(omega)off and R(omega)off were measured simultaneously in the absence of the terahertz field, while R(omega)on and S(omega)on [=I(omega)on] were measured in the presence of it. I(omega)off was obtained by R(omega)onS(omega)off/R(omega)off. Thus, we could measure I(omega)off and I(omega)on simultaneously in the double-beam configuration and derive the terahertz field waveform using Eq. (3). In the present study, the data were acquired by averaging over a few hundred shots to achieve a good signal-to-noise ratio. Further, an optical delay line was employed in order to vary the time interval between the probe pulse and the terahertz field pulse within the EO crystal. This was necessary to obtain the value of the chirp rate experimentally. The terahertz field waveform thus obtained was shifted on the frequency axis by changing the time interval to derive the value of the chirp rate from the magnitude of the shift as a function of the time interval and convert the abscissa for the terahertz field waveform from frequency to time, as shown below.

Figure 1.

THz-TDS was performed by the conventional method; that is, a fs optical probe pulse was used instead of a chirped one, and the spectrometer with a CCD sensor was replaced with a photodiode. The time interval for sampling the terahertz field was about 100 fs.

EXPERIMENTAL RESULTS AND DISCUSSION

A.Dependence of the terahertz field waveform on the chirp rate

Figure 2(a) depicts the terahertz waveform obtained by EODCP at a chirp rate of −0.24  THz2. The waveform is shifted along the abscissa (wavelength) when the delay time between the chirped probe and terahertz field pulses is changed. Figure 2(b) was obtained by tracking corresponding peak positions [arrow in Fig. 2(a)] in the terahertz waveforms at different delay times. We derived the value of the chirp rate from the slope of the line fitted to the data points. It was found that this value agrees with that obtained by a second-harmonic-generation frequency-resolved optical grating measurement and the estimated value from the distance between the two gratings and the number of grooves ruled on the grating per unit length. The upper axis in Fig. 2(a) is the time axis obtained from the wavelength axis using the chirp rate value. The temporal width of the chirped pulse at this chirp rate was determined to be about 15 ps (FWHM) by cross-correlation measurement between the original fs pulse and the chirped one. The terahertz waveform (black line) is compared with that obtained by THz-TDS (thick gray one) in Fig. 2(c). The waveform due to THz-TDS is monocyclic, although small long-lived oscillation due to absorption by water vapor in the terahertz frequency range is involved; all the measurements in the present study were made in the atmosphere. The waveform measured by EODCP is multicyclic. Further, the amplitude of the oscillation in the waveform derived by EODCP increases gradually with time and then decreases, as seen in Fig. 2(c), except for the central part around 2.5 ps. Such multicyclic behavior was not observed in the terahertz field waveform obtained from EODCP by Jiang and Zhang,9 although the temporal width of the chirped probe pulse in their measurement is comparable to that in this case. However, there are a few experimental studies that report distorted waveforms similar to the multicyclic behavior obtained here from EODCP using a several picosecond wide chirped pulse, although the effect of the distortion is not explicitly mentioned.10,11,12 Further, it should be noted that the temporal period of the multicycle in the terahertz waveform obtained by EODCP is roughly the same as that of the monocycle obtained by THz-TDS [Fig. 2(c)], while the former was three times longer than the latter in the measurement by Jiang and Zhang.9 This implies that the temporal resolution of EODCP is comparable to that of THz-TDS in our case. The experimental results at a chirp rate of −1.3  THz2 for different delay times are shown in Fig. 3. It is obvious that the number of oscillations is small compared with that at −0.24  THz2.

Figure 2. Figure 3.

B.Comparison between the experimental and calculated results

The terahertz field waveforms calculated from Eq. (4) for delay times of −10, 0, and 10 ps are illustrated in Fig. 4, together with the original terahertz field pulse (inset) used for the calculation, where Tc=22  ps, a=−0.12  THz2, DeltaT=0.45  ps, and omega0=375  THz (=800  nm). The values of DeltaT and omega0 are fixed for the calculations below. These parameter values are about the same as those used in the measurement for Fig. 2(a). However, the original terahertz pulse is not symmetrically bipolar, with long-lived oscillations due to the absorption by water vapor, as seen in Fig. 2(c), and the temporal envelope of the chirped pulse deviates from a Gaussian curve in the measurement. The waveform thus calculated reproduces the characteristic features observed in Fig. 2(a). Multicyclic behavior is obviously seen in Fig. 4. In both the experimental and theoretical results, the waveform consists of two damping oscillations moving in opposite directions on the abscissa except for the central part. The chirped probe pulse is modulated by the terahertz field pulse around the peak of the spectrum of the probe pulse at a delay time of 0 ps. The two damping oscillations are inversion symmetric about a wavelength of 800 nm at 0 ps in Fig. 4, while they are asymmetric at ±10  ps. That is, the amplitude at the start point of the damping oscillation on the longer-wavelength side is larger than that on the shorter-wavelength side at 10 ps, while it is smaller at −10  ps. This feature is remarkably enhanced in the measured waveform shown in Fig. 2(a). The waveforms are calculated at [a(THz2),Tc(ps)]=(−0.12,22), (−0.65,4), and (−1.5,1.8), with tau=0  ps in Fig. 5. The waveform is found to change from monocyclic to multicyclic behavior with decreasing magnitude of chirp rate. This tendency agrees with that seen in the experimental results. It should be noted that the multicyclic behavior does not emerge in the waveform derived by Sun et al.10 on the assumption that the stationary phase method is valid. Hence, analysis based on the stationary phase method is not applicable to the case where the terahertz field waveform derived by EODCP is multicyclic.

Figure 4. Figure 5.

C.Dependence on spectral resolution of the spectrometer

Figure 6(a) shows the dependence of the terahertz field waveform obtained by EODCP on the spectral resolution of the spectrometer at a chirp rate of −1.9  THz2. The waveform does not exhibit monocyclic behavior at the high spectral resolution (0.2 nm, FWHM), while it becomes monocyclic at the low resolution (3 nm). According to the fact that the temporal resolution is low at the low spectral resolution in Fig. 6(a), the bandwidth of the spectrum obtained from the Fourier transform of the temporal waveform is narrower at a spectral resolution of 3 nm than at 0.2 nm in Fig. 6(b). Further, we calculated the terahertz field waveform from Eq. (3) through numerical integration of Eq. (2) at a chirp rate of −0.24  THz2 in Fig. 7, where the delta-function and Gaussian function with widths of 0.2, 0.6, and 3.4 nm were used for g(omegaomega[prime]) in Eq. (2). It is found that the multicyclic behavior observed at the high spectral resolutions disappears completely at a resolution of 3.4 nm. Thus, the spectral resolution of the spectrometer plays an important role in determining the temporal shape of the terahertz field obtained by EODCP.

Figure 6. Figure 7.

The terahertz field waveform measured by Jiang and Zhang9 does not show multicyclic features despite the fact that the chirp rate and temporal width of the chirped probe pulse are comparable to those in our case at −0.24  THz2. This may be because the spectral resolution of the spectrometer used by them was not high enough to resolve the multicyclic waveform, as seen in Fig. 6(b), although its value is not described in the literature.9

D.Spectral bandwidth and resolution of terahertz spectroscopy using EODCP

The spectral resolution and bandwidth of THz-TDS is determined by the whole time window of observation and the temporal profile of the terahertz field pulse, respectively, according to time-frequency relation based on the Fourier transform. On the other hand, the spectral bandwidth of terahertz spectroscopy using EODCP is determined not only by the time profile of the original terahertz field pulse but also by the temporal resolution of the measurement.

The temporal resolution due to the spectral resolution (Deltaomega) of the spectrometer is given by Deltaomega/(2|a|) for EODCP; this value must ideally be smaller than the temporal width of the original terahertz waveform. For example, in the case that the spectral resolution and chirp rate are 0.2 nm and −0.24  THz2, respectively, the temporal resolution is about 0.4 ps. The terahertz waveform calculated at a 0.2-nm resolution is found to coincide fairly well with that at the infinite resolution in Fig. 7. On the other hand, when the temporal resolution is about 1.2 ps at a 0.6-nm resolution, the terahertz waveform broadens compared with that at the infinite resolution in Fig. 7. Further, it is found from Fig. 6 that the terahertz field waveform measured broadens at a 3-nm resolution compared with that at 0.2 nm, and hence the bandwidth of the frequency spectrum is narrower than that at 0.2 nm.

The temporal period of the multicycle in the terahertz field waveform measured by EODCP at a chirp rate of −0.24  THz2 and a spectral resolution of 0.2 nm is roughly the same as that of the monocycle obtained by THz-TDS [Fig. 2(c)]. The frequency spectrum obtained by EODCP (Fig. 8), accordingly, has a spectral bandwidth (about 2 THz) comparable to that obtained by THz-TDS [Fig. 6(b)], although the spectral shapes are not the same.

Figure 8.

As for the time window of the observation, a range of a few tens of picoseconds is required to obtain spectral resolution as small as 1  cm−1 (30 GHz). EODCP with a small chirp rate leads to a narrow time window of observation under the limited spectral width of the probe pulse (~10  nm in the present study), and hence a high spectral resolution is not achieved. In Fig. 6(b), the frequency spectrum of the terahertz field pulse obtained by THz-TDS shows sharp gaps due to absorption by water vapor, e.g., around 1.2 THz. However that obtained by EODCP at −1.9  THz2 (chirp rate) and 0.2 nm (spectral resolution) does not because of its narrow time window (~5  ps), i.e., low spectral resolution (~0.15  THz). On the other hand, there exist sharp gaps due to absorption by water vapor in the frequency spectrum obtained at −0.24  THz2 (Fig. 8). Therefore, EODCP with a small chirp rate has a crucial drawback for use in terahertz spectroscopy, even if the terahertz field waveform approaches monocyclic behavior. Thus, we must use a chirped probe pulse with a chirp rate that assures a time window for observation of a few tens of picoseconds, although we will possibly encounter multicyclic behavior in the terahertz waveform in such a case.

The multicyclic behavior in the terahertz field waveform calculated from Eq. (4) at a chirp rate of −0.24  THz2 leads to the oscillating behavior in the frequency spectrum shown in Fig. 8, although the envelope connecting the peaks in the spectrum coincides roughly with the spectrum obtained from the Fourier transform of the original terahertz field waveform (inset of Fig. 4). We consider the effect of the oscillating behavior on the derivation of the optical constants of materials by terahertz spectroscopy using EODCP. The quadratic term of the modulation coefficient is neglected in the derivation of Eq. (4), which implies that the EODCP instrument is assumed to exhibit a linear response to the terahertz field for the small modulations considered in this work. On the assumption of such a linear response, the waveform of the terahertz field obtained by EODCP is expressed by

<i>E</i><sub>ref</sub>(<i>t</i>) = [integral]<sub>0</sub><sup>[infinity]</sup><i>phi</i>(<i>t</i> − <i>t</i><sub>1</sub>)<i>E</i><sub>0</sub>(<i>t</i><sub>1</sub>)<i>d</i><i>t</i><sub>1</sub>,5

where phi(t) and E0(t) are the time profiles of the instrumental response and the original terahertz field, respectively. If the terahertz field is transmitted through the material under study and measured by EODCP, its waveform is given by

<i>E</i><sub>sig</sub>(<i>t</i>) = [integral][integral] <i>phi</i>(<i>t</i> − <i>t</i><sub>1</sub>)<i>eta</i>(<i>t</i><sub>1</sub> − <i>t</i><sub>2</sub>)<i>E</i><sub>0</sub>(<i>t</i><sub>2</sub>)<i>d</i><i>t</i><sub>1</sub><i>d</i><i>t</i><sub>2</sub>,6

where the material is assumed to exhibit a linear response to the terahertz field and the response function is denoted as eta(t). Thus, the complex transmission coefficient <i>eta</i>-tilde(omega) of the material is derived by

<i>eta</i>-tilde(<i>omega</i>) = ((<i>E</i>-tilde<sub>sig</sub>(<i>omega</i>))/(<i>E</i>-tilde<sub>ref</sub>(<i>omega</i>))),7

with

<i>E</i>-tilde<sub>sig</sub>(<i>omega</i>) = <i>phi</i>-tilde(<i>omega</i>)<i>eta</i>-tilde(<i>omega</i>)<i>E</i>-tilde<sub>0</sub>(<i>omega</i>)  and  <i>E</i>-tilde<sub>ref</sub>(<i>omega</i>) = <i>phi</i>-tilde(<i>omega</i>)<i>E</i>-tilde<sub>0</sub>(<i>omega</i>),8

where the terms with superscripts, e.g.,<i>E</i>-tildesig(omega), are the Fourier transforms of the corresponding terms without superscripts, e.g., Esig(t). The convolution theorem is used in the derivation of Eq. (8). The Fourier-transformed spectra of the instrument response function and the original terahertz waveform are canceled out in Eq. (7), and so have no effect on derivation of the optical constant of the material. The complex refractive index is derived from the complex transmission coefficient.25

The terahertz field waveform obtained by EODCP depends on the delay time between the terahertz field and chirped probe pulses, as shown in Figs. 2(a)3. This implies that the instrument response of EODCP varies with the delay time. If the refractive index of the material under study is not unity in the terahertz spectral range, the terahertz field pulse transmitted through the material is delayed from that passing through free space by an amount determined by the refractive index. Thus, the two terahertz field pulses temporally modulate different portions of the chirped probe pulse. Therefore, the instrument response is different for the two terahertz fields. A method for removing the effect of variable instrument response is to make the two terahertz field pulses overlap with the same (temporal) portion of the chirped probe pulse using an optical delay line, although the time difference between the two terahertz field pulses must be measured in advance. The chirped probe beam can be modulated by the terahertz field at any temporal portion using an optical delay line, as seen in Figs. 2(a)3. Further investigation on how to correct for the change in the instrument response will be reported elsewhere. Thus, one can determine the optical constant of a material by terahertz spectroscopy using EODCP with spectral resolution and bandwidth values which are comparable to those obtained by THz-TDS if appropriate values are adopted for the chirp rate and the spectral resolution of the spectrometer and if the instrument response <i>phi</i>-tilde(omega) is canceled out through the calculation of the optical constant in EODCP.

CONCLUSION

We have measured the temporal waveform of a terahertz field pulse using EO detection with a chirped optical probe beam and a spectrometer as a function of the chirp rate and spectral resolution of the spectrometer. A double-beam configuration was employed in this measurement to simultaneously measure the probe pulses with and without terahertz field modulation. The waveform thus measured at a chirp rate of −0.24  THz2 exhibited multicyclic behavior and approached monocyclic behavior with increasing magnitude of the chirp rate, although the original terahertz field waveform was monocyclic. Further, it was found that the broadening of the spectral resolution of the spectrometer leads to a change from a multicyclic waveform to a monocyclic waveform. Furthermore, an expression for the terahertz field waveform measured using the chirped probe pulse was derived without using the stationary phase method. The chirp rate and spectral resolution dependences of the terahertz field waveform measured were in agreement with those calculated. Moreover, it was demonstrated that this technique will be applicable to the terahertz spectroscopy if appropriate values are adopted for the chirp rate and the spectral resolution of the spectrometer and if the instrument response is canceled out through the calculation of the complex transmission coefficient.

ACKNOWLEDGMENTS

We would like to thank Professor M. Hosoda, Dr. A. Yokoyama, and Mr. Y. Toyota for the fruitful discussions.

REFERENCES


  1. S. Nashima, O. Morikawa, K. Takata, and M. Hangyo, J. Appl. Phys. 90, 837 (2001). [ISI] first citation in article
  2. S. R. Keiding, J. Phys. Chem. A 101, 5250 (1997). [Inspec] [ISI] first citation in article
  3. B. N. Flanders, R. A. Cheville, D. Grischkowsky, and N. F. Scherer, J. Phys. Chem. 100, 11824 (1996). [Inspec] [ISI] first citation in article
  4. J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996). [Inspec] [ISI] first citation in article
  5. C. Rønne, K. Jensby, B. J. Loughnane, J. Fourkas, O. F. Nielsen, and S. R. Keiding, J. Chem. Phys. 113, 3749 (2000). [ISI] first citation in article
  6. S. E. Whitmire, D. Wolpert, A. G. Markelz, J. R. Hillebrecht, J. Galan, and R. R. Birge, Biophys. J. 85, 1269 (2003). [ISI] [MEDLINE] first citation in article
  7. A. G. Markelz, A. Roitberg, and E. J. Heilweil, Chem. Phys. Lett. 320, 42 (2000). first citation in article
  8. M. Brucherseifer, M. Nagel, P. H. Bolivar, and H. Kurz, Appl. Phys. Lett. 77, 4049 (2000). [ISI] first citation in article
  9. Z. Jiang and X. -C. Zhang, Appl. Phys. Lett. 72, 1945 (1998). [ISI] first citation in article
  10. B. Yellampalle, K. Y. Kim, G. Rodriguez, J. H. Glownia, and A. J. Taylor, Opt. Express 15, 1376 (2007). [MEDLINE] first citation in article
  11. K. Y. Kim, B. Yellampalle, G. Rodriguez, R. D. Averitt, A. J. Taylor, and J. H. Glownia, Appl. Phys. Lett. 88, 041123 (2006). first citation in article
  12. B. Yellampalle, K. Y. Kim, G. Rodriguez, J. H. Glownia, and A. J. Taylor, Appl. Phys. Lett. 87, 211109 (2005). first citation in article
  13. F. G. Sun, Z. Jiang, and X. -C. Zhang, Appl. Phys. Lett. 73, 2233 (1998). first citation in article
  14. J. Shan, A. S. Weling, E. Knoesel, L. Bartels, M. Bonn, A. Nahata, G. A. Reider, and T. F. Heinz, Opt. Lett. 25, 426 (2000). [ISI] [MEDLINE] first citation in article
  15. E. Donth, The Glass Transition: Relaxation Dynamics in Liquids and Disordered Materials (Springer, Berlin, 2001). first citation in article
  16. Relaxations in Complex Systems 5: Proceedings of the Fifth International Discussion Meetings on Relaxations in Complex Systems, edited by F. Affouard, M. Descamps, and K. L. Ngai (Elsevier, New York, 2006). first citation in article
  17. W. Doster, S. Cusack, and W. Petry, Nature (London) 337, 754 (1989). [Inspec] [ISI] [MEDLINE] first citation in article
  18. M. Kataoka, H. Kamikubo, J. Yunoki, F. Tokunaga, T. Kanaya, Y. Izumi, and K. Shibata, J. Phys. Chem. Solids 60, 1285 (1999). [Inspec] [ISI] first citation in article
  19. H. Nakagawa, M. Kataoka, Y. Joti, A. Kitao, K. Shibata, A. Tokuhisa, I. Tsukushi, and N. Go, Physica B 385–386, 871 (2006). [Inspec] first citation in article
  20. S. V. Adichtchev, N. V. Surovtsev, J. Wiedersich, A. Brodin, V. N. Novikov, and E. A. Rössler, J. Non-Cryst. Solids 353, 1491 (2007) (ScienceDirect). [Inspec] first citation in article
  21. R. Mcelroy and K. Wynne, Phys. Rev. Lett. 79, 3078 (1997). [ISI] first citation in article
  22. E. Knoesel, M. Bonn, J. Shan, F. Wang, and T. F. Heinz, J. Chem. Phys. 121, 394 (2004). [MEDLINE] first citation in article
  23. K. Y. Kim, B. Yellampalle, J. H. Glownia, A. J. Taylor, and G. Rodriguez, Phys. Rev. Lett. 100, 135002 (2008). [MEDLINE] first citation in article
  24. K. Y. Kim, B. Yellampalle, G. Rodriguez, and J. H. Glownia, Bull. Am. Phys. Soc. 50, 117 (2005). first citation in article
  25. L. Duvillaret, F. Garet, and J. L. Coutaz, IEEE J. Sel. Top. Quantum Electron. 2, 739 (1996). first citation in article

FIGURES


Full figure (9 kB)

Fig. 1. The schematic diagram of the experimental setup for the EODCP of the terahertz field pulse with a spectrometer equipped with a CCD image sensor. A double-beam configuration was employed to obtain the probe pulses with modulation due to the terahertz field and without it simultaneously. ITO is indium tin oxide. First citation in article


Full figure (24 kB)

Fig. 2. (a) The terahertz field waveforms obtained by EODCP at five delay times between the input terahertz field and chirped probe pulses, at a chirp rate of −0.24  THz2. The ordinate of each waveform is shifted for the sake of clarity. (b) Corresponding peak positions [arrow in (a)] in the terahertz waveform, represented by frequency, plotted as a function of the delay time (open squares). The chirp rate is obtained from the slope of the solid line fitted to the data. The upper abscissa is derived from the lower one using the chirp rate value in (a). (c) The waveform obtained by EODCP (black line) at −0.24  THz2 is compared with that obtained by THz-TDS (thick gray line). First citation in article


Full figure (13 kB)

Fig. 3. Terahertz field waveforms obtained by EODCP at five delay times, at a chirp rate of −1.3  THz2. The ordinate of each waveform is shifted for the sake of clarity. The upper abscissa is derived from the lower one using the chirp rate value. First citation in article


Full figure (20 kB)

Fig. 4. The terahertz field waveforms calculated from Eq. (4) at three delay times and a chirp rate of −0.24  THz2, where Tc=22  ps, DeltaT=0.45  ps, and omega0=375  THz. The ordinate of each waveform is shifted for the sake of clarity. The original terahertz field waveform, ETHz(t) in Eq. (1), is depicted in the inset. The values of DeltaT and omega0 are fixed for the calculations below. First citation in article


Full figure (9 kB)

Fig. 5. Chirp rate dependence of the terahertz field waveform calculated from Eq. (4), where [a(THz2),Tc(ps)]=(−0.12,22), (−0.65,4), and (−1.5,1.8), with tau=0  ps. The chirp rate value is given by 2a. The ordinate of each waveform is shifted for the sake of clarity. First citation in article


Full figure (15 kB)

Fig. 6. (a) Terahertz field waveforms measured by EODCP at two spectral resolutions of the spectrometer [0.2 nm (black line) and 3 nm (thick gray line)] at a chirp rate of −1.9  THz2. (b) Frequency spectra obtained from the Fourier transforms of the two terahertz field waveforms, together with that obtained by THz-TDS. A logarithmic scale is used for the ordinate. First citation in article


Full figure (9 kB)

Fig. 7. Terahertz field waveform calculated from Eq. (3) through numerical integration of Eq. (2) at a chirp rate of −0.24  THz2, where the delta-function and Gaussian functions with widths of 0.2, 0.6, and 3.4 nm are used for g(omegaomega[prime]). The ordinate of each waveform is shifted for the sake of clarity. First citation in article


Full figure (9 kB)

Fig. 8. Frequency spectra obtained from the Fourier transform of the terahertz field waveform calculated from Eq. (4) at a chirp rate of −0.24  THz2 (black line) and the original waveform (thick gray line). The frequency spectrum obtained by EODCP (experiment) at −0.24  THz2 is plotted, where the ordinate is shifted for the sake of clarity. A logarithmic scale is used for the ordinate. First citation in article

FOOTNOTES

*Electronic mail: murakami.hiroshi@jaea.go.jp.


Up: Issue Table of Contents
Go to: Previous Article | Next Article
Other formats: HTML (smaller files) | PDF ( kB)