Abstract
We study the generation of spin-orbit (SO) modes via four-wave mixing (FWM)-based parametric amplification. SO modes carry quantized total angular momentum (TAM), and we show that FWM processes that generate new signals conserve TAM. This is a generalization of prior research which operated in a regime where FWM processes conserved spin and orbital angular momenta independently. We calculate the growth rates of new modes for both degenerate and nondegenerate pump configurations. Our theory is validated against numerical simulations for the cases where the generated signals are in the same SO mode(s) as the pump(s). We also calculate the growth rates of signals in SO modes other than the pumps.
© 2024 Optica Publishing Group under the terms of the Optica Open Access Publishing Agreement
1. Introduction
Light that carries orbital angular momentum (OAM), also referred to as optical vortices or twisted light, is characterized by a helical or twisted wavefront. In contrast to spin angular momentum (SAM), where photons are limited to two states, OAM allows for, in principle, an infinite set of spatially orthogonal states [1]. This feature makes OAM-carrying light a highly attractive platform for spatial division multiplexing and quantum key distribution in free space as well as guided-wave communication systems [2–4]. Vortex beams also have applications in particle trapping and optical tweezers [5], optical metrology [6], quantum optics [7,8], nanoscopy [9,10], and remote sensing [11,12].
For many of these applications, the ability to create light with a desired combination of SAM and OAM at specific frequencies and high powers is key. Generation and amplification of free space OAM beams, which are defined as having an azimuthally varying phase structure exp[iℓφ] where ℓ and φ are the OAM topological charge and azimuthal coordinate respectively, has been explored using Raman amplification in plasma [13] and four-wave mixing (FWM) in 85Rb atomic vapor [14]. In azimuthally symmetric fibers, high refractive index contrasts produce spin-orbit (SO) interactions. In simple terms, this refers to the dependence of the effective index of a fiber mode of a given OAM topological charge ℓ on its SAM [15,16]. Raman amplification of “SO modes” in fibers have been studied and demonstrated recently [17].
Of particular interest here is FWM-based parametric amplification of SO modes in fiber, for which a typical experimental configuration is provided in Fig. 1 [18,19]. High power SO modes are injected into a fiber that supports SO modes; one mode, labeled (+,+2), has a dominant left-circular polarization (LCP) component and OAM ℓ = + 2, while the other, labeled (–,+5), has a dominant right-circular polarization (RCP) component and OAM ℓ = + 5. Within the fiber, FWM processes can amplify an ambient noise field into Stokes and anti-Stokes pairs of new SO modes. Conservation of energy dictates that the Stokes and anti-Stokes modes will have a symmetric detuning frequency Ω when the pumps have the same frequency, while conservation of linear and angular momentum will dictate which of the supported SO modes are excited. An optical spectrum analyzer, combined with bandpass filters and an SO mode imaging scheme, can fully characterize the fields produced in the fiber. We note that Raman scattering will also likely occur within the fiber, producing a broader Stokes-shifted peak which will be in a superposition of the pump SO modes [17,18].
Although FWM-based parametric amplification has been described by applying the multimode nonlinear Schrödinger equations (MM-NLSEs) [18,20] to SO fiber modes, there remains a need for a more complete theoretical description of nonlinearly interacting SO modes that provides insight into how the angular momentum of light is transformed during the process. This is important not only for the aforementioned applications, but also in the context of the recent surge in interest in multimode nonlinear photonics [21–23]. In this work, we present, for the first time to our knowledge, a generalized theory of FWM-based parametric amplification of SO modes in azimuthally symmetric waveguides. We consider degenerate and nondegenerate pump cases for the spontaneous generation of new SO modes. We show that these processes conserve total angular momentum (TAM), but under certain conditions conserve SAM and OAM independently. We also compare our theory to numerical simulations for the cases where the generated Stokes and anti-Stokes signals are in the same SO modes as the pumps – equivalent to scalar modulation instability (SMI) and cross phase modulation instability (XPMI) of the pump beams.
This paper is organized as follows. First, we review the properties of SO modes in azimuthally symmetric waveguides. We then outline the derivation and results for the FWM-based parametric amplification processes and growth rates. These analytical results are compared to Fourier split-step simulations for the SMI and XPMI cases utilizing a fiber designed to accentuate SO interactions. Finally, we consider the generation of signals in SO modes different from the pump(s) by way of intermodal FWM.
2. Spin-orbit modes in azimuthally symmetric waveguides
Fibers with azimuthally symmetric refractive index profiles support modes ${{\mathbf E}_q}({{\mathbf r},t} )= {\tilde{{\mathbf F}}_q}({{{\mathbf r}_ \bot },{\omega_0}} )\exp [{i({{\beta_q}({{\omega_0}} )z - {\omega_0}t} )} ]$, where ${{\mathbf E}_q}({{\mathbf r},t} )$ is the complex envelope of the electric field and ${\tilde{{\mathbf F}}_q}$ is the mode’s electric field profile which has the form
J denotes the total quantized angular momentum. Each photon in the mode carries $J\hbar $ angular momentum in the $\hat{{\mathbf z}}$ direction. An equivalent representation of ${\tilde{{\mathbf F}}_q}$ is produced by transforming to the circular polarization basis,
Both TE and TM modes carry zero TAM. On the other hands, modes with J ≠ 0 do carry a net TAM; in addition, one spin component will be dominant in the sense that, e.g., $\int {{d^2}r\,{{|{\tilde{F}_q^ \pm } |}^2}} > \int {{d^2}r\,{{|{\tilde{F}_q^ \mp } |}^2}}$. It is common to express such spatial modes in terms of the OAM quanta ℓ carried by the dominant spin component rather than J, e.g. for a LCP-dominant mode
3. Four wave mixing of spin-orbit modes: theory
We now consider how SO modes interact in FWM processes. The derivation for the governing equations is based on the generalized multimode nonlinear Schrödinger equation (MMNLSE) as described in Refs [26,27]. We consider the electric and magnetic fields inside the fiber as the sum of spatial modes q, each at a frequency ωq which are in the vicinity of a common carrier frequency ω0,
Note that, by requiring the envelopes to have the form Aq(z), we are assuming continuous wave (CW) signals; we will lift this assumption later in the paper to allow for temporal variations in the envelopes. Proceeding with the derivation of the MMNLSE, we arrive at
In this work, we are concerned with FWM processes driven by the instantaneous third-order nonlinear susceptibility. Thus the nonlinear polarization PNL is given by
Combining Eqs. (9)–(11) and assuming β1(ω1) ≠ β2(ω2), the governing equations for the pumps are given by
Since $Q_q^{\textrm{(SPM)}}$ and $Q_{q,j}^{\textrm{(XPM)}}$ are real-valued, |A1(z)|2 and |A2(z)|2 are z-invariant and Eqs. (12) can be solved directly even though they form a coupled nonlinear system. The equation for an arbitrary signal S will include cross phase modulation by the two pumps and numerous other nonlinear terms. If we only consider signals that can undergo amplification, even from noise, one obtains
With ωS < ωA, S can naturally be identified as the Stokes mode and A the anti-Stokes. The governing equation for mode A is obtained by the substitution S ↔ A. In addition, we can identify FWM1(2) as corresponding to degenerate FWM from pump 1(2), and FWM3 to nondegenerate FWM.
Angular momentum selection rules, embedded within the overlap integrals Q, dictate which SO modes interact with each other via FWM. Considering the forms of the spatial modes provided in Eq. (1), the azimuthal integration for $Q_{\textrm{S,A}}^{\textrm{(FWM1)}}$ is nonzero only if 2J1 = JS + JA. Similarly, $Q_{\textrm{S,A}}^{\textrm{(FWM2)}}$ requires 2J2 = JS + JA, and $Q_{\textrm{S,A}}^{\textrm{(FWM3)}}$ requires J1 + J2 = JS + JA. Thus, since J determines the TAM of the mode, one can say the TAM of the light is conserved during the FWM processes. When the modes involved in a particular FWM process carry nonzero TAM and have negligibly small recessive spin components such that the approximate form ${\tilde{{\mathbf F}}_{({p, \pm ,\ell } )}}({{\mathbf r}_ \bot },{\omega _0}) \mathbf \simeq {\hat{\boldsymbol{\sigma}}_ \pm }{\tilde{F}_{({p, \pm ,\ell } )}}(r,{\omega _0})\exp [{i\ell \varphi } ]$ is accurate, the overlap integrals are nonzero only when OAM and SAM are conserved independently. For example, in this regime, Q(FWM3) ≠ 0 only when ℓ1+ ℓ2 = ℓS + ℓA and s1 + s2 = sS + sA, where s = + 1 for LCP-dominant modes and s = −1 for RCP-dominant modes. This latter case was observed experimentally in Ref. [18].
Assuming the frequency matching conditions (delta functions in Eq. (14)) are satisfied, and that the pumps carry undepleted powers P1 and P2, one can show that each FWM term can lead to the growth of new SO modes. The resultant growth rates for each FWM process are
A Stokes and anti-Stokes mode pair will stably grow when g is real-valued. The growth rate is maximized under the phase matching condition κ = 0, and larger pump powers will increase the bandwidth over which the modes will grow. For SO modes, additional mode selectivity is achieved by the constraint of conservation of TAM. Note that the growth rates in Eqs. (16)–(18) are the power e-folding lengths, e.g. PS(z)=PS(0) exp[gz].
4. Simulation overview and fiber description
Directly solving Eqs. (12) and (14) would require construction of a Nm × Nω system of equations, where Nm is the number of modes and Nω the number of discrete frequency bins. One can reduce this to a Nm system by allowing the mode envelopes to have a time dependence, i.e. Aq(z) → Aq(z,t). Doing so and rederiving Eqs. (12) and (14) yields
Since the Fourier split-step simulation technique requires a finite-duration pulse, we use a T = 0.8 ns length supergaussian pulse, i.e. ${A_q}({z = 0,t} )\propto \,\exp [{ - {{({{{2t} / T}} )}^{16}}} ]$. For each test case we perform 100 simulations with white Gaussian noise added to the beam profiles before propagation. This acts as a seed for signal generation. We then analyze the ensemble average spectra after propagation. When performing the Fourier transforms of the final envelopes, we only utilize the beam content within the pulse width to better approximate CW behavior. Dispersion terms up to fourth order (β4) are included to ensure that the phase matching processes within the simulations are accurate across the bandwidth of interest.
As for the fiber, we consider an air-core ring fiber as shown in Fig. 2(a) [15–19,24]. It consists of an 8 µm radius air-filled core surrounded by a 2 µm thick high refractive index ring of a GeO2-SiO2 mixture, followed by a cladding of pure SiO2. The ring-shaped refractive index profile allows for the stable guiding of SO modes which have characteristic donut intensity profiles. Examples of mode profiles are shown in Fig. 2(b). We observe that, for the (0,+,–2) mode, the recessive spin component (RCP in this case) is comparable to, but still smaller than, the dominant spin component (LCP). On the other hand, the higher-$\ell $ mode (0,+,+8) has a fairly negligible recessive spin component. To characterize this effect across all modes in the fiber, we quantify the relative strengths of each mode’s ${\hat{{\mathbf e}}_j} = \{{{{\hat{{\boldsymbol \sigma }}}_ + },{{\hat{{\boldsymbol \sigma }}}_ - },\hat{\mathbf{z}}} \}$ component as
The results for a selection of modes are shown in Fig. 2(c). As discussed earlier, the two spin components for TE and TM modes have the same strengths. Additionally, as ℓ increases, the recessive spin component does become weak. Generally speaking, modes with higher ℓ can be well approximated by the profile ${\tilde{{\mathbf F}}_{({p, \pm ,\ell } )}}({{\mathbf r}_ \bot },{\omega _0}) \mathbf \simeq {\hat{\boldsymbol{\sigma}}_ \pm }{\tilde{F}_{({p, \pm ,\ell } )}}(r,{\omega _0})\exp [{i\ell \varphi } ]$, and FWM processes amongst these modes will, also to good approximation, conserve OAM and SAM independently. On the other hand, TE, TM, and modes with low ℓ are better described by their full forms Eqs. (3)–(6), and will not have such a restriction. What quantifies as “high” versus “low” ℓ will vary depending on the fiber design and the level of analytical precision desired.
Further details on the fiber design and its modal properties are presented in the Supplement 1.
5. Scalar modulation instability of the (0,–,+8) mode
We first explore the case of a single pump: a 0.8 ns supergaussian pulse in the (0,–,+8) spatial mode at 1.064 µm, with 1 kW peak power propagating over 3 meters. We consider only new Stokes and anti-Stokes signals which are also in the (0,–,+8) spatial mode; this is equivalent to scalar modulation instability (SMI). The relevant dispersion and overlap coefficients are provided in the first column of Table 1. Note that, for SMI, Q(SPM) = Q(FWM1). Since this mode is in the anomalous dispersion regime it will undergo the typical scalar modulation instability and develop broad, symmetric sidebands about the pump frequency with peaks at λS = 1.108 µm and λA = 1.023 µm. In addition, the dispersion curve is phase matched at a further detuning with peaks at λS = 1.305 µm, λA = 0.898 µm. The simulation results are compared to the predictions of Eq. (16) in Fig. (3).
We observe that the locations of the new modes are accurately predicted, including the spectral widths, although the gains are slightly overpredicted. The extra oscillations present are attributed to other FWM processes not considered in the analytic derivations. That is, any phase-matched set of frequencies β(ωa) + β(ωb) + θ(ωa) + θ(ωb) = β(ωc) + β(ωd) + θ(ωc) + θ(ωd) will interact through FWM. The analytic results presented are meant to capture the first-order effects of stably growing modes generated by the pump, i.e. the case where ωa = ωb = ω0, where ω0 is the pump frequency. Furthermore, these and following simulations were performed using a variety of sampling configurations, e.g. z-step size and time-domain sample spacing Δt, to ensure such features are not the result of numerical artifacts.
When considering signals that may be generated in new spatial modes, and the large number of spatial modes that are available in the fiber, it may be desirable to ensure the new signals have smaller bandwidths to prevent spectral overlap and avoid possible interference in signal growth. In Eqs. (16)–(18) we observe that decreasing the pump power will decrease the signal bandwidth, though a longer fiber will be needed to reach the same power levels. To test this, we repeat the prior case but with a 100 W pump and 30 m fiber length. The results are shown in Fig. (4). As predicted, the peak gains remain unchanged while the spectral widths are reduced to ≲1 nm for the further detuned signals.
6. Cross-phase modulation instability of the (0,±,+2) spatial modes
We next explore the case of two pumps: one in the (0, –,+2) spatial mode, and the other in the (0,+,+2) spatial mode. Here, again, both pumps consist of 0.8 ns supergaussian pulses centered at 1.064 µm, have 100 W peak power, and propagate 30 m. We consider only new Stokes and anti-Stokes signals in either (0,±,+2) spatial modes. This is equivalent to cross phase modulation instability (XPMI). The relevant dispersion and overlap coefficients are provided in the second and third columns of Table 1 (note that, for XPMI, Q(XPM) = Q(FWM3)).
The simulation results are presented in Fig. (5). Neither mode has phase-matched SMI processes as the (0,–,+8) spatial mode does; however, there is a XPMI process predicted to occur, based on Eq. (18), with the (0,–,+2) spatial mode at λS = 1.1138 µm and the (0,+,+2) spatial mode at λA = 1.0185 µm. This is consistent with XPMI observed in other fiber types, where the slow mode (higher group index and β1) is Stokes shifted, while the fast mode (lower group index and β1) is anti-Stokes shifted [19,29]. Equation (18) also predicts a set of peaks detuned by ∼0.1 nm from the pump. These correspond to the first-order approximation to the phase matching condition, (β1,S – β1,A) Ω = θ1 + θ2 – θS – θA, where Ω is the detuning. These signals are not present in the simulation: neither at the end of the 30 m fiber, nor any intermediately sampled distance. Since the derivation of the growth rates provided in Eqs. (16)–(18) effectively assume distinct modes which are sufficiently separated in frequency, the very narrow detuning of these modes, and even partial overlap in the gain curves with the pump, may feed other nonlinear processes which suppresses their growth.
7. Generation of new spatial modes
We now study the generation of spatial modes different than the pump(s), continuing with the case of pumps in the (0,±,+2) spatial modes at 1.064 µm. The FWM overlap integrals, Eqs. (15), are evaluated for all mode permutations and, when they are nonzero, Eqs. (16)–(18) are calculated for the peak growth rates. The results for degenerate and nondegenerate FWM processes are shown in Figs. (6) and (7) respectively. The growth rates are normalized to the pump power and are accurate for pumps up to ∼1 kW (changes in the Stokes and anti-Stokes frequencies due to nonlinear phase shifts are at most ∼1 nm for 1 kW pumps for this test case). The x-axis denotes the Stokes and anti-Stokes modes generated. Note that all generated modes have zero radial order, and thus p is omitted from the mode labels to make the figures more legible. In addition, the last numbers in the labels are the wavelengths in µm. To emphasize how angular momentum evolves between the pump(s) and amplified modes, processes which conserve OAM and SAM independently are marked with blue circles while processes which do not are marked with red x’s.
Figure (6) (left) provides the growth rates for modes generated by degenerate FWM of the (0,+,+2) pump. Recall that the TAM of an SO mode is quantized by J in Eqs. (1) and (2), i.e. J = s + ℓ, where s is the dominant spin component and ℓ is the transverse OAM of that dominant spin component, with J = 0 for TE and TM modes. Thus, to obtain a nonzero FWM overlap integral, the selection rule 2J1 = JA + JS = + 6 must be satisfied. All of the modes shown satisfy this condition. The differences in gain are due to variations in the mode profiles and thus overlap integral values. An equivalent selection rule of 2J2 = JA + JS = + 2 is enforced for degenerate FWM from the (0, –,+2) pump, which is shown in Fig. (6) (right). Figure (7) shows the growth rates for modes generated by nondegenerate FWM of the two pumps. Here, the selection rule J1 + J2 = JS + JA = + 4 is enforced. We note that Raman scattering of the pumps will likely produce a peak at ∼1.12 µm, which may mask or disrupt the generation of the modes shown [18].
Under different conditions, such as for a fiber profile with a thicker “ring” section and thus smaller spin-orbit coupling, or when only larger-ℓ modes are involved in the FWM processes, the recessive spin components of the mode profiles can reasonably be discarded, i.e. ${\tilde{{\mathbf F}}_{({p, \pm ,\ell } )}}({{\mathbf r}_ \bot },{\omega _0}) \mathbf \simeq {\hat{\boldsymbol{\sigma}}_ \pm }{\tilde{F}_{({p, \pm ,\ell } )}}(r,{\omega _0})\exp [{i\ell \phi } ]$. In this case the overlap integrals for FWM-based parametric amplification would enforce independent conservation of SAM and OAM, as has been demonstrated in prior experiments [16,18]. However, for the fiber design and pump configuration considered in this paper, such an approximation is unjustified. For example, one of the most efficient degenerate FWM processes of the (0,–,+2) pump will generate the (0,+,–2) and (0,–,+4) spatial modes at λS = 1.492 µm and λA = 0.827 µm, respectively. Such a process would be “forbidden” under the reduced case.
8. Conclusion
In this work, we presented a generalized theory of four-wave mixing (FWM)-based parametric amplification for modes in azimuthally symmetric optical fibers. Strictly speaking, the modes in such fibers are characterized by a total angular momentum (TAM) quanta J, where each photon carries $J\hbar $ angular momentum along the propagation axis. Such modes have nonzero LCP and RCP spin components. We showed that this feature translates to selection rules during FWM-based parametric amplification, wherein TAM must be conserved. This is a generalization of prior research [16,18], where the fiber design and mode subset utilized allowed for the modes to be well-approximated as “pure” spin-orbit (SO) modes rather than TAM modes. In such cases the selection rules are more restrictive and correspond to conservation of orbital angular momentum (OAM) and spin angular momentum (SAM) independently.
In addition, we derived growth rates of new modes via degenerate and nondegenerate parametric amplification. Numerical simulations were performed for the cases where the generated Stokes and anti-Stokes modes are in the same spatial mode as the pump(s) – equivalent to scalar and cross-phase modulation instability of the pump beam(s). Our numerical simulations are in good agreement with our theory. We also calculated the growth rates of all new SO modes from a nondegenerate pump configuration and found that processes which “break” the selection rules for independent conservation of SAM and OAM, but are allowed under conservation of TAM, can be significant.
Funding
Office of Naval Research.
Acknowledgments
The authors would like to thank Prof. Siddharth Ramachandran (Boston University) and Dr. Vladimir Markov (Advanced Systems & Technologies, Inc.) for useful discussions and the Office of Naval Research (Q. Saulter) for funding this UMD research.
Disclosures
The authors declare no conflicts of interest.
Data availability
The data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.
Supplemental document
See Supplement 1 for supporting content.
References
1. L. Allen, M. Beijersbergen, R. Spreeuw, et al., “Orbital angular momentum of light and the transformation of Laguerre-Gaussian laser modes,” Phys. Rev. A 45(11), 8185–8189 (1992). [CrossRef]
2. G. Gibson, J. Courtial, and M. Padgett, “Free-space information transfer using light beams carrying orbital angular momentum,” Opt. Express 12(22), 5448–5456 (2004). [CrossRef]
3. Y. Ren, Z. Wang, P. Liao, et al., “Experimental characterization of a 400 Gbit/s orbital angular momentum multiplexed free-space optical link over 120 m,” Opt. Lett. 41(3), 622–625 (2016). [CrossRef]
4. M. Mirhosseini, M. Malik, Z. Shi, et al., “Efficient separation of the orbital angular momentum eigenstates of light,” Nat. Commun. 4(1), 2781 (2013). [CrossRef]
5. D. G. Grier, “A revolution in optical manipulation,” Nature 424(6950), 810–816 (2003). [CrossRef]
6. E. Otte and C. Denz, “Optical trapping gets structure: structured light for advanced optical manipulation,” Appl. Phys. Rev. 7(4), 041308 (2020). [CrossRef]
7. M. Krenn, M. Malik, M. Erhard, et al., “Orbital angular momentum of photons and the entanglement of Laguerre-Gaussian modes,” Phil. Trans. R. Soc. A. 375(2087), 20150442 (2017). [CrossRef]
8. R. Fickler, G. Campbell, B. Buchler, et al., “Quantum entanglement of angular momentum states with quantum numbers up to 10,010,” Proc. Natl. Acad. Sci. U.S.A. 113(48), 13642–13647 (2016). [CrossRef]
9. L. Yan, P. Gregg, E. Karimi, et al., “Q-plate enabled spectrally diverse orbital-angular-momentum conversion for stimulated emission depletion microscopy,” Optica 2(10), 900–903 (2015). [CrossRef]
10. T. A. Klar, E. Engel, and S. W. Hell, “Breaking Abbe’s diffraction resolution limit in fluorescence microscopy with stimulated emission depletion beams of various shapes,” Phys. Rev. E 64(6), 537 (2001). [CrossRef]
11. M. P. J. Lavery, F. C. Speirits, S. M. Barnett, et al., “Detection of a spinning object using light's orbital angular momentum,” Science 341(6145), 6145 (2013). [CrossRef]
12. N. Cvijestic, G. Milione, E. Ip, et al., “Detecting lateral motion using light's orbital angular momentum,” Sci. Rep. 5(1), 15422 (2015). [CrossRef]
13. J. Vieira, R. M. G. M. Trines, E. P. Alves, et al., “Amplification and generation of ultra-intense twisted laser pulses via stimulated Raman scattering,” Nat. Commun. 7(1), 10371 (2016). [CrossRef]
14. N. S. Mallick and T. N. Dey, “Four-wave mixing-based orbital angular momentum translation,” J. Opt. Soc. Am. B 37(6), 1857–1864 (2020). [CrossRef]
15. P. Gregg, P. Kristensen, A. Rubano, et al., “Enhanced spin orbit interaction of light in highly confining optical fibers for mode division multiplexing,” Nat. Commun. 10(1), 4707 (2019). [CrossRef]
16. S. K. Dacha, W. Zhu, A. Agrawal, et al., “Nonlinear rotation of spin-orbit coupled states in hollow ring-core fibers,” Opt. Express 30(11), 18481–18495 (2022). [CrossRef]
17. K. Rottwitt, J. G. Koefoed, K. Ingerslev, et al., “Inter-modal Raman amplification of OAM fiber modes,” APL Photonics 4(3), 030802 (2019). [CrossRef]
18. X. Liu, E. N. Christensen, K. Rottwitt, et al., “Nonlinear four-wave mixing with enhanced diversity and selectivity via spin and orbital angular momentum conservation,” APL Photonics 5(1), 010802 (2020). [CrossRef]
19. S. K. Dacha, H. F. Elder, W. Zhu, et al., “Observation of conjugate OAM mode generation via modulational instability in a hollow ring-core fiber,” in Optica Nonlinear Optics Topical Meeting 2023, Honolulu (2023). [CrossRef]
20. P. Horak and F. Poletti, “Multimode nonlinear fibre optics: theory and applications,” Recent Progress in Optical Fiber Research, 1–24, (2012).
21. K. Krupa, A. Tonello, B. M. Shalaby, et al., “Spatial beam self-cleaning in multimode fibres,” Nature Photon. 11(4), 237–241 (2017). [CrossRef]
22. K. Krupa, A. Tonello, A. Barthélémy, et al., “Multimode nonlinear fiber optics, a spatiotemporal avenue,” APL Photonics 4(11), 110901 (2019). [CrossRef]
23. L. G. Wright, Z. Liu, D. A. Nolan, et al., “Self-organized instability in graded-index multimode fibres,” Nature Photon. 10(12), 771–776 (2016). [CrossRef]
24. G. Guerra, M. Lonardi, A. Galtarossa, et al., “Analysis of modal coupling due to birefringence and ellipticity in strongly guiding ring-core OAM fibers,” Opt. Express 27(6), 8308–8326 (2019). [CrossRef]
25. R. Bhandari, “Nature of the orbital angular momentum (OAM) fields in a multilayered fiber,” OSA Continuum 4(6), 1859–1874 (2021). [CrossRef]
26. F. Poletti and P. Horak, “Description of ultrashort pulse propagation in multimode optical fibers,” J. Opt. Soc. Am. B 25(10), 1645–1654 (2008). [CrossRef]
27. C. Antonelli, M. Shtaif, and A. Mecozzi, “Modeling of nonlinear propagation in space-division multiplexed fiber-optic transmission,” J. Lightwave Technol. 34(1), 36–54 (2016). [CrossRef]
28. G. Agrawal, Nonlinear Fiber Optics, 5th ed. (Academic, 2013).
29. A. Kudlinski, A. Bendahmane, D. Labat, et al., “Simultaneous scalar and cross-phase modulation instabilities in highly birefringent photonic crystal fiber,” Opt. Express 21(7), 8437–8443 (2013). [CrossRef]
30. L. Li, J. Leng, P. Zhou, et al., “Modulation instability induced by intermodal cross-phase modulation in step-index multimode fiber,” Appl. Opt. 58(16), 4283–4287 (2019). [CrossRef]