Abstract
We analyse theoretically the spatiotemporal dynamics of two incoherent counterpropagating Airy beams interacting in a photorefractive crystal under focusing conditions. For a large enough nonlinearity strength the interaction between the two Airy beams leads to lightinduced waveguiding. The stability of the waveguide is determined by the crystal length, the nonlinearity strength and the beam’s intensities and is improved when comparing to the situation using Gaussian beams. We further identify the threshold above which the waveguide is no longer static but evolves dynamically either timeperiodically or even chaotically. Above the stability threshold, each Airysoliton moves erratically between privileged output positions that correspond to the spatial positions of the lobes of the counterpropagating Airy beam. These results suggest new ways of creating dynamically varying waveguides, optical logic gates and chaosbased computing.
Introduction
Instabilities, selfoscillations and chaos are fundamental processes in nonlinear optics. Multiple beams’ interactions in nonlinear media, even without external feedback, can give rise to beam selftrapping and spatial solitons that may further destabilise to spatiotemporal dynamics and then, eventually, chaos^{1,2,3}. Multiple parameters, such as the optical intensity or the misalignment of the interacting beams^{4,5}, enable to control the sequence of bifurcations from stationary dynamics to deterministic chaos^{6}.
Interestingly the onset of spatiotemporal instabilities observed for various beams’ configurations in different nonlinear media presents the same evolution pattern: initial diffraction, collapse to the soliton shape, then timeperiodic dynamics to chaotic instabilities, where the interacting solitons rotate and twist around each other in an erratic way^{7,8,9,10,11}. Solitons are selftrapped light beams propagating without change in a diffractive nonlinear medium. They are generated by the photoinduced refractive index in media, as in Kerr or photorefractive crystals or in nematic liquid crystals^{10,12,13}. The interacting solitons enable various optically induced waveguiding structures in the nonlinear media. The control of the mutual exchange of energy between interacting solitons enables to create alloptical guiding, dividing and switching devices^{14,15} and even over large distances in different media^{16}. Apart from Gaussian beams, other diffractive beam profiles such as optical vortices present solitonic behaviour under selffocusing conditions and exhibit similar dynamical routes to instabilities^{17}.
Recent works have shown the possibility to induce spatial solitons from selffocusing of ideally nondiffractive beam profiles including optical Airy beams^{18,19,20}. As a truncated solution of the ideal Airy waveform, the optical Airy beam has the advantage of combining the parabolic trajectory and selfhealing properties of the Airy wave solution over a finite distance with the diffractive beam properties for larger propagation distances. Airy beams have been extensively studied in the recent years. First discovered in quantum mechanics^{21}, the Airy wave packet has been suggested as a nonspreading and selfhealing solution of the Schrödinger equation with a parabolic propagation. Due to these unique properties, the optical analogy has been of a great interest. In 2007, Christodoulides et al. has generated the first optical Airy beam using a spatial light modulator (SLM)^{22}. The use of an SLM or other generation methods, as a photonic crystal^{23}, allows for controlling the exact linear propagation trajectory of the beam^{24,25}. Applications of the Airy beams are very large, from optical micromanipulation of particles^{26}, to laser processing^{27} or optical routing^{28}, including plasmonic circuitry and surface tweezers using Airy plasmons^{29}.
As already mentioned, under highfocusing conditions the Airy beam undergoes solitonlike behaviour. The main part of the beam power is focused in a socalled offshooting soliton while the remaining part propagates as in the linear case. Nonstationary dynamics of a single Airysoliton have been shown in Kerr media and called moving solitons^{30}. The collision of two Airy beams also suggests a large variety of interaction schemes. The solitonic interactions have been demonstrated in the spatial domain using copropagating Airy beams^{31,32,33,34} and in the temporal domain using Airy pulses^{35,36}. Recently we have studied the counterpropagating (CP) configuration in photorefractive media^{37}. First results have shown more complex stationary waveguide structures than those induced by interacting Gaussian beams. A single Airy beam leads to waveguiding structures with multiple outputs. The additional interactions induced by a CP beam allows for achieving complex waveguiding structures that would otherwise require the CP interactions of more than two Gaussian beams. But these stationary structures only exist for a limited range of the nonlinearity strength. As will be shown here, by increasing the nonlinear coupling strength, spatiotemporal dynamics appear that result from the interactions of the two counterpropagating (CP) Airy beams and that differ from what is known from CP Gaussian beams’ systems.
In this paper we analyse numerically the spatiotemporal dynamics of two incoherent counterpropagating Airy beams interacting in a photorefractive crystal under focusing conditions. When a positive external electrical field is applied on the crystal, multiple waveguiding structures are photoinduced. If we increase the nonlinearity of the system, we demonstrate the existence of a threshold curve above which non steadystate dynamics appear. The threshold leading to spatiotemporal dynamics can be tuned by the crystal length, the nonlinearity strength and the beams’ intensities. For similar operating conditions this threshold is larger for interacting Airy beams when compared to interacting Gaussian beams, hence demonstrating the larger stability range of the photoinduced Airy waveguides. Above the threshold the position of the offshooting soliton varies periodically in time. The period is of the same order of magnitude than the material nonlinear optics timescale and the amplitude is determined by the Airy properties, hence both properties can be controlled by the beam interaction scheme. When further increasing the nonlinear strength and/or the crystal length, this timeperiodic dynamics bifurcates to a chaoticlike dynamics of the offshooting soliton. The erratic motion of the soliton encompasses privileged positions that match the input positions of the multiple lobes of the counterpropagating Airy beams. The engineering of the Airy beam properties therefore allows to modify the topology of the resulting chaotic motion. These findings not only address the important issue of the Airy waveguiding stability but also suggest new ways of creating dynamically varying waveguides, optical logic gates and chaosbased computing.
Results
Modelisation and simulation scheme of two counterpropagating Airy beams
To study the nonlinear interactions of two CP Airy beams in a photorefractive medium, we consider the propagation of two identical onedimensional Airy beams along the longitudinal zaxis. The two Airy beams are initially injected at each side of the crystal and have both the lobe size and aperture parameter a, as given by the following equations^{8}:
where F_{0} and B_{0} correspond to the electric field amplitudes of respectively the forward beam F and the backward beam B, Ai represents the Airy function, x_{0} an arbitrary transverse scale and a the truncation factor. Counterpropagating beams are herein defined by their opposite propagation directions along the longitudinal zaxis [Fig. 1]. Each Airy beam is formed of successive lobes, the first order (centered at x = −x_{0}) being the main lobe. The nonlinear propagation of these two incoherent CP beams can be expressed as follows:
where is the nonlinear photorefractive coupling strength (r_{eff} is the effective component of the electrooptic tensor and E_{e} the external electric field), E_{0} is the homogeneous part of the xcomponent of the photorefractive spacecharge field. The temporal evolution of E_{0} is calculated using a relaxationtype dynamics given by: , where τ is the relaxation time of the crystal and . A positive external electrical field is applied along the caxis of the crystal (parallel to the xaxis) [Fig. 1(b)]. Through the Pockels effect the optical Airy beams locally photoinduce a refractive index variation in the photorefractive crystal which leads to optically induced complex waveguiding structures in the photorefractive nonlinear material. The induced refractive index distribution is then related to the combination of the multiplexed focused Airy beams F and B.
Under focusing conditions, the Airy beam undergoes nonlinear interactions inside the photorefractive medium and most of the beam turns into an “offshooting soliton”, while a small fraction of the power remains a selfaccelerating linear packet [Fig. 1(b)]^{18,30}. As shown recently^{37}, the refractive index variation structure photoinduced by a single or two CP Airy beams enables to guide optical beams along the crystal similar to systems of two CP conventional beams^{14}. As interactions of two CP Gaussian beams or two CP vortices lead to similar spatiotemporal dynamics, we question whether selfaccelerating Airy beams undergo a similar spatiotemporal behaviour and how the Airy properties influence the dynamics. To describe these dynamics, we analyse the spatiotemporal evolution of the output position of the forward offshooting soliton at the crystal’s output plane z = L [Fig. 1(c)].
Onset of spatiotemporal dynamics
The interaction schemes of two CP beams depend on two main control parameters, that are the nonlinear coupling constant Γ and the crystal length L. Figure 2 shows a stability diagram in the plane of the following parameters: (Γ, L). It depicts the various spatiotemporal dynamics of the forwardpropagating “offshooting” soliton’s intensity at the output I(x, z = L) for two different intensities [Fig. 2(a)] and [Fig. 2(b)]. For an Airy lobe’s waist x_{0} = 10μm (a = 0.01), the parameter range corresponds to 1cm ≤ L ≤ 10cm and, for Γ, an external electrical bias field of a few kV/cm.
For low Γvalues (Γ = 3), the nonlinearity Γ applied on the system is not high enough to create locally a large refractive index variation inside the crystal by the photorefractive effect and therefore to induce an offshooting soliton. Still, the propagation of each Airy beam optically induces a curved waveguide along the deflecting Airy trajectory^{28}. We call this region ‘static waveguide without offshooting soliton’. For a larger nonlinearity strength, each CP Airy beam undergoes selftrapping and a part of the beam’s energy turns into an “offshooting” soliton [Fig. 1(b)]. We define the existence of an offshooting soliton, when at least of the input intensity exits at z = L and can be clearly distinguished from the linear output beam. Since almost half of the energy is stored in the first Airy lobe^{30}, the nonlinearity of the system mostly influences the main lobes and the offshooting solitons. The interaction of the two CP Airy beams then leads to various new static waveguide structures and we call this region ‘static waveguide with offshooting soliton’. As presented in reference^{37}, the photoinduced waveguide structure enables a Gaussian beam to exit the crystal at a single or at two output positions simultaneously. The parabolic trajectory of the CP Airy beams enables waveguiding structures even for transverse shifts of the interacting beams that by far exceed the beam waist. When we still increase the nonlinearity Γ, the waveguide is no longer steady in time but rather shows stable timeperiodic dynamics: the offshooting soliton evolves from a constant transverse output position to an output position that oscillates harmonically in time along the xaxis [Figs 1(c) and 3(c)]. We call this region ‘harmonic oscillations’. Similar to the case of CP Gaussian beams^{8}, the critical nonlinearity strength that delimits the onset of timeperiodic oscillations of the waveguide decreases with the increase of the crystal length L, see the line labelled ‘threshold staticdynamic’ in Fig. 2.
For an even larger Γ and/or crystal length L, the timeperiodic waveguide dynamics is replaced by chaoticlike spatiotemporal dynamics. The position of the offshooting soliton does not vary periodically in time but rather in an erratic way. As will be shown later, while the trajectory is erratic in time, the motion of the offshooting soliton is attracted towards the input positions of the lobes of the counterpropagating Airy beam. We call this parameter region ‘chaotic waveguide’. The critical nonlinearity leading to unstable waveguiding decreases with the increase of the crystal length L, as is also true for CP Gaussian^{8} and vortex beams (see the line labelled ‘threshold dynamicunstable’ in Fig. 2).
Interestingly, we identify two additional regions. In both cases (a) and (b) the timeperiodic dynamically varying waveguide may restabilise to a static waveguide when increasing the nonlinearity. The offshooting soliton stabilises again at a constant output position. The possibility to stabilise again the photoinduced waveguiding by increasing the nonlinearity strength has not been observed earlier with CP Gaussian beams and is related to the multilobe shape of the Airy beams. Therefore this suggests an advantage in using CP Airy beams. We also identify another parameter region where the position of the offshooting soliton varies periodically in time but not in an harmonic way. We have simply called this region ‘timeperiodic waveguide’. This specific dynamics bifurcates from the harmonic waveguide case but is also observed as a bifurcation of the chaotic waveguide case. We shall detail these dynamics and their bifurcations in the next section.
Similarly to other CP beams’ systems, the intensity of the input beams is an important parameter. When increasing the total optical intensity injected in the crystal through the CP beams the refractive index variations increases, hence resulting in more nonlinear interactions; see [Fig. 2(b)]. When we compare the Fig. 2(a,b), the critical nonlinearity that leads to either a timeperiodic waveguide or even chaotic waveguide for a normalized intensity [Fig. 2(a)] is larger than for [Fig. 2(b)]. The stability of the waveguide is therefore reduced by the increase of the optical intensity.
Finally it is worth comparing the critical nonlinearity that leads to dynamically varying waveguide (our dashed line) in the case of CP Airy beams with the one computed for CP Gaussian beams (dotted line). Besides the fact that Airyinduced waveguides have more complex features than Gaussianinduced waveguides, it appears also that, the Airyinduced waveguides are stable in a large range of parameters and in particular for a large range of nonlinearity strength and/or crystal length. This unique property of Airyinduced waveguides is related to the diffractionfree propagation and multilobe shape of Airy beams.
Detailed route to chaos
In this section we analyse the nature and the evolution of the spatiotemporal dynamics of two CP Airy beams for a fixed crystal length L when the nonlinear coupling strength Γ is increased. Physically the nonlinearity is increased through the positive electrical bias field applied on the crystal. Although as mentioned earlier the stability of the photoinduced waveguide depends on both the crystal length L and the beam intensities, we shall restrict ourselves to one case where and L = 5.5L_{d}. For Airy beams with the parameters x_{0} = 10μm, a = 0.1, it corresponds to a crystal length of L = 28mm (see arrow [Fig. 2(a)]). This case illustrates the complexity underlying the sequence of bifurcations to spatiotemporal instabilities of the waveguide [Fig. 3(a)]. A similar sequence of bifurcations occurs when varying the system parameters. For each value, we simulate the propagation of two CP Airy beams over t_{f} = 100τ_{0}, where at each crystal’s side the main lobe of the CP Airy beams is centered around x = −x_{0} for its input position. We then display the spatiotemporal dynamics of the forward offshooting soliton at the crystal’s output side z = L along the transverse xaxis by plotting the offshooting transverse position versus time [Fig. 3(b–g)]. To avoid the transient dynamics, we detect the extreme positions taken by the offshooting soliton within the times t_{1} = 20τ_{0} and t_{f} = 100τ_{0}. The bifurcation diagram on Fig. 3(a) resumes the position of the spatial output of the offshooting soliton during time: for each value, the various dots display the xextrema taken by the offshooting soliton along time.
The diagram on Fig. 3(a) displays the route to instabilities from a system with a weak nonlinearity (Γ = 9) to a highly nonlinear system (Γ > 16). For , the bifurcation diagram displays a steadystate transverse output position of the offshooting soliton during time [Fig. 3(b)]. The steadystate case depicted in Fig. 3(b) corresponds to the waveguide structures demonstrated in^{37}, where the CP Airy beams and their offshooting solitons coexist in the crystal. When , two extrema of the x position of the offshooting soliton appear for a given Γ value. The timetrace of the position of the offshooting soliton displays a sinusoidal evolution [Fig. 3(c)]. We observe a stable oscillating dynamics, where the offshooting soliton rotates periodically around its characteristic position x = 2.5x_{0}. The period of the sinusoidal oscillation is about 3.5τ_{0}, i.e. is of the same order of magnitude than the material nonlinear optics timescale. The amplitude of the oscillation is determined by the Airy properties and in particular their deflection characteristics. Indeed the amplitude of the oscillation is larger for a longer crystal since by increasing the crystal length, the CP Airy beams deflect more before colliding. The Fig. 3(d) shows that this oscillating soliton dynamics restabilises when increasing the nonlinearity leading to a new static waveguide structure. This singular case will be explained in details and illustrated in Fig. 4(b) of the next section.
When increasing the nonlinearity strength above Γ = 13, the position of the offshooting soliton presents an erratic motion along the output plane [Fig. 3(a)]. Figure 3(e) indicates that the soliton tends to follow alternatively a complex timeperiodic, then a chaoticlike behaviour. As depicted on Fig. 3(f), when we increase the nonlinearity (Γ = 14.8), the chaoticlike evolution of the position of the offshooting solution stabilises to a timeperiodic dynamics where oscillations at a slower timescale modulate the dynamics with a higher amplitude than for the harmonic oscillation depicted in Fig. 3(c).
Finally when Γ > 15.3, the timeperiodic dynamically varying waveguide becomes unstable and the position of the offshooting soliton rotates in an erratic way around the single Airy case position x = 0 and the CP main lobe’s position x = −x_{0} [Fig. 3(g)]. It is worth noting that for a very high coupling strength (Γ > 16.3) the erratic motion of the offshooting soliton encompasses additional attractive xvalues at . Interestingly these xvalues correspond to the respective input positions of the second, third and fourth lobe orders of the CP Airy beam [Fig. 3(a)]. The characteristics of this chaotic soliton motion will be further discussed in the next sections.
Restabilisation of the waveguide for large nonlinearity strength
As previously emphasized in Fig. 2, dynamically varying waveguides photoinduced by two counterpropagating Airy beams can restabilise when the nonlinearity Γ increases. So far, the dynamical behaviour of two CP Gaussian beams has only shown an evolution from a steadystate, then timeperiodic to chaotic like regimes. The Fig. 4 compares the typical static waveguides with offshooting soliton that can be observed in the two stability diagrams [Fig. 2(a–b)]. Figure 4(a) corresponds to the waveguide structures presented in^{37} where a Gaussian probe beam can be guided along the crystal to one or two outputs. Figure 4(b) illustrates the new waveguide structure in the restabilization zone (‘static waveguide with off shooting soliton’ above the line labelled ‘threshold staticdynamic’ [Fig. 2]). This waveguide structure offers the same type of photoinduced waveguides as in Fig. 4(a) but with a better coupling efficiency (up to instead of ) in the offshooting solitons due to the stronger focusing nonlinearity. Figure 4(c) also presents a particular steadystate structure in the strong intensity case at two parameter points: (L = 5, Γ = 7.5) and (L = 5.5, Γ = 7). Contrary to the usual steadystate case, where the CP Airy beam induces a transverse shift of the offshooting soliton^{37}, the offshooting soliton changes its output position and merges exactly into the main lobe of the CP Airy beam. In the configurations (b) and (c) where the nonlinearity of the system is increased through or the initial intensity, the spacecharge field photoinduced by the multiple lobe orders of the two CP Airy beams has a significant role in the interaction schemes in the photorefractive crystal. Although the main power is transferred into the offshooting soliton during the nonlinear propagation of the Airy beams^{30}, the secondary lobes of the Airy beams are essential for the restabilization of our system above the conventional steadystate threshold curve.
Chaotic motion of Airyinduced soliton
As previously emphasized, the dynamical behaviour of the photoinduced waveguide significantly depends on the crystal length. In the ‘chaotic waveguide’ region [Fig. 2], we propose therefore to compare the situation of a short crystal (e.g. L = 2.5L_{d} = 13mm) and the one of a long crystal (e.g. L = 5.5Ld = 28mm). It is worth mentioning that due to their parabolic trajectory, the two CP Airy beams intersect at for the short crystal and at for the long crystal. As a result, the photoinduced waveguides originate mostly from the interaction of the first and second lobe orders of the CP Airy beams in the case of a short crystal and from the interaction of the four first lobe orders of the CP Airy beam in the case of a longer crystal. The resulting waveguiding structure in the case of a longer crystal will be larger along the transverse axis ( instead of in the case of a short crystal). Similarly, in the case of a long crystal, the transverse trajectory of the offshooting solitons will shift from its typical transverse position x_{soliton} = 0 towards the + xdirection.
In particular, as illustrated in Fig. 3(a,g) for a long crystal, under high nonlinear conditions, the system of two CP Airy beams shows a peculiar instability pattern: the output positions of the offshooting soliton in the unstable regime appear to be attracted toward very specific output positions, which correspond to the respective input positions of the different lobe orders of the CP Airy beam. Figure 5 depicts the statistical distribution of the output position of the offshooting soliton in (a) a short crystal (L = 2.5L_{d}) and (b) in a long crystal (L = 5.5L_{d}) with . The output positions of the offshooting soliton are not distributed in a continuous way but rather in a discrete way. The privileged output positions match with the input positions of the Airy lobes of the CP backward beam (blue zone). Also the attraction strength, measured by the highest probability in the plotted histograms, decreases for the higher lobe orders, as the spacecharge field related to the energy distribution of the Airy beam decreases along the −xaxis. In the short crystal case (a), the offshooting soliton is also attracted towards the + xaxis, at the output positions of the Airy lobes of the forward beam (green zone).
Such spatially localized instabilities have never been observed in an optical system using CP beams^{8,11}. Our system made of CP Airy beams therefore creates a chaotic motion of the offshooting soliton whose topology can be engineered by both the Airy beam properties and the photorefractive crystal nonlinearity strength and length. Recent years have seen a tremendous interest in applications of optical chaos for alloptical signal processing including optical generation of random numbers. The most conclusive proposals so far have used the temporal chaotic output of semiconductor lasers^{38}. The digital sampling of optical chaos allows to extract random bits at high bit rate^{39}. The extension to massive parallel computing is however limited in that it requires either a large number of such chaotic lasers or the use of uncorrelated emission from individual laser longitudinal or transverse modes. In the present scheme, one has access to a chaotic output (the erratic motion of the offshooting soliton) that is by essence spatially multiplexed at discrete positions that match the locations of Airy beam lobes. Our findings therefore suggest innovative ways of performing multiplexed chaosbased optical computing.
Discussion
To conclude, the interaction of two CP Airy beams in a photorefractive crystal leads to peculiar spatiotemporal dynamics. The system evolves from static to timeperiodic then chaotic waveguides when increasing the nonlinearity strength and the crystal length. We demonstrated the existence of a threshold curve above which non steadystate dynamics appear. By comparison to similar studies using CP Gaussian beams, photoinduced Airy waveguides are stable for a larger range of parameters. Also on the route to instabilities, we identify a singular additional region where dynamical waveguides restabilise to static waveguides with a better coupling efficiency. When the system bifurcates to the chaoticlike dynamics, the offshooting soliton moves in an erratic way with privileged positions that match the input positions of the multiple lobes of the CP Airy beams. Such spatially localized instabilities suggest innovative ways of performing optical computing based on spatiotemporal chaos. The unique properties of static and dynamic Airy waveguide structures motivate experimental demonstration and implementation in different nonlinear optical media.
Additional Information
How to cite this article: Wiersma, N. et al. Spatiotemporal dynamics of counterpropagating Airy beams. Sci. Rep. 5, 13463; doi: 10.1038/srep13463 (2015).
References
Silberberg, Y. & Joseph, I. Instabilities, selfoscillation and chaos in a simple nonlinear optical interaction. Phys. Rev. Lett. 43, 1541–1543 (1982).
Firth, W. & Paré, C. Transverse modulational instabilities for counterpropagating beams in Kerr media. Opt. Lett. 13, 1096–1098 (1988).
Gauthier, D., Malcuit, M., Gaeta, A. & Boyd, R. Polarization bistability of counterpropagating laser beams. Phys. Rev. Lett. 64, 1721–1724 (1990).
Gaeta, A., Boyd, R., Ackerhalt, J. & Milonni, P. Instabilities and chaos in the polarizations of counterpropagating light fields. Phys. Rev. Lett. 58, 2432–2435 (1987).
Henninot, J., Blach, J. & Warenghem, M. Experimental study of the nonlocality of spatial optical solitons excited in nematic liquid crystal. J. Opt. A Pure App. Opt. 9, 20 (2007).
Chang, R. & Meystre, P. Instabilities of counterpropagating beams in a nonlinear medium . Phys. Rev. A 44, 3188–3200 (1991).
Denz, C. et al. Dynamics of formation and interaction of photorefractive screening solitons. Phys. Rev. E 60, 6222–6225 (1999).
Belić, M. et al. Counterpropagating selftrapped beams in photorefractive crystals. J. Opt. B Quantum Semiclassical Opt. 6, S190 (2004).
Jander, P. et al. Spatiotemporal dynamics of counterpropagating photorefractive selftrapped beams. In Photorefractive Effects, Materials and Devices. 522 (Optical Society of America, 2005).
Chen, Z., Segev, M. & Christodoulides, D. N. Optical spatial solitons: historical overview and recent advances. Rep. Prog. Phys. 75, 086401 (2012).
Petrović, M., Belić, M., Denz, C. & Kivshar, Y. Counterpropagating optical beams and solitons. Laser Photon. Rev. 5, 214–233 (2011).
Haelterman, M., Sheppard, A. & Snyder, A. Bimodal counterpropagating spatial solitarywaves. Opt. Comm. 103, 145–152 (1993).
Wolfersberger, D. & Tranca, D. 2D infrared selffocusing in bulk photorefractive SBN. Opt. Mater. Express 1, 1178–1184 (2011).
Petter, J. & Denz, C. Guiding and dividing waves with photorefractive solitons. Opt. Comm. 188, 55–61 (2001).
Peccianti, M., Conti, C., Assanto, G., De Luca, A. & Umeton, C. Alloptical switching and logic gating with spatial solitons in liquid crystals. Appl. Phys. Lett. 81, 3335–3337 (2002).
Rotschild, C., Alfassi, B., Cohen, O. & Segev, M. Longrange interactions between optical solitons. Nature Phys. 2, 769–774 (2006).
Belić, M., Jović, D., Prvanović, S., Arsenović, D. & Petrović, M. Counterpropagating selftrapped beams in optical photonic lattices. Opt. Express 14, 794–799 (2006).
Kaminer, I., Bekenstein, R., Nemirovsky, J. & Segev, M. Nondiffracting accelerating wave packets of Maxwell’s equations. Phys. Rev. Lett. 108, 163901 (2012).
Hu, Y. et al. Reshaping the trajectory and spectrum of nonlinear Airy beams. Opt. Lett. 37, 3201–3203 (2012).
Panagiotopoulos, P. et al. Nonlinear propagation dynamics of finiteenergy Airy beams. Phys. Rev. A 86, 013842 (2012).
Berry, M. & Balazs, N. Nonspreading wave packets. Am. J. Phys 47, 264–267 (1979).
Siviloglou, G. A. & Christodoulides, D. N. Accelerating finite energy Airy beams. Opt. Lett. 32, 979–981 (2007).
Ellenbogen, T., VolochBloch, N., GananyPadowicz, A. & Arie, A. Nonlinear generation and manipulation of Airy beams. Nature Photon. 3, 395–398 (2009).
Hu, Y. et al. Optimal control of the ballistic motion of Airy beams. Opt. Lett. 35, 2260–2262 (2010).
Dolev, I., Ellenbogen, T. & Arie, A. Switching the acceleration direction of Airy beams by a nonlinear optical process . Opt. Lett. 35, 1581–1583 (2010).
Baumgartl, J. et al. Optical path clearing and enhanced transmission through colloidal suspensions. Opt. Express 18, 17130–17140 (2010).
Fahrbach, F. O., Simon, P. & Rohrbach, A. Microscopy with selfreconstructing beams.Nature Photon. 4, 780–785 (2010).
Rose, P., Diebel, F., Boguslawski, M. & Denz, C. Airy beam induced optical routing. App. Phys. Lett. 102, 101101 (2013).
Minovich, A. E. et al. Airy plasmons: nondiffracting optical surface waves. Laser Photon. Rev. 8, 221–232 (2014).
Allayarov, I. & Tsoy, E. Dynamics of Airy beams in nonlinear media. Phys. Rev. A 90, 023852 (2014).
Zhang, Y. et al. Soliton pair generation in the interactions of Airy and nonlinear accelerating beams. Opt. Lett. 38, 4585–4588 (2013).
Shen, M., Gao, J. & Ge, L. Solitons shedding from Airy beams and bound states of breathing Airy solitons in nonlocal nonlinear media. Sci. Rep. 5 (2015).
Zhang, Y. et al. Interactions of Airy beams, nonlinear accelerating beams and induced solitons in Kerr and saturable nonlinear media. Opt. Express 22, 7160–7171 (2014).
Driben, R., Konotop, V. V. & Meier, T. Coupled Airy breathers. Opt. Lett. 39, 5523–5526 (2014).
Fattal, Y., Rudnick, A. & Marom, D. M. Soliton shedding from Airy pulses in Kerr media. Opt. Express 19, 17298–17307 (2011).
Rudnick, A. & Marom, D. M. Airysoliton interactions in Kerr media. Opt. Express 19, 25570–25582 (2011).
Wiersma, N., Marsal, N., Sciamanna, M. & Wolfersberger, D. Alloptical interconnects using Airy beams. Opt. Lett. 39, 5997–6000 (2014).
Sciamanna, M. & Shore, K. A. Physics and applications of laser diode chaos. Nature Photon. 9, 151–162 (2015).
Virte, M., Mercier, E., Thienpont, H., Panajotov, K. & Sciamanna, M. Physical random bit generation from chaotic solitary laser diode. Opt. Express 22, 17271–17280 (2014).
Acknowledgements
The authors acknowledge the support of Conseil Régional de Lorraine, Fondation Supélec, the IAP P7/35 (BELSPO) with the ‘Photonics@be’ project (2012–2017) and Fonds européen de développement régional (FEDER) with the project ‘PHOTON’ (2015).
Author information
Affiliations
Contributions
N.W. performed the numerical simulations. All the authors contributed to the data analysis and to the writing of the manuscript.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Wiersma, N., Marsal, N., Sciamanna, M. et al. Spatiotemporal dynamics of counterpropagating Airy beams. Sci Rep 5, 13463 (2015). https://doi.org/10.1038/srep13463
Received:
Accepted:
Published:
Further reading

Suppression of collapse for twodimensional Airy beam in nonlocal nonlinear media
Scientific Reports (2017)

Experimental study on the propagation characteristics of ring Airy Gaussian vortex beams
Applied Physics B (2017)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.