• Photonics Research
  • Vol. 7, Issue 7, 754 (2019)
Dae Seok Han and Myeong Soo Kang*
Author Affiliations
  • Department of Physics, Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, South Korea
  • show less
    DOI: 10.1364/PRJ.7.000754 Cite this Article Set citation alerts
    Dae Seok Han, Myeong Soo Kang. Reconfigurable generation of optical vortices based on forward stimulated intermodal Brillouin scattering in subwavelength-hole photonic waveguides[J]. Photonics Research, 2019, 7(7): 754 Copy Citation Text show less

    Abstract

    An all-optically reconfigurable generation of optical vortices would be highly beneficial to the implementation of next-generation optical communication and advanced information processing. The previously demonstrated approaches based on the parametric nonlinear optical processes, however, have exhibited limited conversion efficiency due to the group velocity mismatch and nonlinear phase shifts, and require the cumbersome preparation of either the optical element or initial seed beam having a non-zero topological charge. Here, we propose and analyze a novel scheme for highly efficient all-optical generation and control of optical vortices based on the dynamic acoustic vortex grating created by forward stimulated intermodal Brillouin scattering in a subwavelength-hole photonic waveguide. The dual-frequency pump beams in two different hybrid optical modes drive an acoustic vortex mode, which transforms a signal in the fundamental optical mode into an optical vortex mode. This scheme not only eliminates the need for the initial preparation of an angular-momentum-carrying medium or an optical vortex seed but also guarantees high modal purity and nearly 100% conversion efficiency assisted by the energy-momentum conservation. We also investigate the feasibility and practicability of the subwavelength-hole waveguides by examining the intermodal conversion efficiency and robustness of guidance of the optical vortices, taking into account the impact of the Kerr-type nonlinear effects on the intermodal Brillouin interactions based on our rigorous full-vectorial analytical theory.

    1. INTRODUCTION

    Spatially structured light—an optical wave with a customized transverse profile of intensity, polarization, and phase—has been in the limelight as a versatile tool for exploring fundamental science and applications such as micro/nano-fabrication [1], topological statistics [2], astronomy [3], imaging [4], optical tweezers [5], quantum optics [6], and optical communications [7]. In particular, optical vortex beams (OVBs) possessing phase singularity at the beam center and a donut-shaped intensity pattern exhibit extraordinary properties that lead to new types of light–matter interactions, which in some cases enhance the performance of photonic devices and make a technological breakthrough. A notable example is a radical improvement in super-resolution fluorescence microscopy [4], where the conventional resolution limit can be circumvented by adopting the donut-shaped laser beam to cause the fluorescent molecules to be depleted only around the dark hollow beam center and leave a tiny beam spot from the fluorophore. Furthermore, helically phased light beams carry orbital angular momentum (OAM), which has been utilized mainly for exerting a torque precisely on the micro/nano-particles trapped in optical tweezers [5]. The optical OAM also provides intrinsically infinite dimensionality, bringing ongoing attention to high-dimensional quantum cryptography [6]. Mode-division multiplexing based on the OAM modes has also emerged as a scintillating solution to the exponentially growing demand for data capacity in optical communications [7].

    Until very recently, OAM-carrying OVBs have been generated in most cases using specially designed optical media having topological charges, e.g., spiral phase plates [8], liquid crystal Q-plates [9], fiber gratings [10], and helically twisted fibers [11]. In addition, the development of holographic techniques using digital spatial light modulators (SLMs) has allowed for OVB generation based on real-time computer-controlled wavefront shaping [12]. On the other hand, parametric nonlinear optical processes have also attracted ongoing attention as an alternative means to generate and amplify the OVBs, as they enable dynamic all-optical manipulation and offer high damage threshold compared to the SLMs. All-optical creation of the OAM states has been relentlessly investigated using second-harmonic generation [13], intermodal four-wave mixing (FWM) [14], and parametric down-conversion [15], which would be beneficial especially to next-generation ultrahigh-capacity optical communications. However, these schemes have suffered from several disadvantages that hinder efficient OVB generation. The group velocity mismatch between the pump and signal gives rise to their temporal walk-off, which reduces the interaction length. In addition, the undesired nonlinear phase shift induced by the self-phase modulation (SPM) and cross-phase modulation (XPM) also limits the overall efficiency of the OVB generation. A multi-level atomic ensemble with a photonic bandgap structure has also been studied for transferring and controlling the OVBs via FWM and six-wave mixing, which could be modulated by superimposing the nonlinear phase shift in the atomic ensemble [1619]. In this case, the interactions between the atomic ensemble and the electromagnetic fields could enhance the nonlinear optical susceptibility through the induced atomic coherence and the decrease of optical group velocity. However, they are not all-optical in the true sense because an optical seed carrying a topological charge is still required for their operation and therefore should be initially prepared using a delicately structured medium.

    Meanwhile, the interaction between light and acoustic phonons has offered fundamentally different concepts of all-optical signal processing, as it has no all-optical analog among the parametric nonlinear optical processes. For instance, new optical frequencies or spatial modes can be excited efficiently through strong photon–phonon interactions even in the absence of optical resonances, which has been demonstrated in various forms, e.g., cascaded stimulated Raman-like scattering [20], phonon-mediated optical signal transmission [21], and unidirectional interband modulation [22]. The amplification and generation of OAM-carrying light via stimulated Raman scattering and stimulated Brillouin scattering (SBS) has been recently studied in liquid [23], gas [24], and plasmas [25]. On the other hand, in solid-state photonic waveguides such as optical fibers, it has only been demonstrated to date that the angular momentum in the externally applied acoustic wave can be transferred to the guided optical mode via acousto-optic interactions [26]. If the acoustic phonons could be optically generated and controlled in photonic waveguides, they would furnish promising techniques in photonic integrated circuits for all-optical manipulation of OAM light.

    In this paper, we show that nonlinear optoacoustic interactions via forward stimulated intermodal Brillouin scattering (FSIBS) in suitably designed photonic waveguides can be employed for all-optical generation and control of OAM-carrying optical vortex modes (OVMs), which excludes the need for the initial preparation of any vortex seed. As shown in the schematic diagram of Fig. 1(a), when dual-frequency optical pump waves are launched into the waveguide, one frequency (ω1) being in the circularly polarized fundamental mode (i.e., HE11x+iHE11y) and the other (ω2) in the higher-order hybrid mode (e.g., HE21even), their beat produces oscillating optical forces (i.e., electrostrictive force and radiation pressure), which in turn coherently drive two orthogonal acoustic resonance (AR) modes of the frequency Ω=ω1ω2. The phases of the two AR modes differ from each other by 90°, creating an acoustic vortex mode (AVM) [26]. The optical signal of the frequency ω3 in the circularly polarized fundamental mode is then converted into the frequency-downshifted (ω4=ω3Ω) targeted OVM via interaction with the AVM. We note that the intermodal coupling through the optically generated acoustic vortex grating might be related to the conventional backward-SBS-based dynamic grating in polarization-maintaining fibers, where a traveling acoustic wave created by two counterpropagating optical pump waves of one polarization reflects the orthogonally polarized signal [27]. It should be emphasized that in sharp contrast to the previously studied SBS-based mode conversion [28], in our case an optical signal in the circularly polarized fundamental mode is scattered into the associated OVM rather than the higher-order hybrid mode, the same as the higher-order pump mode, due to the total energy-momentum conservation. Hence, our scheme in principle guarantees high-purity OVM generation and can be exploited to switch the OVM state simply by varying the polarization state of the pump wave.

    (a) Schematic diagram of optical vortex (OV) generation based on forward stimulated intermodal Brillouin scattering. A circularly polarized fundamental mode (a1) and a higher-order hybrid mode (a2) drive a coherent acoustic vortex (AV) (b), which in turn transforms the incident signal (a3) in the circularly polarized fundamental mode into the OV beam (a4). (b) Vector diagram describing the angular momentum conservation in this process.

    Figure 1.(a) Schematic diagram of optical vortex (OV) generation based on forward stimulated intermodal Brillouin scattering. A circularly polarized fundamental mode (a1) and a higher-order hybrid mode (a2) drive a coherent acoustic vortex (AV) (b), which in turn transforms the incident signal (a3) in the circularly polarized fundamental mode into the OV beam (a4). (b) Vector diagram describing the angular momentum conservation in this process.

    While such a photon–phonon interaction is generally construed as a nonlinear wave-mixing process, it can also be described as an angular momentum transfer process between photons and phonons [Fig. 1(b)]. For paraxial beams propagating in free space or weakly guiding optical fibers, the spin angular momentum (SAM) and OAM are individually conserved in the photon–phonon interactions [23,25]. The AVM, having no SAM, can then convey to photons with its OAM only, their polarization state being unchanged. On the other hand, in strongly confining optical waveguides (e.g., microfibers, small-solid-core photonic crystal fibers, and on-chip integrated waveguides), the significant electric fields at the waveguide boundaries yield spin-orbit coupling [29]. In this case, the optical modes can no longer be treated as the eigenstates of either the SAM (S) or the OAM (L) operator, and the total angular momentum (AM) J=S+L should be a conserved quantity. In sharp contrast to the OAM transfer processes in the paraxial regime [23,25], it turns out that both the SAM and OAM of the photons can be converted simultaneously in the FSIBS-based OVM generation process. This result would lead to a number of fundamentally intriguing and intricate phenomena in combination with the distinct features of optical AM, e.g., quantum spin Hall effect of light and optical spin-momentum locking [30].

    2. PRINCIPLES

    Here, we consider a silica-glass subwavelength-hole waveguide (SWHW) suspended in the air, as shown in Fig. 2(a), which has been proposed for robust OVM guidance against perturbations (e.g., bends and twists) [31] and can be implemented in the form of photonic crystal fiber having a circular air hole at the core center [32]. We numerically solve the eigenequations for both the optical and acoustic modes of the three-layered circularly symmetric SWHW [33,34]. We obtain the optical intensity profiles and phase distributions for the three lowest-order OVMs and express each of them in terms of the conventional hybrid modes, as summarized in Fig. 2(b). In general, the OVM, OVlmσ, can be represented as a linear combination of the corresponding even and odd hybrid modes with the phase difference of 90°, i.e., where σ (=+, −) and l (=1,2,3,) denote the SAM and OAM quantum numbers, respectively, in the paraxial approximation [26], and m is the radial mode number. For consistency, we designate the circularly polarized fundamental HE11 mode as the OV01± mode, although it is not an OVM but carries only an SAM.

    (a) Silica-glass subwavelength-hole waveguide (SWHW) suspended in the air, together with the Cartesian coordinates. (b) Intensity profiles of the circularly polarized fundamental mode and the three lowest-order optical vortex modes (OVMs), together with the corresponding compositions of the conventional even/odd hybrid modes for λ=1550 nm, Rin=1.2 μm, and Rout=3.0 μm. The light blue arrows represent the electric field distributions for the hybrid modes and the polarization states of the transverse electric fields for the OVMs. The rightmost column corresponds to the phase pattern of the x component of the electric field, arg(Ex).

    Figure 2.(a) Silica-glass subwavelength-hole waveguide (SWHW) suspended in the air, together with the Cartesian coordinates. (b) Intensity profiles of the circularly polarized fundamental mode and the three lowest-order optical vortex modes (OVMs), together with the corresponding compositions of the conventional even/odd hybrid modes for λ=1550  nm, Rin=1.2  μm, and Rout=3.0  μm. The light blue arrows represent the electric field distributions for the hybrid modes and the polarization states of the transverse electric fields for the OVMs. The rightmost column corresponds to the phase pattern of the x component of the electric field, arg(Ex).

    We determine the AVMs that are phase-matched with the optical modes by taking into account the dispersion of several acoustic modes [Fig. 3(a)]. For a given pair of the core radius Rout and hole radius Rin, we find an AVM satisfying the phase-matching condition q=β1β2 for the two pump beams, where q and β1,2 are the axial wavevectors of the AVM and the two pump optical modes, respectively. The AVM, AVjm, carrying a topological charge of j, can be expressed as a linear combination of the two degenerate conventional acoustic modes having an azimuthal mode number of l=j with a 90° phase difference [26]. For example, the AV±1,1 mode consists of the two orthogonal flexural modes F11x and F11y (“C” in Fig. 3) as AV±1,1=F11x±iF11y, while the AV±2,1 mode is comprised of the two axially anti-symmetric torsional-radial modes as AV±2,1=TR21(0°/90°)±iTR21(45°/45°) (“D” in Fig. 3). (See Visualization 1 and Visualization 2 that display the vibrational motions of the AV1,1 and the AV2,1 modes, respectively.) We also consider the axially symmetric radial R01 mode with no azimuthal dependence (l=0), which we designate simply as the AV01 mode (“B” in Fig. 3). We note that we exclude the longitudinal L01 mode (“A” in Fig. 3) in our analysis, as it is not efficiently driven by quasi-transverse optical forces that can yield the forward intermodal coupling.

    (a) Dispersion curves for the acoustic modes with the lowest radial mode number (m=1) in silica-glass SWHWs. The red, purple, and blue curves represent the acoustic modes with the three lowest azimuthal mode numbers, l=0, 1, and 2, respectively. The frequencies and wavevectors of the acoustic modes are in units of VL/Rout and π/Rout, respectively, where VL=5.972 km/s is the longitudinal sound velocity in fused silica. Some examples of the phase-matched acoustic modes are marked with green filled circles. (b) Displacement distributions of the longitudinal L01 (A), radial R01 (B), flexural F11 (C), and torsional-radial TR21 (D) acoustic modes, together with the resulting acoustic vortex modes (AVMs). The solid and dashed curves represent the deformed and undeformed boundaries, respectively, and the color map is used to describe the profile of the magnitude of total displacement. All the calculations in (a) and (b) are carried out at the wavelength of λ=1550 nm and for the SWHW parameters of Rin=1.2 μm and Rout=3.0 μm. See Visualization 1 and Visualization 2, which display the vibrational motions of the AV1,1 and AV2,1 modes, respectively.

    Figure 3.(a) Dispersion curves for the acoustic modes with the lowest radial mode number (m=1) in silica-glass SWHWs. The red, purple, and blue curves represent the acoustic modes with the three lowest azimuthal mode numbers, l=0, 1, and 2, respectively. The frequencies and wavevectors of the acoustic modes are in units of VL/Rout and π/Rout, respectively, where VL=5.972  km/s is the longitudinal sound velocity in fused silica. Some examples of the phase-matched acoustic modes are marked with green filled circles. (b) Displacement distributions of the longitudinal L01 (A), radial R01 (B), flexural F11 (C), and torsional-radial TR21 (D) acoustic modes, together with the resulting acoustic vortex modes (AVMs). The solid and dashed curves represent the deformed and undeformed boundaries, respectively, and the color map is used to describe the profile of the magnitude of total displacement. All the calculations in (a) and (b) are carried out at the wavelength of λ=1550  nm and for the SWHW parameters of Rin=1.2  μm and Rout=3.0  μm. See Visualization 1 and Visualization 2, which display the vibrational motions of the AV1,1 and AV2,1 modes, respectively.

    We model the propagation of the optical and acoustic modes as well as their interactions using the coupled-mode theory, under the viewpoint that each of the two pairs of optical modes—the two pump waves (a1 and a2) and the incident/coupled signals (a3 and a4) is coupled via the FSIBS-induced acoustic vortex grating (b) (Fig. 1). We also take into account the possibility that these interactions can be significantly influenced by the optical Kerr effect, as in the previously reported phase-locked Brillouin frequency comb generation in a nonlinear waveguide cavity that was affected by the Kerr-induced FWM [35]. Even in the absence of an optical cavity, the SPM and XPM can also make a discernible impact on the dynamics of the SBS-based wave-mixing process [36,37]. Furthermore, the nondegenerate intermodal FWM of the two pump waves and the incident signal may not only affect the targeted OVM signal (a4) but also create an unwanted idler wave (a5) in another OVM that is degenerate with the converted signal (ω5=ω4=ω3Ω) but has the opposite total AM quantum number, j5=j4. In our case, this idler wave gives rise to the detriment of the modal purity of the OVM signal, and therefore it is necessary to examine how significantly the idler wave affects the entire process. Taking all these into consideration, we derive the following steady-state (/t=0) coupled-mode equations for the slowly varying amplitudes of the ν-th (ν=1, 2, 3, 4, 5) optical mode, aν(z), and the AVM, b(z): where ων and Ω represent the frequencies of the ν-th optical mode and the AVM, respectively, αb is the acoustic decay rate, Q and R are the optoacoustic overlap integrals [38,39], and γνSPM, γμνXPM, and γνFWM are the nonlinear coefficients associated with the SPM, XPM, and Kerr-induced FWM, respectively [40]. We also introduce the wavevector mismatch κ=q(β3β4), where β3,4 are the axial wavevectors of the incident and coupled signal modes, while keeping the frequency condition ω3ω4=Ω.

    3. RESULTS AND DISCUSSION

    In our scheme, while the pump wave a1 is kept in the OV01+ (=HE11x+iHE11y) mode, the other a2 can be in any higher-order hybrid mode, which determines the generated OVM according to the total AM conservation. For example, the FSIBS process with the OV01+ and the HE21even modes creates the AVM having the total AM jb=j1j2=1, i.e., the AV1,1 mode, where j1=1 and j2=2 are the total AM of the pump beam a1 and one of the two constituent OVMs of the pump beam a2, the OV11+ mode, respectively. The other constituent OVM of a2, the OV1,1 mode, does not participate in the optoacoustic interaction due to the energy-momentum conservation. The AV1,1 mode converts the signal beam in the OV01+ mode with j3=1 into the OV11+ mode with j4=j3jb=2.

    Here, we focus on the AM transfer processes for the three lowest higher-order hybrid pump modes, the HE21, EH11, and HE31 modes, which we denote as PROC1, PROC2, and PROC3, respectively. To figure out which type of AVM can mediate each AM transfer process efficiently, we obtain and analyze the optical force distribution generated by the two pump beams, as illustrated in Fig. 4. In all the three cases, the total optical force is a linear combination of the two orthogonal force components with the phase difference of 90° and thus can drive the two orthogonal acoustic modes, which composes an AVM. In PROC1, since the HE11x and HE21even modes have the same (opposite) symmetry with respect to the y=0 (x=0) plane, the resulting optical force is symmetric (anti-symmetric) with respect to the y=0 (x=0) plane [Figs. 4(a) and 4(b)]. Then, the F11x mode having the same symmetry with the azimuthal mode number l=1 can be excited. Similarly, the optical force induced by the HE11y and HE21even modes can drive the orthogonal F11y mode, where the minus sign indicates that the optical force is anti-aligned with the displacement of the F11y mode [Figs. 4(c) and 4(d)]. This reversal of the optical force direction can be seen from the fact that by rotating the xy coordinates by 90°, one can find the phase difference of 180° between the optical force created by the HE11y and HE21even modes and that by the HE11x and HE21even modes. Since the two orthogonal flexural modes have a phase difference of 90°, the AV1,1 mode can be excited. A similar explanation holds for PROC3. In PROC2, on the other hand, the HE11x (HE11y) and EH11even modes have the same (opposite) symmetries with respect to both the x=0 and y=0 planes. We note that in Fig. 4, we choose a particular global phase in such a way that the longitudinal component of the total electric field in the waveguide is purely real. Then, only the real part of the optical force spatially matches the profile of the AV01 mode, whereas the imaginary part that is anti-symmetric with respect to both planes does not have any effect on the generation of the AV01 mode.

    Real and imaginary parts of the typical optical force distributions for the three types of angular momentum (AM) transfer processes, PROC1, PROC2, and PROC3, which involve different higher-order hybrid optical modes: HE21, EH11, and HE31 modes, respectively. All the calculations are performed at the wavelength of λ=1550 nm and for the SWHW parameters of Rin=1.2 μm and Rout=3.0 μm. In addition, the photoelastic coefficients are set as p11=0.121 and p12=0.270. For the electrostriction body force, the blue color intensities and black arrows depict the magnitude and direction, respectively, of the optical force. The optical boundary force is composed of the electrostriction force f(es) (blue) and radiation pressure f(rp) (red), where their magnitude and direction are plotted by the arrows in the same scale for all the three processes for comparison. We choose a particular global phase in such a way that the longitudinal component of the total electric field in the waveguide is purely real.

    Figure 4.Real and imaginary parts of the typical optical force distributions for the three types of angular momentum (AM) transfer processes, PROC1, PROC2, and PROC3, which involve different higher-order hybrid optical modes: HE21, EH11, and HE31 modes, respectively. All the calculations are performed at the wavelength of λ=1550  nm and for the SWHW parameters of Rin=1.2  μm and Rout=3.0  μm. In addition, the photoelastic coefficients are set as p11=0.121 and p12=0.270. For the electrostriction body force, the blue color intensities and black arrows depict the magnitude and direction, respectively, of the optical force. The optical boundary force is composed of the electrostriction force f(es) (blue) and radiation pressure f(rp) (red), where their magnitude and direction are plotted by the arrows in the same scale for all the three processes for comparison. We choose a particular global phase in such a way that the longitudinal component of the total electric field in the waveguide is purely real.

    We investigate the feasibility of our scheme of the optoacoustic OVM generation at the telecom wavelength (1550  nm) and assess the stability of guidance of the OVMs in silica-glass SWHWs. For typical mode conversion via FSIBS, |R/Q|1 in Eq. (2), because the two pairs of optical modes are set equal to each other (e1e3 and e2e4) [41]. In this case, the input powers of the two pump beams should be strongly unbalanced to obtain the 100% mode conversion from a3 to a4. However, this requires a very long waveguide for building up the acoustic amplitude b(z) significantly. In strong contrast, in our case of the OVM generation, the optoacoustic overlap integral between the two pump beams |Q|=|e1|Δε¯|e2| is exactly half of that involving the OVM |R|=|e3|Δε¯|e4| because the OVM, OVlmσ, is a linear combination of the two hybrid optical modes as Eq. (1). As a result, 100% conversion from a3 to a4 can be achieved without the need for the unbalanced pump powers and impractically long interaction lengths. To quantitatively evaluate the feasibility, we define 3 dB coupling length, ζ, as the propagation length over which half the input signal power is converted into the targeted OVM or the idler wave (i.e., the other constituent OVM) at balanced input powers of the two pump beams, i.e., P10=P20, where Pν0Pν(z=0) is the input power of the ν-th optical wave. To examine the robustness of the OVM transmission in the waveguide, we check the two following requirements. First, it is essential to have a sufficiently large effective index separation between the adjacent optical modes of the same OAM order. In our case, this requirement is ensured by the large index step between the silica ring and air [42]. It is widely accepted that the effective index separation of Δneff>104 is enough for stable OVM propagation, where Δneff is defined as the effective index difference between the adjacent vector modes, e.g., among the TE01, TM01, and HE21 modes for PROC1, and between the EH11 and HE31 modes for PROC2 and PROC3. In addition, the generated OVMs should possess high modal purity and low propagation loss simultaneously. Here, we define the modal purity as the power portion of the desired OVM in the converted signal wave, i.e., P4/(P4+P5). We estimate the modal power fraction confined in the core (p), which is in practice related to the energy leakage into the cladding of the photonic crystal fibers [43]. The larger modal power fraction leads to a lower propagation loss.

    We calculate Δneff, p, and ζ for the three AM transfer processes over a range of the core radius Rout and the radius ratio MRin/Rout. In all the following analyses, for simplicity, we fix the acoustic quality factor as 1000 for fused silica glass [44] and the total pump power as P0=P1+P2=1  W. The larger Rout and smaller M yield a larger power fraction p in the core. It turns out that in the silica-glass SWHWs with Rout3  μm, the Kerr nonlinear coefficients, γνSPM, γμνXPM, and γνFWM, are estimated below 0.05  W1·m1, whereas the FSIBS peak gain coefficient g0 is much larger than them, e.g., g010  W1·m1 for PROC1 and g01  W1·m1 for PROC2 and PROC3, so the optical Kerr effect is negligible in Eq. (2). The modal purity of the generated OVM is then almost 100%, as the FWM-induced generation of the idler wave is suppressed. On the other hand, in the SWHWs made of a higher-index material like silicon, the Kerr nonlinear coefficients become comparable to the FSIBS peak gain coefficient, and thus the optical Kerr effect cannot be ignored, which we will not discuss here in more detail. In PROC1 [Fig. 5(a)], a large effective index separation Δneff>104 can be achieved over a wide range of (Rout,M). As the radius ratio M rises, the radiation pressure on the silica–air boundaries gets enhanced and more dominant in the FSIBS process, and thus ζ gets shorter at the expense of the transmission loss. The 3 dB mode conversion occurs typically at 1  m propagation over a broad range of (Rout,M), which is short enough to implement with the current fabrication technologies. In PROC2 and PROC3, we can also find the range of (Rout,M) suitable for our scheme. However, the range of useful waveguide parameters is not so wide compared to the PROC1 case, because the typical value of ζ is on the order of hundreds of meters. ζ decreases as Rout falls, and Rout should be small enough (<2  μm) in order to obtain sufficiently small values of ζ. Moreover, Δneff between the EH11 and the HE31 modes is relatively small and can be even zero. We note that at certain waveguide parameters, the FSIBS can be completely suppressed due to the counterbalance between the electrostriction and radiation pressure [38,39], as displayed by red striped regions, i.e., Rout>3  μm and M0.8 for PROC2 in Fig. 5(b) and Rout>3  μm and M0.05 for PROC3 in Fig. 5(c). We highlight that we here report for the first time, to our knowledge, the case of FSIBS elimination, while the cancellation of the conventional backward surface acoustic wave Brillouin scattering and forward intramodal (Raman-like) Brillouin scattering has been recently reported in Refs. [45,46] and Ref. [47], respectively.

    Color maps displaying the 3 dB coupling length ζ in silica-glass SWHWs over a range of Rout and M≡Rin/Rout for (a) PROC1, (b) PROC2, and (c) PROC3, where the pink curves indicate the cutoff of the HE21, EH11, and HE31 mode, respectively. The contour plots on top of the color maps show the minimum Δneff between the adjacent vector modes, i.e., among the TE01, TM01, and HE21 modes for (a), and between the EH11 and HE31 modes for (b) and (c). The black curves correspond to the power fraction p=85% in the waveguide core. The OVMs have lower confinement losses over the regions below the black curves.

    Figure 5.Color maps displaying the 3 dB coupling length ζ in silica-glass SWHWs over a range of Rout and MRin/Rout for (a) PROC1, (b) PROC2, and (c) PROC3, where the pink curves indicate the cutoff of the HE21, EH11, and HE31 mode, respectively. The contour plots on top of the color maps show the minimum Δneff between the adjacent vector modes, i.e., among the TE01, TM01, and HE21 modes for (a), and between the EH11 and HE31 modes for (b) and (c). The black curves correspond to the power fraction p=85% in the waveguide core. The OVMs have lower confinement losses over the regions below the black curves.

    We then investigate the spatial evolution of each optical wave. For simplicity, we set the waveguide parameters as (Rout,M)=(3  μm,0.4), at which for PROC1 the effective index separation, power fraction in the core, FSIBS peak gain coefficient, and 3 dB coupling length are calculated as Δneff=4×103, p=0.98, g0=15  W1·m1, and ζ=46  cm, respectively, and the Kerr nonlinear coefficients (γνSPM, γμνXPM, and γνFWM) are in the range of 48  W1  ·km1. Here, we use the long-waveguide approximation, assuming that the acoustic modes are damped more rapidly than the optical pump power changes [38]. This assumption is valid because the acoustic decay length (αb1=16  cm) is smaller than ζ. We calculate the spectral dependence of the mode conversion efficiency, P4/P30, over the entire range of the input pump power P10, while fixing the total pump power at P0=1  W and the propagation length at 3 m. Here, we set the frequency difference Δωijωiωj between the two pump beams (a1 and a2) as well as between the incident and coupled signals (a3 and a4) to exactly match the acoustic frequency Δω12=Δω34=Ω=2π×858  kHz. We determine the wavevector mismatch κ considering the optical dispersion of the silica-glass SWHWs. As previously mentioned, unity conversion efficiency can be obtained at P10=P20 [Fig. 6(a)]. At this balanced pump-power condition, the full width at half maximum (FWHM) operating bandwidth for the signal is determined as 26 GHz (0.21 nm) from Fig. 6(b). We can further optimize the waveguide length by setting the two input pump powers unbalanced (P10P20). We can in general think of two regimes of the input pump power ratio P10:P20, as shown in Fig. 6(c). When P10<P20, the acoustic power is concentrated at the beginning of the waveguide section, which leads to insufficient coupling strength for full mode conversion. On the other hand, when P10>P20, sufficient acoustic power can be sustained over a longer section of the waveguide to even over-couple the signal mode. In this case, it is possible to control the waveguide length to achieve unity conversion. We optimize the waveguide length and recalculate P4/P30 in Fig. 6(d), where we find the maximum FWHM operating bandwidth of 50 GHz (0.40 nm) at P10:P20=3:1.

    Spatial evolution of the optical and acoustic powers in PROC1 in silica-glass SWHW with (Rout,M)=(3 μm,0.4). (a) Optical and acoustic powers as functions of the propagation distance, z, at zero signal frequency detuning Δω13≡ω1−ω3=0. P1 and P2 are normalized to P0=1 W, P3 and P4 to P30=1 mW, and Pb to |b(z=0)|2. (b) Mode conversion efficiency P4(zmax)/P30 after propagation along zmax=3 m, where the optical dispersions of the silica-glass SWHWs are taken into account. (c) Same plots as (a) but at the different signal frequency detuning Δω13=10 GHz with different input pump power ratios of P10:P20=1:3 (dashed curves) and P10:P20=3:1 (solid curves). (d) Mode conversion efficiency at the optimum waveguide length P4(zopt)/P30, where zopt is the propagation length at which P4(z) is maximized.

    Figure 6.Spatial evolution of the optical and acoustic powers in PROC1 in silica-glass SWHW with (Rout,M)=(3  μm,0.4). (a) Optical and acoustic powers as functions of the propagation distance, z, at zero signal frequency detuning Δω13ω1ω3=0. P1 and P2 are normalized to P0=1  W, P3 and P4 to P30=1  mW, and Pb to |b(z=0)|2. (b) Mode conversion efficiency P4(zmax)/P30 after propagation along zmax=3  m, where the optical dispersions of the silica-glass SWHWs are taken into account. (c) Same plots as (a) but at the different signal frequency detuning Δω13=10  GHz with different input pump power ratios of P10:P20=1:3 (dashed curves) and P10:P20=3:1 (solid curves). (d) Mode conversion efficiency at the optimum waveguide length P4(zopt)/P30, where zopt is the propagation length at which P4(z) is maximized.

    The theory we have developed so far can be applied to any types of nonlinear intermodal optoacoustic coupling and arbitrary waveguide cross-sectional geometries. In particular, recent studies suggest that such nonlinear intermodal coupling could also be implemented in on-chip integrated waveguides, which can be fabricated using the standard CMOS-compatible technology [22,28]. While our approach of OVM generation has a unique advantage over other kinds of nonlinear optical methods, i.e., no initial OVM seed is required, there are still some important aspects that should be further addressed to put our scheme into practical use. First, it is essential to ensure that the hybrid modes are well controlled and generated. Mode converters would be needed at the input of the waveguide to couple light into the higher-order hybrid modes. Also, the mode filtering is required to separate the generated OVM (a4) from the residual fundamental mode (a3). We note that various types of multimode waveguide-based photonic devices [48,49] would provide possible solutions for this purpose, which have been experimentally demonstrated in the fiber [50] and integrated waveguide platform [51]. For example, the multimode Y-junctions based on the adiabatic evolution of effective modal indices can be used as mode combiners/splitters and mode converters [49]. Another issue is the narrow operating spectral bandwidth, as shown in Fig. 6. The limited bandwidth is due to the group velocity difference between the two hybrid modes, e.g., the HE11 and the HE21 modes, which increases the wavevector mismatch κ with the detuning of signal wavelength from the pump wavelength. While it is not straightforward to match the dispersions of the two optical modes in simple step-index SWHWs, in photonic crystal fibers, the group velocities of the two different optical modes can be made equal to each other by adequately adjusting the pitch and hole size in the cladding [52]. The operating spectral bandwidth would be then significantly improved, potentially beyond 10  nm [52]. One might think that the strong pump beams may give rise to the temperature increase of the waveguide, which alters the effective indices of the optical modes. However, we determine from the heat diffusion analysis [53] that the influence of such a thermo-optic effect is not significant in practical waveguides such as silica-glass photonic crystal fibers [32], as the temperature increase is far below 1°C at the 1 W pump power. Finally, the robust guidance of the hybrid modes in optical waveguides is impeded by the random linear mode coupling due to the residual defects and random birefringence fluctuations.

    4. CONCLUSION

    In conclusion, nonlinear optoacoustic interactions via FSIBS can be employed for the generation of OAM-carrying optical vortices, which has a unique advantage of excluding the need for any initial optical vortex seed. We have presented numerical results of the intermodal coupling involving the OVM and the AVM, and it turns out that a waveguide with a subwavelength-hole of proper diameter at the core center is suitable for not only the robust guidance but also the efficient optoacoustic generation of the OVMs. We point out that the FSIBS process is the right nonlinear effect for all-optical OVM generation, where nearly unity conversion efficiency and high modal purity are guaranteed with the aid of the energy-momentum conservation. Our scheme can be realized in various optical systems, such as microstructured fibers, microcavities, and integrated optical circuits.

    Our findings open an avenue for providing a key technology for next-generation ultrahigh-bandwidth optical communication systems. For instance, OAM-selective photon–phonon interactions could be employed to fabricate the spatial mode multiplexers [7]. We have shown that coupling with the optoacoustically generated AVM enables the selective excitation of one of the constituent OVMs of the pump beams. Such all-optical control of the AVMs in the waveguides might also be exploited for the application of axial torque to axisymmetric objects and the transport of modulated helicity [54]. Furthermore, introducing AM transfer to the intermodal coupling allows us to generate and manipulate non-classical correlations between OAM-carrying photons and phonons, and to explore new types of quantum phenomena [55].

    References

    [1] K. Toyoda, K. Miyamoto, N. Aoki, R. Morita, T. Omatsu. Using optical vortex to control the chirality of twisted metal nanostructures. Nano Lett., 12, 3645-3649(2012).

    [2] M. R. Dennis, R. P. King, B. Jack, K. O’Holleran, M. J. Padgett. Isolated optical vortex knots. Nat. Phys., 6, 118-121(2010).

    [3] G. Foo, D. M. Palacios, G. A. Swartzlander. Optical vortex coronagraph. Opt. Lett., 30, 3308-3310(2005).

    [4] S. W. Hell, J. Wichmann. Breaking the diffraction resolution limit by stimulated emission: stimulated-emission-depletion fluorescence microscopy. Opt. Lett., 19, 780-782(1994).

    [5] H. He, M. E. J. Friese, N. R. Heckenberg, H. Rubinsztein-Dunlop. Direct observation of transfer of angular momentum to absorptive particles from a laser beam with a phase singularity. Phys. Rev. Lett., 75, 826-829(1995).

    [6] M. Mirhosseini, O. S. Magaña-Loaiza, M. N. O’Sullivan, B. Rodenburg, M. Malik, M. P. J. Lavery, M. J. Padgett, D. J. Gauthier, R. W. Boyd. High-dimensional quantum cryptography with twisted light. New J. Phys., 17, 033033(2015).

    [7] N. Bozinovic, Y. Yue, Y. Ren, M. Tur, P. Kristensen, H. Huang, A. E. Willner, S. Ramachandran. Terabit-scale orbital angular momentum mode division multiplexing in fibers. Science, 340, 1545-1548(2013).

    [8] M. W. Beijersbergen, R. P. C. Coerwinkel, M. Kristensen, J. P. Woerdman. Helical-wavefront laser beams produced with a spiral phaseplate. Opt. Commun., 112, 321-327(1994).

    [9] L. Marrucci, C. Manzo, D. Paparo. Optical spin-to-orbital angular momentum conversion in inhomogeneous anisotropic media. Phys. Rev. Lett., 96, 163905(2006).

    [10] H. Xu, L. Yang. Conversion of orbital angular momentum of light in chiral fiber gratings. Opt. Lett., 38, 1978-1980(2013).

    [11] C. N. Alexeyev, M. A. Yavorsky. Generation and conversion of optical vortices in long-period helical core optical fibers. Phys. Rev. A, 78, 043828(2008).

    [12] A. W. Lohmann, D. P. Paris. Binary Fraunhofer holograms, generated by computer. Appl. Opt., 6, 1739-1748(1967).

    [13] K. Dholakia, N. B. Simpson, M. J. Padgett, L. Allen. Second-harmonic generation and the orbital angular momentum of light. Phys. Rev. A, 54, R3742-R3745(1996).

    [14] W. Jiang, Q. F. Chen, Y. S. Zhang, G. C. Guo. Computation of topological charges of optical vortices via nondegenerate four-wave mixing. Phys. Rev. A, 74, 043811(2006).

    [15] J. Arlt, K. Dholakia, L. Allen, M. J. Padgett. Parametric down-conversion for light beams possessing orbital angular momentum. Phys. Rev. A, 59, 3950-3952(1999).

    [16] Y. Zhang, Z. Nie, Y. Zhao, C. Li, R. Wang, J. Si, M. Xiao. Modulated vortex solitons of four-wave mixing. Opt. Express, 18, 10963-10972(2010).

    [17] Z. Zhang, D. Ma, Y. Zhang, M. Cao, Z. Xu, Y. Zhang. Propagation of optical vortices in a nonlinear atomic medium with a photonic band gap. Opt. Lett., 42, 1059-1062(2017).

    [18] D. Zhang, X. Liu, L. Yang, X. Li, Z. Zhang, Y. Zhang. Modulated vortex six-wave mixing. Opt. Lett., 42, 3097-3100(2017).

    [19] Z. Wang, J. Yang, Y. Sun, Y. Zhang. Interference patterns of vortex beams based on photonic band gap structure. Opt. Lett., 43, 4354-4357(2018).

    [20] M. S. Kang, A. Nazarkin, A. Brenn, P. St.J. Russell. Tightly trapped acoustic phonons in photonic crystal fibres as highly nonlinear artificial Raman oscillators. Nat. Phys., 5, 276-280(2009).

    [21] H. Shin, J. A. Cox, R. Jarecki, A. Starbuck, Z. Wang, P. T. Rakich. Control of coherent information via on-chip photonic-phononic emitter-receivers. Nat. Commun., 6, 6427(2015).

    [22] E. A. Kittlaus, N. T. Otterstrom, P. Kharel, S. Gertler, P. T. Rakich. Non-reciprocal interband Brillouin modulation. Nat. Photonics, 12, 613-619(2018).

    [23] W. Gao, C. Mu, H. Li, Y. Yang, Z. Zhu. Parametric amplification of orbital angular momentum beams based on light-acoustic interaction. Appl. Phys. Lett., 107, 041119(2015).

    [24] M. F. Andersen, C. Ryu, P. Cladé, V. Natarajan, A. Vaziri, K. Helmerson, W. D. Phillips. Quantized rotation of atoms from photons with orbital angular momentum. Phys. Rev. Lett., 97, 170406(2006).

    [25] J. T. Mendonça, B. Thidé, H. Then. Stimulated Raman and Brillouin backscattering of collimated beams carrying orbital angular momentum. Phys. Rev. Lett., 102, 185005(2009).

    [26] P. Z. Dashti, F. Alhassen, H. P. Lee. Observation of orbital angular momentum transfer between acoustic and optical vortices in optical fiber. Phys. Rev. Lett., 96, 043604(2006).

    [27] K. Y. Song, W. Zou, Z. He, K. Hotate. All-optical dynamic grating generation based on Brillouin scattering in polarization-maintaining fiber. Opt. Lett., 33, 926-929(2008).

    [28] E. A. Kittlaus, N. T. Otterstrom, P. T. Rakich. On-chip inter-modal Brillouin scattering. Nat. Commun., 8, 15819(2017).

    [29] J. Petersen, J. Volz, A. Rauschenbeutel. Chiral nanophotonic waveguide interface based on spin-orbit interaction of light. Science, 346, 67-71(2014).

    [30] K. Y. Bliokh, D. Smirnova, F. Nori. Quantum spin Hall effect of light. Science, 348, 1448-1451(2015).

    [31] Y. Yue, L. Zhang, Y. Yan, N. Ahmed, J.-Y. Yang, H. Huang, Y. Ren, S. Dolinar, M. Tur, A. E. Willner. Octave-spanning supercontinuum generation of vortices in an As2S3 ring photonic crystal fiber. Opt. Lett., 37, 1889-1891(2012).

    [32] G. S. Wiederhecker, C. M. B. Cordeiro, F. Couny, F. Benabid, S. A. Maier, J. C. Knight, C. H. B. Cruz, H. L. Fragnito. Field enhancement within an optical fibre with a subwavelength air core. Nat. Photonics, 1, 115-118(2007).

    [33] K. Okamoto. Fundamentals of Optical Waveguides(2006).

    [34] R. A. Waldron. Some problems in the theory of guided microsonic waves. IEEE Trans. Microwave Theory Tech., 17, 893-904(1969).

    [35] T. F. S. Büttner, I. V. Kabakova, D. D. Hudson, R. Pant, C. G. Poulton, A. C. Judge, B. J. Eggleton. Phase-locking and pulse generation in multi-frequency Brillouin oscillator via four wave mixing. Sci. Rep., 4, 5032(2014).

    [36] H. Shin, W. Qiu, R. Jarecki, J. A. Cox, R. H. Olsson, A. Starbuck, Z. Wang, P. T. Rakich. Tailorable stimulated Brillouin scattering in nanoscale silicon waveguides. Nat. Commun., 4, 1944(2013).

    [37] R. Van Laer, B. Kuyken, D. Van Thourhout, R. Baets. Interaction between light and highly confined hypersound in a silicon photonic nanowire. Nat. Photonics, 9, 199-203(2015).

    [38] C. Wolff, M. J. Steel, B. J. Eggleton, C. G. Poulton. Stimulated Brillouin scattering in integrated photonic waveguides: forces, scattering mechanisms, and coupled-mode analysis. Phys. Rev. A, 92, 013836(2015).

    [39] P. T. Rakich, C. Reinke, R. Camacho, P. Davids, Z. Wang. Giant enhancement of stimulated Brillouin scattering in the subwavelength limit. Phys. Rev. X, 2, 011008(2012).

    [40] S. Afshar, T. M. Monro. A full vectorial model for pulse propagation in emerging waveguides with subwavelength structures part I: Kerr nonlinearity. Opt. Express, 17, 2298-2318(2009).

    [41] X. Huang, S. Fan. Complete all-optical silica fiber isolator via stimulated Brillouin scattering. J. Lightwave Technol., 29, 2267-2275(2011).

    [42] S. Ramachandran, P. Kristensen, M. F. Yan. Generation and propagation of radially polarized beams in optical fibers. Opt. Lett., 34, 2525-2527(2009).

    [43] H. Zhang, X. Zhang, H. Li, Y. Deng, X. Zhang, L. Xi, X. Tang, W. Zhang. A design strategy of the circular photonic crystal fiber supporting good quality orbital angular momentum mode transmission. Opt. Commun., 397, 59-66(2017).

    [44] R. M. Shelby, M. D. Levenson, P. W. Bayer. Guided acoustic-wave Brillouin scattering. Phys. Rev. B, 31, 5244-5252(1985).

    [45] O. Florez, P. F. Jarschel, Y. A. V. Espinel, C. M. B. Cordeiro, T. P. Mayer Alegre, G. S. Wiederhecker, P. Dainese. Brillouin scattering self-cancellation. Nat. Commun., 7, 11759(2016).

    [46] D. S. Han, I.-M. Lee, K. H. Park, M. S. Kang. Extremely polarization-sensitive surface acoustic wave Brillouin scattering in subwavelength waveguides. Appl. Phys. Lett., 113, 121108(2018).

    [47] D. S. Han, I.-M. Lee, K. H. Park, M. S. Kang. Polarization-selective control of nonlinear optomechanical interactions in subwavelength elliptical waveguides. Opt. Express, 27, 1718-1726(2019).

    [48] L. B. Soldano, E. C. M. Pennings. Optical multi-mode interference devices based on self-imaging: principles and applications. J. Lightwave Technol., 13, 615-627(1995).

    [49] J. D. Love, N. Riesen. Single-, few-, and multimode Y-junctions. J. Lightwave Technol., 30, 304-309(2012).

    [50] R. Ismaeel, T. Lee, B. Oduro, Y. Jung, G. Brambilla. All-fiber fused directional coupler for highly efficient spatial mode conversion. Opt. Express, 22, 11610-11619(2014).

    [51] L.-W. Luo, N. Ophir, C. P. Chen, L. H. Gabrielli, C. B. Poitras, K. Bergmen, M. Lipson. WDM-compatible mode-division multiplexing on a silicon chip. Nat. Commun., 5, 3069(2014).

    [52] M. Bache, H. Nielsen, J. Lægsgaard, O. Bang. Tuning quadratic nonlinear photonic crystal fibers for zero group-velocity mismatch. Opt. Lett., 31, 1612-1614(2006).

    [53] D. C. Brown, H. J. Hoffman. Thermal, stress, and thermo-optic effects in high average power double-clad silica fiber lasers. IEEE J. Sel. Top. Quantum Electron., 37, 207-217(2001).

    [54] L. Zhang, P. L. Marston. Angular momentum flux of nonparaxial acoustic vortex beams and torques on axisymmetric objects. Phys. Rev. E, 84, 065601(2011).

    [55] R. Riedinger, S. Hong, R. A. Norte, J. A. Slater, J. Shang, A. G. Krause, V. Anant, M. Aspelmeyer, S. Gröblacher. Non-classical correlations between single photons and phonons from a mechanical oscillator. Nature, 530, 313-316(2016).

    Dae Seok Han, Myeong Soo Kang. Reconfigurable generation of optical vortices based on forward stimulated intermodal Brillouin scattering in subwavelength-hole photonic waveguides[J]. Photonics Research, 2019, 7(7): 754
    Download Citation