Hostname: page-component-78c5997874-4rdpn Total loading time: 0 Render date: 2024-11-07T23:29:18.377Z Has data issue: true hasContentIssue false
  • English
  • Français

Spatial distribution of self-seeded air lasers induced by the femtosecond laser filament plasma

Published online by Cambridge University Press:  18 September 2024

Tao Zeng Open the ORCID record for Tao Zeng [Opens in a new window] ,
Nan Li  and
Yuliang Yi
Tao Zeng*
Affiliation:
School of Physical Science and Technology, Chongqing key Laboratory of Micro & Nano Structure Optoelectronics, Southwest University, Chongqing 400715, PR China
Nan Li
Affiliation:
School of Physical Science and Technology, Chongqing key Laboratory of Micro & Nano Structure Optoelectronics, Southwest University, Chongqing 400715, PR China
Yuliang Yi
Affiliation:
School of Physical Science and Technology, Chongqing key Laboratory of Micro & Nano Structure Optoelectronics, Southwest University, Chongqing 400715, PR China
*
†Email address for correspondence: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

The femtosecond laser filament-induced air laser plays a significant role for the remote sensing of air pollutants. The spatial distributions of air laser intensity were investigated experimentally in previous studies. However, the mechanism of the air laser propagation properties inside the filament plasma has not been quite clear yet. Moreover, few studies have been dedicated to the reproduction of the air laser profile from nitrogen molecules propagating in the filament plasma based on the numerical simulation method. In this study, the lasing action of the air laser from the transition of the first negative (0,0) band of nitrogen ions at 391 nm was simulated during the femtosecond laser filamentation. The beam profile of the air laser changes from a Gaussian or super-Gaussian shape to an outer ring structure by increasing the filament length or nitrogen ion density, which is in accord with the previous experimental result. A multiple-diffraction effect has been proposed to clarify the mechanism of the outer rings beam pattern formation, which is induced by the dynamical interaction between the lasing effect and diffraction effect of the air laser propagating inside the filament plasma. In addition, the amplified air laser power as a function of both the filament length and nitrogen ion density was investigated. Our study would pave the way to improve the energy conversion efficiency and directivity of remote air lasers, which would be significant for remote sensing applications.

Keywords

femtosecond laser filamentair lasernitrogen iron densityspatial distribution
Type
Research Article
Information
Journal of Plasma Physics , Volume 90 , Issue 4 , August 2024 , 905900404
Check for updates
Copyright
Copyright © The Author(s), 2024. Published by Cambridge University Press

1. Introduction

Remote sensing of atmosphere pollutants relies on optical techniques for detecting molecules and atoms based on lasers in order to identify and quantify different species. The light emission collection in the backward direction of the pumped laser is restricted by the non-directional characteristic of the spontaneous emission from molecules of interest. Amplified spontaneous emission based on femtosecond laser filamentation has provided the opportunity to solve this problem and appears as an attractive technique to realize the remote sensing of air pollutants (Luo, Liu & Chin Reference Luo, Liu and Chin2003; Dogariu et al. Reference Dogariu, Michael, Scully and Miles2011; Yuan et al. Reference Yuan, Wang, Teranishi, Sridharan, Lin, Zeng and Chin2013, Reference Yuan, Wang, Lu, Chin and Zeng2014; Zhang et al. Reference Zhang, Lu, Zhang, Fan, Zhou, Liang, Yuan, Zhuang, Dorfman and Liu2021; Xue et al. Reference Xue, Zhang, Guo, Zhang, Qi, Sun, Gong, Lin and Liu2022). The laser intensity inside a light filament core could remain at more than 1013 W cm−2 over a very long distance (Liu Reference Liu2014). Various species in the atmosphere will be ionized, resulting in the emission of the fingerprint fluorescence. Propagating through the filament, the fluorescence signal could be enhanced, which is referred to as an air laser (Kartashov et al. Reference Kartashov, Alisauskas, Andriukaitis, Pugžlys, Shneider, Zheltikov, Chin and Baltuska2012; Ni et al. Reference Ni, Chu, Jing, Zhang, Zeng, Yao, Li, Xie, Zhang, Xu, Chin, Cheng and Xu2013; Zhou et al. Reference Zhou, Zhang, Lu, Liang and Liu2020).

The generation of well-collimated air lasers at 845 nm in the backward and forward directions was investigated with an ultraviolet focused pump beam. By the interaction of the pump beam with molecular oxygen, the doughnut-shaped mode of the air laser was observed with a divergence angle of 40 mrad (Dogariu et al. Reference Dogariu, Michael, Scully and Miles2011). A super-Gaussian shape of a backward ultraviolet air laser at 337 nm and 357 nm with a divergence angle of 1.6 mrad was produced. The air lasing action results from nitrogen molecules pumped with a high-power midinfrared laser pulse (Kartashov et al. Reference Kartashov, Alisauskas, Andriukaitis, Pugžlys, Shneider, Zheltikov, Chin and Baltuska2012). A forward harmonic-seeded air laser from nitrogen molecules driven by intense midinfrared femtosecond laser pulses was reported. The lasing effect refers to the amplification of the harmonic radiation with the presence of population inversion induced by ${B^2}\sum _u^ + \to {X^2}\sum _g^ +$ transition between the ground and excited cation states (Yao et al. Reference Yao, Zeng, Xu, Li, Chu, Ni, Zhang, Chin, Cheng and Xu2011). In addition, the influence of gas pressure on the air laser beam intensity and profile was investigated. It indicates that the forward self-seeded air laser from nitrogen molecules changes from a Gaussian shape to an outer ring pattern with the increase of gas pressure (Liu et al. Reference Liu, Brelet, Point, Houard and Mysyrowicz2013). Moreover, our previous study revealed that a multiple ring pattern of the air laser was generated induced by the nitrogen ions inside a femtosecond laser filament (Zeng et al. Reference Zeng, Zhao, Liu and Chin2014). However, the mechanism of the propagation properties and spatial distribution of the air laser during the femtosecond laser filamentation has not been quite clear up to now.

In this work, the spatial distribution of the air laser propagating through the femtosecond laser filament plasma has been investigated based on numerical simulation. The result indicates that the air laser beam pattern after the filament changes from a Gaussian or super-Gaussian shape to a multiple outer rings structure by increasing the filament length or nitrogen ion density. Moreover, the number of the outer rings increases versus the filament length or nitrogen ion density. The mechanism of the air laser intensity spatial distribution has been revealed. In addition, the amplified air laser beam power as a function of the filament length and nitrogen ion density has been investigated.

2. Theoretical method

As the nitrogen molecules were ionized to $\textrm{N}_2^ +$ due to multiphoton or tunnel ionization effects during the femtosecond laser filamentation process, the ground state ${X^2}\sum _g^ +$ of the ion was pumped to the excited state ${B^2}\sum _u^ +$. The fluorescence of $\textrm{N}_2^ +$ with a wavelength of 391 nm was stimulated during the vibrational transition of the $\textrm{N}_2^ +$ first negative (0,0) band $({B^2}\sum _u^ + \to {X^2}\sum _g^ + )$ (Doering & Yang Reference Doering and Yang1996). The photon density ${N_{e\textrm{,}\;\textrm{391}\;\textrm{nm}}}$ of the stimulated fluorescence signal could be derived from the following (Kartashov & Shneider Reference Kartashov and Shneider2017):

(2.1)\begin{equation}{N_{e\textrm{,}\;\textrm{391}\;\textrm{nm}}} = {W_{21}}\tau {N_{e\textrm{,}\;\textrm{N}_2^ + }},\end{equation}

where ${N_{e\textrm{,}\;\textrm{N}_2^ + }}$ represents the nitrogen ion density; ${W_{21}}$ and $\tau$ denote the probability of the fluorescence emission and the population inversion lifetime, respectively. The fluorescence could be produced only in the region where the nitrogen ion density is larger than a threshold value of $\textrm{1}{\textrm{0}^{\textrm{15}}}\;\textrm{c}{\textrm{m}^{ - 3}}$.

In this work, the spatial distribution of the nitrogen ion density in the filament plasma obeys the following equation:

(2.2)\begin{equation}{N_{e\textrm{,}\;\textrm{N}_2^ + }}(r,l) = {N_{e,\max }}\exp \left[ { - \left( {\frac{{{r^2}}}{{r_0^2}} + \frac{{{l^2}}}{{l_0^2}}} \right)} \right].\end{equation}

Here the nitrogen ion density distribution satisfies the cylindrical symmetry. Here r and l refer to the variables of the axis along the radial direction and the longitudinal direction, respectively. Here r 0 is the radius of the filament plasma and is set to be 50 μm in our simulation, which is a common value in the previous studies (Chin et al. Reference Chin, Hosseini, Liu, Luo, Theberge, Akozbek, Becker, Kandidov, Kosareva and Schroeder2005; Berge et al. Reference Berge, Skupin, Nuter, Kasparian and Wolf2007; Couairon & Mysyrowicz Reference Couairon and Mysyrowicz2007; Kasparian & Wolf Reference Kasparian and Wolf2008; Kandidov, Shlenov & Kosareva Reference Kandidov, Shlenov and Kosareva2009). Here l 0 is the half-length of the filament plasma. Here ${N_{e,\max }}$ denotes the peak nitrogen ion density. Both l 0 and ${N_{e,\max }}$ depend on the initial femtosecond laser beam parameters and focusing conditions.

A stationary, homogeneous and isotropic plasma model was selected without considering the temporal variation. The lifetime of the filament plasma is at the nanosecond time scale (Mehr & Biondi Reference Mehr and Biondi1969; Johnson & Fowler Reference Johnson and Fowler1970; Brian & Mitchell Reference Brian and Mitchell1990; Bahati et al. Reference Bahati, Jureta, Belic, Cherkani-Hassani, Abdellahi and Defrance2001; Chauveau et al. Reference Chauveau, Perrin, Riviere and Soufiani2002; Jain, Gupta & Suk Reference Jain, Gupta and Suk2019), which is much slower than the air laser pulse duration at 391 nm in the picosecond range (Xu et al. Reference Xu, Azarm, Bernhardt, Kamali and Chin2009; Yao et al. Reference Yao, Chu, Liu, Chen, Xu and Cheng2018; Wang et al. Reference Wang, Ding, Wilkins, Athanasakis-Kaklamanakis, Zhang, Liu and Hu2023). Therefore, this simplification of the plasma model is reasonable.

The propagation of the fluorescence from nitrogen molecules at 391 nm in the filament plasma could be based on the following wave equation (Chiron et al. Reference Chiron, Lamouroux, Lange, Ripoche, Franco, Prade, Bonnaud, Riazuelo and Mysyrowicz1999):

(2.3)\begin{equation}2\textrm{i}{k_0}\frac{{\partial A}}{{\partial z}} + \left( {\frac{{{\partial^2}}}{{\partial {x^2}}} + \frac{{{\partial^2}}}{{\partial {y^2}}}} \right)A + 2\frac{{k_0^2}}{{{n_0}}}{n_{\textrm{nl}}}A - \textrm{i}{k_0}\alpha A = 0,\end{equation}

where A denotes the air laser amplitude envelope, and ${k_0}$ is the air laser wavenumber. Here z represents the air laser propagation distance; n 0 is the air refractive index; n nl refers to the nonlinear refractive index corresponding to the defocusing effect of the filament plasma. The simulation of the air laser propagating in the filament plasma was based on the finite difference method. FORTRAN software was used for conducting the simulation procedure.

It is worth mentioning that since we focused exclusively on the spatial distribution of the air laser, the temporal aspect of the emitted light field was not considered during the simulation. The validity of this simplified model for the nonlinear propagation of a laser pulse has been demonstrated in previous studies (Dubietis et al. Reference Dubietis, Tamosauskas, Fibich and Ilan2004; Liu & Chin Reference Liu and Chin2007). Thus, temporal effects including self-steeping, group-velocity dispersion effects were neglected in our simulation model. Actually, this kind of simplification would significantly influence the temporal dynamics of the air laser. However, to a lesser extent it would influence the spatial evolution of the air laser in principle, which is clarified as follows. Typically, the air laser radius inside the filament is approximately 50 μm. The corresponding longitudinal scale of spatial transformation $\textrm{2}\pi {a^2}/\lambda$ is approximately 2 cm. The second-order group velocity dispersion coefficient $k^{\prime\prime}$ for air is 0.2 fs2 cm−1 (Akozbek et al. Reference Akozbek, Scalora, Bowden and Chin2001). For an air laser pulse with pulse duration at the picosecond time scale (Xu et al. Reference Xu, Azarm, Bernhardt, Kamali and Chin2009; Yao et al. Reference Yao, Chu, Liu, Chen, Xu and Cheng2018; Wang et al. Reference Wang, Ding, Wilkins, Athanasakis-Kaklamanakis, Zhang, Liu and Hu2023), the characteristic dispersion length scale ${\tau ^\textrm{2}}/k^{\prime\prime}$ of the temporal changes is much longer than the spatial transformation scale. Therefore, the spatial effects develop much faster than the temporal effects and dominate the air laser intensity distribution (Akozbek et al. Reference Akozbek, Scalora, Bowden and Chin2001).

The amplitude of the spontaneous fluorescence emission at 391.4 nm at the l position of the filament plasma was denoted by A 0(l). The relative value of A 0(l) was achieved by the square root of the intensity of the fluorescence emission. The fluorescence intensity was proportional to the photon density, which could be derived based on (2.1). The fluorescence signal with amplitude of A 0(l) at z = l position in the filament plasma would propagate in the filament with a gain effect due to the population inversion of the nitrogen ions, which leads to the air laser with amplitude A 1(l,z) at a propagation distance z position by conducting the numerical simulation based on (2.3). Then the total air laser amplitude A 2(z) at the z position inside the filament plasma could be derived by the integral of A 1(l,z) along the particular filament length from l = 0 to l = z, as follows:

(2.4)\begin{equation}{A_2}(z) = {A_0}(z) + \int_{l = 0}^{l = z} {{A_1}(l,z)} .\end{equation}

Total air laser amplitude A 2(z) at the z position after the filament plasma could be derived as follows:

(2.5)\begin{equation}{A_2}(z) = \int_{l = 0}^{l = {L_0}} {{A_1}(l,z)} .\end{equation}

Here, L 0 corresponds to the total filament length.

The nitrogen ions induced by the filamentation could be served as a gain medium for the air laser. Here $\alpha$ corresponds to the gain coefficient of the air laser propagating through the filament plasma, which could be derived based on the following relationship (Sprangle et al. Reference Sprangle, Penano, Hafizi, Gordon and Scully2011):

(2.6)\begin{equation}\alpha = \sigma {N_{e\textrm{,}\;\textrm{N}_2^ + }},\end{equation}

where $\sigma$ is the stimulated emission cross-section for the first negative (0,0) band vibrational transition at 391 nm with a constant value of $14.8 \times {10^{ - 18}}\;\textrm{c}{\textrm{m}^2}$ according to the previous work (Doering & Yang Reference Doering and Yang1996). Finally, the total air laser power at any propagation distance inside or after the filament plasma could be accumulated by the amplified fluorescence emission at any formal filament position with the nitrogen ion density larger than the threshold value.

Population inversion between the vibrational levels of the excited ${B^2}\sum _u^ +$ and the ground ${X^2}\sum _g^ +$ states of the nitrogen molecular ion $\textrm{N}_2^ +$ was achieved, enabling lasing at the 391 nm emission wavelength. The fraction of the inverted $\textrm{N}_2^ + {B^2}\sum _u^ +$ ions was calculated in the previous study, which corresponds to a value of approximately 9 % (Doering & Yang Reference Doering and Yang1996). The fraction of the inverted molecules would influence the $\textrm{N}_2^ + {B^2}\sum _u^ +$ state cross-section. Therefore, both the photon density ${N_{e\textrm{,}\;\textrm{391}\;\textrm{nm}}}$ of the stimulated fluorescence signal in (2.1) and the gain coefficient of the air laser in the filament plasma in (2.6) would be influenced by the inversion fraction. Therefore, for a relativly low inversion fraction of the nitrogen molecular ion $\textrm{N}_2^ +$, both the air laser intensity and the number of the outer rings through the filament plasma would decrease.

It has been demonstrated that a remote pump laser at 226 nm wavelength could drive a two-photon dissociation of molecular oxygen and a two-photon excitation of one of the resulting oxygen atom fragments (Yuan et al. Reference Yuan, Wang, Lu, Chin and Zeng2014). The two-photon excitation transition from the ground state to the exited state produces fluorescence emission at 845 nm, which leads to the high gain and lasing properties in the backward direction. Actually, the lasing effect could be realized by the excitation of all three main components of air with a particular pump femtosecond laser wavelength, including nitrogen, oxygen and argon. The air lasing based on the ionic nitrogen molecules has attracted particular attention due to the fruitful observed effects and the high concentration in air. Pumped by the femtosecond laser pulses at 800 nm, nitrogen ions give rise to a series of narrow emission lines at 391.4 nm, 427.8 nm and 471.2 nm. The $\textrm{N}_2^ +$ lasing at 391 nm associated with the ${B^2}\sum _u^ + \to {X^2}\sum _g^ +$ emission induced by the femtosecond laser filament is noteworthy among the air lasing emission lines for that the population transferring process from the ground state could be monitored. Moreover, the air lasing at 391 nm has been commonly observed in the experiments under certain conditions. Therefore, in this study, we mainly focus on the air lasing at the 391.4 nm wavelength based on the ionic nitrogen molecules. The air lasing process based on the oxygen molecules was not included in the simulation.

3. Results and discussion

The air laser intensity distributions for various filament plasma lengths were investigated. The peak nitrogen ion density was chosen to be 1017 cm−3, which could be obtained in the laser filament pumped by a tightly focused femtosecond laser. The results are indicated in figure 1. The filament lengths are 2 cm, 4 cm, 10 cm, 20 cm and 30 cm for figures 1(a), 1(b), 1(c), 1(d) and 1(e), respectively. Figures 1(a 1)–1(e 1) depict the spatial distributions of the nitrogen ion density along the propagation direction for different filament lengths. The centre of the filament plasma is at z = 30 cm. Figure 1(a 2–e 2) show the spatial distribution of the air laser beam intensity along the propagation direction. Figure 1(a 3–e 3) correspond to the air laser intensities on the cross-section at the right end of the filament, while figure 1(a4–e4) correspond to the beam intensity patterns on the cross-section at the distance of 30 cm from the right-hand end of the filament. Figure 1(a5–e5) refer to the divergent angular distributions of the air laser intensities.

Figure 1. (a1)–(e1) Spatial distributions of the filament-induced nitrogen ion density versus propagation distance for different filament lengths. (a2)–(e2) Spatial distributions of the air laser intensity along the propagation direction. (a3)–(e3) Transverse patterns of air laser intensity distributions at the right end of the filament. (a4)–(e4) Transverse patterns of air laser intensity distributions at the distance of 30 cm from the right-hand end of the filament. (a5)–(e5) Divergent angular distributions of the air laser intensity. The filament lengths are (a) 2 cm, (b) 4 cm, (c) 10 cm, (d) 20 cm and (e) 30 cm.

Figure 2. (a1)–(e1) Spatial distributions of the air laser intensity along the propagation direction for different nitrogen ion densities. (b2)–(e2) Transverse patterns of the air laser intensity at distance of 70 cm. (b3)–(e3) Divergent angular distributions of the air laser intensity. (b4)–(e4) The air laser intensity at the centre of the beam cross-section as a function of propagation distance. The peak nitrogen ion densities are (a) 1016 cm−3, (b) 5 × 1016 cm−3, (c) 1017 cm−3, (d) 1.2 × 1017 cm−3 and (e) 1.4 × 1017 cm−3.

Figure 3. (a) Air laser beam power along the propagation direction for filament length of 20 cm and peak nitrogen ion density of 1017 cm−3. (b) The laser beam power after the filament plasma as a function of the filament length. (c) The laser beam power after the filament plasma as a function of the peak nitrogen ion density.

For the filament length of 2 cm, figures 1(a4) and 1(a5) indicate that a Gaussian profile of the air laser with a divergence angle of around 2 mrad is formed. When the filament length reaches 4 cm, a doughnut beam pattern of the air laser occurs after propagating through the filament at z = 32 cm as shown in figure 1(b 3). Figures 1(b4) and 1(b5) show that the air laser beam pattern can be regarded as the superposition of an inner Gaussian beam and an outer doughnut ring beam, resulting from the defocusing effect of the filament plasma and diffraction effect in the free space for the air laser. For the filament length of 10 cm, the inner Gaussian beam and the outer doughnut ring beam get apart from each other, as indicated in figure 1(c 3–c5). As shown in figure 1(d), when the filament length is 20 cm, a multiple rings structure occurs on the air laser beam cross-section after passing through the filament. This could be attributed to an interesting multiple-diffraction effect for the air laser inside the filament plasma. At the beginning of the air laser propagating in the filament plasma, since the lasing effect of the stimulated fluorescence overcomes the diffraction effect induced by the filament plasma, the air laser intensity gets increased with the propagating distance. After the air laser passes through the filament for several centimetres, the gain effect and the diffraction effect get balanced, leading to the highest value of the air laser intensity. Then the air laser intensity starts to drop down as the diffraction effect overcomes the lasing effect. During this process, an outer ring structure forms on the cross-section of the air laser beam. For a relatively long filament length, this dynamic process could be repeated for several times, resulting in a multiple rings structure pattern of the air laser. This is the so-called multiple-diffraction phenomenon proposed in this paper. Further, for the filament length of 30 cm, the number of the outer rings gets larger as shown in figure 1(e 3). At the distance z = 75 cm, the outer rings merge together and the outline of the air laser intensity divergent angular distribution becomes smooth as indicated in figures 1(e4) and 1(e5).

For the lasing action around 391 nm self-seeded by the second harmonic around 400 nm generated by the laser pulse in the filament plasma, the backward stimulated emission was relatively weak compared with the forward one, which was demonstrated in the experiment (Liu et al. Reference Liu, Brelet, Point, Houard and Mysyrowicz2013). The second harmonic generation was attributed to the non-uniform plasma distribution driven by the ponderomotive force or an electric-field-induced third-order mixing, leading to the difference of the harmonic intensity between the forward and backward directions. Conversely, it has been demonstrated that the air laser at 845 nm with a high gain can be achieved in the near-infrared region by pumping with a focused ultraviolet laser. Moreover, the lasing is equally strong in the forward and backward directions (Dogariu et al. Reference Dogariu, Michael, Scully and Miles2011). In this study, it is assumed that the filament plasma is uniformly distributed. The lasing action is initiated by fluorescence emission whose intensity depends on the filament plasma distribution, leading to equally strong lasing in both forward and backward directions. Therefore, both backward and forward propagation of air lasing could be assumed in the simulation.

Further, the air laser intensity spatial distributions for different peak nitrogen ion densities were studied. The filament length was set to be 20 cm. The nitrogen ion densities are 1016 cm−3, 5 × 1016 cm−3, 1017 cm−3, 1.2 × 1017 cm−3 and 1.4 × 1017 cm−3 for figures 2(a), 2(b), 2(c), 2(d) and 2(e), respectively. Figure 2(a 1–e 1) show the spatial distributions of the air laser intensity along the propagation direction, while figure 2(b 2–e 2) represent the air laser beam patterns on the cross-section at propagation distance of 70 cm. The corresponding divergent angular distributions of the air laser intensity are depicted in figure 2(a 3–e 3). Figure 2(a4–e4) depict the air laser intensity at the centre of the cross-section as a function of propagation distance.

As shown in figures 2(a 2) and 2(a 3), for a relative low peak nitrogen ion density of 1016 cm−3, a super-Gaussian profile of the air laser occurs after the filament plasma with a divergent angle of around 2 mrad. Inside the filament, the air laser peak intensity increases due to the lasing effect and then decreases induced by the diffraction effect. Moreover, there exists only one peak as shown in figure 2(a4). The outer ring structure of the air laser is formed with the peak nitrogen ion density of 5 × 1016 as indicated in figure 2(b). There are two peaks for the air laser intensity at the centre of the cross-section along the propagation direction as shown in figure 2(b4), resulting from the multiple-diffraction effect induced by a relative high nitrogen ion density distribution. By increasing the peak nitrogen ion density of the filament plasma, both the number of the outer rings for the air laser and the number of the peaks of the central air laser intensity along the propagation direction increase as shown in figure 2(c–e). It is worth mentioning that the peak intensities of the air laser along the propagation direction for different nitrogen ion density distributions almost remain at a constant order of magnitude due to the balance between the lasing effect and the diffraction effect.

In a previous experimental study (Liu et al. Reference Liu, Brelet, Point, Houard and Mysyrowicz2013), the lasing radiation profiles around 391 nm exhibits a Gaussian pattern in the nitrogen gas with a relatively low pressure. Upon increase of nitrogen gas pressure, an outer ring structure surrounding the central spot appears. This result is consistent with the simulation result in figure 2, which shows the same variation of the lasing radiation profile with the increase of the nitrogen ion density. This experimental result can be clarified by the multiple-diffraction effect induced by the relative high nitrogen ion density. In Zeng et al. Reference Zeng, Zhao, Liu and Chin2014, angular distribution of the backscattered nitrogen fluorescence induced by femtosecond laser filamentation in air has been studied in the experiment. Air lasing emission at 357 nm also forms a multiple rings profile.

A representative result for the air laser beam power along the propagation direction is shown in figure 3(a). The filament length is chosen to be 20 cm, while the peak nitrogen ion density is 1017 cm−3. The result indicates that the air laser experiences a gain effect inside the filament plasma and the beam power remains constant after passing through the filament. Figure 3(b) represents the air laser beam power after the filament as a function of the filament length. The amplified air laser beam power increases heavily when the filament length is less than 10 cm. Since the propagation of the air laser inside the filament for one time of diffraction is almost within 10 cm, the gain effect overcomes the diffraction effect during this filament region. Thus, the air laser beam intensity experience increasing continuously inside the filament. When the filament length is larger than 10 cm, the increase slope of the amplified air laser power versus the filament length drops down. This results from the multiple-diffraction effect. Since in some parts of the filament the diffraction effect overcomes the gain effect, the air laser beam intensity cannot increase continuously.

In addition, the laser beam power after the filament plasma as a function of the peak nitrogen ion density has been investigated as shown in figure 3(c). The filament length is set to be 20 cm. Since the lasing effect overcomes the diffraction effect for a relatively low nitrogen ion density, the air laser power increases heavily when the peak nitrogen ion density is less than 5 × 1016 cm−3. Then the air laser power increases slowly with the peak nitrogen ion density larger than 5 × 1016 cm−3. When the peak nitrogen ion density is larger than 1.2 × 1017 cm−3, the air laser power increases almost exponentially versus the peak nitrogen ion density. This is caused by the fact that the lasing effect of the filament plasma plays a dominant role in the propagating of the air laser for a relatively large nitrogen ion density.

When the incident pulse energy was below the energy threshold for harmonic generation, external seeding pulse was required to initiate lasing action of $\textrm{N}_2^ +$. Lasing emission around 391 nm could also self-seeded by the second harmonic around 400 nm generated by the laser pulses in the filament plasma for the case with the incident laser energy exceeding the second harmonic generation threshold energy (Liu et al. Reference Liu, Brelet, Point, Houard and Mysyrowicz2013). Self-induced white light inside the filament covering the ${B^2}\sum _u^ + \to {X^2}\sum _g^ +$ transition was serving as the seed to be amplified, leading to the self-seeded lasing action at 428 nm in air. The white light seed pulse was generated through spectral broadening mechanisms such as self-phase modulation and self-steepening of the femtosecond laser pulse (Chu et al. Reference Chu, Li, Xie, Ni, Yao, Zeng, Zhang, Jing, Xu, Cheng and Xu2013). In this study, the air lasing was self-seeded by the spontaneous fluorescence emission with the field amplitude depending on the nitrogen ion density in the filament plasma, which has been included in the simulation. For external seed laser passing through the filament plasma, the air laser would also form a multiple rings structure. Moreover, the multiple rings size would rely on the nitrogen ion density distribution.

4. Conclusion

In a summary, the lasing action of the fluorescence emitted from the transmission of the first negative (0,0) band of nitrogen ions at 391 nm was simulated during the femtosecond laser filamentation. The air laser intensity spatial distributions show that the beam profile changes from a Gaussian or super-Gaussian shape to an outer rings structure by increasing the filament length or nitrogen ion density. This result meets well the experimental result in a previous study. The multiple-diffraction effect induced by the dynamical interaction between the lasing effect and diffraction effect of the filament plasma causes the outer rings structure pattern of the air laser. This study provides a feasible way to understand the mechanism of the air laser propagating properties and improve the energy conversion efficiency and directivity of remote air lasers, which would be significant for remote sensing applications based on femtosecond laser filamentation.

Acknowledgements

Editor V. Malka thanks the referees for their advice in evaluating this paper.

Funding

This work was supported by the National Natural Science Foundation of China (grant number 12004316).

Declaration of interests

The authors report no conflict of interest.

Data availability

Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

References

Akozbek, N., Scalora, M., Bowden, C.M. & Chin, S.L. 2001 White-light continuum generation and filamentation during the propagation of ultra-short laser pulses in air. Opt. Commun. 191, 353–362.Google Scholar
Bahati, E.M., Jureta, J.J., Belic, D.D., Cherkani-Hassani, H., Abdellahi, M.O. & Defrance, P. 2001 Electron impact dissociation and ionization of $\textrm{N}_2^ +$. J. Phys. B 34, 2963–2973.Google Scholar
Berge, L., Skupin, S., Nuter, R., Kasparian, J. & Wolf, J.-P. 2007 Ultrashort filaments of light in weakly ionized, optically transparent media. Rep. Prog. Phys. 70, 1633.Google Scholar
Brian, J. & Mitchell, A. 1990 The dissociative recombination of molecular ions. Phys. Rep. 186, 215–248.Google Scholar
Chauveau, S., Perrin, M.Y., Riviere, P. & Soufiani, A. 2002 Contributions of diatomic molecular electronic systems to heated air radiation. J. Quant. Spectrosc. Radiat. Transfer 72, 503–530.Google Scholar
Chin, S.L., Hosseini, S.A., Liu, W., Luo, Q., Theberge, F., Akozbek, N., Becker, A., Kandidov, V.P., Kosareva, O.G. & Schroeder, H. 2005 The propagation of powerful femtosecond laser pulses in opticalmedia: physics, applications, and new challenges. Can. J. Phys. 83, 863.Google Scholar
Chiron, A., Lamouroux, B., Lange, R., Ripoche, J.F., Franco, M., Prade, B., Bonnaud, G., Riazuelo, G. & Mysyrowicz, A. 1999 Numerical simulations of the nonlinear propagation of femtosecond optical pulses in gases. Eur. Phys. J. D 6, 383.Google Scholar
Chu, W., Li, G., Xie, H., Ni, J., Yao, J., Zeng, B., Zhang, H., Jing, C., Xu, H., Cheng, Y. & Xu, Z. 2013 A self-induced white light seeding laser in a femtosecond laser filament. Laser Phys. Lett. 11 (1), 015301.Google Scholar
Couairon, A. & Mysyrowicz, A. 2007 Femtosecond filamentation in transparent media. Phys. Rep. 441, 47.Google Scholar
Doering, J.P. & Yang, J. 1996 Comparison of the electron impact cross section for the $\textrm{N}_2^ +$ first negative (0,0) band (λ 3914 Å) measured by optical fluorescence, coincidence electron impact, and photoionization experiments. J. Geophys. Res. 101, 19723.Google Scholar
Dogariu, A., Michael, J.B., Scully, M.O. & Miles, R.B. 2011 High-gain backward lasing in air. Science 331, 442.Google Scholar
Dubietis, A., Tamosauskas, G., Fibich, G. & Ilan, B. 2004 Multiple filamentation induced by input-beam ellipticity. Opt. Lett. 29, 1126–1128.Google Scholar
Jain, A., Gupta, D.N. & Suk, H. 2019 Electron–ion recombination effect on electron acceleration by an intense laser pulse. IEEE Trans. Plasma. Sci. 47, 4891–4897.Google Scholar
Johnson, A.W. & Fowler, R.G. 1970 Measured lifetimes of rotational and vibrational levels of electronic states of N2. J. Chem. Phys. 53, 65–72.Google Scholar
Kandidov, V.P., Shlenov, S.A. & Kosareva, O.G. 2009 Filamentation of high-power femtosecond laser radiation. Quant. Electron. 39, 205.Google Scholar
Kartashov, D., Alisauskas, S., Andriukaitis, G., Pugžlys, A., Shneider, M.N., Zheltikov, A.M., Chin, S.L. & Baltuska, A. 2012 Free-space nitrogen gas laser driven by a femtosecond filament. Phys. Rev. A 86, 033831.Google Scholar
Kartashov, D. & Shneider, M.N. 2017 Femtosecond filament initiated, microwave heated cavity-free nitrogen laser in air. J. Appl. Phys. 121, 113303.Google Scholar
Kasparian, J. & Wolf, J.-P. 2008 Physics and applications of atmospheric nonlinear optics and filamentation. Opt. Express 16, 466.Google Scholar
Liu, W. 2014 Intensity clamping during femtosecond laser filamentation. Chin. J. Phys. 52, 465.Google Scholar
Liu, Y., Brelet, Y., Point, G., Houard, A. & Mysyrowicz, A. 2013 Self-seeded lasing in ionized air pumped by 800 nm femtosecond laser pulses. Opt. Express 21, 22791.Google Scholar
Liu, W. & Chin, S.L. 2007 Abnormal wavelength dependence of the self-cleaning phenomenon during femtosecond-laser-pulse filamentation. Phys. Rev. A 76, 013826.Google Scholar
Luo, Q., Liu, W. & Chin, S.L. 2003 Lasing action in air induced by ultrafast laser filamentation. Appl. Phys. B 76, 337.Google Scholar
Mehr, F.J. & Biondi, M.A. 1969 Electron temperature dependence of recombination of $\textrm{O}_2^ +$ and $\textrm{N}_2^ +$ ions with electrons. Phys. Rev. 181, 264–271.Google Scholar
Ni, J., Chu, W., Jing, C., Zhang, H., Zeng, B., Yao, J., Li, G., Xie, H., Zhang, C., Xu, H., Chin, S.-L., Cheng, Y. & Xu, Z. 2013 Identification of the physical mechanism of generation of coherent $\textrm{N}_2^ +$ emissions in air by femtosecond laser excitation. Opt. Express 21, 8746.Google Scholar
Sprangle, P., Penano, J., Hafizi, B., Gordon, D. & Scully, M. 2011 Remotely induced atmospheric lasing. Appl. Phys. Lett. 98, 211102.Google Scholar
Wang, Q., Ding, P., Wilkins, S.G., Athanasakis-Kaklamanakis, M., Zhang, Y., Liu, Z. & Hu, B. 2023 Theoretical treatment on externally-seeded superradiance from $\textrm{N}_2^ +$ in femtosecond laser filamentation in low-pressure nitrogen. Front. Phys. 10, 1090346.Google Scholar
Xu, H.L., Azarm, A., Bernhardt, J., Kamali, Y. & Chin, S.L. 2009 The mechanism of nitrogen fluorescence inside a femtosecond laser filament in air. Chem. Phys. 360 (1–3), 171–175.Google Scholar
Xue, J., Zhang, N., Guo, L., Zhang, Z., Qi, P., Sun, L., Gong, C., Lin, L. & Liu, W. 2022 Effect of laser repetition rate on the fluorescence characteristic of a long-distance femtosecond laser filament. Opt. Lett. 47, 5676–5679.Google Scholar
Yao, J., Chu, W., Liu, Z., Chen, J., Xu, B. & Cheng, Y. 2018 An anatomy of strong-field ionization-induced air lasing. Appl. Phys. B 124, 1–17.Google Scholar
Yao, J., Zeng, B., Xu, H., Li, G., Chu, W., Ni, J., Zhang, H., Chin, S.L., Cheng, Y. & Xu, Z. 2011 High-brightness switchable multiwavelength remote laser in air. Phys. Rev. A 84 (5), 051802.Google Scholar
Yuan, S., Wang, T., Lu, P., Chin, S.L. & Zeng, H. 2014 Humidity measurement in air using filament-induced nitrogen monohydride fluorescence spectroscopy. Appl. Phys. Lett. 104, 091113.Google Scholar
Yuan, S., Wang, T., Teranishi, Y., Sridharan, A., Lin, S.H., Zeng, H. & Chin, S.L. 2013 Lasing action in water vapor induced by ultrashort laser filamentation. Appl. Phys. Lett. 102, 224102.Google Scholar
Zeng, T., Zhao, J.Y., Liu, W. & Chin, S.L. 2014 Backward angular distribution of air lasing induced by femtosecond laser filamentation. Laser Phys. Lett. 11 (7), 075401.Google Scholar
Zhang, X., Lu, Q., Zhang, Z., Fan, Z., Zhou, D., Liang, Q., Yuan, L., Zhuang, S., Dorfman, K. & Liu, Y. 2021 Coherent control of the multiple wavelength lasing of $\textrm{N}_2^ +$: coherence transfer and beyond. Optica 8 (5), 668–673.Google Scholar
Zhou, D., Zhang, X., Lu, Q., Liang, Q. & Liu, Y. 2020 Time-resolved study of the lasing emission from high vibrational levels of $\textrm{N}_2^ +$. Chin. Opt. Lett. 18, 023201.Google Scholar
Figure 0

Figure 1. (a1)–(e1) Spatial distributions of the filament-induced nitrogen ion density versus propagation distance for different filament lengths. (a2)–(e2) Spatial distributions of the air laser intensity along the propagation direction. (a3)–(e3) Transverse patterns of air laser intensity distributions at the right end of the filament. (a4)–(e4) Transverse patterns of air laser intensity distributions at the distance of 30 cm from the right-hand end of the filament. (a5)–(e5) Divergent angular distributions of the air laser intensity. The filament lengths are (a) 2 cm, (b) 4 cm, (c) 10 cm, (d) 20 cm and (e) 30 cm.

Figure 1

Figure 2. (a1)–(e1) Spatial distributions of the air laser intensity along the propagation direction for different nitrogen ion densities. (b2)–(e2) Transverse patterns of the air laser intensity at distance of 70 cm. (b3)–(e3) Divergent angular distributions of the air laser intensity. (b4)–(e4) The air laser intensity at the centre of the beam cross-section as a function of propagation distance. The peak nitrogen ion densities are (a) 1016 cm−3, (b) 5 × 1016 cm−3, (c) 1017 cm−3, (d) 1.2 × 1017 cm−3 and (e) 1.4 × 1017 cm−3.

Figure 2

Figure 3. (a) Air laser beam power along the propagation direction for filament length of 20 cm and peak nitrogen ion density of 1017 cm−3. (b) The laser beam power after the filament plasma as a function of the filament length. (c) The laser beam power after the filament plasma as a function of the peak nitrogen ion density.