Abstract
The generation of terahertz waves from air-plasma induced by femtosecond three-color harmonic pulses with a frequency ratio of 1:2:m (m is a positive integer), based on the transient photocurrent model and the sawtooth-like electric field formed via multi-color harmonic pulses superposition, has been theoretically investigated. It can be seen that when the air is saturated ionized and the electron density reaches the same maximum, for the same number of harmonic pulses, terahertz conversion efficiency is not always higher when the electric field shape in the composed pulse envelope is closer to a sawtooth waveform and more asymmetric. Besides, the specific wavelength combination schemes of femtosecond three-color harmonic pulses with the frequency ratios of 1:2:3 and 1:2:4 have also been simulated, which can significantly enhance the generation of terahertz waves, and are realized by adding only a set of optical parametric amplifiers on the basis of the conventional two-color laser pulse case at the frequency ratio of 1:2. Our study will be helpful to obtain intense terahertz sources and provide guidance for experimental operations.
Terahertz (THz) waves have attracted significant attention in the fields of physics, chemistry, physiology and medicine
In the paper, based on the sawtooth-like electric field formed via femtosecond multi-color harmonic pulse superposition and the conventional scheme that femtosecond two-color laser pulses at the frequency ratio of 1:2 excite the air-plasma to radiate THz waves experimentally, we theoretically investigated the variation trends of THz wave generation by femtosecond three-color harmonic pulses with the frequency ratio of 1:2:m (3≤m≤20) according to the transient photocurrent model (TPCM), and distinguished the similarity between the electric field shape of the composed pulse and the perfect sawtooth waveform by a simple and efficient method. The results reveal when the electric field shape in the composed pulse envelope is closer to the sawtooth waveform and more asymmetric, the THz conversion efficiency is not always higher for the same number of harmonic pulses. The specific wavelength combination schemes of three-color harmonic pulses at the frequency ratios of 1:2:3 and 1:2:4 have been also studied and presented. These specific wavelength combination schemes can significantly enhance the THz wave generation, and be achieved by adding only a set of optical parametric amplifiers (OPA) on the basis of the conventional two-color laser pulse case at frequency ratio of 1:2. In addition, we demonstrate when the relative phase of the highest harmonic is varied, the maximum THz yield of the three-color pulse scheme at the frequency combination of 1:2:4 is slightly smaller than that of the three-color pulse case at the frequency combination of 1:2:3, but its overall THz yield can be maintained better stability and smaller fluctuation. When the amplitude ratio of the highest harmonic is appropriately increased, the three-color pulse scheme at the frequency combination of 1:2:4 can generate stronger THz waves than that of the three-color pulse case at the frequency combination of 1:2:3.
As a theoretical mechanism of THz wave generation in air-plasma excited via femtosecond two-color laser pulses, the TPCM can also be applied to the analysis of multi-color harmonic pulse schemes
, | (1) |
where is the number of multi-color harmonic pulses, is the electric field amplitude, is the central frequency of each pulse, is the pulse duration. Then, according to the TPCM, the composed pulse ionizes air to form an air-plasma. The type of ionization can be classified by the Keldysh parameter
, | (2) |
where is the electron charge, is the frequency of the laser, is the electron mass, is the laser pulse induced ponderomotive potential, is the ionization potential of the gas. When , the main ionization mechanism is tunnel ionization and the ionization rate can be expressed as
, | (3) |
where is the atomic frequency unit, is the atomic unit of the electric field. By tunnel ionization, bound electrons move away from the neutral gas particle as the Coulomb barrier is suppressed by the instantaneous electric field of the composed pulse, a large number of free electrons will be produced, and the time-varying electron density is given by
, | (4) |
where is the electron density, is the neutral gas density. Then, the free electrons will drift under the electric field of the composed pulse. The drift velocity is
, | (5) |
where is the time to produce free electrons. The drifting motion of free electrons will generate a photocurrent, and the time-varying photocurrent density is presented as
, | (6) |
where is the collision frequency of electron-ion and electron-neutral particle. The pulse duration is much shorter than the average collision time between particles, so the collision frequency can be ignored. Finally, the time-varying photocurrent density generates THz waves. The electric field amplitude of THz waves can be expressed as
. | (7) |
According to
Thus, we calculate the electric field peak-peak amplitude of THz pulse radiated from femtosecond three-color harmonic pulses with the frequency ratio of 1:2:m (3≤m≤20), and compare them with that generated by the two-color pulse case, as shown in

Fig. 1 The THz electric field peak-peak amplitude of three-color harmonic pulse schemes with the frequency ratio of 1:2:m (3≤m≤20), when m=2, it indicates the two-color pulse case at the frequency ratio of 1:2, red dot lines connecting data (blue) are used to guide vision only and do not mean anything
图1 频率比为1:2:m(3≤m≤20)的三色谐波脉冲方案产生的太赫兹波电场峰峰值,当m=2时,表示频率比为1:2的双色脉冲方案,连接数据点(蓝色)的红色点线仅用于视觉引导,不代表任何物理意义
In order to explore the different THz wave output from three-color harmonic pulse cases at various frequency ratios, we select several typical frequency ratios as research objects.

Fig. 2 The time-domain waveform of THz waves generated from three-color harmonic pulse schemes at different frequency ratios and the pulse scheme containing the perfect sawtooth electric field, respectively
图2 不同频率比的三色谐波脉冲方案和包含完美锯齿形电场的脉冲方案产生的太赫兹波时域波形

(a)

(b)
Fig. 3 (a) The pulse envelope composed by three-color harmonic pulses at different frequency ratios, and the pulse envelope containing the perfect sawtooth electric field, respectively, (b) the electric field structure in the time range of -6 fs to 6 fs in the composed pulse envelope, the black dashed line is the grid line
图3 (a)不同频率比的三色谐波脉冲合成的光脉冲包络和包含完美锯齿形电场的光脉冲包络,(b)合成脉冲包络中-6 fs至6 fs时间范围内的电场结构,黑色虚线是网格线
Both the electron drift velocity and the time-varying electron density are important parameters related to the THz yield and determined by the electric field of the composed pulse. The electron drift velocity represents the ability to drift away from their parent ions. The electron drift velocity also contains the information on both the shape and asymmetry of the electric field of the composed pulse

(a)

(b)

(c)

(d)

(e)

(f)
Fig. 4 In the time range of -21 fs to -10 fs, the electric field structure, electron drift velocity, and ionization rate corresponding to the pulse composed by three-color harmonic pulses at different frequency ratios and the pulse containing the perfect sawtooth electric field, (a) 1:2:3, (b) 1:2:4, (c) 1:2:5, (d) 1:2:8, (e) 1:2:10, (f) the perfect sawtooth electric field, respectively
图4 在-21 fs到-10 fs的时间内,不同频率比的三色谐波脉冲叠加合成的光脉冲和包含完美锯齿形电场的光脉冲对应的电场结构、电子漂移速度和电离率,(a) 1:2:3,(b)1:2:4,(c)1:2:5,(d)1:2:8,(e)1:2:10,(f)完美锯齿形电场
For the sawtooth-like electric field of each frequency ratio, the electric field peak and the electron drift velocity peak are out of sync, ionization is interrupted before the velocity reaches the peak, and the velocity during ionization can be obviously divided into two parts according to the positive and negative signs, as shown in
We mark the time corresponding to the electric field peak and the adjacent velocity peak with two black dashed lines perpendicular to the time axis in
Frequency ratio | 1:2:3 | 1:2:4 | 1:2:5 | 1:2:8 | 1:2:10 | sawtooth |
---|---|---|---|---|---|---|
Time interval (fs) | 0.67 | 0.59 | 1.19 | 0.85 | 0.72 | 0 |
The time-varying electron density and ionization rate represent the ability to produce electrons.

Fig. 5 The time-varying electron density of three-color harmonic pulse schemes at different frequency ratios and the pulse scheme containing the perfect sawtooth electric field, respectively
图5 不同频率比的三色谐波脉冲方案和包含完美锯齿形电场的脉冲方案对应的时变电子密度

Fig. 6 The time-varying current density of three-color harmonic pulse schemes at different frequency ratios, and the pulse scheme containing the perfect sawtooth electric field, respectively
图6 不同频率比的三色谐波脉冲方案和包含完美锯齿形电场的脉冲方案对应的时变电流密度
Next, it is considered how to design the wavelength combination scheme of three-color harmonic pulses in the actual operation suitably. At present, for the wavelength combination of the two-color pulse scheme, the basic experimental setup is usually a set of a femtosecond laser amplifier system that can generate 800 nm pulse and a β-barium borate (BBO) crystal that can partially double the frequency of 800 nm pulse to the frequency of 400 nm pulse, then the two-color pulses of 800 nm + 400 nm combination are focused into the air and induced air-plasma to generate THz waves. On the basis of the experimental equipment used in the conventional two-color pulse scheme and the calculation results of 1:2:m, we simulate several specific wavelength combination schemes of femtosecond three-color harmonic pulses at the frequency ratios of 1:2:3 and 1:2:4 which can generate intense THz waves. Subject to the above given pump pulse parameters,

Fig. 7 The time-domain waveform of THz waves generated by different wavelength combination schemes
图7 不同波长组合方案产生的太赫兹波时域波形
Frequency ratio | Wavelength combination (nm) | THz peak-peak electric field amplitude (arb. units) | Approximate THz amplitude ratio for 800+400 nm combination |
---|---|---|---|
1:2 | 800+400 | 1.939 5 | 1 |
1:2:3 | 800+400+267 | 3.055 2 | 1.57 |
1:2:3 | 2 400+1 200+800 | 9.117 9 | 4.70 |
1:2:4 | 800+400+200 | 2.900 7 | 1.50 |
1:2:4 | 1 600+800+400 | 5.742 6 | 2.96 |
1:2:4 | 3 200+1 600+800 | 11.970 9 | 6.17 |
2:3:4 | 800+533+400 | 1.383 3 | 0.71 |
2:3:4 | 1 200+800+600 | 1.909 9 | 0.98 |
By comparison, it can be seen the wavelength combination of 3 200 nm+1 600 nm+800 nm can generate extremely intense THz waves at the corresponding frequency ratio of 1:2:4, the second is from the wavelength combination of 2 400 nm+1 200 nm+800 nm at the frequency ratio of 1:2:3, the third is from the wavelength combination of 1 600 nm+800 nm+400 nm at the frequency ratio of 1:2:4. As for the wavelength combination of 800 nm+400 nm+267 nm at the frequency ratio of 1:2:3, and the wavelength combination of 800 nm+400 nm+200 nm at the frequency ratio of 1:2:4, the output THz waves from them are also stronger than that from the two-color combination of 800 nm+400 nm. When the frequency ratio is 2:3:4, the yields of THz waves generated by both 800 nm+533 nm+400 nm combination and 1 200 nm+800 nm+600 nm combination are weaker than that generated by two-color combination of 800 nm+400 nm, so these two wavelength combinations can be excluded.
In addition, with the number of the harmonic pulses increasing, it is recognized that the control of relevant parameters becomes more complex when the requirements of forming the sawtooth-like electric field are met, for example, proportional control, phase control, etc. Overcoming these difficulties depends on the satisfactory experimental equipment and the excellent experimental skill. We have known the 800 nm + 400 nm + 267 nm combination scheme at the frequency ratio of 1:2:3 is attempted to output intense THz waves via continuously doubling the frequency for 800 nm pulse with the double BBO crystal instead of OPA and modulating the phase using a phase compensation device with attosecond precision
Here, it should be noted that not only the selection of frequency ratio can affect the THz wave generation, but also the use of longer wavelength combination can optimize the THz yield. From the perspective of electric field structure, the longer wavelength can also change the electric field structure in the composed pulse. It is similar to the two-color pulse scheme, with longer wavelengths increasing the THz wave output. We take 800 nm+400 nm+200 nm combination, 1 600 nm+800 nm+400 nm combination and 3 200 nm+1 600 nm+800 nm combination at the frequency ratio of 1:2:4 as examples. With the increase of wavelength,


Fig. 8 (a) The pulse envelope composed of different wavelength combinations at the frequency ratio of 1:2:4, respectively, (b) the electric field structure in the time range of -10 fs to 10 fs in the composed pulse envelope
图8 (a)频率比为1:2:4的不同波长组合方案合成的脉冲包络,(b)合成脉冲包络中-10 fs至10 fs时间范围内的电场结构
In previous studies, the electric field amplitude ratio, relative phase, and duration of three-color pulses are determined according to



Fig. 9 The dependence of the THz electric field peak-peak amplitude on the parameters of the highest harmonic at frequency ratios of 1:2, 1:2:3 and 1:2:4, including (a) the relative phases, (b) the amplitude ratio, (c) the pulse duration, respectively, the circle dots in (b) correspond to the electric field amplitude ratio of the highest harmonic for three-color pulses in Eq. (1)
图9 太赫兹波电场峰峰值随频率比为1:2、1:2:3和1:2:4的光脉冲方案最高次谐波参数变化的关系,包括(a)相对相位,(b)幅值比,(c)脉宽,(b)中的圆点与公式(1)中三色光脉冲最高次谐波的电场幅值比相对应
The THz wave generation from the air-plasma induced via femtosecond three-color harmonic pulses with frequency ratios of 1:2:m (m is a positive integer) has been investigated based on the TPCM and the sawtooth-like electric field formed via the multi-color harmonic pulse superposition. It shows that when the air is saturated ionized and the electron density reaches the same maximum, for the same number of harmonic pulses, the THz conversion efficiency is not always higher when the electric field shape in the composed pulse envelope is closer to the sawtooth waveform and more asymmetric. Furthermore, we present several specific wavelength combination schemes of femtosecond three-color harmonic pulses, which are 800 nm+400 nm+267 nm combination and 2 400 nm+1 200 nm+800 nm combination at the frequency ratio of 1:2:3, 800 nm + 400 nm + 200 nm combination,1 600 nm+800 nm+400 nm combination and 3 200 nm+1 600 nm+800 nm combination at the frequency ratio of 1:2:4, respectively. These specific wavelength combination schemes can significantly enhance the THz wave generation by adding only a set of OPA on the basis of the conventional two-color laser pulse scheme. This optimization of the THz yield is related to both the selection of frequency ratio and the longer wavelength combination. With the progress of laser technology, we believe the three-color scheme and even the multi-color scheme will become the more effective way to generate intense THz waves. Finally, we demonstrate when the relative phase of the highest harmonic is varied, the maximum THz yield of the three-color pulse scheme at the frequency combination of 1:2:4 is slightly smaller than that of the three-color pulse case at the frequency combination of 1:2:3, but its overall THz yield can be maintained the better stability and the smaller fluctuation. When the amplitude ratio of the highest harmonic is appropriately increased, the three-color pulse scheme at the frequency combination of 1:2:4 can generate stronger THz waves than that of the three-color pulse case at the frequency combination of 1:2:3. Our study will be helpful to obtain intense THz sources and provide guidance for practical operations.
References
Tonouchi M. Cutting-edge terahertz technology [J]. Nat. Photonics, 2007, 1(2): 97-105. 10.1038/nphoton.2007.3 [Baidu Scholar]
Davies A G, Burnett A D, Fan W H, et al. Terahertz spectroscopy of explosives and drugs [J]. Mater. Today, 2008, 11(3): 18-26. 10.1016/s1369-7021(08)70016-6 [Baidu Scholar]
Blanchard F, Doi A, Tanaka T, et al. Real-time terahertz near-field microscope [J]. Opt. Express, 2011, 19(9): 8277-8284. 10.1364/oe.19.008277 [Baidu Scholar]
Kampfrath T, Tanaka K, Nelson K A. Resonant and nonresonant control over matter and light by intense terahertz transients [J]. Nat. Photonics, 2013, 7(9): 680-690. 10.1038/nphoton.2013.184 [Baidu Scholar]
Zhang X C, Shkurinov A, Zhang Y. Extreme terahertz science [J]. Nat. Photonics, 2017, 11(1): 16-18. 10.1038/nphoton.2016.249 [Baidu Scholar]
Du H W. Systematical study on the role of laser second harmonic generation in the terahertz air coherent detection [J]. Opt. Laser Technol., 2020, 130: 106312. 10.1016/j.optlastec.2020.106312 [Baidu Scholar]
Zhang D F, Kroh T, Ritzkowsky F, et al. THz-Enhanced DC Ultrafast Electron Diffractometer [J]. Ultrafast Science, 2021, 2021(2): 1-7. 10.34133/2021/9848526 [Baidu Scholar]
Hamster H, Sullivan A, Gordon S, et al. Subpicosecond, electromagnetic pulses from intense laser-plasma interaction [J]. Phys. Rev. Lett., 1993, 71(17): 2725-2728. 10.1103/physrevlett.71.2725 [Baidu Scholar]
Cook D J, Hochstrasser R M. Intense terahertz pulses by four-wave rectification in air [J]. Opt. Lett., 2000, 25(16): 1210-1212. 10.1364/ol.25.001210 [Baidu Scholar]
Xie X, Dai J M, Zhang X C. Coherent Control of THz Wave Generation in Ambient Air [J]. Phys. Rev. Lett., 2006, 96(7): 075005. 10.1103/physrevlett.96.075005 [Baidu Scholar]
Kim K Y, Glownia J H, Taylor A J, et al. Terahertz emission from ultrafast ionizing air in symmetry-broken laser fields [J]. Opt. Express, 2007, 15(8): 4577-4584. 10.1364/oe.15.004577 [Baidu Scholar]
Kim K Y, Taylor A J, Glownia J H, et al. Coherent control of terahertz supercontinuum generation in ultrafast laser–gas interactions [J]. Nat. Photonics, 2008, 2(10): 605-609. 10.1038/nphoton.2008.153 [Baidu Scholar]
Kim K Y. Generation of coherent terahertz radiation in ultrafast laser-gas interactions [J]. Phys. Plasmas, 2009, 16(5): 056706. 10.1063/1.3134422 [Baidu Scholar]
Oh T I, You Y S, Kim K Y. Two-dimensional plasma current and optimized terahertz generation in two-color photoionization [J]. Opt. Express, 2012, 20(18): 19778. 10.1364/oe.20.019778 [Baidu Scholar]
Nguyen A, Kaltenecker K J, Delagnes J C, et al. Wavelength scaling of terahertz pulse energies delivered by two-color air plasmas [J]. Opt. Lett., 2019, 44(6): 1488-1491. 10.1364/ol.44.001488 [Baidu Scholar]
Zhang L L, Wang W M, Wu T, et al. Observation of Terahertz Radiation via the Two-Color Laser Scheme with Uncommon Frequency Ratios [J]. Phys. Rev. Lett., 2017, 119(23): 235001. 10.1103/physrevlett.119.235001 [Baidu Scholar]
Yoo Y J, Jang D, Kim K Y. Highly enhanced terahertz conversion by two-color laser filamentation at low gas pressures [J]. Opt. Express, 2019, 27(16): 22663-22673. 10.1364/oe.27.022663 [Baidu Scholar]
Kim K Y, Glownia J H, Taylor A J, et al. High-power broadband terahertz generation via two-color photoionization in gases [J]. IEEE J. Quantum Electron., 2012, 48(6): 797-805. 10.1109/jqe.2012.2190586 [Baidu Scholar]
Kuk D, Yoo Y J, Rosenthal E W, et al. Generation of scalable terahertz radiation from cylindrically focused two-color laser pulses in air [J]. Appl. Phys. Lett., 2016, 108(12): 121106. 10.1063/1.4944843 [Baidu Scholar]
Zhao J Y, Guo L J, Chu W, et al. Simple method to enhance terahertz radiation from femtosecond laser filament array with a step phase plate [J]. Opt. Lett., 2015, 40(16): 3838-3841. 10.1364/ol.40.003838 [Baidu Scholar]
Houard A, Liu Y, Prade B, et al. Strong Enhancement of Terahertz Radiation from Laser Filaments in Air by a Static Electric Field [J]. Phys. Rev. Lett., 2008, 100(25): 255006. 10.1103/physrevlett.100.255006 [Baidu Scholar]
Wang T J, Ju J J, Liu Y X, et al. Waveform control of enhanced THz radiation from femtosecond laser filament in air [J]. Appl. Phys. Lett., 2017, 110(22): 221102. 10.1063/1.4984599 [Baidu Scholar]
Martínez P G d A, Babushkin I, Bergé L, et al. Boosting Terahertz Generation in Laser-Field Ionized Gases Using a Sawtooth Wave Shape [J]. Phys. Rev. Lett., 2015, 114(18): 183901. 10.1103/physrevlett.114.183901 [Baidu Scholar]
Petersen P B, Tokmakoff A. Source for ultrafast continuum infrared and terahertz radiation [J]. Opt. Lett., 2010, 35(12): 1962-1964. 10.1364/ol.35.001962 [Baidu Scholar]
Lu C H, Zhang C Y, Zhang L Q, et al. Modulation of terahertz-spectrum generation from an air plasma by tunable three-color laser pulses [J]. Phys. Rev. A, 2017, 96(5): 053402. 10.1103/physreva.96.053402 [Baidu Scholar]
Bagley J D, Moss C D, Sorenson S A, et al. Laser-induced plasma generation of terahertz radiation using three incommensurate wavelengths [J]. J. Phys. B At. Mol. Opt. Phys., 2018, 51(14): 144004. 10.1088/1361-6455/aac6ef [Baidu Scholar]
Vaičaitis V, Balachninaite O, Morgner U, et al. Terahertz radiation generation by three-color laser pulses in air filament [J]. J. Appl. Phys., 2019, 125: 173103. 10.1063/1.5078683 [Baidu Scholar]
Ma D N, Dong L Q, Zhang R, et al. Enhancement of terahertz wave emission from air plasma excited by harmonic three-color laser fields [J]. Opt. Commun., 2021, 481: 126533. 10.1016/j.optcom.2020.126533 [Baidu Scholar]
Wang S F, Xiao H C, Peng Y. Mechanism study of terahertz radiation regulation in a multi-color laser field [J]. J. Opt. Soc. Am. B, 2020, 37(11): 3325-3334. 10.1364/josab.399626 [Baidu Scholar]
Alirezaee H, Sharifian M, Darbani S M R, et al. Terahertz radiation emission from three-color laser-induced air plasma [J]. Eur. Phys. J. Plus, 2020, 135(4): 342. 10.1140/epjp/s13360-020-00320-z [Baidu Scholar]
Pei M J, Lu C H, Wang X W, et al. Enhancement and modulation of terahertz radiation by multi-color laser pulses [J]. Chin. Phys. B, 2018, 27(8): 084209. 10.1088/1674-1056/27/8/084209 [Baidu Scholar]
Zhang L, Wang G L, Zhao S F, et al. Controlling of strong tunable THz emission with optimal incommensurate multi-color laser field [J]. Phys. Plasmas, 2017, 24(2): 023116. 10.1063/1.4976549 [Baidu Scholar]
Zhou Z Y, lv Z H, Zhang D W, et al. Analysis of THz generation by multicolor laser pulses with various frequency ratios [J]. Phys. Rev. A, 2020, 101(4): 043422. 10.1103/physreva.101.043422 [Baidu Scholar]
Rodriguez G, Siders C W, Chunlei G, et al. Coherent ultrafast MI-FROG spectroscopy of optical field ionization in molecular H2, N2, and O2 [J]. IEEE J. Sel. Top. Quant. Electron., 2001, 7(4): 579-591. 10.1109/2944.974229 [Baidu Scholar]
Wu P H, Yang C S, Chen P H, et al. Optimization of Terahertz Source via an Ambient Air-Based Multi-Color Photoionization [C]. 2019 44th International Conference on Infrared, Millimeter, and Terahertz Waves (IRMMW-THz), IEEE, 2019, 1-2. 10.1109/irmmw-thz.2019.8873885 [Baidu Scholar]
Liu S J, Fan Z Q, Lu C H, et al. Coherent control of boosted terahertz radiation from air plasma pumped by a femtosecond three-color sawtooth field [J]. Phys. Rev. A, 2020, 102(6): 063522. 10.1103/physreva.102.063522 [Baidu Scholar]