Optimal chirped probe pulse length for terahertz pulse measurement

A detailed analysis of the relationship between the duration of the chirped probe pulse and the bipolar terahertz (THz) pulse length in the spectral encoding technique is carried out. We prove that there is an optimal chirped probe pulse length (or an optimal chirp rate of the chirped probe pulse) matched to the input THz pulse length and derive a rigorous relationship between them. We find that only under this restricted condition the THz signal can be correctly retrieved. ©2008 Optical Society of America OCIS codes: (040. 2235) Far infrared or terahertz; (110. 0110) Imaging systems References and links 1. S.P. Mickan and X.-C. Zhang, "T-ray sensing and imaging," Interrnational Journal of High Speed Electronics and Systems, 12, 601(2003). 2. Q. Wu and X.-C.Zhang, “7 terahertz broadband GaP electro-optic sensor,”Appl. Phys. Lett. 70, 1784-1786 (1997). 3. Zhiping Jiang and X.-C. Zhang, “Electro-optic measurement of THz field pulses with a chirped optical beam,” Appl. Phys. Lett. 72, 1945-1947 (1998). 4. Jie. Shan, Aniruddha S. Weling, Ernst Knoesel, Ludwig Bartels, Mischa Bonn, Ajay Nahata, Georg A. Reider, and Tony F. Heinz, “single-shot measurement of terahertz electromagnetic pulses by use of electrooptic sampling,”Opt. Lett. 25, 426-428 (2000). 5. Steven P. Jamison, Jingling Shen, A. M. Macleod and W. A. Gillespie, and D. A. Jaroszynski, “Hightemporal-resolution, single-shot characterization of terahertz pulses,”Opt. Lett. 28, 1710-1712 (2003). 6. B. Yellampalle, K. Y. Kim, G. Rodriguez, J. H. Glownina, and A. J.Taylor, Appl. Phys. Lett. 87, 211109 (2005). 7. K. Y. Kim, B. Yellampalle, G. Rodriguez, R. D. Averitt, A. J. Taylor, and J. H. Glownia, “Single-shot, interferometric, high-resolution, terahertz field diagnostic,” Appl. Phys. Lett. 88, 041123 (2006). 8. Jeroen van Tilborg, “Coherent terahertz radiation from laser-wakefield-accelerated electron beams”, Doctoral Dissertation, p.76 (2006). 9. F. G. Sun, Zhiping Jiang, and X.-C. Zhang, “Analysis of terahertz pulse measurement with a chirped probe beam,”Appl. Phys. Lett. 73, 2233-2235 (1998). 10. J. R. Fletcher, “Distortion and uncertainty in chirped pulse THz spectrometers,” Opt. Express 10, 14251430 (2002). 11. Yellampalle, K. Y. Kim, G. Rodriguez, J. H. Glownia, and A. J. Taylor, “Details of electro-optic terahertz detection with a chirped probe pulse,” Opt. Express 15, 1376-1383 (2007). For past two decades THz pulse or T-ray detection techniques have been developed due to the requirements to numerous applications of THz science including material characterization, biomedical imaging, and tomography. One of the typical T-ray detection techniques is electrooptic sampling (EOS) , a scanning technique. Through EOS, the electric field profile including phase and amplitude of a THz pulse can be measured. The main benefit of this technique is of its high temporal resolution depending on the short probe pulse length. However, the measured THz profiles represent averaged waveforms and a long time is required to obtain each profile. In addition, EOS is not suitable to be used in some cases such as the T-rays with low repetition rate or with strong shot-to-shot fluctuations and experiments #95893 $15.00 USD Received 13 May 2008; revised 2 Jul 2008; accepted 2 Jul 2008; published 1 Aug 2008 (C) 2008 OSA 4 August 2008 / Vol. 16, No. 16 / OPTICS EXPRESS 12342 that may have low duty cycles. To solve these problems three different single-shot detection techniques 3, 4, 5 were developed. Among these techniques are the popular and practical spectral encoding 3 and the cross-correlation 5 techniques which normally are used for different cases. Compared to each other, the spectral encoding technique is more convenient to use due to its simple optical arrangement, capability of measuring the THz signal in realtime and providing THz spatiotemporal imaging with its disadvantage of limited temporal resolution. The cross-correlation technique has higher resolution which only depends on the duration of the short probe pulse while it has some disadvantages such as its more complicated optical arrangement and its absence of the ability of spatiotemporal imaging. To enhance the temporal resolution of the spectral encoding technique, an interferometric retrieval algorithm was proposed 6 and was applied in experiments 7 recently. This technique can provide transform-limited temporal resolution which is mainly limited by the spectral bandwidth of the optical probe pulse, regardless of its chirp. However, with this technique a complicated matrix inversion equation needs to be resolved numerically to retrieve the THz waveform . Furthermore, this algorithm needs to be improved due to the poor signal-to-noise ratio and the dominant presence of algorithm artifacts such as the fast oscillations when it was applied in some experiments . The spectral encoding technique not only has been applied to many experiments but also has been analyzed theoretically by some authors 9, 10 after it was proposed. In their detailed analysis the temporal resolution Tmin of the detection system was obtained as Tmin≈ (T0Tc) 1/2 where T0 is the original probe pulse length and Tc is the chirped probe pulse length. It was proposed that the THz signal could be retrieved without distortion if the pulse length T of the detected THz field satisfies the condition T ≥ Tmin≈ (T0Tc) . Though their theory can explain some experimental results, we found that this condition is a relatively inaccurate description. According to Sun’s analysis , to reduce the distortion, Tc must be reduced so that the temporal resolution could be improved while the smallest Tc should be larger than the THz duration. Hence a compromise between the temporal resolution and the chirp rate is desirable. In other words, there would be an optimal chirp rate or an optimal chirped probe pulse length that can be applied to measure the THz signal without distortion if the THz pulse length is given. In this paper we present a detailed analysis and a rigorous deduction of the relationship between the chirped probe pulse length Tc and the THz pulse length T. We prove that there is an optimal chirped probe pulse length Tco (or an optimal chirp rate) that matches with the input THz pulse length T and only under this restricted condition the THz signal can be retrieved properly. For the sake of argument, we use the same definitions of the electric field of the original probe pulse, the chirped probe pulse, and the THz pulse as described in reference 9 in following deduction. The electric field component of the original probe pulse with a central frequency ω0 and an envelope Gaussian function can be written as ( ) 2 2 0 0 0 ( ) exp E t t T i t ω − = − − . Here T0 is the pulse length, which is related to the laser spectral bandwidth Δω0 through T0=2/Δω0. After stretched by a grating pair or a dispersive glass, the electric field component of the probe pulse can be written as ( ) 2 2 2 0 ( ) exp c c E t t T i t i t α ω − = − − − where Tc is the pulse length after chirping, and 2α is the chirp rate which is approximated as the laser bandwidth divided by the chirped laser pulse length, i.e. 0 2 / c T α ω ≈ Δ . Assuming a bipolar THz waveform as ( ) 1 2 2 ( ) exp THz E t tT t T − − = − − , here we define a characteristic time 2T which is the interval between the maximum and the minimum and define T as the THz pulse length for simplicity. When the chirped probe pulse and a THz pulse co-propagate in an electro-optic crystal, the chirped pulse is modulated by the pulsed THz field through the Pockels effect. Considering #95893 $15.00 USD Received 13 May 2008; revised 2 Jul 2008; accepted 2 Jul 2008; published 1 Aug 2008 (C) 2008 OSA 4 August 2008 / Vol. 16, No. 16 / OPTICS EXPRESS 12343 the case that the polarizer and the analyzer are crossed to each other, the electric field component of the chirped probe beam modulated by the THz field can be written as 3, : ( ) ( )[1 ( )] ( ) ( ) ( ) m c THz c THz c E t E t kE t E t kE t E t τ τ = + − = + − with ETHz(t) the electric field of the THz waveform, τ is the relative time delay between the probe pulse and the THz pulse, and k is the modulation constant. Here we assume τ = 0 for simplicity due to the fact that the probe pulse can be synchronized with the THz pulse by adjusting the time delay in principle. Thus we have ( ) ( ) ( ) ( ) ( ) ( ) m c THz c c s E t E t kE t E t E t E t = + = + . The spectral modulation is spatially separated on the CCD of the spectrometer. Assuming that the measured signal on a CCD pixel corresponds to the optical frequency ω1 and the spectral resolution function of the spectrometer is g(ω1-ω), the spectrum of the background chirped pulse Ic in case of no THz field modulation and the modulated spectrum Im by THz pulse can be expressed as the convolution of the spectral function and the square of the Fourier transform of Ec(t) and Em(t), respectively. They are: 2 1 ( ) ( ) c c I g E ω ω ω ∝ − ∗

For past two decades THz pulse or T-ray detection techniques have been developed due to the requirements to numerous applications 1 of THz science including material characterization, biomedical imaging, and tomography.One of the typical T-ray detection techniques is electrooptic sampling (EOS) 2 , a scanning technique.Through EOS, the electric field profile including phase and amplitude of a THz pulse can be measured.The main benefit of this technique is of its high temporal resolution depending on the short probe pulse length.However, the measured THz profiles represent averaged waveforms and a long time is required to obtain each profile.In addition, EOS is not suitable to be used in some cases such as the T-rays with low repetition rate or with strong shot-to-shot fluctuations and experiments that may have low duty cycles.To solve these problems three different single-shot detection techniques 3,4,5 were developed.Among these techniques are the popular and practical spectral encoding 3 and the cross-correlation 5 techniques which normally are used for different cases.Compared to each other, the spectral encoding technique is more convenient to use due to its simple optical arrangement, capability of measuring the THz signal in realtime and providing THz spatiotemporal imaging with its disadvantage of limited temporal resolution.The cross-correlation technique has higher resolution which only depends on the duration of the short probe pulse while it has some disadvantages such as its more complicated optical arrangement and its absence of the ability of spatiotemporal imaging.To enhance the temporal resolution of the spectral encoding technique, an interferometric retrieval algorithm was proposed 6 and was applied in experiments 7 recently.This technique can provide transform-limited temporal resolution which is mainly limited by the spectral bandwidth of the optical probe pulse, regardless of its chirp.However, with this technique a complicated matrix inversion equation needs to be resolved numerically to retrieve the THz waveform 6 .Furthermore, this algorithm needs to be improved due to the poor signal-to-noise ratio and the dominant presence of algorithm artifacts such as the fast oscillations when it was applied in some experiments 8 .
The spectral encoding technique not only has been applied to many experiments but also has been analyzed theoretically by some authors 9,10 after it was proposed.In their detailed analysis the temporal resolution T min of the detection system was obtained as where T 0 is the original probe pulse length and T c is the chirped probe pulse length.It was proposed that the THz signal could be retrieved without distortion if the pulse length T of the detected THz field satisfies the condition T ≥ T min ≈ (T 0 T c ) 1/2 .Though their theory can explain some experimental results, we found that this condition is a relatively inaccurate description.According to Sun's analysis 9 , to reduce the distortion, T c must be reduced so that the temporal resolution could be improved while the smallest T c should be larger than the THz duration.Hence a compromise between the temporal resolution and the chirp rate is desirable.In other words, there would be an optimal chirp rate or an optimal chirped probe pulse length that can be applied to measure the THz signal without distortion if the THz pulse length is given.In this paper we present a detailed analysis and a rigorous deduction of the relationship between the chirped probe pulse length T c and the THz pulse length T. We prove that there is an optimal chirped probe pulse length T co (or an optimal chirp rate) that matches with the input THz pulse length T and only under this restricted condition the THz signal can be retrieved properly.
For the sake of argument, we use the same definitions of the electric field of the original probe pulse, the chirped probe pulse, and the THz pulse as described in reference 9 in following deduction.The electric field component of the original probe pulse with a central frequency ω 0 and an envelope Gaussian function can be written as ( ) . Here T 0 is the pulse length, which is related to the laser spectral bandwidth Δω 0 through T 0 =2/Δω 0 .After stretched by a grating pair or a dispersive glass, the electric field component of the probe pulse can be written as ( ) where T c is the pulse length after chirping, and 2α is the chirp rate which is approximated as the laser bandwidth divided by the chirped laser pulse length, i.e.
. Assuming a bipolar THz waveform as ( ) , here we define a characteristic time 2T which is the interval between the maximum and the minimum and define T as the THz pulse length for simplicity.
When the chirped probe pulse and a THz pulse co-propagate in an electro-optic crystal, the chirped pulse is modulated by the pulsed THz field through the Pockels effect.Considering the case that the polarizer and the analyzer are crossed to each other, the electric field component of the chirped probe beam modulated by the THz field can be written as 3,11 : with E THz (t) the electric field of the THz waveform, τ is the relative time delay between the probe pulse and the THz pulse, and k is the modulation constant.Here we assume τ = 0 for simplicity due to the fact that the probe pulse can be synchronized with the THz pulse by adjusting the time delay in principle.Thus we have . The spectral modulation is spatially separated on the CCD of the spectrometer.Assuming that the measured signal on a CCD pixel corresponds to the optical frequency ω 1 and the spectral resolution function of the spectrometer is g(ω 1 -ω), the spectrum of the background chirped pulse I c in case of no THz field modulation and the modulated spectrum I m by THz pulse can be expressed as the convolution of the spectral function and the square of the Fourier transform of E c (t) and E m (t), respectively.They are: The difference between I m and I c then can be expressed as , where ) , respectively.
For a small amplitude of E THz (t), k<<1, |E s (ω 1 )| 2 is negligible and thus can be ignored.Then we have . Defining a function S (ω 1 ) which is related to the THz signal as follows: For a linear chirp, the probe instantaneous frequency ω 1 measured on the CCD is proportional to the time t′ .In other words, the THz field profile can be directly obtained through measurements of the spectral modulation , and introducing two dimensionless pulse lengths m=T/T 0 and n=T c /T 0 , one obtains Given T and T c , γ is a constant which does not affect the profile of the reconstructed THz field.Due to the fact that it is in the exponent which is related to t′ , χ is the main factor which affects the profile of the reconstructed THz field.Though the function 0 ( / ) f t T ′ has a plateau around the center of the time window and it oscillates at both sides of the time   One can see clearly from Eq. ( 6) that the function ( ) ( ') In other words, the real THz profile can be retrieved if χ =1 while severe distortion would be introduced if χ >1 or χ <1.When m>1.4 (this condition can be easily be fulfilled due to the fact that T >>T 0 ) then the only reasonable solution of equation χ =1 follows as .
Here n is the dimensionless optimal probe pulse length.Thus we obtain the optimal duration T co =nT 0 and the optimal chirp rate . This means that it is possible to retrieve the THz signal without distortion if one uses a chirped probe pulse with the optimal duration T co matched to a input THz pulse length T = mT 0 according to Eq. (7).The relationship between m and n according to Eq. ( 7) is shown in Fig. 2. One can see that the dimensionless optimal chirped probe pulse length n increases nonlinearly with the input dimensionless THz pulse length m.Fig. 2. Relationship (solid line) between the dimensionless optimal chirped pulse length n and the input dimensionless THz pulse length m.The restricted condition T c ≤T 2 /T 0 or T ≥(T 0 T c ) 1/2 described by Sun's theory is also shown as the area under the curve n=m 2 (or T c =T 2 /T 0 , dashed line) and the curve itself compared with the optimal T c curve.
By comparing between the condition described by Sun 9 and Fletcher 10 and the relationship of Eq. ( 7), one can find that their description is different from ours.In their theory,  6) under the conditions of input THz pulse length T = 0.5ps (m =20 when T 0 =25fs and m =50 when T 0 =10fs) and their corresponding optimal chirped pulse lengths T co =2.6437ps when T 0 =25fs, T co = 4.1943ps when T 0 =10fs, respectively, which are calculated from equation (7).The retrieved THz fields calculated according to the equation ( 6  the condition to retrieve the original THz field without distortion was described as T ≥ (T 0 T c ) 1/2 or T c ≤ T 2 /T 0 (i.e.n≤ m 2 if we introduce the dimensionless pulse length).The curve of n=m 2 (marked "T c =T 2 /T 0 ") and the curve of the optimal T co (marked "optimal T c ") obtained from Eq. ( 7) are shown in Fig. 2. One can find that the condition "T c ≤ T 2 /T 0 " is inaccurate as this condition includes the whole area under the curve n=m 2 and the curve itself, while our analysis indicates that only those chirped probe pulses with duration T co satisfy the curve obtained from equation ( 7) can be applied to reconstruct the original THz signal.
To verify the validity of our theory, we calculated a representative example with the original THz pulse length T =0.5ps.We assumed in our calculation that the Fourier-limited pulse length T 0 =25fs and T 0 =10fs.According to Eq. ( 7) the optimal chirped probe pulse lengths T co matched to above T can be calculated to be as 2.6437ps when T 0 =25fs and 4.1943ps, when T 0 =10fs, respectively.With these T and their T co , we calculated the THz fields according to Eq. ( 6).The results are shown in Fig. 3 in which T 0 =25fs in (a) and T 0 =10fs in (b).In Fig. 3, the dotted black lines represent the original THz fields, while the solid red lines are the retrieved THz fields calculated with the optimal probe pulse length T co .One can find that the retrieved THz fields with the T co match the original ones very well.
To see the behavior of the retrieved THz field due to the discrepancy ΔT c =T c -T co , we calculated four cases with T c = T co /2, T c = T co /4, T c = 2T co , and T c = 4T co which fulfill the condition T ≥ (T 0 T c ) 1/2 of Sun's theory 9 with the same T(=0.5ps).The results are also shown in Fig. 3.One can clearly see that all these cases introduce distortions and the larger the difference of T c to T co is, the more severe the distortion of the retrieved THz field will be.In detail, if T c <T co , the retrieved THz waveform will be compressed while its spectrum would be broaden compared to the input original THz field.Whereas, if T c >T co , the retrieved THz waveform will be broaden while its spectrum would be compressed.
In addition, one can find that the larger the T 0 is, the greater the difference of the retrieved THz waveform to the same discrepancy c T Δ will be if one compares Fig. 3(a In other words, if an input THz pulse is measured using different probe beams, the shorter the original probe pulse length (or the broader the spectrum of the original probe pulse) is, with the same discrepancy to their optimal chirped pulse length, the smaller the distortion of the retrieved THz signal will be.
We also simulated the whole retrieving process of the THz waveforms under the condition of k<<1 (for example k = 0.02).We find that the simulation results are the same to the ones shown in Fig. 3. Our simulations confirm our deduction process and prove the validity of Eq. (7).The detail of these simulation results will be published in the future.It should be noted that our analysis is based on a bipolar THz field.For other THz pulses such as with multiple cycles, we prove that there still is an optimal chirped probe pulse duration T co by means of our simulation, while the characteristic time of the THz field should to be redefined and the expression of Eq. (7) needs to be modified.The detailed theoretical deduction will be performed in the future work.
In summary, we conduct a rigorous deduction and a detailed analysis of the relationship between the chirped pulse length T c and the input bipolar THz pulse length T in the spectral-encoding technique.We prove that there is an optimal chirped pulse duration T co or an optimal chirp rate matched to the input THz pulse length T. We derive a relationship between the dimensionless optimal chirped pulse length n and the input dimensionless THz pulse length m.With this relationship the optimal duration and the optimal chirp rate of the chirped probe pulse can be calculated according to T co =nT 0 and respectively.We find that only under this restricted condition the measured THz signal can be retrieved without distortion.

Fig. 3 .
Fig. 3. (Color on line) Retrieved THz waveforms vs the original THz waveforms (black dotted lines) with T 0 =25fs in (a) and T 0 =10fs in (b).The red solid lines are the retrieved THz waveforms calculated according equation (6) under the conditions of input THz pulse length T = 0.5ps (m =20 when T 0 =25fs and m =50 when T 0 =10fs) and their corresponding optimal chirped pulse lengths T co =2.6437ps when T 0 =25fs, T co = 4.1943ps when T 0 =10fs, respectively, which are calculated from equation(7).The retrieved THz fields calculated according to the equation (6) with T c = T co /2, T c = T co /4, T c = 2T co , and T c = 4T co are also shown in each figure.
Fig. 3. (Color on line) Retrieved THz waveforms vs the original THz waveforms (black dotted lines) with T 0 =25fs in (a) and T 0 =10fs in (b).The red solid lines are the retrieved THz waveforms calculated according equation (6) under the conditions of input THz pulse length T = 0.5ps (m =20 when T 0 =25fs and m =50 when T 0 =10fs) and their corresponding optimal chirped pulse lengths T co =2.6437ps when T 0 =25fs, T co = 4.1943ps when T 0 =10fs, respectively, which are calculated from equation(7).The retrieved THz fields calculated according to the equation (6) with T c = T co /2, T c = T co /4, T c = 2T co , and T c = 4T co are also shown in each figure.