• Photonics Research
  • Vol. 13, Issue 6, 1736 (2025)
Maximilian Lenz* and Pietro Musumeci
Author Affiliations
  • Department of Physics and Astronomy, University of California at Los Angeles, Los Angeles 90095-1547, USA
  • show less
    DOI: 10.1364/PRJ.558770 Cite this Article Set citation alerts
    Maximilian Lenz, Pietro Musumeci, "Single-shot electro-optic sampling with arbitrary terahertz polarization," Photonics Res. 13, 1736 (2025) Copy Citation Text show less

    Abstract

    With the recent development of diversity electro-optic sampling (DEOS), significant progress has been made in the range of applicability of single-shot EOS measurements, allowing broadband THz waveforms to be captured in a single shot over large temporal windows. In addition to the decrease in acquisition time compared to standard multishot data acquisition, this technique allows measurements on systems far from equilibrium with large shot-to-shot noise or with irreversible or poorly repeatable dynamics. Although DEOS has been demonstrated and verified for linearly polarized THz waveforms, we investigate the effects resulting from the presence of a secondary polarization component. This imposes new challenges for accurate waveform reconstruction, and opens the opportunity to measure out complex polarization states such as arbitrary elliptically polarized THz field. We demonstrate a single-shot diversity-electro-optic-sampling-based approach to capture both x- and y-THz fields simultaneously with a single (110)-cut EO crystal for THz polarimetry and ellipsometry over a wide range of frequencies.

    1. INTRODUCTION

    The ability to capture terahertz (THz) waveforms has led to the development of many important research tools such as THz time-domain spectroscopy (TDS), material characterization, biomedical imaging, electron beam diagnostics, and scanning tunneling microscopy [14]. Electro-optic sampling (EOS) is a commonly used technique that measures out the field-induced birefringence in an electro-optic crystal using a probe infrared laser pulse, allowing direct measurements of the temporal electric field waveform in the THz frequency range [58].

    Traditional multishot schemes where the time evolution of the birefringence is mapped by scanning the delay of the ultrashort laser pulse over several shots face two main challenges. First, these are not applicable in low repetition rate cases or in setups with large shot-to-shot time-of-arrival fluctuations between the THz and the probe laser. These importantly include most electron-beam-based field sources, such as THz free electron lasers, coherent transition and synchrotron radiation, and near-fields of the beam itself [915]. Even when the THz source and probe laser are perfectly synchronized, multishot TDS systems still fail when analyzing samples that are dynamically changing with poor repeatability [1622], highlighting the necessity of single-shot detection schemes. Secondly, in its typical implementation, multishot EOS measures the electric field along one axis, not allowing arbitrary polarization states to be measured without reorienting the optical setup and performing a second independent scan.

    The first challenge can be addressed using single-shot setups where the waveform is encoded either spatially or spectrally in the IR pulse. Spatial encoding requires non-collinear propagation in the crystal and typically can be useful only over short time windows, limiting the application range. In spectral encoding approaches, the probe pulses are stretched in time by adding a frequency chirp so that the THz waveform information can be collected over very long (100s of ps) time windows in a single shot by simply recording the probe pulse spectrum [2325]. However, this technique has been severely limited by distortions that appear when trying to capture wide frequency ranges over very long time windows [26]. With the new advancement of diversity electro-optic sampling (DEOS) by Roussel et al. [27], these limitations can be overcome and accurate waveform reconstruction of long THz waveforms with broad spectral ranges is now achievable. The approach essentially relies on the use of a two-channel encoding scheme with a shifted frequency response so that combining the information from the two channels makes it possible to reconstruct the full waveform.

    The second challenge related to the THz field polarization reconstruction has been addressed in multishot EOS schemes by mechanically rotating the optical components of the readout [28,29]. For single-shot polarimetry, a general formalism based on probe pulse stretching was developed in the visible range, avoiding the need for mechanical intervention in the setup [30]. This concept was recently adapted to the THz range [31] allowing the capture of a THz waveform with polarization information in a single-shot EOS scheme, but without addressing the limitations associated with spectral-temporal encoding [26,27].

    In this paper, we will demonstrate how the DEOS formalism can be extended to arbitrary polarization states, allowing THz polarization measurements while still addressing the waveform distortions that appear in the spectral-temporal encoding scheme. We present an analytical model for the measured signals and a reconstruction method for THz waveforms that are polarized at an arbitrary angle or ellipticity. We also highlight additional effects in DEOS measurements that alter the reconstruction results. Finally, using an optical-rectification-based THz source, we demonstrate the use of a stretched pulse EOS setup to measure the THz wave polarization state, allowing THz time-domain ellipsometry to be performed in a single shot over a wide frequency range.

    2. REVIEW OF DIVERSITY ELECTRO-OPTIC SAMPLING

    We start the discussion with a quick review of the EOS technique and its phase-diversity implementation to establish the formalism, which we will later use to retrieve the polarization of the THz electric field. For an electro-optic crystal cut in the (110) plane, a THz field incident at angle α with respect to the [110] axis results in a birefringent index-ellipsoid with slow and fast axes in the (110) plane [see Fig. 1(a)] having refractive index [27,32] n1,2(α)=n0+η1,2(α)ETHz,where η1,2(α)=n03r414(sinα±1+3cos2α), n0 is the refractive index in the absence of an electric field, r41 is the EO coefficient of the crystal, and ETHz is the applied THz electric field. The orientation of the index ellipsoid is determined by the angle Ψ between the slow axis and the [110] axis, related to α by tan2Ψ=2cotα [32] and portrayed in Fig. 1. The birefringence Δn(α)=n1(α)n2(α) is maximized for α=0 as in the standard configuration for EOS [see Fig. 1(b)].

    ZnTe crystal orientation for arbitrary angle (a), standard/balanced configuration (b), and DEOS configuration (c).

    Figure 1.ZnTe crystal orientation for arbitrary angle (a), standard/balanced configuration (b), and DEOS configuration (c).

    The induced birefringence can be read out measuring the polarization state of the probe pulse, which can be modeled using Jones-matrix formalism [33]. In the crystal frame, the EO-induced matrix is diagonal and can be written as MXtal=[eiϕ1(α,ETHz)00eiϕ2(α,ETHz)],where ϕ1,2(α,ETHz)=k0L(n0+η1,2(α)ETHz), where L is the length of the crystal and k0 is the wavenumber of the probe pulse. Typically, after the crystal an EOS setup includes a quarter-waveplate (QWP, η=π/2) followed by a half-waveplate (HWP, η=π) with corresponding Jones matrices W(η,θ)=eiη/2[cos2θ+eiηsin2θ(1eiη)sinθcosθ(1eiη)sinθcosθsin2θ+eiηcos2θ],so that the full change in polarization state of the probe pulse is described by the single matrix M=W(π,θH)W(π/2,θQ)R1(Ψ)MXtalR(Ψ),in which R is the rotation matrix and θH, θQ are the angles of the waveplates’ fast axes with respect to the [110] axis of the EO crystal [33]. The Jones transport operates on the incoming probe electric field vector resulting in an output signal: Es=ME0(cosθIRsinθIR),where the incident field has magnitude E0 and is polarized at angle θIR with respect to the [110] axis. For the standard configuration we set α=0 [i.e., Ψ=π/4; see Fig. 1(b)] and use a linearly polarized IR probe pulse with θIR=0, a QWP at θQ=π/4, and an HWP at θH=0, so that the modulated probe pulse becomes (Es,1Es,2)=E02(cos(Γ2)sin(Γ2)i(cos(Γ2)+sin(Γ2)))12(1Γ2i(1+Γ2))for small phase shifts. Γ is related to the applied THz field by Γ=ϕ1ϕ2=k0Ln03r41TEO=:βETHz,where TEO is a correction term that accounts for Fresnel losses, absorption, and dephasing effects in the crystal [34,35]. We define β as the calibration factor between the electric field and the induced phase shift. By separately measuring the intensity of the two polarizations Is,i|Es,i|2, in the limit of small Γ, we can retrieve the phase shift as Γ=Is,2Is,1Is,2+Is,1and hence the applied electric field by inverting Eq. (7). The standard configuration results in the strongest modulation of the probe pulse [32,34], but as discussed below, in a single-shot spectrally encoded scheme where the probe pulse is stretched, temporal distortions appear and other configurations might be favored [27].

    A. Single-Shot Phase Diversity Electro-Optic Sampling

    To better understand the bandwidth limitations in spectrally encoded setups as well as the phase diversity solution, we review the principles of photonic time stretching in more detail [36]. A chirped probe pulse can be modeled as E0(t)=A(t)exp(iC2t2iω0t),with A(t) being the pulse envelope, C=ω/t the chirp parameter, and ω0 the center frequency. For a linear chirp, the spectral-temporal correlation is expressed as t=(ωω0)/C. The corresponding spectral amplitude is obtained by Fourier transform for a slowly varying A(t) as [27] E˜0(ω)i2πCA(ωω0C)exp(i(ωω0)22C).

    The modulation of the probe pulse by the EO effect in the crystal can be expressed as Es(t)=T(t)E0(t), with T(t) being the THz-induced modulation function, which is written in a general form for small phase retardation Γ(t) as T(t)=T(1+meiψΓ(t)).

    For example, for the two polarizations in the standard configuration case the transfer function coefficients can be extracted from Eq. (6) as T1=1/2,T2=i/2,m1,2=1/2,ψ1,2=0. The signal spectrum is given by the convolution integral E˜s(ω)=T˜(ω)*E˜0(ω) and after background subtraction and normalization the recorded signal can be written as Y(ω)=|E˜s(ω)|2|TE˜0(ω)|2|TE˜0(ω)|2=2e[E˜0*(ω)meiψΓ˜E˜0(ω)]|E˜0(ω)|2.

    Using the expression from Eq. (10) and assuming that the slowly varying spectral amplitude A can be carried out of the integral and simplified away by the normalization (see Appendix A), we have Y(ω)=2eΓ˜(Ω)meiΩ2/C+ψeiΩ(ωω0)/CdΩF(H(Ω)Γ˜(Ω)),where Γ˜(Ω) is the Fourier transform of the applied THz field and the Fourier transform F is intended to be taken with respect to the stretched time t, correlated to ω by the frequency chirp. In this picture, H(Ω) represents the frequency response transfer function of the measurement technique. In fact, if H(Ω)=1, then the measured spectrum signal would be an exact copy of Γ(t) and we would have perfect spectral-temporal encoding. In practice, the chirp C affects our ability to retrieve the signal Γ as, after using the Hermitian property of Γ(Ω) to simplify the real part operator, we have H(Ω)=2mcos(Ω22C+ψ),and all spectral information is lost for frequencies in which H(Ω)=0. A detailed derivation of the transfer function is presented in Appendix A.1. In the standard EOS configuration H1,2(Ω)=cos(Ω22C) and the signals retrieved from each probe pulse polarization have the same zeros so that no information can be retrieved for those frequencies. As we increase the measurable time window by stretching the probe pulse, we decrease C, causing the zeros of the transfer function to move to lower frequencies as seen in Fig. 2 where we plot the absolute value of the transfer function H(Ω) for different chirp parameters C. As an example, for C=1, we would have no way to retrieve electromagnetic field information at 282 GHz, 488 GHz, 631 GHz, and so on.

    Transfer function values in absolute units for different chirp parameters and different frequencies.

    Figure 2.Transfer function values in absolute units for different chirp parameters and different frequencies.

    Phase diversity EOS makes use of the parameter ψ in Eq. (14) to enable full spectral reconstruction. By adjusting the polarization readout away from the standard configuration, we can tune ψ1 and ψ2 for the two separate polarizations Es,1, Es,1 so that the zeros of the transfer functions H1(Ω) and H2(Ω) are offset from each other. For example, by orienting the THz field to be perpendicular to the [001] axis (i.e., α=π/2, Ψ=0) and then setting our incoming probe pulse polarization to θIR=π/4 and adjusting θQ=0 and θH=π/8 as in Fig. 6(c) (shown later), it can be shown that the two transfer functions take on the form of [27] H1,2PD=12cos(Ω22Cπ4).(Our definition of Γ leads to a factor of two difference in the DEOS transfer function compared to that of Roussel et al.). This polarization readout is referred to as the DEOS configuration and allows the full waveform to be recovered, as Y1 and Y2 now contain complementary information. The spectral profile of the THz signal is then recovered via Γ˜(Ω)=H1(Ω)Y˜1(Ω)+H2(Ω)Y˜2(Ω)H12(Ω)+H22(Ω).

    The temporal profile is then obtained with an inverse Fourier transform and has been shown to result in accurate waveform reconstruction in tabletop experiments where multishot and single-shot EOS measurements can be easily compared [27].

    It is important to note that the spectral intensity of the unmodulated probe pulse |E˜0(ω)|2 in the DEOS case cannot be recorded as easily as in the standard configuration where |E˜0(ω)|2=|E˜s,1std(ω)|2+|E˜s,2std(ω)|2. For DEOS, a reference line is necessary to capture this normalization spectrum and retrieve Y1,Y2 separately. Another drawback of the proposed DEOS polarization readout is that the induced response is weaker, and thus signal-to-noise levels are lower.

    3. THEORY OF SINGLE-SHOT POLARIMETRY

    Let us now investigate the cases of arbitrary THz polarization, studying how adding an independent orthogonal polarization affects the measured EOS signal. We assume that the fields in two perpendicular directions independently contribute to the induced birefringence. This can be modeled by introducing a secondary phase shift Γy=βETHzy, where β is the calibration factor from Eq. (7). Referring to the Jones calculus described in Eq. (4), the second field polarization can be taken into account by applying the transport matrix M after an additional transport element R(Ψ(αy))MXtal(αy)R(Ψ(αy)), where αy=α+π/2. In the standard EOS configuration, when assuming small phase shifts Γx,Γy1, the modulated probe pulse is then described by (see Appendix A.2 for full expression) EsstdE02(1Γx2+iΓy2i(1+Γx2+iΓy2)),and using the same reasoning that led to Eqs. (13) and (14), the Fourier transforms of the observed signals result in Y˜1,2std=cos(Ω22C)=:HxΓ˜xsin(Ω22C)=:HyΓ˜y.

    It is then immediate to note that in the standard EOS analysis using Eq. (8) the contribution of the y-polarization is canceled and Eq. (7) only allows us to retrieve Γx. Second, a spectral signature of the y-polarization is only present in the signal in the case of a chirped pulse, as in the limit of a short (or compressed) pulse the transfer functions Hx(Ω)C1 and Hy(Ω)C0.

    In the case of linear polarization incident at an angle of 45° with respect to the x-axis (α=π4), Γx=Γy=12Γ and we note that we can rewrite Y˜1,2std as cos(Ω22C±π4)Γ˜ as we would expect when computing the transfer functions in the standard configuration for α=π/4 instead of α=0. This proves that our extended model is in agreement with the DEOS formalism, and this particular case shows that diversity arises in the standard configuration when rotating the THz polarization from α=0 to α=π/4.

    In a DEOS setup, in the presence of arbitrary THz polarization, we calculate the modulated probe pulse to be (see Appendix A.2 for full expression) EsPDE02(1(1i)Γx4+iΓy2i(1+(1+i)Γx4+iΓy2)),which results in the relationship Y˜1,2PD=12cos(Ω22Cπ4)Γ˜xsin(Ω22C)Γ˜y.

    It is worth noting here that, in the presence of ΓyEyTHz, the reconstruction from Eq. (16) fails, as it is contaminated by the presence of Γy that requires separate extraction. In fact, in the presence of a vertical polarization field EyTHz, DEOS-based reconstruction of Γx is not possible without prior knowledge of Γy.

    This can be seen in general by writing the relationship between the applied THz field and the measured signals as (Y˜1Y˜2)=(H1xH1yH2xH2y)(Γ˜xΓ˜y).

    In order to extract Γ˜x and Γ˜y from the observed signals Y˜1, Y˜2, the matrix needs to be invertible, so this can be done only for those frequencies within the observed signal range for which the matrix determinant is non-zero. Note that in both the standard and DEOS configurations, this is not the case and the retrieval of the two fields is not possible. In fact, full reconstruction of both Γx and Γy simultaneously and independently is an open challenge, as there is no known polarization readout in which the transfer function matrix of Eq. (21) has a non-zero determinant over the entire spectral range.

    A. Known Arbitrary Polarization

    In case the polarization state of the incoming THz wave is known, full waveform reconstruction is possible. The incident field in this case can be described without loss of generality by Γ˜x=rΓ˜ and Γ˜y=1r2eiϕΓ˜, with r[0,1] and ϕ[π,π). For example, if r=1, we retrieve linear polarization along x and for r=12, ϕ=±π2, we have a left- or right-handed circular polarization.

    In the standard EOS configuration, the measured observables take on the form Y˜1std=(rcos(Ω22C)1r2eiϕsin(Ω22C))=:H1,mstdΓ˜,Y˜2std=(+rcos(Ω22C)1r2eiϕsin(Ω22C))=:H2,mstdΓ˜.

    Since these newly defined transfer functions H1,m,H2,m have offset zeros, we can reconstruct the unknown spectral amplitude Γ using an equivalent of Eq. (16) by substituting H1,H2 with H1,m,H2,m. Note that if some ellipticity is present (i.e., r0,1 and ϕ0,π), we obtain complex-valued transfer functions with |H1,m|,|H2,m|0 over the entire spectral range, and the waveform can be retrieved independently from each channel simply by dividing out the transfer functions. The reconstruction would work in a similar way in the DEOS configuration after substituting the correct transfer functions from Eq. (20) in H1,mPD and H2,mPD. Interestingly, this suggests the possibility of reconstructing the full THz waveform in a single shot using spectral-temporal encoding simply by using a polarizer and a properly oriented well-characterized THz waveplate on the THz path.

    A numerical example is presented in Fig. 3 in which we demonstrate reconstruction with different configuration readouts. We start from the spectrum of a single-cycle THz pulse Γ˜act. with a bandwidth extending to 1.5 THz. The probe pulse is an 800 nm, 40 fs pulse, linearly stretched to 18 ps (FWHM), resulting in a chirp parameter C=3.413ps2. We simulate the modulation of the probe pulse using Jones matrices and obtain the measured signals Y˜1,Y˜2 from Eq. (13). The original Γ˜act. is then reconstructed using the appropriate choice of transfer functions (i.e., standard or phase diversity) depending on the EOS configuration readout. In Fig. 3 top left, we show the simulation of a standard configuration polarization readout with the THz polarization adjusted by 45° with respect to the [110] axis (α=π/4) so that both Γ˜x,y are present. In the bottom left, we show the simulation of the DEOS polarization readout with the THz polarization perpendicular to the [110] (α=π/2) as prescribed for the DEOS setup. The Fourier transforms of the measured signals Y˜1,Y˜2 are plotted together with the corresponding transfer function H1,H2 to highlight the lack of response in the measurement at certain frequencies. Note that the curves in the standard 45° and DEOS 0° cases have the same zeros (highlighted by dashed vertical lines), but differ in amplitude, allowing higher signal-to-noise ratios to be achieved in the standard 45° configuration. In both cases, complementary information is present in the two channels and it is possible to fully reconstruct the original signal. On the right, we present a simulated dataset in the DEOS polarization readout with the THz polarization rotated 45° away from the expected THz polarization, causing both Γx and Γy to be present. The zeros of the two transfer functions are closer together, which is undesirable for reconstruction, as information is lost over wider regions of the spectrum. We also show the accurately reconstructed signal (enabled by prior knowledge of the THz polarization), but highlight the reconstruction error that could occur by neglecting the additional polarization component Γy and applying the classical DEOS transfer functions from Eq. (15) for reconstruction.

    Simulated spectral datasets for different polarization configurations. The top left shows the expected signals under the standard configuration with α=π4 and the bottom left shows the classic DEOS configuration. The zeros of the transfer functions (equal in these two configurations) are indicated with dotted lines. On the right the simulated DEOS polarization readout with THz polarized at α=π4 is used to showcase the difference between a correct reconstruction resulting from the knowledge of the incoming THz polarization and a traditional DEOS reconstruction.

    Figure 3.Simulated spectral datasets for different polarization configurations. The top left shows the expected signals under the standard configuration with α=π4 and the bottom left shows the classic DEOS configuration. The zeros of the transfer functions (equal in these two configurations) are indicated with dotted lines. On the right the simulated DEOS polarization readout with THz polarized at α=π4 is used to showcase the difference between a correct reconstruction resulting from the knowledge of the incoming THz polarization and a traditional DEOS reconstruction.

    B. Ellipsometry

    Let us now study in detail the case in which the polarization state is unknown and needs to be retrieved as in the case for THz-ellipsometry application. For simplicity, we assume a birefringent sample of thickness d with a complex index ellipsoid defined by nx*, ny* in which a THz pulse propagates at 45° between the two main axes, so that the incident THz pulse is described by Γ˜x,yin(Ω)=Γ˜(Ω)2. The transmitted THz pulse is attenuated due to the absorption in the material as well as Fresnel losses at the dielectric interfaces and acquires a phase shift due to the refractive indices of the sample, resulting in an output EOS waveform described by (Γ˜x(Ω)Γ˜y(Ω))=Γ˜2(τx  exp(iΦx(Ω))τy  exp(iΦy(Ω))).

    Note that in general, the radiation pulse could experience multiple reflections in the sample, resulting in a somewhat complicated dependency of the transmission coefficient from the real nx,y and imaginary κx,y parts of the index of refraction. In these cases, the reconstruction of the dielectric response becomes a more involved task [37]. Often in single-shot EOS setups, the reflections are not even recorded as the induced time delay can be larger than the temporal window of the measurement. Here we just assume that absorption is small and neglect the Fabry-Perot effect so that we can write τx,y=4nx,y(nx,y+1)2eΩκx,yd/c as the transmission coefficient and Φx,y(Ω)=Ωc(nx,y1)d.

    By comparing the measured polarization components before and after inserting the sample we can retrieve the induced attenuation and phase shifts on the two polarizations and hence the sample birefringence. In practice, we need to be careful of the transfer functions in this process.

    In the case of standard configuration single-shot EOS, we rearrange the measured waveforms Hx,yΓ˜x,y=(Y˜2Y˜1)/2 to separate the orthogonal polarization contributions with Hx,y being the transfer functions defined in Eq. (18). Then multiplying by the respective filter functions Hx,y, we can leverage Hx,y20 to reliably extract the phase Φx,y(Ω)=arg(Hx,y(Y˜1Y˜2)2Γ˜)=arg(2Hx,y2Γ˜x,yΓ˜)after normalizing by the spectrum measured without the sample. This works as long as Hx,y(Ω)0, which defines the range of frequencies over which we can perform sample ellipsometry. It is important to note that Φ(π,π) due to the definition of the arg-function on the complex plane. Therefore, in order to retrieve nx,y by inverting Eq. (25) we need to first unwrap the phase so that it is a continuous function across the entire spectrum. Numerically, unwrapping the phase for a given dataset is a delicate task, as experimental noise can cause unwanted extrapolations.

    Once the real part of the index of refraction is recovered, we can use κx,y(Ω)=cΩdRe(lnHx,y(Y˜1Y˜2)2Γ˜lnHx,y2ln4nx,y(nx,y+1)2)to retrieve the absorption coefficient.

    For common materials that have a steady and continuous refractive index in the spectral region of interest, it is possible to extrapolate the refractive index in the regions where the spectral response of the measurement is minimal. Instead of extracting the refractive indices analytically, in fact, the problem can be approached numerically, modeling nx,y as a second-order polynomial in Ω and κx,y as a constant in the frequency domain. The unknown model parameter (the polynomial coefficients of the real index of refraction and its imaginary part in x and y, respectively) can be found by minimizing the difference between the computed and measured spectral profiles.

    A simulated dataset is presented in Fig. 4 for a single-cycle THz waveform Γ=tσtexp(t2σt2), with σt=0.5  ps. The EOS signals are collected after numerically propagating the waveform through a birefringent crystal that has real and complex indices of refraction different in x and y and with a parabolic frequency dependence, as shown by the solid lines in Fig. 5. For this particular example, we used nx=50.3f1.3f2+0.7f3, nx=60.7f+1.3f20.7f3, κx=0.1+0.01f+0.03f20.02f3, κy=0.05+0.01f+0.003f2. The simulated EOS waveforms are obtained with Jones-matrix calculations. The resulting Γ˜x,y are modulated by the standard configuration transfer functions Hx,y from Eq. (18) as shown in Fig. 4. We follow the procedures outlined above to extract nx,y and κx,y both by analytical inversion and numerical optimization overlaying the results in Fig. 5. The real part of the index ellipsoid is retrieved with great accuracy, only showing minor deviations from the actual values close to the zeros of the transfer function up to the upper limit of the spectral bandwidth of the input THz pulse. The imaginary part responsible for absorption within the sample is less accurate and shows greater deviations, due to the small values extracted by a logarithmic function [38]. The results are portrayed in Fig. 5 and show some significant advantages compared to the analytical extraction. The modeled refractive index is forced to be continuous and extrapolates well in the regions where the transfer functions are zero, especially when retrieving the imaginary part responsible for absorption.

    Simulated datasets with the actual temporal waveforms displayed on the top and noted with A and measured waveforms that are modulated by the transfer functions on the bottom noted with M.

    Figure 4.Simulated datasets with the actual temporal waveforms displayed on the top and noted with A and measured waveforms that are modulated by the transfer functions on the bottom noted with M.

    Simulated refractive index retrieval with real part on the left and imaginary part on the right. Ext. represents the analytical extraction and Fit is the numerically fitted values. The zeros of the transfer functions are presented in vertical dashed and dotted lines to highlight the frequencies where information cannot be accurately retrieved.

    Figure 5.Simulated refractive index retrieval with real part on the left and imaginary part on the right. Ext. represents the analytical extraction and Fit is the numerically fitted values. The zeros of the transfer functions are presented in vertical dashed and dotted lines to highlight the frequencies where information cannot be accurately retrieved.

    4. EXPERIMENTAL DEMONSTRATION

    To validate reconstruction and polarization detection techniques, we demonstrate single-shot EOS measurements using an optical-rectification-based THz source presented in Fig. 6 in which a Mg-doped LiNbO3 crystal is pumped with 40 fs 0.6 mJ 800 nm pulses with a spectral bandwidth Δλ=25  nm. The strong mismatch of refractive indices of IR and THz in Mg:LiNbO3 requires the use of tilted pulse fronts and a grating-lens optical system to achieve phase matching and efficient conversion from IR to THz [3941].

    THz source experimental setup. HWP and QWP refer to a half- or quarter-waveplate and Pel-BS is shortened for pellicle beamsplitter. The red path represents the IR pump path that is pulse front tilted by the grating and two lenses noted as “f” that can be adjusted in position when varying the pulse compression. The Mg:LiNbO3 crystal converts IR to THz that is then collimated by a Teflon lens and focused down to the ZnTe crystal by a TPX lens. The space in between the two lenses is used for a wiregrid polarizer and/or a sample. The orange path represents the probe path for EOS that is split off and recombined with the THz radiation by pellicle beamsplitters. After the THz-induced modulation, the probe is propagated through polarization optics before being split into horizontal and vertical polarization and then spectrally analyzed. A portion of the probe is split off for a reference line that is not modulated by THz.

    Figure 6.THz source experimental setup. HWP and QWP refer to a half- or quarter-waveplate and Pel-BS is shortened for pellicle beamsplitter. The red path represents the IR pump path that is pulse front tilted by the grating and two lenses noted as “f” that can be adjusted in position when varying the pulse compression. The Mg:LiNbO3 crystal converts IR to THz that is then collimated by a Teflon lens and focused down to the ZnTe crystal by a TPX lens. The space in between the two lenses is used for a wiregrid polarizer and/or a sample. The orange path represents the probe path for EOS that is split off and recombined with the THz radiation by pellicle beamsplitters. After the THz-induced modulation, the probe is propagated through polarization optics before being split into horizontal and vertical polarization and then spectrally analyzed. A portion of the probe is split off for a reference line that is not modulated by THz.

    In our setup, the compressor stage of the pump laser system is adjusted so that stretched pulses of 10–20 ps can be generated. To recover maximum efficiency in the THz optical rectification process, we utilize a two-lens optical system (f1=500  mm and f2=75  mm) from the grating to the LN crystal and adjust the second lens position away from the perfect imaging condition. HWPs before and after the grating are used to maximize the THz generation efficiency. One of the main benefits of this arrangement is that it is possible to perform stretched pulse single-shot EOS diagnostics at different pulse lengths (or chirp parameters C) with minimal variation (<10%) in conversion efficiency and without losing time synchronization between pump and probe.

    The generated THz is collected by a Teflon lens and refocused with a TPX lens down to the EO crystal. In between the two THz lenses we can insert a wiregrid polarizer or/and a birefringent sample. For EOS diagnostics, a small portion of the pump (8%) is split off by a pellicle beamsplitter, propagated through a delay stage to a second pellicle beamsplitter that overlaps THz and IR with minimal loss in THz transmission. The IR transmitted through the pellicle is used as a reference line, bypassing the THz-induced modulation in the EO crystal. The overlapped probe and THz are propagated through a 0.5 mm thick (110)-cut ZnTe crystal and then through a QWP and an HWP before being split into vertical and horizontal polarization with a Wollaston prism. A third pellicle beamsplitter is then used to re-overlap the reference line with the two split polarizations, which are then analyzed by a spectrometer. The spectrometer is home-built and consists of a 1200 lines/mm grating, a 100 mm lens, and a CMOS sensor with 1440×1080 resolution and 3.45  μm×3.45  μm pixel size. It is set up in a 1f1f configuration, which means that the distance between the grating and the lens, as well as the lens and the CMOS, is 100 mm. This results in collimated rays in x and focused in y, forming three spectral lineouts on the CMOS sensor. Two of them are the modulated signals |E˜1,s(ω)|2 and |E˜2,s(ω)|2, while the third is from the reference line.

    To measure Y1,Y2, we must obtain |T1E˜0(ω)|2, |T2E˜0(ω)|2 for each shot. This is performed by measuring all three spectral lines without a THz signal and determining the intensity ratio between the two split polarizations and the reference line. |T1,2E˜0(ω)|2 are then obtained for each individual shot from the calibrated reference line, allowing us to determine Y1,Y2 using Eq. (13).

    The time calibration as well as the spectral-temporal calibration factor C is obtained by adjusting the delay stage by a known amount and measuring the signal position on the spectrometer. A linear fit of the signal position with respect to time then yields the time calibration as tpx=(2.25±0.04)  fspixel. Higher order terms are also considered in the fit but contribute around 1% over the range of the spectrometer and are therefore neglected. With the known spectral bandwidth of our laser pulse, we obtain the chirp parameter as C=ωpx(tpx)1=3.5  ps2. These parameters are important as they inform the transfer functions to be used in the reconstruction of the temporal profiles.

    We first investigate the case in which the THz polarization is reoriented by a wiregrid polarizer. The source emits THz polarized in y^. The detected polarization state transmitted through the polarizer at angle θ in respect to the y-axis is described as ETHzθ(t)=(sinθcosθcos2θ)ETHz(t).

    The wiregrid is then rotated from 0° to 90° in increments of 15° and measurements are performed at each polarizer angle with the captured single-shot temporal and spectral profiles for Γx and Γy presented in Fig. 3. Although these temporal profiles are representative of the THz fields Ex and Ey, it is important to note that the signal is modulated spectrally by the transfer functions, as shown in the spectrum in Fig. 3.

    Given that the THz polarization that illuminates the crystal is known, we can reconstruct the temporal profile accurately, as discussed earlier. In the case of the wiregrid polarizer, the relation between Γx and Γy is given as r=tanθ. Using the transfer functions in Eqs. (22) and (23) we can then reconstruct the temporal profile with Eq. (16). The reconstructed temporal profile and spectrum are presented in Fig. 3 up to 45° with a waveform captured and reconstructed with DEOS for reference.

    To demonstrate THz ellipsometry, we place a y-cut LiNbO3 crystal between the wiregrid polarizer and the TPX lens. This material has shown strong birefringence in the THz range and is well characterized, making it a well-suited sample to demonstrate single-shot-based ellipsometry [42]. First, we perform two separate measurements in which the sample’s ordinary and extraordinary axes are parallel to the THz polarization and the transmitted THz waveform is captured with DEOS-based reconstruction methods. This provides us with a benchmark measurement of the sample index ellipsoid without distortions from the transfer functions and establishes reference values to be compared to and noted as Ref. in Fig. 8.

    Measured and reconstructed signals for different wiregrid polarizer orientations. In the top left are the temporal profiles under the standard configuration polarization readout. Γx and Γy refer to the measured polarization components that are modulated by the transfer functions. Γ refers to the reconstructed waveforms with a DEOS-based measurement and reconstruction as comparison. We note that for the 0° polarization angle in the standard configuration, no accurate reconstruction can be performed as the Γx component is not captured. The top right shows the spectral amplitudes of the individual polarizations for different wiregrid orientations with the transfer functions highlighting the loss of information. The bottom left shows the peak amplitude of Γx and Γy with theoretical expectations for different wiregrid orientations. The bottom right shows the reconstructed spectra with a DEOS dataset as a comparison. The change in amplitude here is due to the filtering effect of the wiregrid polarizer when it is oriented at non-zero angles in respect to the y-axis.

    Figure 7.Measured and reconstructed signals for different wiregrid polarizer orientations. In the top left are the temporal profiles under the standard configuration polarization readout. Γx and Γy refer to the measured polarization components that are modulated by the transfer functions. Γ refers to the reconstructed waveforms with a DEOS-based measurement and reconstruction as comparison. We note that for the 0° polarization angle in the standard configuration, no accurate reconstruction can be performed as the Γx component is not captured. The top right shows the spectral amplitudes of the individual polarizations for different wiregrid orientations with the transfer functions highlighting the loss of information. The bottom left shows the peak amplitude of Γx and Γy with theoretical expectations for different wiregrid orientations. The bottom right shows the reconstructed spectra with a DEOS dataset as a comparison. The change in amplitude here is due to the filtering effect of the wiregrid polarizer when it is oriented at non-zero angles in respect to the y-axis.

    Real part of the ordinary and extraordinary refractive indices no and ne from y-cut LiNbO3. Three different types of measurements are portrayed: Ref. denotes the refractive indices retrieved through two separate measurements by orienting the main axes of the sample to the THz polarization and then fit the index ellipsoid; Ext. refers to the single-measurement technique in which the sample is oriented 45° to the THz polarization, allowing both refractive indices to be sampled simultaneously with the cost of losing information at the zeros of the transfer function; Fit refers to the same setup as Ext. but uses a fitting algorithm to reconstruct the index ellipsoid. The vertical dotted lines are the zeros of the transfer function to visualize where information cannot be analytically retrieved under the single-measurement configuration. The spectral amplitude of the THz signal without the sample is shown in the bottom in green to show how refractive index extraction errors are correlated with low spectral amplitude.

    Figure 8.Real part of the ordinary and extraordinary refractive indices no and ne from y-cut LiNbO3. Three different types of measurements are portrayed: Ref. denotes the refractive indices retrieved through two separate measurements by orienting the main axes of the sample to the THz polarization and then fit the index ellipsoid; Ext. refers to the single-measurement technique in which the sample is oriented 45° to the THz polarization, allowing both refractive indices to be sampled simultaneously with the cost of losing information at the zeros of the transfer function; Fit refers to the same setup as Ext. but uses a fitting algorithm to reconstruct the index ellipsoid. The vertical dotted lines are the zeros of the transfer function to visualize where information cannot be analytically retrieved under the single-measurement configuration. The spectral amplitude of the THz signal without the sample is shown in the bottom in green to show how refractive index extraction errors are correlated with low spectral amplitude.

    The main axes of the sample are then oriented 45° to the incident THz polarization, and the refractive indices are retrieved through analytical extraction and numerical fit. Both are in good agreement with the reference and literature values [43]; however, the analytical extraction deviates slightly at the zeros of the transfer functions. The numerical model is constrained to a third-order polynomial, and therefore does not present the same deviations and is in better agreement with the reference. At frequencies where the spectral amplitude is small, both the analytical extraction and the numerical fit stray more from the reference values. This demonstrates that in a single-shot measurement, information on the entire index ellipsoid can be retrieved over a broad spectral bandwidth with minor limitations due to the transfer functions. This allows THz ellipsometry measurements to be performed on samples that vary in time or where the THz source has significant shot-to-shot fluctuations or low repetition rates.

    5. CONCLUSION

    In summary, in this paper we discuss effects that arise for arbitrary THz polarization in single-shot, spectral-temporal encoding EOS schemes. The transfer-function-based formalism established in DEOS is modified to consider an additional polarization component, allowing complex polarized THz waveforms to be captured and characterized. We highlight the importance of knowing the THz polarization state when performing DEOS with arbitrary polarization and present the necessary modifications to the reconstruction method to avoid inaccurate results. In this framework, we have presented methods to extract THz polarization information over a large spectral range, allowing THz ellipsometry to be performed on an unknown sample. These methods are first demonstrated numerically and then experimentally supported with a table-top optical-rectification-based THz source.

    Using the THz polarization sensitive single-shot methods introduced in this work, we envision the possibility to analyze and study samples that exhibit dynamically changing birefringence in the THz range or exhibit poor repeatability in a single shot and without the need to mechanically intervene. Furthermore, THz sources with significant shot-to-shot jitter can be characterized not only temporally and spectrally but also in their time-dependent polarization state, which could be of particular importance for the characterization of electron-beam-based THz sources such as THz-FELs using helical undulator geometries.

    APPENDIX A

    Transfer Function Derivation

    In this section we dive into more detail on deriving the transfer functions. Starting from the definition of Y(ω), let us elaborate on the steps performed in Eq. (12).

    With E˜s(ω)=E˜0(ω)*T(1+meiΨΓ˜(ω)) we can rewrite E˜s(ω) up to first order in Γ˜(ω) as |E˜s(ω)|2=|T|2|E˜0(ω)|2+|T|22e((E˜0(ω)*Γ˜(ω))meiΨE˜0*(ω))+O(Γ˜2),leading to Y(ω) expression of Eq. (12). However, to fully understand the appearance of the transfer functions, let us expand the second summand: e((E˜0(ω)*Γ˜(ω))meiψE˜0*(ω))=e(2πCA*(ωω0C)exp(i(ωω0)22C)meiψ×Γ˜(Ω)A(ωΩω0C)exp(i(ωΩω0)22C)dΩ)=e(2πC|A(ωω0C)|2meiψ×Γ˜(Ω)exp(iΩ22C)exp(iΩ(ωω0)C)dΩ)=2πC|A(ωω0C)|2me(F(Γ˜(Ω)exp(i(Ω22C+ψ)))).

    We assume A(ωω0C) is slowly varying, allowing us to pull it out of the integral. F represents the Fourier transform from Ω to t. After background subtraction and normalization, we then have Y(ω)=2me(F(Γ˜(Ω)exp(i(Ω22C+ψ)))).

    To further evaluate this expression, note e(F(f(ω)))=(F(f(ω))+F(f*(ω)))/2. And given that Γ˜(Ω) is a Hermitian function this allows to express Y(ω) as Y(ω)=F(Γ˜(Ω)2mcos(Ω22C+ψ)).

    THz Polarimetry

    In this section we give more details on the derivation of the modulation functions for single-shot THz polarimetry. Our goal is to obtain an expression for the modulated signal field of the form Es,i=TiE0(1+mixeiψixΓx+miyeiψiyΓy),so that using the same steps of the derivation leading to Eq. (A4) we can express the Fourier transform of the measured waveforms as Y˜i=2mixcos(Ω22C+ψix)Γ˜x+2miycos(Ω22C+ψiy)Γ˜y.

    We start with the Jones calculus from Eq. (4) adding an additional matrix element that describes the modulation caused by the orthogonal polarization THz field EyTHzΓy. This results in a total transport matrix: M=W(π,θH)W(π/2,θQ)R(Ψ(αy))MXtal(αy)R(Ψ(αy))R(Ψ(α))MXtal(α)R(Ψ(α)),where αy=α+π/2 and θH, θQ are the angles of the HWP and QWP fast axes with respect to the [110] axis of the EO crystal. In the standard configuration α=0, θH=0, θQ=π/4, θIR=0 and the resulting electric field modulation is expressed as Esstd=ME0(10)=E02(exp(iΓx2)(exp(iΓy2)iexp(iΓy2))exp(iΓx2)(exp(iΓy2)+iexp(iΓy2)))E02(1Γx2+iΓy2i+iΓx2Γy2)+O(Γ2)after using small angle approximation.

    From this expression we can read the transfer function coefficients directly as T1=1/2, T2=i/2, m1x=1/2, m2x=1/2, m1y=m2y=1/2, ψ1x=ψ2x=0, and ψ1y=ψ2y=π/2, which are then put into Eq. (A6), leading directly to Eq. (18). In the case of the DEOS configuration, θIR=π/4α=π/2, θH=π/8, and θQ=0, an additional polarization that leads to a phase shift Γy results in the modulation EsPD=E022(exp(iΓy2)(exp(iΓx2)+1)+iexp(iΓy2)(exp(iΓx2)1)iexp(iΓy2)(exp(iΓx2)+1)+exp(iΓy2)(exp(iΓx2)1))E02(1(1i)Γx4+iΓy2i(1+(1+i)Γx4+iΓy2))+O(Γ2),from which we can read off T1=1/2, T2=i/2, m1x=1/(22), m2x=1/(22), m1y=m2y=1/2, ψ1x=π/4, ψ2x=π/4, and ψ1y=ψ2y=π/2, which leads to Eq. (20).

    References

    [1] M. Koch, D. M. Mittleman, J. Ornik. Terahertz time-domain spectroscopy. Nat. Rev. Methods Primers, 3, 48(2023).

    [2] J. van Tilborg, C. Schroeder, C. Filip. Terahertz radiation as a bunch diagnostic for laser-wakefield-accelerated electron bunches. Phys. Plasma, 13, 056704(2006).

    [3] D. Mittleman, M. Gupta, R. Neelamani. Recent advances in terahertz imaging. Appl. Phys. B, 68, 1085-1094(1999).

    [4] T. L. Cocker, V. Jelic, M. Gupta. An ultrafast terahertz scanning tunnelling microscope. Nat. Photonics, 7, 620-625(2013).

    [5] X.-C. Zhang. Free-space electro-optic sampling of terahertz beams. Appl. Phys. Lett., 67, 3523-3525(1995).

    [6] A. Nahata, A. S. Weling, T. F. Heinz. A wideband coherent terahertz spectroscopy system using optical rectification and electro-optic sampling. Appl. Phys. Lett., 69, 2321-2323(1996).

    [7] G. Gallot, D. Grischkowsky. Electro-optic detection of terahertz radiation. J. Opt. Soc. Am. B, 16, 1204-1212(1999).

    [8] A. Nahata, D. H. Auston, T. F. Heinz. Coherent detection of freely propagating terahertz radiation by electro-optic sampling. Appl. Phys. Lett., 68, 150-152(1996).

    [9] M. Lenz, A. Fisher, A. Ody. Electro-optic sampling based characterization of broad-band high efficiency THz-FEL. Opt. Express, 30, 33804-33816(2022).

    [10] A. Fisher, M. Lenz, A. Ody. Towards higher frequencies in a compact prebunched waveguide THz-FEL. Nat. Commun., 15, 7582(2024).

    [11] M. C. Hoffmann, S. Schulz, S. Wesch. Coherent single-cycle pulses with MV/cm field strengths from a relativistic transition radiation light source. Opt. Lett., 36, 4473-4475(2011).

    [12] R. Pompili, A. Cianchi, D. Alesini. First single-shot and non-intercepting longitudinal bunch diagnostics for comb-like beam by means of electro-optic sampling. Nucl. Instrum. Methods Phys. Res. Sect. A, 740, 216-221(2014).

    [13] C. Scoby, P. Musumeci, J. Moody. Electro-optic sampling at 90 degree interaction geometry for time-of-arrival stamping of ultrafast relativistic electron diffraction. Phys. Rev. Spec. Top. Accel Beams, 13, 022801(2010).

    [14] M. Ota, K. Kan, S. Komada. Longitudinal and transverse spatial beam profile measurement of relativistic electron bunch by electro-optic sampling. Appl. Phys. Express, 14, 026503(2021).

    [15] G. Berden, S. Jamison, E.-A. Knabbe. Benchmarking of electro-optic monitors for femtosecond electron bunches. Phys. Rev. Lett., 99, 164801(2007).

    [16] H. Hamster, A. Sullivan, S. Gordon. Subpicosecond, electromagnetic pulses from intense laser-plasma interaction. Phys. Rev. Lett., 71, 2725-2728(1993).

    [17] K. Y. Kim, B. Yellampalle, J. Glownia. Measurements of terahertz electrical conductivity of intense laser-heated dense aluminum plasmas. Phys. Rev. Lett., 100, 135002(2008).

    [18] G. T. Noe, I. Katayama, F. Katsutani. Single-shot terahertz time-domain spectroscopy in pulsed high magnetic fields. Opt. Express, 24, 30328-30337(2016).

    [19] S. W. Teitelbaum, B. K. Ofori-Okai, Y.-H. Cheng. Dynamics of a persistent insulator-to-metal transition in strained manganite films. Phys. Rev. Lett., 123, 267201(2019).

    [20] S. M. Teo, B. K. Ofori-Okai, C. A. Werley. Invited article: single-shot THz detection techniques optimized for multidimensional THz spectroscopy. Rev. Sci. Instrum., 86, 051301(2015).

    [21] Z. Zhang, F. Y. Gao, J. B. Curtis. Terahertz field-induced nonlinear coupling of two magnon modes in an antiferromagnet. Nat. Phys., 20, 801-806(2024).

    [22] J. Zhang, X. Tan, M. Liu. Cooperative photoinduced metastable phase control in strained manganite films. Nat. Mater., 15, 956-960(2016).

    [23] Z. Jiang, X.-C. Zhang. Electro-optic measurement of THz field pulses with a chirped optical beam. Appl. Phys. Lett., 72, 1945-1947(1998).

    [24] I. Wilke, A. M. MacLeod, W. A. Gillespie. Single-shot electron-beam bunch length measurements. Phys. Rev. Lett., 88, 124801(2002).

    [25] B. Yellampalle, K. Kim, G. Rodriguez. Algorithm for high-resolution single-shot THz measurement using in-line spectral interferometry with chirped pulses. Appl. Phys. Lett., 87, 211109(2005).

    [26] F. Sun, Z. Jiang, X.-C. Zhang. Analysis of terahertz pulse measurement with a chirped probe beam. Appl. Phys. Lett., 73, 2233-2235(1998).

    [27] E. Roussel, C. Szwaj, C. Evain. Phase diversity electro-optic sampling: a new approach to single-shot terahertz waveform recording. Light Sci. Appl., 11(2022).

    [28] N. Nemoto, T. Higuchi, N. Kanda. Highly precise and accurate terahertz polarization measurements based on electro-optic sampling with polarization modulation of probe pulses. Opt. Express, 22, 17915-17929(2014).

    [29] N. Yasumatsu, S. Watanabe. Precise real-time polarization measurement of terahertz electromagnetic waves by a spinning electro-optic sensor. Rev. Sci. Instrum., 83, 023104(2012).

    [30] K. Oka, T. Kato. Spectroscopic polarimetry with a channeled spectrum. Opt. Lett., 24, 1475-1477(1999).

    [31] Q. Lin, X. Yuan, X. Zeng. Single-shot terahertz polarization detection based on terahertz time-domain spectroscopy. Photonics Res., 10, 1374-1379(2022).

    [32] Q. Chen, M. Tani, Z. Jiang. Electro-optic transceivers for terahertz-wave applications. J. Opt. Soc. Am. B, 18, 823-831(2001).

    [33] E. Collett. Field Guide to Polarization(2005).

    [34] J. Van Tilborg, C. Schroeder, C. Filip. Temporal characterization of femtosecond laser-plasma-accelerated electron bunches using terahertz radiation. Phys. Rev. Lett., 96, 014801(2006).

    [35] H. Bakker, G. Cho, H. Kurz. Distortion of terahertz pulses in electro-optic sampling. J. Opt. Soc. Am. B, 15, 1795-1801(1998).

    [36] F. Coppinger, A. Bhushan, B. Jalali. Photonic time stretch and its application to analog-to-digital conversion. IEEE Trans. Microw. Theory Tech., 47, 1309-1314(1999).

    [37] J.-L. Coutaz, F. Garet, V. Wallace. Principles of Terahertz Time-domain Spectroscopy(2018).

    [38] K. Cheung, D. Auston. A novel technique for measuring far-infrared absorption and dispersion. Infrared Phys., 26, 23-27(1986).

    [39] J. Hebling, G. Almasi, I. Z. Kozma. Velocity matching by pulse front tilting for large-area THz-pulse generation. Opt. Express, 10, 1161-1166(2002).

    [40] J. Fülöp, L. Pálfalvi, G. Almási. Design of high-energy terahertz sources based on optical rectification. Opt. Express, 18, 12311-12327(2010).

    [41] T. Kroh, T. Rohwer, D. Zhang. Parameter sensitivities in tilted-pulse-front based terahertz setups and their implications for high-energy terahertz source design and optimization. Opt. Express, 30, 24186-24206(2022).

    [42] Y. M. Sun, Z. L. Mao, B. H. Hou. Giant birefringence of lithium niobate crystals in the terahertz region. Chin. Phys. Lett., 24, 414(2007).

    [43] X. Wu, C. Zhou, W. R. Huang. Temperature dependent refractive index and absorption coefficient of congruent lithium niobate crystals in the terahertz range. Opt. Express, 23, 29729-29737(2015).