Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Nonellipsometric electro-optic sampling of terahertz waves in GaAs

Open Access Open Access

Abstract

We report efficient high-resolution electro-optic sampling detection of broadband terahertz radiation by measuring direct intensity modulation of a femtosecond fiber laser beam induced by terahertz field in a 1-cm thick GaAs crystal.

© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Free space electro-optic (EO) sampling of terahertz waveforms by femtosecond laser pulses [1] is a technique that is commonly used in terahertz time-domain spectroscopy [2], along with photoconductive sampling [3]. In the standard EO sampling scheme, the probe optical pulse co-propagates collinearly with the terahertz pulse in an EO crystal and accumulates a variation of polarization caused by the Pockels effect. By measuring the polarization variation as a function of delay between the optical and terahertz pulses, the terahertz electric field can be mapped. To achieve efficient EO sampling for a given wavelength λ of the probe optical beam, one should use a specific EO crystal that provides optical-terahertz velocity matching. For example, ZnTe is routinely used to perform EO sampling with a Ti:sapphire laser (λ ≈ 0.8 μm). Unfortunately, no crystal is known to provide velocity matching at the wavelength (λ ≈ 1.55 μm) of femtosecond fiber lasers, which are an appropriate light source for low-cost compact terahertz spectrometers.

In noncollinear extensions of the ellipsometric technique [4,5], the probe optical beam propagates in an EO crystal at the Cherenkov angle to the terahertz beam. This allows one to achieve synchronous propagation of the probe optical pulse and a terahertz wavefront for an arbitrary optical wavelength in crystals with an optical-terahertz velocity mismatch, such as, LiNbO3 [4] and GaAs (at the wavelengths larger than 1.33 μm) [5].

Recently, a nonellipsometric method of EO sampling has been proposed [6]. In this method, direct modulation of the probe optical beam intensity, induced by terahertz field in an EO crystal, is measured without any polarization optics. The mechanism of the nonellipsometric sampling has been elucidated as spatial separation of the sum-frequency-generation (SFG) and difference-frequency-generation (DFG) contributions of opposite polarities to the optical intensity modulation [7]. The separation is achieved due to noncollinear propagation of the optical and terahertz beams. By using a Si-prism-coupled LiNbO3 crystal, it was shown that noncollinear nonellipsometric sampling could provide efficiency comparable to that of the standard ellipsometric technique [6]. Later, spectral filtering was proposed as another way to frustrate mutual compensation of the SFG and DFG contributions in the collinear scheme [8].

In this paper, we demonstrate experimentally noncollinear nonellipsometric EO sampling with a femtosecond fiber laser (λ = 1.55 μm) in a GaAs crystal. Due to a smaller, than in LiNbO3, velocities mismatch between optical and terahertz waves in GaAs, we implemented noncollinear nonellipsometric scheme without a Si-prism coupler, unlike Ref. [6]. Such simplification of the scheme facilitates its practical applications in compact terahertz spectrometers. We obtained detection efficiency of the same order of magnitude as that of the noncollinear ellipsometric technique [5]. Using a 1-cm thick crystal ensured higher, as compared to mm- or sub-mm-thick crystals [9,10], spectral resolution (as high as a few GHz [5]). We derived theoretically and confirmed experimentally optimal orientation of the terahertz and optical polarizations with respect to the crystallographic axes of the crystal. This orientation provides a reasonable tradeoff between high modulation depth and practicality of implementation of the noncollinear nonellipsometric EO sampling in GaAs. We also verified experimentally the theory developed in Ref. [7].

2. Experimental

The experimental setup is shown in Fig. 1(a). A femtosecond Er3+–doped fiber laser with 1.55 μm central wavelength, 70 fs pulse duration (FWHM), and 100 MHz repetition rate was used as an optical source. The optical beam was split into the pump beam (35 mW average power), which triggered a photoconductive antenna (PCA) on an InGaAs/InAlAs substrate, and probe beam (30 mW) whose intensity was varied by using a half-wave plate (λ/2) and Glan prism (GP). PCA was biased with a ± 20 V, 10 kHz square wave voltage. By using a TPX lens and a parabolic mirror, the terahertz beam was collimated and focused onto the 1-cm thick <110>-cut GaAs crystal (with 5 × 5 mm2 cross-section). From knife-edge measurements with use of a Golay cell, the 1/e beam width was estimated to be ~1.2 mm. The corresponding Rayleigh length is as high as ~1.5 cm (at the frequency ~0.5 THz of the spectrum maximum). Therefore, diffraction broadening of the terahertz beam was negligible.

 figure: Fig. 1

Fig. 1 (a) Schematic of the experimental setup. (b) The polarizations of the optical and terahertz beams that optimize the efficiency of the nonellipsometric sampling.

Download Full Size | PDF

The probe beam was collimated by a 25.4 mm focal length lens f1 and focused onto the GaAs crystal by a lens f2, whose focal length was varied to provide different beam widths in the crystal. We used the 1/e-widths of 80, 120, 180, and 500 µm measured at the entrance surface of the crystal. The corresponding Rayleigh lengths (>2 cm) substantially exceed the crystal thickness. The incident angle of the probe beam on the entrance surface of the crystal was taken equal to 45° to ensure the propagation of the probe beam in the crystal at the Cherenkov angle β ≈ 12° to the terahertz beam [5].

The setup could be switched from the nonellipsometric scheme to the ellipsometric one [5] by using a flipped mirror (it also required rotating the GaAs crystal by 90°). In the ellipsometric scheme, a standard combination of a quarter-wave plate (λ/4) and Wollaston prism (WP) was used to separate orthogonal polarizations of the probe beam. To improve the signal-to-noise ratio in the nonellipsometric scheme, we used a D-shape mirror, which divided the probe beam into two halves, which correspond to the SFG and DFG processes, and directed them to a balanced detector. Due to opposite polarities of the SFG and DFG signals they are added constructively whereas the background signals from the two photodiodes are subtracted thus improving the signal-to-noise ratio. To minimize water vapor absorption, the measurements were made in dry air conditions.

In the ellipsometric operation mode, the [001] crystallographic axis of the crystal was aligned along the optical beam polarization and orthogonal to the terahertz beam polarization, i.e., at θ = −90° [Fig. 1(b)]. This orientation is optimal for ellipsometric detection in <110>-cut zinc-blende crystals [11–13]. In the nonellipsometric mode, we rotated the crystal by 90° to the position θ = 0° shown in Fig. 1(b). This orientation, as we explain below, may be optimal in terms of high optical intensity modulation depth and practical convenience. To verify experimentally a theoretically predicted dependence of the nonellipsometric EO signal on angle θ, we rotated the crystal around the [110] direction [Fig. 1(b)].

3. Results and discussion

Figure 2 shows the EO signals obtained by using the nonellipsometric method with different 1/e-widths (80, 120, 180, and 500 µm) of the probe beam in the crystal. The signal obtained by the ellipsometric method with an 80-µm wide probe beam is shown for reference. From Fig. 2, it is seen that the efficiencies of both methods are comparable but the shapes of the output signals differ significantly. For a narrow (80-µm wide) beam, the shape of the nonellipsometric signal is close to the time derivative of the ellipsometric one (inset in Fig. 2). With increasing the beam width, the nonellipsometric signal deviates from the time derivative shape and becomes smoother.

 figure: Fig. 2

Fig. 2 EO signals obtained by the nonellipsometric (NEL) method with different widths of the probe beam. EO signal obtained by the ellipsometric (EL) method is shown for reference. The inset compares the normalized shapes of the NEL 80 μm signal (solid) and the time derivative of the EL signal (dashed).

Download Full Size | PDF

Figure 3 shows the corresponding amplitude spectra. The experimental results in Fig. 3 confirm the theoretical prediction [7] that noncollinear nonellipsometric scheme operates as a bandpass filter, whose transmission band is defined by the probe beam width. According to Ref. [7], the output signal of the nonellipsometric scheme can be evaluated as

s(τ)0dΩH(Ω)Im[ETHz(Ω)eiΩτ],
where ETHz(Ω) is the Fourier transform of the input terahertz field and H(Ω)=exp(Ω2/Ωc2)erf(Ω/Ωc) is the transfer function, with erf(ξ) the ordinary error function and Ωc the cutoff frequency. The cutoff frequency is Ωc=22c/(anTHzsinβ), where a is the 1/e2 half-width of the probe beam (2 times smaller than the 1/e-widths we use here), nTHz = 3.6 is the terahertz refractive index of GaAs [5], and c is the speed of light. Figure 4 shows H(Ω) for the probe beam widths used in the experiment. The spectrum of the EO signal obtained by the ellipsometric method with an 80-µm wide beam is shown for reference.

 figure: Fig. 3

Fig. 3 Amplitude spectra of the EO signals shown in Fig. 2.

Download Full Size | PDF

 figure: Fig. 4

Fig. 4 Transfer function of the nonellipsometric scheme for different widths of the probe beam. The spectrum of the ellipsometric EO signal is shown for reference.

Download Full Size | PDF

According to Figs. 3 and 4, for the 500-µm wide beam, the scheme operates as a low-pass filter that attenuates high frequencies. Physically, it can be explained by the fact that in noncollinear geometry the pancake-shaped probe pulse is tilted to the terahertz wavefronts, its different parts interact with different phases of the terahertz field. This reduces the EO response, the effect is more pronounced for a wider probe beam and smaller terahertz wavelengths. For the 80-µm wide beam, high frequencies are boosted and low frequencies are filtered out (Figs. 3 and 4). The weakening of the response to low frequencies can be explained by a smaller angular separation of the SFG and DFG contributions to the optical intensity modulation from low frequencies relative to the divergence angle of the probe beam [7]. The optimal width of the probe beam for correctly measuring the spectrum of the used terahertz source (PCA) is ~120 µm (Fig. 3).

The shapes of the nonellipsometric EO signals (Fig. 2) can be explained by using Eq. (1) and Fig. 4. For the 80-µm wide probe beam, the terahertz spectrum lies almost completely at the frequencies significantly lower than the maximum frequency of the transfer function (Fig. 4), where H(Ω) ≈ 2π-1/2Ω/Ωc. Taking into account the phase shift of all spectral components of the input spectrum by π/2 in Eq. (1) (because of Im), the output spectrum is proportional to the input spectrum multiplied by iΩ. This corresponds to differentiating the input terahertz signal with respect to time. With increasing the probe beam width, the output signal deviates from the time derivative shape and becomes smoother.

In the standard collinear geometry, a short coherence length of laser pulses at 1.55 μm and terahertz waves in GaAs leads to the necessity of using crystals thinner than ~1 mm. For thin crystals, an echo, originating from internal reflections of the terahertz pulse in the crystal, is separated from the main EO signal by a short time interval – about 12 ps for a 0.5 mm thick crystal [2]. The necessity to filter out the signal echo puts the upper limit on the time window of EO sampling and, as a result, significantly restricts the spectral resolution – by 40 GHz for a 1 mm thick crystal. In the noncollinear geometry, using cm-thick crystals allows one to achieve an order of magnitude increase in the spectral resolution [5]. In particular, the spectral resolution in Fig. 3 is as high as ~3 GHz.

To provide a highest efficiency of the nonellipsometric scheme, the terahertz and optical beam polarizations should be different from those in the ellipsometric scheme. The optimal polarizations, shown in Fig. 1(b), can be found from following consideration. Fourier transform (with respect to time) of the optical electric field after the GaAs crystal can be written as

Eopt(ω)=E(0)(ω)+ΔE(ω),
where E(0)(ω) is the optical field transform in the absence of the terahertz pulse and ΔE(ω) is a small perturbation arising from the interaction of the optical and terahertz fields in the crystal, |ΔE(ω)||E(0)(ω)|. The perturbation ΔE(ω) can be found by solving Fourier-transformed Maxwell’s equations with the nonlinear polarization
PiNL(ω)=2ε00dΩχijk(2)Ej(0)(ωΩ)EkTHz(±Ω)
included as a source (χijk(2) is the second-order nonlinear susceptibility tensor and ε0 is the vacuum permittivity) [7]. It should be noted that the SFG and DFG contributions to the sum in Eq. (3) are spatially separated and, therefore, do not compensate each other [7]. Vector ΔE(ω) is directed along PNL(ω). By splitting ΔE(ω) into two orthogonal parts ΔE(ω)=ΔE(ω)+ΔE(ω) with ΔE(ω)E(0)(ω) and ΔEE(0)(ω), we write the optical intensity at the frequency ω as
I(ω)|E(0)(ω)+ΔE(ω)|2+|ΔE(ω)|2|E(0)(ω)|2+2Re[E(0)(ω)ΔE*(ω)],
where small quadratic terms are neglected. The optical intensity modulation depth
ΔI(ω)=I(ω)I(0)(ω)2Re[E(0)(ω)ΔE*(ω)],
with I(0)(ω)|E(0)(ω)|2, is proportional to ΔE(ω) and practically independent of ΔE(ω). Thus, to achieve a highest modulation depth ΔI, one should choose such polarizations of the optical and terahertz fields that maximize the component of PNL(ω), which is directed along E(0)(ω), and therefore maximize ΔE(ω). In zinc-blende crystals, such as GaAs, where the tensor χijk(2) has three nonzero components d14 = d25 = d36 [11], the above formulated requirement can be fulfilled for parallel polarizations of the optical and terahertz beams directed at ≈55° to the [001] axis [14]. This efficiency-maximizing configuration is not, however, convenient for practical use. For a standard rectangular shape of a <110>-cut crystal, such as in Fig. 1(b), it requires a special holder to rotate the crystal and adjust its [001] axis at 55° to the beam polarizations. Instead, we propose to use the configuration shown in Fig. 1(b). This configuration provides only slightly (≈15%) lower optical intensity modulation depth but at the same time does not require any adjustment. Thus, it may be a reasonable tradeoff between high modulation depth and practicality, in particular, for compact terahertz spectrometers.

To corroborate the consideration above, we derived the dependence of the output nonellipsometric signal ΔI on the crystal rotation angle θ [Fig. 1(b)] for fixed polarizations of the optical and terahertz fields [Fig. 1(b)]. This dependence coincides with the θ-dependence of the PNL(ω) component parallel to E(0)(ω). By using the tensor χijk(2) [11], we obtained the dependence ΔIcosθ(23cos2θ) that was confirmed experimentally (Fig. 5).

 figure: Fig. 5

Fig. 5 The output nonellipsometric EO signal as a function of θ: experiment (crosses) and theory (solid line).

Download Full Size | PDF

4. Conclusion

We have demonstrated efficient broadband EO sampling of terahertz pulses in a GaAs crystal without any polarization and coupling optics. The optimal, in terms of high modulation depth and practical convenience, polarizations of the terahertz and probe optical beams have been theoretically found and experimentally verified. Using a cm-thick crystal provided spectral resolution as high as a few GHz. The measured dependencies of the output terahertz waveform and spectrum on the probe beam width agree well with the theoretical predictions [7].

Funding

Ministry of Education and Science of the Russian Federation (3.3854.2017/4.6, 3.6252.2017/9.10); Japan Science and Technology Agency (20100220).

References and links

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

2. A. Tomasino, A. Parisi, S. Stivala, P. Livreri, A. C. Cino, A. C. Busacca, M. Peccianti, and R. Morandotti, “Wideband THz time domain spectroscopy based on optical rectification and electro-optic sampling,” Sci. Rep. 3(1), 3116 (2013). [CrossRef]   [PubMed]  

3. D. H. Auston, K. P. Cheung, and P. R. Smith, “Picosecond photoconducting Hertzian dipoles,” Appl. Phys. Lett. 45(3), 284–286 (1984). [CrossRef]  

4. M. Tani, K. Horita, T. Kinoshita, C. T. Que, E. Estacio, K. Yamamoto, and M. I. Bakunov, “Efficient electro-optic sampling detection of terahertz radiation via Cherenkov phase matching,” Opt. Express 19(21), 19901–19906 (2011). [CrossRef]   [PubMed]  

5. E. A. Mashkovich, A. I. Shugurov, S. Ozawa, E. Estacio, M. Tani, and M. I. Bakunov, “Noncollinear electro-optic sampling of terahertz waves in a thick GaAs crystal,” IEEE Trans. THz Sci. Technol. 5(5), 732–736 (2015).

6. M. Tani, T. Kinoshita, T. Nagase, K. Horita, C. T. Que, E. Estacio, K. Yamamoto, and M. I. Bakunov, “Non-ellipsometric detection of terahertz radiation using heterodyne EO sampling in the Cherenkov velocity matching scheme,” Opt. Express 21(8), 9277–9288 (2013). [CrossRef]   [PubMed]  

7. M. I. Bakunov, S. D. Gorelov, and M. Tani, “Nonellipsometric noncollinear electrooptic sampling of terahertz waves: A comprehensive theory,” IEEE Trans. THz Sci. Technol. 6(3), 473–479 (2016).

8. I. E. Ilyakov, G. Kh. Kitaeva, B. V. Shishkin, and R. A. Akhmedzhanov, “Terahertz wave electro-optic measurements with optical spectral filtering,” Appl. Phys. Lett. 106(12), 121101 (2015). [CrossRef]  

9. M. Nagai, K. Tanaka, H. Ohtake, T. Bessho, T. Sugiura, T. Hirosumi, and M. Yoshida, “Generation and detection of terahertz radiation by electro-optical process in GaAs using 1.56 μm fiber laser pulses,” Appl. Phys. Lett. 85(18), 3974–3976 (2004). [CrossRef]  

10. Z. Zhao, A. Schwagmann, F. Ospald, D. C. Driscoll, H. Lu, A. C. Gossard, and J. H. Smet, “Thickness dependence of the terahertz response in <110>-oriented GaAs crystals for electro-optic sampling at 1.55 microm,” Opt. Express 18(15), 15956–15963 (2010). [CrossRef]   [PubMed]  

11. Q. Chen, M. Tani, Z. Jiang, and X.-C. Zhang, “Electro-optic transceivers for terahertz-wave applications,” J. Opt. Soc. Am. B 18(6), 823–831 (2001). [CrossRef]  

12. N. C. J. van der Valk, T. Wenckebach, and P. C. M. Planken, “Full mathematical description of electro-optic detection in optically isotropic crystals,” J. Opt. Soc. Am. B 21(3), 622–631 (2004). [CrossRef]  

13. R. D. Santos, S. Ozawa, V. Mag-Usara, S. Azuma, A. Tuico, V. Copa, A. Salvador, K. Yamamoto, A. Somintac, K. Kurihara, H. Kitahara, M. Tani, and E. Estacio, “Cherenkov-phase-matched nonlinear optical detection and generation of terahertz radiation via GaAs with metal-coating,” Opt. Express 24(22), 24980–24988 (2016). [CrossRef]   [PubMed]  

14. S. P. Kovalev and G. Kh. Kitaeva, “Terahertz electro-optical detection: optical phase or energy measurements,” J. Opt. Soc. Am. B 30(10), 2650–2656 (2013). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1
Fig. 1 (a) Schematic of the experimental setup. (b) The polarizations of the optical and terahertz beams that optimize the efficiency of the nonellipsometric sampling.
Fig. 2
Fig. 2 EO signals obtained by the nonellipsometric (NEL) method with different widths of the probe beam. EO signal obtained by the ellipsometric (EL) method is shown for reference. The inset compares the normalized shapes of the NEL 80 μm signal (solid) and the time derivative of the EL signal (dashed).
Fig. 3
Fig. 3 Amplitude spectra of the EO signals shown in Fig. 2.
Fig. 4
Fig. 4 Transfer function of the nonellipsometric scheme for different widths of the probe beam. The spectrum of the ellipsometric EO signal is shown for reference.
Fig. 5
Fig. 5 The output nonellipsometric EO signal as a function of θ: experiment (crosses) and theory (solid line).

Equations (5)

Equations on this page are rendered with MathJax. Learn more.

s( τ ) 0 dΩ H( Ω )Im[ E THz ( Ω ) e iΩτ ],
E opt (ω)= E (0) (ω)+ΔE(ω),
P i NL ( ω )=2 ε 0 0 dΩ χ ijk ( 2 ) E j ( 0 ) ( ωΩ ) E k THz ( ±Ω )
I(ω) | E (0) (ω)+Δ E (ω) | 2 + | Δ E (ω) | 2 | E (0) (ω) | 2 +2Re[ E (0) (ω)Δ E * (ω) ],
ΔI(ω)=I(ω) I (0) (ω)2Re[ E (0) (ω)Δ E * (ω) ],
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.