Competition mechanism of multiple four-wave mixing in highly nonlinear fiber: spatial instability and satellite characteristics
Liang ZHAO
, 1
,
Junqiang SUN
2
,
Xinliang ZHANG
2
,
Cong CHEN
3
Expand
1. Wuhan Foreign Languages School, Wuhan 430022, China
2. Wuhan National Laboratory for Optoelectronics, School of Optoelectronic Science and Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
3. Naval University of Engineering, Wuhan 430033, China
Received date: 20 Jul 2012
Accepted date: 30 Aug 2012
Published date: 05 Dec 2012
Copyright
2014 Higher Education Press and Springer-Verlag Berlin Heidelberg
Competition mechanism in multiple four-wave mixing (MFWM) processes is demonstrated theoretically. Provided considering only two waves injected into a highly nonlinear fiber (HNLF), there are three modes displaying comprehensive dynamic behaviors, such as fixed points, periodic motion, and chaotic motion. Especially, Mode C of MFWM is emphasized by analyzing its phase-space trajectory to demonstrate nonlinear wave-wave interactions. The study shows that, when the phase-space trajectory approaches or gets through a saddle point, a dramatic power depletion for the injected wave can be realized, with the representative point moving chaotically, but when phase-space trajectories are distributed around a center point, the power for the injected wave is retained almost invariable, with the representative point moving periodically. Finally, the evolvement of satellite wave over an optical fiber is investigated by comparing it with the interference pattern in Young’s double-slit experiment.
Liang ZHAO, Junqiang SUN, Xinliang ZHANG, Cong CHEN. Competition mechanism of multiple four-wave mixing in highly nonlinear fiber: spatial instability and satellite characteristics[J]. Frontiers of Optoelectronics, 2012, 5(4): 414-428. DOI: 10.1007/s12200-012-0280-z
Introduction
Four-wave mixing (FWM), the most pervasive and persuasive modern technology in optical communication, marked by rapid change and growth, is moving into a new era, an era of extraordinary sophistication and versatility. It enjoys the obvious advantage of generating new frequencies, giving rise to a potential breakthrough in the applicative regime, such as broadcast all-optical wavelength conversion [1,2], fiber optical parametric amplification (FOPA) [3,4], all-optical switching [5,6], and 2R regeneration [7,8].
With the rapid development of multiple four-wave mixing (MFWM) theory, it has attracted increasing attention from experimentalists and theorists alike. Liu et al. have developed MFWM in theory and experiment when pump waves have equal frequency spacing, and have found the self-stability effect in optical fibers under the undepleted approximation [9,10]. Moreover, Cappellini and Trillo have investigated the degenerate FWM effect in a strong-interaction regime [11], with the aid of Hamilton system and phase-space trajectory. Similar investigations can be found in Refs. [12-14]. In fact, when taking account of two waves injected into a piece of highly nonlinear fiber (HNLF) and the occurrence of the high conversion efficiency, there are three sub-FWM processes, corresponding to three modes, existing in HNLF, as seen in Fig. 1. Moreover, the above three modes will interact and compete with each other, leading to energy switching from one wave to another in the largest extent.
An important phenomenon in MFWM is the power-depletion effect, which can find many applications in optical communications system, such as all-optical switching and notch filter. In Ref. [15], we theoretically demonstrated a novel notch filter using the dramatic power-depletion effect in HNLF, which can filter out the undesired wave through switching its energy to other waves. In experiment, Marhic et al. realized 92% pump depletion in a continuous-wave FOPA [16], with a 200-mW pump at 1560 nm in an 11-km dispersion-shifted fiber. In our paper, we put forward the competition-mechanism conception in MFWM processes, and we investigate Mode C of MFWM in detail by yielding the phase-space equation and the fixed points analytically from the coupled-wave equations to analyze the stability of the system. Furthermore, we discuss the stability of the conserved system (referred to the MFWM system) by deriving the equilibrium points from the potential well equation analytically. At last, the evolvement of satellite wave over an optical fiber is compared with Young’s double-slit experiment [17,18], and the analytical expressions, describing periodic and amplitude characteristics of satellite wave, are yielded, which can explain the physical implication of the dynamical process.
Competition mechanisms in MFWM processes
Figure 1 illustrates three modes in HNLF, which indicates that MFWM processes can be decomposed into three sub-FWM processes, namely Mode A, Mode B and Mode C, whose angular frequency relations satisfy 2ω1 = ω2 + ω3, 2ω2 = ω1 + ω4, ω1 + ω2 = ω3 + ω4, respectively. One can easily obtain the following conclusion based on Ref. [9]: when the center wavelength in terms of the injected wavelength is just located at around the zero-dispersion wavelength, Mode A will play a dominate role in MFWM processes when P1>>P2, and Mode B will occupy the absolute advantage when P2>>P1. Last but not least, Mode A and Mode B will be in dynamical balance when P1 = P2 (P1≈P2), and then Mode C will dominate in the whole process. When neglecting the walk-off, absorption, and dispersion effects in HNLF, the coupled-wave equations can be robustly expressed as
Here, Aj (j = 1, 2, 3, 4) is the complex amplitude for each wave, and Δβj is the linear phase mismatch profile, which can be analytically derived by expanding the propagation constant in Taylor series [19]. Qualitative and semiquantitative method here is adopted to analyze the energy switching among the waves. Under the circumstances of slowly-varied amplitude approximation and undepleted approximation, we can integrate the first FWM term of Eq. (1) with the result ofwhere represents the amount of power-varying for wave 1 induced by the first FWM term in Eq. (1), and the subscript “1” represents wave 1 while the superscript “1-FWM” denotes the first FWM term, Pj (j = 1, 2, 3, 4) denotes the power of each wave. Similarly, one can yield all the involved FWM terms for (i = 1, 2, 3, 4 and j = 1, 2, 3), which are listed in Table 1. Notably, in all the FWM terms, there is a common factor, namely sinc function, featuring oscillation and damping characteristics, which conforms to the experimental results demonstrated by Hart et al. [20].
Fig.1 Schematic diagram for the decomposition of MFWM
Tab.1 Contribution of all the FWM terms to the injected and sideband waves
1-FWM
2-FWM
3-FWM
1
2
3
4
Remarkably, Mode A and Mode B have been extensively investigated in the previous literatures [3,10,13] for the degenerate FWM case, but Mode C, corresponding to the non-degenerate FWM, is seldom mentioned before. Differing from the non-degenerate FWM process in a FOPA model, in which three different waves are prerequisite to excite the system, only two different waves are needed for Mode C in our system. Under the conditions of P1 ≈ P2 and Δβ1>>0, Δβ2>>0, Δβ3→0 [9], Eqs. (1)-(4) are modified to
Spatial instability
This section focuses on the characteristics of the solutions for Eqs. (6)-(9), and a deeper insight into the spatial instability for nonlinear eigenmodes of interaction is investigated extensively. In fact, as a consequence of the existence of spatially unstable eigensolutions, small variations in certain input conditions lead to quite different spatial evolutions and therefore to dramatically different output power redistributions among the waves. By inserting Aj = ξjexp(iϕj), (j = 1, 2, 3, 4), Г = γ·z, η = and κ3 = Δβ3/γ into Eqs. (6)-(9), one can easily obtainwith the Manley-Rowe relations [9]where the constant x1, x2 and x3 are determined by the initial wave amplitude ξj (j = 1, 2, 3, 4). The relative phase difference θ3 [9] is determined by
The system governed by Eqs. (6)-(9), besides the conserving total power (see Eq. (11)), admits of the following additional invariant, namely Hamiltonian,
Hence, Hamilton system equations explicitly read as
Equations (6)-(9) reduce to a simplified form, namely one-dimensional Hamilton system, as shown in Eqs. (14) and (15). The advantage of this simplification, concerning on the topology of the phase-space orbit, gives us an insight into the properties of energy exchange among waves, which controls the stability characteristic of a nonlinear system. Furthermore, the property of equilibrium points (also called singular points or fixed points) can be readily investigated by exploiting the well-developed theory for Hamilton system.
Phase-space analysis for equal-amplitude sidebands
In this case, the injected wave power is regulated by Mode C, and the condition of P1 = P2 should be strictly satisfied, which gives rise to the generation of equal-power sidebands. To simplify the calculation, supposing x1 = x2 = 1 and x3 = 0, Eqs. (14) and (15) are transformed into
To best unveil the stability characteristic of Hamilton system, one can transform Eqs. (16) and (17) into the form under the rectangular coordinate,where x = ηcosθ3, y = ηsinθ3. Moreover, the Jacobi matrix is available from Eqs. (18) and (19):
From Eqs. (18) and (19), one can also explicitly yield a valuable equation describing the phase-space characteristics,where C is a constant, dependent seriously on the boundary condition of the system, and Eq. (21) depicts the characteristic of phase-space portraits. The topology of phase-space orbits is significantly affected by the characteristic of singular point, which makes a significant influence to the stabilization of the system. From Eqs. (18) and (19), one can easily obtain singular points for the following four different types:
1) (x, y) = ((κ3 + 6)/12, 0). Remarkably, the liner-phase-mismatch profile should satisfy the condition of κ3>-6. Furthermore, one can get the eigenvalue equation: λ2 = (κ3 + 6)(κ3-6)/3. If κ3 lies in the interval of (-6, 6), the singular point will turn out to be a center point; otherwise, if κ3 is in the range of (6, +∞), the singular point will be a saddle point, and they are both controlled by the characteristic of eigenvalue [10].
2) (x, y) = ((κ3-2)/4, 0). Based on the similar principal, the eigenvalue equation can be yielded with λ2 = (κ3 + 2)(κ3-2), under the condition of κ3<2. Notably, the singular point will be a center point as long as -2<κ3<2; otherwise, the singular point will be turned into a saddle point if κ3<-2.
3) (x, y) = (0, 0). This specially fixed point resides in the liner-phase-mismatch condition of κ3∈(-∞, -6)∪(2, +∞), characterized by the eigenvalue equation λ2 = -(κ3 + 2)2, and it will be definitely the center point.
4) . Interesting is that the condition for the existence of this type singular point is -2<κ3<6, and the fixed point will be nothing but a saddle point. This conclusion is drawn by analyzing the characteristics of the eigenvalue equation of λ2 = - (κ3 + 2)(κ3 - 6).
Consider the anomalous-dispersion regime first. Figures 2(a) and 2(b) display phase-space portraits with κ3 = -7 and κ3 = -3, describing FWM dynamics, and one can clearly see that phase-space portraits contain only one center point, around which the representative point of the field moves on stable periodic trajectories in the physical limits regime, and thus the sideband wave varies in a stable mode. The fact that, κ3 with a relatively large value leads to a small periodic exchange of energy between the injected and sideband waves, is demonstrated in our case. Moreover, out of the physical limits regime, the unstable saddle point, which would be present when two eigensolutions exchange their stability, does not assume the physical significance.
Figure 2(c) shows a completely different scenario when κ3 is increased up to κ3 = -1. In fact, when κ3 = -2 (referred to Fig. 3(a)), the exchange of stability between eigensolutions is a standard characteristic exhibiting the transcritical bifurcation, in which the original stable mode vanishes, accompanied with the emergence of two newly central points (stable mode) and two newly saddle points (unstable mode). Moreover, if the representative point moves around the saddle point, the energy can be significantly coupled from the injected to sideband wave. As for two saddle points in Fig. 2(c), it is clearly found that they simultaneously have both the stable and unstable characteristics. On the one hand, the arrow on the trajectory points to the equilibrium point marked with yellow color, featuring the stability characteristic; on the other hand, the arrow, deviating from the above equilibrium point, points to another equilibrium point, exhibiting the unstable characteristic for the first equilibrium point. The trajectory that connects two saddle points is named “heteroclinic orbit” and the corresponding equilibrium point is called “heteroclinic point.” As seen in Fig. 2(c), the spatial instability of sidebands entails that power exchange from the injected to sideband waves is greatly enhanced. Remarkably, the spatial period of the conversion diverges to infinity on the loop heteroclinic separatrix that stems from the saddle points at ().
Going on increasing the value of κ3 up to the normal-dispersion regime, it is found that the topology of phase-space orbits retains invariable, until κ3 is increased up to κ3 = 6 (see Fig. 3(a)). As shown in Fig. 2(e), phase-space trajectories are also divided into four separated regions, which is similar with Fig. 2(c). Differently, the saddle points are displaced from the left to right half-plane. Interestingly, as alluded earlier, when κ3≥6, the topology of phase-space orbits is dramatically changed and comprises one stable center point (in the physical limits regime) and one unstable saddle point (out of the physical limits regime), and this characteristic is illustrated in Fig. 2(d) with κ3 = 10. In fact, transcritical bifurcation takes place under the condition of κ3 = 6.
Fig.2 Phase-space portraits with abnormal dispersion (see (a), (b) and (c)) and normal dispersion (see (d), (e) and (f)) regimes. (a) κ3 = -7; (b) κ3 = -3; (c) κ3 = -1; (d) κ3 = 10; (e) κ3 = 5; (f) κ3 = 1. Physical limits require that -1≤ηcosθ3≤1 and -1≤ηsinθ3≤1
As illustrated in Fig. 3(a), distribution of singular point I and II implies that cosθ3 is proportional to κ3, and the center and saddle points are marked with circle and asterisk, respectively. Remarkably, the critical point 1 (6, 1) and point 2 (-2, -1) simultaneously have both the stable and unstable characteristics, and the singular point IV, expressed as , varies with κ3 in a nonlinear manner, and it travels along with two different routes, marked with different colors in Fig. 3(b). Remarkably, the projection of 3D plot at x-y plane is just a perfect circle, which implies that the representative point moves on a circle orbit. Finally, the equilibrium point (0, 0) will survive under the condition of κ3∈(-∞, -6)∪(2, +∞), and one can safely obtain the conclusion that if the representative point moves on a circle orbit around the point (0, 0), energy exchange between the injected and sideband waves is small.
Fig.3 Illustration of distribution of singular points with different values of κ3. (a) Plot for singular point I (x, y) = ((κ3 + 6)/12, 0), singular point II (x, y) = ((κ3-2)/4, 0) and singular point III (x, y) = (0, 0); (b) plot for singular point IV
Phase-space analysis for approximate-equal-amplitude sidebands
An assumption of x1 = 1 and x2 ≈ 1 with x3 = α ≈ 0 is made in this section, referred to the case for approximate-equal-amplitude sidebands. Strictly speaking, Mode C survives under the boundary condition of x1 = x2 = 1 and x3 = 0, but for x2 ≈ 1, Mode C also plays a dominate role in MFWM processes. One can easily get the following equations from Eqs. (14) and (15) with
Equations (22) and (23) can be adopted to analyze a conserved Hamilton system describing Mode C. For the sake of simplicity, we can analyze the Hamilton system under the polar coordinate with the Jocabi matrix given by
The relation of x and y can be explicitly yielded from Eqs. (22) and (23):where C, x and y have been defined in Eq. (21). Remarkably, Eq. (25) can be degenerated into Eq. (21) if α = 0. Similarly, the stable points for Eqs. (22) and (23) can be classified as the following two aspects:
1) or . The fixed point, given by under the rectangle coordinate, will exist under the condition of . Furthermore, one can easily obtain the eigenvalue with , which implies that the fixed point will be definitely an unstable saddle point.
2) (η, θ3) = (η0, π) or (η0, 2π). η0 is the root of the following equation:where a = (9α + κ3-6)/6, b = [27α2 + 10(κ3-6)α + 12-4κ3-]/48, c = [-α3 + 2(κ3-10) α2 + (12-4κ3-)/48], d = α2/12. The eigenvalue equation for this case is
By substituting the root η0 of Eq. (26) into Eq. (27), the eigenvalue λ can be thus obtained, and we can further judge the type of a singular point that if 0<η0<1, it will be a center point; otherwise, it will turn to be a saddle point. Under the rectangle coordinate, the corresponsive fixed points are (x, y) = (±η0, 0), which implies that the fixed points appear at x-axis and are located symmetrically in terms of y-axis.
Figure 4 depicts distribution of singular points when taking account of α = 0.01. As mentioned before, firstly, the point satisfies the relation of x2 + y2 = 1, which represents that it moves on a circle trajectory, as shown in Fig. 4(a). It is worthy noting that when κ3 = 3.5, two saddle points P3 (0.3725, 0.928) and P4 (0.3725, -0.928) can be obtained, corresponding to yellow points in Fig. 5(b). Secondly, singular points (x, y) = (±η0, 0) are illustrated in Fig. 4(b), and interesting is that κ3 varies with η in a linear manner, but in different directions. Moreover, two critical points, provided with both the stable and unstable characteristics, are identified to (1, 6.03) and (1, -2.01).
Fig.4 Illustration of distribution of singular points with different values of κ3. (a) Plot for singular point ; (b) plot for singular point (x, y) = (±η0, 0)
Figure 5 illustrates phase-space portraits with κ3 = -3, κ3 = 3.5 and κ3 = 7, respectively. From what has been discussed above, we know that κ3 = -2.01 and κ3 = 6.03 are two critical linear-phase-mismatch values that dramatically change the topology structure of phase-space trajectories, and that further influence energy exchange among the waves. As shown in Figs. 5 (a) and 5(c), the saddle point shifts from the left to right half-plane for phase-space orbits, with respect to the center point, during which the number of saddle point is also increased up to two, as seen in Fig. 5(b).
Fig.5 Illustration for phase-space portraits with (a) κ3 = -3; (b) κ3 = 3.5 and (c) κ3 = 7. Physical limits require that -1≤ηcosθ3≤1 and -1≤ηsinθ3≤1
Figure 6(a) illustrates sideband power as a function of Γ when κ3 varies from -5 to 5. As long as κ3 increases approximately from -2.5 to 0.5, intense interaction between the injected and sideband wave takes place, and the sideband wave can obtain high energy from the injected wave. In addition, projection of Fig. 6(a) at η-κ3 plane is also displayed in the inset of Fig. 6(b), and the red curve in Fig. 6(b) indicates fractional power loss (FPL) for the injected wave as a function of κ3, which can be also regarded as the envelope of the inset. Here, FPL is defined as η/1 × 100%, reflecting power attenuation of the injected wave. It is expected that FPL will exceed 50% when κ3 ranges from -2.4 to 0.4, implying the occurance of a dramatic energy exchange between the injected and sideband waves. Furthermore, a FPL of about 100% can be attained if κ3 ≈ -0.4, and then almost all the injected wave power is delivered to the sideband wave.
Fig.6 (a) Sideband power as a function of Γ with different values of κ3; (b) fractional power loss as a function of κ3
To further unveil the relation between phase-space orbit and power exchange among the waves, we explicitly mimic and compare phase-space orbits, η vs. Г curve and cosθ3 vs. Г curve, under the condition of κ3 = -1 and κ3 = 4, corresponding to unstable and stable modes, as seen in Fig. 7. Especially, Fig. 7(a) illustrates unstable phase-space orbits, on which the representative point moves in a larger range than that shown in Fig. 7(d). η vs. Г curve in Fig. 7(b) indicates that the sideband wave can obtain approximate 80% injected-wave power, which is much larger than that shown in Fig. 7(e) and the sideband wave power in Fig. 7(b) varies with Г periodically with its period almost four times larger than that shown in Fig. 7(e). On another perspective, a large period for η vs. Г curve is benefit to energy accumulation when the sideband wave propagating along an optical fiber, and the amount of period is significantly dependent on cosθ3 vs. Г curve, which in facet reflects the varying rate of power for the sideband wave. The varying rate for cosθ3 vs. Г curve in Fig. 7(c) is much slower than that in Fig. 7(f), and this characteristic can be adopted to explain the phenomena in Figs. 7(b) and 7(e).
Fig.7 Illustration of phase-space trajectories representing (a) an unstable saddle point with κ3 = -1 and (d) a stable center point with κ3 = 4. η vs. Г with (b) κ3 = -1 and (e) κ3 = 4. cosθ3 vs. Г with (c) κ3 = -1 and (f) κ3 = 4
Through using simultaneous Eqs. (14) and (15), and taking account of the initial conditions of x1 = x2 = 1 and x3 = 0, one can obtainwith and . Equation (28) can be well explained as an evolution equation for the velocity of an ideal point moving in a potential well V(η) with a totally conserved energy E, which may be integrated to yield the solution in an implicit form,with the potential well explicitly reading asand the totally conserved energyis uniquely characterized by the initial settings for H and κ3. Γ can be explicitly yielded from Eq. (29) with the form of periodic Jacobi elliptic and hyperbolic function [10]. The form of the solution for Γ presumably depends on the root of the fourth-order polynomial f(η).
Potential well is a key parameter that can estimate the capability of power exchange among waves in a conserved Hamilton system, and it also influences the characteristic of a singular point and the topology of phase-space orbits. Figure 8(a) illustrates potential well energy V(η) vs. normalized power η with H = 1 and different values of κ3, and all the concave points, namely stable equilibrium points, are assembled to form curve 1 in Fig. 8(a). What’s more, stable eigensolutions are available from curve 1, which are increased with the decrement of κ3. Conversely, all the salient points in Fig. 8(a), assembled to form curve 2, corresponds to unstable equilibrium points, and they are decreased with the decrement of κ3. Figure 8(b) further illustrates potential well energy V(η) vs. normalized power η with κ3 = 1 and different values of H. Notably, different values of H denote different initial conditions for injected waves, and stable equilibrium points with η = 0.5 will be turned into unstable equilibrium points if H increases from<2.5 to>2.5 (see curve 2), where H = 2.5 is the critical point marked with five-pointed star in Fig. 8(b). Curve 1 in Fig. 8(b) is the aggregation of salient points, which does not assume physical significance. Simulation for potential well energy as a function of η and κ3 with H = 1 has been implemented, as shown in Fig. 8(c), which intuitively illustrates the evolution of stable (green line) and unstable (black and red lines) equilibrium evolves with the increment of κ3.
Fig.8 (a) Potential well energy V(η) vs. normalized power η with H = 1 and different values of κ3; (b) potential well energy V(η) vs. normalized power η with κ3 = 1 and different values of Hamiltonian parameter H; (c) 3D plot for potential well energy V(η) as function of normalized power η and linear phase mismatch κ3 when H = 1. Physical limits require that 0≤η≤1
To obtain the extreme points, one should find the roots from the following third-order polynomial equationand the roots can be analytically obtained as follows:with
Here, e = (κ3-6)/4, D = (q/2)2 + (p/3)3, f = [12+ 4H- (κ3 + 2)2]/24, g = (κ3 + 2)(1-H)/24, T = -q/2, p = f-e2/3, q = 2e3/27-ef/3+ g, , . The complex root implies that there is no extreme point in potential well curve. Therefore, we can discuss the property of the roots in Eq. (34) from the following three aspects:
1) If D>0, there is only one real root (y1) and two conjugated complex roots (y2 and y3), and thus there is only one equilibrium point in V(η) vs. η curve.
2) If D<0, the roots obtained are all real number, and hence there are three equilibrium points in V(η) vs. η curve.
3) If D = 0 and p = q = 0, the three roots are identified with η1,2,3 = -e/3, corresponding to only one equilibrium point in V(η) vs. η curve; if D = 0 and (q/2)2 = - (p/3)3≠0, the roots are , and , corresponding to two equilibrium points in V(η) vs. η curve.
From what has been discussed above, the equilibrium point in potential well curve can be robustly governed by κ3 and H, which would further influence energy exchange between the injected and sideband waves.
Characteristic of satellite
Satellite wave is one of the sideband waves with angular frequency ω4, which is adjacent to the injected wave with angular frequency ω2 and is phase conjugated to the injected wave with angular frequency ω1, as seen in Fig. 1. Figure 9(a) displays the well known Young’s double-slit experiment based on two coherent waves interfering with each other, in which the two waves experience different optical paths and thus form alternatively light (green line) and dark (black line) fringes. A similar phenomenon also takes place when satellite wave travels over an optical fiber, originated from the MFWM effect, as shown in Fig. 9(b), and the sameness and difference of the two phenomena are demonstrated in detail in Table 2.
Fig.9 (a) Comparison of Young's double-slit experiment and (b) evolvement of satellite wave over an optical fiber
Tab.2 Table 2 Comparison of Young’s double-slit experiment and evolvement of satellite wave over an optical fiber
Young’s double-slit experiment
evolvement of satellite wave
sameness
alternatively bright and dark fringes
need of two injected waves
related to phase difference of two injected waves
polarization dependence of two injected waves
difference
need of two coherent waves
free of coherent characteristic for the two waves
injected waves with a uniform wavelength
injected waves with different wavelengths
distribution of fringes depending on the injected wavelength and dimensional structure of the experimental installation
distribution of fringes depending on the injected wavelengths and powers
No new component generation
new component generation
Non-essential of optical fiber
need of optical fiber
Supposing |A1|>|A2|>>|A3|>|A4|, |A1(z)| ≈ |A1(0)| and |A2(z)| ≈ |A2(0)|, the satellite-wave amplitude can be explicitly derived from Eqs. (1)-(4),
The solution indicated in Eq. (35) is improved a lot compared to the result reported in Ref. [9], whose authors do not consider the impact of idler wave on satellite wave. Since satellite wave is a second-order FWM product, the impact of idler wave on satellite wave cannot be ignored. The big bracket in Eq. (35) indicates that pump, signal and idler waves together contribute to the amplitude of satellite wave. The satellite-wave power can be yielded from Eq. (35),with
Note that Δβ1 + Δβ2 + Δβ3 = 0, and it leads to the result of a = c. Thus, Eq. (36) can be further simplified with
The extreme points yielded from Eq. (38) are
To make clear the physical implication of Eq. (39), we define Δϕ = Δk × z = (k1 - k2)z = (a - b)z/2, where Δk represents the difference of different net phase mismatch profiles, and Δϕ denotes accumulation of phase-mismatch difference along an optical fiber, which determines the intense of satellite wave. In this case, an assumption that Δk retains invariable and z varies at will is made. At first, when z = 4lπ/(a - b), P4 ≈ 4{[(m + n)/a]sin[2alπ/(a - b)] - (m/b)sin[2blπ/(a - b)]}2, and it is remarkable that 2alπ/(a - b) -2blπ/(a - b) = 2lπ, which shows that the two terms [(m + n)/a]sin[2alπ/(a - b)] and (m/b)sin[2blπ/(a - b)] are in-phase. Hence, the whole process is identified to a interfered-destructive case, giving rise to a dark fringe emerging at the point of z = 4lπ/(a - b). Secondly, when z = 2(2l + 1)π/(a - b), P4 ≈ 4{[(m + n)/a]sin[a(2l + 1)π/(a - b)] + (m/b)sin[b(2l + 1)π/(a - b)]}2. We note that a(2l + 1)π/(a - b)-b (2l + 1)π/(a - b) = (2l + 1)π, so the two terms [(m + n)/a]sin[a(2l + 1) π/(a - b)] and (m/b)sin[b(2l + 1)π/(a - b)] are out of phase. This issue, however, also equals to a interfered-destructive case, because of the “+” sign located between the two terms. When z = (2l + 1)π/(a - b), P4 ≈ 4(m + n)2/a2sin2[a(2l + 1)π/(2a - 2b)] + 4m2/b2sin2[b(2l + 1)π/(2a - 2b)], which can be further simplified as: P4 ≈ 2(m + n)2/a2 + 2m2/b2, and it is identified to a interfered-constructive term.
Furthermore, Eq. (38) can be transformed into another form, that iswith
Mathematically, Eqs. (40) and (41) point out that the expression for the power of satellite wave can be approximately regarded as a complex of three cosine functions, and the irregular period T for satellite power vs. fiber length curve is controlled by |Ta|, |Tb| and |Ta - b|, respectively. Note that, Ta, Tb and Ta - b are related to a and b, which means that one can adjust the injected wavelength detune to change the irregular period T. Further investigation shows that the irregular period can be yielded through the manipulation of T = [|Ta|, |Tb|, |Ta - b|], while the middle bracket denotes the operation of searching the least common multiple in an array.
Figure 10(a) shows the evolvement of satellite wave over an optical fiber, under the condition of λ1 = 1555 nm, λ2 = 1545 nm, P10 = 20 mW, P20 = 10 mW. The red dash line and green dots represent the numerical simulation and the analytical solution, respectively, and they are in high agreement. It can be clearly seen that satellite power varies with fiber length in a sine-like trajectory, in line with the analysis for Eqs. (38) and (39). By changing the boundary condition with λ1 = 1555 nm, λ2 = 1548 nm, P10 = 10 mW, P20 = 5 mW and L = 1000 m, Fig. 10(b) again displays satellite power as a function of fiber length, and the interfered-constructive dots marked with blue color meet the condition of z = (2l + 1)π/(a - b). Moreover, interfered-destructive points marked with green and red dots satisfy the condition of z = 4lπ/(a - b) and z = 2(2l + 1)π/(a - b), respectively. Also, the numerical simulation is in line with the analytical solution in this case. In fact, there have been some common points between the property of satellite wave evolvement and the coherent-superimposed principal. First, Eq. (38) for satellite power is similar with the expression for the intensity of interference pattern [17,18]. Second, for the coherent-superimposed principal, constructive and destructive points depend critically on the optical-path difference (identified with the phase difference) of two injected waves when they reach some point in space. Similarly, when the satellite wave travels to a specific location in an optical fiber, difference of different sets of phase-mismatch profiles in MFWM processes determines whether the fringe is bright or dark, corresponding to the constructive or destructive point.
Fig.10 Simulations for satellite power as a function of fiber length. (a) The red dash line represents the numerical result while the green dots denote the analytical solution. The related parameters are λ1 = 1555 nm, λ2 = 1545 nm, P10 = 20 mW, P20 = 10 mW and L = 500 m; (b) The green and red dots correspond to interfered-destructive points satisfying the condition of z = 4lπ/(a - b) and z = (4l+ 2)π/(a - b) respectively, the blue dots correspond to interfered-constructive points satisfying the condition of z = (2l+ 1)π/(a - b). The related parameters are λ1 = 1555 nm, λ2 = 1548 nm, P10 = 10 mW, P20 = 5 mW and L = 1000 m
Figure 11 illustrates satellite power as a function of fiber length with different injected wavelengths to check how the injected wavelength influences the periodic characteristic, as described by T = [|Ta|, |Tb|, |Ta - b|]. When λ1 = 1555 nm and λ2 = 1548 nm, one can easily obtain Ta = -176.9358 m, Tb = -88.4493 m, Tc = 176.8615 m, and thus T ≈ 176.8615 m, which matches well with the numerical simulation shown in Fig. 11(a). Moreover, Figs. 11(b), 11(c), 11(d) further demonstrate that, the period changes dramatically when setting different injected wavelengths, which implies that the values of Ta, Tb and Tc significantly rely on the injected wavelength. This characteristic is very useful in the applicative regime, especially for the optical switch and the logic gate [21], that we can obtain the satellite power to a largest extent by prudently selecting the injected wavelength provided the fiber length is fixed.
Fig.11 Satellite power as function of fiber length with different injected wavelengths: (a) λ1 = 1555 nm, λ2 = 1548 nm; (b) λ1 = 1560 nm, λ2 = 1548 nm; (c) λ1 = 1550 nm, λ2 = 1545 nm; (d) λ1 = 1550 nm, λ2 = 1540 nm. The other parameters are P10 = 10 mW, P20 = 5 mW
The nonlinear interaction between four optical fields involves a single non-degenerate (Mode C) and two degenerate FWM (Mode A and B) processes in the third-order nonlinear mixings, and we emphasize Mode C by simplifying the coupled-wave equations into the eigenequations, as well as analyzing the phase-space trajectory. The existence of unstable eigensolutions has important consequences for evolution of each wave and the power exchange. Our solutions permit one to predict high-efficiency wavelength conversion for initially mismatched waves, and approximate 100% injected-wave power could be presumably coupled into the sideband wave, as long as we properly set the boundary condition. Furthermore, under the condition of equalized and approximately-equalized injected powers, we have investigated phase-space orbits intensively by analyzing the type of singular point to judge the amount of power exchange. Finally, we demonstrate the characteristic of satellite wave by yielding the solution for its power analytically from the coupled-wave equations under certain approximations, and compare the evolvement of satellite wave with the interference pattern in Young’s double-slit experiment. The above investigations can be found potential applications in 2R regeneration [7,8], optical switching [5], multi-channel all optical wavelength conversion [2], and UWB generation and transmission [22].
Acknowledgements
This paper was supported by the National Basic Research Program of China (No. 2011CB301704), the National Natural Science Foundation of China (Grant Nos. 60901006 and 61106046), and the Youth Science Funds (No. 51109215).
Geraghty D F, Lee R B, Verdiell M, Ziari M, Mathur A, Vahala K J. Wavelength conversion for WDM communication systems using four-wave mixing in semiconductor optical amplifiers. IEEE Journal of Quantum Electronics, 1997, 3(5): 1146-1155
Hansryd J, Andrekson P A, Westlund M, Li J, Hedekvist P O. Fiber-based optical parametric amplifiers and their applications. IEEE Journal of Quantum Electronics, 2002, 8(3): 506-520
Hedekvist P O, Karlsson M, Andrekson P A. Fiber four-wave mixing demultiplexing with inherent parametric amplification. Journal of Lightwave Technology, 1997, 15(11): 2051-2058
Nie Z Q, Zheng H B, Zhang Y, Zhao Y, Zuo C, Li C, Chang H, Xiao M. Experimental demonstration of optical switching and routing via four-wave mixing spatial shift. Optics Express, 2010, 18(2): 899-905
Ciaramella E, Curti F, Trillo S. All-optical signal reshaping by means of four-wave mixing in optical fibers. IEEE Photonics Technology Letters, 2001, 13(2): 142-144
Cappellini G, Trillo S. Third-order three-wave mixing in single-mode fibers: exact solutions and spatial instability effects. Journal of the Optical Society of America. B, Optical Physics, 1991, 8(4): 824-838
Trillo S, Wabnitz S, Kennedy T A. Nonlinear dynamics of dual-frequency-pumped multiwave mixing in optical fibers. Physical Review A, 1994, 50(2): 1732-1747
Zhao L, Sun J Q. Investigation of tunable trap filter utilizing intense signal four wave mixing model in highly non-linear fiber. Optics Communications, 2009, 282(14): 2975-2982
Marhic M E, Wong K K Y, Ho M C, Kazovsky L G. 92% pump depletion in a continuous-wave one-pump fiber optical parametric amplifier. Optics Letters, 2001, 26(9): 620-622
Gachet D, Brustlein S, Rigneault H. Revisiting the Young’s double slit experiment for background-free nonlinear Raman spectroscopy and microscopy. Physical Review Letters, 2010, 104(21): 213905
18
Intonti F, Riboli F, Caselli N, Abbarchi M, Vignolini S, Wiersma D S, Vinattieri A, Gerace D, Balet L, Li L H, Francardi M, Gerardino A, Fiore A, Gurioli M. Young’s type interference for probing the mode symmetry in photonic structures. Physical Review Letters, 2011, 106(14): 143901
19
Liu X M. Theory and experiments for multiple four-wave-mixing processes with multifrequency pumps in optical fibers. Physical Review A, 2008, 77: 043818
20
Hart D L, Judy A F, Roy R, Beletic J. Dynamical evolution of multiple four-wave-mixing processes in an optical fiber. Physical Review E: Statistical Physics, Plasmas, Fluids, and Related Interdisciplinary Topics, 1998, 57(4): 4757-4774
Wu B B, Fu S N, Wu J, Shum P, Ngo N Q, Xu K, Hong X, Lin J. Simultaneous implementation of all-optical OR and AND logic gates for NRZ/RZ/CSRZ ON-OFF-keying signals. Optics Communications, 2010, 283(3): 349-354
Zhao L, Sun J Q, Huang D. Photonic generation of ultrawideband signals by exploiting gain saturation of dark pump pulse with double undershoots in a highly nonlinear fiber. Optics Communications, 2011, 284(6): 1669-1676