# The coupled nonlinear Schrodinger equations describing power and phase for modeling phase-sensitive parametric amplification in silicon waveguides.

1. IntroductionNowadays, 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 [[omega].sub.p] passing their energy to a signal wave at frequency [[omega].sub.s] and an idler wave at frequency [[omega].sub.i] as the relation 2[[omega].sub.p] = [[omega].sub.s] + [[omega].sub.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 [DELTA]k = [DELTA][beta] + 2[[gamma].sub.p] [P.sub.pump] = 0. Here [DELTA][beta] = [k.sub.s] + [k.sub.i] - 2[k.sub.p] is the phase mismatch due to the linear dispersion, [P.sub.pump] is the pump power, [[gamma].sub.p] = [[omega].sub.p] [n.sub.2]/c[A.sub.eff] is the nonlinear waveguide parameter, [n.sub.2] = 12[[pi].sup.2] [[chi].sup.(3)]/[n.sub.0]c is the nonlinear index coefficient, c is the speed of light, [n.sub.0] is the linear refractive index, and [A.sub.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

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)

where [A.sub.j] is the slowly varying amplitude (j = p, s, i), z is the propagation distance, and [[beta].sub.2j] is the GVD coefficient. Time T = t - z/[v.sub.gp] is measured in a reference frame moving with pump pulse traveling at speed [v.sub.gp]. The two walk-off parameters of the signal and idler are defined as [d.sub.s] = [[beta].sub.1s] - [[beta].sub.1p] and [d.sub.i] = [[beta].sub.1i] - [[beta].sub.1p], respectively, where [[beta].sub.1j] is the inverse of the group velocity. The nonlinear coefficient [[gamma].sub.j] = [[omega].sub.j] [n.sub.2]/c[A.sub.eff] with [n.sub.2] = 6 x [10.sup.-18] [m.sup.2] [W.sup.-1] and [[beta].sub.TPA] = 5 x [10.sup.-12] m [W.sup.-1] is the coefficient of TPA at the wavelength of 1550 nm [2].

In (1), [[alpha].sub.j] accounts for the linear loss and [[alpha].sub.fj] = [[sigma].sub.j] [N.sub.c] represents FCA, where [[sigma].sub.j] is the FCA coefficient and [N.sub.c] is the free-carrier density generated by pump, signal, and idler pulses. [[delta].sub.nfj] = [[zeta].sub.j] [N.sub.c] is the free-carrier induced index change. These free-carrier parameters are obtained by solving [3, 19]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (2)

where [[lambda].sub.j] is the wavelength, [[lambda].sub.ref] = 1550 nm, h is Planck's constant, and [[tau].sub.c] [approximately equal to] 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.sub.j] = [absolute value of [A.sub.j]] exp (i[[phi].sub.j]), [P.sub.j] = [[absolute value of [A.sub.j]].sup.2]; (1) can be rewritten in terms of optical powers and phase:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (3)

where [theta] = [DELTA][beta]z + [[phi].sub.s] + [[phi].sub.i] - 2[[phi].sub.p] is the relative phase and [[phi].sub.p], [[phi].sub.s], and [[phi].sub.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 length h, 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:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (4)

The above ordinary differential equations can be solved by Runge-Kutta method, and the results can be expressed as [P.sub.p] (z + h)', [P.sub.s] (z + h)', [P.sub.i] (z + h)', [[phi].sub.p] (z + h)', [[phi].sub.s] (z + h)', and [[phi].sub.i] (z + h)', respectively. In the second step, dispersion acts alone mathematically,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (5)

where [F.sub.T] represents the Fourier-transform operation and w 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 Figure 1. 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.sub.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 [[theta].sup.Out-PIA] = [DELTA][[beta].sub.1] [L.sub.1] + [[phi].sub.s] ([L.sub.1]) + [[phi].sub.i] ([L.sub.1]), where [DELTA][[beta].sub.l] is the linear phase mismatch of the silicon waveguide with width of [W.sub.1] and [L.sub.1] is the length of the first segment. The second segment is channel waveguide with width of [W.sub.2], which has a normal dispersion at the pump wavelength. The relative phase at the input of the PSA is [[theta].sup.In-PSA] = [DELTA][[beta].sub.1] [L.sub.1] + [DELTA][[beta].sub.2] [L.sub.2] + [[phi].sub.s] ([L.sub.1] + [L.sub.2]) + [[phi].sub.l] ([L.sub.1] + [L.sub.2]) - 2[[phi].sub.p] ([L.sub.1] + [L.sub.2]), where [DELTA][[beta].sub.2] is the linear phase mismatch of the silicon waveguide with width of [W.sub.2] and [L.sub.2] is the length of the second segment. By changing the dispersion and length of the second segment, the relative phase [[theta].sup.In-PSA] can be set to an arbitrary value. The third silicon waveguide with width of W3 acts as a PSA, which can amplify or deamplify the signal depending on the relative phase at the input of PSA [[theta].sup.In-PSA].

The widths [W.sub.1], [W.sub.2], and [W.sub.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.2dB/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[[omega].sub.p] = [[omega].sub.s] [[omega].sub.i], the wavelength of the idler is 1801.7 nm. The dispersion parameters of our model are listed in Table 1.

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 1mW From Figure 2(a), it is found that the relative phase [theta] quickly increases to [pi]/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 [DELTA][[beta].sub.l] = -1016 [m.sup.-1]. It is clear that the signal peak power increases with the increase of the propagation length. That is because [theta] > 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 [theta] < 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 0Out-PIA [[theta].sup.Out-PIA] is 0.127[pi] 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 [DELTA][[beta].sub.2] can compensate the negative linear phase mismatch of the first segment [DELTA][[beta].sub.l]. The relative phase is tuned to 0.63[pi] 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 [DELTA][[beta].sub.3] = -389.5 [m.sup.-1]. Therefore, the second segment of the width-modulated waveguide can tune the relative phase to an appropriate value to realize an effective PSA.

With the increase of the propagation length of the second segment, the relative phase increases from 0.127[pi] to 2[pi] 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 [[theta].sup.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 [[theta].sup.In-PSA] is tuned to be 0.63[pi], while the minimal gain is 1.3 dB for [[theta].sup.In-PSA] = 1.64[pi]. From Figure 3, it is clear that the [[theta].sup.In-PSA] should be tuned between 0.127[pi] and 0.98[pi] 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.

The linear phase mismatch [DELTA][beta] 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.sub.2] of the second segment when the length of the second segment is 2 mm. The relative phase [[theta].sup.In-PSA] can be tuned from -0.11[pi] to 0.95[pi] by tailoring [W.sub.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[pi] < [[theta].sup.In-PSA] < 0.7[pi], which means that [W.sub.2] should be tailored to satisfy the condition: 610 nm < [W.sub.2] < 660 nm.

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.

http://dx.doi.org/10.1155/2014/621751

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This work was supported by the National Natural Science Foundation of China under Grant no. 61275134.

References

[1] J. Leuthold, C. Koos, and W. Freude, "Nonlinear silicon photonics," Nature Photonics, vol. 4, no. 8, pp. 535-544, 2010.

[2] L. Yin and G. P. Agrawal, "Impact of two-photon absorption on self-phase modulation in silicon waveguides," Optics Letters, vol. 32, no. 14, pp. 2031-2033, 2007.

[3] Q. Lin, J. Zhan, P. M. Fauchet, and G. P. Agrawal, "Ultrabroadband parametric generation and wavelength conversion in silicon waveguides," Optics Express, vol. 14, no. 11, pp. 4786-4799, 2006.

[4] Z. Wang, H. Liu, N. Huang, Q. Sun, and J. Wen, "Impact of dispersion profiles of silicon waveguides on optical parametric amplification in the femtosecond regime," Optics Express, vol. 19, no. 24, pp. 24730-24737, 2011.

[5] X. Li, Z. Wang, and H. Liu, "Optimizing initial chirp for efficient femtosecond wavelength conversion in silicon waveguide by split-step Fourier method," Applied Mathematics and Computation, vol. 218, no. 24, pp. 11970-11975, 2012.

[6] Q. Zhang and M. I. Hayee, "Symmetrized split-step fourier scheme to control global simulation accuracy in fiber-optic communication systems," Journal of Lightwave Technology, vol. 26, no. 2, pp. 302-316, 2008.

[7] H. Wang, "Numerical studies on the split-step finite difference method for nonlinear Schrodinger equations," Applied Mathematics and Computation, vol. 170, no. 1, pp. 17-35, 2005.

[8] X. Xu and T. Taha, "Parallel split-step Fourier methods for nonlinear Schrodinger-type equations," Journal of Mathematical Modelling and Algorithms, vol. 2, no. 3, pp. 185-201, 2003.

[9] G. M. Muslu and H. A. Erbay, "Higher-order split-step Fourier schemes for the generalized nonlinear Schrodinger equation," Mathematics and Computers in Simulation, vol. 67, no. 6, pp. 581-595, 2005.

[10] Z. Tong, C. Lundstrm, P. A. Andrekson, M. Karlsson, and A. Bogris, "Ultralow noise, broadband phase-sensitive optical amplifiers, and their applications," IEEE Journal on Selected Topics in Quantum Electronics, vol. 18, no. 2, pp. 1016-1032, 2012.

[11] J. Hansryd, P. A. Andrekson, M. Westlund, J. Li, and P. Hedekvist, "Fiber-based optical parametric amplifiers and their applications," IEEE Journal on Selected Topics in Quantum Electronics, vol. 8, no. 3, pp. 506-520, 2002.

[12] R. Hang, J. Lasri, P. S. Devgan, V. Grigoryan, P. Kumar, and M. Vasilyev, "Gain characteristics of a frequency nondegenerate phase-sensitive fiber-optic parametric amplifier with phase self-stabilized input," Optics Express, vol. 13, no. 26, pp. 10483-10493, 2005.

[13] S. Zhang, Z. Deng, and W. Li, "A precise Runge-Kutta integration and its application for solving nonlinear dynamical systems," Applied Mathematics and Computation, vol. 184, no. 2, pp. 496-502, 2007.

[14] M. Z. Liu, S. F. Ma, and Z. W. Yang, "Stability analysis of Runge-Kutta methods for unbounded retarded differential equations with piecewise continuous arguments," Applied Mathematics and Computation, vol. 191, no. 1, pp. 57-66, 2007.

[15] B. S. Attili, K. Furati, and M. I. Syam, "An efficient implicit Runge-Kutta method for second order systems," Applied Mathematics and Computation, vol. 178, no. 2, pp. 229-238, 2006.

[16] Z. Wang, H. Liu, N. Huang, Q. Sun, and J. Wen, "Influence of spectral broadening on femtosecond wavelength conversion based on four-wave mixing in silicon waveguides," Applied Optics, vol. 50, no. 28, pp. 5430-5436, 2011.

[17] R. L. Espinola, J. I. Dadap, R. M. Osgood Jr., S. J. McNab, and Y. A. Vlasov, "C-band wavelength conversion in silicon photonic wire waveguides," Optics Express, vol. 13, no. 11, pp. 4341-4349, 2005.

[18] G. P. Agrawal, Nonlinear Fiber Optics, 4th edition, 2007.

[19] L. Yin, Q. Lin, and G. P. Agrawal, "Soliton fission and supercontinuum generation in silicon waveguides," Optics Letters, vol. 32, no. 4, pp. 391-393, 2007.

[20] B. Jin, J. Yuan, C. Yu et al., "Efficient and broadband parametric wavelength conversion in a vertically etched silicon grating without dispersion engineering," Optics Express, vol. 22, pp. 6257-6268, 2014.

[21] H. Rong, Y. Kuo, A. Liu, M. Paniccia, and O. Cohen, "High efficiency wavelength conversion of 10 Gb/s data in silicon waveguides," Optics Express, vol. 14, no. 3, pp. 1182-1188, 2006.

Xuefeng Li, (1) Zhaolu Wang, (2) and Hongjun Liu (2)

(1) School of Science, Xi'an University of Post & Telecommunications, Xi'an 710121, China

(2) Xi'an Institute of Optics and Precision Mechanics, Chinese Academy of Sciences, Xi'an 710119, China

Correspondence should be addressed to Xuefeng Li; lixfpost@163.com

Received 16 June 2014; Revised 13 August 2014; Accepted 14 August 2014; Published 24 August 2014

Academic Editor: Ray K.L. Su

TABLE 1 Width [A.sub.eff] [alpha][beta] [[beta].sub.1p] ([mu]/[m.sup.2]) (1/m) (ps/m) [W.sub.1] 0.25 -1016 12912 [W.sub.2] 0.31 384.1 12791 [W.sub.3] 0.27 -389.5 12843 Width [[beta].sub.1s] [[beta].sub.1f] [[beta].sub.2p (ps/m) (ps/m) ([ps.sup.2]/m) [W.sub.1] 12951 12972 -0.04603 [W.sub.2] 12842 12850 0.004552 [W.sub.3] 12891 12901 -0.02276 Width [[beta].sub.2s [[beta].sub.2i ([ps.sup.2]/m) ([ps.sup.2]/m) [W.sub.1] 0.5444 -0.7975 [W.sub.2] 0.6014 -0.782 [W.sub.3] 0.5749 -0.8014

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Li, Xuefeng; Wang, Zhaolu; Liu, Hongjun |

Publication: | Journal of Applied Mathematics |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 3367 |

Previous Article: | On the domination number of Cartesian product of two directed cycles. |

Next Article: | An environmental-economic dispatch method for smart microgrids using VSS_QGA. |

Topics: |