CC BY
26
INTEGRATED FIBER-BASED TRANSVERSE MODE CONVERTER
A.V. Gavrilov12, V.S. Pavelyev12 1 Image Processing Systems Institute of RAS - Branch of the FSRC "Crystallography and Photonics " RAS, Samara, Russia,
2 Samara National Research University, Samara, Russia
Abstract
A transverse mode converter based on a binary microrelief implemented directly on the end-face of a few-mode fiber was numerically investigated. The results of numerical simulation demonstrated the converter to form LP-11 and LP-21 modes with high efficiency, providing a more-than 92 % mode purity. Transformations of modes excited by a fiber microbending were also numerically investigated. The excited beams were shown to save their mode purity even in a strong bending as the arising parasitical modes were mostly unguided by the fiber. The resulting beam power and mode content were also demonstrated to depend on the beam and bending mutual orientation for beams with strong rotational symmetry.
Keywords: transverse mode, binary microrelief, few-mode fiber, diffractive optical element.
Citation: Gavrilov AV, Pavelyev VS. Integrated fiber-based transverse mode converter. Computer Optics 2017; 41(4): 510-514. DOI: 10.18287/2412-6179-2017-41-4-510-514.
Acknowledgements: This work was partly supported by the Ministry of Education and Science of the Russian Federation (project 16.7894.2017/6.7).
Introduction
The highly efficient fiber transverse mode conversion is an actual task for many practical applications [1 - 6]. In [2, 3] it was suggested to use diffractive optical elements (DOEs) for mode conversion and excitation in graded-index fiber. The use of DOEs for optical waveguide radiation control is considered also in [7, 8].
Later the same approach was used for excitation of selected mode in a few-mode step-index fiber [4, 9]. However, use of a separated DOE (phase or amplitude one) with transmission function calculated via amplitude-phase coding methods [2, 3] for selected fiber mode forming leads to energy losses due to appearing of parasitical diffractive orders. Use of a separated phase binary DOE [4, 9] does not allow forming selected modal distribution with required purity due to neglecting mode amplitude distribution.
Another disadvantage of fiber mode excitation by DOEs lies in necessity of DOE adjustments regarding fiber optical axis [2 - 4]. In [10] it was suggested to control the angle of incidence of launching beam for selective excitation of coupled modes in a homogeneous coupled multi-core fiber. The approach [10] also requires complicated adjustment procedure.
The realization of micro relief directly on the end face of fiber is considered, for example, in [11]. The present paper is devoted to investigation of integrated fiber-based mode converter based on binary diffractive microrelief implemented directly on an end-face of few-mode stepindex fiber. The few-mode fiber acts as a "filter" of un-guided modes and therefore there is no need to use amplitude-phase coding [2, 3] to increase purity of the mode.
1. Investigation of mode excitation using binary on-fiber microrelief The height function h(x, y) for binary microrelief can be obtained from the diffractive element phase function f(x, y) using the following formula [3]:
f( x, y )1
h ( x, y ) =
2p(ncore -1)
(1)
where X is the wavelength, ncore is the refractive index of fiber core. The formula is applicable if refractive index of the cladding does not differ significantly from the one of the core.
To check these assumptions and to explore beam propagation in fiber a numerical simulations using beam propagation method [12, 13] was made. Parameters of Corning Glass SMF-28 fiber were used (ncore = 1.4619, ncladding = 1.457, core diameter d = 8.3 mm). For the wavelength l = 0.63 mm this fiber acts as a few-mode one. To cover fiber core the calculation window was a square with 16.6^m side, discretization step was 0.2l.
The selective excitation of modes LP-11 and LP-21 was investigated. Their excitation requires phase functions, which are represented in Fig. 1a, b (white stands for 0 and black stands for p). Using Eq. (1) one can obtain corresponding microreliefs, which are schematically shown in Fig. 1 c, d.
Fig. 1. LP-11 and LP-21 modes excitation. Phase functions (a) and (b) were implemented as microreliefs (c) and (d)
In the first simulation fiber end was illuminated by Gaussian beam (with its waist conformed to the LP-01
mode of the fiber). Propagation was simulated for the distance of 5 mm. Resulting beam amplitude distributions and beam phase distributions are shown in Fig. 2.
y ten) y (m)
~0 x (\un) jj -8.3 0 x (|xm) Fig. 2. LP-01, LP-11 and LP-21 modes excitation with Gaussian beam. Cross-section amplitude distributions (a), (c) and (e) and phase distributions (b), (d) and f) are formed using no microrelief, and microreliefs shown in Fig. 1 (c) and (d), respectively
The amplitude and phase distributions looked promising, but it is essential to know how fast beam becomes stable and how much energy is carried by the beam and the excited mode. Fig. 3 shows dependencies of beam power in the fiber core on propagation distance from fiber's end for modes LP-11 and LP-21 (obviously, there is no reason to investigate such a dependence for the LP-01 mode in case of excitation using conformed Gaussian beam).
Relative power 1.0
0.8
0.6
0.4
0.2
----LP-11
-LP-21
0
0.2
0.<
0.6 0.8 1.0 Propagation distance (mm)
Fig. 3. Beam power in the fiber core depending on propagation distance for LP-11 (upper line) and LP-21 (lower line)
Stabilization distances were found to be approximately 0.45 mm for LP-11 and 0.25 mm for LP-21. Beam powers were approximately 73 % and 55 % of initial
powers, and about 92.5 % and 93.5 % of power were carried by the desired modes, respectively.
So, excitation effectiveness was sufficient and most of the beam power was carried by desired modes. Moreover, incoming beam was transformed into the mode on a very short distance (less than 1 mm). It can be explained by the following. Right after diffractive microrelief a beam is a combination of guided fiber modes (with some powers) and some unguided modes, which leave the fiber. Consequently, fiber fragment located behind a microrelief acts as an effective cut-off filter.
The second simulation was made to investigate modes excitation using different source of light. Semiconductor diodes are one of popular sources of coherent beams. They produce beams with complicated enough wave-fronts, but using collimators one can obtain a beam with a wavefront close enough to a planar wave. Also, as core diameter of the fiber in study is small, one can assume that that in fiber core region a beam distribution is very close to a constant one. So, the following results represent illumination of fiber end by a planar wave (with other simulation parameters being the same as in the previous simulation). Resulting beam amplitude distributions and beam phase distributions are shown in Fig. 4.
y (m) y ten)
16.6
__ -16.6_
b) ~16-6 y ten) i6.6
0 x ten)
-16.6 ) -16.6
y (m) i6.6
__ -16.61__
f)
Fig. 4. LP-01, LP-11 and LP-21 modes excitation with Gaussian beam. Cross-section amplitude distributions (a), (c) and (e) and phase distributions (b), (d) and f) are formed using no microrelief, and microreliefs shown in Fig. 1 (c) and (d), respectively
Obviously, results differ from previous ones (especially it can be seen for the LP-01 mode). In the same way it is essential to know how fast modes stabilize and what percent of power is carried by desired modes. Fig. 5
shows dependencies of beam power in the fiber core on propagation distance from fiber's end for three modes. Relative power
1.0 0.8 0.6 0.4 0.2
0
1
3 4 5 6 7 8 9 10 Propagation distance {mm)
Fig. 5. Beam power in the fiber core depending on propagation distance forLP-01 (upper line), LP-11 (middle line) and LP-21 (lower line)
The first significant difference is that power outcome is much lower in comparison to the case of illumination by Gaussian beam. Beam powers were approximately 12 %, 8 % and 8 % of initial powers. But it is explained easily by initial distribution, which carries most of the power in the cladding area, so this part of the beam just leaves the fiber.
The second important result is that the distance, required for beam stabilization, is much longer (about 4 mm against 0.5 mm in the previous simulation). The reason for this is that planar wave creates more parasitical unguided modes and also transfers some power to cladding, so longer distance is required to filter all this unguided powers out.
The third essential point is that mode content is obviously different even for stabilized beams (see, e.g. Fig. 4a). About 49.3 %, 71.8 % and 93.8 % of power were carried by the LP-01, LP-11 and LP-21 modes, respectively.
So, planar waves also can be used to excite higher modes, but stabilization distance is longer and the higher is the mode the better it can be excited using a planar wave. Fiber fragment located behind a microrelief still acts as an effective cut-off filter, but there is more to cut-off.
2. Investigation of modes transformation in fiber bending
It is known that making distortions to a fiber affects the guided beam power and its mode content (for multimode fibers) [2, 3]. From this point it is interesting to know how higher modes survive and preserve their purity in such distortions. Here one simple but revealing distortion, namely fiber bending, was studied.
For demonstration the case was simulated, where fiber lies on two supports with 1 mm between them and the bending occurs due to a force applied to the fiber straight between the supports (yet we did not take photoelasticity effects into account).
When a beam consists of the main mode mostly, it is obvious that some power will leave the fiber and some power will transfer to higher modes (see Fig. 6). But if a beam itself is a higher mode, an outcome is more complicated.
When the LP-01 mode propagates through a fiber bending it has no sense to take the mutual orientation of the mode and bending direction into account (as the mode has radial symmetry). But as the LP-11 mode has a strong
rotational symmetry, its mutual orientation with bending direction can be essential. Two contrary orientations were simulated. In the first one (orientation one) the mode symmetry axis lied in the same plane as the bended fiber, in the second one (orientation two) the mode symmetry axis was orthogonal to the fiber plane.
21.65' 0-21.65 \
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 Fig. 6. Amplitude distribution in fiber section for LP-01 mode propagating through fiber microbending (with bending depth of 70 /Mn)
Deeper investigation of beam profiles during their propagation showed that in orientation one both spots in beam propagated through bending with power loss, but saved the beam structure (see Fig. 7). y (\un) y {wn)
7.5
f)
Fig. 7. Amplitude distribution in cross-section (orthogonal to the fiber axis) for the LP-11 mode in orientation one with bending depth of 74 /mm for different propagation distances: 0 mm (a), 0.2 mm (b), 0.4 mm (c), 0.6 mm (d), 0.8 mm (e) and 1 mm (f)
In contrast, in orientation two the outer spot (regarding the fiber curvature center) leaved the fiber core while the inner one preserved and mutated almost to the main mode with loss of power (see Fig. 8).
Obviously, these two orientations should show different characteristics in propagation through a fiber bending.
The LP-21 mode has a rotational symmetry also, but due to its amplitude distribution structure the effect of
different fiber bending and mode mutual orientations should not be so significant. For the sake of completeness it is necessary to add that simulation for LP-21 mode was made for the case when one symmetry axis lied in the fiber plane and the other was orthogonal to the plane. y(\im) 7.5
Fig. 8. Amplitude distribution in cross-section (orthogonal to
the fiber axis) for the LP-11 mode in orientation two with bending depth of 74 jum for different propagation distances: 0 mm (a), 0.2 mm (b), 0.4 mm (c), 0.6 mm (d), 0.8 mm (e) and 1 mm (f)
Fig. 9 shows dependencies of beam and mode powers remained in the fiber core on the bending depth.
Relative power 1.01—*-
60 75 90 Bending depth, [im
Fig. 9. Beam power in the fiber core depending on microbenging depth for LP-01 (the first from the top), two orientations of LP-11 (the second and the third from the top) and LP-21 (the first from the bottom). Thick lines represent full beam powers and thin lines represent mode powers
It can be seen that the main mode LP-01 is the least sensitive while LP-11 is more sensitive and LP-21 is even more sensitive. So, modes with higher numbers have higher sensibility. It can be explained by the fact that the higher is the mode number, the greater part of mode
power propagates closer to fiber cladding. Consequently, bending affects them more than lower modes.
Also it can be seen that two orientations of LP-11 mode showed significantly different sensitivity to bending. Obviously, other orientations will give results between the obtained ones.
Another effect is revealed by gaps between beam power and mode power graphs. The gap for LP-01 increases as some power transfers to higher modes. But gaps for higher modes decrease because there are less higher modes to transfer power into. Also the product of higher mode propagation through a bending is scarcely to be similar to lower modes, especially the main one. The exception is the situation when bending is so strong, so most of the beam power leaves the fiber and the leftovers can be considered as a mode mixture noise (see graphs for LP-11 orientation two and LP-21). So, here full beam power is close to the power of the mode.
Conclusion
In summary, transverse mode converter based on binary microrelief implemented directly on an end-face of a few-mode fiber was investigated. Results of numerical simulation demonstrate high efficiency of LP-11 and LP-21 mode forming by the converter. Difference between illuminating beam amplitude and amplitude of mode to be formed results mainly in excitation of un-guided modes and, consequently, amplitude-phase coding is not required to obtain high mode purity (more than 92 %) of a beam in case of illumination by a Gaussian beam. In case of illumination by a planar wave effectiveness is much lower, but mode purity is higher the higher is an excited mode.
Transformations of excited modes by fiber microbending were numerically investigated. It was shown that excited beams save their mode purity even in a strong bending, as arising parasitical modes are mostly unguided by the fiber. It was also demonstrated that resulting beam power and mode content depend on beam and bending mutual orientation for beams with strong rotational symmetry.
References
[1] Kersey AD. A review of recent developments in fiber optic sensor technology. Optical Fiber Technology 1996; 2(3): 291-317. DOI: 10.1006/ofte.1996.0036.
[2] Soifer V, Golub M. Laser beam mode selection by computer generated holograms. Boca Raton: CRC Press; 1994. ISBN 0-8493-2476-9.
[3] Soifer VA, ed. Methods for Computer Design of Diffrac-tive Optical Elements. New York: John Wiley & Sons; 2002. ISBN: 978-0-471-09533-0.
[4] Karpeev SV, Pavelyev VS, Khonina SN, Kazanskiy NL, Gavrilov AV, Eropolov VA. Fibre sensors based on transverse mode selection. J Mod Opt 2007; 54(6): 833-844. DOI: 10.1080/09500340601066125.
[5] Carpenter J, Wilkinson TD. Characterization of multimode fiber by selective mode excitation. Journal of Lightwave Technology 2012; 30(10): 1386-1392. DOI: 10.1109/JLT.2012.2189756.
[6] Soifer VA, Karpeev SV, Pavelyev VS, Kazanskiy NL, Gavrilov AV. Fibre-optic device for measuring transverse
deformation vector [In Russian]. Pat RF of Invent N 2386105 C1 of April 10, 2010, Russian Bull of Inventions N10, 2010.
[7] Alferov SV, Khonina SN, Karpeev SV. Study of polarization properties of fiber-optics probes with use of a binary phase plate. J Opt Soc Am A 2014; 31(4): 802-807. DOI: 10.1364/JOSAA.31.000802.
[8] Karpeev SV, Khonina SN. Experimental excitation and detection of angular harmonics in a step-index optical fiber. Opt Mem Neural Networks 2007; 16(4): 295-300. DOI: 10.3103/S1060992X07040133.
[9] Gavrilov AV, Karpeev SV, Kazanskiy NL, Pavelyev VS, Duparre M, Luedge B, Schroeter S. Selective excitation of step-index fiber modes. Proc SPIE 2006; 6605: 660508. DOI: 10.1117/12.728461.
[10] Kokubun Y, Komo T, Takenaga K, Tanigawa S, Matsuo S. Selective mode excitation and discrimination of four-core homogeneous coupled multi-core fiber. Opt Express 2011; 19(26): B905-B914. DOI: 0.1364/OE.19.00B905.
[11] Pavelyev VS, Moiseev OYu, Volkov AV, Eropolov VA, Dmitriev SV, Karpeev SV, Artyushenko VG, Kashin VV. Realization and characterization of diffraction microrelief on the end faces of silver-halide waveguide. Proc SPIE 2008; 6994: 69940Q. DOI: 10.1117/12.780695.
[12] Chiou Y-P, Chang H. Efficient beam-propagation method based on Pad'e approximants in the propagation direction. Opt Lett 1997; 22(13): 949-951. DOI: 10.1364/OL.22.000949.
[13] Ma Ch, van Keuren E. A simple three dimensional wide-angle beam propagation method. Opt Express 2006; 14(11): 4668-4674. DOI: 10.1364/OE.14.004668.
Authors' information
Andrey Vadimovich Gavrilov (b. 1982), graduated with honors (2005) from the S.P. Korolyov Samara State Aerospace University (SSAU, presently, Samara National Research University). He has a candidate degree in Physical and Mathematical Science (2009). At present he is an Associate Professor at Technical Cybernetics department of Samara University. The area of interests includes simulation of light propagation in waveguides as well as parallel and distributed computing. He's list of publications contains more than 30 scientific papers. E-mail: andrey.w.gavrilov@gmail.com .
Vladimir Sergeevich Pavelyev, Doctor of Physical and Mathematical Sciences; Head of Nanoengineering department of the Samara National Research University, chief researcher of the Image Processing Systems Institute of the RAS - Branch of the FSRC "Crystallography and Photonics" RAS. Research interests: diffractive microoptics, optical waveguides, nanophotonics, optical metamaterials synthesis, technologies for forming of micro- and nanostructures. Email: pavelyev10@mail. ru .
Code of State Categories Scientific and Technical Information (in Russian - GRNTI): 29.33.39.
Received April 02, 2017. The final version - August 17, 2017.
