We investigate numerically the interactions of two in-phase and out-of-phase Airy beams and nonlinear (NL) accelerating beams in Kerr and saturable NL media, in one transverse dimension. We find that bound and unbound soliton pairs, as well as single solitons, can form in such interactions. If the interval between two incident beams is large relative to the width of their first lobes, the generated soliton pairs just propagate individually and do not interact. However, if the interval is comparable to the widths of the maximum lobes, the pairs interact and display varied behavior. In the in-phase case, they attract each other and exhibit stable bound, oscillating, and unbound states, after shedding some radiation initially. In the out-of-phase case, they repel each other and, after an initial interaction, fly away as individual solitons. While the incident beams display acceleration, the solitons or soliton pairs generated from those beams do not.

We investigate numerically interactions between two in-phase or out-of-phase Airy beams and nonlinear accelerating beams in Kerr and saturable nonlinear media in one transverse dimension. We discuss different cases in which the beams with different intensities are launched into the medium, but accelerate in opposite directions. Since both the Airy beams and nonlinear accelerating beams possess infinite oscillating tails, we discuss interactions between truncated beams, with finite energies. During interactions we see solitons and soliton pairs generated that are not accelerating. In general, the higher the intensities of interacting beams, the easier to form solitons; when the intensities are small enough, no solitons are generated. Upon adjusting the interval between the launched beams, their interaction exhibits different properties. If the interval is large relative to the width of the first lobes, the generated soliton pairs just propagate individually and do not interact much. However, if the interval is comparable to the widths of the maximum lobes, the pairs strongly interact and display varied behavior.

We study periodic inversion and phase transition of normal, displaced, and chirped finite energy Airy beams propagating in a parabolic potential. This propagation leads to an unusual oscillation: for half of the oscillation period the Airy beam accelerates in one transverse direction, with the main Airy beam lobe leading the train of pulses, whereas in the other half of the period it accelerates in the opposite direction, with the main lobe still leading – but now the whole beam is inverted. The inversion happens at a critical point, at which the beam profile changes from an Airy profile to a Gaussian one. Thus, there are two distinct phases in the propagation of an Airy beam in the parabolic potential – the normal Airy and the single-peak Gaussian phase. The length of the single-peak phase is determined by the size of the decay parameter: the smaller the decay, the smaller the length. A linear chirp introduces a transverse displacement of the beam at the phase transition point, but does not change the location of the point. A quadratic chirp moves the phase transition point, but does not affect the beam profile. The two-dimensional case is discussed briefly, being equivalent to a product of two one-dimensional cases.

We investigate the propagation of one-dimensional and two-dimensional (1D, 2D) Gaussian beams in the fractional Schrödinger equation (FSE) without a potential, analytically and numerically. Without chirp, a 1D Gaussian beam splits into two nondiffracting Gaussian beams during propagation, while a 2D Gaussian beam undergoes conical diffraction. When a Gaussian beam carries linear chirp, the 1D beam deflects along the trajectories z = ±2(x − x₀), which are independent of the chirp. In the case of 2D Gaussian beam, the propagation is also deflected, but the trajectories align along the diffraction cone "Equation missing"<!-- image only, no MathML or LaTex --> and the direction is determined by the chirp. Both 1D and 2D Gaussian beams are diffractionless and display uniform propagation. The nondiffracting property discovered in this model applies to other beams as well. Based on the nondiffracting and splitting properties, we introduce the Talbot effect of diffractionless beams in FSE.

We analytically and numerically investigate an anharmonic propagation of two-dimensional beams in a harmonic potential. We pick noncentrosymmetric beams of common interest that carry orbital angular momentum. The examples studied include superposed Bessel–Gauss (BG), Laguerre–Gauss (LG), and circular Airy (CA) beams. For the BG beams, periodic inversion, phase transition, and rotation with periodic angular velocity are demonstrated during propagation. For the LG and CA beams, periodic inversion and variable rotation are still there but not the phase transition. On the whole, the “center of mass” and the orbital angular momentum of a beam exhibit harmonic motion, but the motion of the beam intensity distribution in detail is subject to external and internal torques and forces, causing it to be anharmonic. Our results are applicable to other superpositions of finite circularly asymmetric beams.

We investigate controllable spatial modulation of circular autofocusing Airy beams, under action of different dynamic linear potentials, both theoretically and numerically. We introduce a novel treatment method in which the circular Airy beam is represented as a superposition of narrow azimuthally-modulated one-dimensional Airy beams that can be analytically treated. The dynamic linear potentials are appropriately designed, so that the autofocusing effect can either be weakened or even eliminated when the linear potential exerts a “pulling” effect on the beam, or if the linear potential exerts a “pushing” effect, the autofocusing effect can be greatly strengthened. Numerical simulations agree with the theoretical results very well.