Advances in Optical TechnologiesVolume 2013 (2013), Article ID 168923, 8 pageshttp://dx.doi.org/10.1155/2013/168923

Research Article

Splice Loss of Graded-Index Fibers: Accurate Semianalytical Descriptions Using Nelder-Mead Nonlinear Unconstrained Optimization with Three-Parameter Fundamental Modal Field

Raja Roy Choudhury,¹ Arundhati Roy Choudhury,² and Mrinal Kanti Ghose³

¹Applied Electronics and Instrumentation Department, Sikkim Manipal Institute of Technology, Majitar, Sikkim 737136, India²Physics Department, Sikkim Manipal Institute of Technology, Majitar, Sikkim 737136, India³Computer Science Department, Sikkim Manipal Institute of Technology, Majitar, Sikkim 737136, India

Received 5 March 2013; Accepted 4 June 2013

Academic Editor: Yong Zhao

Copyright © 2013 Raja Roy Choudhury et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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 [2–8], 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 [18–20]. 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 [1–3] 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.

Table 1: Values of optimizing parameters with different normalized frequencies for trapezoidal index profile.

1.839506
1.487480
0.006634
1.6000

1.854451
1.152568
0.040590
1.8000

1.878781
1.006541
0.097612
2.0000

1.908874
0.923490
0.171137
2.2000

1.941083
0.867514
0.253638
2.4000

1.954907
0.814689
0.289126
2.6000

1.987346
0.781444
0.277328
2.8000

1.967303
0.722072
0.256257
3.0000

2.207612
0.797867
0.283156
3.2000

1.930533
0.637208
0.226140
3.4000

Table 2: Values of optimizing parameters with different normalized frequencies for triangular index profile.

1.836461
1.962401
−0.000309
1.7000

1.843169
1.323589
0.014877
1.9000

1.858852
1.102224
0.050720
2.1000

1.881018
0.986778
0.103042
2.3000

1.907224
0.914502
0.167036
2.5000

1.934891
0.862840
0.237307
2.7000

1.954757
0.817735
0.290207
2.9000

1.947828
0.764272
0.271233
3.1000

1.941593
0.718887
0.255127
3.3000

2.150373
0.784010
0.278234
3.5000

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].

Figure 1: Variation of power transmission coefficients with the normalized angular offset p for splicing of two identical single mode triangular index fibers with (exact numerical results [10, 13]; results by our approximation; results based on Gaussian approximation [10, 13]).

Figure 2: Variation of power transmission coefficients with the normalized transverse offset Δ for splicing of two identical single mode triangular index fibers with .7 (exact numerical results [10, 13]; results by our approximation; results based on Gaussian approximation [10, 13]).

Figure 3: Variation of power transmission coefficients with the normalized angular offset for splicing of two identical single mode trapezoidal index fibers with (exact numerical results [10, 13]; results by our approximation; results based on Gaussian approximation [10, 13]).

Figure 4: Variation of power transmission coefficients with the normalized transverse offset Δ for splicing of two identical single mode trapezoidal index fibers with (exact numerical results [10, 13]; results by our approximation; results based on Gaussian approximation [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.

References

R. Tewari, S. I. Hosain, and K. Thyagarajan, “Scalar variational analysis of single mode fibers with Gaussian and smoothed-out profiles,” Optics Communications, vol. 48, no. 3, pp. 176–180, 1983.

A. Sharma, S. I. Hosain, and A. K. Ghatak, “The fundamental mode of graded-index fibres: simple and accurate variational methods,” Optical and Quantum Electronics, vol. 14, no. 1, pp. 7–15, 1982.

A. Sharma and A. K. Ghatak, “A variational analysis of single mode graded-index fibers,” Optics Communications, vol. 36, no. 1, pp. 22–24, 1981.

A. Ankiewicz and G.-D. Peng, “Generalized Gaussian approximation for single-mode fibers,” Journal of Lightwave Technology, vol. 10, no. 1, pp. 22–27, 1992.

S. Chieh, W.-H. Tsai, and M.-S. Wu, “Extended Gaussian approximation for single-mode graded-index fibers,” Journal of Lightwave Technology, vol. 12, no. 3, pp. 392–395, 1994.

G. De Angelis, G. Panariello, and A. Scaglione, “Variational method to approximate the field of weakly guiding optical fibers by Laguerre-Gauss/Bessel expansion,” Journal of Lightwave Technology, vol. 17, no. 12, pp. 2665–2674, 1999.

F. Chiadini, G. Panariello, and A. Scaglione, “Variational analysis of matched-clad optical fibers,” Journal of Lightwave Technology, vol. 21, no. 1, pp. 96–105, 2003.

Q. Cao and S. Chi, “Approximate analytical description for fundamental-mode fields of graded-index fibers: beyond the Gaussian approximation,” Journal of Lightwave Technology, vol. 19, no. 1, pp. 54–59, 2001.

A. K. Ghatak, R. Srivastava, I. F. Faria, K. Thyagarajan, and R. Tiwari, “Accurate method for characterising single-mode fibres: theory and experiment,” Electronics Letters, vol. 19, no. 3, pp. 97–99, 1983.

J. P. Meunier and S. I. Hosain, “An efficient model for splice loss evaluation in single-mode graded-index fibers,” Journal of Lightwave Technology, vol. 9, no. 11, pp. 1457–1463, 1991.

S. I. Hosain, A. Sharma, and A. K. Ghatak, “Splice-loss evaluation for single-mode graded-index fibers,” Applied Optics, vol. 21, no. 15, pp. 2716–2720, 1982.

U. C. Paek, “Dispersionless single-mode fibers with trapezoidal-index profiles in the wavelength region near 1.5 m,” Applied Optics, vol. 22, no. 15, pp. 2363–2369, 1983.

K. Kawano and T. Kitoh, Introduction to Optical Waveguide Analysis, John Willy and Sons, New York, NY, USA, 2001.

J. C. Lagarias, J. A. Reeds, M. H. Wright, and P. E. Wright, “Convergence properties of the Nelder-Mead simplex method in low dimensions,” SIAM Journal on Optimization, vol. 9, no. 1, pp. 112–147, 1998.

J. M. Parkinson and D. Hutchinson, “An investigation into the efficiency of variants on the simplex method,” in Numerical Methods for Nonlinear Optimization, F. A. Lootsma, Ed., pp. 115–135, Academic Press, New York, NY, USA, 1972.

M. H. Wright, “Direct search methods: once scorned, now respectable,” in Proceedings of the Dundee Biennial Conference in Numerical Analysis, D. F. Griffiths and G. A. Watson, Eds., pp. 191–208, Addison Wesley; Longman, Harlow, UK, 1996.

J. A. Nelder and R. Mead, “A simplex method for function minimization,” Computer Journal, vol. 7, pp. 308–313, 1965.

T. H. Rowan, Functional stability analysis of numerical algorithms [Ph.D. thesis], University of Texas, Austin, Tex, USA, 1990.

S. Singer and S. Singer, “Complexity analysis of Nelder-Mead search iterations,” in Proceedings of the 1st Conference on Applied Mathematics and Computation, Dubrovnik, Croatia, 1999, M. Rogina, V. Hari, N. Limić, and Z. Tutek, Eds., pp. 185–196, PMF-Matematički odjel, Zagreb, 2001.

S. Singer and S. Singer, “Efficient implementation of the Nelder-Mead search algorithm,” Applied Numerical Analysis and Computational Mathematics, vol. 1, no. 3, pp. 524–534, 2004.

R. Roychoudhury and A. Roychoudhury, “Accurate semi analytical model of an optical fiber having Kerr nonlinearity using a robust nonlinear unconstrained optimization method,” Optics Communications, vol. 284, no. 4, pp. 1038–1044, 2011.

I. S. Gradshteyn and I. M. Ryzhik, Table of Integrals, Series and Products, Academic Press, New York, NY, USA, 1980.

M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover, New York, NY, USA, 1981.