# The Coupled Nonlinear Schrödinger Equations Describing Power and Phase for Modeling Phase-Sensitive Parametric Amplification in Silicon Waveguides

Authors: Xuefeng Li; Zhaolu Wang; Hongjun Liu
Journal: Journal of Applied Mathematics
 (2014)

---

## Abstract

The coupled nonlinear Schrödinger (NLS) equations describing power and phase of the optical waves are used to model phase-sensitive (PS) parametric amplification in a width-modulated silicon-on-insulator (SOI) channel waveguide. Through solving the coupled NLS equations by the split-step Fourier and Runge-Kutta integration methods, the numerical results show that the coupled NLS equations can perfectly describe and character the PS amplification process in silicon waveguides.


---

## Body

## 1. Introduction

Nowadays, silicon has emerged as a highly attractive material for nonlinear photonic integration [1]. Compared with highly nonlinear fiber, the SOI platform has inherent advantages due to the large values of Kerr parameter and Raman gain coefficient, the tight confinement of the optical mode, and the mature and low-cost fabrication process [2]. Optical parametric amplifications based on four-wave-mixing (FWM) in SOI waveguide have been theoretically investigated with the model of the coupled NLS equations describing the slowly varying amplitude of the optical waves [3–5]. The most commonly used numerical scheme for solving the NLSE is the split-step Fourier (SSF) method due to its simplicity for implementation and high computational efficiency [6–9]. Since phase-sensitive amplifiers (PSA) have the potential applications in optical communication, optical processing, photon detection, and optical spectroscopy and sensing [10], it is significant and crucial to model and investigate phase-sensitive amplification in SOI waveguide by using the coupled NLS equations describing power and phase of the optical waves for nonlinear photonic integration. The coupled NLS equations describing power and phase have been used to analyze the parametric process in fibers [11, 12]. Particularly, Hansryd et al. used the coupled NLS equations to analyze fiber-based optical parametric amplifiers [11]. Compared with optical fibers, the silicon waveguide has some additional complications, such as two-photon absorption (TPA), free-carrier absorption (FCA), and free-carrier dispersion (FCD). Therefore, we should consider TPA, FCA, and FCD in silicon waveguide to develop the coupled NLS equations.In this paper, through analyzing the coupled NLS equations describing slowly varying amplitude of the optical waves in FWM process, we develop the model of PS amplification based on coupled NLS equations describing power and phase of the optical waves in a width-modulated SOI channel waveguide. The model describes the power and phase variation in the silicon waveguide, which can be solved by the split-step Fourier and Runge-Kutta integration methods [5, 7–9, 13–15]. The calculation process using the solving methods has been discussed. The numerical results show that the coupled NLS equations can perfectly describe and character the PS amplification process in silicon waveguides.


## 2. Theory

The PS parametric amplification can be realized using FWM effect. Here, we focus on the degenerate FWM, which typically involves two pump photons at frequencyω
p passing their energy to a signal wave at frequency ω
s and an idler wave at frequency ω
i as the relation 2
ω
p
=
ω
s
+
ω
i holds. The signal wave is amplified and the idler wave is generated during the FWM process. Moreover, the phase-matching among the interacting waves is required in the FWM process, which is achieved when the mismatch in the propagation constants of the pump, signal, and idler waves is compensated by the phase shift due to SPM and XPM, such that Δ
k
=
Δ
β
+
2
γ
p
P
pump
=
0. Here Δ
β
=
k
s
+
k
i
-
2
k
p is the phase mismatch due to the linear dispersion, P
pump is the pump power, γ
p
=
ω
p
n
2
/
c
A
eff is the nonlinear waveguide parameter, n
2
=
12
π
2
χ
3
/
n
0
c is the nonlinear index coefficient, c is the speed of light, n
0 is the linear refractive index, and A
eff is the effective area of the propagating mode, respectively.The pump and signal waves are identically polarized in the fundamental quasi-TM mode. To depict the nonlinear optical interaction of the pump, signal, and idler in the waveguide, we use the formulism described in [16–18] and take into account the effects of TPA, FCA, and FCD. The coupled NLS equations describing slowly varying amplitude of the different optical waves read as
(1)
∂
A
p
∂
z
+
i
β
2
p
2
∂
2
A
p
∂
T
2
=
-
1
2
α
p
+
α
f
p
A
p
+
i
γ
p
+
i
β
TPA
2
A
eff
A
p
2
A
p
+
2
i
γ
p
+
i
β
TPA
2
A
eff
A
s
2
+
A
i
2
A
p
+
i
2
π
λ
p
δ
n
f
p
A
p
+
2
i
γ
p
A
s
A
i
A
p
*
exp
⁡
i
Δ
β
z
,
∂
A
s
∂
z
+
d
s
∂
A
s
∂
T
+
i
β
2
s
2
∂
2
A
s
∂
T
2
=
-
1
2
α
s
+
α
f
s
A
s
+
i
γ
s
+
i
β
TPA
2
A
eff
A
s
2
A
s
+
2
i
γ
s
+
i
β
TPA
2
A
eff
A
p
2
+
A
i
2
A
s
+
i
2
π
λ
s
δ
n
f
s
A
s
+
i
γ
s
A
p
2
A
i
*
exp
⁡
-
i
Δ
β
z
,
∂
A
i
∂
z
+
d
i
∂
A
i
∂
T
+
i
β
2
i
2
∂
2
A
i
∂
T
2
=
-
1
2
α
i
+
α
f
i
A
i
+
i
γ
i
+
i
β
TPA
2
A
eff
A
i
2
A
i
+
2
i
γ
i
+
i
β
TPA
2
A
eff
A
p
2
+
A
s
2
A
i
+
i
2
π
λ
i
δ
n
f
i
A
i
+
i
γ
i
A
p
2
A
s
*
exp
⁡
-
i
Δ
β
z
,
where A
j is the slowly varying amplitude (j
=
p
,
s
,
i), z is the propagation distance, and β
2
j is the GVD coefficient. Time T
=
t
-
z
/
v
g
p is measured in a reference frame moving with pump pulse traveling at speed v
g
p. The two walk-off parameters of the signal and idler are defined as d
s
=
β
1
s
-
β
1
p and d
i
=
β
1
i
-
β
1
p, respectively, where β
1
j is the inverse of the group velocity. The nonlinear coefficient γ
j
=
ω
j
n
2
/
c
A
eff with n
2
=
6
×
10
-
18 m2W−1 and β
TPA
=
5
×
10
-
12 mW−1 is the coefficient of TPA at the wavelength of 1550 nm [2].In (1), α
j accounts for the linear loss and α
f
j
=
σ
j
N
c represents FCA, where σ
j is the FCA coefficient and N
c is the free-carrier density generated by pump, signal, and idler pulses. δ
n
f
j
=
ζ
j
N
c is the free-carrier induced index change. These free-carrier parameters are obtained by solving [3, 19]
(2)
σ
j
=
1.45
×
1
0
-
21
λ
j
λ
ref
2
m
2
,
ζ
j
=
-
1.35
×
1
0
-
27
λ
j
λ
ref
2
m
3
,
∂
N
c
z
,
t
∂
t
=
π
β
TPA
2
h
ω
p
A
eff
2
A
p
z
,
t
4
-
N
c
z
,
t
τ
c
,
where λ
j is the wavelength, λ
ref
=
1550 nm, h is Planck’s constant, and τ
c
≈
1 ns is the carrier lifetime. In order to describe the power and phase of the different waves, we should find the coupled NLS equations describing power and phase. Let A
j
=
A
j
exp
⁡
i
φ
j, P
j
=
A
j
2; (1) can be rewritten in terms of optical powers and phase:
(3)
∂
P
p
∂
z
=
β
2
p
∂
P
p
∂
T
∂
φ
p
∂
T
+
β
2
p
P
p
∂
2
φ
p
∂
T
2
-
α
p
+
α
f
p
P
p
W
-
β
TPA
A
eff
P
p
2
-
2
β
TPA
A
eff
P
s
+
P
i
P
p
W
-
4
γ
p
P
s
P
i
P
p
2
1
/
2
sin
θ
,
∂
P
s
∂
z
=
-
d
s
∂
P
s
∂
T
+
β
2
s
∂
P
s
∂
T
∂
φ
s
∂
T
+
β
2
s
P
s
∂
2
φ
s
∂
T
2
-
α
s
+
α
f
s
P
s
W
-
β
TPA
A
eff
P
s
2
-
2
β
TPA
A
eff
P
p
+
P
i
P
s
W
+
2
γ
s
P
s
P
i
P
p
2
1
/
2
sin
θ
,
∂
P
i
∂
z
=
-
d
i
∂
P
i
∂
T
+
β
2
i
∂
P
i
∂
T
∂
φ
i
∂
T
+
β
2
i
P
i
∂
2
φ
i
∂
T
2
-
α
i
+
α
f
i
P
i
W
-
β
TPA
A
eff
P
i
2
-
2
β
TPA
A
eff
P
p
+
P
s
P
i
W
+
2
γ
i
P
s
P
i
P
p
2
1
/
2
sin
θ
,
∂
φ
p
∂
z
=
-
β
2
p
2
P
p
∂
2
P
p
∂
T
2
+
β
2
p
2
∂
φ
p
∂
T
2
+
γ
p
P
p
+
2
π
λ
p
δ
n
f
p
W
+
2
γ
p
P
s
+
P
i
+
2
γ
p
P
s
P
i
1
/
2
cos
⁡
θ
,
∂
φ
s
∂
z
=
-
d
s
∂
φ
s
∂
T
-
β
2
s
2
P
s
∂
2
P
s
∂
T
2
+
β
2
s
2
∂
φ
s
∂
T
2
+
γ
s
P
s
W
+
2
π
λ
s
δ
n
f
s
+
2
γ
s
P
p
+
P
i
+
γ
s
P
p
2
P
i
P
s
1
/
2
cos
⁡
θ
,
∂
φ
i
∂
z
=
-
d
i
∂
φ
i
∂
T
-
β
2
i
2
P
i
∂
2
P
i
∂
T
2
+
β
2
i
2
∂
φ
i
∂
T
2
+
γ
i
P
i
+
2
π
λ
i
δ
n
f
i
W
+
2
γ
i
P
p
+
P
s
+
γ
i
P
p
2
P
s
P
i
1
/
2
cos
⁡
θ
,
where θ
=
Δ
β
z
+
φ
s
+
φ
i
-
2
φ
p is the relative phase and φ
p, φ
s, and φ
i are the phases of the pump, signal, and idler. Equations (3) can be solved by split-step Fourier and Runge-Kutta integration methods, which are discussed in the following part.First, we assume the propagation length is divided into a large number of segments. Each segment length is defined as step lengthh, which is small enough that the nonlinear effects and dispersion can be treated independently. Two steps should be carried out for the propagation length from z to z
+
h. In the first step, the nonlinearity acts alone, and (3) are simplified as ordinary differential equations:
(4)
d
P
p
d
z
=
-
α
p
+
α
f
p
P
p
-
β
TPA
A
eff
P
p
2
-
2
β
TPA
A
eff
P
s
+
P
i
P
p
-
4
γ
p
P
s
P
i
P
p
2
1
/
2
sin
θ
,
d
P
s
d
z
=
-
α
s
+
α
f
s
P
s
-
β
TPA
A
eff
P
s
2
-
2
β
TPA
A
eff
P
p
+
P
i
P
s
+
2
γ
s
P
s
P
i
P
p
2
1
/
2
sin
θ
,
d
P
i
d
z
=
-
α
i
+
α
f
i
P
i
-
β
TPA
A
eff
P
i
2
-
2
β
TPA
A
eff
P
p
+
P
s
P
i
+
2
γ
i
P
s
P
i
P
p
2
1
/
2
sin
θ
,
d
φ
p
d
z
=
γ
p
P
p
+
2
π
λ
p
δ
n
f
p
+
2
γ
p
P
s
+
P
i
+
2
γ
p
P
s
P
i
1
/
2
cos
⁡
θ
,
d
φ
s
d
z
=
γ
s
P
s
+
2
π
λ
s
δ
n
f
s
+
2
γ
s
P
p
+
P
i
+
γ
s
P
p
2
P
i
P
s
1
/
2
cos
⁡
θ
,
d
φ
i
d
z
=
γ
i
P
i
+
2
π
λ
i
δ
n
f
i
+
2
γ
i
P
p
+
P
s
+
γ
i
P
p
2
P
s
P
i
1
/
2
cos
⁡
θ
.
The above ordinary differential equations can be solved by Runge-Kutta method, and the results can be expressed as P
P
z
+
h
′, P
s
z
+
h
′, P
i
z
+
h
′, φ
p
z
+
h
′, φ
s
z
+
h
′, and φ
i
z
+
h
′, respectively. In the second step, dispersion acts alone mathematically,
(5)
P
p
z
+
h
=
F
T
-
1
F
T
P
p
z
+
h
′
exp
⁡
-
2
β
2
p
ω
2
φ
p
z
h
,
P
s
z
+
h
=
F
T
-
1
F
T
P
s
z
+
h
′
exp
⁡
⁡
i
ω
d
s
h
-
2
β
2
s
ω
2
φ
s
z
h
,
P
i
z
+
h
=
F
T
-
1
F
T
P
i
z
+
h
′
exp
⁡
⁡
i
ω
d
i
h
-
2
β
2
i
ω
2
φ
i
z
h
,
φ
p
z
+
h
=
F
T
-
1
F
T

φ
p
z
+
h
′
exp
⁡
⁡
ω
2
β
2
p
2
φ
p
z
h
-
ω
2
β
2
p
φ
p
z
2
h

,
φ
s
z
+
h
=
F
T
-
1
F
T
φ
s
z
+
h
′
i
ω
d
s
h
+
ω
2
β
2
s
2
φ
s
z
h
-
ω
2
β
2
s
φ
s
z
2
h
W
W
W
W
W
×
exp
⁡
⁡
i
ω
d
s
h
+
ω
2
β
2
s
2
φ
s
z
h
-
ω
2
β
2
s
φ
s
z
2
h
,
φ
i
z
+
h
=
F
T
-
1
F
T
φ
i
z
+
h
′
i
ω
d
i
h
+
ω
2
β
2
i
2
φ
i
z
h
-
ω
2
β
2
i
φ
i
z
2
h
W
W
W
W
W
×
exp
⁡
⁡
i
ω
d
i
h
+
ω
2
β
2
i
2
φ
i
z
h
-
ω
2
β
2
i
φ
i
z
2
h
,
where F
T represents the Fourier-transform operation and ω is the frequency in the Fourier domain [18]. Therefore, the coupled NLS equations describing powers and phase of different optical waves are solved by split-step Fourier and Runge-Kutta integration methods. The detailed numerical results in the third part of the paper will show that the split-step Fourier and Runge-Kutta integration methods provide accurate and stable solution for the coupled NLS equations describing power and phase.Here the coupled NLS equations describing power and phase are used to model and simulate the PS parametric amplification in a width-modulated SOI waveguide, which is comprised of three segments of channel waveguides with different widths and identical height as shown in Figure1. Tapers are used to connect the three segments to avoid the mode mismatch induced by the variation of width [20]. In the first segment, the SOI waveguide with width of W
1 has an anomalous dispersion at the pump wavelength of 1550 nm, which acts as a PIA to amplify the signal and generate an idler. The relative phase at the output of the PIA is θ
Out-PIA
=
Δ
β
1
L
1
+
φ
s
(
L
1
)
+
φ
i
(
L
1
)
-
2
φ
p
(
L
1
), where Δ
β
1 is the linear phase mismatch of the silicon waveguide with width of W
1 and L
1 is the length of the first segment. The second segment is channel waveguide with width of W
2, which has a normal dispersion at the pump wavelength. The relative phase at the input of the PSA is θ
In-PSA
=
Δ
β
1
L
1
+
Δ
β
2
L
2
+
φ
s
(
L
1
+
L
2
)
+
φ
i
(
L
1
+
L
2
)
-
2
φ
p
(
L
1
+
L
2
), where Δ
β
2 is the linear phase mismatch of the silicon waveguide with width of W
2 and L
2 is the length of the second segment. By changing the dispersion and length of the second segment, the relative phase θ
In-PSA can be set to an arbitrary value. The third silicon waveguide with width of W
3 acts as a PSA, which can amplify or deamplify the signal depending on the relative phase at the input of PSA θ
In-PSA.Figure 1
Illustration of the phase-sensitive parametric amplification in a width-modulated SOI channel waveguide with identical height.The widthsW
1, W
2, and W
3 of the width-modulated SOI waveguide are assumed as 500 nm, 650 nm, and 580 nm, respectively, while the height is 800 nm. The linear propagation losses of the three segments are set to be 0.3 dB/cm, 0.2 dB/cm, and 0.25 dB/cm, respectively [21]. The PSA process is theoretically investigated with pump pulse of 20 ps at the wavelength of 1550 nm and continuous-wave signal at the wavelength of 1360 nm. According to the relation 2
ω
p
=
ω
s
+
ω
i, the wavelength of the idler is 1801.7 nm. The dispersion parameters of our model are listed in Table 1.Table 1
Width
A
eff (μm2)
Δ
β (1/m)
β
1
p (ps/m)
β
1
s (ps/m)
β
1
i (ps/m)
β
2
p (ps2/m)
β
2
s (ps2/m)
β
2
i (ps2/m)
W
1
0.25
−1016
12912
12951
12972
−0.04603
0.5444
−0.7975
W
2
0.31
384.1
12791
12842
12850
0.004552
0.6014
−0.782
W
3
0.27
−389.5
12843
12891
12901
−0.02276
0.5749
−0.8014


## 3. Results and Discussion

The coupled NLS equations of (3) are solved by using split-step Fourier and Runge-Kutta integration methods to investigate the PS parametric amplification process in the width-modulated SOI channel waveguide as shown in Figure 2. The initial pump peak power is set to be 5 W, while the initial signal power is set to be 1 mW. From Figure 2(a), it is found that the relative phase θ quickly increases to π
/
2 when the initial phases of pump and idler are zero due to the generation of the idler and then decreases as the propagation length increases because of the large negative linear phase mismatch Δ
β
1
=
-
1016 m−1. It is clear that the signal peak power increases with the increase of the propagation length. That is because θ
>
0 along the propagation length, and the energy of the pump is transferred to signal and idler according to (3). It is expected that the signal peak power will decrease when further increasing the propagation length leads to θ
<
0, because the energy will be transferred back to pump. The output signal peak power is up to 5.07 mW, and relative phase at the output of the PIA θ
Out-PIA is 0.127
π when the length of the first segment is 8 mm. The second segment of the waveguide stops the decrease of the relative phase, because the positive linear phase mismatch of the second segment Δ
β
2 can compensate the negative linear phase mismatch of the first segment Δ
β
1. The relative phase is tuned to 0.63
π and the signal peak power is increased to 7.57 mW through the second segment of the width-modulated waveguide as shown in Figure 2(b). The phase-sensitive amplification occurs in the third segment as shown in Figure 2(c). It exhibits exponential gain, and the signal peak power is amplified to 56.8 mW when the propagation length of the third segment is 10 mm. The relative phase in the third segment decreases as the propagation length increases due to a negative linear phase mismatch Δ
β
3
=
-
389.5 m−1. Therefore, the second segment of the width-modulated waveguide can tune the relative phase to an appropriate value to realize an effective PSA.The signal peak power and relative phase as a function of the propagation length for the first segment (a), the second segment (b), and the third segment (c), respectively.
(a)
(b)
(c)With the increase of the propagation length of the second segment, the relative phase increases from0.127
π to 2
π as shown in Figure 3. However, the signal gain has a peak and a valley along the propagation length, which depends on the relative phase θ
In-PSA. Here, we define the signal gain as the ratio of the output signal power to the input signal power for the width-modulated SOI waveguide. The maximum signal gain of 17.5 dB is obtained when θ
In-PSA is tuned to be 0.63
π, while the minimal gain is 1.3 dB for θ
In-PSA
=
1.64
π. From Figure 3, it is clear that the θ
In-PSA should be tuned between 0.127
π and 0.98
π for the PSA to obtain a gain larger than 15 dB, which means the length of the second segment should be less than 4.6 mm. From Figures 2 and 3, it is clear that the coupled NLS equations describing power and phase of the optical waves solved by split-step Fourier and Runge-Kutta integration methods provide an accurate description of the PS amplification process.Figure 3
The relative phase at the input of PSAθ
In-PSA and signal gain as a function of the propagation length for the second segment.The linear phase mismatchΔ
β is simulated for different waveguide widths as shown in Figure 4(a), which can be used to simulate the phase-sensitive amplification process for different width of the second segment. The relative phase and signal gain are investigated by tailoring the width W
2 of the second segment when the length of the second segment is 2 mm. The relative phase θ
In-PSA can be tuned from -
0.11
π to 0.95
π by tailoring W
2 ranging from 500 nm to 700 nm as shown in Figure 4(b). It is found that the relative higher gain over 17 dB can be obtained for 0.41
π
<
θ
In-PSA
<
0.7
π, which means that W
2 should be tailored to satisfy the condition: 610
nm
<
W
2
<
660
nm.The linear phase mismatchΔ
β (a), the relative phase at the input of PSA θ
In-PSA and signal gain (b) as a function of the waveguide width for the second segment.
(a)
(b)


## 4. Conclusion

The phase-sensitive parametric amplification process in a width-modulated silicon waveguide is described by the model of coupled NLS equations describing power and phase of the optical waves, which can be solved by the split-step Fourier and Runge-Kutta integration methods. Numerical results show that the split-step Fourier and Runge-Kutta integration methods provide accurate and stable solution for coupled NLS equations, which can perfectly describe the PS amplification process in silicon waveguides, and the PS amplification can be achieved in our designed width-modulated silicon waveguide.



---

*Source: 621751-2014-08-24.xml*