Nonlinear cone optical fiber —
propagation characteristic
Ad a m Ma je w sk i, Ar t u r Ka r c z e w sk i
Warsaw University of Technology, Institute o f Electronic Systems, u l Nowowiejska 15/19, 0 0 — 665 Warszawa, Poland.
Numerical analysis of the nonlinear cone optical fiber using the coupled mode theory and the split-step Fourier method is described. The nonlinear cone optical fiber is divided into a large number o f appropriate cylindrical segments. In the analysis of light propagation between two segments of various diameters, the coupled mode theory is used. Propagation of a pulse, in the cylindrical nonlinear segment, is analyzed using the split-step Fourier method. Some results of computation for the selected cone conlinear optical fibers are also presented.
1. Introduction
A cone waveguide is a guiding structure for electromagnetic waves that gradually becomes narrower towards one of the ends [1]. The analysis of such guides can be carried out in two steps. In the first one, the fiber is approximated by a large number of segments of appropriate diameter (the staircase model of the optical fiber), Fig. 1.
1
Fig. 1. Staircase model of optical fiber
The field distribution on the boundary of two segments of different diameter is described using the coupled mode theory [2] — [6]. In the case of the staircase of the nonlinear optical fiber (Fig. I), the exchange of power between two segments is described by the equation [3]
A
2W ) =
v/p? > № 12|¿¡?W )
where E0 is the electric field, # 12 is the amplitude coupling coefficient
(1)
JJ
E ^ E ^ d x d y <P12J l
\ E ^ \ 2dxdy JJ \EP\2dxdyÊn is the normalized electric field
(2)
(3)
where P0 represents the power of the mode in the segment
In the second step, the nonlinear Schrodinger equation (NLSE) of the form [7],
[
8]
j8^ - W a 4 +y]A]lA = 0
(4)*NL'
is solved for each segment. In Eq. (4), A represents the envelope function, y
*cf j2 o
denotes the nonlinearity parameter, = — =■ is the dispersion parameter of the
dco
fiber, is the nonlinear refractive coefficient, which for S i0 2 is equal to 1.2* 10“ 22 m2/V2, k = A ti is the mode effective area [7], [8], /? is the phase
A
constant, X is the wavelength. Equation (4) is valid for pulses wider than 100 fs. If a
losses are taken into consideration, the component of j - A must be added to
JL
Eq. (4), where a represents attenuation.
The idea of the split-step Fourier method is based on the mutual physical interaction of nonlinearity and dispersion. The solution of the NLSE by the split-step Fourier method is carried out in two steps. In the first one, the influence of nonlinearity on the pulse shape is analysed, in the second step, it is just the dispersion effect that is analysed. After some transformations, Eq. (4) takes the form
[9] ^ = {D + N )A oz where d2 a ST2 2 (5) (
6
)is the differential operator which describes the dispersion and the losses in a nonlinear guide,
N = jy \A \2 (7)
is a nonlinear differential operator describing the nonlinear effects occurring during the propagation of the optical pulses. After some mathematics, we obtain
A (z+ h,T) = exp(hD) exp(hN)A(z, T). (8)
To implement the method, the guide is divided into a large number of segments of width h. The pulse is carried twice through each segment [9], [10], once as a nonlinear segment, where D = 0, and then as a dispersive segment, where TV = 0 (Fig. 2\
N o n lin e a r ity se g m e n t
D isp ersion se g m e n t
Fig. 2. Scheme of calculations
In the case of the nonlinear segment, Eq. (8) has the form
A i (z+ hyT) = exp(hN)A(z,T). (9)
The exponential operator exp(hD), which represents the dispersive segment, is convenient for calculations in the frequency domain
A (z+ h,T) = F -^expthD ticD K F lA^z+ h'T)-]} (10)
where F denotes the Fourier transform operation, D(jco) is given by Eq. (6) when the differential operator d/dT is replaced by (jco)y namely
0 = 7^ 2 - “ · (11)
2. N um erical results
To prove how the nonlinearity y and dispersion p2 vary, calculations of the cone nonlinear optical fiber for shifted dispersion characteristics have been carried out. Results of the computation for A = 0.9%, where A = ——— is the relative
n2
core-cladding difference, have been presented (Figs. 3 — 5).
Fig. 3. Amplitude o f the electric field vs core radius (A = 0.9%, X = 1.S5 pm, r [pm])
Fig. 5. Nonlinearity y vs. core radius (d = 0 .9 % , X = 1.55 pm, r [pm ], y [W -1 m " 1] )
Fig. 6. Dispersion P2 vs. core radius (d = 0 .9 % , X = 1.55 pm, r [pm ], P 2 [ps2/k m ])
Using the nonlinearity y and dispersion /?2, calculation of the electric field distribution, for shifted dispersion characteristic, has been carried out (Fig. 6). In Figure 7, fundamental solitons in the cylindrical and the cone nonlinear optical fiber are presented, where A = 0.95%.
3. Conclusions
Owing to the cone structure of the optical guide the values of nonlinear and dispersion parameters are influenced in a such way that dispersion due to losses is partially compensated. As the diameter of the guide is getting smaller the dispersion diminishes and the intensity of light becomes greater in respect to its initial
t/t0
Fig. 7. Fundamental soliton in the cone nonlinear optical fiber {A = 0.95%, A = 1.55 pm, a = 0.3 dB/km, z = 5 km). Dashed line — input shape. Dotted line — the tapered nonlinear optical fiber where initial value o f radius a is 2.27 pm, final value of a is 2.17 pm. Solid fine — the nonlinear optical fiber where a = 2 2 1 pm
magnitude, so mutual compensation of the nonlinear and dispersion effects is preserved over longer distance than in cylindrical structure. Both numerical methods, coupled mode theory and split-step Fourier method, have proved to be successful; high accuracy (better than 1%) is achieved when the segments are 2 m long. The staircase model of the cone guide enables a simulation of imperferctions, which may occur due to technological process.
References
[1] Spo r l e d e r F., Un g e r H. G , IEE Electromagnet Wave Series 6 (1979). [2] Sn y d e r A. W., IEEE Trans. Microwave Theory Tech. 18 (1970), 383.
[3] Eb e l in g K . J., In teg rierte O pto elek tro n ik , W ellenleiteroptik, P h oton ik, H a lb leiter, Springer-Verlag, Berlin 1989.
[4] Yajima H , IEEE J. Quantum Electron. 14 (1978), 749.
[5] Petykiewicz J., In te g ra te d O p tics: P h ysica l F un dam entals, (in Polish), [Ed.] WNT, Warszawa 1989. [6] Ma r c u s e D., T h e o ry o f D ielectric W aveguides, Academic Press, New York 1974.
[7] Ma jew ski A . S oliton s in O p tic a l G uides, (in Polish), [Ed.] Wydawnictwa Politechniki Warszawskiej, Warszawa 1991.
[ 8 ] Ma jew sk i A , N on lin ear F ib er O ptics, (in Polish), [Ed.] Wydawnictwa Politechniki Warszawskiej, Warszawa 1993.
[9] Ag r a w a l G. P., N o n lin ea r F ib e r O ptics, Academic Press 1989.
[10] Ha g e sa w a A., Ta p p e r t F., Appl. Phys. L ett 23 (1973), 142.