Self-focusing/Defocusing of Hermite-Sinh-Gaussian Laser Beam in Underdense Inhomogeneous Plasmas

Abstract

The self-focusing/defocusing of Hermite-sinh-Gaussian (HshG) laser beam in underdense inhomogeneous plasmas is studied by using higher-order approximation theory. It is found that Hermite mode index and the fluctuation of the periodic plasma density have a significant effect on the dielectric constant and laser beam self-focusing/self-defocusing. With the increase of mode index, the high-order HshG laser beam is beneficial to suppress self-focusing and enhance self-defocusing. In addition, the effects of decentered parameters, beam intensity, and plasma non-uniformity on self-focusing/self-defocusing are discussed.

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

There are a lot of non-linear phenomena in the interaction between intense laser and plasma, such as self-focusing/self-defocusing [1], filamentation [2], stimulated Raman scattering [3] and stimulated Brillouin scattering [4]. These nonlinearities have important effects on inertial confinement fusion [5] and laser particle accelerator [6]. In particular, plasma self-focusing can improve the yield of these applications by increasing the interaction length between the laser and the medium, but also increase the success of these applications by balancing the natural diffraction of the laser beam to transmit more Rayleigh lengths [7]. Therefore, the self-focusing generated by the interaction of laser pulses with plasma is a very worthy topic of study [8].

In interaction intense laser and plasms, the mechanism of self-focusing/defocusing mainly include collisional nonlinearity [9], ponderomotive nonlinearity [10, 11], relativistic nonlinearity [12]. Nanda and Kant [13] researched self-focusing of cosh-Gaussian laser beams in collisionless non-uniform magnetic plasmas. By comparing extraordinary mode and ordinary mode, it is found that increasing the plasma density ramp can make the former stronger self-focusing. Asgharnejad et al. [7] investigated self-focusing a Gaussian laser pulse in a weakly relativistic and ponderomotive regime inside a collisionless warm quantum plasma. Zhang et al. [14] analyzed relativistic self-focusing of Gaussian laser beam, the results show that relativistic effect can enhance the oscillation of dielectric function and increase the self-focusing. Abedi-Varaki and Jafari [15] have present that a decrease in plasma temperature and an increase in oblique magnetic field increases defocusing for left-handed polarization. Extending the High-harmonic generation phase matching cutoff by controlling the rapid self-defocusing effect of driving laser is experimentally and theoretically demonstrated by Sun et al. [16].

Lots of research of self-focusing/self-defocusing have focused on interaction of Gaussian beam and homogeneous/inhomogeneous plasmas. Studies by Gupta et al. [17] and Kant et al. [18] have shown that the existence of plasma density fluctuations can enhance self-focusing. Wani et al. [19] have found that electron temperature can reduce defocusing and plasma density transition increases the range of self-focusing effect. Recently, many researchers have analyzed the effects of different types of laser beams and plasma on self-focusing/self-defocusing. Laser beam extends from Gaussian type to the others type of beam, such as Hermite-Gaussian beam [20], cosh-Gaussian beam [21], super Gaussian beam [22]. Various types plasmas have been researched, such as cold quantum plasma with exponential density variation [23]. Effects of magnetic field on relativistic self-focusing, self-phase modulation and self-trapping of cosh-Gaussian laser beam is exploited by Gill et al. [24]. Kant et al. [25] observed relativistic self-focusing of Hermite-cosh-Gaussian laser pulse in plasma with exponential plasma density ramp. Devi and Malik [26] studied the relativistic self-focusing of super-Gaussian laser beams, they have found that the external magnetic field and the cyclotron motion of electrons is important for self-focusing. Based on the moment theory, Malik and Devi [27] developed a formulation to study the self-defocusing of super-Gaussian laser beam in tunnel ionized plasmas, and found that laser beam can get more defocused for higher intensity and smaller spot size. Kaur et al. [28] used paraxial approximation method to study propagation characteristics of Hermite-cosh-Gaussian laser beam in rippled density plasmas and strong oscillatory self-focusing and defocusing are observed in relativistic case.

In this paper, applying higher-order paraxial theory and WKB approximation, the self-focusing/self-defocusing are investigated in the interaction of Hermite-Sinh-Gaussian (HshG) laser beam and the non-uniform underdense plasma with periodic density fluctuations.

2. Propagation Theory of Laser Beam in Plasmas

In laser-plasma interaction, according to Maxwell’s equation theory, the wave equation of laser propagation in plasma can be expressed as,$$\begin{array}{}\text{(1)}& {\nabla }^{2}E-\frac{{\omega }^2}{c^2}\epsilon E+\nabla \frac{E\cdot \nabla \epsilon }{\epsilon }=0.\end{array}$$

For $k^{-2}{\nabla }^2\ln \epsilon <<1$ and $k^2=\epsilon {\omega }^2/c^2$, the wave equation can be converted to $$\begin{array}{}\text{(2)}& \frac{{\partial }^{2}E}{\partial z^2}+\frac{{\partial }^{2}E}{\partial r^2}+\frac{1}{r}\frac{\partial E}{\partial r}+\frac{{\omega }^2}{c^2}\epsilon E=0.\end{array}$$

When the laser propagates in the plasma along the z-direction, the specific form of the field can be expressed as,$$\begin{array}{}\text{(3)}& E=Ar,z{e}^{-i\omega t-kz}.\end{array}$$

In this paper, the interaction between HshG laser beam and inhomogeneous plasma is studied. For equation (3), the general solution of the field amplitude is HshG beam.$$\begin{array}{}\text{(4)}& Er,z=\frac{E_0}{fz}{H}_n\frac{\sqrt{2}r}{r_{0}fz}{e}^{b^2/4}{e}^{-{r/r_{0}f+b/2}^2}-e^{-{r/r_{0}f-b/2}^2}.\end{array}$$$E_0$ is the amplitude of HshG laser beam at $r=z=0$, Hermite polynomial $H_n\sqrt{2}r/r_{0}fz$, when $n=0$, $H_n\sqrt{2}r/r_{0}fz=1$; when $n=1$, $H_n\sqrt{2}r/r_{0}fz=2r$; when $n=2$, $H_n\sqrt{2}r/r_{0}fz=4r^2-2$. $r$ is the radial component in the cylindrical coordinate system, $r_0$ is the initial beam width, $b$ represents decentred parameter, the dimensionless beam width parameter of laser beam is $fz$.

Non-uniform plasmas showing periodic variations have been found experimentally and theoretically by researchers in the process of laser-plasma interaction. The plasma system studied in this paper is based on the experimental results of Ouahid et al. [29]; electron density $N_{e}z$ is sinusoidal function. The expression is as follows,$$\begin{array}{}\text{(5)}& N_{e}z=N_{0e}1+\alpha \sin qz,\end{array}$$where $N_{e}z$, $N_{0e}$, $\alpha$, $q$ represents electron density in equilibrium steady state, initial electron density, parameter of plasma fluctuation, wave number of plasma density fluctuation.

In laser-plasma interaction, when plasma density fluctuation exists, the dielectric function of plasma composed of linear part ${\epsilon }_0$ and non-linear part $\varphi E{E}^{\ast }$ can be expressed as $\epsilon r,z={\epsilon }_0-\varphi E{E}^{\ast }$. In the limit $c^2/r_0^2{\omega }_{p0}^2<1$, the linear and non-linear part can be presented as follows ${\epsilon }_{0}z=1-{\omega }_{p0}^{2}1+\alpha \sin qz/\gamma {\omega }^2$, $\varphi E{E}^{\ast }={\omega }_{p0}^{2}1+\alpha \sin qz1-{1/1+\beta E{E}^{\ast }}^{1/2}/\gamma {\omega }^2$. ${\omega }_{p0}=\sqrt{4\pi n_0{e}^2/m}$ is the plasma frequency, $e$ is the electronic charge and $m$ shows the rest mass. Here the relativistic factor can be written as $\gamma =\sqrt{1+\beta E{E}^{\ast }/2}$, where $\beta ={e/m\omega c}^2$. By applying higher-order paraxial theory and extending the dielectric function up to the fourth higher-order term.$$\begin{array}{}\text{(6)}& \epsilon r,z={\epsilon }_{0}z-{\epsilon }_{2}z\frac{r^2}{r_0^2}-{\epsilon }_{4}z\frac{r^4}{r_0^4},\end{array}$$where ${\epsilon }_{0}z$ is the linear part, ${\epsilon }_{2}z$ and ${\epsilon }_{4}z$ are expansion coefficients of $\epsilon r,z$. Based on the research results of Pathak et al. [30]; using WKB approximation and higher-order paraxial theory, dielectric constant components of different Hermite mode indexes can be obtained in interaction of HshG beams and plasma.For mode index $n=0$$$\begin{array}{}\text{(7)}& \begin{array}{c}{\epsilon }_0=1-t+\frac{3t}{4}\frac{B}{f^2}{e}^{b^2/2}D-\frac{27t}{32}\frac{B^2}{f^4}{e}^{b^2}{D}^2,\hfill \\ {\epsilon }_2=t\frac{B^2}{f^6}{e}^{b^2}\frac{3}{32}4b^2-36D+\frac{27}{16}{a}_2{D}^2-\frac{tB}{4f^4}{e}^{b^2/2}7b^2-6+3a_{2}D,\hfill \\ \begin{array}{c}{\epsilon }_4=\frac{t{B}^2}{f^8}{e}^{b^2}\frac{3}{32}72+\frac{107}{2}{b}^4-102b^2+\frac{27D^2}{32}{a}_2^2+2a_4+\frac{3D}{16}{a}_{2}42b^2-36\\ -\frac{tB}{4f^6}{e}^{b^2/2}6+a_{2}7b^2-6+3a_{4}D.\end{array}\hfill \end{array}\end{array}$$

For mode index $n=1$$$\begin{array}{}\text{(8)}& \begin{array}{c}{\epsilon }_0=1-t,\hfill \\ {\epsilon }_2=-\frac{6tB}{f^4}{e}^{b^2/2}D,\hfill \\ {\epsilon }_4=\frac{54t{B}^2}{f^8}{e}^{b^2}{D}^2-\frac{2tB}{f^6}{e}^{b^2/2}7b^2-6+3a_{2}D.\hfill \end{array}\end{array}$$

For mode index $n=2$$$\begin{array}{}\text{(9)}& \begin{array}{c}{\epsilon }_0=1-t+\frac{3tB}{f^2}{e}^{b^2/2}D-\frac{27t}{2}\frac{B^2}{f^4}{e}^{b^2}{D}^2,\hfill \\ {\epsilon }_2=\frac{9t}{4}\frac{B^2}{f^6}{e}^{b^2}47b^2-6-96D+12a_{2}D-\frac{tBD}{4f^4}{e}^{b^2/2}47b^2-6-96D+12a_{2}D,\hfill \\ \begin{array}{c}{\epsilon }_4=\frac{3t{B}^2}{32f^8}{e}^{b^2}{47b^2-6-96D+12a_{2}D}^2+\frac{9t{B}^{2}D}{4f^8}{e}^{b^2}24+192-96a_2+12a_{4}D+4a_2-327b^2-6\\ -\frac{tB}{4f^6}{e}^{b^2/2}24+192-96a_2+12a_{4}D+4a_2-327b^2-6,\end{array}\hfill \end{array}\end{array}$$where in equations (7)–(9), $t={\omega }_{p0}^2/{\omega }^{2}1+\alpha \sin q\xi$ denotes the parameters related to plasma density, $B=\beta E_0^2$, $D=1-b^2/2+b^4/8$ denote intensity dependence and intensity distribution dependence of dielectric constant. For wave equation (2), substitute (3) into equation (2), ${\partial }^{2}A/\partial z^2$ be ignored, can get,$$\begin{array}{}\text{(10)}& -2ik\frac{\partial A}{\partial z}-iA\frac{\partial k}{\partial z}+\frac{{\partial }^{2}A}{\partial r^2}+\frac{1}{r}\frac{\partial A}{\partial r}+\frac{{\omega }^2}{c^2}\epsilon r,z-{\epsilon }_{0}zA=0.\end{array}$$

Following Nanda et al. [31]; complex amplitude A can be written as follows,$$\begin{array}{}\text{(11)}& Ar,z=A_{0}r,z{e}^{-ikSr,z},\end{array}$$$A_{0}r,z$ and $Sr,z$ are the real functions of space, according to Habibi and Ghamari [32]; they can be expressed this way$$\begin{array}{}\text{(12)}& A_0^2=\frac{E_0^2}{f^2}{e}^{b^2/2}{H_n\frac{\sqrt{2}r}{r_{0}f}}^{2}1+\frac{a_2{r}^2}{r_0^2{f}^2}+\frac{a_4{r}^4}{r_0^4{f}^4}{e}^{-2r/r_{0}f+b^2/2}-e^{-2r^2/r_0^2{f}^2+b^2/2},\text{(13)}& Sr,z=S_{0}z+\frac{r^2}{r_0^2}{S}_{2}z+\frac{r^4}{r_0^4}{S}_{4}z.\end{array}$$

$a_2,a_4$, $S_0,S_2,S_4$ are the functions of $z$, $S_0,S_2,S_4$ are the components of the eikonal. $S_2=r_0^2/2fdf/d\xi$ corresponds to the spherical curvature of the wave front, $S_4$ indicates deviation from spherical curvature. $a_2,a_4$ are indicative of the departure of the beam from the HshG nature. Can get a complicated differential equation by bringing in equation (10) with equation (11), the following equations can be obtained if the real and imaginary parts are equal.$$\begin{array}{}\text{(14)}& 2\frac{\partial S}{\partial z}+\frac{2S}{k}\frac{\partial k}{\partial z}+{\frac{\partial S}{\partial r}}^2=\frac{1}{k^2{A}_0}\frac{{\partial }^2{A}_0}{\partial r^2}+\frac{1}{r}\frac{\partial A_0}{\partial r}-\frac{r^2}{r_0^2}\frac{{\epsilon }_{2}z}{{\epsilon }_{0}z}-\frac{r^4}{r_0^4}\frac{{\epsilon }_{4}z}{{\epsilon }_{0}z},\text{(15)}& \frac{\partial A_0^2}{\partial z}+A_0^2\frac{{\partial }^{2}S}{\partial r^2}+\frac{1}{r}\frac{\partial S}{\partial r}+\frac{\partial A_0^2}{\partial r}\frac{\partial S}{\partial r}+\frac{A_0^2}{k}\frac{dk}{dz}=0.\end{array}$$

By substituting equations (12) and (13) in equation (15) for different Hermite mode indexes, and equating the coefficients of $r^2$ and $r^4$ on both sides of equation, the equation relation of parameter $a_2$ and $a_4$ is represented as.

For mode index $n=0$$$\begin{array}{}\text{(16a)}& \frac{d{a}_2}{d\xi }=-16S_4^{\prime }{f}^2,\text{(16b)}& \frac{d{a}_4}{d\xi }=-24a_2+44-2b^2{S}_4^{\prime }{f}^2.\end{array}$$

For mode index $n=1$$$\begin{array}{}\text{(17a)}& \frac{d{a}_2}{d\xi }=-24S_4^{\prime }{f}^2,\text{(17b)}& \frac{d{a}_4}{d\xi }=-32a_2+44-2b^2{S}_4^{\prime }{f}^2.\end{array}$$

For mode index $n=2$$$\begin{array}{}\text{(18a)}& \frac{d{a}_2}{d\xi }=-16S_4^{\prime }{f}^2,\text{(18b)}& \frac{d{a}_4}{d\xi }=-24a_2+64+44-2b^2{S}_4^{\prime }{f}^2,\end{array}$$where $S_4^{\prime }=S_4\omega /c$. Same method, equations (12) and (13) in equation (14) to get the following equations governing the beam width $f$ and $S_4^{\prime }$.

For mode index $n=0$$$\begin{array}{c}\text{(19a)}& \frac{d^{2}f}{d{\xi }^2}=\frac{-2a_2^2+8-8a_2-46b^4/3+4a_2{b}^2-8b^2+8a_4}{{\epsilon }_0{f}^3}-\frac{{\epsilon }_{2}f{\rho }^2}{{\epsilon }_0}-\frac{1}{2{\epsilon }_0}\frac{d{\epsilon }_0}{d\xi }\frac{df}{d\xi },\text{(19b)}& \frac{d{S}_4^{\prime }}{d\xi }=\frac{8a_2{a}_4-4a_2^3-52a_2{b}^4/3-16a_4+8a_2^2+104b^4/3+8b^2{a}_4-4a_2^2{b}^2-52b^6/3}{2{\epsilon }_0{f}^6}-\frac{4S_4^{\prime }}{f}\frac{df}{d\xi }-\frac{S_4^{\prime }}{2{\epsilon }_0}\frac{d{\epsilon }_0}{d\xi }-\frac{{\epsilon }_4{\rho }^2}{2{\epsilon }_0}.\end{array}$$

For mode index $n=1$$$\begin{array}{}\text{(20a)}& \frac{d^{2}f}{d{\xi }^2}=\frac{8-6a_2^2-8a_2-98b^4/3+4a_2{b}^2-8b^2+16a_4}{{\epsilon }_0{f}^3}-\frac{{\epsilon }_{2}f{\rho }^2}{{\epsilon }_0}-\frac{1}{2{\epsilon }_0}\frac{d{\epsilon }_0}{d\xi }\frac{df}{d\xi },\text{(20b)}& \frac{d{S}_4^{\prime }}{d\xi }=\frac{2a_2-4+2b^{2}4a_4-2a_2^2-26b^4/3}{2{\epsilon }_0{f}^6}-\frac{4S_4^{\prime }}{f}\frac{df}{d\xi }-\frac{S_4^{\prime }}{2{\epsilon }_0}\frac{d{\epsilon }_0}{d\xi }-\frac{{\epsilon }_4{\rho }^2}{2{\epsilon }_0}.\end{array}$$

For mode index $n=2$$$\begin{array}{}\text{(21a)}& \frac{d^{2}f}{d{\xi }^2}=\frac{-2a_2^2+72-40a_2-46b^4/3+4a_2{b}^2-40b^2+8a_4}{{\epsilon }_0{f}^3}-\frac{{\epsilon }_{2}f{\rho }^2}{{\epsilon }_0}-\frac{1}{2{\epsilon }_0}\frac{d{\epsilon }_0}{d\xi }\frac{df}{d\xi },\text{(21b)}& \frac{d{S}_4^{\prime }}{d\xi }=\frac{2a_2-20+2b^{2}4a_4-64-2a_2^2-26b^4/3}{2{\epsilon }_0{f}^6}-\frac{4S_4^{\prime }}{f}\frac{df}{d\xi }-\frac{S_4^{\prime }}{2{\epsilon }_0}\frac{d{\epsilon }_0}{d\xi }-\frac{{\epsilon }_4{\rho }^2}{2{\epsilon }_0}.\end{array}$$

Normalize the above equations, $\xi =cz/\omega r_0^2=z/R_d$ is the dimensionless propagation distance, $\rho =\omega r_0/c$ represents the dimensionless original beam width.

3. Discussion

In the paper, intensity of incident laser beam is $1.23\times {10}^{17}\text{W}/{\text{cm}}^2$, $\beta E_0^2=0.01,0.02,0.1$ and 0.2 (corresponding intensities are $1.23\times {10}^{15}\text{W}/{\text{cm}}^2$, $2.46\times {10}^{15}\text{W}/{\text{cm}}^2$, $1.23\times {10}^{16}\text{W}/{\text{cm}}^2$ and $2.46\times {10}^{16}\text{W}/{\text{cm}}^2$ respectively), the frequency and spot size of incident laser beam as $\omega =1.78\times {10}^{15}\text{rad}/\text{s}$, $r_0=1.68\times {10}^{-5}\text{m}$ in interaction system. Appling fourth-order Runge-Kutta method to solve equations (16a)–(21b), the phenomena of self-focusing/self-defocusing in the interaction of intense laser and the inhomogeneous plasma are discussed. The initial conditions are $f=1$, $df/d\xi =0$, $a_2=1$, as well as $d{a}_2/d\xi =0$ at $\xi =0$, plasmas density is $N_{e}z=N_{0e}1+\alpha \sin qz$.

Figures 1(a)–1(c) shows the linear component of the dielectric function of the HshG laser beam interaction of plasma with periodic density fluctuation at mode index n = 0, 1 and 2. The linear part of dielectric function shows sinusoidal periodic variations similar to oscillation in plasma of periodic density fluctuation. The peak value of the linear part decreases slightly with the increase of mode index, and the oscillation has intensified the tendency. For the non-linear part ${\epsilon }_2$, ${\epsilon }_4$ at different mode indexes, the value is much smaller compared to the linear part, and mode index has a more obvious influence on the non-linear part. In figure (d), (e) and (f), it can be seen that mode index obvious affect the peak value of non-linear part, the peak position drifts to the left as the mode increases, and finally all tend to stable. From the above linear and non-linear parts of the dielectric function, it shows an approximate periodic oscillatory variation due to the plasma periodicity, and its oscillatory variation trend to increase with the increase of mode index of Hermite polynomial. In addition, the mode index have more obvious effect on the non-linear part comparison linear part.

[figure(s) omitted; refer to PDF]

Figure 2 presents beam width variation for different Hermite mode indexes of HshG laser beams. In the periodic underdense inhomogeneous plasma, the beam width all show oscillatory changes, firstly, beam width becomes small and then increases. During the propagation of the laser beam, ponderomotive force influence on laser beam is obvious at the start, and beam width reduces and appears to convergence, resulting in self-focusing. Then, the influence ponderomotive force of laser beam decrease and beam width increase continuously with distance, and beam width presents diverge obviously, resulting in beam self-defocusing. However, at different mode indexes, the size and position of the minimum value of beam width are different. Increase of mode index, the size slightly increase and position appear early, and the beam diverges more quickly, which means the self-defocusing is more intense. For $n=0$, reaches the minimum value 0.55 at $\xi \approx 0.14$. For $n=1$, reaches the minimum value 0.6 at $\xi \approx 0.1$. For $n=2$, reaches the minimum value 0.998 at $\xi \approx 0.02$. This is due to the combined effects of HshG laser beam and periodic underdense non-uniform plasma, causing a similar periodic oscillatory variation of dielectric function, which leads to self-focusing and self-defocusing. In addition, with the increase of mode index, the oscillation of the dielectric function is enhanced, indicating that the instability is enhanced, which will affect the coupling of the laser plasma, so that the quality of Gaussian beam decreases, self-focusing reduce and self-defocusing enhance.

[figure(s) omitted; refer to PDF]

Figure 3 presents the effect of decentered parameters b on beam width $f$ in different modes indexes. It is shown that the decentered parameter has obvious influence on the beam width. For $n=0$, the beam width slowly increases in $b=0$, and in $b\ne 0$, it first decreases, then increases, and with the increase of $b$, the variation of beam width is strengthened. $n=1$ is similar to $n=0$. For (e) $n=2$, beam width increases slightly at first, and then almost linearly increases, moreover, beam width decreases slightly for the same propagation distance with increase of $b$. (a) (c) correspond to (b) (d), the rules are consistent, the purpose is to study the effect of smaller changes in b on the self-focusing, and optimize the value of the decentered parameter. (f) depicts the variation of the beam width when $\beta E_0^2=0.3$, the beam width shows an oscillation trend. Compared with (e), it not only shows the influence of decentered parameters on self-focusing, but also reflects the importance of laser intensity to self-focusing. In summary, (b) (d) (f) show that the decentered parameter is very sensitive to the self-focusing, which is consistent with that described in the paper of Nanda et al. [33]. For the zero-order and 1-order HshG beams, under the influence of the decentered parameter, there is a significant beam self-focusing at $b\ne 0$, however, for the 2-order HshG beam, it presents obvious beam self-defocusing, in this case, the effect of decentered parameter on self-defocusing is relatively slight. It can be clearly seen that the effect of decentered parameters on self-focusing of the HshG beam is significantly stronger than self-defocusing.

[figure(s) omitted; refer to PDF]

Figure 4 depicts the variation of the beam width for different light intensities. The figure is shown that, beam width increases at the same propagation distance with the increase of light intensity, it implies that self-focusing decrease and self-defocusing increase. In addition, it is easier to form self-focusing with the decrease of laser intensity, while self-focusing becomes weak when the light intensity increases to 0.2, it may be because laser intensity is close or exceed the threshold value of forming self-focusing, leading to self-focusing reduce or disappear. In Figure 4(c), for $n=2$ , when laser intensity increases, beam width continues to increase, self-focusing is very weak and produces obvious self-defocusing. Hence, it is shown that the high-order HshG beam and high beam intensity are benefit to restrain self-focusing and enhance self-defocusing.

[figure(s) omitted; refer to PDF]

Figure 5 depicts the variation of the beam width for different plasma density parameter. With the increase of plasma inhomogeneity, the beam width increases at $n=0$ and $n=1$. At $n=2$, as the plasma inhomogeneity increases, beam width slowly decreases first and then increases. It implies that the inhomogeneity of plasma has a significant effect on the beam width because the plasma density fluctuation affects the dielectric function, which directly affects laser propagation process in plasma, lead to the change of beam width, and increase non-uniformity of plasmas, it can strength self-focusing/self-defocusing.

[figure(s) omitted; refer to PDF]

4. Conclusion

In the paper, apply high-order paraxial theory to analyze the self-focusing/self-defocusing of HshG laser beam in underdense inhomogeneous plasma. Self-focusing and self-defocusing laws are explored by analyzing the variation of the dielectric constant and the influence of the decentered parameters, laser intensity and plasma density on beam width. The results show that, due to the combined effects of HshG laser beam and periodic underdense non-uniform plasma, lead to dielectric function present a similar periodic oscillation variation. Under the influence of ponderomotive force and diffraction effect, lead to beam self-focusing/self-defocusing. With increase of Hermite mode index, causing the position of self-focusing/self-defocusing is bring forward, and beam defocusing is strengthened, it implies that the higher mode index HshG beam is benefit to reduce self-focusing and enhance defocusing. In addition, the decentered parameter, beam intensity and non-uniformity of the plasma has a significant influence on the beam width. Influence of the decentered parameter on beam self-focusing is stronger than self-defocusing. In the high-order HshG beam, high beam intensity is benefit to restrain self-focusing and enhance self-defocusing, and increase non-uniformity of plasmas, it can strength self-focusing/self-defocusing.

Acknowledgments

The work is supported by the National Natural Science Foundation of China (Grant no. 11447169) and the Natural Science Foundation of Guangxi province (2018GXNSFAA138180, 2016GXNSFAA380071, and 2016GXNSFBA380204).

References

[1] V. F. Kovalev, V. Y. Bychenkov, "Analytic theory of relativistic self-focusing for a Gaussian light beam entering a plasma: renormalization-group approach," Physical Review, vol. 99,DOI: 10.1103/PhysRevE.99.043201, 2019.

[2] M. B. Hassan, I. J. Abd-Ali, A. O. Soary, "The filamentation instability of nonparaxial laser beam inside magnetized plasma," Results in Physics, vol. 14,DOI: 10.1016/j.rinp.2019.102386, 2019.

[3] H. Peng, Y. L. Zuo, H. Y. Zhu, J. Q. Su, "Forward Raman scattering of the seed pulse in strongly coupled stimulated Brillouin amplification in plasma," Physics of Plasmas, vol. 25 no. 1,DOI: 10.1063/1.5020327, 2018.

[4] Z. Li, Y. Zuo, J. Su, S. Yang, "The filamentation effect in short pulse amplification by strong-coupling stimulated Brillouin scattering," Physics of Plasmas, vol. 26 no. 9,DOI: 10.1063/1.5094513, 2019.

[5] C. Lu, V. Tikhonchuk, S. Weber, "Analytic solutions for delocalized heat transport," Plasma Physics and Controlled Fusion, vol. 63 no. 7,DOI: 10.1088/1361-6587/abf766, 2021.

[6] K. V. Lotov, A. P. Sosedkin, A. V. Petrenko, L. D. Amorim, J. Vieira, R. A. Fonseca, L. O. Silva, E. Gschwendtner, P. Muggli, "Electron trapping and acceleration by the plasma wakefield of a self-modulating proton beam," Physics of Plasmas, vol. 21 no. 12,DOI: 10.1063/1.4904365, 2014.

[7] D. Asgharnejad, T. Mohsenpour, S. Mirzanejhad, "Investigation of self-focusing a Gaussian laser pulse in a weakly relativistic and ponderomotive regime inside a collisionless warm quantum plasma," Chinese Journal of Physics, vol. 73, pp. 304-313, DOI: 10.1016/j.cjph.2021.06.018, 2021.

[8] X. Xia, B. Xu, "Off-axial contribution of beam self-focusing in plasma with density ripple," Optik, vol. 125 no. 19, pp. 5899-5903, DOI: 10.1016/j.ijleo.2014.07.044, 2014.

[9] M. Singh, D. N. Gupta, "Laser-pulse compression in a collisional plasma under weak-relativistic ponderomotive nonlinearity," Physics of Plasmas, vol. 23 no. 5,DOI: 10.1063/1.4951722, 2016.

[10] S. D. Patil, P. P. Chikode, M. V. Takale, "Turning point temperature of self-focusing at laser-plasma interaction with weak relativistic-ponderomotive nonlinearity: effect of light absorption," Journal of Optics, vol. 47 no. 2, pp. 174-179, DOI: 10.1007/s12596-018-0448-z, 2018.

[11] P. Rawat, G. Purohit, "Self-focusing of a cosh-Gaussian laser beam in magnetized plasma under relativistic-ponderomotive regime," Contributions to Plasma Physics, vol. 59 no. 2, pp. 226-235, DOI: 10.1002/ctpp.201800066, 2019.

[12] A. Sharma, V. K. Tripathi, "Relativistic and ponderomotive self-focusing of a laser pulse in magnetized plasma," Laser and Particle Beams, vol. 30 no. 4, pp. 659-664, DOI: 10.1017/s0263034612000481, 2012.

[13] V. Nanda, N. Kant, "Strong self-focusing of a cosh-Gaussian laser beam in collisionless magneto-plasma under plasma density ramp," Physics of Plasmas, vol. 21 no. 7,DOI: 10.1063/1.4889862, 2014.

[14] G. Zhang, Q. Liang, X. Xia, "Relativistic self-focusing in the interaction of laser beam and plasma with periodical density ripple," Laser and Particle Beams, vol. 38 no. 4, pp. 244-250, DOI: 10.1017/s0263034620000300, 2020.

[15] M. Abedi-Varaki, S. Jafari, "Nonlinear interaction of intense left- and right-hand polarized laser pulse with hot magnetized plasma," Journal of Plasma Physics, vol. 83 no. 4,DOI: 10.1017/s0022377817000460, 2017.

[16] H.-W. Sun, P.-C. Huang, Y.-H. Tzeng, J.-T. Huang, C. D. Lin, C. Jin, M.-C. Chen, "Extended phase matching of high harmonic generation by plasma-induced defocusing," Optica, vol. 4 no. 8, pp. 976-981, DOI: 10.1364/optica.4.000976, 2017.

[17] D. N. Gupta, M. S. Hur, H. Suk, "Additional focusing of a high-intensity laser beam in a plasma with a density ramp and a magnetic field," Applied Physics Letters, vol. 91 no. 8,DOI: 10.1063/1.2773943, 2007.

[18] N. Kant, M. A. Wani, A. Kumar, "Self-focusing of Hermite-Gaussian laser beams in plasma under plasma density ramp," Optics Communications, vol. 285 no. 21-22, pp. 4483-4487, DOI: 10.1016/j.optcom.2012.05.065, 2012.

[19] M. A. Wani, V. Thakur, H. S. Ghotra, N. Kant, "Effect of axial electron temperature and plasma density ramp on self-focusing/defocusing of a laser beam in plasma," Optik, vol. 192,DOI: 10.1016/j.ijleo.2019.162963, 2019.

[20] J. Wadhwa, A. Singh, "Second harmonic generation of self-focused Hermite-Gaussian laser beam in collisional plasma," Optik, vol. 202,DOI: 10.1016/j.ijleo.2019.01.116, 2020.

[21] T. U. Urunkar, S. D. Patil, A. T. Valkunde, B. D. Vhanmore, K. M. Gavade, M. V. Takale, "Effect of critical beam radius on self-focusing of cosh-Gaussian laser beams in collisionless magnetized plasma," Communications in Theoretical Physics, vol. 70 no. 2, pp. 220-224, DOI: 10.1088/0253-6102/70/2/220, 2018.

[22] H. K. Malik, L. Devi, "Relativistic self focusing and frequency shift of super-Gaussian laser beam in plasma," Results in Physics, vol. 17,DOI: 10.1016/j.rinp.2020.103070, 2020.

[23] V. Thakur, N. Kant, "Stronger self-focusing of a chirped pulse laser with exponential density ramp profile in cold quantum magnetoplasma," Optik, vol. 172, pp. 191-196, DOI: 10.1016/j.ijleo.2018.07.027, 2018.

[24] T. S. Gill, R. Kaur, R. Mahajan, "Relativistic self-focusing and self-phase modulation of cosh-Gaussian laser beam in magnetoplasma," Laser and Particle Beams, vol. 29 no. 2, pp. 183-191, DOI: 10.1017/s0263034611000152, 2011.

[25] N. Kant, S. Vij, S. K. Chakravarti, J. P. Kushwaha, V. Thakur, "Relativistic self-focusing of Hermite-cosh-Gaussian laser beam in magnetoplasma with exponential plasma density ramp," Communications in Theoretical Physics, vol. 71 no. 12, pp. 1469-1474, DOI: 10.1088/0253-6102/71/12/1469, 2019.

[26] L. Devi, H. K. Malik, "Role of magnetic field on self focusing of super-Gaussian laser beam under relativistic effect," Optik, vol. 207,DOI: 10.1016/j.ijleo.2020.164439, 2020.

[27] H. K. Malik, L. Devi, "Self-defocusing of super-Gaussian laser beam in tunnel ionized plasmas," Optik, vol. 222,DOI: 10.1016/j.ijleo.2020.165357, 2020b.

[28] R. Kaur, T. Singh Gill, R. Mahajan, "Relativistic effects on evolution of a q-Gaussian laser beam in magnetoplasma: application of higher order corrections," Physics of Plasmas, vol. 24 no. 5,DOI: 10.1063/1.4983309, 2017.

[29] L. Ouahid, L. Dalil-Essakali, A. Belafhal, "Evolution of the beam-width parameter of zeroth-order Bessel-Gaussian beams in collisional plasma with density ripple," Optical and Quantum Electronics, vol. 51 no. 4,DOI: 10.1007/s11082-019-1818-8, 2019.

[30] N. Pathak, M. Kaur, S. Kaur, T. S. Gill, "Non-paraxial theory of self-focusing/defocusing of Hermite cosh Gaussian laser beam in rippled density plasmas," Contributions to Plasma Physics, vol. 59 no. 10,DOI: 10.1002/ctpp.201900026, 2019.

[31] V. Nanda, H. S. Ghotra, N. Kant, "Early and strong relativistic self-focusing of cosh-Gaussian laser beam in cold quantum plasma," Optik, vol. 156, pp. 191-196, DOI: 10.1016/j.ijleo.2017.10.147, 2018.

[32] M. Habibi, F. Ghamari, "Improved focusing of a cosh-Gaussian laser beam in quantum plasma: higher order paraxial theory," IEEE Transactions on Plasma Science, vol. 43 no. 7, pp. 2160-2165, DOI: 10.1109/tps.2015.2440319, 2015.

[33] V. Nanda, N. Kant, M. A. Wani, "Sensitiveness of decentered parameter for relativistic self-focusing of Hermite-cosh-Gaussian laser beam in plasma," IEEE Transactions on Plasma Science, vol. 41 no. 8, pp. 2251-2256, DOI: 10.1109/tps.2013.2268164, 2013.