1. Introduction
The number of applications that require intense terahertz (THz) pulses has been rapidly increasing in recent years [
1,
2]. Optical rectification (OR) of the femtosecond laser pulses in lithium niobate (LN) crystals has emerged as the most powerful way to generate broadband, nearly-single-cycle THz pulses with spectra centered at frequencies below 1 THz. Single-cycle THz pulses with record energy of 0.4 mJ [
3] and pump-to-THz conversion efficiency η = 3.8% [
4] have been reported. To satisfy phase matching condition, the pump beam is reflected off a grating to acquire a tilted-pulse-front (TPF), which is then subsequently imaged onto the crystal using a lens or a telescope. However, imaging errors limit the effective length of the crystal [
5,
6]; therefore, further increase in THz generation performance is challenging. Attempts to mitigate this problem by using contacting grating [
7], stair-step echelons [
8,
9], and a combination of stair-step echelon with reflecting grating [
10] still do not result in THz generation with energies or efficiencies close to the abovementioned values. It is well established [
11,
12] that the effective interaction length in semiconductor materials such as ZnTe, GaP, and GaAs can be significantly larger due to operation ability at smaller tilt angles < 30° and lower THz-wave absorption at room temperature. However, to escape 2nd- and 3rd-order multiphoton pump absorption (MPA), these materials have to be pumped at longer wavelengths
λp ≥ 1.7 μm, where it is still relatively challenging to obtain femtosecond laser pulses with the required high power.
It should be noted that MgO-doped LN crystals have the advantages of a large nonlinear coefficient (especially among the inorganic crystals), sufficiently wide bandgap (the lowest order effective MPA is four-photon absorption at the 1.03 µm pump wavelength), high damage threshold and, finally, the ability for spatial sign-modulation of the nonlinear coefficient. The latter is important, as nowadays, high-power pulses can be used to pump periodically-poled lithium niobate (PPLN) crystals. The large area PPLN crystals with a cross section 3 mm × 3 mm are commercially available, and even a crystal with dimensions 1 cm × 1.5 cm was recently used for THz generation [
13]. The ability to significantly reduce THz-waves absorption in PPLN crystal by cryogenic cooling makes it possible to use crystals with a length of about 2–3 cm [
13,
14]. Therefore, investigations of alternative methods for nearly single-cycle broadband THz pulse generation in domain-engineered PPLN crystals are interesting.
It has been well established that the use of aperiodically-poled lithium niobate (APPLN) leads to the generation of broadband THz pulses [
15,
16,
17,
18]. However, its temporal form is not a nearly single-cycle THz pulse (i.e., it is not determined by the spectral bandwidth) because different spectral components are radiated from different positions of the crystal. Like in the optical region, the technique of broadband pulse compression using Bragg grating, chirped mirror, and pairs of gratings or prisms, can be used. However, strong THz beam divergence limits its applications.
The opportunity of nearly single-cycle THz pulse generation in APPLN crystal in conjunction with a THz chirped mirror has been analyzed in detail by Yahaghi and co-worker in Ref. [
19]. However, the difficulties of THz chirped mirror fabrication and its application in THz generator schemes complicates practical implementation. To avoid this problem, a nearly single-cycle THz pulse generation in chirped APPLN crystal using a pair of chirped optical pump pulses with various chirp rates was recently proposed and analyzed [
20]. The relative chirp between the pump pulses compensates the temporal shift between THz frequencies generated at different locations in the crystal. As a consequence, the THz pulse emerges compressed upon exiting the crystal.
Taking advantage of the above idea, here, we introduce a new scheme to generate THz pulses with controllable number of THz field oscillations, from nearly single- to multi-cycles and, correspondingly, on the spectral bandwidth from broadband to narrowband. In contrast to the laser pulse format used in [
20], here, the sequence of transform-limited laser pulses with chirped delay between adjacent pulses (so-called pulse-position modulated (PPM) signal in optical communication) is considered as the pump wave. Such sequences of laser pulses can be obtained by using a birefringent crystal array [
21] or modern pulse shaping techniques [
22].
The operation principle is similar to that presented in [
20]; it is physically based on the fact that the temporary overlap of a sufficiently large number of phase-locked spectral components leads to the formation of an ultrashort pulse. It is well-known [
16,
17,
18] that in the case of APPLN crystal pumping by single transform-limited laser pulse the generated THz waveform corresponds to domain structure of the crystal. More exactly, it corresponds to spatial reversal replica (or so-called phase-conjugate replica) of the APPLN domain lengths distribution, considering that the exciting laser pulse propagates faster than the radiated THz-wave. The set of these temporally-shifted replicas is originated when the pump wave is the sequence of laser pulses with delays Δ
tm (
m = 1, 2, 3…
M). By appropriately choosing the Δ
tm, it is possible to obtain a temporary overlap of all generated frequencies (in the range of Δ
ωTHz) that results in the formation of a nearly single-cycle THz pulse in the case of sufficiently large
M and Δ
ωTHz, where Δ
ωTHz is the bandwidth of THz generation in APPLN by separate laser pulse.
This paper provides a simple theoretical model for the efficient generation of a nearly single-cycle THz pulse in the chirped APPLN crystal pumped by a sequence of laser pulses with chirped delay times. Closed-form expression for the spectrum and waveform of the generated THz pulse is obtained. A peak THz electric field strength of 0.3 MV/cm is predicted for the APPLN crystal pumped by the sequence of laser pulses with peak intensity of the separate pulse in the sequence about 20 GW/cm2. By focusing the THz beam and using higher pump power, the electric field strength can be increased by an order of magnitude.
2. Theoretical Model
Let us consider a transform-limited optical pulse with Gaussian temporal and spatial profiles that propagates along the
x-axis of APPLN crystal having a nonuniform distribution of the nonlinear coefficient
d33 =
d33(
x). The APPLN crystal consists of
N number domains having lengths of
lk (
k = 1, 2…
N), which are reduced along the direction of the pump beam propagation. In every adjacent domain the sign of the nonlinear coefficient
d33 is opposite (
Figure 1).
To calculate the electric field of THz generation, the radiating antenna model [
23,
24] is used. For simplicity, pump pulse depletion and its spatial and temporal distortions during propagation in the crystal are neglected. To avoid optical and THz wave reflections at the interfaces of the crystal, we assume that it is incorporated in a linear medium having the same refractive index. With these assumptions, the instantaneous intensity of the pump pulse is given by
where
a is the beam waist radius,
ng is the group refractive index,
c is the velocity of light, 2
σ is the pulse duration at the 1
/e level that is related to the full width at half maximum by
σFWHM = 2(ln2)
1/2σ, and
I0 is the peak on-axis intensity at entrance of the crystal
x = 0, which is related to pulse energy ℑ by
I0 = ℑ/π
3/2a2σ.
It is known [
25] that in the frequency domain low-frequency nonlinear polarization
PNL is related to the pump intensity by the following expression:
where
,
is the Fourier transform of the on-axis intensity
at
x = 0,
n0 is the refractive index at the laser frequency, and
kg =
ωng/
c.
To calculate the THz far-field produced by the nonlinear polarization, the APPLN crystal is considered as an antenna fed by the
z-oriented current given by
In antenna theory, the calculation of the radiation field for a given distribution of excitation currents leads to the integration in the so-called radiation integral over spatial coordinates
[
26]. To simplify the integration over transverse coordinates
y and
z, it is assumed that the pump beam is completely inside the crystal. Finally, by presenting THz field as the superposition of fields radiated by the separate crystal domains, we determined that the THz field is given by
where
µ0 is the permeability of vacuum,
S =
πa2 is the pump beam area,
lk is the length of separate
k-domain of the APPLN crystal,
nTHz is the refractive index of the crystal at THz frequencies, Δ
n =
nTHz −
ng is the refractive indexes mismatch,
L is the length of the crystal, and
R is the distance between the point of THz field observation and the exit surface of the APPLN crystal.
From the above, it follows that the temporal form of the radiated THz pulse can be easily calculated. Indeed, the factor
iωI(
ω) is the Fourier transform of temporal derivative of the intensity
the exponential factors in the sum represent temporal shifts of
on the corresponding values, and the exponential factor before the summation symbol is related to the time of THz-wave propagation from the entrance surface of the crystal to a point of THz field observation. Hence, the THz electric field in the time domain is given by
where
tr =
t −
t0 is the retarded time.
The above equation can be easily modified for the case, when APPLN crystal is pumped by a sequence of the laser pulses with delays Δ
tm =
τ2m-1 −
τ1 and
m = 1, 2, 3…
N/2. Using the superposition principle, the generated THz field can be expressed by
Thus, the THz pulse form and spectrum can be easily calculated for a given distribution of the domain length lk and pump pulse duration σ.
3. Results and Discussions
Let’s consider linearly-chirped APPLN, in which the domain lengths are decreasing from the front to the end of the crystal as lk = [l1 − (k − 1)δ], where δ = lk − lk+1 is the length difference of neighboring domains and k = 1, 2, …32. Note that central frequency radiated from the k-domain is inversely proportional to its length lk and, therefore, higher THz frequencies are generated closer to the exit surface of the crystal. Because THz absorption at higher frequencies is larger, such a design of the crystal is favorable to mitigating THz-wave damping.
Using Equation (5), the calculated THz waveform in the crystal with
l1 = 90 µm,
δ = 1.6 µm, and
σ = 0.4 ps is shown in
Figure 2. Also, in the same figure, we present the THz waveforms obtained for pump pulses delayed at Δ
tm =
τ2m-1 −
τ1 for
m = 2 and
m = 3, respectively. It is seen that there is a temporary overlap of the fields generated by the first positive domain
k = 1 and the third (the fifth) domain, if they are excited with delays Δ
t2 (Δ
t3), respectively.
Therefore, temporary overlap of the fields generated by all positive domains will occur if the APPLN crystal is pumped by a sequence of laser pulses having delays Δtm with m = 1, 2, 3…N/2. As a result of the fields overlapping, a nearly single-cycle THz pulse can be formed, if the number of domains N in the crystal and value of δ (characterizing the domain’s length gradient) are sufficiently large. These conditions are necessary to obtain a broadband THz generation in APPLN crystal pumped by the single laser pulse. The use of a sequence of laser pulses with suitable delays Δtm is needed to satisfy the certain phase relationship for the spectral components of the broadband radiation. Obviously, in the case of relatively narrowband radiation (or small value of δ), the temporary overlap of THz fields will lead to the formation of a few cycle THz pulses, with the number of the cycles depending on δ.
The above statements are illustrated in
Figure 3a, where using Equation (6), we present the temporal forms of THz pulse
EΣ(
tr), excited by the sequence of pulses having delays Δ
tm (
m = 1, 2, 3…16) for cases of
δ = 1.6 µm,
δ = 0.8 µm, and
δ = 0.5 µm, respectively. The required sequence of laser pulses is schematically illustrated in
Figure 3b for cases of
δ = 1.6 µm and
δ = 0.8 µm, respectively.
For many applications, especially for THz-driven particle acceleration, the low-frequency part of the THz spectrum (0.3 to 1.5 THz) is most useful, because the longer wavelength compares well with typical sizes of particle bunches [
10]. In our case, the THz spectrum is given by the product of the spectrum
E(
ω) generated by the first pump pulse and the factor
F = [1 +
exp(−
iωΔ
t2) +
exp(−
iωΔ
t3) + … +
exp(−
iωΔ
tM)] representing application of the sequence of pump pulses with the delays Δ
tm. The dependences of spectral amplitude
EΣ(
f) on frequency in the cases of
δ = 1.6 µm and
δ = 0.8 µm are presented in
Figure 4. It is seen that the spectrum of the nearly single-cycle THz pulse (
δ = 1.6 µm) is broadband and centered mainly at frequencies lower than 1 THz. The narrowband multicycle THz radiation corresponds to the
δ = 0.8 µm.
It should be noted that all graphs presented in
Figure 2,
Figure 3 and
Figure 4 are built with constant factor accuracy, i.e., arbitrary units are used for the values of functions. However, the absolute value of the field strength is very important parameter. By including the numerical value of
A0 in the calculations using Equations (5) and (6), we determine that the maximum electric field strength
EΣmax is estimated to be about 300 kV/cm at the distance
R = 20 mm, if the beam area
S and the peak intensity
I0 of every pulse in the sequence are 7 mm
2 and 20 GW/cm
2, respectively. Obviously, the estimated value can be increased by enhancing the pump power and by focusing the THz beam with a short focus lens. The choice of the numerical value for
I0 is related to the established fact that the generated THz field strength linearly depends on the pump intensity for up to ~20 GW/cm
2 [
3]. As for the distance
R = 20 mm, it satisfies the well-known condition of the far-field approximation [
26].
It is useful to note that there is a physically more transparent way to estimate peak electric field strength. It was already mentioned that THz pulses radiated by the first domain
E1(
tr) and by all other positive domains are temporarily overlapped if the sequence of laser pulses with delays Δ
tm =
τ2m-1 −
τ1 (
m = 2, 3…
N/2) is used. Therefore, the generated electric field can be approximately presented as
EΣ(
tr) ≈
NE1(
tr)/2. Using expansion into a Taylor series, the first term of Equation (6) can be rewritten as
where
is the second order derivative of the pump intensity.
Using above approximation, the maximum field strength can be given by EΣmax ≈ B0NI0, where B0 = µ0d33Sl1/πσ2Rn0c. By substituting the numerical values used in the previous field estimation, we obtain EΣmax ≈ 320 kV/cm, which is in good agreement with the earlier obtained result.