Abstract
Analytically discussed is the dynamical properties of the spiraling elliptic beams in nonlocal nonlinear media. This class of spiraling elliptic beams carry the orbital angular momentum (OAM), and can rotate on the cross section perpendicular to the propagation direction during the propagations. The optical intensity, the beam width, and specially the angular velocity are both analytically and numerically discussed in details. We shown that both the deviations from the critical power and the deviations from the critical OAM can make the spiraling elliptic beams breathe. The decrease (increase) of the OAM or the increase (decrease) of the power can both make the spiraling elliptic breathers contract (diffract), however, there still exist differences between them. The rotating speed can be changed by the input optical power or the input OAM, which may have potential applications in the controlling of the optical beams.
© 2015 Optical Society of America
1. Introduction
The technologies associated with orbital angular momentum (OAM), including spatial light modulators and hologram design, have found their own applications ranging from optical tweezers to microscopy [1]. A kind of the spiraling elliptic solitons carrying the OAM was introduced in saturable nonlinear media [2], and was systematically discussed in media of the nonlocal nonlinear media without anisotropy [3], where it was claimed that the OAM can result in the effective anisotropic diffraction for the spiraling elliptic beams. In the nonlinear media with linear anisotropy, the OAM is not conserved. Depending on the linear anisotropy of the media, two kinds of evolution behaviors for the dynamic breathers, rotations and molecule-like librations were predicated analytically and confirmed in numerical simulations [4]. Very recently, the polarized vector spiraling elliptic solitons were reported to stably exist in nonlocal nonlinear media [5]. It was shown that the spiraling elliptic solitons can be formed and can propagate stably in the waveguide induced by orthogonally polarized Gaussian beams.
It has been claimed that [3] only when the input optical power is equal to the critical power and the input OAM is equal to the critical OAM the spiraling elliptic beams can evolve as the spiraling elliptic solitons. In fact, external perturbations will easily make the input power and the input OAM deviate from their critical values, in this general case the evolution of the spiraling elliptic beams is still an open question. For the special case of the strongly nonlocal nonlinear (SNN) media in which the characteristic length of the material response function is much larger than the beam width, the propagation equation can be linearized to the well-known Snyder-Mitchell model (SMM) [6]. The propagation of optical beams in SNN media can be simply regarded as a self-induced fractional Fourier transform [7], then there exists a one-to-one correspondence between free propagation and propagation in SNN media [8]. In this paper we use the established one-to-one correspondence to analytically investigate the dynamic properties of the spiraling elliptic beams. We find that both the deviations from the critical power and the deviations from the critical OAM can make the spiraling elliptic beams breathe. The decrease (increase) of the OAM or the increase (decrease) of the power can both make the spiraling elliptic breathers contract (diffract), however, there still exist differences between them. In addition, the rotating speed can be changed by the input optical power and the input OAM, it may have potential applications in the controlling of the optical beams.
2. Theoretical model and analytical solution
The propagation of the (1+2)-dimensional optical beam in nonlocal nonlinear media is described by the nonlocal nonlinear Schrödinger equation (NNLSE) [6, 9, 10]
where ψ(r, z) is the complex amplitude envelope, Δn = ∬R(r − r′)|ψ(r′)|2d2r′ is the light-induced nonlinear refractive index, z is the longitudinal coordinate, r is the transverse coordinate vector with r = xex + yey, is the transverse Laplacian operator, and R is normalized symmetrical real spatial response of the media such that Δn = ∬R(r′)d2r′ = 1. In the case of strongly nonlocal nonlinear (SNN) media we need only keep the first two terms of the expansion of Δn and the NNLSE is simplified to the Snyder-Mitchell model [6, 11] where , and P0 = ∬|ψ(r′)|2d2r′ is the input optical power.With an established one-to-one correspondence between free and SNN propagation [7, 8, 12], it would be easy to deal with the propagation of an arbitrary field in SNN media by using the well-developed theory of free propagation. The expression of paraxial optical beams at any propagation distance z in SNN media can be obtained by the following integral formula [7, 12]:
where ψ(x0, y0, 0) is the input optical beam. The input optical beam considered here has the form of which has an elliptic shape with b and c being the major semiaxis and the minor semiaxis, and has a cross term Θx0y0 in the phase. We have shown that the cross term in the phase will result in the rotation of the optical beam [3]. The optical beam expressed as Eq. (4) owns the OAM It is known that the NNLSE, i.e. Eq. (1) conserves the OAM [13]. Substituting Eq. (4) into Eq. (3) yields the expression of the optical beam at any propagation distance z: where withThe analytic solution, i.e. Eq. (6), the key result of the paper, is an exact one to the Snyder-Mitchell model, i.e. Eq. (2). The spiraling elliptic beam expressed as Eq. (6) carries the OAM and rotates on the cross section perpendicular to the propagation direction z during the propagations.
3. Discussion of solutions
From Eq. (6), it can be easily shown that the shape of the spiraling elliptic beam is an inclined ellipse, whose major axis and minor axis are not parallel to x-axis and y-axis. The elliptic optical beam, i.e. Eq. (6), can be transformed to its standard elliptic form by a coordinate rotation through angle ϑ, which satisfies the following equation
where One of the semi-axes of the standard elliptic optical spot is the other is The angular velocity of the optical beam in the xyz-coordinate frame can be obtained from Eq. (18)3.1. Spiraling elliptic solitons
Generally speaking, the semi-axes wb and wc are the functions of the propagation distance z. As the special case, the spiraling elliptic solitons can be obtained by letting and . For the spiraling elliptic solitons, from Eqs. (20) and (21) we have
The critical OAM can be obtained by combining Eq. (5) and Eq. (24) as Substitution of Eqs. (23) and (25) into Eq. (22) yields which shows that the spiraling elliptic solitons make constant-angular rotations. Then the period of rotation can be obtained as .We suppose the material response to be the Gaussian function [14, 15]
where wm is the width of the response function. The ratio of the width of the response function to the scale in the transverse dimension occupied by the optical beam determines the degree of nonlocality, which reads . The larger is δ, the stronger is the degree of nonlocality. For the response function of Eq. (27), we obtain that . Figs. 1 present the propagation dynamics of the spiraling elliptic solitons. The excellent agreement between the exact analytical spiraling elliptic soliton solutions of the Snyder-Mitchell model, i.e. Eq. (2), and the numerical simulation of the NNLSE for the case of strong nonlocality δ = 10 is shown in Figs. 1. When the degree of nonlocality becomes weaker (δ = 5), the approximations of the analytical results of the Snyder-Mitchell model to the exact ones of the NNLSE are a little bit worse. When the degree of nonlocality is weak enough (δ = 1), the dynamics of the spiraling elliptic solitons will be unstable.3.2. Spiraling elliptic breathers
In general, for the case that either the input optical power or the input OAM is not equal to the critical power and the critical OAM, all the spiraling elliptic beams will evolve as the spiraling elliptic breathers, following Eq. (6). In the following, we will consider the breathing behavior of the spiraling elliptic breathers, and discuss the evolutions of the optical intensity, the beam width, and the angular velocity in details.
3.2.1. Breathing behavior resulting from input optical power
It is well known that [6] for P0 > Pc, self-focusing initially overcomes diffraction and the beam initially contracts shown in Figs. 2, whereas for P0 < Pc, the reverse happens and the beam initially expands shown in Figs. 3. The excellent agreement between the exact analytical spiraling elliptic breathers of the Snyder-Mitchell model and the numerical simulation of the NNLSE for the case of strong nonlocality δ = 10 are also shown in Figs. 2 and Figs. 3, where we have assumed that the input optical power P0 = 1.5Pc and P0 = 0.5Pc, respectively, but meanwhile kept the input OAM equal to the critical OAM.
We start our considerations with the evolution of the on-axis optical intensity for the spiraling elliptic breathers. Note, the intensities of the spiraling elliptic beams are periodic functions of the propagation distance z, the period can be obtained from Eq. (7) as , which can be rewritten as
for the Gaussian response function of Eq. (27). It is shown from Eq. (28) that the evolutional periods of the spiraling elliptic breathers decrease as the input optical powers increase. Taking parameters as P0 = 1.5Pc and P0 = 0.5Pc, the evolutional periods can be calculated to be 4.1 and 7.1, respectively. The analytical results are confirmed by the numerical simulation shown in Fig. 4. Under the assumption that the input OAM equals to the critical OAM, the on-axis optical intensity will increase initially and then decrease if P0 > Pc, and will decrease initially and then increase if P0 < Pc. Both the maximum on-axis optical intensity for the case that P0 > Pc and the minimum on-axis optical intensity for the case that P0 < Pc can be obtained from Eq. (7) as , the ratio of which to the critical optical intensity is Im/Ic = P0/Pc, with the critical optical intensity being . Therefore, for the spiraling elliptic breathers the ratio of the on-axis optical intensity to the input optical intensity vibrates between 1 and P0/Pc, which also agrees well with the numerical simulation shown in Fig. 4.And then, we investigate the evolution of the major semiaxis wb and the minor semiaxis wc of the spiraling elliptic breathers. By using Eqs. (20) and (21), after somewhat complicated calculations we obtain that the ratio of the major semiaxis wb to b and the ratio of the minor semiaxis wc to c vibrate between 1 and . Fig. 5 shows the evolution of the semi-axes wb and wc of the spiraling elliptic breathers, and the analytical solution agrees well with the numerical simulation. When we take P0 = 1.5Pc, the major semiaxis wb vibrates between 1.6 (i.e., ) and 2, and the minor semiaxis wc vibrates between 0.8 (i.e., ) and 1. And when we take P0 = 0.5Pc, the major semiaxis wb vibrates between 2 and 2.8 (i.e., ), and the minor semiaxis wc vibrates between 1 and 1.4 (i.e., ).
In the end, we discuss the angular velocity of the spiraling elliptic breathers. For the case that the input optical power is not equal to the critical power but the input OAM is equal to the critical OAM, the angular velocity can be obtained by using Eqs. (22) and (24)
where the critical optical power Pc and the critical angular velocity ωc are expressed as Eqs. (23) and (26), respectively. It can be found that the spiraling elliptic breathers perform nonuniform rotations, the angular velocity is not a constant but is a periodic function of the propagation distance z. From Eq. (29) with the aid of Eqs. (23) and (26), after some simple calculations, the angular velocity can be found vibrating between ωc and , as shown in Fig. 6. It is interesting to note that the angular velocity of the spiraling elliptic breathers depends on the optical power P0. Consequently, the angular velocity of the spiraling elliptic beams can be controlled by the input optical power.3.2.2. Breathing behavior resulting from input OAM
In this section, we consider the evolutions of the spiraling elliptic breathers for the case that the input OAM is not equal to the critical OAM but keeping P0 = Pc. It is shown that [2, 3] the OAM can introduce an effective diffraction into the spiraling elliptic beams. Therefore, when we increase the input OAM, the spiraling elliptic beams will diffract in the whole, as shown in Fig. 7. Conversely, when we decrease the input OAM, the spiraling elliptic beams will contract in the whole, as shown in Fig. 8. Comparing Figs. 2 and Figs. 8 (or comparing Figs. 3 and Figs. 7), it seems that both the decrease (increase) of the OAM and the increase (decrease) of the power have the similar effects on the evolutions of the spiraling elliptic breathers in the sense that they both make the spiraling elliptic breathers contract (diffract). But we will show that there still exist differences between the effects of the decrease (increase) of the OAM and the increase (decrease) of the power on the evolutions of the spiraling elliptic breathers.
Firstly, we investigate the effect of the OAM on the evolution of the on-axis optical intensity for the spiraling elliptic breathers. To this end, we assume that the input OAM departs from the critical OAM but the input optical power is equal to the critical power, in this case we can obtain from Eq. (7) that the ratio of the on-axis optical intensity I(z) to the input optical intensity I(0) has such three extrema as
with κ denoting the extrema of I(z)/I(0). Fig. 9 gives the evolution of the on-axis optical intensity when σ ≠ σc and P0 = Pc. From Eqs. (30), (31) and (32), the three extrema of the intensity can be obtained as 1, 1.05 and 0.69 when σ = 1.5σc, and obtained as 1, 0.97 and 1.37 when σ = 0.5σc, which are confirmed by the numerical simulation in Fig. 9. Comparing Fig. 9 with Fig. 4, it is observed that although both the decrease (increase) of the OAM and the increase (decrease) of the power can make the spiraling elliptic breathers contract (diffract), the effect of the OAM is obviously more complicated. The decrease (increase) of the OAM makes spiraling elliptic beams contract (diffract) on the whole, but in some small propagation ranges the spiraling elliptic beams diffract (contract).Subsequently, we consider the effects of the OAM on the evolutions of the semi-axes for spiraling elliptic beams. When we increases the input OAM, the OAM will enhance the diffraction, and the semi-axes will be increased to their maxima
respectively, as shown in Fig. 10. Differently form the case in Fig. 5 that it can be observed that the minor semiaxis initially contracts and the major semiaxis initially expands in Fig. 10. On the contrary, when the input OAM decreases, the effective diffraction of the spiraling elliptic beams will be weakened, then the semi-axes will be decreased to their minimums wbm and wcm, respectively, but the major semiaxis will initially contract and the minor semiaxis will initially expand, as shown in Fig. 10.In addition, the angular velocity of the spiraling elliptic breathers can be obtained from Eq. (22) when σ ≠ σc and P0 = Pc, as shown in Fig. 11. It is shown from Fig. 11 that the angular velocity of the spiraling elliptic breathers also depends on the input OAM. It was claimed that [2] when the elliptic beams are focused by a cylindrical lens, and the major axis of the elliptically shaped beam makes an angle with the cylindrical lens, the elliptic beams will carry the OAM. Therefore, the angular velocity of the spiraling elliptic beams can also be changed by the angle made between the major axis of the elliptically shaped beam makes and the cylindrical lens.
Before concluding, we give a short discussion of the applicability of the Snyder-Mitchell model of Eq. (2), and perform a comparison between the Snyder-Mitchell model and the actual physical systems. The Snyder-Mitchell model [6] is a phenomenological model, and is also the simplified model of the NNLSE for the case of the strong nonlocality (i.e. the the degree of nonlocality δ ≫ 1) and the condition that the response function R(x, y) is symmetrical and regular (or at least twice differentiable) at x = y = 0, such as the extremely instructive Gaussian shaped response function of Eq. (27). The derivation from the NNLSE to the Snyder-Mitchell model is shown for the strong nonlocality in [11]. On the other hand, the response functions of the actual physical systems, such as the SNN media nematic liquid crystal (NLC) and the lead glass, are always singular [11]. The singularity in the response functions results in the actual physical systems such as the NLC and the lead glass can not be simplified to the Snyder-Mitchell model [11]. However, the Snyder-Mitchell model keeps the main features of the SNN media, and the physical properties do not depend strongly on the shapes of the response functions for the SNN media [16, 17]. For example, the theoretical predictions by the Snyder-Mitchell model [6], such as the accessible solitons and the attraction of spatial solitons, have been observed in experiments in the NLC [18, 19, 20] and the lead glass [21]. It is worth mentioning that [7] the fractional Fourier transform existing the SNN media was predicted by the Snyder-Mitchell model, and was also observed in the lead glass. Therefore, the dynamics of elliptic beams in highly nonlocal media discussed in the paper by the Snyder-Mitchell model pave the way to the observation in the NLC and the lead glass in experiments.
4. Conclusion
We analytically discussed the dynamical properties of the spiraling elliptic beams in nonlocal nonlinear media. This class of spiraling elliptic beams carry the orbital angular momentum (OAM), and can rotate on the cross section perpendicular to the propagation direction during the propagations. We obtained the analytical solutions for the evolution period,the beam width, and specially the angular velocity of the spiraling elliptic beams in strongly nonlocal nonlinear media. Although both the deviations from the critical power and the deviations from the critical OAM can make the spiraling elliptic beams breathe, there still exist differences between them. The breathing behaviors resulting from the OAM are more complicated. The rotating speed can be changed by the input optical power and the input OAM, which may have potential applications in the controlling of the optical beams.
Acknowledgments
This research was supported by the National Natural Science Foundation of China, Grant Nos. 11274125 and 11474109.
References and links
1. S. Franke-Arnold, L. Allen, and M. Padgett, “Advances in optical angular momentum,” Laser Photon. Rev. 2(4), 299–313 (2008). [CrossRef]
2. A. S. Desyatnikov, D. Buccoliero, M. R. Dennis, and Y. S. Kivshar, “Suppression of Collapse for Spiraling Elliptic Solitons,” Phys. Rev. Lett. 104, 053902 (2010). [CrossRef] [PubMed]
3. G. Liang and Q. Guo, “Spiraling elliptic solitons in nonlocal nonlinear media without anisotropy,” Phys. Rev. A 88, 043825 (2013). [CrossRef]
4. G. Liang, Q. Guo, Q. Shou, and Z.M. Ren, “Dynamics of elliptic breathers in saturable nonlinear media with linear anisotropy,” J. Opt. 16, 085205 (2014). [CrossRef]
5. G. Liang and H. G. Li, “Polarized vector spiraling elliptic solitons in nonlocal nonlinear media,” Opt.Commun. 352, 39 (2015).
6. A. W. Snyder and D. J. Mitchell, “Accessible Solitons,” Science 276, 1538 (1997). [CrossRef]
7. D.Q. Lu, W. Hu, Y. J. Zheng, Y. B. Liang, L. G. Cao, S. Lan, and Q. Guo, Phys. Rev. A “Self-induced fractional Fourier transform and revivable higher-order spatial solitons in strongly nonlocal nonlinear media,” Phys. Rev. A 78, 043815 (2008). [CrossRef]
8. D. Q. Lu, W. Hu, and Q. Guo, “The relation between optical beam propagation in free space and in strongly nonlocal nonlinear media,” Eur. Phys. Lett. 86, 44004 (2009). [CrossRef]
9. D.J. Mitchell and A.W. Snyder, “Soliton dynamics in a nonlocal medium,” J. Opt. Soc. Am. B 16, 236 (1999). [CrossRef]
10. W. Krolikowski, O. Bang, J. J. Rasmussen, and J. Wyller, “Modulational instability in nonlocal nonlinear Kerr media,” Phys. Rev. E 64, 016612 (2001). [CrossRef]
11. Q. Guo, D. Q. Lu, and D. M. Deng, “Nonlocal spatial optical solitons,” in Advances in Nonlinear Optics, X. Chen, Q. Guo, W. She, H. Zeng, and G. Zhang, eds. (De Gruyter, 2015), pp. 277–306.
12. Z. P. Dai, Z. J. Yang, S. M. Zhang, and Z. G. Pang, “Propagation of anomalous vortex beams in strongly nonlocal nonlinear media,” Opt.Commun. 350, 19 (2015).
13. A. I. Yakimenko, V. M. Lashkin, and O. O. Prikhodko, “Dynamics of two-dimensional coherent structures in nonlocal nonlinear media,” Phys. Rev. E 73, 066605 (2006). [CrossRef]
14. W. Krolikowski, O. Bang, N. I. Nikolov, D. Neshev, J. Wyller, J. J. Rasmussen, and D. Edmundson, “Modulational instablity, solitons and beam propagation in spatially nonlocal nonlinear media,” J. Opt. B 6, S288 (2004). [CrossRef]
15. Q. Guo, B. Luo, F. Yi, S. Chi, and Y. Xie, “Large phase shift of nonlocal optical spatial solitons,” Phys. Rev. E 69, 016602 (2004). [CrossRef]
16. J. Wyller, W. Krolikowski, O. Bang, and J.J. Rasmussen, “Generic features of modulational instability in nonlocal Kerr media,” Phys. Rev. E 66, 066615 (2002). [CrossRef]
17. O. Bang, W. Krolikowski, J. Wyller, and J.J. Rasmussen, “Collapse arrest and soliton stabilization in nonlocal nonlinear media,” Phys. Rev. E 66, 046619 (2002). [CrossRef]
18. C. Conti, M. Peccianti, and G. Assanto, “Route to nonlocality and observation of accessible solitons,” Phys. Rev. Lett. 91, 073901 (2003). [CrossRef] [PubMed]
19. C. Conti, M. Peccianti, and G. Assanto, “Observation of optical spatial solitons in highly nonlocal medium,” Phys. Rev. Lett. 91, 073901 (2003). [CrossRef]
20. M. Peccianti, K.A. Brzdkiewicz, and G. Assanto, “Nonlocal spatial soliton interactions in nematic liquid crystals,” Opt. Express 27(6), 1460-1462 (2002). [CrossRef]
21. C. Rotschild, O. Cohen, O. Manela, M. Segev, and T. Carmon, “Solitons in Nonlinear Media with an Infinite Range of Nonlocality: First Observation of Coherent Elliptic Solitons and of Vortex-Ring Solitons,” Phys. Rev. Lett. 95, 213904 (2005). [CrossRef] [PubMed]