Artificial Immune Network Design of Optical Multiplexers / Demultiplexers

An efficient design approach of directional couplers based multiplexers/demultiplexers for optical communication applications, by using an adapted artificial immune network algorithm for optimization (opt-AiNet), is presented and validated by using the beam propagation method. Two key multiplexers/demultiplexers based on planar waveguides and optical fiber, directional couplers, have been optimized in order to validate the efficiency and usefulness of the opt-AiNet.

The artificial immune network algorithm for optimization and its adaptations is presented in section II, couplers and multiplexers based in, are presented in section III, divided into planar and optical fiber based items.Optimization results are given in section IV and finally the conclusions and work in progress are presented.

II. IMMUNE NETWORK ALGORITHM FOR OPTIMIZATION -OPT-AINET
The Artificial Immune System (AIS) researchers and techniques have recently been expanded by the benchmarking and tackling real-world applications [4]- [10].According to [5], its story starts on 90's, when Forrest et al. [6] proposed a selection algorithm based on the generation of T-cells to discriminate the antigen-antibody, called negative selection method.Therefore, the AIS are methaheuristics techniques being widely explored in different applications, in economy analysis [7], [8] and environmental techniques and processes improvements [9], [10].
The artificial immune network algorithm is one of the computational representation from the vertebrate immune system, being a general schematic between the antigenic and antibody relations [11], [12], [13].Herein, is considered the number of antibodies fabricated in B-cells (Bone Marrow) to identify their geometrically complementarity, also called affinity [13].
The affinity level implies in the B-cell expansions (Figure 1), consequently, in the fabrication of new antibodies by cloning, which its can suffer some random changes (mutations) in their characteristics.It is an important concept adapted in the algorithm as important immune operator to improve the optimization performance.These new antibodies and antigenic population will be aggregated in the previous population to re-measure the affinities fees again.
Considering this full population, some macrophage or T-cell are signalized to suppress some antibodies/antigen, keeping in account a pre-determined region to evaluate the individuals affinities and maintain only best antibodies [11].This complete procedure of antigenic interconnection, clonal expansion, mutation and cell selections is called clonal selection [12].
From this macro vertebrate immune system description, this work presents an adapted artificial immune system, as shown in Fig. 2, under the opt-AiNet version, combining two different mutation operators and an updated antibody suppression process to improve the convergence performance for computational electromagnetic applications.The herein Opt-AiNet algorithm presented has been developed under the Object Orientation definitions by the C++ language, allowing its reuse in further different applications, and follows the main theoretical aspects described in [12], [13], [17], [18].Furthermore, this computational implementation is fast enough to be disregarded in computational powerful requirements analysis, being intrinsically related only to the objective functions.For examples, when the numerical methods, such as FDTD or FEM, are adopted to model any electromagnetic devices as presented in 0 [16], [19], [20], [21].

III. COUPLER BASED MULTIPLEXER
The coupler based multiplexer/demultiplexers, as shown in Fig. 3, separate or combine two wavelength by using the directional coupler, the coupling lengths, Lc1 and Lc2, corresponding to λ1 and λ2, respectively, must satisfy the condition L = n Lc1 = m Lc2, where n and m must be odd and even integer numbers (or vice-versa).This ensures that each power signal will be present only in one of the outputs.The operation principle can be seen in Fig. 4, where n = 5 (odd) and m = 4 (even).where, P1 and P2 are the waveguide output power.Pin is the input power, C is the coupling factor obtained by the coupling mode theory and z is the distance measured since the input port.Therefore, the total power transfer to the other waveguide occurs when C.z = π/2.Consequently, the coupling distance can be expressed by (3).
For planar waveguide DCs, the coupling factor C is given by (4), where d is the waveguide width, S the distance between the waveguides, ) 0.5 and k 0 = 2π/λ is the wave number in free space, λ is the wavelength, n 1 and n 2 is the guide and substrate refractive indexes, respectively.The n eff is the effective index obtained by solving the transcendental equation (5).
The coupling length in optical fiber DCs shown in Fig. 2. (b) can be obtained in an approximate way by using the analytical solution given in [1], where, W= ak0(neff2-n12)1/2, U = ak0 (n22-neff2)1/2, a is the optical fiber radius, d is the distances between the fibers from their centers, k0 is the free space wavenumber given by 2π/λ, λ is the free space wavelength, n1 and n2 are the refractive indexes of the core and the cladding, respectively.K0 and K1, are the Bessel functions.
The n eff is the effective index obtained by solving the transcendental equation (5).
where J 0 and J 1 are Bessel functions.We can observe the nonlinear relation among all the parameters.

IV. NUMERICAL RESULTS
The first multiplexer/demultiplexer device designed is the planar waveguide operating in wavelengths of 1.41 µm and 1.55 µm.The device length has been imposed as a restriction, in this way L, should be 200 µm.The AIS algorithm presented satisfactory results for the convergence and for the device behavior.The parameters obtained after the optimization attended the imposed restrictions.In In Fig. 5. it is possible to observe the propagation of the input signal which is placed at the left port, and after the propagation it can reach the left or the right port, depending on the operating wavelength.The results obtained for the planar mux/demux are summarized in Table I.The second application is the design of multiplexer/demultiplexer based on optical fiber waveguides.In this case, two CDWM wavelengths were chosen, being λ1 = 1.29 μm and λ2 = 1.59 μm.The maximum device length was fixed to be L = 1 mm.The core and cladding of the optical fiber were considered to be 13.5% GeO2−86.5% SiO2 and quenched SiO2, respectively.
The wavelength dependence of the refractive indexes has been taken into account by using the Sellmeier coefficients given in [22].The values of a and d were computed by using the proposed AIS.
The objective functions, for both applications, have penalties and bonuses in accordance to the limits of Lc1, Lc2, P1 and P2.For Lc1 and Lc2 the objective functions must be higher than 0 and lower than Lc limits, when these conditions are not satisfied, the final objective value is given by 12 .fit Lc Lc  From the initial conditions, power outputs P1 and P2 promote the final objective value.When P1 and P2 have values higher than 0.98 or lower than 0.002, the final value is given by   In order to validate these results, the beam propagation method in frequency domain [23] has been used to simulate the propagation of the fundamental mode along this coupler.The results are shown in Fig. 7 and 8.It can be observed that the field corresponding to λ1 = 1.29 μm appears at the same input fiber after propagating 941 μm (Fig. 7), while for λ2 = 1.59 μm the field appears at the second fiber after propagating 941 μm (Fig. 8).Another observation is that overlapping field in both cases can be The results obtained for the optical fiber mux/demux are summarized in Table II.

V. CONCLUSIONS
The use of efficient parallel platforms and algorithms may certainly help to reduce substantially the CPU time.In this work, the AIS in conjunction with analytical formulae have been used to design planar and optical fiber couplers based CDWM multiplexers/demultiplexers in an efficient way.The optimizations presented here showed a satisfactory convergence, optimized results with the constraints imposed and short runtime.The results were validated by using the FEM-BPM.

Fig. 7 Fig. 8
Fig. 7 Spatial distribution of the field intensity at the input and output of the coupler for λ 1 = 1.29 μm (a) z = 0 and (b) z = 941 μm.