Abstract
Using the power flow equation (PFE), this article investigates mode coupling in step-index (SI) multimode (MM) polymer optical fiber (POF). This equation’s coupling coefficient was initially fine-tuned so that it could appropriately reconstruct previously recorded far-field (FF) power distributions. The equilibrium mode distribution (EMD) and steady-state distribution (SSD) in the SI MM POF were found to be obtained at lengths Lc = 15 m and zs = 41 m, respectively. These lengths are substantially shorter than their glass optical fiber counterparts. Such characterization of the investigated POF can be used in its employment as a part of the communication or sensory system. Namely, the POF’s bandwidth is inverse linear function of fiber length (z−1) below the coupling length Lc. However, it has a z−1/2 dependence beyond this equilibrium length. Thus, the shorter the coupling length Lc, the sooner transition to the regime of slower bandwidth decrease occurs. It is also important to be able to determine a modal distribution at a certain length of the POF employed as a part of optical fiber sensory system.
Introduction
Silica optical fiber has a lot of advantages such as a low loss, lightweight, and high bandwidth [1–4]. Although silica optical fiber has a lot of advantages, is not appropriate for short-distance applications such as automotive [5], local area networks [6], and visible light communication (VLC) applications [7]. On the other hand, POFs are made of poly (methylmethacrylate) (PMMA) [8, 9], polycarbonate (PC) [10, 11], ZEONEX ® [12, 13], TOPAS ® [14, 15], poly dimethyl siloxane (PDMS) [16, 17] or hydrogel [18, 19].
Apart from optical fiber material, optical fibers usually have step-index or graded-index refractive index distribution, and can operate in a single-mode or multimode regime. There are several typical commercial POFs, such as Mitsubishi Rayon’s Eska Extra EH 4001 SI MM POF [20] and CYTOP® graded index (GI) POF [21], both of them could transmit a large number of modes, and have been used for short-range communications and sensing area in fiber optic sensory systems [22]. Additionally, micro-structured polymer optical fiber which was first produced by Argyros’ group, is one very specific type of POF [23]. However, due to the transmission loss and but-coupling with commercial silica optical fiber, only one Start-up Company supplied the commercial product [24]. PMMA SI MM POF is still the most common one with lots of commercial products on the market.
Mode coupling has a significant impact on the transmission performance of PMMA SI MM POF. The transfer of power between nearby modes is represented by such coupling. It causes the launched light’s angular power distribution to be disrupted [25, 26]. It is caused by the optical fiber’s intrinsic perturbation effects (for example, variations of the refractive index distribution and microscopic bends). The angular power distribution is predicted to be influenced by the launch conditions and mode coupling characteristics. Thus, a Gaussian beam launched at an angle θ0 in respect to the optical fiber axis, at the output end of a short piece of optical fiber, is seen as a sharp ring pattern (the ring diameter depends on θ0). Mode coupling, on the other hand, causes such a ring pattern to be distorted in longer fibers and finally transformed into disk. The coupling length Lc marks the fiber length at which the ring pattern of the highest order guiding mode evolved into disk, indicating that an EMD is established. Coupling can fully complete at optical fiber length zs (zs> Lc), which is referred to as an SSD.
The light transmission performance of PMMA SI MM POF is investigated in this work. The modal power diffusion model is used to obtain the coupling coefficient D. In this way, we computed the Lc required to accomplish the EMD and the length zs to achieve the SSD. By such a characterization of mode coupling process in an optical fiber one can predict at which length z = Lc one can expect that bandwidth decrease with length would start to decelerate, as explained in more details in Results and Discussion. It is also important to know the state of mode coupling in an optical fiber employed as a part of optical fiber sensory system, especially in terms of the modal distribution at certain fiber length.
Power Flow Equation
The Gloge’s PFE has the following form [27]:where P(θ, z) is the angular power distribution, θ is the angle of propagation, z is the distance of the propagation, D is the coupling coefficient (assumed constant [27, 28]) and is the modal attenuation. Eq. 1 can be reduced to [29]:Steady-state solution of Eq. 2 is given as [28]:where J0 is the Bessel function of the first and zero-order, γ0 [m−1] = 2.4052D/θc2 is the attenuation coefficient. In this work, solved Eq. 2 using the explicit finite difference method [29].
In our earlier published work, we proposed a method which enables that the coupling coefficient D can be obtained from just two output angular power distributions P(θ, z) in the case of centrally launched beam [30]. As an alternative, in this work we propose a method for calculating the coupling coefficient D on the basis of the measured FFPs p(x, z), illustrated in Figure 1.
FIGURE 1

Schematics of measuring FF power distributions.
Results and Discussion
We used the PFE (2) to calculate the lengths Lc and zs for the POF that Ribeiro et al. [31] studied experimentally. This POF (Mitsubishi Rayon’s Eska Extra EH 4001) had NA = 0.47, the inner critical angle = 18o ( = 27.4o measured in air), core refractive index ncore = 1.4897, numerical aperture NA = 0.46 and fiber diameter d = 1 mm. Ribeiro et al. focused a He-Ne laser beam at 633 nm onto the input end of the fiber in their experiment (Figure 1). They used a CCD to detect the FF patterns of the fiber output.
The coupling coefficient D is obtained by its fine-tuning in the PFE, thus recreating the Ribeiro et al.’s measured FF patterns in this fiber for lengths of z = 1 m and z = 19 m (Figure 2). This necessitated numerically calculating the PFE for various values of D. The value of D = 9.2 × 10−4 rad2/m enabled the best fit between the calculated and measured output power distributions.
FIGURE 2

The measured FF pattern from (A) ∼1 m and (B) ∼19 m PMMA SI MM POF lengths. The light was launched at θ0 = 0o [31].
A Gaussian launch beam with (FWHM)z=0 = 7o is used in the computations. Figure 3 shows the output power distribution for a z = 1 and 19 m length POF with an input angle of θ0 = 0o obtained as numerical solution of Eq. 2. The distributions P (x, z = 1 and 19 m, L = 110 mm) are generated from the distributions P (θ, z = 1 and 19 m, L = 110 mm), where x = L⋅tgθ and L is the receiving distance.
FIGURE 3

Normalized output power distribution P (x, z = 1 and 19 m, L = 110 mm) at the end of (A) 1 m and (B) 19 m long PMMA SI MM POF, obtained by solving the PFE for Gaussian launch distribution with input angle θ0 = 0°, obtained from the P (x, z = 1 and 19 m, L = 110 mm), for x = L⋅tgθ, where L is the receiving distance.
In Figure 4, our numerical results for the normalized output angular power distribution for input angles θ0 = 0, 5, 10, and 15o are shown. At short fiber lengths, a significant mode coupling is observed for low order modes, as can be seen in Figure 4A. The EMD is achieved at length z = Lc = 15 m (Figure 4C), while the SSD is established at zs = 41 m in Figure 4D.
FIGURE 4

The normalized angular power distribution at the end of PMMA SI MM POF, obtained by solving the PFE for four Gaussian launch distributions with input angles θ0 = 0o (—), θ0 = 5o (− −), θ0 = 10o (•••) and θ0 = 15o (−•−), with (FWHM)z=0 = 7o (g represent the analytical SSD).
The coupling coefficient D for the PMMA SI MM POF evaluated in this work is similar to those which we obtained for other investigated POFs (∼10−4 rad2/m), so the characteristic lengths for coupling, Lc and zs, are similar (Lc is between ≃15–35 m and zs is between ≃40 and 100 m [25, 29, 30]). We have previously reported that glass optical fibers show the weakest strength of mode coupling (D ∼ 10−7 to 10−6 rad2/m), so their SSD lengths zs are between ≃ 1–10 km [32].
It is worth noting that the length dependence of the bandwidth of POF is determined by mode coupling behavior. The bandwidth is inverse linear function of fiber length (z−1) below the coupling length Lc. However, it has a z−1/2 dependence beyond this equilibrium length. Thus, the shorter the coupling length Lc, the sooner transition to the regime of slower bandwidth decrease occurs [26]. It is obvious that mode coupling has a beneficial influence on bandwidth. Therefore, it is of great interest to characterize mode coupling in an optical fiber in order to predict its transmission characteristics, especially to determine at which coupling length Lc one can expect that bandwidth would start to improve. It is also important to know the state of mode coupling in an optical fiber employed as a part of optical fiber sensory system, especially in terms of the modal distribution at certain fiber length.
Conclusion
We investigate a mode coupling along a PMMA SI MM POF previously investigated experimentally by Ribeiro et al. [31]. To appropriately recreate the measured FF patterns reported before, the coupling coefficient D in the PFE is tweaked. As a result, the lengths Lc and zs that characterize the coupling process are obtained. Such characterization of the investigated PMMA SI MM POF can be used in its employment as a part of communication or sensory system. In practice, by characterization of mode coupling in an optical fiber one can predict its transmission characteristics, especially to determine length-dependent bandwidth behavior. Since POF’s bandwidth is inverse linear function of fiber length (z−1) below the coupling length Lc, while it has a z−1/2 dependence beyond this equilibrium length, it is obvious that the shorter the coupling length Lc leads to the sooner transition to the regime of slower bandwidth decrease. It is also important to be able to determine a modal distribution at a certain length of the fiber employed as a part of optical fiber sensory system.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
SS: methodology, conceptualization, formal analysis, validation, writing—original draft. RM: conceptualization, formal analysis, funding acquisition, writing—review and editing. AD: formal analysis, validation, writing—review and editing. BD: validation, writing—review and editing. AS: validation, writing—review and editing.
Funding
This research was funded by the National Natural Science Foundation of China (62003046, 6211101138); the Strategic Research Grant of City University of Hong Kong (Project No. CityU 7004600); Serbian Ministry of Education, Science and Technological Development (Agreement No. 451-03-68/2020-14/200122); Science Fund of the Republic Serbia (Agreement No. CTPCF-6379382); Guangdong Basic and Applied Basic Research Foundation (2021A1515011997); The Innovation Team Project of Guangdong Provincial Department of Education (2021KCXTD014); Special project in the key field of Guangdong Provincial Department of Education (2021ZDZX1050).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1.
MinRLiuZPereiraLYangCSuiQMarquesC. Optical Fiber Sensing for marine Environment and marine Structural Health Monitoring: A Review. Opt Laser Techn (2021) 140:107082. 10.1016/j.optlastec.2021.107082
2.
XiangMFuSXuOLiJPengDGaoZet alAdvanced DSP Enabled C-Band 112 Gbit/s/λ PAM-4 Transmissions with Severe Bandwidth-Constraint. J Lightwave Technol (2022) 40(4):987–96. 10.1109/jlt.2021.3125336
3.
MaedaHSaitoKSasaiTHamaokaFKawaharaHSekiTet alReal-time 400 Gbps/carrier WDM Transmission over 2,000 Km of Field-Installed G654E Fiber. Opt Express (2020) 28(2):1640–6. 10.1364/oe.383471
4.
LiYFanHZhangLWangKWuGLiuZ. A bend-tolerant BOTDR Distributed Fiber Sensor. Opt Commun (2022) 514:128110. 10.1016/j.optcom.2022.128110
5.
POF automotive application. Automotive (2022). Available at: https://www.kdpof.com/automotive/(Access on 04 02, 2022).
6.
ForniFShiYTranN-CVan Den BoomHTangdionggaEKoonenT. Multi-Format Wired and Wireless Signals over Large-Core Plastic Fibers for in-home Network. J Lightwave Technol (2018) 36:1. 10.1109/JLT.2018.2839029
7.
ApoloJAOrtegaBAlmenarVAlmenarV. Hybrid POF/VLC Links Based on a Single LED for Indoor Communications. Photonics (2021) 8(7):254. 10.3390/photonics8070254
8.
LiuLZhengJDengSYuanLTengC. Parallel Polished Plastic Optical Fiber-Based SPR Sensor for Simultaneous Measurement of RI and Temperature. IEEE Trans Instrum Meas (2021) 70:1–8. 10.1109/TIM.2021.3072136
9.
HuXWoyessaGKinetDJantingJNielsenKBangOet alBDK-doped Core Microstructured PMMA Optical Fiber for Effective Bragg Grating Photo-Inscription. Opt Lett (2017) 42(11):2209–12. 10.1364/ol.42.002209
10.
ZubelMGFasanoAWoyessaGTMinRLeal-JuniorAGTheodosiouAet alBragg Gratings Inscribed in Solid-Core Microstructured Single-Mode Polymer Optical Fiber Drawn from a 3D-Printed Polycarbonate Preform. IEEE Sensors J (2020) 20(21):12744–57. 10.1109/JSEN.2020.3003469
11.
MoslanMSOthmanMHDSamavatiASalimMAMRahmanMAIsmailAFet alFabrication of Polycarbonate Polymer Optical Fibre Core via Extrusion Method: The Role of Temperature Gradient and Collector Speed on its Characteristics. Opt Fiber Techn (2020) 55:102162. 10.1016/j.yofte.2020.102162
12.
WoyessaGRasmussenHKBangO. Zeonex - a Route towards Low Loss Humidity Insensitive Single-Mode Step-index Polymer Optical Fibre. Opt Fiber Techn (2020) 57:102231. 10.1016/j.yofte.2020.102231
13.
DashJNChengXGunawardenaDSTamH-Y. Rectangular Single-Mode Polymer Optical Fiber for Femtosecond Laser Inscription of FBGs. Photon Res (2021) 9(10):1931–8. 10.1364/prj.434252
14.
MarquesCAFMinRJuniorALAntunesPFasanoAWoyessaGet alFast and Stable Gratings Inscription in POFs Made of Different Materials with Pulsed 248 Nm KrF Laser. Opt Express (2018) 26(2):2013–22. 10.1364/oe.26.002013
15.
MarkosCStefaniANielsenKRasmussenHKYuanWBangO. High-T_g TOPAS Microstructured Polymer Optical Fiber for Fiber Bragg Grating Strain Sensing at 110 Degrees. Opt Express (2013) 21(4):4758–65. 10.1364/oe.21.004758
16.
Leal-JuniorAGuoJMinRFernandesAJFrizeraAMarquesC. Photonic Smart Bandage for Wound Healing Assessment. Photon Res (2021) 9(3):272–80. 10.1364/prj.410168
17.
ShentuZKangJZhuZWangLGuoYXuTet alNo-core PDMS Fiber for Large-Scale Strain and Concentration Measurement. Opt Fiber Techn (2021) 63:102531. 10.1016/j.yofte.2021.102531
18.
ElsherifMHassanMUYetisenAKButtH. Hydrogel Optical Fibers for Continuous Glucose Monitoring. Biosens Bioelectron (2019) 137(15):25–32. 10.1016/j.bios.2019.05.002
19.
GuoJLiuXJiangNYetisenAKYukHYangCet alHighly Stretchable, Strain Sensing Hydrogel Optical Fibers. Adv Mater (2016) 28(46):10244–9. 10.1002/adma.201603160
20.
Mitsubishi POF SH4001. Ibsstore (2022). Available at: https://www.ibselectronics.com/ibsstore/sh4001-mitsubishi-rayon-eska-plastic-optical-fiber-cable.html (Access on 04 02, 2022).
21.
WoyessaGTheodosiouAMarkosCKalliKBangO. Single Peak Fiber Bragg Grating Sensors in Tapered Multimode Polymer Optical Fibers. J Lightwave Technol (2021) 39(21):6934–41. 10.1109/JLT.2021.3103284
22.
HeRTengCKumarSMarquesCMinR. Polymer Optical Fiber Liquid Level Sensor: A Review. IEEE Sensors J (2022) 22(2):1081–91. 10.1109/jsen.2021.3132098
23.
van EijkelenborgMLargeMArgyrosAZagariJManosSIssaNet alMicrostructured Polymer Optical Fibre. Opt Express (2001) 9(7):319–27. 10.1364/oe.9.000319
24.
Shute Company. Shute (2022). Available at: https://shute.dk/(Access on 04 02, 2022).
25.
SavovićSDjordjevichA. Mode Coupling in Strained and Unstrained Step-index Plastic Optical Fibers. Appl Opt (2006) 45:6775–80. 10.1364/AO.45.006775
26.
MateoJLosadaMAZubiaJ. Frequency Response in Step index Plastic Optical Fibers Obtained from the Generalized Power Flow Equation. Opt Express (2009) 17:2850–60. 10.1364/oe.17.002850
27.
GamblingWAPayneDNMatsumuraH. Mode Conversion Coefficients in Optical Fibers. Appl Opt (1975) 14:1538–42. 10.1364/ao.14.001538
28.
RousseauMJeunhommeL. Numerical Solution of the Coupled-Power Equation in Step-Index Optical Fibers. IEEE Trans Microwave Theor Techn. (1977) 25:577–85. 10.1109/TMTT.1977.1129162
29.
DjordjevichASavovicS. Investigation of Mode Coupling in Step index Plastic Optical Fibers Using the Power Flow Equation. IEEE Photon Technol Lett (2000) 12:1489–91. 10.1109/68.887704
30.
SavovićSDjordjevichA. Method for Calculating the Coupling Coefficient in Step-index Optical Fibers. Appl Opt (2007) 46:1477–81. 10.1364/AO.46.001477
31.
RibeiroRSilvaVBalodYBarberoAPGermanoSSantosPd. Profile Scanning Measurements of Far-Field from Plastic Optical Fibre (POF) Passive Devices for Bandwidth Characterization in Datacom Networks. Rio de Janeiro: XXVI Simposio Brasileiro de Telecomunicaciones-BrT´ (2008). p. 08. 10.14209/sbrt.2008.42632
32.
DjordjevichASavovićSTsePWDrljačaBSimovićA. Mode Coupling in Strained and Unstrained Step-index Glass Optical Fibers. Appl Opt (2010) 49:5076–80. 10.1364/AO.49.005076
Summary
Keywords
polymer optical fiber, modal distribution, coupling length, bandwidth, distance dependence, sensory system
Citation
Savović S, Djordjevich A, Drljača B, Simović A and Min R (2022) Calculation of the Coupling Coefficient in Step-Index Multimode Polymer Optical Fibers Based on the Far-Field Measurements. Front. Phys. 10:927907. doi: 10.3389/fphy.2022.927907
Received
25 April 2022
Accepted
09 May 2022
Published
25 May 2022
Volume
10 - 2022
Edited by
Sushank Chaudhary, Quanzhou Institute of Equipment Manufacturing (CAS), China
Reviewed by
Lei Liu, Harbin Engineering University, China
Yuyang Peng, Macau University of Science and Technology, Macao SAR, China
Updates
Copyright
© 2022 Savović, Djordjevich, Drljača, Simović and Min.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rui Min, ruimin@bnu.edu.cn
This article was submitted to Optics and Photonics, a section of the journal Frontiers in Physics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.