Control of the polarization direction of isolated attosecond pulses using inhomogeneous two-color fields

We theoretically demonstrate the control of the polarization direction of isolated attosecond pulses (IAPs) with inhomogeneous two-color fields synthesized by an 800-nm fundamental pulse and a 2000-nm control pulse having crossed linear polarizations. The results show that by using the temporally and spatially shaped field, the high-order harmonic generation (HHG) process can be efficiently controlled. An ultra-broad supercontinuum ranging from 150th to 400th harmonics which covers the water window region is generated. Such a supercontinuum supports the generation of a 64-as linearly polarized IAP, whose polarization direction is at about 45° with respect to the x axis. Moreover, we analyze the influence of the inhomogeneity parameters and the relative angle of the fundamental and control pulses on the IAP generation. It is shown that the polarization direction of the IAP can rotate in a wide range approximately from 8° to 90° relative to the x axis when the inhomogeneity parameters and the relative angle vary.

In this paper, we theoretically demonstrate control of the polarization direction of IAPs with the cross-LP inhomogeneous two-color field, which is synthesized by an 800-nm fundamental pulse and a 2000-nm control pulse. Based on the quantum results obtained by solving the two-dimensional time-dependent Schrödinger Eq. (2D-TDSE), we show that the HHG can be efficiently controlled by using the inhomogeneous two-color field. Then a 250-th (150-400th) supercontinuum covering the water window region and a 64-as IAP with polarization direction at about 45° with respect to the x axis are produced. Moreover, we have also discussed the influence of the inhomogeneity parameters and the relative angle of the two pulses on the polarization direction of the generated IAPs. The results show that the polarization direction of IAPs can vary in a wide range approximately from 8° to 90° relative to the x axis .

theoretical Model
In our simulations, the harmonic spectra are obtained by solving 2D-TDSE, which is given by (atomic units are used throughout this paper) is the soft-core potential. The soft core parameter α is chosen to be 0.1195 to match the ground ionization potential of neon atom which is 0.7925 a.u (21.6 eV).
is the potential due to laser-electron interaction. The x, y components of the inhomogeneous field are given by Here the parameters ε x and ε y define the strengths of the spatial inhomogeneity of the laser field along x, y directions, respectively 44,48,50 . In our calculations, an 800-nm fundamental pulse and a 2000-nm control pulse with crossed linear polarizations are adopted to synthesize the two-color field. The intensities of these two laser pulses are chosen to be . × W cm 3 0 10 / 14 2 and . × W cm 6 0 10 / 13 2 , respectively. The electric fields of the x, y components of the two-color fields can be written as Here E 0 , E 1 , ω 0 , ω 1 φ 0 and φ 1 are the amplitudes, frequencies and carrier-envelope phases (CEPs) of the fundamental and control pulses, respectively. θ and Δφ are the relative angle and the relative phase between the two pulses, respectively. π = f t t T ( ) sin ( / ) 2 is the laser envelope. T = 10T 0 , T 0 is the optical cycle of the fundamental pulse.
We use the split-operator method to solve Eq. (1) 58 . To avoid spurious reflections from the boundaries, the electron wavefunction is multiplied by a mask function at each time step 59 . The neon atom is in the initial (ground) state before the laser is turned on. The ground state is obtained by imaginary time propagation with the soft-core potential. Once the electron wavefunction ψ x y t ( , , ) is obtained, the time-dependent dipole acceleration along x and y direction is calculated by the Ehrenfest theorem 60

a t x y t H x y t H x y t x x y t a t x y t H x y t H x y t y x y t
x y The HHG spectrum is obtained by Fourier transforming time-dependent dipole acceleration, Here q is the harmonic order. Figure 1(a-c) show the x, y components of the calculated harmonic spectrum and corresponding time-frequency distributions in the homogeneous (ε x = ε y = 0) two-color laser field. The ellipticity of harmonics is also presented (the yellow dot) in Fig. 1(a). Here, the fundamental and control pulses are cross-polarized with the relative angle θ of 0.3π. The relative phase and CEPs of these two laser pulses are set to be 0. The electric fields of the fundamental (the black line) and control (the green line) pulses are shown in the inset of Fig. 1(a). As can be seen from Fig. 1(b), the harmonics above 70th are contributed by two emission peaks, P 1 and P 2 . Among them, the peak P 2 has much higher intensity than P 1 , it dominates the harmonic emission, therefore the harmonics above 70th become continuous. The maximum energy of P 2 is at the 90th harmonic, corresponding to the harmonic cutoff. Besides, the peak P 2 has two branches: the short and long quantum paths. The interference of these quantum paths leads to the modulation structure of the x component of the supercontinuum (70th-90th) (the red line), as shown in Fig. 1(a). Similar results have also been found in the y component of the harmonic spectrum. Due to the interference of the short and long quantum paths of the peak P′ 2 which has much higher intensity than P′ 1 [see Fig. 1(c)], the y component of the supercontinuum above 70th presents modulation structure [see the blue line in Fig. 1(a)]. Moreover, as shown in Fig. 1(a), the harmonic ellipticity of the supercontinuum is 0 (the yellow dot). By superposing the harmonics from 75th to 90th around in the supercontinuum, we can get a LP attosecond double pulse due to the presence of the short and long quantum paths. The polarization direction of the generated LP attosecond double pulse is approximately 33° relative to the x axis [see Fig. 1

Results and Discussion
To generate an IAP, we consider modulating the electron dynamics by using a cross-LP inhomogeneous two-color field. Here, the inhomogeneity parameters of the fundamental and control pulses are taken as ε x = 0.003 and ε y = 0.01, respectively. Other parameters are the same as in Fig. 1(a). Figure 2(a) shows the x, y components of the generated harmonic spectrum in the two-color inhomogeneous field. One can see that the intensities of the x, y components of the harmonic spectrum are comparable. But the harmonic cutoffs are extended to the 400th harmonic and a smooth supercontinuum with a bandwidth of 250th harmonics [from 150th to 400th] covering the water window region is successfully produced. For a deeper insight, we present the corresponding time-frequency distributions in Fig. 2(b-c). One can see that the harmonics above 150th are contributed by only one emission peak P 1 and P′ 1 for the x and y harmonic components, respectively. The maximum energy of the highest emission peaks P 1 and P′ 1 reaches 400th harmonic, corresponding to the harmonic cutoff. Moreover, only the short quantum paths survive for the peaks P 1 and P′ 1 . Therefore a broadband 250th (from 150th to 400th) supercontinuum which spans the water window region, is generated [see Fig. 2(a)]. Furthermore, the harmonic ellipticity of the supercontinuum is still 0 (the yellow dot). By superposing harmonics from 250th to 300th in the supercontinuum, a 64-as LP IAP is obtained, as shown in Fig. 2(d). The polarization direction of the LP IAP is at about 45° relative to x axis [see Fig. 2(d)]. Note that we choose the harmonics in a window of 50 harmonic orders in width to synthesize the attosecond pulse throughout the paper in order to guarantee that an IAP can be generated in the inhomogeneous field when the inhomogeneity parameters and relative angle vary.  Except for the inhomogeneity parameters, other parameters are the same as those in Fig. 1(a). Figure 3. (a,d) The x, y components of the generated harmonic spectrum and the 3D plot of the electric field of the IAP generated by superposing 250th to 300th harmonics in the inhomogeneous two-color field with the inhomogeneity parameter along y direction ε y of 0.005. (b,e,c,f) are same as (a,d), but for the cases of ε y = 0.015 and ε y = 0.02, respectively. Except for ε y , other parameters are the same as those in Fig. 2(a) www.nature.com/scientificreports www.nature.com/scientificreports/ In Fig. 3(a-c), we discuss the influence of the inhomogeneity along y direction on the generated harmonic spectra in the two-color inhomogeneous field. The generated IAPs by superposing 250th to 300th harmonics are also presented in Fig. 3(d-f). Here, the inhomogeneity parameter ε y of the y component of the electric field is chosen to be 0.005, 0.015 and 0.02, respectively. Other parameters are the same as those in Fig. 2(a). One can see that the cutoff extension and smooth supercontinuum generation are robust against the variation of the inhomogeneity. However, the intensities of the x, y components of the generated harmonic spectrum are no longer the same in comparison with the case of ε y = 0.01 [shown in Fig. 2(a)]. For the case of ε y = 0.005 shown in Fig. 3(a), the intensity of the x harmonic component is higher than that of the y harmonic component. By superposing harmonics from 250th to 300th, an IAP with the duration of 77as, whose polarization direction is about 27° relative to the x axis, is generated [as shown in Fig. 3(d)]. While in the case of ε y = 0.015, the y harmonic component has higher intensity than the x harmonic component [as shown in Fig. 3(b)]. An IAP with the duration of 64as and the polarization direction of about 54° relative to the x axis is generated [see Fig. 3(e)]. As ε y increases to 0.02, the magnitude difference between the x, y components of the generated harmonic spectrum is further enlarged. Therefore, the polarization direction of the generated IAP also changes, which is at about 68° relative to the x axis. The duration of the generated IAP is 67as. It is worth noting that the generation of the broadband supercontinuum and ultra-short IAP with different inhomogeneity parameters in the cross-LP inhomogeneous two-color field is similar to that in the commonly used parallel-polarized inhomogeneous two-color field 49,50 . However, the polarization direction of the IAP can change with the variation of the inhomogeneity parameters in the cross-polarized configuration while the variation of inhomogeneity parameters may not change the polarization direction of the IAP in the parallel-polarized configuration.
In Fig. 4, we also investigate the effect of the relative angle θ between the fundamental and control pulses on the generation of the IAP. In our simulations, the relative angle θ is chosen to be 0.1π, 0.2π, and 0.4π, respectively. Other parameters are the same as those in Fig. 2(a). The IAPs are generated by superposing the harmonics from 250th to 300th in the water window region. One can see that LP IAPs with the duration below 80as are generated Figure 4. The 3D plot of the electric field of the IAP generated by superposing 250th to 300th harmonics for the cases of the relative angle θ between the fundamental and control pulses of (a) 0.1π, (b) 0.2π and (c) 0.4π, respectively. Except for θ, other parameters are the same as those in Fig. 2(a). www.nature.com/scientificreports www.nature.com/scientificreports/ for all the relative angle values. Moreover, the polarization directions of the generated IAPs are at about 8°, 23° and 57° with respect to the x axis, when the relative angle θ is chosen to be 0.1π, 0.2π, and 0.4π, respectively.
Finally, in Fig. 5, we present variation of the polarization direction of the generated IAPs in the inhomogeneous two-color field when the inhomogeneity parameters ε x , ε y and the polarization angle θ vary. In Fig. 5(a) the IAPs are generated by superposing harmonics from 200th to 250th harmonics while in Fig. 5(b) the IAPs are obtained by superposing harmonics from 250th to 300th harmonics. It is obvious that the polarization directions of the generated IAPs change with the variations of the inhomogeneity parameters and the polarization angle of the two-color fields, which can rotate in a wide range approximately from 8° to 90° with respect to the x axis.

conclusion
In conclusion, we have theoretically demonstrated the control of the polarization direction of the LP IAP using the two-color laser field having the crossing linear polarization. Based on the numerical solution of the 2D-TDSE, we show that HHG process can be efficiently controlled with the temporally and spatially shaped laser fields. An ultra-broad supercontinuum ranging from 150th to 400th which covers the water window region is obtained. By filtering out harmonics from 250th to 300th, a LP IAP in the water-window region with the duration of 64as is generated, whose polarization direction is at approximately 45° relative to the x axis. Moreover, the influence of the inhomogeneity parameters of the fundamental and control pulses and the relative angle between these two pulses on the generated IAP has also been investigated. It shows that the polarization directions of the generated IAPs can rotate in a wide range approximately from 8° to 90° with respect to the x axis when the inhomogeneity parameters and the relative angle vary.