All-optical interconnects using Airy beams

We analyze numerically optical waveguide structures in photorefractive media induced by one or two incoherent counter-propagating (CP) Airy beams. Under nonlinear focusing conditions, we show that for a single Airy beam or for two CP beams with various input positions, multiple waveguiding structures can be photo-induced in the medium. Optical Gaussian beams can therefore be guided with a deflecting trajectory and/or even split into several output beams. These results enable new configurations for all-optical interconnections. © 2014 Optical Society of America OCIS codes: (190.4420) Nonlinear optics, transverse effects in; (190.6135) Spatial solitons. http://dx.doi.org/10.1364/OL.39.005997

The Airy wave function has been initially suggested by Berry and Balazs [1] as a solution for the Schrödinger equation with a non-spreading wave propagation along a parabolic trajectory. First experimentally observed in free space in 2007 [2] by means of a spatial light modulator (SLM), Airy beams have been also generated using photonic crystals [3]. Later, a main interest has been to control their trajectories [4,5] and to maintain their Airy shape [6] in nonlinear media. Their unique properties make Airy beams good candidates for micro manipulation of particles [7], imaging in scattering media [8], and waveguiding for optical routing communications [9].
Recently, the self trapping character of Airy beams in biased nonlinear media has suggested interesting dynamics such as soliton-like behaviors and interactions of co-propagating Airy functions [10,11]. By contrast to the co-propagating case, a counter-propagating (CP) configuration combined with the large transverse dimension of the Airy shape and their curved trajectory allow several interaction schemes where one or several lobes of the CP beams can, for example, face each other. The interaction of CP Airy beams opens therefore new interesting fields for optical interconnections.
In this Letter, we analyze numerically optical waveguide structures created in photorefractive (PR) media by one or two incoherent CP Airy beams under nonlinear self-focusing conditions. The scheme of incoherent CP beams allows strong interactions because of crosscoupling of the beams via the combined refractive index structure that is not present in the co-propagating case [12]. When a positive external electrical field is applied on the PR crystal, we show that for a single Airy beam or for two CP beams even strongly misaligned, multiple waveguiding structures are photo-induced. In that way, an optical Gaussian beam can be linearly guided along a deflecting trajectory or split into several output beams.
Our typical interaction scheme is depicted on Fig. 1.It illustrates the linear propagation of two 1D CP Airy beams in an unbiased PR crystal. Airy profiles, Fx; z and Bx; z, are given by the following equations: where F 0 and B 0 correspond to the wave amplitudes of respectively the forward beam F and the backward beam B, Ai represents the Airy function, x 0 an arbitrary transverse scale, a the truncation factor, and d the additional normalized shifting distance between the initial input beams F and B. CP beams are herein defined by their opposite propagation directions along the longitudinal z axis [ Fig. 1]. Each Airy beam is formed of successive lobes, the first order being the main lobe. The nonlinear propagation of these two incoherent CP beams can be expressed as follows [13]: where Γ kn 0 x 0 2 r eff E e is the nonlinear photorefractive coupling strength (r eff is the effective component of the electro-optic tensor and E e the external electric field), and E 0 is the homogeneous part of the x-component of the photorefractive space-charge field. The temporal evolution of E 0 is calculated using a relaxation-type dynamics given by τ∂ t E 0 E 0 −I 0 ∕1 I 0 , where τ is the relaxation time of the crystal and I 0 jFj 2 jBj 2 . A positive external electrical field is applied along the c axis of the crystal (parallel to the x axis [ Fig. 1]). This type of interaction enables us to optically induce complex waveguiding structures in the PR nonlinear material through the Pockels effect [14]. The induced refractive index distribution is then related to the combination of the multiplexed focused Airy beams F and B.
The induced refractive index profile is numerically simulated via a fast Fourier transform beam propagation method. We fix the crystal length L 1 cm, Γ 9, a 0.09, x 0 7.5 μm, the normalized input energy of a beam F 0 2.5 p , and keep B 0 and d as free parameters. We first consider optical waveguides created by a single Airy beam F (B 0 0) propagating from z 0 to z L (z 0 → L) with a focusing bias electric field. Figure 2 displays the intensity of the self-focused Airy beam. The results concern the stationary state that is reached after a transient duration equal to 6τ.Atz 0, we indicate only two positions of high index variations corresponding to the two first lobe orders of the input Airy function (1,2 in Fig. 2). Interestingly and similar to the references [10,15], the initial Airy beam with its curved trajectory turns into an "off-shooting" soliton as it propagates along the z axis (1 0 at z L in Fig. 2). This simple configuration created by a single Airy beam under high-focusing conditions raises the question whether the photo-induced structure can be used for various linear guidings of Gaussian beams. To test the waveguide structure in Fig. 2, we now inject at the different input positions 1; 2; 1 0 a Gaussian probe beam with its size matching with the photo-induced waveguide at each input position 1; 2; 1 0 . We consider the linear propagation by simulating Eq. (3) or Eq. (4) and keeping the space-charge field E 0 equal to the stationary E 0 x; z resulting from the nonlinear propagation of the Airy beam. The crystal length L 1 cm corresponds to 3.3L d , with L d 2kx 2 0 the characteristic diffraction length of the Gaussian beam (n 0 2.3, λ 532 nm). The results show two effective waveguides with different output structures: The first one [Figs. 3(a1)-3(b1)] guides 66% of the initial energy along the crystal [red line, Fig. 3(b1)]. Based on Fig. 2(a), we stress that the waveguide 1 → 1 0 is first formed by the self-focused main Airy lobe, then by the "off-shooting" soliton which no longer satisfies the Airy beam's properties but rather the ones of a self-focused Gaussian beam [15]. This guiding differs from the situation analyzed in [9] in that, under high-focusing conditions, the output position of the waveguide 1 → 1 0 does not relate to the deflection of the Airy beam but is rather determined by the position of the "off-shooting" soliton. The scenario is different when the Gaussian probe beam is injected in 1 0 and propagates in the other direction (z L → 0) [Figs. 3(a2)-3(b2)]. We show that the energy of the initial Gaussian beam splits now into two outputs: 30% in 1 and 20% in 2. The demultiplex structure (1 0 → 1; 2) allows for a large spacing between the two outputs (1 and 2): ≈3x 0 [ Fig. 3(b2)]. It is worth emphasizing that such a waveguide with one input and two outputs cannot be achieved with a single Gaussian beam and therefore utilizes the unique properties of the Airy beam. The waveguide 2 → 1 0 does not present strong guiding of the Gaussian probe beam along the crystal [Figs. 3(a3)-3(b3)]. This low guiding efficiency can be explained by the lower refractive index structure photo-induced by the second lobe order of the Airy beam at z 0 and by the asymmetrical waveguide structure between z 0 and z L. Those results demonstrate that, although the waveguiding is bidirectional (i.e., a Gaussian probe beam on both crystal sides is guided through the structure), the asymmetry of the waveguide leads to different output positions and guiding efficiencies depending on the input position of the probe beam.
To enlarge the possible interconnection schemes, we now study the interactions of two self-focused CP Airy beams in the PR crystal. We restrict our investigation to two configurations: the first one where the two CP Airy beams (F and B with Results for F and B initially aligned (d 0): First, it is worth noting that the backward beam modifies the trajectory of the forward beam F [Figs. 4(a1)-4(b1)]. Indeed, at z L, the position of the "off-shooting" soliton 1 0 (respectively 3 0 at z 0) is not the same compared to the situation with only one Airy beam (Fig. 2). We notice an additional drift of 2x 0 . This configuration enables therefore larger transverse distances for optical guiding   (b1)]. The strong interaction due to cross-coupling of the beams via the combined refractive index structure and energy flows between not only the main CP lobes but also the neighbors [ Fig. 4(a1)], [4] can explain the possible energy losses on the depicted structures. As a main conclusion we can stress that the use of two CP Airy beams allows for achieving complex waveguiding structures that would otherwise require the CP interactions of more than two Gaussian beams.
Results for F and B initially misaligned In this configuration, higher energy is exchanged within the CP main and secondary lobes. If a probe beam is injected into the PR crystal, different output cases are observed [ Fig. 6]. Similar to the structure generated when d 0, if the probe beam is injected into inputs formed by "off-shooting" solitons (1 0 or 3 0 ), it will split into two beams at the Airy lobes positions on the other side of the crystal with input-output shifts up to 6x 0 (3 0 → 3; 4 [Figs. 6(a2)-6(b2)]). However, if the Gaussian beam is first guided by a main Airy lobe (1 or 3), we not only observe, as in Figs. 5(a1)-5(b1), an output beam at the soliton position (1 0 or 3 0 ), but also a second one appears at the position of the main lobe of the CP beam (3 or 1) [Figs. 6(a1)-6(b1)]. Such a configuration where interlobe interaction areas inside the crystal are possible, enables creating 1-to-2-demultiplexers guiding up to half the energy of the input beam (outputs: 33%-15%) with an important transverse shift x. The resulting waveguide structures show interesting features. First the waveguide occurs although the positions of the two input beams are largely shifted (2x 0 ). The same simulation, but using shifted Gaussian beams at the entrances 1 and 3, does not create any waveguide. Second the interaction scheme allows for controlling the spatial positions of the waveguide outputs by varying the transverse shift d between the CP input beams.
The high-focusing conditions of the propagation medium has led to a soliton-like behavior of the Airy beam, and our interaction studies have shown first coupling effects between the CP Airy beams resulting in different waveguiding possibilities. It is worth stressing that all previous cases correspond to stationary situations in that the waveguide remains fixed when time increases above a transient duration of the order of the relaxation time of the crystal [ Fig. 7(a)]. Although not discussed in detail, we do find parameters where the waveguiding is not stationary anymore. In Figs. 7(b) and 7(c), when increasing the nonlinear parameter Γ, we observe a situation where the intensity profile at the crystal output varies periodically [ Fig. 7(b)] or irregularly [ Fig. 7(c)] over time.
To conclude, we have demonstrated waveguiding achieved by the nonlinear interactions of CP Airy beams.   In comparison with waveguide structures created by Gaussian beams with the same dimensions and focusing conditions [12], the interaction schemes between two Airy beams offer larger waveguide possibilities thanks to their unique Airy shape and properties: (1) A single Airy beam photo-induces a waveguide with two outputs; (2) two CP Airy beams induce multiple output waveguiding structures that cannot be achieved with only two CP Gaussian beams; and (3) the resulting waveguiding remains possible even for transverse shifts of the interacting beams that by far exceed the beam waist. The numerical simulations presented in this Letter motivate an experimental implementation of CP Airy beams induced optical routing. In addition the first spatiotemporal dynamics observed in a CP Airy beam system encourage further investigation into the relations between the Airy properties and the solitonic interactions.