Abstract
We investigate both analytically and numerically the propagation dynamic of on-axis and off-axis cosine-Gaussian (CG) beams in a linear medium with quadratic external potential. CG beam propagation evolves periodically with a period depended on the potential depth (α) and whether the beam shape is symmetrical with respect to optical axis. In each period, the CG beam first splits into two sub-beams with different accelerated direction; they then reverse the accelerated direction owing to the quadratic external potential and finally merge again to reproduce its initial shape, and the whole process repeats periodically. The intensity oscillation period of the off-axis CG beam is double times than that of the on-axis one. At the special position, the beam (or spectral) shape is strongly related to the initial spectral (beam) shape. The corresponding scaled relationship is that the spatial intensity Ix (or spatial frequency axis k) is α times the spectral intensity Ik (or space axis x). The interaction of two spatially separated CG beams still exhibit periodic evolution with complex structure in the regime of focal point. The propagation dynamics of two-dimensional CG beams are also presented. When the propagation distance is exactly an integer multiple of half period, there are four focal points in the diagonal position.
© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement
1. Introduction
The spreading of wave packets originated from dispersive or/and diffraction is a universal phenomenon, which is a limiting factor in many practical applications [1]. Therefore, it becomes important to develop methods to overcome or reduce such phenomenon. Non-diffracting waves, also known as localized waves are able to resist diffraction and have attracted significant and increasing research interest [2–5]. For instance, the different kinds of solitons are of the well-known example of non-diffracting waves in nonlinear regimes [6–8]. On the other hand, in the linear approximation, non-diffracting beams (or non-dispersive pulses) have also theoretically proposed and experimentally demonstrated [9–15]. Non-diffracting wave packets have also been expanded from optics to surface plasmonic, acoustics, and electronics [16].
Durnin et al. first found non-diffracting Bessel beams in 1987 [17,18], which stimulated much attention for finding another kinds of non-diffracting beams. Cosine beams (CBs), a type of non-diffracting beam, are exact solution with the form of a cosine function to the Helmholtz equation, whose intensity shape is independent of the propagation distance [19,20]. However, it is impossible to realize CBs physically because its total energy is infinite at any transverse plane. To overcome such problem, then finite-energy CBs were introduced by imposing a truncated factor on the CBs [21]. Cosine-Gaussian (CG) beam is a typical one of the truncated version [22], which can be considered as special cases of the Hermite-sinusoidal-Gaussian beams as the order of Hermite-polynomial function is zero [23,24]. Compared to the ideal CBs, CG beam still exhibits unique ability of diffraction-free and self-healing in a finite propagation distance and beyond that region they are inevitably dispersive [25]. Because of their self-healing and non-diffracting natures, there are much increasing attention towards the propagation dynamics of the CG beam in different system, including free space [26], turbulent atmosphere [27–30], uniaxial crystals [31], correlated beams [32–37], plasmonic [38–41] and even water waves in hydrodynamics [42,43], and so on.
As an effective method to control beams propagation, the external potentials, such as linear and parabolic potentials, have been introduced. Linear and dynamic linear potentials have been used to manipulate the beam trajectory during propagation [44–47]. As for parabolic potential, propagation dynamics of different kinds of beams, such as symmetric and asymmetric beams, as well as beams carried orbital angular momentum have been investigated [48–60]. There are some novel and interesting features during propagation process, including periodic inversion and focusing [48–54], and self-induced Fourier transformation [55]. Besides, the manipulation of the propagation dynamics of beams or pulses via its spectral phase modulation have also been disclosed [61,62].
The CG beams can be treated as a superposition of two Gaussian beams with different accelerating direction. The impact of linear potential on its propagation dynamic has been reported [44]. In this work, we investigated the propagation dynamics of on-axis and off-axis CG beams in linear media with external quadratic potential. We obtained the analytical expression of CG beams propagation in the presence of external quadratic potential by Fourier analysis, which is quite different from previous researches related the CG beams based on ABCD matrix. The analytical results were confirmed by directly numerical calculation of the Schrödinger equation with external quadratic potential. The propagation dynamics of dual CG beams as well as two-dimensional CG beams propagations are also discussed.
2. Theoretical model
In the paraxial approximation, beam propagation in a linear media with quadratic external potential can be described by the following dimensionless Schrödinger equation [48–54]
The analytical expression of general solution to Eq. (1) can be written as [50]
whereNotably seen in Eq. (4), the integral can be seen as the Fourier transform of $\phi ({\xi ,0} )\exp ({ib{\xi^2}} )$ with the spatial frequency K, that is if we chooses an exact input beam $\phi ({\xi ,0} )$, the analytical evolution solution will be obtained by computing the Fourier transform of $\phi ({\xi ,0} )\exp ({ib{\xi^2}} )$. In simple terms, we can describe the propagation of beams in a parabolic potential in terms of the periodic conversion from the initial beam to the Fourier transform with quadratic chirp: which process can be named self-Fourier transform [54].
3. Analytical and numerical results
3.1 Single CG beam
The initial beam was chosen as a general CG beam,
By using analytical expression of Eq. (7), we plot the propagation dynamics of CG beam as a function of propagation distance. Figure 1 shows the spatial evolution of the on-axis and off-axis CG beams with $\omega = 1$ and $\sigma = 0.1$ in a linear media with different potential depth $\alpha = 0.2$ and $\alpha = 0.4$. It can be clearly seen from Fig. 1 that the CG beam exhibits periodic evolution pattern with initial shape recurring at special propagation distance. The CG beam first splits into two branches and then they are focused into spots individually. With a further increase in propagation distance, such two sub-beams are combined again. The process of separation-focus-combination is reproduced periodically. The period of beam pattern and the distance between two focal spots are independent of the transverse displacement (${x_0}$). While they are very sensitive to the depth of potential ($\alpha $), which decreases with an increasing $\alpha $. This is apparently disclosed by comparing left column of Fig. 1 ($\alpha = 0.2$) with right one ($\alpha = 0.4$). By changing the values of transverse displacement (${x_0}$) from zero to nonzero, the CG beam is switched from on-axis state to off-axis states. As a result, the beam shape changes from symmetric to asymmetric with respect to $x = 0$. So, it is easy to conclude that the intensity period of symmetric beam evolution is double times that of asymmetric one in linear media with quadratic potential.
In order to deep the understanding of periodic evolution presented in Fig. 1 and quantify the impact of $\alpha $ and ${x_0}$ on the cosine beam propagation, we further analyze Eq. (7). Within mind, CG beam described by Eq. (5) can be decomposed into two Gaussian beams with opposite linear chirp. Their trajectories can be deduced from Eq. (6)
To confirm analytical results shown in Fig. 1, we directly solve the Eq. (1) with input of CG beams by using split-step Fourier method. Figure 2 displays numerical results of spatial and spectral evolution of CG beams versus propagation distance. All parameters used in numerical simulations are same as Fig. 1. By comparing Figs. 2(a) and 2(c) to Figs. 1(c) and 1(a), they exhibit perfect agreement. In order to quantitatively compare numerical and analytical results, the spatial intensities at two special propagation distances obtained from numerical simulation and analytical expression are shown in Figs. 3(a) and 3(b). Here the on-axis CG beam was selected for comparison. When the propagation distance is the integer multiple of the intensity period ${T _s} = {{{T_{as}}} / 2}$, that is $z = m{T _s}$ with $m = 1,2,3 \cdots$, the corresponding beam's shape is same as the initial one $\phi ({x,z = m{T_s}} )= \phi ({x,0} )$; when the propagation distance is the odd integer multiple of the half intensity period, that is ${z_h} = {{({2m + 1} ){T _s}} / 2}$ with $m = 1,2,3 \cdots$, the beam's shape can be expressed as
A comparison of the spatial and spectral evolution plotted in Fig. 2 shows that the spectral pattern resembles the spatial one, there is only a phase shift of one quarter period (${\pi / {2\alpha }}$) and different scales. To understand such interesting feature of Fig. 2, by use of Fourier transformer, Eq. (1) can be written as in moment space
Next, we carry out a series of numerical simulations to confirm the period ${T _s} = {{2\pi } / \alpha }$ and distance between two focus spots $\Delta {\kern 1pt} x = {{2\omega } / \alpha }$ obtained from analytical expression. Both of them depend on the potential depth. The period scales inversely with the potential depth. The separation between two focal spots varies linearly with modulation frequency and inversely with the potential depth. Figure 4 shows the dependence of ${T _s}$ and $\varDelta x$ with the potential depth or the modulation frequency. It can be clearly seen from Fig. 4 that the agreement between the analytical results and numerical simulations holds very well. Once again, the analytical predictions have been confirmed by numerical experiments.
3.2 Dual CG beams
Next, we investigate the two CG beams interaction governed by Eq. (1). We first consider the case of two beams with in-phase. So, the incident beam can be written as
As seen in Eq. (14), the propagation of dual CG beams is made of four Gaussian beams evolution, which accelerate along the four branches of trajectories. Correspondingly, their trajectories can be expressed as
Let focus our attention on the focal regime. Compared to the Gaussian shape of focal spots for the case of single input CG beam, dual CG beams will form two multi-peaked structures at the propagation distance ${z_F} = {{({2m + 1} ){T_s}} / 2}$. While the locations of other two deviate from ${z_F}$. One location is $m{T_s} + {z_1}$, the other one is $m{T_s} + {z_2}$. From Eq. (16), ${z_{1,2}}$ can be obtained
With further increase propagate distance to ${z_2}$, the bean intensity profile is same as that at ${z_1}$, the expression changes only in phase, which can be written as
Figures 6(b)–6(c) display the analytical and numerical results of beam's intensity distribution at three propagation distances ${z_1}$, ${z_2}$ and ${{{T_s}} / 2}$. For comparison, the case of incident beam profile was shown in Fig. 6(a). The numerical results show an excellent agreement with analytical predictions. As the beams propagate along the media, sub-beams are created after splitting process, and then come close to each other. At the propagation distance ${z_1}$ shown in Fig. 6(b), the shape is constructed by an oscillatory structure covering the center range and other two smooth peaks located on its wings. The peaks on both sides are Gaussian profiles, but the innermost two sub-beams merge together to form a multiple peaks structure, resulting from its coherent superposition. Because their relative phase difference is equal to that at input, it still a CG beam with the same modulation frequency and different width compared with that of input beam. This feature can be directly observed by comparing the numbers of multiple peaks. As the propagation distance is increased to halve period $z = {{{T_s}} / 2}$, the inner and outer beams with farthest separation fully overlap each other, leading to form two new CG beams with small modulation frequency compared to that of incident Gaussian beam. Because the accumulated phase shifts of inner beams are different from that of outer beams. As a result, the modulation frequency changes. But the new CG beams are symmetric with respect to $x = 0$. At the propagation distance ${z_2}$, the beam profile is mirror symmetric to that at ${z_1}$. Although, the both shapes are the same, the left and right parts are exchanged.
3.3 Two-dimensional CG beam
In the former sections, we have showed the harmonic oscillation and self-Fourier transform of the one-dimensional CG beam propagating in linear media with parabolic potential. The question is, will the same phenomena still appear when the beam extends to two dimensions? Next, we analyze the two-dimensional CG beam propagation, which was governed by
As can be seen from Fig. 7, we can clearly see recurrence of incident beam at the integer multiple of the period ${T_s}$, and when the propagation distance is an odd integer multiple of the half period ${z_{odd}} = {{({2m + 1} ){T_s}} / 2}$, the beam will be divided into four spots, whose position can be computed from the analytical solution expression:
According to the properties of Gaussian beam, one can see that the beam divided into four spots, as shown in Fig. 7(c), whose intensity profiles are of Gaussian distributions. Their center coordinates are $({ - {\omega / \alpha }, - {\omega / \alpha }} )$, $({ - {\omega / \alpha },{\omega / \alpha }} )$, $({{\omega / \alpha }, - {\omega / \alpha }} )$ and $({{\omega / \alpha },{\omega / \alpha }} )$, which can be directly seen from the Eq. (26). These four focal spots located at two diagonals of square, not at center perpendicular lines. The reason is that x and y dimensions are indeed found to be correlated during beam propagation. As seen in Fig. 7, the modulation frequencies in two dimensions are the same, the focal spots are going to appear on the diagonal. It is also easy to understand by use of the relation between sum and difference of angles in trigonometry functions as $2\cos ({{\omega_1}x} )\cos ({{\omega_2}y} )= \cos ({{\omega_1}x + {\omega_2}y} )+ \cos ({{\omega_1}x - {\omega_2}y} )$. Therefore, the focal spots are satisfied relationship ${\omega _2}x \pm {\omega _1}y = 0$, and the corresponding coordinates are $({{{ \pm {\omega_1}} / \alpha },\textrm{ }{{ \pm {\omega_2}} / \alpha }} )$. Moreover, the positions of focal spots can be controlled by varying the modulation frequency of CG beam. This relationship also indicates the focal spots do not appear at center perpendicular lines in two dimensional cases. It is the prerogative of a one-dimensional CG beam.
4. Conclusion
In summary, we have analytically and numerically investigated the propagation dynamics of on-axis and off-axis CG beams in a medium with an external quadratic potential. The CG beams are the interference pattern originated from two Gaussian beams with different accelerating directions. The analytical expression of the CG beams propagation in a quadratic potential showed their periodic evolution described by the Fourier transform of the initial beam with a quadratic chirp. Due to the presence of external quadratic potential, the CG beams experience a periodic oscillation process, which can be characterized as splitting, focus and coalescence. The oscillated period and the distance between two focal spots are inversely proportional to the potential depth ($\alpha $). They remain the same for both on-axis and off-axis CG beams. While the trajectory of sub beams made a transition from sine curves to cosine-like curves as the incident on-axis CG beams were substituted by the off-axis one. Compared to the on-axis CG beams, the off-axis CG beams are asymmetric with respect to $x = 0$, Consequently, The period of intensity evolution of on-axis CG beam is half of that of on-axis one.
The dual CG beams, can be treated as a combination of two CG beams with transverse displacements, will perform complex interference and form four multi-peaks structure near the phase transition point in a period. Two-dimensional CG beams have the similar phenomenon to the one-dimensional case. But at the focal point, there are four focal spots located at four corner, their locations and interval between focal points are determined by transverse modulation frequencies and the depth of quadratic potential. The present investigation will be helpful to deepen understanding of CG beams with quadratic external potential, and also point out the potential applications of beam shaping technology.
Funding
National Natural Science Foundation of China (61975130).
Disclosures
The authors declare no conflicts of interest.
References
1. M. Born, E. Wolf, and A. B. Bhatia, Principles of Optics: Electromagnetic Theory of Propagation, Interference and Diffraction of Light, 7th ed. (Cambridge University, 1999).
2. M. V. Berry and N. L. Balazs, “Nonspreading wave packets,” Am. J. Phys. 47(3), 264–267 (1979). [CrossRef]
3. M. Mazilu, D. J. Stevenson, F. Gunn-Moore, and K. Dholakia, “Light beats the spread: non-diffracting beams,” Laser Photon. Rev. 4(4), 529–547 (2010). [CrossRef]
4. D. Mihalache, “Linear and nonlinear light bullets: recent theoretical and experimental studies,” Rom. J. Phys. 57, 88820J (2012). [CrossRef]
5. H. E. Hernández-Figueroa, E. Recami, and M. Zamboni-Rached, Non-Diffracting Waves, (Wiley, 2013).
6. Y. S. Kivshar and G. P. Agrawal, Optical Solitons: From Fibers to Photonic Crystals, (Academic, 2003).
7. F. Lederer, G. I. Stegeman, D. N. Christodoulides, G. Assanto, M. Segev, and Y. Silberberge, “Discrete solitons in optics,” Phys. Rep. 463(1-3), 1–126 (2008). [CrossRef]
8. Z. Chen, M. Segev, and D. N. Christodoulides, “Optical spatial solitons: historical overview and recent advances,” Rep. Prog. Phys. 75(8), 086401 (2012). [CrossRef]
9. A. Chong, W. Renninger, D. N. Christodoulides, and F. W. Wise, “Airy-Bessel wave packets as versatile linear light bullets,” Nat. Photonics 4(2), 103–106 (2010). [CrossRef]
10. D. Abdollahpour, S. Suntsov, D. G. Papazoglou, and S. Tzortzakis, “Spatiotemporal airy light bullets in the linear and nonlinear regimes,” Phys. Rev. Lett. 105(25), 253901 (2010). [CrossRef]
11. H. Li, X. Huang, Q. Cao, Y. Zhao, P. Li, C. Wan, and A. Chong, “Generation of three-dimensional versatile vortex linear light bullets (Invited Paper),” Chin. Opt. Lett. 15(3), 030009 (2017). [CrossRef]
12. W. P. Zhong, B. Milivoj, and Y. Q. Zhang, “Three-dimensional localized Airy-Laguerre-Gaussian wave packets in free space,” Opt. Express 23(18), 23867–23876 (2015). [CrossRef]
13. D. Bongiovanni, B. Wetzel, Y. Hu, Z. Chen, and R. Morandotti, “Optimal compression and energy confinement of optical Airy bullets,” Opt. Express 24(23), 26454–26463 (2016). [CrossRef]
14. Y. L. Peng, B. Chen, P. Xi, M. L. Zhou, L. P. Zhang, D. D. Li, and D. M. Deng, “Self-accelerating Airy-Ince-Gaussian and Airy-Helical-Ince-Gaussian light bullets in free space,” Opt. Express 24(17), 18973–18985 (2016). [CrossRef]
15. J. F. Wu, F. Q. Huang, Y. Q. Chen, L. Cheng, L. F. Li, and D. Deng, “Three-dimensional localized Airy elegant Laguerre–Gaussian light bullets in free space,” J. Opt. Soc. Am. B 34(4), 808–813 (2017). [CrossRef]
16. N. K. Efremidis, Z. Chen, M. Segev, and D. N. Christodoulides, “Airy beams and accelerating waves: an overview of recent advances,” Optica 6(5), 686–701 (2019). [CrossRef]
17. J. Durnin, “Exact solutions for nondiffracting beams. I. The scalar theory,” J. Opt. Soc. Am. A 4(4), 651–654 (1987). [CrossRef]
18. J. Durnin, J. J. Miceli, and J. H. Eberly, “Diffraction-free beams,” Phys. Rev. Lett. 58(15), 1499–1501 (1987). [CrossRef]
19. S. M. Wang, Q. Lin, and X. H. Lu, “Super-diffraction limiting propagation of cos beams,” Appl. Laser 14(2), 49–51 (1994).
20. S. Wang, Q. Lin, and X. Lu, “Realization of super-diffraction limiting propagation by cos beams,” Optik 100(1), 8–10 (1995).
21. Z. Jiang, “Truncation of a two-dimensional nondiffracting cos beam,” J. Opt. Soc. Am. A 14(7), 1478–1481 (1997). [CrossRef]
22. S. M. Wang, Q. Lin, and X. Q. Jiang, “Cosine-Gauss beam,” Acta. Photonica Sinica 28(4), 367–370 (1999).
23. L. W. Casperson and A. T. Anthony, “Hermite–sinusoidal-Gaussian beams in complex optical systems,” J. Opt. Soc. Am. A 15(4), 954–961 (1998). [CrossRef]
24. C. G. V. Julio and A. B. Miguel, “Helmholtz–Gauss waves,” J. Opt. Soc. Am. A 22(2), 289–298 (2005). [CrossRef]
25. A. Bencheikh, S. Chabou, O. C. Boumeddine, H. Bekkis, A. Benstiti, L. Beddiaf, and W. Moussaoui, “Cosine beam: diffraction-free propagation and self-healing,” J. Opt. Soc. Am. A 37(11), C7–C14 (2020). [CrossRef]
26. X. Q. Wang, Z. P. Ke, and B. D. Lü, “Propagation properties of cosine-Gaussian beams,” Laser Technol. 25(2), 121–125 (2001).
27. T. E. Halil and B. Yahya, “Analysis of reciprocity of cos-Gaussian and cosh-Gaussian laser beams in a turbulent atmosphere,” Opt. Express 12(20), 4659–4674 (2004). [CrossRef]
28. T. E. Halil, “Hermite-cosine-Gaussian laser beam and its propagation characteristics in turbulent atmosphere,” J. Opt. Soc. Am. A 22(8), 1527–1535 (2005). [CrossRef]
29. G. Q. Zhou and X. X. Chu, “Propagation of a partially coherent cosine-Gaussian beam through an ABCD optical system in turbulent atmosphere,” Opt. Express 17(13), 10529–10534 (2009). [CrossRef]
30. W. J. Wang and Z. S. Wu, “Propagation of annular cos-Gaussian beams through turbulence,” J. Opt. Soc. Am. A 35(7), 1165–1172 (2018). [CrossRef]
31. T. Bin, “Hermite-cosine-Gaussian beams propagating in uniaxial crystals orthogonal to the optical axis,” J. Opt. Soc. Am. A 26(12), 2480–2487 (2009). [CrossRef]
32. C. L. Ding, K. Olga, Y. T. Zhang, and L. Z. Pan, “Cosine-Gaussian correlated Schell-model pulsed beams,” Opt. Express 22(1), 931–942 (2014). [CrossRef]
33. Z. R. Mei and O. Korotkova, “Cosine-Gaussian Schell-model sources,” Opt. Lett. 38(14), 2578–2580 (2013). [CrossRef]
34. H. F. Xu, Z. Zhang, J. Qu, and W. Huang, “Propagation factors of cosine-Gaussian-correlated Schell-model beams in non-Kolmogorov turbulence,” Opt. Express 22(19), 22479–22489 (2014). [CrossRef]
35. L. Z. Pan, C. L. Ding, and H. X. Wang, “Diffraction of cosine-Gaussian-correlated Schell-model beams,” Opt. Express 22(10), 11670–11679 (2014). [CrossRef]
36. Z. R. Mei and O. Korotkova, “Electromagnetic cosine-Gaussian Schell-model beams in free space and atmospheric turbulence,” Opt. Express 21(22), 27246–27259 (2013). [CrossRef]
37. Z. R. Mei, E. Shchepakina, and O. Korotkova, “Propagation of cosine-Gaussian-correlated Schell-model beams in atmospheric turbulence,” Opt. Express 21(15), 17512–17519 (2013). [CrossRef]
38. J. Lin, J. Dellinger, P. Genevet, B. Cluzel, F. de Fornel, and F. Capasso, “Cosine-Gauss plasmon beam: a localized long-range nondiffracting surface wave,” Phys. Rev. Lett. 109(9), 093904 (2012). [CrossRef]
39. X. S. Hao, X. Feng, P. Zhao, F. Liu, K. Y. Cui, W. Zhang, and Y. D. Huang, “Polarization-controllably launching localized cosine-Gauss beam with spatially varied metallic nano-apertures,” Opt. Express 27(16), 22053–22073 (2019). [CrossRef]
40. X. Q. He, T. G. Ning, R. J. Li, P. Li, J. J. Zheng, and J. Li, “Dynamical manipulation of cosine-Gauss beams in a graphene plasmonic waveguide,” Opt. Express 25(12), 13923–13932 (2017). [CrossRef]
41. K. Xiao, S. B. Wei, C. J. Min, G. H. Yuan, S. W. Zhu, T. Lei, and X. C. Yuan, “Dynamic cosine-Gauss plasmonic beam through phase control,” Opt. Express 22(11), 13541–13546 (2014). [CrossRef]
42. S. H. Fu, Y. Tsur, J. Y. Zhou, L. Shemer, and A. Arie, “Propagation dynamics of nonspreading cosine-Gauss water-wave pulses,” Phys. Rev. Lett. 115(25), 254501 (2015). [CrossRef]
43. S. H. Fu, Y. Tsur, J. Y. Zhou, L. Shemer, and A. Arie, “Self-similar propagation of Hermite-Gauss water-wave pulses,” Phys. Rev. E 93(1), 013127 (2016). [CrossRef]
44. N. K. Efremidis, “Airy trajectory engineering in dynamic linear index potentials,” Opt. Lett. 36(15), 3006–3008 (2011). [CrossRef]
45. W. Liu, D. N. Neshev, I. V. Shadrivov, A. E. Miroshnichenko, and Y. S. Kivshar, “Plasmonic Airy beam manipulation in linear optical potentials,” Opt. Lett. 36(7), 1164–1166 (2011). [CrossRef]
46. H. Zhong, Y. Zhang, M. R. Belic, C. Li, F. Wen, Z. Zhang, and Y. Zhang, “Controllable circular Airy beams via dynamic linear potential,” Opt. Express 24(7), 7495–7506 (2016). [CrossRef]
47. X. Huang, Z. Deng, and X. Fu, “Dynamics of finite energy Airy beams modeled by the fractional Schrödinger equation with a linear potential,” J. Opt. Soc. Am. B 34(5), 976–982 (2017). [CrossRef]
48. M. Newstein and K. Lin, “Laguerre-Gaussian periodically focusing beams in a quadratic index medium,” IEEE J. Quantum Electron. 23(5), 481–482 (1987). [CrossRef]
49. J. C. Gutiérrez-Vega and M. A. Bandres, “Ince-Gaussian beams in a quadratic-index medium,” J. Opt. Soc. Am. A 22(2), 306–309 (2005). [CrossRef]
50. Y. Q. Zhang, M. R. Belić, L. Zhang, W. P. Zhong, D. Y. Zhu, R. M. Wang, and Y. P. Zhang, “Periodic inversion and phase transition of finite energy Airy beams in a medium with parabolic potential,” Opt. Express 23(8), 10467–10480 (2015). [CrossRef]
51. X. Peng, Y. Peng, L. Zhang, D. Li, and D. Deng, “Reversed Airy Gaussian and Airy Gaussian vortex light bullets in harmonic potential,” Laser Phys. Lett. 14(5), 055002 (2017). [CrossRef]
52. Y. Zhang, X. Liu, M. R. Belic, W. Zhong, Y. Zhang, and M. Xiao, “Propagation dynamics of a light beam in a fractional Schrödinger equation,” Phys. Rev. Lett. 115(18), 180403 (2015). [CrossRef]
53. Y. Q. Zhang, X. Liu, M. R. Belić, W. P. Zhong, S. M. Petrović, and Y. Zhang, “Automatic Fourier transform and self-Fourier beams due to parabolic potential,” Ann. Phys. 363, 305–315 (2015). [CrossRef]
54. Y. Q. Zhang, X. Liu, M. R. Belić, W. P. Zhong, F. Wen, and Y. P. Zhang, “Anharmonic propagation of two-dimensional beams carrying orbital angular momentum in a harmonic potential,” Opt. Lett. 40(16), 3786–3789 (2015). [CrossRef]
55. R. Zhao, F. Deng, W. Yu, J. Huang, and D. Deng, “Propagation properties of Airy-Gaussian vortex beams through the gradient index medium,” J. Opt. Soc. Am. A 33(6), 1025–1031 (2016). [CrossRef]
56. F. Liu, J. W. Zhang, W. P. Zhong, M. R. Belić, Y. Zhang, Y. P. Zhang, F. L. Li, and Y. Q. Zhang, “Manipulation of Airy beams in dynamic parabolic potentials,” Ann. Phys. 532(4), 1900584 (2020). [CrossRef]
57. Z. Pang and D. Deng, “Propagation properties and radiation forces of the Airy Gaussian vortex beams in a harmonic potential,” Opt. Express 25(12), 13635–13647 (2017). [CrossRef]
58. X. Peng, Y. Peng, D. Li, L. Zhang, J. Zhuang, F. Zhao, X. Chen, X. Yang, and D. Deng, “Propagation properties of spatiotemporal chirped Airy Gaussian vortex wave packets in a quadratic index medium,” Opt. Express 25(12), 13527–13538 (2017). [CrossRef]
59. L. Feng, J. Zhang, Z. Pang, L. Wang, T. Zhong, X. Yang, and D. Deng, “Propagation properties of the chirped Airy beams through the gradient-index medium,” Opt. Commun. 402, 60–65 (2017). [CrossRef]
60. Y. Qian and S. Zhang, “Manipulation of accelerating beam in gradient potentials,” Opt. Commun. 435, 409–412 (2019). [CrossRef]
61. Y. Lan, F. Hu, and Y. Qian, “Generation of spirally accelerating optical beams,” Opt. Lett. 44(8), 1968–1971 (2019). [CrossRef]
62. L. Zhang, K. Liu, H. Zhong, J. Zhang, Y. Li, and D. Fan, “Effect of initial frequency chirp on Airy pulse propagation in an optical fiber,” Opt. Express 23(3), 2566–2576 (2015). [CrossRef]
63. D. Mendlovic and H. M. Ozaktas, “Fractional Fourier transforms and their optical implementation: I,” J. Opt. Soc. Am. A 10(9), 1875–1881 (1993). [CrossRef]
64. D. Lu, W. Hu, Y. Zheng, Y. Liang, L. Cao, S. Lan, and Q. Guo, “Self-induced fractional Fourier transform and revivable higher-order spatial solitons in strongly nonlocal nonlinear media,” Phys. Rev. A 78(4), 043815 (2008). [CrossRef]