Abstract

A faster and accurate semianalytical formulation with a robust optimization solution for estimating the splice loss of graded-index fibers has been proposed. The semianalytical optimization of modal parameters has been carried out by Nelder-Mead method of nonlinear unconstrained minimization suitable for functions which are uncertain, noisy, or even discontinuous. Instead of normally used Gaussian function, as the trial field for the fundamental mode of graded-index optical fiber a novel sinc function with exponentially and ( is the normalized radius of the optical fiber) decaying trailing edge has been used. Due to inclusion of three parameters in the optimization of fundamental modal solution and application of an efficient optimization technique with simple analytical expressions for various modal parameters, the results are found to be accurate and computationally easier to find than the standard numerical method solution.

1. Introduction

Single mode fiber is considered as the most important broadband transmission media for optical communication system. Achieving accurate values of modal field distribution in such fiber is very essential, as it can provide basic solutions for wave equation and many useful properties like splice loss, microbending loss, fiber coupling, and the prediction of intramodal dispersion [1]. However, the various expressions for the fundamental modal field that have been reported so far are not able to predict propagation constant and modal parameters exactly in all regions of single mode operation [2]. The Gaussian approximation shows poor accuracy for lower normalized frequency region although this region may involve single mode fiber operation [2]; however, it can perform satisfactorily only for higher normalized frequency region and give good result near the cut-off frequency of next higher mode [3]. Besides, it is also equally important that the approximation should describe the field in the cladding accurately, as it is useful in the study of evanescent coupling problem. To overcome these inefficiencies, an exponentially and decaying trailing edge fundamental modal field solution in core-cladding interface region has been considered.

To achieve higher accuracy compared to Gaussian function, the Gaussian-Hankel [2], the generalized Gaussian [4], the extended Gaussian [5], and the Laguerre-Gauss/Bessel expansion approximation [6, 7] have been proposed so far. An approximate analytical description with no requirement for optimization has also been presented [8]. But such analytical expression may not work for all specifications of an optical fiber. In the proposed formulation, Nelder-Mead method of nonlinear unconstrained minimization and the process of minimization of core parameter () for all specific requirements have been used to achieve an accurate and computationally appropriate result.

Unlike the existing reported fundamental modal solution with one or two parameters [28], an attempt has been made to propose a three-parameter fundamental modal field solution for graded-index fiber to introduce more flexibility to solve the fundamental modal solution more accurately, especially in core-cladding interface region wherein the solution has different form (exponentially and decaying trailing edge). Ghatak et al. [9] had arrived at simple analytical expressions to describe different optical fiber characteristics by implementing variational technique. Again, the optimization process requires expressions for propagation constant  and core parameter . The analytical expressions for and used for the present study involve many fiber parameters, such as core radius (), refractive indices of core and cladding ( and ), aspect ratio (), and wavelength (), where is the free space wave number. Hence, any desired specification can be incorporated by varying these parameters. Now, the task of optimization can be carried out by using Nelder-Mead method of nonlinear unconstrained minimization, to meet a particular design.

For graded-index optical fiber at the splices, the power transmission coefficients with transverse and angular mismatch have been estimated by using the methods given by Meunier and Hosain [10] and Hosain et al. [11]. For arbitrarily graded-index fiber, the Gaussian approximation does not give accurate result at lower normalized frequency or in cases where the power law profile deviates from its simplest form [12]. Further, the numerical solution requires rigorous computations and specialized numerical techniques [13]. However, using the proposed three-parameter fundamental modal solution coupled with Nelder-Mead method of nonlinear unconstrained minimization, the algorithm becomes comparatively easier to be implemented on an ordinary personal computer, which provides computationally more efficient result [14, 15] than standard numerical method and yields excellent agreement with exact solutions. This is achieved due to the fact that requisite analytical formulae are deduced beforehand and then parameters of those analytical expressions are found by optimization using Nelder-Mead simplex method for nonlinear unconstrained minimization. Furthermore, Nelder-Mead simplex method for nonlinear unconstrained minimization is a direct search method [16, 17] which does not require any derivative information, so it can optimize nonstationary functions, as needed for the problems under study [1820]. The proposed semianalytical model can also be used in the study of nonlinear fiber [21].

2. Formulation of the Problem

2.1. Theory

Splice loss can be evaluated analytically with the help of the following equations [22]: where and are given in (26). , , and are the three variational parameters present in the fundamental modal solution.

is exponential integral given by [23] as follows: is the cosine integral function, defined by [23], as follows: where is Euler's constant 0.5772.

are the Whittaker functions which are solutions to the Whittaker differential equation [23].

is the generalized hypergeometric function , where [23] with, where is the gamma function [23].

2.2. Basic Formulations

The refractive index profile for a weakly guiding fiber is given by where the normalized profile functions for the trapezoidal and triangular index profiles () are given by Here, is the aspect ratio, is the normalized radius (=), is the core radius, is the actual radius of the optical fiber, and and are, respectively, the refractive indices of the core axis and cladding.

For the present study, the following approximations for the fundamental mode as the trial field have been proposed: where , , and are the three variational parameters present in the fundamental modal solution.

To employ variational technique, first the scalar variational expression for the propagation constant  as given by (27) has been considered and is shown in equations through (28) to (30) as follows: for and , for and , for , where Now, the core parameter is given by

Now for a fixed value of normalized frequency, the core parameter is minimized with respect to the variational parameters , , and . Once the optimized values of these three parameters are obtained, the propagation constant and other design parameters can be obtained as explained in the next section.

2.3. Splice Loss

For small angular misalignment () at the splice of two optical fibers, following Hosain et al. [11], the well-known overlap integral can be represented as where , being refractive index of the index matching fluid joining the fibers and being the angular misalignment.

The transmission coefficient at the splice with angular mismatch can then be expressed as Expanding the exponential term, (33) can be written as and from (33), is given by

According to Hosain et al. [11], only the first four terms in (35a) are enough to obtain sufficient accuracy for misalignment up to , which corresponds to for an optical fiber with and μm working at a wavelength μm. Here, up to the fifth term of (35a) have been calculated and the required expressions are given in (9)-(10).

The transmission coefficient at the splice for a transverse offset is expressed as follows: where is the normalized transverse offset, and in practice for Δ ≤ 0.8, one can approximately write [11] where Integrals given in (38)–(41) can be evaluated by using the expression of fundamental modal field given by (26). Hence, (34) and (36) can be evaluated with the help of (35a), (35b), and (37).

2.4. Evaluation of Integrals

Evaluation of integrals to determine propagation constant and splice loss is presented in (1)–(23). Substituting (1)–(6) into (28), (29), and (30), analytical expression of propagation constant can be obtained with the help of (26) and (27). The transmission coefficient (34) at the splice with angular mismatch can be obtained by substituting (9)-(10) into (35a) and (35b). Using (31), (1), (2), (5), and (6), (38) can be obtained. Equations (40) and (41) can be evaluated using (11). Once (38)–(41) are evaluated, (see (37)) can be calculated. Then, the analytical expression of transmission coefficient at the splice between two identical optical fibers having a transverse offset can be evaluated using (36).

3. Results and Discussions

Detailed comparison between the proposed formulation and available exact numerical results [10, 13] has been carried out in terms of accuracy assessment. It has been justified by many authors [13] that two-parameter approximations are more accurate than single-parameter approximation. The proposed approximation of fundamental field involving three optimizing parameters incorporates more flexibility to modify the fundamental modal solution of optical fibers having different specifications. Optimized values of these parameters for different normalized frequencies are given in Tables 1 and 2 for a particular specification of optical fiber having trapezoidal and triangular index profiles, respectively. Values for other normalized frequencies having different specification of optical fiber can also be obtained by using Nelder-Mead method of nonlinear unconstrained minimization.

In order to verify the feasibility of the proposed approximation, the outcomes of the proposed study have been compared with the earlier reported numerical results [10, 13]. In the present study, and are considered for trapezoidal index profile, which corresponds to a typical dispersion shifted silica fiber with μm, , and zero dispersion wavelength at μm [12]. For triangular index profile, and have been chosen, taking μm, , and zero dispersion wavelength at μm [12].

For evaluation of splice loss, the applicability of the proposed formulations in case of power transmission coefficients and at splices between two identical optical fibers has been considered. Gaussian approximation gives accurate result for the evaluation of transmission coefficient only in the region near the cutoff of single mode operation, but it leads to considerable error throughout the single mode region [11]. The variation of with and with , in case of splicing of two identical triangular index fibers, has been plotted in Figures 1 and 2, respectively. Similarly, the variations of these power transmission coefficients for the case of splicing of two identical trapezoidal index fibers are illustrated in Figures 3 and 4. For the practical range of and , the results obtained by the proposed approximation are identically matching with the exact available and numerical results [10, 13].

4. Conclusions

An accurate three-parameter approximation of fundamental modal field solution of an optical fiber has been presented, which can effectively be used to estimate the power transmission coefficients in case of splicing of two identical single mode graded-index fibers in presence of both transverse and angular misalignments. Taking trapezoidal and triangular index fibers as examples, it has been shown that the results obtained with our function are excellently matching with the exact available and numerical results [10, 13]. Besides providing values of optimizing parameters involved in the approximate field obtained by Nelder-Mead method of nonlinear unconstrained minimization, all related simplified analytical expressions have also been presented, which can be used directly by optical fiber designer while predicting splice losses of an optical fiber, having triangular and trapezoidal index profiles for a wide range of normalized frequencies. The salient features of the proposed solution are easy computation on an ordinary personal computer and a robust algorithm for nonlinear unconstrained optimization being applied in an optical fiber having triangular and trapezoidal index profiles.