1. Introduction
During the past decades, investigation of Airy beams [
1,
2,
3,
4,
5,
6,
7,
8,
9] have been one of the hot topics in the area of optics due to their novel propagation properties, e.g., ballistic motion [
10,
11], self-healing [
12], etc. In linear regime, mechanisms of Airy beams are widely used in particle clearing [
13], light bullets [
14], surface plasmons [
15,
16,
17], and electron beams [
18,
19]. On the contrary, nonlinear manipulation of Airy beams is another novel problem, including nonlinear generation of Airy beams [
20,
21], propagation of Airy beams in nonlinear media [
22,
23,
24,
25,
26,
27,
28,
29], solitons generated form Airy beams [
30] and spatiotemporal Airy light bullets [
31,
32,
33].
In nonlinear media, interactions between Airy beams have also been extensively studied in recent years. In particular, interactions of Airy beams were studied in photorefractive crystal [
34,
35,
36,
37,
38,
39], Kerr media [
40,
41,
42], optical fiber [
43,
44,
45], quadratic nonlinear medium [
46], fractional media [
47,
48], photonic lattices [
49], as well as nonlocal nonlinear media [
50,
51,
52,
53]. Compared with interactions of conventional beams [
54], many anomalous interactions of Airy beams have been reported [
51,
52,
53].
In nonlinear media, besides nonlinear self-focusing effect, linear diffraction is another effect induced by optical beams themselves during propagation [
55]. Apart from quadratic diffraction, high-order diffraction should not be neglected [
56]. In particular, fourth-order diffraction has deep impacts on nonlinear optics, e.g., self-similar propagation of optical pulses [
57], modulation instability [
58], pure-quartic soliton laser [
59], and soliton states [
60,
61,
62,
63,
64,
65], etc.
The effects of third-order [
66,
67] and fourth-order [
68,
69,
70] diffraction/dispersion on the propagation of Airy pulses/beams have been investigated in detail. In our previous work, we have also studied interactions of Airy beams in local Kerr, saturable, and nonlocal media with fourth-order diffraction [
71]. We demonstrated that stable bound states of Airy solitons pairs can be obtained in nonlocal medium with fourth-order diffraction which are always repulsive in local media [
71].
In this paper, based on our previous works [
52,
71], we demonstrate that the interactions of solitons generated from two in-phase Airy beams in Kerr nonlinear media can be controlled by a fundamental Gaussian beam and fourth-order diffraction. Numerical results with split-step Fourier transform method reveal the interactions of Airy beam are sensitive to the beam intervals, coupling constant of fourth-order diffraction as well as the relative phase and amplitude of the fundamental Gaussian beam. In particular, many anomalous interactions are obtained and physical mechanisms are explained in detail. Stable bound states of breathing Airy soliton pairs can be formed with appropriate beam parameters and coupling constants of fourth-order diffraction.
2. Physical model and basic equations
Considering the propagation of an optical beams in nonlinear media with fourth-order diffraction, the slowly varying envelope of the beam
is described by the following normalized nonlinear Schrödinger equation[
61,
62],
here
x and
z denote transverse and longitude spatial coordinates [
40,
41], and
represents the strength of fourth-order diffraction [
61,
62,
63], respectively. When
, the fourth-order diffraction is normal, whereas, it is anomalous when
[
61,
62,
63]. The nonlinear refractive index change of the Kerr media
can be written as
In this work, we concentrate on the manipulation of interactions of Airy beams with the help of a fundamental Gaussian beam and fourth-order diffraction. Thus we assume the incident beam is composed of two in-phase shifted counter propagating Airy beams, but with an additional fundamental Gaussian beam located in the middle with a relative phase between them [
52],
where
A represents the amplitude of the Airy beams and we set
to ensure solitons can be generated from Airy beams,
B is the parameter controlling the beam separations (intervals),
C denotes the amplitude of the fundamental Gaussian beam [
52], and
is the decaying factor of the Airy beams [
1], respectively.
is the parameter controlling the phase shift, in particular,
(
) describes Airy beams and fundamental Gaussian beam are in-phase (out of phase) [
52].
The interactions of the beams can be investigated numerically by integrating Eq. (1) directly with split-step Fourier transform method. In the following, with fourth-order diffraction, we will demonstrated the interactions of Airy beams modulated by in-phase and out-of-phase fundamental Gaussian beams, respectively.
3. Interactions of Airy beams with in-phase Gaussian beam and fourth-order diffraction
Firstly, we consider the interactions of Airy beams and in-phase (
) Gaussian beam with fourth-order diffraction. In
Figure 1, we display the numerical results of the interactions with some different beam intervals
and fourth-order diffraction coupling constants
in the case of weak Gaussian beam (
). Previous work has shown that dynamics of interactions depend crucially on the amplitude of the Gaussian beam [
52].
For comparison, we also re-do some previous results. When
and
, it represents the interactions of two in-phase Airy beams in Kerr media with only quadratic diffraction [
40,
41,
51,
52,
71] and attractions between beams are always dominant. It is obvious that bound states of breathing soliton pairs with different breathing periods are obtained [
Figure 1(a1-c1)]. The strongest attraction occurs at
[
Figure 1(b1)] due to strong self-trapping effect [
40,
41,
51,
52]. With a in-phase (
) weak fundamental Gaussian beam (
), the attraction between Airy beams are enhanced and the corresponding breathing periods are reduced significantly [
Figure 1(a2-c2)]. The physical mechanisms for above results have been explained extensively in previous works [
52].
The dynamics of bound states can be modulated by fourth-order diffraction. In particular, without fundamental Gaussian beam (C=0), the interactions between both in-phase and out-of-phase Airy beams under fourth-order diffraction have been studied previously [
71]. Thus in this paper we focus on the condition of
when fourth-order diffraction is taken into account (
). For normal fourth-order diffraction (
), with an in-phase Gaussian beam, as shown in
Figure 1(a3-c3), the interactions of Airy beams exhibit repulsions for all beam intervals due to the fact that the diffraction effect of the beams is enhanced. The nonlinearity is unable to overcome the strong diffraction induced by normal fourth-order diffraction, and therefore no stable bound state can be formed [
71]. On the other hand, when fourth-order diffraction is anomalous (
), situations become different. Comparing with case without fourth-order diffraction
[
Figure 1(a2-c2)], the attractions between Airy beams are significantly strengthened with breathing periods of all the bound states become smaller [
Figure 1(a4-c4)]. The reason is that the diffraction effect of the beams is weakened with anomalous fourth-order diffraction, so nonlinearity will reinforce the interactions between Airy beams [
71].
In general, the beam interval
B is an important parameter in controlling the interactions of Airy beams. When
B is larger, e.g.,
, only two parallel solitons can be generated from two well-separated Airy beams [
40,
41]. Apart from the results obtained above, anomalous interactions happen when
. As illustrated in
Figure 1, an in-phase fundamental Gaussian beam always enhances the attractions of two in-phase Airy beams when
. When increasing the amplitude
C gradually, attraction between the Airy beams with
decreases firstly and subsequently increases [
52], as shown in
Figure 2(a2-c2). This interesting phenomena can be explained as follow: the nonlinearity induced by the intensity of the incident beam [Eq. (3)] decreases firstly and then increases subsequently, leading to the weakness and then enhancement of the attractions [
52]. Similar with
Figure 1, normal (
) and anomalous (
) fourth-order diffraction will weaken and strengthen the attractions of the Airy beams, result in the repulsions [
Figure 2(a1-c1)] and bound states [
Figure 2(a3-c3)] with smaller breathing periods, respectively.
When the Gaussian beam is strong light, the amplitude of the fundamental Gaussian beam
C is large enough. Similarly, we show in
Figure 3 interactions of airy beams with different beam intervals
and fourth-order diffraction coupling constants
. When
, a single solitons generated from fundamental Gaussian beam is located in the cental of two solitons shedding from Airy beams [
Figure 3(a2-d2)]. Different from the case of enhancement of the attractions with in-phase weak light Gaussian beam [
Figure 1(a2-c2)], repulsions appear for all beam intervals [
Figure 3(a2-d2)]. For smaller beam intervals
[
Figure 3(a2)] and
[
Figure 3(b2)], repulsions appear after a very short propagation distance. Although bound states are formed initially, for larger beam intervals
[
Figure 3(c2)] and
[
Figure 3(d2)], repulsions overcome attractions eventually [
41,
52].
Surprisingly, fourth-order diffraction plays an unique role in the interactions of Airy beam with an in-phase strong light Gaussian beam. As shown in
Figure 3, repulsions are weakened for
[
Figure 3(a1-d1)] and strengthened for
[
Figure 3(a3-d3)]. This interesting result is quite different with the case of weak light Gaussian beam [
Figure 1 and 2], where normal (anomalous) fourth-order diffraction always strengthen (weaken) the repulsions. The anomalous phenomena is very similar with the interactions of two in-phase Airy beams in saturable nonlinear media with fourth-order diffraction [
71].
4. Interactions of Airy beams with out-of-phase Gaussian beam and fourth-order diffraction
Next, we study the interactions of Airy beams and out-of-phase (
) fundamental Gaussian beam in nonlinear Kerr media with fourth-order diffraction. Comparing
Figure 1(a1-c1), when beam intervals are smaller (
) and without fourth-order diffraction (
), weak light out-of-phase fundamental Gaussian beam weakens the attractions of Airy beams with breathing periods increase obviously [
52], as shown in
Figure 4(a2-c2). This is the general case that out-of-phase Gaussian beam provides a repulsive force between Airy beams [
52].
Similar with results of
Figure 1 and 3, dynamics of interactions of Airy beams are shown in
Figure 4(a1-c1) and 4(a3-c3) when normal and anomalous fourth-order diffraction is used, respectively. Specifically, with normal fourth-order diffraction (
), attractions become weaker and repulsions appear for all the beam intervals. On the contrary, attractions become stronger and breathing periods of the bound states decrease for anomalous fourth-order diffraction (
).
In the case of the beam interval
, without fourth-order diffraction (
) and fundamental Gaussian beam (
), stationary bound state of two in-phase Airy beams is formed [
Figure 2(a2)]. When an out-of-phase Gaussian beam with amplitude
is added, repulsion appears between the two Airy beams and one fundamental solitons generated from Gaussian beam is located in the center, as shown in
Figure 5(a2). However, repulsion will be weakened [
Figure 5(b2)] when the amplitude of Gaussian beam
C increases, and bound state of breathing mode can be even formed [
Figure 5(c2)] with appropriate intensity of Gaussian beam. We emphasize that the anomalous interaction demonstrated in
Figure 5(c2) is different from conventional interactions between out-of-phase beams, where repulsions are always dominant [
54].
Similar with the interactions of Airy beam with an in-phase strong light Gaussian beam shown in
Figure 3 or interactions of two in-phase Airy beams in saturable nonlinear media [
71], fourth-order diffraction plays the same role in the interactions between Airy beam and out-of-phase weak light Gaussian beam when the beam interval is
. As shown in
Figure 5(a1-b1), repulsions are also weakened (or attractions are enhanced) and breathing period decreases for
[
Figure 5(c1)] when fourth-order diffraction is normal. On the other hand, repulsions are enhanced for all the interactions with anomalous fourth-order diffraction [
Figure 5(a3-c3)].
In the case of out-of-phase strong light Gaussian beam, for beam intervals
and
, similar with
Figure 4(a2-b2), we numerically check that the breathing periods of the bound states become larger (not shown). When
, as shown in
Figure 6(a2), repulsion appears instead of attraction [
Figure 4(c2)] when
. Similar result happens to the case of
, repulsion emerges again [
Figure 6(b2)] despite of bound state can be formed with appropriate out-of-phase fundamental Gaussian beam [
Figure 5(c2)].
Similar with
Figure 3 and 5, anomalous fourth-order diffraction (
) strengthens the repulsions of the interactions, as shown in
Figure 6(a3) and 6(b3). However, normal fourth-order diffraction (
) weakens the repulsions for
[
Figure 6(b1)], and stationary bound state can even be formed for
[
Figure 6(a1)]. These impacts of fourth-order diffraction are different with outcomes obtained in
Figure 1, 2 and 4.
Figure 1.
(Color online) The intensity trajectories of interactions between Airy beams and in-phase () weak light Gaussian beams. The beam intervals are B=0 (a1-a4), B=1 (b1-b4), and B=2 (c1-c4). The fourth-order diffraction coupling constants are , , and for (a3-c3), (a1-c2), and (a4-c4), respectively. The amplitudes are .
Figure 1.
(Color online) The intensity trajectories of interactions between Airy beams and in-phase () weak light Gaussian beams. The beam intervals are B=0 (a1-a4), B=1 (b1-b4), and B=2 (c1-c4). The fourth-order diffraction coupling constants are , , and for (a3-c3), (a1-c2), and (a4-c4), respectively. The amplitudes are .
Figure 2.
(Color online) The intensity trajectories of interactions between Airy beams and in-phase () weak light Gaussian beams. The beam intervals are B=3. The fourth-order diffraction coupling constants are , , and for (a1-c1), (a2-c2), and (a3-c3), respectively. The amplitudes are .
Figure 2.
(Color online) The intensity trajectories of interactions between Airy beams and in-phase () weak light Gaussian beams. The beam intervals are B=3. The fourth-order diffraction coupling constants are , , and for (a1-c1), (a2-c2), and (a3-c3), respectively. The amplitudes are .
Figure 3.
(Color online) The intensity trajectories of interactions between Airy beams and in-phase () strong light Gaussian beams. The beam intervals are B=0 (a1-a3), B=1 (b1-b3), B=2 (c1-c3) and B=3 (d1-d3). The fourth-order diffraction coupling constants are , , and for (a1-d1), (a2-d2), and (a3-d3), respectively. The amplitudes are .
Figure 3.
(Color online) The intensity trajectories of interactions between Airy beams and in-phase () strong light Gaussian beams. The beam intervals are B=0 (a1-a3), B=1 (b1-b3), B=2 (c1-c3) and B=3 (d1-d3). The fourth-order diffraction coupling constants are , , and for (a1-d1), (a2-d2), and (a3-d3), respectively. The amplitudes are .
Figure 4.
(Color online) The intensity trajectories of interactions between Airy beams and out-of-phase () weak light Gaussian beams. The beam intervals are B=0 (a1-a3), B=1 (b1-b3), and B=2 (c1-c3). The fourth-order diffraction coupling constants are , , and for (a1-c1), (a2-c2), and (a3-c3), respectively. The amplitudes are .
Figure 4.
(Color online) The intensity trajectories of interactions between Airy beams and out-of-phase () weak light Gaussian beams. The beam intervals are B=0 (a1-a3), B=1 (b1-b3), and B=2 (c1-c3). The fourth-order diffraction coupling constants are , , and for (a1-c1), (a2-c2), and (a3-c3), respectively. The amplitudes are .
Figure 5.
(Color online) The intensity trajectories of interactions between Airy beams and out-of-phase () weak light Gaussian beams. The beam intervals are B=3. The fourth-order diffraction coupling constants are , , and for (a1-c1), (a2-c2), and (a3-c3), respectively. The amplitudes are .
Figure 5.
(Color online) The intensity trajectories of interactions between Airy beams and out-of-phase () weak light Gaussian beams. The beam intervals are B=3. The fourth-order diffraction coupling constants are , , and for (a1-c1), (a2-c2), and (a3-c3), respectively. The amplitudes are .
Figure 6.
(Color online) The intensity trajectories of interactions between Airy beams and out-of-phase () strong light Gaussian beams. The beam intervals are B=2 (a1-a3) and B=3 (b1-b3). The fourth-order diffraction coupling constants are , , and for (a1-b1), (a2-b2), and (a3-b3), respectively. The amplitudes are .
Figure 6.
(Color online) The intensity trajectories of interactions between Airy beams and out-of-phase () strong light Gaussian beams. The beam intervals are B=2 (a1-a3) and B=3 (b1-b3). The fourth-order diffraction coupling constants are , , and for (a1-b1), (a2-b2), and (a3-b3), respectively. The amplitudes are .