Design of Silica Multimode Optical Fibers with Extremely Enlarged Core Diameter for Laser-Based Multi-Gigabit Short-Range Optical Networks †

This work presents an alternative fast and simple method for the design of a refractive index profile of silica multimode optical fibers (MMFs) with extremely enlarged core diameters of up to 100 μm for laser-based multi-gigabit short-range optical networks. We demonstrate some results of 100 μm core MMF graded index profile optimization performed by a proposed solution, which provides a selected mode staff differential mode delay (DMD) reduction over the “O”-band under particular launching conditions. Earlier on, a developed alternative model for a piecewise regular multimode fiber optic link operating in a few-mode regime for the computation of laser-excited optical pulse dynamics during its propagation over an irregular silica graded-index MMF with an extremely large core diameter, is utilized to estimate the potentiality of fiber optic links with the described MMFs. Here, we also present the comparison results of the simulation of 10GBase-LX optical signal transmission over 100 μm core MMFs with conventional and optimized graded-index refractive index profiles.


Introduction
Nowadays, laser-based optical signal transmission techniques over silica optical fibers with core diameters enlarged in comparison with standard telecommunication single-mode fibers are widely used in various applications of high-bit-rate optical networks.This technique, combined with special launching conditions, provides a few-mode regime, where laser-excited optical emission is transferred over large core optical fibers not by total-mode staff, but only by a limited number of mode components.Since the IEEE 802.3z standard was ratified by the Institute of Electrical and Electronics Engineers (IEEE) in 1998 [1,2], this technique began to be widely used for short-range in-premise multi-gigabit networks, such as cable systems of data and computation centers, storage and local area networks, etc.
Since 1998, almost every five years, a new-generation of MMFs of Cat.OM2 + -OM5 (these categories are outlined and declared in the ISO/IEC 11801 standard ratified by the International Organization for Standardization (ISO) and the International Electrotechnical Commission (IEC)) known as "LOMFs" (laser-optimized multimode-fibers) with a core diameter of 50 µm were specifically developed, designed, and ratified for multi-gigabit optical networks with laser-based transceivers and multimode fiber optic links operating in a few-mode regime [2][3][4]6].
Nowadays, this technique has become in demand for on-board and industrial network applications requiring 1 Gbps and greater bit rates [7][8][9], where silica MMFs with enlarged core diameters of up to 100 µm are the most preferable solution.Here, cable systems operate in harsh environmental conditions: high-density mounting spaces, vibrations, dust/suspensions in the air, etc.Therefore, the greater the fiber core diameter, the greater the total network reliability: (a) Optical fibers with a large core diameter of 100 µm (100/125 and 100/140) have almost the same immunity to bending loss as well as specially designed single-mode fibers (SMFs) of ITU-T G.657 Recommendations declared by the International Telecommunication Union-Telecommunication sector (ITU-T) as a solution for indoor "last/first mile" cable systems.(b) Unlike the SMFs of ITU-T G.657 Recommendations, the fiber optic connections of MMFs with a 100 µm core are insensitive to vibrations, and (c) the dust on connector ferrules, due to their great core size.(d) Additionally, the silica core and silica cladding provide lower losses in comparison to the polymer or plastic core/cladding optical fibers.However, the most commercially available MMFs with 100 µm core diameters are targeted only to multimode regimes and are based on multimode light emitting diode (LED) transceivers providing low bit rates of 10-100 Mbps data transmission, while their potentiality for laser-based data transmission in a few-mode regime still requires additional research.
The rest of this paper is structured as follows.Section 2 introduces the features of the DMD effect occurring during laser-excited optical signal transmission over MMFs, presents a brief overview of related works devoted to the design of MMFs with improved bandwidths, and describes in detail the proposed alternative fast and simple method for the optimization of the refractive index profile for MMFs with a reduced DMD, based on an earlier-developed math tool that has been recently modified for the design of 100 µm core MMFs operating under laser-based multi-gigabit short-range optical networks.Section 3 presents an optimized sample of a 100 µm core MMF refractive index profile that provides reduction of DMD over the "O"-band, and some simulation results of 10GBase-LX optical signal transmissions over irregular 100 µm core MMFs with both typical and optimized graded-index refractive index profiles.Section 4 is concerned with discussion that is focused on the analysis of 10GBase-LX channel quality parameters, based on a recalculation of the simulation results, while Section 5 presents our conclusions.

Related Works
Laser-excited optical signal transfer over MMF by a limited number of mode components primarily differs from a multimode regime by the appearance of the DMD effect [1][2][3][4][5][6].DMD is a main issue for short-range multi-gigabit optical networks based on MMFs.Generated by a coherent source, optical emission is launched to an MMF core with a much greater diameter then the laser spot size.This great difference between the fiber core diameter and the laser spot size leads to the excitation of a limited number of particular order modes from the total MMF guided mode staff.Moreover, during propagation in MMF, the amplitudes of these selected mode components transferring the optical signal are non-equal.The description of the amplitude difference is mainly defined by its launching conditions: the typical fiber optic adaptor, radial or angular misalignments, a special matching device, etc.Therefore, the optical pulse at the transmitter end of the fiber optic link contains several modes with non-equal amplitudes and differing mode velocities (e.g., mode delays).As a result, these modes come to the receiver end with the spread of both amplitudes and delays relative each other, which becomes critical at multi-gigabit rates.DMD strongly distorts the pulse form; it splits at the receiver, and this generates a "glove" effect.Moreover, DMD for each particular "laser-MMF" pair may strong differ, due to the launching uncertainty provided by the variations of sources, connections, and fiber parameters.
A detailed overview of the known solutions concerned with the development of MMF with improved bandwidth is represented in published monographs [10].Some of them were designed as LOMFs, so that they are targeted not only for modal dispersion reduction, but also for DMD suppression.Most of the known approaches for LOMFs are based on DMD monitoring during fiber preform manufacturing.As proposed in [11][12][13], the objective function estimates the desired refractive index profile proximity to a profile, which is optimal for the smallest deviation of delays between the total MMF mode staff that satisfies the cut-off condition that corresponds to a multimode regime.In [14][15][16], the object function is the difference between the integral and the local value of the profile grade parameter α i for the preform refractive index profile approximated by a simple power function and some optimal α-parameter that is also the best for a multimode regime.Here, local profile-grade parameter correction is performed by reproducing the optimal for the total mode staff DMD diagram-the DMD value distribution on the corresponding mode orders.It is formed by the corresponding modification of the initial DMD diagram of the basic model of optical fibers under overfilled launching conditions, which is the "worst case" for bandwidth, while centralized launching is well known to be the simplest method of implementation under field conditions for bandwidth improvement of a fiber optic link with MMFs operating in a few-mode regime [1,5,10].
Therefore, we suppose that the following outlined factors should be taken into account during the design of MMFs with reduced DMD, or so-called low differential mode delay fibers (LDMDFs) that are optimized for operation under a few-mode regime: (1) Real and approved commercially available LOMFs of Cat.OM2 + -OM5 structure, geometry, refractive index profile, and parameters of data should be used initially/approximated, but not the first generation MMFs of Cat.OM1 and OM2, as well as the model of optical fibers with the ideal graded refractive index profiles as represented by one or a set of simple power-law α-functions.(2) The development or modification of a fiber structure/refractive index profile has to be based on comparison with the DMD diagram.(3) The laser source type and its parameters, including the generated emission initial transverse mode staff and the launching conditions, should be taken into account.(4) Reduced DMD should be provided, not only at particular wavelengths but over the spectral range (for example, conventional telecommunication "O"-band with a reference wavelength λ = 1300-1310 nm).

Design of an LDMDF-Graded Refractive Index Profile
The proposed solution is based on the design of a specialized form of refractive index profile providing a selected guided mode staff delay equalization in relation to some reference delay value, t REF .The designed LDMDF structure weakly satisfies the guiding optical waveguide approximation.It contains a fused silica core doped by germanium and fluorine, bounded by a pure fused silica outer solid cladding.Here, unlike the known solutions, we utilize a stratification method [17] approach to describe the desired refractive index profile.As a result, the designed weakly guiding optical fiber with an arbitrary axially symmetric refractive index profile is represented in the form of a multilayered optical fiber in the core region.It is considered as a finite set on N layers, where the refractive index value stays constant: and the profile function f(R) describing the refractive index profile by a known relation: is written as follows: where N is the local profile parameter; n k is the refractive index of k layers (k = 0 . . .N); n max is the maximal core refractive index; n N is the outer cladding refractive index; ∆ = n 2 max − n 2 N /2n 2 max is the profile height parameter; R k = r k /a is the normalized radial coordinate of the k-th layer; r k is the radial coordinate of the k-th layer; a is the designed LDMDF core radius.
As a result, the desired refractive index profile form in the LDMDF core region is selected in such a way that it should satisfy the minimization of some objective function F described by the proposed simple formula: where t d(j) is the desired value of delay for the j-th guided mode LP lm(j) computed at the corresponding synthesis iteration; t REF is some reference value of mode delay that is applied for DMD diagram equalization; M is the total number of the mode components transferring a laser source excited few-mode optical signal with a normalized amplitude that is not less than A j > 0.1, and with a core power (known also as an optical confinement factor) that is not less then P (j) co ≥ 0.5.Here, the total number of modes M taken into account depends on the following factors: (1) the designed fiber basic geometry parameters (e.g., core diameter and profile height parameter), (2) the launching conditions, (3) the emission of the initial transverse mode staff at the laser output defined by the type of source-vertical surface-emitting laser (VCSEL) or single-mode laser, e.g., Fabry-Perot laser diode (LD)/distributed feedback laser (DFB-laser), and (4) the prediction of a new guided mode of excitation, with A j > 0.1 during the following optical signal propagation over MMF, due to mode mixing and power diffusion effects provided by real fiber irregularities and its micro-/macro-bends and tensions/stress occurring under fiber optic cable manufacturing, installation, and maintenance.Unlike the known solutions, we propose setting the reference delay t REF from the range of values containing the DMD diagram formed for only M selected guided modes, with particular orders propagating in the new generation of LOMFs of Cat.OM2 + -OM4.The objective function F in Equation ( 4) is minimized by the Nelder-Mead simplex method, whose efficiency was demonstrated in [18,19] concerned with the design of optical fibers.

Extension of Modified Gaussian Approximation
During objective function (4) minimization, a direct problem of the optical fiber with the current refractive index profile should be solved during each iteration.Here, objective function arguments are an array of local parameters h k that completely describe the optical fiber refractive index profile.Therefore, a fast and simple method for the computation of both the fundamental and higher-order mode parameters propagating in MMF is required.We propose utilizing earlier on the developed modification of a Gaussian approximation [20] that is generalized and extended [21] for the evaluation of arbitrary order-guided mode dispersion parameters propagating in a weakly guiding optical fiber, with an arbitrary axially symmetric refractive index profile in the core region, bounded by a single solid outer cladding.This extended modified Gaussian approximation (EMGA) is based on a combination of "classical" Gaussian approximation [22] and a stratification method [17] for the researched MMF complicated graded refractive index profile representation.The proposed approach permits the derivation of an analytical formula for the square core mode parameter U 2 in the form of finite nested sums [21] from a well-known integral variational expression [22]: where where R 0 = ρ 0 /a is the equivalent (as a result of Gaussian approximation) normalized mode field radius (MFR); ρ 0 is the equivalent MFR; a is the MMF core radius; l is the azimuthal mode; m is the radial mode of LP lm orders.V = k 0 an max √ 2∆ is the normalized frequency; k 0 = 2π/λ is the wavenumber; λ is the operating wavelength; b (l,m) p is the expansion factor of the Laguerre polynomial representation in the form of a finite power series [23,24]: The characteristic equation of the equivalent MFR that has the generalized form according to the Gaussian approximation [22] is and it also leads to the following analytical expression: where Here, the normalized equivalent MFR R 0 is the result of the numerical solution of the characteristic Equation (7) after substitution of the researched optical fiber geometrical parameters and the particular azimuthal and radial orders of the analyzed mode.In EMGA (as well as in the "classical" Gaussian approximation), R 0 is a basic single variational parameter, which completely describes the mode dispersion characteristics.Following the substitution of R 0 to the variational expression expressed in Equation (5), this permits the evaluation of the core mode parameter U, which relates with the propagation constant β by the following well-known ratio [17,20]: The solution of the characteristic Equation ( 7) by taking into account the further substitution to Equation ( 5) and then to Equation (8) should satisfy the guided mode cut-off condition [1,17,20]: We also propose considering the optical confinement factor P core as the second criterion for the identification of the "ghost" solutions: The last one is also described by the following analytical expression earlier on, which is derived within the Gaussian approximation approach [20,21]: The proposed EMGA provides low degrees of error, with a reduction of total computational time, especially for the calculations of higher-order mode parameters under low error, even during the analysis of large core optical fibers with complicated forms refractive index profiles, including real commercially available MMF graded-index profiles with refractive index fluctuations and technological defects.The results of EMGA verification by the rigorous mixed finite element method are represented by details in earlier published studies [25].

Mode Delay
The mode delay is inversely related to the propagation constant β [1,17,22], and, respectively, to the core mode parameter derivative ∂U 2 /∂λ [21]: where while the equivalent MFR derivative analytical expressions are also written as finite nested sums with the following form [21]: where S (1)

Mode of Chromatic Dispersion Parameter
The same technique is applied to derive the analytical expressions for the mode of the chromatic dispersion parameter D, which is related to the propagation constant β by the well-known formula [1,17,22]: Therefore, by differentiating Equation ( 13), we shall obtain a second partial derivative of the square of the propagation constant: and the expression for chromatic dispersion parameter can be rewritten in the following form [21]: After differentiating Equations ( 14) and ( 15) with respect to wavelength, the desired second derivatives of the square core mode parameter U 2 and the normalized equivalent mode field radius R 0 will be obtained: where and (1) where

Material Dispersion and Refractive Index Profile Parameters
The well-known Sellmeier equation [1,17] is utilized to take into account the material dispersion: where A i and B i are Sellmeier's coefficients (B i is also denoted as the resonance wavelength), which has been empirically measured for GeO 2 -SiO 2 glass [26,27] under several particular dopant concentrations.
Here, we shall apply the method described in detail in the published work [28], to estimate the Sellmeier coefficients at the graded-index profile points.
A passage from the differentiation operator to its square is utilized analogously to the propagation of constant derivatives.This passage will greatly simplify the expressions for the first and second derivatives of the refractive index written via the Sellmeier equation: Therefore, the first and second derivatives of the profile height parameter ∆ can be expressed as follows: ∂∆ ∂λ while derivatives of the local profile parameter h k lead to the following expressions: Finally, by applying Equations ( 22) and ( 23), the first and second derivatives of the normalized frequency V can be expressed by the following formulas: We propose utilizing a previously developed model of a piecewise regular multimode fiber optic link operating in a few-mode regime to simulate laser-excited optical pulse propagation, and to estimate the potentiality of a 100 µm core MMF, as well as an optimized 100 µm core LDMDF for laser-based multi-gigabit data transmission.This model is described in detail in the published work [21], where its experimental verification is also demonstrated.The proposed solution is based on piecewise regular representation, combined with the general approach of the split-step method [29] to simulate the processes of mode mixing and power diffusion occurring due to mode coupling (Figure 1).Here, single silica weakly guiding circular MMF with an arbitrary axially symmetric refractive index profile with single continuous outer cladding is considered.According to the piecewise regular representation, the fiber is divided into regular spans with length ∆z.Inside of the span, the fiber geometry parameters are considered as a constant, and the modes propagate independently without interaction and mixing.It is supposed that each excited guided mode with a propagation constant, varying from one regular span to another span, satisfies a cut-off condition for whole regular spans composing the fiber.Additionally, it is assumed that at the link transmitter end, each excited guided mode begins transferring a single optical pulse with a particular form that is identical to the input signal (for example, Gaussian).During pulse propagation over a regular span, its amplitude decreases, due to mode attenuation.The signal is mainly distorted due to the difference between the group velocities and the amplitudes of modes, i.e., the DMD effect.Additionally, the transfer by the pulse spreading of each mode due to chromatic dispersion for a particular mode is taken into account.
The boundaries of regular spans can be represented generally as ideal axial alignment splices of two almost similar optical fibers with mismatching parameters.However, it is correct only for "straight" fibers.Therefore, we propose simulating fiber bends by the representation of boundaries in the form of splices of two mismatched fibers with some low angular misalignment [21].The mode power redistribution between the amplitudes of signal components as a result of mode interaction is estimated by mode coupling of coefficient matrix computing at the joints of regular spans.Here, only guided modes are considered as the main issue under pulse dynamics research during the propagation over a multimode fiber link in a few-mode regime.However, power loss due to component transformation from guided to leaky mode and reflections are also indirectly taken into account.
independently without interaction and mixing.It is supposed that each excited guided mode with a propagation constant, varying from one regular span to another span, satisfies a cut-off condition for whole regular spans composing the fiber.Additionally, it is assumed that at the link transmitter end, each excited guided mode begins transferring a single optical pulse with a particular form that is identical to the input signal (for example, Gaussian).During pulse propagation over a regular span, its amplitude decreases, due to mode attenuation.The signal is mainly distorted due to the difference between the group velocities and the amplitudes of modes, i.e., the DMD effect.Additionally, the transfer by the pulse spreading of each mode due to chromatic dispersion for a particular mode is taken into account.
The boundaries of regular spans can be represented generally as ideal axial alignment splices of two almost similar optical fibers with mismatching parameters.However, it is correct only for "straight" fibers.Therefore, we propose simulating fiber bends by the representation of boundaries in the form of splices of two mismatched fibers with some low angular misalignment [21].The mode power redistribution between the amplitudes of signal components as a result of mode interaction is estimated by mode coupling of coefficient matrix computing at the joints of regular spans.Here, only guided modes are considered as the main issue under pulse dynamics research during the propagation over a multimode fiber link in a few-mode regime.However, power loss due to component transformation from guided to leaky mode and reflections are also indirectly taken into account.At the receiver end, the resulting pulse envelope is considered as a superposition of all mode components of the signal.Here, it is proposed that a well-known expression is applied for the frequency response of the signal transferred by the M mode components LPlm over a regular multimode fiber with length z [1]: At the receiver end, the resulting pulse envelope is considered as a superposition of all mode components of the signal.Here, it is proposed that a well-known expression is applied for the frequency response of the signal transferred by the M mode components LP lm over a regular multimode fiber with length z [1]: where F is the direct Fourier transform; h Tx (t) is the initial pulse at the transmitter end; A are first-and second-order dispersion parameters.These dispersion parameters are elements of the well-known Taylor series expansion approximation of the propagation constant frequency dependence β(ω) [1,17,29]: where Here, is the mode delay and 2 is the group velocity dispersion associated with the chromatic dispersion of the p-th guided mode, LP lm .
Therefore, according to introduced piecewise regular representation of the irregular multimode fiber link, the frequency response of a few-mode optical signal, transferred by M guided modes over an irregular MMF with length z under a given particular length of regular span ∆z by taking Equation (30) into account, can be written in the following form [21]: where N z = E(z/∆z); E(x) is the integer part of the real number x.
The resulting pulse response at the receiver end of the irregular multimode link is computed by the following simple expression: where F −1 is the inverse Fourier transform; [x]* is the complex conjugate of x.
Therefore, EMGA is applied for computing the dispersion parameters of the desired selected order guided modes at each regular span of the researched irregular MMF, including mode delays (or group velocity) and chromatic dispersion parameters.Differential mode attenuation is estimated by the known empirical expression proposed by Yabre in [30], which is based on experimental data represented by the same author in [31]: where µ = 2m + l + 1 is principal mode number; α 0 (λ) is the attenuation of lower-order modes (it is supposed to be equal to the attenuation at the correspondence wavelength mentioned in fiber specification); M 0 is total number of modes satisfying the cutoff condition for the analyzed fiber: where g is the gradient factor of the smoothed α-profile.
The model presented passed experimental approbation: a good agreement between its simulation and direct experimental measurements of pulse response was obtained, which was described in detail in the published work [21].

Mode Coupling
Research on the irregularity of MMF is concerned with core diameter variations.Therefore, we propose setting it directly via an array from the reports of optical fiber diameter measurements that were produced during fiber drawing.Micro-and macro-bends are simulated by random equivalent low (θ = 2.0 • . . .4.0 • ) angular misalignments [21] at the boundaries of regular spans, while mode coupling coefficient redistribution at the span boundaries as well as at the optical interconnections of the transmitter and receiver ends may be estimated by the well-known mode field overlap integral method [1,17,22], taking into account the particular angular misalignment inserted.Here, we utilized the overlap integral method in combination with the introduced EMGA, which takes into account the local features of the real silica optical fiber refractive index profile and decreases the computational error.A passage from the generalized form of the overlap, integral to the analytical expression for the arbitrary order mode coupling coefficient estimation at the centralized splice of the optical fibers with mismatched parameters without any misalignments, was proposed in [32].The analytical expression derivation is described in detail in [33].Finally, the formula for the arbitrary order mode coupling coefficient at the central splice is written as follows: where mode coupling occurs only for modes with the same azimuthal order l; Γ is the gamma function; ρ m and ρ n are the injected LP lm and excited LP ln mode field radiuses, respectively.The analytical expression for the arbitrary order mode coupling coefficient at the optical fiber splice under a low angle misalignment θ < 10 • derived and represented in [33] has the following form: where 1 F 1 is the confluent hypergeometric function of the first kind [23,24]: where b (l m ,m−1) q and b (l n ,n−1) q are Laguerre polynomial expansion factors of Equation ( 6); n θ is the refractive index of the launching medium (air gap, core of adjusting/exiting fiber, etc.).

Low Differential Mode Delay Fiber with a Large 100 µm Core Diameter
In the first step, we utilized a LOMF 50/125 refractive index profile scaled to a 100 µm core diameter as the reference for further computation of the selected guided mode staff DMD diagram, to define the boundaries of the reference mode delay t REF range.The reference graded refractive index profile was set by a data protocol of measurements, performed by a certified optical fiber analyzer for a sample of a commercially available LOMF of ISO/IEC Cat.OM2+/OM3 [34], and scaled up to a 100 µm core.This is represented in Figure 2. for a sample of a commercially available LOMF of ISO/IEC Cat.OM2+/OM3 [34], and scaled up to a 100 μm core.This is represented in Figure 2.This profile is quite smooth, and it may be considered as an almost graded power law profile.However, it is not a pure α-profile, and it differs from the last one by the "non-power-law" envelope and by weak local fluctuations of the refractive index improving close to the fiber core center area.By using EMGA, we performed an analysis of a reference 100 μm core MMF with a refractive index over the "O"-band.According to the computational results, it supports the propagation of 145 modes This profile is quite smooth, and it may be considered as an almost graded power law profile.However, it is not a pure α-profile, and it differs from the last one by the "non-power-law" envelope and by weak local fluctuations of the refractive index improving close to the fiber core center area.By using EMGA, we performed an analysis of a reference 100 µm core MMF with a refractive index over the "O"-band.According to the computational results, it supports the propagation of 145 modes LP lm (l = 0, . . ., 24; m = 1, . . ., 11) satisfying the cutoff conditions with an optical confinement factor P co > 0.5, whose distribution over the total mode staff at wavelength λ = 1310 nm is shown in Figure 3.This profile is quite smooth, and it may be considered as an almost graded power law profile.However, it is not a pure α-profile, and it differs from the last one by the "non-power-law" envelope and by weak local fluctuations of the refractive index improving close to the fiber core center area.By using EMGA, we performed an analysis of a reference 100 μm core MMF with a refractive index over the "O"-band.According to the computational results, it supports the propagation of 145 modes LPlm (l = 0, ... ,24; m = 1, … ,11) satisfying the cutoff conditions with an optical confinement factor Pco > 0.5, whose distribution over the total mode staff at wavelength λ = 1310 nm is shown in Figure 3.During the design, with an MMF with an improved bandwidth under a multimode regime, all guided mode staff delays should be equalized, while according to the presented method, only guided modes with normalized amplitudes of more than 0.1 should be considered under particular launching conditions.Here, we propose the centralized launching of the optical signal generated by conventional Fabry-Perot laser diode-generating emission at λ = 1310 to the MMF core end via a standard single mode optical fiber (ITU-T Recommendations G.652), which in the "worst case" may support the propagation of two modes, LP01 and LP11, at the mentioned wavelength.According to the model described above for the simulation of a piece-wise regular multimode fiber optic link operating in a few-mode regime, mode mixing and power diffusion due to micro-and macro-bends are During the design, with an MMF with an improved bandwidth under a multimode regime, all guided mode staff delays should be equalized, while according to the presented method, only guided modes with normalized amplitudes of more than 0.1 should be considered under particular launching conditions.Here, we propose the centralized launching of the optical signal generated by conventional Fabry-Perot laser diode-generating emission at λ = 1310 to the MMF core end via a standard single mode optical fiber (ITU-T Recommendations G.652), which in the "worst case" may support the propagation of two modes, LP 01 and LP 11 , at the mentioned wavelength.According to the model described above for the simulation of a piece-wise regular multimode fiber optic link operating in a few-mode regime, mode mixing and power diffusion due to micro-and macro-bends are simulated as connections of the researched MMF regular spans performed with weak angular misalignment φ.This supposition was confirmed by experimental verification of the model [21], while it was found that for MMF coiled on a typical fiber spool, this angular misalignment can vary as much as φ = 2.0-3.5 • .Therefore, we performed the computation of normalized amplitude distribution dynamics over a guided mode staff under centralized mode LP 01 and LP 11 launching conditions, and further mode power redistribution due to bends for 11 regular spans of piece-wise regular representation.The results of the mode staff normalized amplitude dynamics calculation are shown in Figure 4. Here, we noticed that only 21 guided modes LP lm (l = 0, . . ., 7; m = 1, . . ., 5) transfer the majority of the optical signal during its propagation over the reference 100 µm core MMF, and that their delays vary from 4.922 up to 4.924 ns/km, defining the range for the reference value t REF setting, while normalized amplitudes of other components are less than 0.1.
We performed the synthesis of LDMDF 100/125 desired refractive index profiles for various combinations of t BASE and α-parameters of a first iteration graded profile.The result of 100 µm core LDMDF refractive index profile optimization under t REF = 4923.08ns/km and α = 1.900 for centralized launching conditions is shown in Figure 5, while a comparison between the DMD diagrams for the reference MMF and a sample of LDMDF is shown in Figure 6.
conditions, and further mode power redistribution due to bends for 11 regular spans of piece-wise regular representation.The results of the mode staff normalized amplitude dynamics calculation shown in Figure 4. Here, we noticed that only 21 guided modes LPlm (l = 0, … ,7; m = 1, … ,5) transfer the majority of the optical signal during its propagation over the reference 100 μm core MMF, and that their delays vary from 4.922 up to 4.924 ns/km, defining the range for the reference value tREF setting, while normalized amplitudes of other components are less than 0.1.We performed the synthesis of LDMDF 100/125 desired refractive index profiles for various combinations of tBASE and α-parameters of a first iteration graded profile.The result of 100 μm core LDMDF refractive index profile optimization under tREF = 4923.08ns/km and α = 1.900 for centralized launching conditions is shown in Figure 5, while a comparison between the DMD diagrams for the reference MMF and a sample of LDMDF is shown in Figure 6.Spectral curves of DMD regarding to the fundamental mode DMD LP01 (λ) computed over the "O"-band for both the reference 100 µm core MMF and LDMDF samples are represented in Figure 7: (a) corresponds to the total "O" band, while (b) shows the LDMDF and DMD curve fragment over a λ = 1310 nm wavelength region.For the reference MMF, DMD varies from 1589.7 ps/km up to 1950.5 ps/km over the "O"-band, while the optimized refractive index profile provides a decrease of DMD from 5.0-6.0 times at the boundaries of the "O"-band, and up to 9.0 and more times at a λ = 1310 nm region in comparison with the reference MMF.Here, the total value of DMD LP01 is not more than 300 ps/km over all of the "O"-band; the lowest DMD LP01 is 176.2 ps/km, and it corresponds to 1285 nm; for λ = 1310 nm, DMD LP01 255.7 ps/km.Spectral curves of DMD regarding to the fundamental mode DMDLP01(λ) computed over the "O"-band for both the reference 100 μm core MMF and LDMDF samples are represented in Figure 7: (a) corresponds to the total "O" band, while (b) shows the LDMDF and DMD curve fragment over a λ = 1310 nm wavelength region.For the reference MMF, DMD varies from 1589.7 ps/km up to 1950.5 ps/km over the "O"-band, while the optimized refractive index profile provides a decrease of DMD from 5.0-6.0 times at the boundaries of the "O"-band, and up to 9.0 and more times at a λ = 1310 nm region in comparison with the reference MMF.Here, the total value of DMDLP01 is not more than 300 ps/km over all of the "O"-band; the lowest DMDLP01 is 176.2 ps/km, and it corresponds to 1285 nm; for λ = 1310 nm, DMDLP01 = 255.7 ps/km.

Simulation of 10GBase-LX Fiber Optic System Optical Pulse Propagation
In the first step, we simulated a 10GBase-LX fiber optic system optical pulse propagation over a 100 μm core MMF with a total length of 2 km, and described the above refractive index profile presented in Figure 2, which corresponds to the measured real commercially available sample of the LOMF ISO/IEC 11801 OM2+/OM3 Cat.(MMF 50/125) profile [34], which was scaled up to a 100 μm core region.According to the previous paragraph results, this MMF provides the DMD with about 1800 ps/km at an operating wavelength of the 10GBase-LX system-λ = 1310 nm.
Fiber irregularity is considered as the core diameter variations that are assigned by the array, taken from reports of commercially available MMF 50/125 outer diameter measurements performed during fiber drawing [35] under a mean deviation of ±0.271 μm from the average value of 125.022 μm.Its fragment of over 2 km in length regarding the reference value of the core radius and statistics are shown in Figure 8.The length of a regular span is set to be equal to the fiber diameter measurement workstation spatial resolution, ∆z = 8 m.Micro-and macro-bends of the researched fiber are simulated by random equivalent low angular misalignments [21] at the boundaries of regular spans.We consider a "stressed and bended" MMF that is simulated by normal distribution under a mean of 3.0° and a deviation of 0.5°.Its distribution over 2 km of MMF length and its statistics are presented in Figure 9.

Simulation of 10GBase-LX Fiber Optic System Optical Pulse Propagation
In the first step, we simulated a 10GBase-LX fiber optic system optical pulse propagation over a 100 µm core MMF with a total length of 2 km, and described the above refractive index profile presented in Figure 2, which corresponds to the measured real commercially available sample of the LOMF ISO/IEC 11801 OM2+/OM3 Cat.(MMF 50/125) profile [34], which was scaled up to a 100 µm core region.According to the previous paragraph results, this MMF provides the DMD with about 1800 ps/km at an operating wavelength of the 10GBase-LX system-λ = 1310 nm.
Fiber irregularity is considered as the core diameter variations that are assigned by the array, taken from reports of commercially available MMF 50/125 outer diameter measurements performed during fiber drawing [35] under a mean deviation of ±0.271 µm from the average value of 125.022 µm.Its fragment of over 2 km in length regarding the reference value of the core radius and statistics are shown in Figure 8.The length of a regular span is set to be equal to the fiber diameter measurement workstation spatial resolution, ∆z = 8 m.Micro-and macro-bends of the researched fiber are simulated by random equivalent low angular misalignments [21] at the boundaries of regular spans.We consider a "stressed and bended" MMF that is simulated by normal distribution under a mean of 3.0 • and a deviation of 0.5 • .Its distribution over km of MMF length and its statistics are presented in Figure 9.
taken from reports of commercially available MMF 50/125 outer diameter measurements performed during fiber drawing [35] under a mean deviation of ±0.271 μm from the average value of 125.022 μm.Its fragment of over 2 km in length regarding the reference value of the core radius and statistics are shown in Figure 8.The length of a regular span is set to be equal to the fiber diameter measurement workstation spatial resolution, ∆z = 8 m.Micro-and macro-bends of the researched fiber are simulated by random equivalent low angular misalignments [21] at the boundaries of regular spans.We consider a "stressed and bended" MMF that is simulated by normal distribution under a mean of 3.0° and a deviation of 0.5°.Its distribution over 2 km of MMF length and its statistics are presented in Figure 9.We simulated the propagation of a Gaussian-shaped optical pulse with a full width of half maximum FHWM = 90.91 ps that corresponds to a 10 Gbps bit rate excited by laser diodes operating at a wavelength λ = 1310 nm, generating two LP-mode (linear polarized mode) optical emissions containing the fundamental LP01 and the higher-order mode LP11.Here, the "worst case" of launching conditions was considered, which corresponds to a connection via a conventional fiber optic adapter simulated as a connection with an angular misalignment θ = 4.20° [36] between a standard singlemode optical fiber (ITU-T Recommendations G.652) of a laser source pigtail and a link 100 μm core MMF.
The results of the simulation of 10GBase-LX optical pulse dynamics during propagation over a 100 μm core MMF with a total length of 2 km without special launching conditions is represented in Figure 10a,b (pulse dynamics and diagram of pulse propagation).It was unexpected, but the optical signal kept its envelope without strong DMD distortions, up to a distance of almost 0.45 km from the transmitter end, even without special launching conditions, which is quite sufficient for on-board cable system fiber optic links that are characterized by extremely short ranges.However, further strong DMD occurred with increasing distance.
In the second step, we researched the potentiality of the designed 100 μm core LDMDF sample potentiality for laser-based 10 Gbps data transmission under the same "the worst case" launching conditions, and variations of both the core diameter and the equivalent angular misalignment at the boundaries of the regular spans.The results of the simulation of 10GBase-LX optical pulse dynamics during propagation over a 100 μm core LDMDF with a total length of 2 km without special launching conditions is represented in Figure 11.The LDMDF sample shows a great potentiality for applications in multi-gigabit network fiber optic links: here, DMD distortions did not occur over the entire monitored 2 km length.We simulated the propagation of a Gaussian-shaped optical pulse with a full width of half maximum FHWM = 90.91 ps that corresponds to a 10 Gbps bit rate excited by laser diodes operating at a wavelength λ = 1310 nm, generating two LP-mode (linear polarized mode) optical emissions containing the fundamental LP 01 and the higher-order mode LP 11 .Here, the "worst case" of launching conditions was considered, which corresponds to a connection via a conventional fiber optic adapter simulated as a connection with an angular misalignment θ = 4.20 • [36] between a standard single-mode optical fiber (ITU-T Recommendations G.652) of a laser source pigtail and a link 100 µm core MMF.
The results of the simulation of 10GBase-LX optical pulse dynamics during propagation over a 100 µm core MMF with a total length of 2 km without special launching conditions is represented in Figure 10a,b (pulse dynamics and diagram of pulse propagation).It was unexpected, but the optical signal kept its envelope without strong DMD distortions, up to a distance of almost 0.45 km from the transmitter end, even without special launching conditions, which is quite sufficient for on-board cable system fiber optic links that are characterized by extremely short ranges.However, further strong DMD occurred with increasing distance.
In the second step, we researched the potentiality of the designed 100 µm core LDMDF sample potentiality for laser-based 10 Gbps data transmission under the same "the worst case" launching conditions, and variations of both the core diameter and the equivalent angular misalignment at the boundaries of the regular spans.The results of the simulation of 10GBase-LX optical pulse dynamics during propagation over a 100 µm core LDMDF with a total length of 2 km without special launching conditions is represented in Figure 11.The LDMDF shows a great potentiality for applications in multi-gigabit network fiber optic links: here, DMD distortions did not occur over the entire monitored 2 km length.

Discussion
The results of optical signal propagation simulation over a 100 μm core MMF and LDMDF with a length of 2 km presented in the previous section allow for the localization of a maximal distance where the optical pulse keeps its envelope, and where a strong DMD "glove" effect does not arise, while estimation of the fiber optic system channel quality requires the evaluation of the bit error ratio (BER) or the Q-factor, by taking into account the dispersion distortions.Therefore, we propose computing the optical signal dispersion at the receiver end under a particular distance by a known method for the processing of pulses with arbitrary forms, as described in detail in the monograph [37], and other issues, while, for the first time, the estimation of the 10GBase-LX channel Q-factor, performed by using "the worst case" link model developed by the IEEE 802.3zTask Force, declared and ratified the ability for laser-based data transmission to proceed over MMFs with a 1 Gbps bit rate [38,39].
By analyzing the computation results of pulse dynamics according to the above-mentioned method, we found that pulse dispersion at the receiver end under a maximal 0.45 km length of the

Discussion
The results of optical signal propagation simulation over a 100 μm core MMF and LDMDF with a length of 2 km presented in the previous section allow for the localization of a maximal distance where the optical pulse keeps its envelope, and where a strong DMD "glove" effect does not arise, while estimation of the fiber optic system channel quality requires the evaluation of the bit error ratio (BER) or the Q-factor, by taking into account the dispersion distortions.Therefore, we propose computing the optical signal dispersion at the receiver end under a particular distance by a known method for the processing of pulses with arbitrary forms, as described in detail in the monograph [37], and other issues, while, for the first time, the estimation of the 10GBase-LX channel Q-factor, performed by using "the worst case" link model developed by the IEEE 802.3zTask Force, declared and ratified the ability for laser-based data transmission to proceed over MMFs with a 1 Gbps bit rate [38,39].
By analyzing the computation results of pulse dynamics according to the above-mentioned method, we found that pulse dispersion at the receiver end under a maximal 0.45 km length of the

Discussion
The results of optical signal propagation simulation over a 100 µm core MMF and LDMDF with a length of 2 km presented in the previous section allow for the localization of a maximal distance where the optical pulse keeps its envelope, and where a strong DMD "glove" effect does not arise, while estimation of the fiber optic system channel quality requires the evaluation of the bit error ratio (BER) or the Q-factor, by taking into account the dispersion distortions.Therefore, we propose computing the optical signal dispersion at the receiver end under a particular distance by a known method for the processing of pulses with arbitrary forms, as described in detail in the monograph [37], and other issues, while, for the first time, the estimation of the 10GBase-LX channel Q-factor, performed by using "the worst case" link model developed by the IEEE 802.3zTask Force, declared and ratified the ability for laser-based data transmission to proceed over MMFs with a 1 Gbps bit rate [38,39].
By analyzing the computation results of pulse dynamics according to the above-mentioned method, we found that pulse dispersion at the receiver end under a maximal 0.45 km length of 100 µm core MMF, providing a negligible DMD effect, is D = 83.583ps.The pulse response is shown in Figure 12a.During the next step, we performed an estimation of the Q-factor for the conventional 9.2 dB budget 10GBase-LX transceiver, by taking into account the pulse dispersion, as well as other distortion factors/penalties according to the model [38,39].The computed eye-diagram envelope is represented in Figure 12b: Here, the Q-factor is 4.968, which is unacceptably lower in comparison with the reference value Q REF = 7.04, which corresponds to the "worst case" BER = 10 −12 required by a 10GBase-LX specification without a forward error correction (FEC) technique [40].As a result, a maximal 100 µm core MMF length L = 0.286 km was localized: here, the total pulse dispersion is D = 62.329 ps (Figure 12a), which provides a desired Q-factor of 8.236 (Figure 12b).12a.During the next step, we performed an estimation of the Q-factor for the conventional 9.2 dB budget 10GBase-LX transceiver, by taking into account the pulse dispersion, as well as other distortion factors/penalties according to the model [38,39].The computed eye-diagram envelope is represented in Figure 12b: Here, the Q-factor is 4.968, which is unacceptably lower in comparison with the reference value QREF = 7.04, which corresponds to the "worst case" BER = 10 −12 required by a 10GBase-LX specification without a forward error correction (FEC) technique [40].As a result, a maximal 100 μm core MMF length L = 0.286 km was localized: here, the total pulse dispersion is D = 62.329 ps (Figure 12a), which provides a desired Q-factor of 8.236 (Figure 12b).Regarding the 100 μm core LDMDF, the optical pulse maintains its envelope over the total fiber length of 2 km, and its dispersion at the receiver end is D = 92.113ps.As expected, this is an unacceptably high value that leads to a low Q-factor of 3.656.For this reason, we also localized the maximal distance for a 100 μm core LDMDF link that provides the required 10GBase-LX specification Q-factor value: here, it is L = 1.020 km with a pulse dispersion D = 43.669ps and Q = 8.660 (Figure 13a  Regarding the 100 µm core LDMDF, the optical pulse maintains its envelope over the total fiber length of 2 km, and its dispersion at the receiver end is D = 92.113ps.As expected, this is an unacceptably high value that leads to a low Q-factor of 3.656.For this reason, we also localized the maximal distance for a 100 µm core LDMDF link that provides the required 10GBase-LX specification Q-factor value: here, it is L = 1.020 km with a pulse dispersion D = 43.669ps and Q = 8.660 (Figure 13a,b).
Photonics 2018, 5, x FOR PEER REVIEW 19 of 22 100 μm core MMF, providing a negligible DMD effect, is D = 83.583ps.The pulse response is shown in Figure 12a.During the next step, we performed an estimation of the Q-factor for the conventional 9.2 dB budget 10GBase-LX transceiver, by taking into account the pulse dispersion, as well as other distortion factors/penalties according to the model [38,39].The computed eye-diagram envelope is represented in Figure 12b: Here, the Q-factor is 4.968, which is unacceptably lower in comparison with the reference value QREF = 7.04, which corresponds to the "worst case" BER = 10 −12 required by a 10GBase-LX specification without a forward error correction (FEC) technique [40].As a result, a maximal 100 μm core MMF length L = 0.286 km was localized: here, the total pulse dispersion is D = 62.329 ps (Figure 12a), which provides a desired Q-factor of 8.236 (Figure 12b).Regarding the 100 μm core LDMDF, the optical pulse maintains its envelope over the total fiber length of 2 km, and its dispersion at the receiver end is D = 92.113ps.As expected, this is an unacceptably high value that leads to a low Q-factor of 3.656.For this reason, we also localized the maximal distance for a 100 μm core LDMDF link that provides the required 10GBase-LX specification Q-factor value: here, it is L = 1.020 km with a pulse dispersion D = 43.669ps and Q = 8.660 (Figure 13a

Conclusions
We here present an alternative fast and simple method for a design refractive index profile of silica LDMDFs with an extremely enlarged core diameter of to 100 µm, for laser-based multi-gigabit short-range optical networks.The results of the 100 µm core LDMDF graded index profile sample optimization performed by the proposed solution, and the following computation of the DMD spectral curve over the "O"-band, were demonstrated and compared with a "reference" 100 µm core MMF, where the refractive index profile is similar to the conventional commercially available LOMF, scaled up to a 100 µm core.For the reference MMF, DMD varies from 1589.7 ps/km up to 1950.5 ps/km over the above-mentioned wavelength band.The optimized refractive index profile provides a decreasing of DMD from 5.0-6.0 times at the boundaries of the "O"-band, up to 9.0 and more times at the λ = 1310 nm region, in comparison with the reference MMF.Here, the DMD is not more than 300 ps/km over the entire "O"-band, while the least DMD LP01 = 176.2ps/km, and it corresponds to 1285 nm, and DMD LP01 = 255.7 ps/km at λ = 1310 nm.We utilized a previously developed alternative method for the simulation of piecewise regular multimode fiber optic links operating in a few-mode regime for the computation of laser-excited 10GBase-LX optical pulse dynamics during its propagation over irregular 100 µm core optical fibers-the reference MMF and the optimized LDMDF.While most of the commercially available MMFs with a 100 µm core diameter are targeted to only multimode regimes, and the utilized transceivers based on multimode LEDs provide low bit rates of 10-100 Mbps of data transmission, our simulation results demonstrate the implementation of 10GBase-LX channels over both short-length MMFs with a LOMF-graded refractive index profile, and scaled up to a 100 µm core and extended distance links provided by LDMDFs with an optimized graded refractive index profile.Here, in spite of uncontrolled launching conditions corresponding to connection via a conventional fiber optic adapter simulated as a connection with an angular misalignment θ = 4.20 • between a standard single-mode optical fiber (ITU-T Recommendations G.652) of a laser source pigtail and a link optical fiber, even the 100 µm core MMF keeps a pulse envelope for a distance of up to 0.5 km, with a total dispersion of about 85 ps, while LDMDF blocks up the DMD over the entire 2 km length with a total pulse dispersion of not more than D = 32 ps under a length L = 0.5 km, and D = 92.50ps under a length L = 2.0 km.According to "the worst case link model" developed by the IEEE 802.3zTask Force, we estimate the Q-factor for the 10GBase-LX channel based on 9.2 dB budget IEEE 802.3aeLX-transceivers and a 100 µm core reference MMF and LDMDF without any special launching conditions under a particular fiber length.As a result, the maximal lengths of both MMF and LDMDF were localized to provide the least reference Q = 7.04 for BER = 10 −12 : for a 100 µm core MMF, this distance was L = 0.286 km under a total pulse dispersion D = 62.329 ps and Q = 8.236; for the 100 µm core LDMDF with a designed special refractive index profile, even under the length L = 1.020 km, the Q-factor is a desirable reference value (Q = 8.660) under a total pulse dispersion of D = 43.669ps.

Figure 1 .
Figure 1.Piecewise regular representation of 200 m of irregular optical fiber with varying core diameters.

Figure 1 .
Figure 1.Piecewise regular representation of 200 m of irregular optical fiber with varying core diameters.

p
and α p are the starting amplitude and the mode attenuation of the p-th guided mode LP lm (p = 1...M); β

Figure 2 .
Figure 2. Reference graded refractive index profile reconstructed by a report of real laser-optimized multimode fiber (LOMF) sample profile measurements and scaled up to a 100 μm core.

Figure 2 .
Figure 2. Reference graded refractive index profile reconstructed by a report of real laser-optimized multimode fiber (LOMF) sample profile measurements and scaled up to a 100 µm core.

Figure 2 .
Figure 2. graded refractive index profile reconstructed by a report of real laser-optimized multimode fiber (LOMF) sample profile measurements and scaled up to a 100 μm core.

Figure 3 .
Figure 3. Optical confinement factor distribution over the mode staff of a reference 100 μm core multimode optical fiber (MMF).

Figure 3 .
Figure 3. Optical confinement factor distribution over the mode staff of a reference 100 µm core multimode optical fiber (MMF).

Figure 4 .
Figure 4. Normalized amplitude distribution dynamics over the guided mode staff of the reference MMF under a centralized mode LP01 and LP11 launching, and following mode mixing and power redistribution due to the fiber bends being represented as a connection of regular spans with weak angular misalignment, which is empirically defined as a range of ϕ = 2.0° … 3.5°.

Figure 4 .
Figure 4. Normalized amplitude distribution dynamics over the guided mode staff of the reference MMF under a centralized mode LP 01 and LP 11 launching, and following mode mixing and power redistribution due to the fiber bends being represented as a connection of regular spans with weak angular misalignment, which is empirically defined as a range of φ = 2.0 • . . .3.5 • .Photonics 2018, 5, x FOR PEER REVIEW 15 of 22

Figure 5 .
Figure 5. Optimized graded refractive index profile for a 100 μm core low differential mode delay fiber (LDMDF).

Figure 5 .
Figure 5. Optimized graded refractive index profile for a 100 µm core low differential mode delay fiber (LDMDF).

Figure 5 .
Figure 5. Optimized graded refractive index profile for a 100 μm core low differential mode delay fiber (LDMDF).

Figure 6 .
Figure 6.DMD diagrams computed for selection under a centralized launching condition mode staff for the reference MMF and LDMDF.

Figure 6 .Figure 7 .
Figure 6.DMD diagrams computed for selection under a centralized launching condition mode staff for the reference MMF and LDMDF.Photonics 2018, 5, x FOR PEER REVIEW 16 of 22

Figure 8 .
Figure 8. Variation of the core diameter along the researched optical fiber length.Figure 8. Variation of the core diameter along the researched optical fiber length.

Figure 8 . 22 Figure 9 .
Figure 8. Variation of the core diameter along the researched optical fiber length.Figure 8. Variation of the core diameter along the researched optical fiber length.Photonics 2018, 5, x FOR PEER REVIEW 17 of 22

Figure 9 .
Figure 9. Variation of the equivalent angular misalignment at the splice of regular spans along the fiber length simulating fiber micro-and macro-bends.

Figure 10 .Figure 10 .Figure 11 .
Figure 10.10GBase-LX optical pulse dynamics during propagation over a 100 µm core MMF with a total length of 2 km without special launching conditions: (a) pulse dynamics; (b) diagram of pulse propagation.

Figure 11 .
Figure 11.10GBase-LX optical pulse dynamics during propagation over a 100 µm core LDMDF with a total length of 2 km without special launching conditions: (a) pulse dynamics; (b) diagram of pulse propagation.

Figure 13 .
Figure 13.10GBase-LX data transmission over a 100 μm core LDMDF without special launching conditions: (a) pulse response; (b) eye diagram envelope.

Figure 13 .
Figure 13.10GBase-LX data transmission over a 100 μm core LDMDF without special launching conditions: (a) pulse response; (b) eye diagram envelope.

Figure 13 .
Figure 13.10GBase-LX data transmission over a 100 µm core LDMDF without special launching conditions: (a) pulse response; (b) eye diagram envelope.
Photonics 2018, 5, x FOR PEER REVIEW 19 of 22 100 μm core MMF, providing a negligible DMD effect, is D = 83.583ps.The pulse response is shown in Figure