Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

Analysis of Combination Frequencies Arising from the Ionization of Gases by Multicolor Femtosecond Pulses

Analysis of Combination Frequencies Arising from the Ionization of Gases by Multicolor... hv photonics Communication Analysis of Combination Frequencies Arising from the Ionization of Gases by Multicolor Femtosecond Pulses 1,2, 1,2 Ivan Laryushin * and Alexander Romanov Institute of Applied Physics (IAP RAS), 603950 Nizhny Novgorod, Russia; mailto:romanoval@ipfran.ru Department of Electrodynamics, Lobachevsky State University of Nizhny Novgorod, 603022 Nizhny Novgorod, Russia * Correspondence: ilaryushin@mail.ru Abstract: Ionization-induced multiwave mixing is attracting much interest nowadays due to the possibility of generating short pulses of secondary radiation over a very wide spectral range, from terahertz to far ultraviolet. This paper presents an analytical method for calculating the amplitudes of arbitrary spectral components of free electron currents arising under the action of multicolor ionizing laser pulses. We show that this method can be used to obtain the dependences of characteristics of a frequency-tunable third harmonic of the intense component of a three-color pulse obtained in an optical parametric generator. The obtained results are in good agreement with quantum mechanical calculations. Keywords: ionization; nonlinear optics; wavemixing; TDSE 1. Introduction At present, the processes associated with the generation of secondary radiation in gases under the action of ionizing laser pulses are widely studied in connection with the possibility of creating compact sources of frequency-tunable super-broadband coherent radiation in wide ranges of central wavelengths (100–10,000 nm) [1–5]. Within this range, Citation: Laryushin, I.; Romanov, A. several different subranges can be distinguished, each of which is relevant for different ap- Analysis of Combination Frequencies plied areas; as a rule, such ultrashort pulses are in demand for ultra-precise measurements Arising from the Ionization of Gases and spectroscopy [6,7], control of various quantum or chemical processes [8,9], biomed- by Multicolor Femtosecond Pulses. ical imaging [10], while the possibility of frequency tuning further expands the range of Photonics 2022, 9, 444. https:// applicability of such sources. This can be done using the frequency tuning of the ionizing doi.org/10.3390/photonics9070444 pulse components in optical parametric generators [11–13]. Recently, the possibility of Received: 6 April 2022 using multi-color (three or more) pulses instead of traditional two-color ionizing pulses Accepted: 20 June 2022 has been increasingly investigated [13–16]. The solution of optimization problems for such Published: 23 June 2022 systems is complicated by the fact that the extremely nonlinear nature of ionization-induced multiwave mixing (IMWM) manifests itself in a strong interdependence of the optimal Publisher’s Note: MDPI stays neutral values of each pulse parameter on the values of all others [14,17–19]. with regard to jurisdictional claims in In this work, we propose a new method for calculating the combination frequencies published maps and institutional affil- iations. in the spectrum of the free electron current generated under the action of a multicolor laser pulse. This method does not rely on the smallness of the intensities of any compo- nents compared to others, or on a large number of mixing waves, and is also suitable for arbitrary elliptical polarizations of quasimonochromatic components. We show that our Copyright: © 2022 by the authors. results are in good agreement with the solutions of three-dimensional time-dependent Licensee MDPI, Basel, Switzerland. Schrödinger equation. This article is an open access article Section 2 of this paper contains a description of the semiclassical approach and the distributed under the terms and derivation of closed expressions for combination frequencies in the free electron current conditions of the Creative Commons spectrum. Section 3 presents the results of semiclassical and quantum-mechanical numeri- Attribution (CC BY) license (https:// cal calculations for two- and three-color ionizing pulses with different polarizations of the creativecommons.org/licenses/by/ quasimonochromatic components and the discussion of the obtained results. 4.0/). Photonics 2022, 9, 444. https://doi.org/10.3390/photonics9070444 https://www.mdpi.com/journal/photonics Photonics 2022, 9, 444 2 of 9 2. Analytical Model In our analytical model, we use the so-called semiclassical approach [2,4,5,14,17,19]. To calculate the current density of free electrons generated by the laser field, we use the system of equations ¶N = (N N)w(jEj), (1) ¶t ¶j e = NE, (2) ¶t m with the initial conditions N(t ! ¥) = 0, j(t ! ¥) = 0. Here, N is free-electron density, N is initial density of neutral particles, w(jEj) is tunneling ionization probability per unit time under the action of electric field E, j is the free-electron current density, e and m are electron charge and mass. The electric field is given by ia w t E = f (t) E e + c.c., (3) å k k=1 where f (t) is a slow envelope, a are some natural numbers, l is the total number of quasi- monochromatic components of the pulse, E are the complex amplitudes of components with frequencies a w k 0 if if kx ky E = E e x + E e y , (4) k kx 0 ky 0 x and y are unit vectors along X and Y axes, E and f are the amplitudes of 0 0 kx,y kx,y projections and phase shifts. In this case, we do not assume the field as a successive series of harmonics, w has the meaning of the least common divisor of frequencies (or the least common multiple of wavelengths) for example, for a pulse of two components at 800 nm and 1200 nm, such a frequency will correspond to a wavelength of 2400 nm (and such a pulse itself, within the framework of the proposed theory, should be understood as a combination of the second and third harmonics of this fundamental frequency). We do not assume that a are small, and the following calculations are also applicable to incommensurate multicolor laser pulses. Since the carrier of laser field is periodic with a period 2p/w , we can also represent the dependence of the ionization probability on time as a superposition of harmonics of the fundamental frequency: isw t w(jE(t)j) = w e , (5) å s where w are slow amplitudes. Moreover, the derivative of the free electron current density isw t can also be expanded in a similar series dj/dt = å dj /dte , with coefficients dj e dN E w E w s+a s k sa k k k + , (6) dt imw w dt s a s + a 0 0 k k k=1 0 0 where N(t) = N [1 exp ( w (t )dt )]. These amplitudes of the spectral components of dj/dt determine the energy spectrum of the secondary radiation generated during the ionization of a gas target. So far, the only unknown quantities in expressions (5) and (6) are the amplitudes w , the rest of this subsection is devoted to finding their values. Traditionally, the polynomial dependence of the ionization probability on the strength of the ionizing field is associated with the multiphoton ionization regime. In this paper, we propose the use of the polynomial approximation also for the probability of tunneling ionization, for example, given by the expression [20] Photonics 2022, 9, 444 3 of 9 E 2E 12jEj a a w(jEj) = 4W exp , (7) jEj 3jEj E where E and W are atomic field and frequency. For practical purposes, it is sufficient a a to carry out such an approximation for a segment limited by the value of E , equal lim to the sum of the amplitudes of all quasi-monochromatic components of the ionizing pulse (since the field strength modulus cannot exceed this sum). The use of polynomial approximation makes it possible, as will be shown below, to express the combination components amplitudes as finite sums of combinations of the ionizing pulse parameters. There are several different ways to expand a given function on a segment into a poly- nomial series, for example, using Legendre polynomials that form an orthogonal system on the segment [1, 1] [21]. Since the argument of w is strictly positive, it is convenient for us to extend it in an even way to the region of negative values. For linearly polarized fields, this corresponds to using the field strength projection E as an argument of this function. 1 d 2 n The n-th Legendre polynomial is defined by a formula P (z) = (z 1) (in our n n n 2 n! dz case z = E/E ), and the expansion coefficients for decomposition w(z) = å a P (z) for n n lim R R 1 1 an arbitrary function can be calculated as a = ( w(z)P (z)dz)/( P (z)dz) (since w n n 1 1 is a function of the modulus of the field strength, we can use only even polynomials and integrate from 0 to 1). This series converges to the original smooth function at each point, and the remainder term for the final expansion satisfies the condition of the least square de- viation from the original function. There are other ways of approximating a given function by a polynomial; we will carry out further reasoning for an arbitrary polynomial, without relying on specific methods of expanding the function. As a suggestion, for simplicity one can use a binomial approximation of the form n +2 n 1 1 w (E) = w(E ) (n n )(E/E ) + (n + 2 n )(E/E ) /2, (8) bin lim 0 1 lim 1 0 lim where n = E w (E )/w(E ) is the effective exponent of w at the point E , and 0 lim lim lim lim n = 2bn /2c. The comparison of the original function (7) with this approximation for 1 0 14 2 E = 0.1E (which for a one-color pulse corresponds to intensity of 3.5 10 W/cm ) lim is shown in the Figure 1. This approximation is convenient primarily because it allows you to reduce the number of necessary calculations to a minimum. However, the error in such a binomial approximation is rather difficult to calculate accurately. On the contrary, when using expansions in Legendre or Chebyshev polynomials, we obtain a guaranteed smallness of the remainder term, from which, due to the linearity of the Fourier transform, follows the smallness of the error in estimating any combination frequency w . Let us represent the square of the field modulus (3) in the form 2 2 i2a w t k 0 jE(t)j / f (t) = I + I e + 0 2a i(a +a )w t i(a a )w t n m 0 n m 0 + I e + I e + c.c. , a +a a a å n m å n m n6=m n6=m (9) I = E E , 0 å k I = E , I = 2E E , I = 2E E . 2a a +a n m a a n k n m n m m Here, vector multiplication is a scalar product without complex conjugation, XY = (XY ). To calculate the harmonics of an even polynomial function, we use the following. Let us ist ist ist have the product of two Fourier series ( A e )( B e ) = C e , then C = A B . å å å å s s s s k k sk By iterating this procedure, we can obtain an expression for the Fourier components of any even power of the field, (n) 2n 2n isw t jE(t)j / f (t) = I e , (10) å s Photonics 2022, 9, 444 4 of 9 where (n) I = I I  ...I I , (11) s s s ss s ...s s å 1 2 n1 1 2 n1 s ,s ,...,s 1 2 n1 summation should be carried out over all non-zero components ofjE(t)j in the expression (9). In certain cases, some of the terms in Equation (9) vanish, for example, for two circular components, either only the sum frequency remains (for opposite directions of polarization rotation) or only the difference frequency (for codirectional rotation), while for orthogonal linear polarizations, on the contrary, both the sum and difference frequencies vanish (see Figure 2). Figure 1. Comparison of the dependence of tunneling ionization probability per unit time on the field strength given by Equation (7) (solid red curve) and the binomial approximation given by Equation (8) with n  4.5 and n = 4 (dashed black curve) on a segment (0, 0.1E ). 0 1 a Finally, for a polynomial 2n w(jE(t)j) = a jE(t)j , (12) å n we obtain (n) 2n w = f (t)a I . (13) s n s Substituting (13) into (6) we obtain a closed analytical formula for the pulse shape of an arbitrary combination frequency of the generated secondary radiation. The obtained analytical expressions, although not always compact, can be used by algorithms for the analytical calculation of the objective function and its derivatives in solving the problems of optimizing the generation of a given combination frequency in the secondary radiation spectrum [22]. The complexity of direct calculations of such optima can be quite large, since each quasi-monochromatic component, even with a given envelope shape, is characterized by a frequency and four more parameters that determine its polarization. In the next section of this work, we consider the simplest illustrative properties of this model, which are observed in the results of numerical calculations. Photonics 2022, 9, 444 5 of 9 Figure 2. Normalized to the maximum value, natural logarithm of the squared spectrum of dj/dt (SCE) from a two-color pulse in hydrogen with varying frequency ratio with l = 800 nm and fixed 14 2 14 2 intensities of components I = 2 10 W/cm , I = 10 W/cm , FWHM duration t = 50 fs 1 2 p and different polarizations of components: (a) collinear linear polarizations; (b) orthogonal linear polarizations; (c) corotating circular polarizations; (d) counterrotating circular polarizations. Noise is suppressed at the 20 exponential level. 3. Numerical Calculations and Discussion The above model may seem complicated to use: iterations in the expressions (6), (11) and (13) look unwieldy. However, in many cases these expressions can be simplified. Figure 2 shows the results of semiclassical calculations of the spectral components of dj/dt for two-color pulses with fixed intensities and durations of the components with a change in the frequency of one of the components, different subplots correspond to different polarizations of the components. The envelopes of the pulse components in the numerical calculations had Gaussian shapes, t ln 4 f (t) = exp , (14) with the pulse FWHM duration t = 50 fs and w corresponding to the wavelength p 1 l = 800 nm. Ionization probability for a hydrogen atom was given by expression (7). Collinear linear polarizations, which are most often used in calculations and in the construction of analytical models, give the most complex picture compared to other op- tions. The fact is that the spectrum of the squared modulus of the field of such two-color pulses contains, in addition to the constant component, the doubled frequencies of both components, as well as the difference and sum frequencies, and all their combinations generated by a strongly nonlinear IMWM seem to be diverse. However, even in this case it is quite easy to analyze Figure 2a and understand from what combinations each of the Photonics 2022, 9, 444 6 of 9 pronounced spectral lines is obtained. At w /w = 1 (this is simply the one-color field) 2 1 the spectrum of dj/dt contains just odd harmonics. In other cases the spectrum contains all the combinations aw + bw where jaj +jbj is a small odd natural number, and the 1 2 bigger this sum, the less efficient is the generation. We note here that setting the field in the form of expression (3) allows us, without additional constructions, to consider the cases of superposition of several channels of multiwave mixing (overlaps of bright straight lines on each of the subplots of Figure 2). For other polarizations of the components, the analysis becomes much simpler. With orthogonal linear polarizations, the spectrum of jE(t)j contains only doubled frequencies of both field components, with codirectional circular polarizations only their difference frequency, and with oppositely rotating only the sum frequency. Accordingly, for bicircular pulses there are much fewer combinations of frequencies in dj/dt spectrum, which can be seen in Figure 2c,d. In several recent experiments were demonstrated the generation of terahertz radiation under the action of two-color and three-color pulses, consisting of the main field at w fre- quency and one or two additional waves from a parametric light generator, the frequencies of which can be detuned from w/2 by some small Dw. In this work we theoretically studied another effect for a three-color pulse with the frequencies of additional fields w/2 Dw and w/2 + Dw, which allows to generate a frequency-tunable third harmonic 3w + 2Dw. For such a pulse we can use the proposed model to determine the interdependence of dj /dt on field amplitudes. The smallest order of wave mixing (number of mixing 3w+2Dw waves) in which the given component appears is 6. If we assign m = E /E and 1 w/2+Dw w m = E /E , and set all three phases to zero, we obtain 2 w w/2Dw 2 2 dj (1 + m ) 2 + 6m 2 2m 3w+2Dw 1 1 1 µ m m + + + /w , (15) 1 2 0 dt w + Dw 2.5w + 3Dw 3.5w + Dw 2.5w + 3Dw 4 4 2 2 2 2 where w µ (1 + m + m + 4[m + m + m m ]). 1 2 1 2 1 2 The obtained analytical results were also compared with the numerical solutions of the three-dimensional time-dependent Schrödinger equation (3D TDSE) " # ¶ h ¯ ih ¯ y(r, t) = r erE(t) + V(r) y(r, t), (16) ¶t 2m where y(t, t) is the electron wavefunction, h ¯ is the reduced Plank constant, V(r) is the parent ion potential. The used potential V(r) for He atom is found based on the density functional theory [23]. The numerical 3D TDSE solution was performed similarly to [23] using expansion of the wavefunction in spherical harmonics with maximum orbital momentum l = 256 and with radial boundary at r = 200r (where r is Bohr radius). The max max B B iwt generated radiation is characterized by Fourier spectrum a = a(t)e dt of the atom dipole acceleration a(t), which determines the time derivative of the electron current density, ¶j/¶t: ¶j e = eN a, a = [EhyjrVjyi]. (17) ¶t m The semiclassical model uses the ionization probability rate for He atom [20], 2/k1 3 3 1 2k 2k ajEj w (jEj) = C exp , (18) He 2/k1 3 2k jEj 3jEj I k which, for simplicity, is written in the atomic units and where I = 0.904 a.u. is ionization potential of He atom, k = 2I , a = 7.0, C = 3.13. Numerically we have investigated the properties of the combination components 3w 2Dw, the results are shown on Figure 3 (the spectral density near the third harmonic Photonics 2022, 9, 444 7 of 9 for pulses with two different values of Dw) and Figure 4 (the dependence of intensities of 3w 2Dw on the intensity of one of the components with fixed intensities of two other in a three-color pulse and Dw = 0.07w). Expression (15) does not correctly describe the entire form of the dependence of dj /dt on m , since to do this, it is necessary to take into 3w+2Dw 1,2 account all the significant higher orders of the IMWM. However, already from it one can see some properties of the generation of a given combination frequency, such as linearity in both amplitudes at low values, and a slight asymmetry: if we put m = m and change the sign of Dw, we see that when the intensities of both additional components are equal, the low-frequency satellite of the third harmonic should have a slightly larger amplitude, which is confirmed by both semiclassical and quantum mechanical calculations. In higher orders there will be different (and more complex) analogous terms in large brackets, but the discussed properties will hold: this component grows linearly and equally with both of additional waves amplitudes, when one of them is small (due to the term m m ) 1 2 and then the growth with m becomes slightly faster than with m . With the use of such an 1 2 expression one can analytically solve various optimization problems, such as searching for the maxima of j /dt under a fixed condition on m and m , e.g., fixed total energy. 3w+2Dw 1 Another interesting example found by us in numerical calculations is that when several channels are superimposed, which occurs when w is a multiple of Dw, the effect of “aligning” the spectral comb occurs (Figure 3b), in our case with a distance between neighboring spectral intensity peaks of 2Dw = 0.2w. Moreover, this happens even if the intensities of the three components are relatively arbitrary, but comparable to each other. 2 2 10 10 TDSE SCE (a) ∆ ω = 0.07ω (b) ∆ ω = 0.1ω I = 0 I = I 0.5ω+∆ ω 0.5ω+∆ ω 0.5ω−∆ ω 0 0 I0.5ω+∆ ω = I0.5ω−∆ ω 10 10 −2 −2 10 10 −4 −4 10 10 −6 −6 10 10 2.0 2.5 3.0 3.5 4.0 2.0 2.5 3.0 3.5 4.0 Harmonicorder Harmonicorder Figure 3. Spectral density S(w) of radiation generated during the interaction of a He atom with three-color laser pulse with sin envelop f (t) and FWHM duration t = 100 fs. The frequencies of the components of the laser pulse w, 0.5w Dw, and 0.5w + Dw correspond to the wavelengths 800 nm, 1860 nm, 1400 nm [panel (a) corresponds to Dw = 0.07w] and 800 nm, 2000 nm, 1333 nm [panel 15 2 (b) corresponds to Dw = 0.1w], respectively. The intensities of the components are I = 10 W/cm , 14 2 14 2 I = 2 10 W/cm , and I = 0 (solid blue lines), 2 10 W/cm (solid orange and 0.5wDw 0.5w+Dw dash red lines). Solid blue and orange lines are 3D TDSE results, dashed red lines are direct numerical solution of semiclassical equations. The vertical dashed lines correspond to the frequencies 3w 2Dw, 3w, 3w + 2Dw. Spectraldensity, S(ω) Photonics 2022, 9, 444 8 of 9 I I 0.5ω+∆ ω 0.5ω−∆ ω 1.00 (a) (b) −1 −1 10 10 0.95 ∝ I ∝ I 0.5ω−∆ ω 0.5ω+∆ ω −2 −2 10 10 n=2.86 0.90 n=3 n=3.14 0.85 −3 −3 10 10 12 13 14 15 12 13 14 15 10 10 10 10 10 10 10 10 2 2 Intensity, I (W/cm ) Intensity, I (W/cm ) 0.5ω−∆ ω 0.5ω+∆ ω Figure 4. The intensity of the generated radiation at frequencies 2.86w (blue lines), 3w (solid black lines), and 3.14w (red lines) as a function of the intensity of laser pulse components 0.5w Dw (a) and 0.5w + Dw (b). The parameters of the laser pulse correspond to Figure 3a. The dashed-dotted gray line in panel (b) is the dependence of the final degree of ionization of He atom on the intensity of the laser pulse component 0.5w + Dw. The vertical dashed lines correspond to the intensity 14 2 2 10 W/cm . 4. Conclusions We have developed an analytical model for calculating the amplitudes of arbitrary spectral components of free electron currents arising under the action of multicolor ionizing pulses. We show that this method can be used to obtain the dependences of characteristics of a frequency-tunable third harmonic of the intense component of a three-color pulse obtained in an optical parametric generator. Our theoretical conclusions based on the analysis of semiclassical equations for free electron currents are supported by the results of quantum mechanical calculations using time-dependent Schrödinger equation in the framework single-active electron approximation for the helium atom. The proposed model can be useful in designing circuits utilizing multicolor ionizing laser radiation for generating short pulses in a wide frequency range from mid-IR to XUV. Author Contributions: Analytical model, I.L.; numerical calculations, A.R. and I.L. All authors have read and agreed to the published version of the manuscript. Funding: This research was funded by the Russian Science Foundation (Grant No. 21-72-00034). Data Availability Statement: The data that support the findings of this study are available from the corresponding author upon reasonable request. Conflicts of Interest: The authors declare no conflict of interest. Abbreviations The following abbreviations are used in this manuscript: IMWM ionization-induced multiwave mixing FWHM full width at half maximum TDSE time-dependent Schrödinger equation SCE semiclassical equations References 1. Brunel, F. Harmonic generation due to plasma effects in a gas undergoing multiphoton ionization in the high-intensity limit. JOSA B 1990, 7, 521–526. [CrossRef] 2. Kostin, V.A.; Vvedenskii, N.V. Generation of few-and subcycle radiation in midinfrared-to-deep-ultraviolet range during plasma production by multicolor femtosecond pulses. Phys. Rev. Lett. 2018, 120, 065002. [CrossRef] [PubMed] Spectraldensity, S(nω) Finaldegreeof ionization, N /N t→∞ m Photonics 2022, 9, 444 9 of 9 3. Tamuliene, ˙ V.; Juškevici ˇ ut ¯ e, ˙ G.; Buožius, D.; Vaicaitis, ˇ V.; Babushkin, I.; Morgner, U. Influence of tunnel ionization to third- harmonic generation of infrared femtosecond laser pulses in air. Sci. Rep. 2020, 10, 17437. [CrossRef] [PubMed] 4. Silaev, A.A.; Romanov, A.A.; Vvedenskii, N.V. Generation of tunable mid-and far-infrared pulses during gas ionization by a chirped two-color laser field. Opt. Lett. 2020, 45, 4527–4530. [CrossRef] [PubMed] 5. Silaev, A.A.; Kostin, V.A.; Laryushin, I.D.; Vvedenskii, N.V. Ionization mechanism of the generation of tunable ultrashort pulses in the mid-infrared range. JETP Lett. 2018, 107, 151–156. [CrossRef] 6. Pazourek, R.; Reduzzi, M.; Carpeggiani, P.A.; Sansone, G.; Gaarde, M.; Schafer, K. Ionization delays in few-cycle-pulse multipho- ton quantum-beat spectroscopy in helium. Phys. Rev. A 2016, 93, 023420. [CrossRef] 7. Blank, V.; Thomson, M.D.; Roskos, H.G. Spatio-spectral characteristics of ultra-broadband THz emission from two-colour photoexcited gas plasmas and their impact for nonlinear spectroscopy. New J. Phys. 2013, 15, 075023. [CrossRef] 8. Kampfrath, T.; Tanaka, K.; Nelson, K.A. Resonant and nonresonant control over matter and light by intense terahertz transients. Nat. Photonics 2013, 7, 680. [CrossRef] 9. Kübel, M.; Siemering, R.; Burger, C.; Kling, N.G.; Li, H.; Alnaser, A.S.; Bergues, B.; Zherebtsov, S.; Azzeer, A.M.; Ben-Itzhak, I.; et al. Steering proton migration in hydrocarbons using intense few-cycle laser fields. Phys. Rev. Lett. 2016, 116, 193001. [CrossRef] 10. Maibohm, C.; Silva, F.; Figueiras, E.; Guerreiro, P.T.; Brito, M.; Romero, R.; Crespo, H.; Nieder, J.B. SyncRGB-FLIM: synchronous fluorescence imaging of red, green and blue dyes enabled by ultra-broadband few-cycle laser excitation and fluorescence lifetime detection. Biomed. Opt. Express 2019, 10, 1891–1904. [CrossRef] 11. Sorenson, S.A.; Moss, C.D.; Kauwe, S.K.; Bagley, J.D.; Johnson, J.A. Enhancing terahertz generation from a two-color plasma using optical parametric amplifier waste light. Appl. Phys. Lett. 2019, 114, 011106. [CrossRef] 12. Vvedenskii, N.V.; Korytin, A.I.; Kostin, V.A.; Murzanev, A.A.; Silaev, A.A.; Stepanov, A.N. Two-color laser-plasma generation of terahertz radiation using a frequency-tunable half harmonic of a femtosecond pulse. Phys. Rev. Lett. 2014, 112, 055004. [CrossRef] [PubMed] 13. Vaicaitis, ˇ V.; Balachninaite, O.; Morgner, U.; Babushkin, I. Terahertz radiation generation by three-color laser pulses in air filament. J. Appl. Phys. 2019, 125, 173103. [CrossRef] 14. Kostin, V.A.; Laryushin, I.D.; Vvedenskii, N.V. Generation of Terahertz radiation by multicolor ionizing pulses. JETP Lett. 2020, 112, 77–83. [CrossRef] 15. Pei, M.J.; Lu, C.H.; Wang, X.W.; Sun, Z.R.; Zhang, S.A. Enhancement and modulation of terahertz radiation by multi-color laser pulses. Chin. Phys. B 2018, 27, 084209. [CrossRef] 16. Liu, S.; Fan, Z.; Lu, C.; Gui, J.; Luo, C.; Wang, S.; Liang, Q.; Zhou, B.; Houard, A.; Mysyrowicz, A.; et al. Coherent control of boosted terahertz radiation from air plasma pumped by a femtosecond three-color sawtooth field. Phys. Rev. A 2020, 102, 063522. [CrossRef] 17. Kostin, V.A.; Laryushin, I.D.; Silaev, A.A.; Vvedenskii, N.V. Ionization-induced multiwave mixing: Terahertz generation with two-color laser pulses of various frequency ratios. Phys. Rev. Lett. 2016, 117, 035003. [CrossRef] 18. Strelkov, V.V. High-order optical processes in intense laser field: Towards nonperturbative nonlinear optics. Phys. Rev. A 2016, 93, 053812. [CrossRef] 19. Laryushin, I. Generation of terahertz radiation by two-color ionizing laser pulses with arbitrary intensity ratios. Photonics 2020, 7, 38. [CrossRef] 20. Tong, X.M.; Lin, C.D. Empirical formula for static field ionization rates of atoms and molecules by lasers in the barrier-suppression regime. J. Phys. B 2005, 38, 2593. [CrossRef] 21. Olver, F.W.; Lozier, D.W.; Boisvert, R.F.; Clark, C.W. NIST Handbook of Mathematical Functions Hardback and CD-ROM; Cambridge University Press: Cambridge, UK, 2010; pp. 435–485. 22. Wright, S.; Nocedal, J. Numerical Optimization; Springer: New York, NY, USA, 1999; pp. 220–242. 23. Romanov, A.A.; Silaev, A.A.; Sarantseva, T.S.; Frolov, M.V.; Vvedenskii, N.V. Study of high-order harmonic generation in xenon based on time-dependent density-functional theory. New J. Phys. 2021, 23, 043014. [CrossRef] http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Photonics Multidisciplinary Digital Publishing Institute

Analysis of Combination Frequencies Arising from the Ionization of Gases by Multicolor Femtosecond Pulses

Laryushin, Ivan; Romanov, Alexander
Photonics , Volume 9 (7) – Jun 23, 2022
Download PDF
9 pages

Loading next page...
 
/lp/multidisciplinary-digital-publishing-institute/analysis-of-combination-frequencies-arising-from-the-ionization-of-9Vi9MikdgY

References (9)

Publisher
Multidisciplinary Digital Publishing Institute
Copyright
© 1996-2022 MDPI (Basel, Switzerland) unless otherwise stated Disclaimer The statements, opinions and data contained in the journals are solely those of the individual authors and contributors and not of the publisher and the editor(s). MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Terms and Conditions Privacy Policy
ISSN
2304-6732
DOI
10.3390/photonics9070444
Publisher site
See Article on Publisher Site

Abstract

hv photonics Communication Analysis of Combination Frequencies Arising from the Ionization of Gases by Multicolor Femtosecond Pulses 1,2, 1,2 Ivan Laryushin * and Alexander Romanov Institute of Applied Physics (IAP RAS), 603950 Nizhny Novgorod, Russia; mailto:romanoval@ipfran.ru Department of Electrodynamics, Lobachevsky State University of Nizhny Novgorod, 603022 Nizhny Novgorod, Russia * Correspondence: ilaryushin@mail.ru Abstract: Ionization-induced multiwave mixing is attracting much interest nowadays due to the possibility of generating short pulses of secondary radiation over a very wide spectral range, from terahertz to far ultraviolet. This paper presents an analytical method for calculating the amplitudes of arbitrary spectral components of free electron currents arising under the action of multicolor ionizing laser pulses. We show that this method can be used to obtain the dependences of characteristics of a frequency-tunable third harmonic of the intense component of a three-color pulse obtained in an optical parametric generator. The obtained results are in good agreement with quantum mechanical calculations. Keywords: ionization; nonlinear optics; wavemixing; TDSE 1. Introduction At present, the processes associated with the generation of secondary radiation in gases under the action of ionizing laser pulses are widely studied in connection with the possibility of creating compact sources of frequency-tunable super-broadband coherent radiation in wide ranges of central wavelengths (100–10,000 nm) [1–5]. Within this range, Citation: Laryushin, I.; Romanov, A. several different subranges can be distinguished, each of which is relevant for different ap- Analysis of Combination Frequencies plied areas; as a rule, such ultrashort pulses are in demand for ultra-precise measurements Arising from the Ionization of Gases and spectroscopy [6,7], control of various quantum or chemical processes [8,9], biomed- by Multicolor Femtosecond Pulses. ical imaging [10], while the possibility of frequency tuning further expands the range of Photonics 2022, 9, 444. https:// applicability of such sources. This can be done using the frequency tuning of the ionizing doi.org/10.3390/photonics9070444 pulse components in optical parametric generators [11–13]. Recently, the possibility of Received: 6 April 2022 using multi-color (three or more) pulses instead of traditional two-color ionizing pulses Accepted: 20 June 2022 has been increasingly investigated [13–16]. The solution of optimization problems for such Published: 23 June 2022 systems is complicated by the fact that the extremely nonlinear nature of ionization-induced multiwave mixing (IMWM) manifests itself in a strong interdependence of the optimal Publisher’s Note: MDPI stays neutral values of each pulse parameter on the values of all others [14,17–19]. with regard to jurisdictional claims in In this work, we propose a new method for calculating the combination frequencies published maps and institutional affil- iations. in the spectrum of the free electron current generated under the action of a multicolor laser pulse. This method does not rely on the smallness of the intensities of any compo- nents compared to others, or on a large number of mixing waves, and is also suitable for arbitrary elliptical polarizations of quasimonochromatic components. We show that our Copyright: © 2022 by the authors. results are in good agreement with the solutions of three-dimensional time-dependent Licensee MDPI, Basel, Switzerland. Schrödinger equation. This article is an open access article Section 2 of this paper contains a description of the semiclassical approach and the distributed under the terms and derivation of closed expressions for combination frequencies in the free electron current conditions of the Creative Commons spectrum. Section 3 presents the results of semiclassical and quantum-mechanical numeri- Attribution (CC BY) license (https:// cal calculations for two- and three-color ionizing pulses with different polarizations of the creativecommons.org/licenses/by/ quasimonochromatic components and the discussion of the obtained results. 4.0/). Photonics 2022, 9, 444. https://doi.org/10.3390/photonics9070444 https://www.mdpi.com/journal/photonics Photonics 2022, 9, 444 2 of 9 2. Analytical Model In our analytical model, we use the so-called semiclassical approach [2,4,5,14,17,19]. To calculate the current density of free electrons generated by the laser field, we use the system of equations ¶N = (N N)w(jEj), (1) ¶t ¶j e = NE, (2) ¶t m with the initial conditions N(t ! ¥) = 0, j(t ! ¥) = 0. Here, N is free-electron density, N is initial density of neutral particles, w(jEj) is tunneling ionization probability per unit time under the action of electric field E, j is the free-electron current density, e and m are electron charge and mass. The electric field is given by ia w t E = f (t) E e + c.c., (3) å k k=1 where f (t) is a slow envelope, a are some natural numbers, l is the total number of quasi- monochromatic components of the pulse, E are the complex amplitudes of components with frequencies a w k 0 if if kx ky E = E e x + E e y , (4) k kx 0 ky 0 x and y are unit vectors along X and Y axes, E and f are the amplitudes of 0 0 kx,y kx,y projections and phase shifts. In this case, we do not assume the field as a successive series of harmonics, w has the meaning of the least common divisor of frequencies (or the least common multiple of wavelengths) for example, for a pulse of two components at 800 nm and 1200 nm, such a frequency will correspond to a wavelength of 2400 nm (and such a pulse itself, within the framework of the proposed theory, should be understood as a combination of the second and third harmonics of this fundamental frequency). We do not assume that a are small, and the following calculations are also applicable to incommensurate multicolor laser pulses. Since the carrier of laser field is periodic with a period 2p/w , we can also represent the dependence of the ionization probability on time as a superposition of harmonics of the fundamental frequency: isw t w(jE(t)j) = w e , (5) å s where w are slow amplitudes. Moreover, the derivative of the free electron current density isw t can also be expanded in a similar series dj/dt = å dj /dte , with coefficients dj e dN E w E w s+a s k sa k k k + , (6) dt imw w dt s a s + a 0 0 k k k=1 0 0 where N(t) = N [1 exp ( w (t )dt )]. These amplitudes of the spectral components of dj/dt determine the energy spectrum of the secondary radiation generated during the ionization of a gas target. So far, the only unknown quantities in expressions (5) and (6) are the amplitudes w , the rest of this subsection is devoted to finding their values. Traditionally, the polynomial dependence of the ionization probability on the strength of the ionizing field is associated with the multiphoton ionization regime. In this paper, we propose the use of the polynomial approximation also for the probability of tunneling ionization, for example, given by the expression [20] Photonics 2022, 9, 444 3 of 9 E 2E 12jEj a a w(jEj) = 4W exp , (7) jEj 3jEj E where E and W are atomic field and frequency. For practical purposes, it is sufficient a a to carry out such an approximation for a segment limited by the value of E , equal lim to the sum of the amplitudes of all quasi-monochromatic components of the ionizing pulse (since the field strength modulus cannot exceed this sum). The use of polynomial approximation makes it possible, as will be shown below, to express the combination components amplitudes as finite sums of combinations of the ionizing pulse parameters. There are several different ways to expand a given function on a segment into a poly- nomial series, for example, using Legendre polynomials that form an orthogonal system on the segment [1, 1] [21]. Since the argument of w is strictly positive, it is convenient for us to extend it in an even way to the region of negative values. For linearly polarized fields, this corresponds to using the field strength projection E as an argument of this function. 1 d 2 n The n-th Legendre polynomial is defined by a formula P (z) = (z 1) (in our n n n 2 n! dz case z = E/E ), and the expansion coefficients for decomposition w(z) = å a P (z) for n n lim R R 1 1 an arbitrary function can be calculated as a = ( w(z)P (z)dz)/( P (z)dz) (since w n n 1 1 is a function of the modulus of the field strength, we can use only even polynomials and integrate from 0 to 1). This series converges to the original smooth function at each point, and the remainder term for the final expansion satisfies the condition of the least square de- viation from the original function. There are other ways of approximating a given function by a polynomial; we will carry out further reasoning for an arbitrary polynomial, without relying on specific methods of expanding the function. As a suggestion, for simplicity one can use a binomial approximation of the form n +2 n 1 1 w (E) = w(E ) (n n )(E/E ) + (n + 2 n )(E/E ) /2, (8) bin lim 0 1 lim 1 0 lim where n = E w (E )/w(E ) is the effective exponent of w at the point E , and 0 lim lim lim lim n = 2bn /2c. The comparison of the original function (7) with this approximation for 1 0 14 2 E = 0.1E (which for a one-color pulse corresponds to intensity of 3.5 10 W/cm ) lim is shown in the Figure 1. This approximation is convenient primarily because it allows you to reduce the number of necessary calculations to a minimum. However, the error in such a binomial approximation is rather difficult to calculate accurately. On the contrary, when using expansions in Legendre or Chebyshev polynomials, we obtain a guaranteed smallness of the remainder term, from which, due to the linearity of the Fourier transform, follows the smallness of the error in estimating any combination frequency w . Let us represent the square of the field modulus (3) in the form 2 2 i2a w t k 0 jE(t)j / f (t) = I + I e + 0 2a i(a +a )w t i(a a )w t n m 0 n m 0 + I e + I e + c.c. , a +a a a å n m å n m n6=m n6=m (9) I = E E , 0 å k I = E , I = 2E E , I = 2E E . 2a a +a n m a a n k n m n m m Here, vector multiplication is a scalar product without complex conjugation, XY = (XY ). To calculate the harmonics of an even polynomial function, we use the following. Let us ist ist ist have the product of two Fourier series ( A e )( B e ) = C e , then C = A B . å å å å s s s s k k sk By iterating this procedure, we can obtain an expression for the Fourier components of any even power of the field, (n) 2n 2n isw t jE(t)j / f (t) = I e , (10) å s Photonics 2022, 9, 444 4 of 9 where (n) I = I I  ...I I , (11) s s s ss s ...s s å 1 2 n1 1 2 n1 s ,s ,...,s 1 2 n1 summation should be carried out over all non-zero components ofjE(t)j in the expression (9). In certain cases, some of the terms in Equation (9) vanish, for example, for two circular components, either only the sum frequency remains (for opposite directions of polarization rotation) or only the difference frequency (for codirectional rotation), while for orthogonal linear polarizations, on the contrary, both the sum and difference frequencies vanish (see Figure 2). Figure 1. Comparison of the dependence of tunneling ionization probability per unit time on the field strength given by Equation (7) (solid red curve) and the binomial approximation given by Equation (8) with n  4.5 and n = 4 (dashed black curve) on a segment (0, 0.1E ). 0 1 a Finally, for a polynomial 2n w(jE(t)j) = a jE(t)j , (12) å n we obtain (n) 2n w = f (t)a I . (13) s n s Substituting (13) into (6) we obtain a closed analytical formula for the pulse shape of an arbitrary combination frequency of the generated secondary radiation. The obtained analytical expressions, although not always compact, can be used by algorithms for the analytical calculation of the objective function and its derivatives in solving the problems of optimizing the generation of a given combination frequency in the secondary radiation spectrum [22]. The complexity of direct calculations of such optima can be quite large, since each quasi-monochromatic component, even with a given envelope shape, is characterized by a frequency and four more parameters that determine its polarization. In the next section of this work, we consider the simplest illustrative properties of this model, which are observed in the results of numerical calculations. Photonics 2022, 9, 444 5 of 9 Figure 2. Normalized to the maximum value, natural logarithm of the squared spectrum of dj/dt (SCE) from a two-color pulse in hydrogen with varying frequency ratio with l = 800 nm and fixed 14 2 14 2 intensities of components I = 2 10 W/cm , I = 10 W/cm , FWHM duration t = 50 fs 1 2 p and different polarizations of components: (a) collinear linear polarizations; (b) orthogonal linear polarizations; (c) corotating circular polarizations; (d) counterrotating circular polarizations. Noise is suppressed at the 20 exponential level. 3. Numerical Calculations and Discussion The above model may seem complicated to use: iterations in the expressions (6), (11) and (13) look unwieldy. However, in many cases these expressions can be simplified. Figure 2 shows the results of semiclassical calculations of the spectral components of dj/dt for two-color pulses with fixed intensities and durations of the components with a change in the frequency of one of the components, different subplots correspond to different polarizations of the components. The envelopes of the pulse components in the numerical calculations had Gaussian shapes, t ln 4 f (t) = exp , (14) with the pulse FWHM duration t = 50 fs and w corresponding to the wavelength p 1 l = 800 nm. Ionization probability for a hydrogen atom was given by expression (7). Collinear linear polarizations, which are most often used in calculations and in the construction of analytical models, give the most complex picture compared to other op- tions. The fact is that the spectrum of the squared modulus of the field of such two-color pulses contains, in addition to the constant component, the doubled frequencies of both components, as well as the difference and sum frequencies, and all their combinations generated by a strongly nonlinear IMWM seem to be diverse. However, even in this case it is quite easy to analyze Figure 2a and understand from what combinations each of the Photonics 2022, 9, 444 6 of 9 pronounced spectral lines is obtained. At w /w = 1 (this is simply the one-color field) 2 1 the spectrum of dj/dt contains just odd harmonics. In other cases the spectrum contains all the combinations aw + bw where jaj +jbj is a small odd natural number, and the 1 2 bigger this sum, the less efficient is the generation. We note here that setting the field in the form of expression (3) allows us, without additional constructions, to consider the cases of superposition of several channels of multiwave mixing (overlaps of bright straight lines on each of the subplots of Figure 2). For other polarizations of the components, the analysis becomes much simpler. With orthogonal linear polarizations, the spectrum of jE(t)j contains only doubled frequencies of both field components, with codirectional circular polarizations only their difference frequency, and with oppositely rotating only the sum frequency. Accordingly, for bicircular pulses there are much fewer combinations of frequencies in dj/dt spectrum, which can be seen in Figure 2c,d. In several recent experiments were demonstrated the generation of terahertz radiation under the action of two-color and three-color pulses, consisting of the main field at w fre- quency and one or two additional waves from a parametric light generator, the frequencies of which can be detuned from w/2 by some small Dw. In this work we theoretically studied another effect for a three-color pulse with the frequencies of additional fields w/2 Dw and w/2 + Dw, which allows to generate a frequency-tunable third harmonic 3w + 2Dw. For such a pulse we can use the proposed model to determine the interdependence of dj /dt on field amplitudes. The smallest order of wave mixing (number of mixing 3w+2Dw waves) in which the given component appears is 6. If we assign m = E /E and 1 w/2+Dw w m = E /E , and set all three phases to zero, we obtain 2 w w/2Dw 2 2 dj (1 + m ) 2 + 6m 2 2m 3w+2Dw 1 1 1 µ m m + + + /w , (15) 1 2 0 dt w + Dw 2.5w + 3Dw 3.5w + Dw 2.5w + 3Dw 4 4 2 2 2 2 where w µ (1 + m + m + 4[m + m + m m ]). 1 2 1 2 1 2 The obtained analytical results were also compared with the numerical solutions of the three-dimensional time-dependent Schrödinger equation (3D TDSE) " # ¶ h ¯ ih ¯ y(r, t) = r erE(t) + V(r) y(r, t), (16) ¶t 2m where y(t, t) is the electron wavefunction, h ¯ is the reduced Plank constant, V(r) is the parent ion potential. The used potential V(r) for He atom is found based on the density functional theory [23]. The numerical 3D TDSE solution was performed similarly to [23] using expansion of the wavefunction in spherical harmonics with maximum orbital momentum l = 256 and with radial boundary at r = 200r (where r is Bohr radius). The max max B B iwt generated radiation is characterized by Fourier spectrum a = a(t)e dt of the atom dipole acceleration a(t), which determines the time derivative of the electron current density, ¶j/¶t: ¶j e = eN a, a = [EhyjrVjyi]. (17) ¶t m The semiclassical model uses the ionization probability rate for He atom [20], 2/k1 3 3 1 2k 2k ajEj w (jEj) = C exp , (18) He 2/k1 3 2k jEj 3jEj I k which, for simplicity, is written in the atomic units and where I = 0.904 a.u. is ionization potential of He atom, k = 2I , a = 7.0, C = 3.13. Numerically we have investigated the properties of the combination components 3w 2Dw, the results are shown on Figure 3 (the spectral density near the third harmonic Photonics 2022, 9, 444 7 of 9 for pulses with two different values of Dw) and Figure 4 (the dependence of intensities of 3w 2Dw on the intensity of one of the components with fixed intensities of two other in a three-color pulse and Dw = 0.07w). Expression (15) does not correctly describe the entire form of the dependence of dj /dt on m , since to do this, it is necessary to take into 3w+2Dw 1,2 account all the significant higher orders of the IMWM. However, already from it one can see some properties of the generation of a given combination frequency, such as linearity in both amplitudes at low values, and a slight asymmetry: if we put m = m and change the sign of Dw, we see that when the intensities of both additional components are equal, the low-frequency satellite of the third harmonic should have a slightly larger amplitude, which is confirmed by both semiclassical and quantum mechanical calculations. In higher orders there will be different (and more complex) analogous terms in large brackets, but the discussed properties will hold: this component grows linearly and equally with both of additional waves amplitudes, when one of them is small (due to the term m m ) 1 2 and then the growth with m becomes slightly faster than with m . With the use of such an 1 2 expression one can analytically solve various optimization problems, such as searching for the maxima of j /dt under a fixed condition on m and m , e.g., fixed total energy. 3w+2Dw 1 Another interesting example found by us in numerical calculations is that when several channels are superimposed, which occurs when w is a multiple of Dw, the effect of “aligning” the spectral comb occurs (Figure 3b), in our case with a distance between neighboring spectral intensity peaks of 2Dw = 0.2w. Moreover, this happens even if the intensities of the three components are relatively arbitrary, but comparable to each other. 2 2 10 10 TDSE SCE (a) ∆ ω = 0.07ω (b) ∆ ω = 0.1ω I = 0 I = I 0.5ω+∆ ω 0.5ω+∆ ω 0.5ω−∆ ω 0 0 I0.5ω+∆ ω = I0.5ω−∆ ω 10 10 −2 −2 10 10 −4 −4 10 10 −6 −6 10 10 2.0 2.5 3.0 3.5 4.0 2.0 2.5 3.0 3.5 4.0 Harmonicorder Harmonicorder Figure 3. Spectral density S(w) of radiation generated during the interaction of a He atom with three-color laser pulse with sin envelop f (t) and FWHM duration t = 100 fs. The frequencies of the components of the laser pulse w, 0.5w Dw, and 0.5w + Dw correspond to the wavelengths 800 nm, 1860 nm, 1400 nm [panel (a) corresponds to Dw = 0.07w] and 800 nm, 2000 nm, 1333 nm [panel 15 2 (b) corresponds to Dw = 0.1w], respectively. The intensities of the components are I = 10 W/cm , 14 2 14 2 I = 2 10 W/cm , and I = 0 (solid blue lines), 2 10 W/cm (solid orange and 0.5wDw 0.5w+Dw dash red lines). Solid blue and orange lines are 3D TDSE results, dashed red lines are direct numerical solution of semiclassical equations. The vertical dashed lines correspond to the frequencies 3w 2Dw, 3w, 3w + 2Dw. Spectraldensity, S(ω) Photonics 2022, 9, 444 8 of 9 I I 0.5ω+∆ ω 0.5ω−∆ ω 1.00 (a) (b) −1 −1 10 10 0.95 ∝ I ∝ I 0.5ω−∆ ω 0.5ω+∆ ω −2 −2 10 10 n=2.86 0.90 n=3 n=3.14 0.85 −3 −3 10 10 12 13 14 15 12 13 14 15 10 10 10 10 10 10 10 10 2 2 Intensity, I (W/cm ) Intensity, I (W/cm ) 0.5ω−∆ ω 0.5ω+∆ ω Figure 4. The intensity of the generated radiation at frequencies 2.86w (blue lines), 3w (solid black lines), and 3.14w (red lines) as a function of the intensity of laser pulse components 0.5w Dw (a) and 0.5w + Dw (b). The parameters of the laser pulse correspond to Figure 3a. The dashed-dotted gray line in panel (b) is the dependence of the final degree of ionization of He atom on the intensity of the laser pulse component 0.5w + Dw. The vertical dashed lines correspond to the intensity 14 2 2 10 W/cm . 4. Conclusions We have developed an analytical model for calculating the amplitudes of arbitrary spectral components of free electron currents arising under the action of multicolor ionizing pulses. We show that this method can be used to obtain the dependences of characteristics of a frequency-tunable third harmonic of the intense component of a three-color pulse obtained in an optical parametric generator. Our theoretical conclusions based on the analysis of semiclassical equations for free electron currents are supported by the results of quantum mechanical calculations using time-dependent Schrödinger equation in the framework single-active electron approximation for the helium atom. The proposed model can be useful in designing circuits utilizing multicolor ionizing laser radiation for generating short pulses in a wide frequency range from mid-IR to XUV. Author Contributions: Analytical model, I.L.; numerical calculations, A.R. and I.L. All authors have read and agreed to the published version of the manuscript. Funding: This research was funded by the Russian Science Foundation (Grant No. 21-72-00034). Data Availability Statement: The data that support the findings of this study are available from the corresponding author upon reasonable request. Conflicts of Interest: The authors declare no conflict of interest. Abbreviations The following abbreviations are used in this manuscript: IMWM ionization-induced multiwave mixing FWHM full width at half maximum TDSE time-dependent Schrödinger equation SCE semiclassical equations References 1. Brunel, F. Harmonic generation due to plasma effects in a gas undergoing multiphoton ionization in the high-intensity limit. JOSA B 1990, 7, 521–526. [CrossRef] 2. Kostin, V.A.; Vvedenskii, N.V. Generation of few-and subcycle radiation in midinfrared-to-deep-ultraviolet range during plasma production by multicolor femtosecond pulses. Phys. Rev. Lett. 2018, 120, 065002. [CrossRef] [PubMed] Spectraldensity, S(nω) Finaldegreeof ionization, N /N t→∞ m Photonics 2022, 9, 444 9 of 9 3. Tamuliene, ˙ V.; Juškevici ˇ ut ¯ e, ˙ G.; Buožius, D.; Vaicaitis, ˇ V.; Babushkin, I.; Morgner, U. Influence of tunnel ionization to third- harmonic generation of infrared femtosecond laser pulses in air. Sci. Rep. 2020, 10, 17437. [CrossRef] [PubMed] 4. Silaev, A.A.; Romanov, A.A.; Vvedenskii, N.V. Generation of tunable mid-and far-infrared pulses during gas ionization by a chirped two-color laser field. Opt. Lett. 2020, 45, 4527–4530. [CrossRef] [PubMed] 5. Silaev, A.A.; Kostin, V.A.; Laryushin, I.D.; Vvedenskii, N.V. Ionization mechanism of the generation of tunable ultrashort pulses in the mid-infrared range. JETP Lett. 2018, 107, 151–156. [CrossRef] 6. Pazourek, R.; Reduzzi, M.; Carpeggiani, P.A.; Sansone, G.; Gaarde, M.; Schafer, K. Ionization delays in few-cycle-pulse multipho- ton quantum-beat spectroscopy in helium. Phys. Rev. A 2016, 93, 023420. [CrossRef] 7. Blank, V.; Thomson, M.D.; Roskos, H.G. Spatio-spectral characteristics of ultra-broadband THz emission from two-colour photoexcited gas plasmas and their impact for nonlinear spectroscopy. New J. Phys. 2013, 15, 075023. [CrossRef] 8. Kampfrath, T.; Tanaka, K.; Nelson, K.A. Resonant and nonresonant control over matter and light by intense terahertz transients. Nat. Photonics 2013, 7, 680. [CrossRef] 9. Kübel, M.; Siemering, R.; Burger, C.; Kling, N.G.; Li, H.; Alnaser, A.S.; Bergues, B.; Zherebtsov, S.; Azzeer, A.M.; Ben-Itzhak, I.; et al. Steering proton migration in hydrocarbons using intense few-cycle laser fields. Phys. Rev. Lett. 2016, 116, 193001. [CrossRef] 10. Maibohm, C.; Silva, F.; Figueiras, E.; Guerreiro, P.T.; Brito, M.; Romero, R.; Crespo, H.; Nieder, J.B. SyncRGB-FLIM: synchronous fluorescence imaging of red, green and blue dyes enabled by ultra-broadband few-cycle laser excitation and fluorescence lifetime detection. Biomed. Opt. Express 2019, 10, 1891–1904. [CrossRef] 11. Sorenson, S.A.; Moss, C.D.; Kauwe, S.K.; Bagley, J.D.; Johnson, J.A. Enhancing terahertz generation from a two-color plasma using optical parametric amplifier waste light. Appl. Phys. Lett. 2019, 114, 011106. [CrossRef] 12. Vvedenskii, N.V.; Korytin, A.I.; Kostin, V.A.; Murzanev, A.A.; Silaev, A.A.; Stepanov, A.N. Two-color laser-plasma generation of terahertz radiation using a frequency-tunable half harmonic of a femtosecond pulse. Phys. Rev. Lett. 2014, 112, 055004. [CrossRef] [PubMed] 13. Vaicaitis, ˇ V.; Balachninaite, O.; Morgner, U.; Babushkin, I. Terahertz radiation generation by three-color laser pulses in air filament. J. Appl. Phys. 2019, 125, 173103. [CrossRef] 14. Kostin, V.A.; Laryushin, I.D.; Vvedenskii, N.V. Generation of Terahertz radiation by multicolor ionizing pulses. JETP Lett. 2020, 112, 77–83. [CrossRef] 15. Pei, M.J.; Lu, C.H.; Wang, X.W.; Sun, Z.R.; Zhang, S.A. Enhancement and modulation of terahertz radiation by multi-color laser pulses. Chin. Phys. B 2018, 27, 084209. [CrossRef] 16. Liu, S.; Fan, Z.; Lu, C.; Gui, J.; Luo, C.; Wang, S.; Liang, Q.; Zhou, B.; Houard, A.; Mysyrowicz, A.; et al. Coherent control of boosted terahertz radiation from air plasma pumped by a femtosecond three-color sawtooth field. Phys. Rev. A 2020, 102, 063522. [CrossRef] 17. Kostin, V.A.; Laryushin, I.D.; Silaev, A.A.; Vvedenskii, N.V. Ionization-induced multiwave mixing: Terahertz generation with two-color laser pulses of various frequency ratios. Phys. Rev. Lett. 2016, 117, 035003. [CrossRef] 18. Strelkov, V.V. High-order optical processes in intense laser field: Towards nonperturbative nonlinear optics. Phys. Rev. A 2016, 93, 053812. [CrossRef] 19. Laryushin, I. Generation of terahertz radiation by two-color ionizing laser pulses with arbitrary intensity ratios. Photonics 2020, 7, 38. [CrossRef] 20. Tong, X.M.; Lin, C.D. Empirical formula for static field ionization rates of atoms and molecules by lasers in the barrier-suppression regime. J. Phys. B 2005, 38, 2593. [CrossRef] 21. Olver, F.W.; Lozier, D.W.; Boisvert, R.F.; Clark, C.W. NIST Handbook of Mathematical Functions Hardback and CD-ROM; Cambridge University Press: Cambridge, UK, 2010; pp. 435–485. 22. Wright, S.; Nocedal, J. Numerical Optimization; Springer: New York, NY, USA, 1999; pp. 220–242. 23. Romanov, A.A.; Silaev, A.A.; Sarantseva, T.S.; Frolov, M.V.; Vvedenskii, N.V. Study of high-order harmonic generation in xenon based on time-dependent density-functional theory. New J. Phys. 2021, 23, 043014. [CrossRef]

Journal

Photonics – Multidisciplinary Digital Publishing Institute

Published: Jun 23, 2022

Keywords: ionization; nonlinear optics; wavemixing; TDSE

There are no references for this article.