Co-propagation of femtosecond vortex beam and the generated third harmonic in air*

ZHAN Lindi， ZHONG Zhiping， XU Litong， XI Tingting?， HAO Zuoqiang

(1 School of Physical Sciences, University of Chinese Academy of Sciences，Beijing 101408, China； 2 Shandong Provincial Engineering and Technical Center of Light Manipulations & Shandong Provincial Key Laboratory of Optics and Photonic Device, School of Physics and Electronics, Shandong Normal University，Jinan 250358, China)

Abstract The third harmonic (TH) vortex generated by the filamentation of femtosecond vortex beam in air is simulated. The influences of perturbation distribution and peak power of the initial fundamental wave (FW) on the evolution of the two-color vortex are investigated. For the initial FW with symmetrical perturbation, the intensity annulus of the two-color vortex is not influenced by the initial peak power and can be kept in the filamentation regime. But the spiral phase distribution is influenced by multiple filamentation under the high peak power condition. In this case, the spiral phase distribution of the two-color vortex is destroyed due to the phase accumulation of the multiple filaments. For the initial FW with asymmetrical perturbation, both the spatial annulus and the spiral phase distribution of the two-color vortex are destroyed when multiple filaments are formed.

Keywords filamentation; third harmonic; femtosecond vortex beam; air; perturbation; peak power

Third-harmonic (TH) generated by femtosecond filamentation has attracted great attention due to the involved nonlinear processes[1]and its potential applications in remote sensing and supercontinuum generation[2-5]. During the filamentation of fundamental wave (FW) in air, the peak intensity is clamped at 1013～1014W/cm2. Correspondingly, the peak intensity of TH generated in this process can reach the value that higher than 1012W/cm2[6]. Due to the strong drive by the filamentation of the FW, the TH pulse keeps synchronization with the FW without the walk-off phenomenon[7], and the characteristics of the TH follow the filamentation dynamics of the FW[8]. It is important for the potential applications to investigate the co-propagation of FW and TH when the initial FW has complicated distribution in space and time. Up to date, several works on the co-filamentation of FW and TH have been reported for the FW with different initial conditions, such as different central wavelengths[9-11], different pulse durations[12], modulation by amplitude mask[13], and introducing plasma grating[14]. In recent years, the nonlinear propagation of femtosecond laser pulses with complex spatial and phase structures, such as Airy beam[15], Bessel beam[16]and optical vortex[17], has been intensively studied. Especially, the optical vortex which has a spiral phase structure and an annular intensity distribution with a singularity, can be used in many potential applications such as optical tweezers[18], optical communication[19]and astronomical observations[20]. The study of the TH generated by the filamentation of femtosecond vortex beam is important not only to the understanding of the nonlinear process of the optical vortex but also to the discovery of new potential applications. However, there have been very few reports on the TH generated by the filamentation of the optical vortex. Although the TH vortex is experimentally reported to be generated by the filamentation of the femtosecond vortex beam in air[21], the maintaining conditions of the two-color vortex, which are important to the applications, are still unknown during this nonlinear co-propagation.

In this paper, we perform simulations on the nonlinear co-propagation of the femtosecond vortex beam and the generated TH vortex in air. The influences of the perturbation distribution and peak power of the initial FW on the evolution of two-color vortex are investigated. Considering the real condition of the femtosecond laser beam, we introduce symmetrical and asymmetrical distributions of the perturbation to the FW in the transverse plane. The keeping conditions of the vortex characteristics, including the annular intensity and the spiral phase distributions, are investigated.

1 Numerical simulation

The equations which are used to simulate the TH generated by the filamentation of femtosecond vortex beam in air[6]are written as

(1)

(2)

(3)

The initial distributions of the FW and the TH are written as

(4)

E3ω(r,t,z=0)=0,

(5)

A=1+Ap,

(6)

whereApis the amplitude of the initial perturbation. To study the influence of the initial perturbation on the evolution of the two-color vortex, we assume that the perturbation has symmetrical or asymmetrical distribution, respectively. For the symmetrical case, the amplitude of the perturbation is written as

(7)

where the peak amplitudeε=0.02 is evaluated according to the experimental condition. The beam waist of the perturbation is assumed to ber1=0.15 mm. For the asymmetrical case, the amplitude of the perturbations is written as

(8)

The fluence and phase distributions of the initial FW are plotted in Fig.1.

2 Results and discussion

First, we simulate the co-propagation of the FW vortex beam and the generated TH for the initial FW with symmetrical perturbation. The initial distribution of the perturbation is described by Eq. (7), and the initial fluence distribution of the FW is shown in Fig.1 (a). Figure 2 (a) and 2(b) show the peak intensities of the FW vortex beam and the generated TH for two input peak powers. For both cases, the evolution of the TH intensity follows that of the FW. The TH intensity increases when the filamentation starts, and decreases when the filamentation ends. Correspondingly, the energy and the conversion efficiency of the TH have a similar evolution, as shown in Fig.2 (c) and 2(d). The energy of the TH increases with the formation of filamentation and is given back to the FW with the termination of filamentation. With the increase of the FW peak power from 31.3 GW to 62.7 GW, both the maximal intensities of the FW and TH have a big increase, from 35 TW/cm2to 70 TW/cm2, and from 0.25 TW/cm2to 1.3 TW/cm2, respectively. Correspondingly, the maximal energy of TH increases from 1 μJ to 4.7 μJ, and the maximal conversion efficiency increases from 0.1% to 0.24%, as shown in Fig.2 (c) and 2(d). On the other hand, with the increase of the power, the onset of the two-color filament becomes earlier, and the filament becomes longer.

(a) Fluence distribution of the FW with symmetrical perturbation; (b) Fluence distribution of the FW with asymmetrical perturbation; (c) Phase distribution of the FW.Fig.1 Fluence distribution and phase distribution of the initial FW

The symmetrical perturbation is initially introduced into the FW vortex beam.Fig.2 Peak intensities of the FW and the TH, and the energy and conversion efficiency of the TH as a function of propagation distance for two input peak powers

Figure 3 (a) and 3(b) show the fluence distributions of the FW and the TH at several typical propagation distances when the initial FW has symmetrical perturbation and peak power 31.3 GW. From Fig.3 (a) and 3(b), we can see that during the filamentation (beforez=3 m), the annular distributions of the FW and the TH in the transverse plane are kept. The evolution of the TH follows that of the FW. When the FW begins to self-focus, four peaks are formed symmetrically in the annulus (see Fig.3 (a)z=1.0 m). Then four filaments are formed with the decrease of the radius of the annulus (see Fig.3 (a)z=2.0 m). The positions of the four filaments are symmetrical and rotational. With the further propagation of the laser, the radius of the annulus increases and the filamentation disperses gradually and tends to termination. Because the TH is mainly generated by the intense FW filaments. The annulus of the TH is changed after the formation of the four filaments (see Fig.3 (b)z=2.0 m). Most TH energy is concentrated in the four filaments. Moreover, the rotation of the TH filaments keeps synchronous with that of the FW. These results suggest that both the FW and the generated TH can keep the spatial distribution of the vortex in the filamentation regime.

The symmetrical perturbation is initially introduced into the FW vortex beam which has a peak power of 31.3 GW.Fig.3 Fluence distributions and phase distributions of the FW and the TH at several typical propagation distances

To make sure whether the FW and the TH are vortices during the co-filamentation, we plot the phase distributions of the two fields when the initial FW has symmetrical perturbation and peak power 31.3 GW, as shown in Fig.3 (c) and 3(d). The initial phase distribution of the FW can be seen from the Fig.1 (c). Because the initial topological charge of the FWmωequals 1, the change of phase in a circle is 2π, and the phase only depends on the azimuthal angle. When the FW vortex beam begins to self-focus and the intensity is not so high (z=1.0 m), the spiral distribution of the phase is kept. Because the intensity distribution relies on the radius for the FW, the phase modulation by Kerr effect results in that the phase distribution not only depends on the azimuthal angle but also depends on the radius. During this process, the FW keeps the characteristics of vortex. Correspondingly, for the TH, the spiral phase distribution is also observed atz=1.0 m. The change of the phase in a circle is 6π, three times that of the FW phase. It suggests that the topological charge of the THm3ωequals 3. The conservation of the topological charge is in accordance with that in Ref. [21]. In the filamentation region (z=2.0 m), the phase distributions of the FW and TH have a slight change. Because the FW filaments intensity is much higher than that in other regions, the phase accumulation due to Kerr effect is much larger. The phase accumulation in the filaments region results in a distortion in the filamentation annulus. The spiral phase in the annulus is destroyed slightly. Correspondingly, the phase distribution of the TH is also distorted in the annulus. From Fig.3 (c) and 3(d), we can also see that during the dispersion of the filaments, the most parts of the FW and TH keeps the spiral phase with a slight distortion (z=3.0 m andz=4.0 m). Therefore, during the co-filamentation, the most parts of the two fields keep the phase distribution of vortices. The fluence evolution and the phase distribution suggest that the topological charge of the TH equals 3, and both the FW and the TH keep the vortex characteristics in the filamentation regime.

The evolution of the two-color vortex is also investigated when the initial FW vortex beam has a higher power of 62.7 GW and symmetrical perturbation. In this case, the fluence and phase distributions of the FW and the TH are shown in Fig.4.

The symmetrical perturbation is initially introduced into the FW vortex beam which has a peak power of 62.7 GW. Out of the black circle the phase distribution of TH is distorted in (d) at z=2.0 m.Fig.4 Fluence distributions and phase distributions of the FW and the TH at several typical propagation distances

During the co-propagation, the evolution of the TH also follows that of the FW. From Fig.4 (a) and 4(b), we can see that for the FW and the TH, four filaments are formed atz=0.76 m. Correspondingly, the spiral phase distributions of the two fields are destroyed slightly in the annulus, as shown in Fig.4 (c) and 4(d). Then, with the interaction of the four filaments, new filaments are generated. The filaments are distributed as a necklace. With the formation of more filaments, the spiral phase distributions are destroyed in most part, as shown in Fig.4 (c) and 4(d) atz=1.0 m. After the termination of the filamentation, the annulus distributions of the intensity are distorted greatly (see Fig.4 (a) and 4(b)z=2.0 m). The phase distributions in most parts of the two fields are also distorted. These results suggest that for the FW with symmetrical perturbation and much higher peak power, both the FW and the TH can maintain the annular intensity distribution of the vortex in the filamentation regime. But the phase distribution of the vortex cannot be maintained.

We also investigate the co-propagation of the FW and the generated TH when the initial FW vortex beam has an asymmetrical perturbation. The initial distribution of the perturbation is described by the Eq. (8), and the initial fluence distribution of the FW is shown in the Fig.1 (b). Figure 5 (a) shows that for the FW with peak power of 31.3 GW, the evolutions of the peak intensity for both the FW and the TH are similar to those in the symmetrical perturbation condition. But the filament intensity is higher than that in the symmetrical perturbation condition, as shown in the Fig.2 (a). Correspondingly, the evolutions of the energy and the conversion efficiency of the TH are also similar to those in the symmetrical perturbation condition. The maximal energy of TH is 1.4 μJ, and the maximal conversion efficiency is 0.14%, higher than those in the symmetrical perturbation condition, as shown in Fig.5 (b) and Fig.2 (c).

The asymmetrical perturbation is initially introduced into the FW vortex beam which has a peak power of 31.3 GW.Fig.5 Peak intensities of the FW and the TH, and the energy and conversion efficiency of the TH as a function of propagation distance

In this case, one hot spot is formed due to the self-focus of the asymmetrical perturbation atz=1.5 m, as shown in Fig.6 (a). Then the multiple filaments are generated in the annulus atz=2.0 m. With the evolution of the multiple filaments, the intensity annulus is destroyed, as shown in Fig.6 (a) atz=3.0 m andz=4.0 m. Figure 6 (b) shows that the evolution of the generated TH follows that of the FW. The annulus of TH is also destroyed with the generation of the multiple filaments. Therefore, the asymmetrical perturbation of the initial FW leads to the spatial breakdown of the two-color vortex. This is due to the asymmetrical energy replenishment during the filamentation, which is caused by the modulation instability[22]. The corresponding phase distributions of the FW and the TH are also plotted in Fig.6 (c) and 6(d). Atz=1.5 m, one spot is formed for the FW and the TH, and the intensities of the FW and the TH are lower. The phase accumulation of this spot is not large. Therefore, the phase distortion of the two fields is not obvious. The spiral phase is also observed for the two fields. After that, the filament intensity becomes higher and higher. With the formation of the multiple filaments, the phase accumulation of these filaments results in the obvious distortion in the phase distributions of the two fields. In this case, the spiral phase of the vortex for both fields is not observed. Therefore, the initial asymmetrical perturbation of the FW results in the damage of the two-color vortex, including the annular intensity distribution and the spiral phase.

3 Conclusion

In this paper, we have simulated the co-propagation of the femtosecond vortex beam and the generated TH vortex in air. The influences of the power and the perturbation distribution on the evolution of the two-color vortex have been investigated. It is shown that for the initial FW with symmetrical perturbation, the spatial annulus of the two-color vortex can be maintained in the filamentation regime. For the low peak power, the spiral phase in most parts of the FW and the TH can also be maintained. But with the increase of the initial peak power, the spiral phase of the two-color vortex is destroyed due to the phase accumulation of the multiple filaments. For the initial FW with asymmetrical perturbation, both the spatial annulus and the spiral phase of the two-color vortex are destroyed when multiple filaments are formed. This study is helpful to understand the nonlinear co-propagation of the two-color vortex and important for the potential applications.

The asymmetrical perturbation is initially introduced into the FW vortex beam which has a peak power of 31.3 GW.Fig.6 Fluence distributions and phase distributions of the FW and the TH at several typical propagation distances