PIERS ONLINE, VOL. 3, NO. 4, 27 462 The Effect of Radiation Coupling in Higher Order Fiber Bragg Gratings Li Yang 1, Wei-Ping Huang 2, and Xi-Jia Gu 3 1 Department EEIS, University of Science and Technology of China, China 2 Department ECE, McMaster University, Canada 3 Department ECE, Ryerson University, Canada Abstract The effect of coupling in second and higher order (marked as N-th order) fiber Bragg gratings is investigated by using coupled mode theory, where the field is solved and expressed in terms of guided waves by using Green s function. Different from the direct index coupling between two contra-propagating guided modes due to the N-th order diffraction, the coupling lead to the indirect self or mutual complex coupling between the two guided waves via the diffractions lower than N-th order. It will become dominant and result in some interesting properties potentially useful in the DFB fiber laser at some special cases such as a rectangular index modification with the duty cycle equals.5. For this reason, the dependences of coupling coefficients, spectrum characteristics and the bandwidth of Bragg peak on the parameters of gratings with rectangular profiles, such as duty cycle, index change and grating length, are simulated and discussed in this paper. DOI: 1.2529/PIERS61613536 1. INTRODUCTION Over the past decades, much theoretical and experimental work has been devoted to the first order fiber Bragg Grating (FBG), while there are fewer reports on the higher order FBGs possibly because of the limitation of the fabrication techniques. However, the distinct properties of higher order FBGs indicate potential applications in laser, sensor and etc, which draw more and more interests 1 5]. One notable property in higher order FBGs is the complex coupling via the diffraction lower than N-th order which is the order of FBG. It means, in the total coupling coefficient, in addition to index coupling coefficient due to the direct coupling between two contra-propagating guided waves similar to the first order gratings, there is coupling coefficient indicating the indirect coupling between the two guided waves via the excitation and resonant of the wave. Since the index coupling coefficient is proportional to the index change if it does not equal zero, while coupling coefficient is proportional to the square of the index change, the effect of coupling will be prominent only when the former one is around zero at some special cases such as a rectangular index modification with the duty cycle equaling.5 3]. At such cases, the complex coupling coefficient provides possibility to control the total coupling coefficient and then the spectrum characteristics in a wider range, which has been proved to be useful in DFB laser based on multilayer planar waveguides 6, 7]. On the other hand, in experiments, usually, the perturbed index profile in FBG is almost sinusoidal; and it will be more similar to rising-cosine after the saturation of the index modification 5]. During the growth of the index change, the higher order gratings attributed to the effect of the higher order Fourier expansion coefficients show up; meanwhile, the coupling enhances gradually, which offers the possibility and necessity to study the effect of coupling in higher order FBGs. In this paper, the spectrum characteristics of higher order FBGs with rectangular index modification in a theoretical model are investigated by using coupled mode theory. Hereinto, the fields are determined in terms of guided waves by using Green s function method. The influence of grating parameters such as duty cycle, index change and grating length on the coupling coefficients and the bandwidth of Bragg peak are simulated and discussed.
PIERS ONLINE, VOL. 3, NO. 4, 27 463 2. THEORETICAL MODEL AND BASIC FORMULATIONS A schematic diagram of a cylindrical symmetric single mode FBG in the theoretical model is shown in Fig. 1. In cylindrical coordinates based on the exp(jωt) time variation, the refractive index in the perturbed optical fiber can be expressed as, ñ 2 ñ 2 co = n 2 co + η m exp ( j 2mπ Λ z) = n 2 + η m exp ( j 2mπ Λ z), < r r (r, z)= m= m= (1) m n 2 cl, r > r where r, n co and n cl are respective the core radius, the refractive index of the core and the cladding in the unperturbed fiber, while n is the refractive index of the core in the reference fiber. The Fourier expansion coefficients of the periodic index modification is η m = 2n co Λ Λ/2 Λ/2 ( f(z) exp j 2mπ ) Λ z dz (2) f(z) is the profile of the periodic modification. Following the previous analysis for a N-th order grating, the electrical field is expressed in terms of two contra-propagation fundamental guided waves (named as LP 1 mode in fiber) and waves due to the resonant via the diffraction lower than N-th 3, 6, 7]. For a fiber with a small index difference, the waves En rad can be solved and expressed in terms of the two guided waves by Green s function in a uniform region 3]. Then from the partial differential equation the electrical field of LP 1 mode satisfies, the following coupled mode equations are obtained as follows, ) ( ) a(z) j da(z) dz db(z) dz = j ( δ + = j ( δ + κ N + m(ab) ) ( b(z) + j κ N + m(ba) sa b(z) ) a(z) (3) In (3), δ is the detune factor, δ = β β = β Nπ Λ κ ±N is the index coupling coefficient analytically expressed as (4) κ ±N = πy A2 r2 k2 η ±N 4β J 2 (k 1tr ) J 2 (k 1t r ) + J 2 1(k 1t r ) ] = k 2 η ±NΓ (5) where Y and A are respectively the mode admittance and normalized amplitude of LP 1 mode in the reference fiber. J (k 1t r ) and J 1 (k 1t r ) are Bessel Functions, where k 1t is the transverse propagation constant of LP 1 mode in the core. More briefly, κ ±N is expressed as the product of the half of the propagation constant in free space, k, the ±N-th order Fourier expansion coefficient, η ±N, and the confinement factor of LP 1 mode, Γ. When N > 1, the analytical expression of self coupling coefficients sa and, mutual coupling coefficients m(ab) and κrad m(ba) are obtained by using Green s function method and analytically expressed as follows, sa = n, n, = n N,N n, m(ab) = n,n n, m(ba) = n N, n (6) i,j jπ 2 Y k 4 = A2 η iη j 4β J 2 (k 1tr ) ( ) k 2 k1t 2 r C(r)J (k 1t r)rdr, n = 1, 2 (7) ] C(r) = krj (k 1t r) H (2) (kr)j 1(kr) H (2) 1 (kr)j (kr) ] +r kh (2) 1 (kr )J (k 1t r ) k 1t H (2) (kr )J 1 (k 1t r ) J (kr) (8)
PIERS ONLINE, VOL. 3, NO. 4, 27 464 where k is propagation constant in uniform region. For symmetrical structure, there are κ N = κ N and κ = rad κrad sa = = κrad m(ab) = κrad m(ba). The total coupling coefficient is defined as κ total = κ N +. The, transmission spectrums are obtained by solving the coupled-mode equations via transfer matrix method. The power is obtained by subtracting the and transmission power from the total power. 3. NUMERICAL RESULTS The effect of coupling in second order FBG with theoretically rectangular index modification is concentrated on. The designed Bragg wavelength is set at 1.55 µm. r (z)...... z z Radiation Coupling Coefficient (1/m) 5 25-25 -5 =.1 Designed λb:1.55um..2.4.6.8 1. 6 4 2-2 -4-6 Index Coupling Coefficient 1/m Figure 1: Schematic diagram of FBG with an arbitrary profile. Figure 2: Dependences of coupling coefficients on duty cycle. The effect of duty cycle is firstly studied. Seen from Fig. 2, the index coupling coefficient sinusoidally varies, obtaining the maximum values when duty cycle equals.25 and.75, equaling zero and changing sign when duty cycle equals.5. The coupling coefficient varies with the square of sinusoidal function, obtaining the maximum value when duty cycle equals.5. These are all determined by the corresponding analytical expressions of coupling coefficients and the Fourier expression properties of the periodic rectangular function. It is thus not difficult to understand the wave becomes prominent with the decrease of the power only when the duty cycle is about.5, seen from Fig. 3. The prominent occurrence of the wave corresponds to the narrowest 3 db bandwidth of the Bragg peak seen from Fig. 4. It also leads to a blue shift of the Bragg wavelength, seen from Fig. 5. Conclusively, at the special case of duty cycle around.5, the dominant complex coupling results in a decrease of the total coupling coefficient and offer a possibility to control the total coupling coefficient in a wider range, which thus leads to some properties different from the first order FBG, such as a lower power, a narrower and blue-shifted Bragg peak. 1.1.1 1E-3 1E-4 1E-5 1E-6 Radiation =.1 Grating Length:2mm Designed λ B :1.55um..2.4.6.8 1. 3.5 rectangle, =.1,Grating Length:2mm 3. Designed λ B :1.55um 2.5 2. 1.5 1..5...2.4.6.8 1. Figure 3: Dependence of and power on duty cycle. Figure 4: Dependence of 3 db bandwidth on duty cycle. Moreover, taking the advantage of other grating parameters, the effect of coupling mentioned above can be further controlled. Seen from Figs. 6 8, due the enhancement of coupling with the increase of the index modification, the power will increase and the Bragg peak will be broader. On the other hand, the increase of the grating length will
PIERS ONLINE, VOL. 3, NO. 4, 27 465 result in the increase of the power and the decrease of the Bragg bandwidth, as shown in Figs. 9 1. 1 transmission =.1 1-1 Grating Length:2mm =.5 1-2 Designed λ B :1.55µm 1-3 1.5496 1.5498 1.55 1.552 1.554 λ (um) Radiation Coupling Coefficient (1/m) -5-1 =.5 Designed λ B =1.55 m Real part Imaginary part -15..1.2.3.4.5 Figure 5: Spectrums at duty cycle equaling.5. Figure 6: Dependences of coupling coefficients on index change. 1 1-1 1-2 1-3 1-4 =.5 Grating Length:2mm Designed λ B :1.55um transmission..1.2.3.4.5.6.5.4.3.2.1. rectangle =.5 Designed λ B =1.55µm Grating Length:2mm..1.2.3.4.5 Figure 7: Dependence of powers on index change. Figure 8: Dependence of 3dB bandwidth on index change. 1 1-1 1-2 transmission =.1 =.5 Designed λ B :1.55µm 1 1 1 1 2 Grating Length (mm).8.7.6.5.4.3.2.1. =.1 =.5 Designed λ B =1.55µm 1 1 1 1 2 Grating Length (mm) Figure 9: Dependence of powers on grating length. Figure 1: Dependence of 3 db bandwidth on index change. 4. CONCLUSION By using Coupled mode theory and Green s function method, the effect of coupling in higher order FBG is modeled and simulated. For a second order FBG with rectangular index modification, when the duty cycle is around.5, the wave becomes prominent. Meanwhile, with the elimination of the usually dominant index coupling, the complex coupling is dominant in the total coupling coefficient. It results in the decrease and a wider controllable range for the total coupling coefficient, and then the decrease of power and Bragg bandwidth. These properties are simultaneously influenced by the index change and grating length. The study will be useful for the applications of high order FBGs, such as the design of DFB fiber laser.
PIERS ONLINE, VOL. 3, NO. 4, 27 466 ACKNOWLEDGMENT Thank the Chinese Academy of Sciences for providing Li YANG with the fellowship and the Ontario Photonics Consortium for financial support. REFERENCES 1. Xie, W. X., et al., Second order diffraction efficiency of Bragg gratings written within germanosilicate fibres, Opt. Commun., Vol. 11, No. 1 2, 85 91, 1993. 2. Jia, H. Z. and Y. L. Li, First- and second-order diffraction characteristics of fiber Bragg gratings, Opt. Commun., Vol. 178, 339 343, 2. 3. Soltani, T., Second order Bragg grating in slab waveguide and optical fiber, Master thesis, McMaster University, 2. 4. Echevarria, J., et al., Uniform fiber Bragg grating first-and second-order diffraction wavelength experimental characterization for strain-temperature discrimination, IEEE, Phot. Tech. Lett., Vol. 13, No. 7, 696 698, 21. 5. Carballar, A., et al., Growth modeling of fiber gratings: a numerical investigation, Fiber and Integrated Optics, Vol. 21, 451 461, 24. 6. Hardy, A., et al., Analysis of second-order gratings, IEEE J. QE, Vol. 25, No. 1, 296 215, 1989. 7. Shams-Zadeh-Amiri, A. M., et al., Second- and higher order resonant gratings with gain or loss-part I: Greens function analysis, IEEE J. QE, Vol. 36, No. 12, 1421 143, 2.