Propagation Property of an Astigmatic sin – Gaussian Beam in a Strongly Nonlocal Nonlinear Media

Based on the Snyder and Mitchell model, a closed-form propagation expression of astigmatic sin-Gaussian beams through strongly nonlocal nonlinear media (SNNM) is derived. The evolutions of the intensity distributions and the corresponding wave front dislocations are discussed analytically and numerically. It is generally proved that the light field distribution varies periodically with the propagation distance. Furthermore, it is demonstrated that the astigmatism and edge dislocation nested in the initial sin-Gaussian beams greatly influence the pattern configurations and phase singularities during propagation. In particular, it is found that, when the beam parameters are properly selected, a vortex beam with perfect doughnut-shaped profile can be obtained for astigmatic sin-Gaussian beams with two-lobe pattern propagating in SNNM.


Introduction
In the past decades, self-trapped optical beams in nonlocal nonlinear media, especially in strongly nonlocal nonlinear media (SNNM), have been a topic of considerable interest due to their theoretical importance and many experimental observations [1,2].SNNM is a novel kind of media in which the characteristic length of the nonlocal response is much larger than the beam width.Distinct from the conventional local nonlinearity, the nonlocality allows the refractive index of the media at a particular point to be related to the beam intensity at all other points.Note that the nonlocal nonlinearity which can support a variety of nonlocal spatial optical solitons exhibits in many physical systems, and some of them have been observed experimentally [3][4][5][6][7][8][9][10].Moreover, it has been reported that a great number of optical beams can steadily propagate in SNNM under sufficient conditions, including Gaussian beams and higher-order Gaussian beams, four-petal Gaussian beams, Lorentz-Gaussian beams, the beams carrying wave front dislocations such as Hermite-, Hermite-cosh-or Laguerre-Gaussian beams, and so on [11][12][13][14][15][16][17][18][19][20][21][22][23][24][25][26].As well known, pure wave front dislocations in a monochromatic wave are divided into two types: one is the longitudinal screw dislocation which is also known as the optical vortex with spiral phase, and the other is the transverse edge dislocation with π-phase shift located along a line in the transverse plane.During the past years, most of the previous investigations on propagation of beams in SNNM have been focused on those beams whose pattern structure remains invariant although the pattern size may change during propagation.Due to the complexity of evolution of beams in the nonlocal domain, it is interesting to find that the changes of wave front dislocations are possible, which turns out to be very important for the propagation characteristics and will modify the pattern structure for a laser beam in SNNM [24,27,28].
Like Hermite-Gaussian beams, sin-Gaussian beams carrying edge dislocations are the special case of the Hermite-sinusoidal-Gaussian beams whose propagation dynamics through various optical systems had been studied widely [29][30][31][32][33][34][35][36][37].However, their propagation in SNNM remains unexplored.In this work, we aim to study the propagation of sin-Gaussian beams with astigmatism in SNNM.A closed-form propagation equation of astigmatic sin-Gaussian beams in SNNM is derived and illustrated with numerical examples.Our results show some novel variations of wave front dislocations occur for astigmatic sin-Gaussian beams propagating in SNNM.The remainder parts of this paper are organized as follows: Section 2 represents the theoretical formulation for astatic sin-Gaussian beams propagating in SNNM.In Section 3 main characters associated with intensity and phase distributions of a propagating sin-Gaussian beam are analytically investigated and a large number of numerical calculations are performed based on the obtained formulae in Section 2. Finally, in Section 4 the main results obtained in this paper are summarized.

Theoretical Formulation
In the Cartesian coordinate system, a sin-Gaussian beam with astigmatism in the source plane z = 0 takes the form as where β x and β y represent the parameters associated with the sin part, k = 2π/λ and λ is the wavelength, w 0 is the spot size of the beam, and R 0x and R 0y are respectively the curvatures along the x and y directions at z = 0.The astigmatism is represented by the curvature difference between x and y directions and can be easily realized in experiment [38,39].
Because the nonlocal nonlinear Schrodinger equation governing the evolution dynamics of optical beams propagating in nonlocal nonlinear medium remains mathematically complicated, in this paper the medium is assumed to be SNNM and, then, Snyder and Mitchell model is appropriate.In fact, it has been found that the analytical solutions obtained from Snyder and Mitchell model and the numerical simulations based on the nonlocal nonlinear Schrodinger equation are in good agreement in the case of strong nonlocality [12,16,18].The relation between optical beams propagating in SNNM and in quadratic-index media or through optical fractional Fourier transform systems has also been determined [40,41].
The propagation of a sin-Gaussian beam in SNNM is governed by the following equation: where γ is the material constant relating to the response function, and p 0 is the input power in the source plane.The solution of Equation ( 2) can be expressed as [15,21,22] where the ABCD transfer matrix for SNNM has been obtained with the following: with z p = 1/ γ 2 p 0 and z s = z/ πz p being the scaled propagation distance.
Substituting Equation (1) into Equation ( 3) and performing some mathematical manipulations we have with u w = u/w 0 (u = x and y) being scaled transversal coordinates and z R = kw 2 0 /2 = πw 2 0 /λ being the Rayleigh distance.Making use of the Gaussian integral formula One obtains where is an axial symmetry factor on a specified observation plane.Equation ( 8) is the main result obtained in this paper and is also the starting equation to discuss the evolution property of the propagating beam in the next section.

Evolutions of Wave Front Dislocations and Intensity Patterns Occurring in an Astigmatic Sin-Gaussian Beam
In this section, we mainly study the evolution of the intensity pattern as well as the wave front dislocation for a sin-Gaussian beam with astigmatism propagating in SNNM using the analytical formulae derived in previous section.According to Equation ( 8) and the relations of A, B, D, α x and α y with the scaled propagation distance z s , one can generally prove Equation (10) reveals that the evolution of the field distribution with the scaled propagation distance z s is periodic and the period of the field distribution evolution equals to 1 corresponding to the propagation distance πz p = πz R p Gc /p 0 where p Gc = 1/ γ 2 z 2 R is the critical power of Gaussian or higher-order Gaussian beams.Obviously, z p = z R corresponds to p 0 = p Gc while z p < z R corresponds to p 0 > p Gc .
It is known that a sin-Gaussian beam carries an edge dislocation along the straight line β x x + β y y = 0 where β x and β y determines the orientation of the edge dislocation and its intensity pattern is a two-lobe configuration for sufficiently small β x and β y .For R 0x = R 0y (i.e., α x = α y ) corresponding to beams without astigmatism, Equation (8) turns to Regardless of astigmatism, for β x or β y being 0, to be specific, assuming β y = 0, Equation ( 8) becomes Therefore, Equations ( 11) and ( 12) clearly demonstrate that the initial edge dislocation remains unchanged for beams without astigmatism propagating in SNNM.
Figure 1 shows typical evolutions of the intensity distributions of sin-Gaussian beams without astigmatism propagating in SNNM at various scaled propagation distance z s for different values of R 0x (R 0y ) and z p (z R ) with sufficiently small β x and β y .Note that only for z p = z R (i.e., p 0 = p Gc ) and R 0x = R 0y = ∞, sin-Gaussian beam can roughly keep its structure and size unchanged during propagation in SNNM, seeing Row 2 in Figure 1.Otherwise, the pattern size varies with a period of πz p although the intensity pattern continues to be a two-lobe structure during propagation.This is consistent with those of Gaussian and higher-order Gaussian beams in SNNM [1].Moreover, the calculations also indicate that the edge dislocation nested in the input sin-Gaussian beam always remains during propagation.However, for sufficiently large values of β x and β y , the intensity pattern changes periodically and is independent from the input power p 0 , which can be seen from Figure 2. In the following we only consider the case of small β x and β y .On the other hand, for β x = 0, β y = 0 and R 0x = R 0y (i.e., α x =α y ) corresponding to beams with astigmatism, Equation (8) turns into: Since the necessary condition of wave front singularity occurrence is E(x, y, z) = 0, it must simultaneously require Solving Equation ( 14) one obtains These are discrete points at which the intensity is null on the observation plane.Especially, n = 0 corresponds to a solitary zero intensity point at (x w , y w ) = (0, 0).This implies that, due to the simultaneous existence of an edge dislocation and astigmatism in an input beam, the wave front dislocation structure varies and the zero-intensity line of an input edge dislocation can evolve into some discrete null intensity points during propagation.Obviously, an astigmatism (R 0x = R 0y ) as well as the orientation of the initial edge dislocation determined by β x (β y ) are very key factors for the variation of wave front dislocation embedded in a propagating sin-Gaussian beam.
Based on the analytical discussions previously mentioned, a large number of numerical calculations are performed and typical numerical examples are chosen to illustrate the general propagation features for the astigmatic sin-Gaussian beam in SNNM.For the case of R 0x R 0y > 0 or R 0x and R 0y having the same sign as shown in Figure 3, the input astigmatic sin-Gaussian beam is a two-lobe pattern only carrying an edge dislocation.Such propagating beam evolves into one with a dark-hollow pattern and a screw dislocation (vortex) because the lines of zero real-and imaginary-part associated with the propagating field distribution intersect at the origin.It can be seen that the pattern configuration changes with period πz p during propagation, and so do the phase singularities varying between an edge-dislocation and a screw dislocation.When the signs of R 0x and R 0y are different or one of them is ∞, the input sin-Gaussian beam carries a mixed dislocation consisting of a vortex and an edge dislocation at the origin because there is an intersection point at (0,0) of both the straight lines β x x + β y y = 0 and x 2 |R 0x | − y 2 |R0y| = 0 which respectively relate to the zero real-and imaginary-part of the input beam.Figure 4 shows a typical example for β x = β y = 1/3, R 0x = 3z R /4 and R 0y = ∞, revealing that the mixed dislocations in the source plane turns into a vortex with period πz p during propagation.The calculations further confirm that the size of the patterns also depends on the values of the input power p 0 .For specified choices of R 0x = −R 0y and smaller values of β x = ± β y , interestingly, it can be observed that an astigmatic sin-Gaussian beam evolves into a doughnut-shaped pattern with a uniform bright enclosure.Figures 5 and 6 show evolutions of the intensity distributions of astigmatic sin-Gaussian beams in SNNM for different values of z p with R 0x = −R 0y = z R , β x = β y = 1/3.Obviously, the propagating beam forms a perfect doughnut-shaped pattern across the transverse plane at z s = 1/2.Note that the size of the preferred doughnut-shaped pattern at z s = 1/2 sensitively depends on the values of the parameter z p , or to be more specific, the size of the doughnut-shaped pattern decreases with an increase of the input power p 0 .Moreover, we also investigate the wave front dislocation evolutions of the propagation field, demonstrating that the canonical vortex is invariant for z p = z R but evolves into noncanonical for other values of z p , which can be seen from these figures.

Figure 2 .
Figure 2. Evolution of intensity distributions of sin-Gaussian beams with the scaled propagation distance z s for β x = 4, β y = 3.The other parameters are R 0x = R 0y = ∞ and z p = z R .