# Modeling Nonlinear Acoustooptic Coupling in Fiber Optics Based on Refractive Index Variation due to Local Bending

Authors: Catalina Hurtado Castano; Rodrigo Acuna Herrera; Pedro I. Torres
Journal: Advances in Optical Technologies
 (2016)

---

## Abstract

A detailed procedure is presented to compute analytically the acoustooptic coupling coefficient between copropagating core and lowest-order cladding modes in tapered fiber optics. Based on the effect of the local bending, the linear and nonlinear variations in the refractive index are modeled. A set of equations and parameters are presented in order to calculate the influence of acoustooptic effect in nonlinear pulse propagation. We will show that as the tapered fiber diameter decreases more energy can be transferred to the cladding and the nonlinear phenomena can compensate the coupling coefficients effects.


---

## Body

## 1. Introduction

The use of light to solve specific problems has led to the development of new knowledge in materials and optical systems which have found applications in technologies such as photonics [1]. The development of fiber optics (FO) has increased the capacity of processing light by offering systems with ultra high speed such as telecommunication networks and fiber laser [2, 3]. Furthermore, FO can be used as a multiplexers, modulators, and tunable filters based on acoustooptic (AO) coupling effect.Although the AO coupling mechanism is well known, the study and design of AO devices remain as a topic that deserves attention, especially if one wants to achieve high frequency processing and reproducibility. Studies of AO effect in fibers have been carried out both theoretically and experimentally, where acoustic waves are considered long period fiber grating (LPFG) which are capable of coupling coherent copropagating modes traveling in neighboring wave guide that interact through the evanescent optical field [4]. The phase matching condition between the modes needs the overlapping among the interacting modes, allowing energy transfer between them. LPFGs offer the possibility of nonlinear coupling interaction between copropagating guided modes [5]. LPFGs have been manufactured with UV lasers to study the propagation of nonlinear pulse in coupled modes [6].Herein, using a local bending model, variations in the refractive index due to the presence of the acoustic wave are calculated. A detailed theory and methodology are presented in order to calculate analytically the nonlinear coupling pulse propagation between copropagating core and linearly polarized (LP) cladding modes in tapered fibers. This paper is structured in the following way: Section2 presents the fundamental theory of local bending model for linear acoustooptic coupling. Section 3 describes the expressions to calculate field profiles and propagation constants of cladding modes. In Section 4, based on local bending, novel equations for the nonlinear coupling between copropagating core-cladding modes are presented. Section 5 presents the results of numerical simulations. Finally, the conclusions are presented in last section.


## 2. Local Bending Model for AO Coupling

In this section, the AO coupling interactions between the fundamental core mode and copropagating cladding modes are analyzed. We see that acoustooptic coupling equations described changes in the amplitude of the mode propagated through an optical fiber, either in core region or in cladding region.The transference of energy between fiber modes can be achieved by the presence of acoustic waves [7]; we assume that a single mode fiber is modulated in the x
-
z plane and the modulation locally bends it with radius of R
(
z
). The effect of the local bending, as a consequence of acoustic wave, can be modeled as a variation in the refractive index [8] of the fiber which is given by the expression(1)
n
2
=
n
c
o
2
1
-
1
+
χ
2
x
R
z
,where n
c
o is the core refractive index for the fundamental fiber mode without bending, x is the local coordinate on the straight fiber, and χ is the silica elastooptic coefficient whose value is -
0.22. Assuming that the amplitude of oscillation is much smaller than the acoustic wavelength, one can approximate the local curvature κ
=
1
/
R
(
z
) as(2)
κ
≈
A
Ω
2
cos
⁡
Ω
z
,where A is acoustic amplitude, Ω
=
2
π
/
Λ is the acoustic wave number, and Λ is acoustic wave period. The coupled-mode equations that describe linear coupling interactions in a long-period grating are(3)
∂
a
c
o
∂
z
=
i
C
c
o
-
c
l
e
i
β
0
ν
z
a
ν
c
l
,
(4)
∂
a
ν
c
l
∂
z
=
i
C
c
l
-
c
o
e
i
β
0
ν
z
a
c
o
,where a
c
o is the core mode amplitude, a
ν
c
l the cladding ν-mode amplitude, and β
0
ν the detuning parameter. In (4) the longitudinal part of the coupling coefficient is neglected since it is substantially much more smaller than the transverse part [9]. The maximum transference of energy is given at the phase condition β
0
ν
(
λ
r
)
=
0, where λ
r is the resonant wavelength. Then, we can write(5)
β
0
ν
λ
r
=
β
c
o
λ
r
-
β
ν
c
l
λ
r
-
Ω
.The simplified expression of the coupling coefficient between the fundamental core mode and any cladding mode is given byC
μ
ν and described in [8] for local bending model; an extensive expression can be seen in [10, Eq. (
34
)]. In this paper, we consider the interaction between the fundamental core mode (L
P
01
c
o) and the first three lowest cladding modes (L
P
1
ν
c
l, ν
=
1,2
,
3) since C
μ
ν is only nonzero for modes that differ by an odd integer in angular momentum [8]. Using the local bending model, C
μ
ν is written as(6)
C
μ
ν
=
A
k
o
2
n
μ
2
Ω
2
1
+
χ
2
β
μ
ξ
μ
t
x
ξ
υ
t
,where subscripts μ and υ reference the core and cladding mode, respectively; ξ
μ
t is the transverse mode profile, normalized according to ξ
μ
t
∣
ξ
υ
t
=
δ
μ
υ. The normalized power carried by the L
P
01
c
o and L
P
1
ν
c
l modes can be written as a function of the propagation distance z as(7)
P
c
o
z
P
c
o
0
=
c
o
s
2
γ
0
ν
z
+
β
0
ν
2
γ
0
ν
2
s
i
n
2
γ
0
ν
z
,
(8)
P
ν
c
l
z
P
c
o
0
=
C
0
ν
γ
0
ν
2
s
i
n
2
γ
0
ν
z
,where γ
0
ν
2
=
β
0
ν
/
2
2
+
C
2, C
2
=
C
0
ν
C
ν
0, and P
ν
c
l
=
|
a
ν
c
l
|
2. The initial power is pumped through the fundamental core region P
c
o
(
0
)
=
|
a
c
o
|
2, and the initial power is pumped through the cladding region P
ν
c
l
(
0
)
=
0.


## 3. LP Cladding Modes

This section is devoted to describe the expressions to calculate field profiles and propagation constants of cladding modes using exact analytical equations for a step index fiber optics.The cladding is considered a multimode wave guide; for this to exist, it is necessary to remove the polymer coating to satisfy the condition of total internal reflection at the external boundary [11]; therefore these modes are very sensitive to change because of the surrounding medium. Cladding modes can be excited by using LPFGs which can couple light from core to cladding modes [11]. This happens due to the fact that the condition of total internal reflection is violated when the optical fiber is bending.We are interested mainly in the interaction between the fundamental core mode (L
P
01) and the lowest (L
P
1
ν) cladding modes in a tapered fiber optics. In the case of the mode guided by the fiber core, field profile and propagation constant are calculated using the expressions given in [12]. For cladding modes, we used the dispersion expression given in [10, 13] for the three-layer fiber model in order to compute the constant propagations. The cladding mode dispersion equation is given by the equations(9)
ζ
0
=
ζ
0
′
,where(10)
ζ
0
=
1
σ
2
u
2
J
K
+
σ
1
σ
2
u
21
u
32
/
n
2
2
a
1
a
2
p
l
a
2
-
1
/
u
2
s
l
a
2
-
K
q
l
a
2
+
J
r
l
a
2
-
u
2
u
32
/
n
2
2
a
2
J
-
u
21
/
n
1
2
a
1
K
p
l
a
2
-
u
32
/
n
1
2
a
2
q
l
a
2
+
u
21
/
n
1
2
a
1
r
l
a
2
,
ζ
0
′
=
σ
1
u
2
u
32
/
a
2
J
-
n
3
2
u
21
/
n
2
2
a
1
K
p
l
a
2
-
u
32
/
a
2
q
l
a
2
-
u
21
/
a
1
r
l
a
2
u
2
n
3
2
/
n
2
2
J
K
+
σ
1
σ
2
u
21
u
32
/
n
1
2
a
1
a
2
p
l
a
2
-
n
3
2
/
n
1
2
K
q
l
a
2
+
J
r
l
a
2
-
n
2
2
/
n
1
2
u
2
s
l
a
2
.The variables and functions given in (10) are the following:(11)
σ
1
=
i
l
n
e
f
f
Z
0
,
σ
2
=
i
l
n
e
f
f
Z
0
,
u
21
=
1
u
2
2
-
1
u
1
2
,
u
32
=
1
w
3
2
-
1
u
2
2
,
u
j
2
=
2
π
λ
2
n
j
2
-
n
e
f
f
2
j
∈
1,2
,
w
3
2
=
2
π
λ
2
n
e
f
f
2
-
n
3
2
,
J
=
J
l
′
u
1
a
1
u
1
J
l
u
1
a
1
,
K
=
K
l
′
w
3
a
2
w
3
K
l
w
3
a
2
,
p
l
r
=
J
l
u
2
r
N
l
u
2
a
1
-
J
l
u
2
a
1
N
l
u
2
r
,
q
l
r
=
J
l
u
2
r
N
l
′
u
2
a
1
-
J
l
′
u
2
a
1
N
l
u
2
r
,
r
l
r
=
J
l
′
u
2
r
N
l
u
2
a
1
-
J
l
u
2
a
1
N
l
′
u
2
r
,
s
l
r
=
J
l
′
u
2
r
N
l
′
u
2
a
1
-
J
l
′
u
2
a
1
N
l
′
u
2
r
,where a
1 and a
2 are the core and cladding radius, respectively. l is the azimuthal number. J
l
(
x
) denoted the Bessel functions of the first kind, N
l
(
x
) are the Bessel functions of the second kind, and K
l
(
x
) are the modified Bessel functions of the second kind; also, J
l
′
(
x
), N
l
′
(
x
), and K
l
′
(
x
) indicate differentiation of Bessel functions with respect to the total argument x. Based on the three-layer model, the electric field components for cladding modes are [10, 14]ρ
≤
a
1
(12)
E
ρ
c
l
=
E
l
ν
c
l
u
1
2
J
l
+
1
u
1
ρ
+
J
l
-
1
u
1
ρ
-
σ
2
ζ
0
n
1
2
J
l
+
1
u
1
ρ
-
J
l
-
1
u
1
ρ
cos
⁡
l
ϕ
,
E
θ
c
l
=
E
l
ν
c
l
u
1
2
J
l
+
1
u
1
ρ
-
J
l
-
1
u
1
ρ
-
σ
2
ζ
0
n
1
2
J
l
+
1
u
1
ρ
+
J
l
-
1
u
1
ρ
sin
⁡
l
ϕ
,
E
z
c
l
=
i
E
l
ν
c
l
π
a
1
u
1
2
u
2
2
σ
2
J
l
u
1
ρ
2
n
2
2
β
l
G
2
p
l
ρ
-
n
2
2
ξ
0
n
1
2
u
2
q
l
ρ
cos
⁡
ν
ϕ
,
F
2
=
J
-
u
21
σ
2
ζ
0
n
1
2
a
1
,
G
2
=
ζ
0
J
+
u
21
σ
a
1
.TheH
E
2
,
ν
c
l conventional mode is equivalent to linearly polarized mode L
P
1
ν
c
l [12], which is our mode of interest. The cladding modes occur in a strict alternating sequence of H
E and E
H, except from the fundamental mode H
E
1,1
c
l [13]. In our case, we set l
=
2 and ν
=
1,3
,
5 in conventional modes which are the first three lowest cladding modes which are the ones we are interested in. Then, Figure 1 shows the E
x-profile for these modes in a standard tapered single mode fiber optics.Figure 1
E
x-profiles for the first three lowest cladding modes L
P
11
c
l, L
P
12
c
l, and L
P
13
c
l in a standard tapered single mode fiber optics.
(a)
ModeL
P
11
c
l
(b)
ModeL
P
12
c
l
(c)
ModeL
P
13
c
l


## 4. Nonlinear Pulse Propagation in Coupled Copropagating Modes

Nonlinear pulse propagation in a linear periodic medium has been studied using coupled-mode theory for the copropagating core and cladding modes in a LPFG [6, 15], where it is necessary to include both the nonlinear effect and the periodic variation of the refractive index by using the equation(13)
n
ω
,
z
=
n
o
+
δ
n
g
z
+
η
2
E
2
=
n
L
+
n
N
L
,where n
o is the nonperturbed linear refractive index, η
2 is the silica nonlinear refractive index, and δ
n
g
(
z
) stands for periodic index variations as a consequence of the grating. If we rewrite the linear part of the refractive index n
L
2
(
ω
,
z
)
=
1
+
f
(
x
,
y
,
z
), where f
(
x
,
y
,
z
) accounts for the refractive index variation associated with the long period grating, one can show that(14)
f
x
,
y
,
z
≃
n
o
2
-
1
×
1
+
2
n
o
1
+
χ
A
x
Ω
2
cos
⁡
Ω
z
.Taking (14), using Maxwell’s equations, and writing the electric field as a sum of the fundamental core and a LP cladding mode in a fiber optic, one can obtain the nonlinear amplitudes evolution [6, 16]:(15)
∂
A
1
∂
z
-
i
∑
n
=
2
∞
i
n
β
1
n
n
!
∂
n
A
1
∂
T
n
-
i
C
12
A
2
-
i
δ
2
A
1
=
i
d
11
A
1
2
+
2
d
12
A
2
2
A
1
,
∂
A
2
∂
z
+
V
g
∂
A
2
∂
T
-
i
∑
n
=
2
∞
i
n
β
2
n
n
!
∂
n
A
2
∂
T
n
-
i
C
12
A
1
+
i
δ
2
A
2
=
i
d
22
A
2
2
+
2
d
12
A
1
2
A
2
;here, A
1 is the amplitude of fundamental core mode and A
2 the amplitude of a LP cladding mode (L
P
11
c
l, L
P
12
c
l, or L
P
13
c
l). The equation β
m
n
=
∂
n
β
m
/
∂
T
n
|
ω
o is n-order of the Taylor expansion series of β
m (m
=
1,2), where ω
o is the pump frequency; V
g
=
β
21
-
β
11 and T
=
t
-
β
11
z. C
12 is the coupling coefficient given in (6) and δ
(
λ
)
=
β
0
c
o
(
λ
)
-
β
ν
c
l
(
λ
)
-
Ω is the detuning, which is the same equation given in (5). The acoustic wave number (
Ω
) and acoustic wave period (
Λ
) was described in Section 2. d
i
js are the overlapped field coefficients associated with the nonlinear phenomena, given by(16)
d
i
j
=
2
π
η
2
λ
A
e
f
f
E
t
,
i
2
∣
E
t
,
j
2
E
t
,
1
2
∣
E
t
,
1
2
,
E
t
,
i is the normalized transverse electric field for the i mode and A
e
f
f is effective fundamental core mode area defined as [16](17)
A
e
f
f
=
E
t
,
1
∣
E
t
,
1
2
E
t
,
1
2
∣
E
t
,
1
2
.The coefficientsd
11 and d
22 are related to self-phase modulation (SPM) for L
P
01
c
o and L
P
1
ν
c
l, respectively. d
12 accounts for cross-phase modulation (XPM) [16].


## 5. Numerical Simulations

The interactions considered in this paper occur between the fundamental core modeL
P
01
c
o and the first three lowest cladding modes L
P
1
m
c
l, with m
=
1,2
,
3. A tapered standard single mode fiber is used for the calculations. Two tapered core diameters are taken into account; they are d
c
o
=
2.5
μ
m and d
c
o
=
3.5
μ
m. Fiber referenceSMF ITU-T. G.652.B was used. The cladding diameter is reduced in the same proportion as the core diameters. Additionally, we set Δ
=
(
n
c
o
-
n
c
l
)
/
n
c
o
=
0.0036, where n
c
o stand for the core refractive index and n
c
l for the cladding. Sellmeier equation for silica is included.


### 5.1. Linear Simulations

To computeC
12 (Eq. (6)) it is necessary to know Λ that can be established by the dispersion relationship [8, 17]:(18)
f
a
=
π
C
e
x
t
d
c
l
/
2
Λ
2
,where f
a is the acoustic wave frequency, C
e
x
t
=
5760
m
/
s is the speed of the acoustic wave in silica, and d
c
l is the cladding diameter. It is assumed that axial strain is not applied on the fiber. Figure 2 shows f
a as a function of wavelength at δ
(
λ
r
)
=
β
0
ν
(
λ
r
)
=
0 described in (5) and coupling between the fundamental core mode and the L
P
1
ν
c
ls modes for two tapered fiber diameters. Figure 2 allows us to identify which frequency needs the fiber to oscillate to get a perfect coupling among the modes for a specific wavelength (λ). We see that no high AO frequency is required to transfer efficiently energy from the core to cladding modes when λ shifts to infrared. Additionally, the smaller the tapered fiber diameter is, the more the required AO frequency decreases, which is a good feature because it could be difficult to find AO drivers at high frequencies.Figure 2
AO frequencyf
a as a function of wavelength at δ
(
λ
)
=
β
0
ν
(
λ
)
=
0 and coupling between the fundamental core mode and the L
P
1
ν
c
ls modes for (a) d
c
o
=
2.5
μ
m and (b) d
c
o
=
3.5
μ
m.
(a)
d
c
o
=
2.5
μ
m
(b)
d
c
o
=
3.5
μ
mFixing AO frequency atf
a
=
1.0
M
H
z, we compute the detuning δ
(
λ
) (5). The results can be seen in Figure 3. As the tapered diameter decreases the detuning δ is blue-shifted. δ-zero crossing moves toward blue side around 250
n
m per 1
μ
m. At fixed wavelength, δ increases for higher cladding order modes. From Figures 2 and 3, one can conclude that, by changing the tapered fiber diameter, AO frequency and detuning can be easily controlled, which makes this type of fiber a good candidate to design AO fiber filter in a large range of wavelengths.Figure 3
Detuning parameterδ as a function of wavelength for coupling between the fundamental core mode and the L
P
1
ν
c
ls modes when (a) d
c
o
=
2.5
μ
m and (b) d
c
o
=
3.5
μ
m.
(a)
d
c
o
=
2.5
μ
m
(b)
d
c
o
=
3.5
μ
mUsing (5) and (6), the resulting normalized transmitted power equation (7) in the core mode at the end of the fiber is shown in Figure 4 [18–20], where P
c
o
(
L
)
/
P
c
o
(
0
) in dB units is plotted. It is set: fiber length of L
=
15
c
m, AO wave amplitude of A
=
63
n
m, and f
a
=
1
M
H
z. According to (18), Λ
=
587.20
μ
m for d
c
o
=
2.5
μ
m and Λ
=
694.79
μ
m for d
c
o
=
3.5
μ
m. The coupling between L
P
01
c
o and L
P
12
c
l exhibits the biggest deep which means that the incident light in the core is partially transferred to the cladding mode at λ
=
680
n
m for d
c
o
=
2.5
μ
m. Similar behavior happens for d
c
o
=
3.5
μ
m at λ
=
923
n
m.Figure 4
Normalized transmitted power in the core mode at the end of the fiber. Setting: fiber length ofL
=
15
c
m, AO wave amplitude of A
=
63
n
m, f
a
=
1
M
H
z, d
c
o
=
2.5
μ
m, and d
c
o
=
3.5
μ
m.


### 5.2. Nonlinear Simulations

For nonlinear pulse simulations, (15) are numerically solved by using the split-step Fourier method [16]. Fiber length of 15
c
m and initial pulse width of T
o
=
50
f
s are considered. We analyze two scenarios: (i) d
c
o
=
2.5
μ
m: pump wavelength λ
o
=
1323
n
m, initial peak power P
o
=
30
k
W, and (ii) d
c
o
=
3.5
μ
m: pump wavelength λ
o
=
1472
n
m, initial peak power P
o
=
3
k
W. We choose the above pump wavelengths because we want to generate solitons at both diameters which allow us to get more insight into nonlinear effects.The linear coupling parameterC
12
(
λ
) is shown in graph in Figure 5 for the two diameters under consideration. C
12 almost scales linearly with wavelength, except for L
P
01
c
o-L
P
11
c
l at d
c
o
=
2.5
μ
m which presents some oscillation. The coupling between L
P
01
c
o and L
P
13
c
l exhibits the weakest strength at both cores diameters since the mode profiles do not overlap too much.Figure 5
AO coupling coefficientC
12 as a function of wavelength for the coupling between L
P
01
c
o and L
P
11
c
l, L
P
12
c
l, and L
P
13
c
l in a tapered standard fiber. Two core diameters are considered, d
c
o
=
2.5
μ
m and d
c
o
=
3.5
μ
m.The plots in Figures6, 7, and 8 present the overlap nonlinear field coefficients as a function of wavelength between the fundamental core mode and the L
P’s cladding modes. In Figure 6, we can see that this parameter that governs the SPM phenomenon in the core region is superior for tapered fiber with greater cross-sectional area for most of the analyzed range of wavelength. This could be due to the fact that more energy goes to the cladding as the core diameter is reduced. From Figure 7, d
22 almost linearly decreases as wavelength increases and has the same value for any L
P’s cladding mode in the same diameter. d
22 has small value since most of the light keeps in the single mode core, which makes SPM not efficiently excited at low power.Figure 6
Overlap field nonlinear coefficientd
11 as a function of wavelength in a tapered standard fiber. Two core diameters are considered, d
c
o
=
2.5
μ
m and d
c
o
=
3.5
μ
m.Figure 7
Overlap field nonlinear coefficientd
22 as a function of wavelength in a tapered standard fiber. Two core diameters are considered, (a) d
c
o
=
2.5
μ
m and (b) d
c
o
=
3.5
μ
m.
(a)
d
c
o
=
2.5
μ
m
(b)
d
c
o
=
3.5
μ
mFigure 8
Overlap field nonlinear coefficientd
12 as a function of wavelength in a tapered standard fiber. Two core diameters are considered, (a) d
c
o
=
2.5
μ
m and (b) d
c
o
=
3.5
μ
m.
(a)
d
c
o
=
2.5
μ
m
(b)
d
c
o
=
3.5
μ
mThe behavior of the overlap coefficientsd
12 is shown in Figure 8; they also decrease for higher wavelengths. The coefficient L
P
01
c
o-L
P
13
c
l is the strongest one, which means that more energy can be nonlinearly transferred from the core to the cladding weather; L
P
13
c
l is efficiently coupled to L
P
01
c
o through the acoustic wave.To complete the needed parameters to simulate (15), the dispersion coefficients β’s are calculated. These can be seen in Table 1.Table 1
Dispersion coefficientsL
P
01
c
o and LP cladding modes.
(a)
d
c
o
=
2.5 μm and λ
0
=
1323 nm
Mode
β
2 [ps2m−1]
β
3 [ps3m−1]
β
4 [ps4m−1]
β
5 [ps5m−1]
L
P
01
c
o
-
1.81
×
1
0
-
4
1.05
×
1
0
-
4
-
1.66
×
1
0
-
7
1.58
×
1
0
-
10
L
P
11
c
l
-
4.67
×
1
0
-
3
7.39
×
1
0
-
5
-
9.56
×
1
0
-
8
9.69
×
1
0
-
11
L
P
12
c
l
-
5.66
×
1
0
-
3
7.57
×
1
0
-
5
-
9.85
×
1
0
-
8
9.98
×
1
0
-
11
L
P
13
c
l
-
6.94
×
1
0
-
3
7.78
×
1
0
-
5
-
1.01
×
1
0
-
7
1.02
×
1
0
-
10
(b)
d
c
o
=
3.5 μm and λ
0
=
1472 nm
Mode
β
2 [ps2m−1]
β
3 [ps3m−1]
β
4 [ps4m−1]
β
5 [ps5m−1]
L
P
01
c
o
-
3.01
×
1
0
-
4
1.01
×
1
0
-
4
-
1.89
×
1
0
-
7
2.06
×
1
0
-
10
L
P
11
c
l
-
1.86
×
1
0
-
2
9.94
×
1
0
-
5
-
1.39
×
1
0
-
7
1.33
×
1
0
-
10
L
P
12
c
l
-
2.00
×
1
0
-
2
1.02
×
1
0
-
4
-
1.43
×
1
0
-
7
1.38
×
1
0
-
10
L
P
13
c
l
-
2.18
×
1
0
-
2
1.06
×
1
0
-
4
-
1.50
×
1
0
-
7
1.44
×
1
0
-
10Using (15), Figures 9 and 10 display the spectral response of all modes under consideration for both tapered fiber diameters. They evidence the generation of Supercontinuum (SC) both at the core and at the cladding. One can see that SC broadens more at the cladding modes when d
c
o
=
3.5
μ
m, even with less power compared to the d
c
o
=
2.5
μ
m case. This makes senses because d
12 is larger for d
c
o
=
3.5
μ
m (see Figure 8).Figure 9
Nonlinear spectral response ford
c
o
=
2.5
μ
m, with pump power of P
0
=
30
k
W, initial width of T
0
=
50
f
s, and pump wavelength of λ
0
=
1323
n
m.Figure 10
Nonlinear spectral response ford
c
o
=
3.5
μ
m, with pump power of P
0
=
3
k
W, initial width of T
0
=
50
f
s, and pump wavelength of λ
0
=
1472
n
m.As last step, we analyze how much energy can be transferred from the core to the LP cladding modes as a function of acoustic frequencyf
a for different pulse width T
0
=
50
f
s, 100
f
s, and 150
f
s; the results are plotted in Figure 11. The energy is normalized according to E
=
P
T
/
P
0
T
0, where P
T is computed as ∫
-
∞
∞
‍
P
(
T
)
d
T. From Figure 11, the transferred energy behaves linearly increasing as f
a and T
0 increase. This shows that although C
12 is proportional to f
a
2, the nonlinear effects linearize the energy and the slope can be controlled by T
0.Figure 11
Normalized transferred energy from the core to the LP cladding modes as a function of acoustic frequencyf
a for different pulse width T
0
=
50
f
s, 100 fs, and 150 fs. (a) For core diameter 2.5 μm and (b) core diameter 3.5 μm.
(a)
d
c
o
=
2.5
μ
m
(b)
d
c
o
=
3.5
μ
m


## 5.1. Linear Simulations

To computeC
12 (Eq. (6)) it is necessary to know Λ that can be established by the dispersion relationship [8, 17]:(18)
f
a
=
π
C
e
x
t
d
c
l
/
2
Λ
2
,where f
a is the acoustic wave frequency, C
e
x
t
=
5760
m
/
s is the speed of the acoustic wave in silica, and d
c
l is the cladding diameter. It is assumed that axial strain is not applied on the fiber. Figure 2 shows f
a as a function of wavelength at δ
(
λ
r
)
=
β
0
ν
(
λ
r
)
=
0 described in (5) and coupling between the fundamental core mode and the L
P
1
ν
c
ls modes for two tapered fiber diameters. Figure 2 allows us to identify which frequency needs the fiber to oscillate to get a perfect coupling among the modes for a specific wavelength (λ). We see that no high AO frequency is required to transfer efficiently energy from the core to cladding modes when λ shifts to infrared. Additionally, the smaller the tapered fiber diameter is, the more the required AO frequency decreases, which is a good feature because it could be difficult to find AO drivers at high frequencies.Figure 2
AO frequencyf
a as a function of wavelength at δ
(
λ
)
=
β
0
ν
(
λ
)
=
0 and coupling between the fundamental core mode and the L
P
1
ν
c
ls modes for (a) d
c
o
=
2.5
μ
m and (b) d
c
o
=
3.5
μ
m.
(a)
d
c
o
=
2.5
μ
m
(b)
d
c
o
=
3.5
μ
mFixing AO frequency atf
a
=
1.0
M
H
z, we compute the detuning δ
(
λ
) (5). The results can be seen in Figure 3. As the tapered diameter decreases the detuning δ is blue-shifted. δ-zero crossing moves toward blue side around 250
n
m per 1
μ
m. At fixed wavelength, δ increases for higher cladding order modes. From Figures 2 and 3, one can conclude that, by changing the tapered fiber diameter, AO frequency and detuning can be easily controlled, which makes this type of fiber a good candidate to design AO fiber filter in a large range of wavelengths.Figure 3
Detuning parameterδ as a function of wavelength for coupling between the fundamental core mode and the L
P
1
ν
c
ls modes when (a) d
c
o
=
2.5
μ
m and (b) d
c
o
=
3.5
μ
m.
(a)
d
c
o
=
2.5
μ
m
(b)
d
c
o
=
3.5
μ
mUsing (5) and (6), the resulting normalized transmitted power equation (7) in the core mode at the end of the fiber is shown in Figure 4 [18–20], where P
c
o
(
L
)
/
P
c
o
(
0
) in dB units is plotted. It is set: fiber length of L
=
15
c
m, AO wave amplitude of A
=
63
n
m, and f
a
=
1
M
H
z. According to (18), Λ
=
587.20
μ
m for d
c
o
=
2.5
μ
m and Λ
=
694.79
μ
m for d
c
o
=
3.5
μ
m. The coupling between L
P
01
c
o and L
P
12
c
l exhibits the biggest deep which means that the incident light in the core is partially transferred to the cladding mode at λ
=
680
n
m for d
c
o
=
2.5
μ
m. Similar behavior happens for d
c
o
=
3.5
μ
m at λ
=
923
n
m.Figure 4
Normalized transmitted power in the core mode at the end of the fiber. Setting: fiber length ofL
=
15
c
m, AO wave amplitude of A
=
63
n
m, f
a
=
1
M
H
z, d
c
o
=
2.5
μ
m, and d
c
o
=
3.5
μ
m.


## 5.2. Nonlinear Simulations

For nonlinear pulse simulations, (15) are numerically solved by using the split-step Fourier method [16]. Fiber length of 15
c
m and initial pulse width of T
o
=
50
f
s are considered. We analyze two scenarios: (i) d
c
o
=
2.5
μ
m: pump wavelength λ
o
=
1323
n
m, initial peak power P
o
=
30
k
W, and (ii) d
c
o
=
3.5
μ
m: pump wavelength λ
o
=
1472
n
m, initial peak power P
o
=
3
k
W. We choose the above pump wavelengths because we want to generate solitons at both diameters which allow us to get more insight into nonlinear effects.The linear coupling parameterC
12
(
λ
) is shown in graph in Figure 5 for the two diameters under consideration. C
12 almost scales linearly with wavelength, except for L
P
01
c
o-L
P
11
c
l at d
c
o
=
2.5
μ
m which presents some oscillation. The coupling between L
P
01
c
o and L
P
13
c
l exhibits the weakest strength at both cores diameters since the mode profiles do not overlap too much.Figure 5
AO coupling coefficientC
12 as a function of wavelength for the coupling between L
P
01
c
o and L
P
11
c
l, L
P
12
c
l, and L
P
13
c
l in a tapered standard fiber. Two core diameters are considered, d
c
o
=
2.5
μ
m and d
c
o
=
3.5
μ
m.The plots in Figures6, 7, and 8 present the overlap nonlinear field coefficients as a function of wavelength between the fundamental core mode and the L
P’s cladding modes. In Figure 6, we can see that this parameter that governs the SPM phenomenon in the core region is superior for tapered fiber with greater cross-sectional area for most of the analyzed range of wavelength. This could be due to the fact that more energy goes to the cladding as the core diameter is reduced. From Figure 7, d
22 almost linearly decreases as wavelength increases and has the same value for any L
P’s cladding mode in the same diameter. d
22 has small value since most of the light keeps in the single mode core, which makes SPM not efficiently excited at low power.Figure 6
Overlap field nonlinear coefficientd
11 as a function of wavelength in a tapered standard fiber. Two core diameters are considered, d
c
o
=
2.5
μ
m and d
c
o
=
3.5
μ
m.Figure 7
Overlap field nonlinear coefficientd
22 as a function of wavelength in a tapered standard fiber. Two core diameters are considered, (a) d
c
o
=
2.5
μ
m and (b) d
c
o
=
3.5
μ
m.
(a)
d
c
o
=
2.5
μ
m
(b)
d
c
o
=
3.5
μ
mFigure 8
Overlap field nonlinear coefficientd
12 as a function of wavelength in a tapered standard fiber. Two core diameters are considered, (a) d
c
o
=
2.5
μ
m and (b) d
c
o
=
3.5
μ
m.
(a)
d
c
o
=
2.5
μ
m
(b)
d
c
o
=
3.5
μ
mThe behavior of the overlap coefficientsd
12 is shown in Figure 8; they also decrease for higher wavelengths. The coefficient L
P
01
c
o-L
P
13
c
l is the strongest one, which means that more energy can be nonlinearly transferred from the core to the cladding weather; L
P
13
c
l is efficiently coupled to L
P
01
c
o through the acoustic wave.To complete the needed parameters to simulate (15), the dispersion coefficients β’s are calculated. These can be seen in Table 1.Table 1
Dispersion coefficientsL
P
01
c
o and LP cladding modes.
(a)
d
c
o
=
2.5 μm and λ
0
=
1323 nm
Mode
β
2 [ps2m−1]
β
3 [ps3m−1]
β
4 [ps4m−1]
β
5 [ps5m−1]
L
P
01
c
o
-
1.81
×
1
0
-
4
1.05
×
1
0
-
4
-
1.66
×
1
0
-
7
1.58
×
1
0
-
10
L
P
11
c
l
-
4.67
×
1
0
-
3
7.39
×
1
0
-
5
-
9.56
×
1
0
-
8
9.69
×
1
0
-
11
L
P
12
c
l
-
5.66
×
1
0
-
3
7.57
×
1
0
-
5
-
9.85
×
1
0
-
8
9.98
×
1
0
-
11
L
P
13
c
l
-
6.94
×
1
0
-
3
7.78
×
1
0
-
5
-
1.01
×
1
0
-
7
1.02
×
1
0
-
10
(b)
d
c
o
=
3.5 μm and λ
0
=
1472 nm
Mode
β
2 [ps2m−1]
β
3 [ps3m−1]
β
4 [ps4m−1]
β
5 [ps5m−1]
L
P
01
c
o
-
3.01
×
1
0
-
4
1.01
×
1
0
-
4
-
1.89
×
1
0
-
7
2.06
×
1
0
-
10
L
P
11
c
l
-
1.86
×
1
0
-
2
9.94
×
1
0
-
5
-
1.39
×
1
0
-
7
1.33
×
1
0
-
10
L
P
12
c
l
-
2.00
×
1
0
-
2
1.02
×
1
0
-
4
-
1.43
×
1
0
-
7
1.38
×
1
0
-
10
L
P
13
c
l
-
2.18
×
1
0
-
2
1.06
×
1
0
-
4
-
1.50
×
1
0
-
7
1.44
×
1
0
-
10Using (15), Figures 9 and 10 display the spectral response of all modes under consideration for both tapered fiber diameters. They evidence the generation of Supercontinuum (SC) both at the core and at the cladding. One can see that SC broadens more at the cladding modes when d
c
o
=
3.5
μ
m, even with less power compared to the d
c
o
=
2.5
μ
m case. This makes senses because d
12 is larger for d
c
o
=
3.5
μ
m (see Figure 8).Figure 9
Nonlinear spectral response ford
c
o
=
2.5
μ
m, with pump power of P
0
=
30
k
W, initial width of T
0
=
50
f
s, and pump wavelength of λ
0
=
1323
n
m.Figure 10
Nonlinear spectral response ford
c
o
=
3.5
μ
m, with pump power of P
0
=
3
k
W, initial width of T
0
=
50
f
s, and pump wavelength of λ
0
=
1472
n
m.As last step, we analyze how much energy can be transferred from the core to the LP cladding modes as a function of acoustic frequencyf
a for different pulse width T
0
=
50
f
s, 100
f
s, and 150
f
s; the results are plotted in Figure 11. The energy is normalized according to E
=
P
T
/
P
0
T
0, where P
T is computed as ∫
-
∞
∞
‍
P
(
T
)
d
T. From Figure 11, the transferred energy behaves linearly increasing as f
a and T
0 increase. This shows that although C
12 is proportional to f
a
2, the nonlinear effects linearize the energy and the slope can be controlled by T
0.Figure 11
Normalized transferred energy from the core to the LP cladding modes as a function of acoustic frequencyf
a for different pulse width T
0
=
50
f
s, 100 fs, and 150 fs. (a) For core diameter 2.5 μm and (b) core diameter 3.5 μm.
(a)
d
c
o
=
2.5
μ
m
(b)
d
c
o
=
3.5
μ
m


## 6. Conclusions

We have shown a complete procedure to study nonlinear pulse propagation in fiber optics when light is coupled from core mode to LP cladding modes as a consequence of acoustic wave. Based on local bending model, we introduce the concept of total variation of the refractive index as a sum of linear and nonlinear parts (Eqs. (13) and (14)). The procedure is applied to a tapered standard single mode fiber. We show that more energy can be nonlinearly transferred to the cladding if the tapered fiber diameter decreases. In addition, we have proven that the nonlinear effects compensate the coupling coefficients (C
12) which is proportional to f
a
2; this makes the transferred energy to cladding be linearized with respect to acoustic frequency.



---

*Source: 3696094-2016-09-18.xml*