Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Polarization-independent amplification and frequency conversion in strongly-birefringent fibers

Open Access Open Access

Abstract

The inverse modulation interaction is a degenerate four-wave mixing process in which two strong pumps drive a weak signal, whose frequency is the average of the pump frequencies. Theoretical analyses and numerical simulations of this process are made for wave frequencies that are near the zero-dispersion frequency of a fiber, in which case dispersion is unimportant, and wave frequencies that are far from the zero-dispersion frequency, in which case dispersion is important. The results show that the inverse modulation interaction in a strongly-birefringent fiber amplifies a linearly-polarized signal by an amount that depends on its phase angle, but not its polarization angle. Phase conjugation and Bragg scattering are nondegenerate four-wave mixing processes in which two strong pumps drive a weak signal and a weak idler. Studies show that phase conjugation and Bragg scattering in strongly-birefringent fibers produce polarization-independent phase-insensitive amplification and frequency conversion, respectively.

©2008 Optical Society of America

1. Introduction

Parametric devices (PDs) based on four-wave mixing (FWM) in fibers can amplify, frequency convert, phase conjugate, regenerate and sample optical signals [1, 2, 3]. However, because transmission fibers do not preserve the signal polarizations, practical PDs must operate in signal-polarization independent manners. The polarization properties of degenerate and nondegenerate FWM, in strongly-birefringent, randomly-birefringent and rapidly-spun fibers were reviewed recently [4, 5, 6]. Most of these processes depend sensitively on the signal polarization. One such polarization-dependent process is the modulation interaction (MI), in which a strong pump amplifies a weak signal and generates a weak idler. The difference between the idler and pump frequencies equals the difference between the pump and signal frequencies.

Devices based on MI (degenerate FWM) can be made polarization independent by the use of polarization-diversity schemes. In one scheme [7, 8], a polarization beam-splitter is used to split the pump and signal into components (x and y), which propagate through a fiber loop in opposite directions. The x-component of the signal experiences MI driven by the x-component of the pump, whereas the y-component of the signal experiences MI driven by the y-component of the pump. Because the x and y components of the waves counterpropagate, the MIs develop independently. If the pump is linearly-polarized at 45 ° to the axes of the beam splitter (so the pump components have equal powers), the signal components are amplified by the same amount, then recombined by the beam splitter. Hence, the output-signal power does not depend on the input-signal polarization.

In the other scheme [9, 10], the pump and signal propagate through a strongly-birefringent fiber (SBF), in the same direction. If the pump and signal frequencies are similar, the interactions between the x-component of the pump and y-component of the signal, and the y-component of the pump and x-component of the signal, are suppressed by birefringence [11]. Hence, the MIs of x-polarized and y-polarized waves develop independently. If the pump is linearly-polarized at 45° to the birefringence axes, the signal components are amplified by the same amount. Hence, the output-signal power does not depend on the input-signal polarization.

A parametric process is termed phase sensitive if the output-signal power depends on the input-signal phase. Otherwise it is termed phase insensitive. The aforementioned realizations of MI are all phase insensitive. However, there is is a degenerate FWM process in which two strong pumps drive a weak signal, whose frequency is the average of the pump frequencies. This process, which is called the inverse MI, is phase sensitive [12, 13]. At present, there is considerable interest in phase-sensitive processes, because of their ability to amplify signals without adding excess noise, and remove phase noise from differential-phase-shift-keyed signals.

In this report, detailed studies are made of the phase and polarization dependence of the inverse MI in a SBF. In the first study (Sec. 3.2), all of the wave frequencies are comparable to the zero-dispersion frequency (ZDF) of the fiber, so the effects of dispersion are unimportant. If the pumps are linearly polarized at 45° to the birefringence axes, a linearly-polarized signal experiences phase-sensitive amplification that does not depend on its polarization angle. However, because dispersion is weak, secondary FWM processes produce a cascade of secondary pumps and secondary signals (idlers) [14, 15, 16]. This cascade consumes bandwidth and depletes the primary pumps. In the second study (Sec. 4.1), most of the wave frequencies are not comparable to the ZDF, so the effects of dispersion are important. Dispersion reduces the FWM cascades, but does not eliminate them completely. Two alternative configurations are considered. In the first configuration, the pump frequencies are far from the ZDF, but the signal frequency is near the ZDF. The presence of a weak cascade does not affect the signal evolution significantly. If the pumps are linearly polarized at 45 ° to the birefringence axes, the signal still experiences polarization-independent phase-sensitive amplification. In the second configuration, the pump and signal frequencies are all far from the ZDF, but close to each other. Although the signal evolves in a phase-sensitive manner, the presence of a moderate cascade limits the signal amplification significantly.

Studies are also made of nondegenerate FWM processes in which two strong pumps drive weak (signal and idler) sidebands (Sec. 4.2). In phase conjugation (PC), the sum of the sideband frequencies equals the sum of the pump frequencies, whereas in Bragg scattering (BS), the difference between the sideband frequencies equals the difference between the pump frequencies [17]. PC in a SBF provides phase-insensitive amplification, whereas BS in a SBF provides phase-insensitive frequency conversion. Dispersion affects both processes. If the pumps are polarized at 45° to the birefringence axes, signals experience polarization-independent amplification and frequency conversion. The weak cascades that occur do not degrade significantly the properties of PC and BS.

2. Wave propagation in a strongly-birefringent fiber

Light-wave propagation in a strongly-birefringent fiber (SBF) is modeled by the coupled Schrödinger equations (CSEs)

zX=iβx(iτ)X+i(γsX2+γcY2)X,
zY=iβy(iτ)Y+i(γcX2+γsY2)Y,

where z is distance, z=/∂z, X and Y are the amplitude (polarization) components of the wave and βx and βy are the dispersion functions of the fiber. In the frequency domain β(ω)=Σn≥1βn(ωc)ωn/n!, where ωc is the carrier frequency of the wave and ω is the difference between the actual and carrier frequencies. To convert from the frequency domain to the time domain one replaces ω by i∂τ, where τ=t-βaz is the retarded time and βa=[β 1x(ωc)+β 1y(ωc)]/2 is the average of the group slownesses. The self-nonlinearity coefficient γs=γK, where γK is the Kerr coefficient, and the cross-nonlinearity coefficient γc=2γK/3 [18, 19].

3. Adjacent waves

Suppose that the wave frequencies are all close to the zero-dispersion frequency (ZDF, ω 0) of the fiber. Then the effects of second- and higher-order dispersion can be neglected. In this (adjacent-wave) case, the CSEs (1) and (2) reduce to the coupled-component equations (CCEs)

(z+βdt)X=i(γsX2+γcY2)X,
(zβdt)Y=i(γcX2+γsY2)Y,

where βd=[β 1x(ωc)-β 1y(ωc)]/2 is the difference between the group slownesses. It follows from Eqs. (3) and (4) that

(z+βdt)Px=0,
(zβdt)Py=0,

where the powers Px=|X|2 and Py=|Y|2. The power of each component is constant, in a frame moving with the group speed of that component. No power is exchanged between the components. Define the retarded times τx=τ-βdz and τy=τ+βdz. Then, by using the method of characteristics [20], one finds that

X(τx,z)=X(τx,0)exp{iγsPx(τx,0)z+iγc0zPy(τx+2βdz,0)dz},
Y(τy,z)=Y(τy,0)exp{iγc0zPx(τy2βdz,0)dz+iγsPy(τy,0)z}.

Notice that the formula for Y can be deduced from the formula for X by interchanging the subscripts x and y, and changing the sign of the walk-off parameter βd. For multiple-frequency inputs, the self- and cross-nonlinearities in Eqs. (7) and (8) have time-independent parts, which produce self-phase modulation (SPM) and cross-phase modulation (CPM), respectively, and time-dependent parts, which produce scalar and vector FWM. (SPM and scalar FWM involve X and Px, or Y and Py, whereas CPM and vector FWM involve X and Py, or Y and Px.) The effects of SPM and scalar FWM accumulate with distance, as do the effects of CPM. However, the effects of vector FWM, which depend on the integral of a periodic function, are bounded and can be neglected for long distances. Apart from CPM, which only produces time-independent phase shifts, the components evolve nearly independently.

3.1. Phase-insensitive cascade

This general conclusion is reinforced by a specific example. Consider the two-frequency boundary (initial) conditions

X(τ,0)=ρ0xexp(iϕ0x)+ρ1xexp(iϕ1x),
Y(τ,0)=ρ0yexp(iϕ0y)+ρ1yexp(iϕ1y),

where the phases ϕjk=-ωjτ+ϕjk(0). These conditions correspond to two waves (0 and 1) with arbitrary input powers, polarizations and phases. It follows from Eqs. (9) and (10) that

Px(τx,0)=ρ0x2+ρ1x2+2ρ0xρ1xcos(ϕ0xϕ1x).
Py(τy,0)=ρ0y2+ρ1y2+2ρ0yρ1ycos(ϕ0yϕ1y).

Each input power has the time average ρ 2 0+ρ 2 1, which produces phase shifts, and the contribution 2ρ 0 ρ 1 cos(ϕ 0-ϕ 1), which oscillates at the difference frequency ωd=ω 1-ω 0 and produces other frequency components (harmonics). It also follows that

0zPx(τy2βdz,0)dz=(ρ0x2+ρ1x2)zρ0xρ1x{sin[ωd(τy2βdz)ϕdx]
sin[ωdτyϕdx]}/βdωd.
0zPy(τx+2βdz,0)dz=(ρ0y2+ρ1y2)z+ρ0yρ1y{sin[ωd(τx+2βdz)ϕdy]
sin[ωdτxϕdy]}/βdωd,

where the phase differences ϕdj=ϕ 1j(0)-ϕ 0j(0). The time-independent terms in Eqs. (13) and (14) grow linearly with z, whereas the time-dependent terms are bounded. The relative size of these terms is βdωdz. The parameters βd=50 ps/Km (which corresponds to the index difference δn=3×10-5), ωd=0.3 Tr/s (which corresponds to the real frequency difference fd=50 GHz) and z ~1 Km are representative of current experiments. For these parameters, βdωdz ~15, so the time-dependent terms (which produce vector FWM) can be neglected. Similar conclusions apply to multiple-frequency inputs, because the time-dependent contributions to Px and Py oscillate at the (pairwise) beat frequencies ωj-ωk. Notice that the remaining terms do not depend on βd, so the formulas for X and Y are related by the aforementioned interchange of indices. Consequently, only formulas for X will be stated explicitly in this report.

FWM cascades produced by two inputs were analyzed in [14, 15]. By proceeding as described therein, one can write the output amplitude in the form

X(τx,z)=nXn(z)exp(iϕnx+iψx),

where the time-dependent phases

ϕnx(τx)=ϕ0x(τx)+n[ϕ1x(τx)ϕ0x(τx)]

and the time-independent phase

ψx(z)=γs(ρ0x2+ρ1x2)z+γc(ρ0y2+ρ1y2)z.

The harmonics

Xn(z)=inρ0xJn(ζx)+in1ρ1xJn1(ζx),

where the distance variable ζx=2γsρ 0x ρ 1x z. The phase of the nth harmonic depends on the (common) phase reference ϕ 0x and the phase difference ϕ 1x-ϕ 0x. However, the power of the harmonic does not depend on either phase. In this sense, each two-input cascade is phase independent. Equations (15)(18), and their analogs for the y-components, show that the components of the nth harmonic have different time-independent phases and their phase modulations arrive at different times.

In general, the powers of the output harmonics depend sensitively on the polarizations of the inputs. However, suppose that the pump amplitudes ρ 0x=ρ 0y=ρ 0 (so the x- and y-polarized processes are driven equally), and the signal amplitudes ρ 1x,ρ 1yρ 0. Then, in the linear regime (ζ j≪1), the signal and primary-idler components

X1(1+iγsρ02z)ρ1x,X1(iγsρ02z)ρ1x.

It follows from Eqs. (19), and their analogs for Y 1 and Y -1, that the signal and primary-idler powers do not depend on the signal polarization (or the pump and signal phases). In this sense, the primary FWM process is signal-polarization independent.

Because solutions (18) and (19) depend on the dimensionless parameters ρ/ρ0 and γsρ 0 z, it is customary to measure the mode powers relative to pump power ρ 2 0 and distance relative to the interaction distance 1/γsρ 2 0. However, to facilitate comparisons between theoretical predictions and numerical results, we chose to measure the pump and signal powers relative to a reference power of 1 mW. The cascade produced by a strong pump and a weak signal is illustrated in Fig. 1, for the initial conditions ρ 0=22.4, ρ 1x=0.283 and ρ 1y=0.141. The input spectrum (γsρ 2 0 z=0) is shown in Fig. 1(a). The pump is polarized at 45° to the slow (x) axis of the fiber, whereas the signal is polarized at about 27°. The output spectrum (γsρ 2 0 z=10) is shown in Fig. 1(b). Not only has the signal been amplified by 20 dB, but a primary idler has also been generated, whose power is comparable to the signal power. The 6-dB difference between the x and y components of the output signal and primary idler equals the difference between the components of the input signal. This result shows that both polarizations are amplified by the same amount: In the linear regime signal amplification and idler generation are polarization independent. However, the differences between the components of the neighboring idlers are 12 dB, which shows that the generation of secondary idlers is polarization dependent.

 figure: Fig. 1.

Fig. 1. Normalized mode powers [Eq. (18)] plotted as functions of mode number for the case in which ρ0=22.4, ρ 1x=0.283 and ρ 1y=0.141: (a) γsρ 2 0 z=0.0 and (b) γsρ 2 0 z=10.0.

Download Full Size | PDF

Red bars denote x components, whereas blue bars denote y components.

To test the theoretical predictions, we made numerical simulations based on the CSEs (1) and (2), for β 3=0.03 ps3/Km, β 4=-3×10-4 ps4/Km and γs=10/Km-W. The pump powers Px=Py=0.5 W and the signal powers Psx=0.08 and Psy=0.02 mW (which correspond to the pump and signal amplitudes of Fig. 1). The pump frequency was the ZDF and the signal-pump frequency difference was 0.31 Tr/s (50 GHz), which is typical of 10-Gb/s communication systems. The simulation results, which are shown in Fig. 2, confirm the theoretical predictions.

 figure: Fig. 2.

Fig. 2. (a) Simulated mode powers [Eqs. (1) and (2)] plotted as functions of frequency for the case in which the pump powers Px=Py=0.5 W and the signal powers Psx=0.08 and Psy=0.02 mW: (a) z=0 Km and (b) z=2.0 Km. Red bars denote x components, whereas blue bars denote y components.

Download Full Size | PDF

In Fig. 3, the simulated signal and idlers gains are plotted as functions of the (input) signal polarization [tan-1(ρ 1y/ρ 1x)], for the same parameters as Fig. 2. The gain variations are about ±1%, so the gain is nearly polarization independent, as predicted.

 figure: Fig. 3.

Fig. 3. Simulated gain [Eqs. (1) and (2)] of the signal (solid curve) and idler (dashed curve) plotted as functions of the signal polarization, for the same parameters as Fig. 2.

Download Full Size | PDF

3.2. Phase-sensitive cascade

Now consider the three-frequency input conditions

X(τ,0)=ρexp(iϕ1)+ρ0xexp(iϕ0)+ρexp(iϕ1).

These conditions correspond to two strong pumps (-1 and 1), which are linearly polarized (LP) at 45° to the birefringence axes, and a weak signal (0), which is LP at an arbitrary angle. The pump powers are equal and the signal frequency is the average of the pump frequencies [ω 0=(ω 1+ω -1)/2]. This process is called the inverse modulation interaction (MI). Define the frequency difference ωd=(ω 1-ω -1)/2, the phase average ϕa=[ϕ 1(0)+ϕ -1(0)]/2 and the phase difference ϕd=[ϕ 1(0)-ϕ -1(0)]/2. Then, by measuring frequency relative to the average frequency ω 0, phase relative to the reference phase ϕa, and time relative to the reference time ϕd/ωd, one can rewrite condition (20) in the simpler form

X(τ,0)=ρexp(iωdτx)+ρ0xexp(iϕ0)+ρexp(iωdτx),

where ϕ 0 is the input phase.

FWM cascades produced by three inputs were also analyzed in [14, 15]. By proceeding as described therein, one can write the output amplitude in the form

X(τx,z)=nXn(z)exp(inωdτx+iψx),

where the time-independent phase

ψx=2(γs+γc)ρ2z+(γsρ0x2+γcρ0y2)z,

and the harmonics

Xn(ζ)=minmJm(ζ)[ρ0xexp(iϕ0)Jn2m(εxζ)2iρJn2m(εxζ)],

where ζ=2 ρ 2 z, εx=2(ρ 0x/ρ)cosϕ 0 and Jl(y)=dJl(y)/dy. Equation (24) and its analog forYn show that the output harmonics depend on the polarization of the input signal in two ways: First, because γsγc, their components experience different phase shifts. Second, their components depend nonlinearly on εj(ρ 0j). However, in the linear regime (εjζ≪1), the signal-induced phase shifts are negligible, the pumps (odd harmonics) do not depend on ρ 0j, and the signal and idlers (even harmonics) are all proportional to ρ 0j. In this regime, the exact solution (24) reduces to the approximate solution

Xn(ζ)in2ρ0x{Jn2(ζ)[exp(iϕ0)+2iζcosϕ0]+Jn2(ζ)2ζcosϕ0],

where n is even. Every contribution to Xn is proportional to ρ 0x, so every contribution to Yn is proportional to ρ 0y. Both components have the same dependence on ϕ 0. Hence, the output signal and idler powers depend on the input signal phase, but not the input signal polarization: In the linear regime, polarization-independent phase-sensitive amplification is possible.

For the signal mode (0), Eq. (25) can be written in the input-output form

X0(ζ)=ρ0[μ(ζ)exp(iϕ0)+ν(ζ)exp(iϕ0)],

where the transfer functions µ(ζ)=J 0(ζ)+ζ J0(ζ)+iζJ 0(ζ) and ν(ζ)=ζJ0(ζ)+iζJ 0(ζ). The signal power attains its extremal values when 2 ϕm=tan-1[(µrνi-µiνr)/(µrνr+µiνi)], where the subscripts r and i denote real and imaginary parts, respectively. The first-quadrant value of 2 φm corresponds to the maximal gain (|µ|+|ν|)2, whereas the third-quadrant value corresponds to the minimal gain (|µ|-|ν|)2. By using the aforementioned formulas for the transfer functions, one finds that

2ϕm(ζ)=tan1{J02(ζ)ζ[J02(ζ)+J12(ζ)]J0(ζ)J1(ζ)}.

The signal gain is plotted as a function of (input) signal phase and distance in Fig. 4. For short distances the gain is maximal when the input phase ϕ 0π/4, whereas it is minimal when ϕ 0≈3π/4. As the distance increases, the phases required for maximal and minimal gain tend nonmonotonically to 0 and π/2, respectively [Fig. 4(a)]. For γsρ 2 z=2 the aforementioned phases are 0.124 and 1.69, respectively [Fig. 4(b)].

 figure: Fig. 4.

Fig. 4. (a) Signal gain [Eq. (26)] plotted as a function of signal phase and distance (γsρ 2 z). Light and dark regions correspond to high and low gains, respectively. The contour spacing is 0.5. (b) Gain plotted as a function of phase for the case in which γsρ 2 z=2.0.

Download Full Size | PDF

The cascades produced by two strong pumps and a weak signal are illustrated in Figs. 5 and 6, for cases in which ρ=22.4 and ρ 0=0.316. The input spectrum is shown in Fig. 5(a). The output spectra for ϕ 0=0.124 and 1.69 are shown in Figs. 6(a) and 6(b) respectively. In the former case the signal was amplified by 10 dB, whereas in the latter it was attenuated by 35 dB. Despite the presence of many secondary pumps (even modes) and idlers (odd modes), the signal experiences significant phase-sensitive amplification and attenuation [14].

 figure: Fig. 5.

Fig. 5. Input power plotted as a function of mode number (frequency) for the case in which there are two pumps and one signal. (a) ρ=22.4 and ρ 0=0.316. (b) P=0.5 W and P 0=0.1 mW. Red bars represent the odd harmonics (pumps), whereas blue bars represent the even harmonic (signal).

Download Full Size | PDF

 figure: Fig. 6.

Fig. 6. Normalized mode power [Eq. (24)] plotted as a function of mode number for cases in which ρ=22.4, ρ 0=0.316 and γsρ 2 z=2.0. (a) For ϕ 0=0.124 the signal is amplified by 10 dB. (b) For ϕ 0=1.69 the signal is attenuated by 35 dB. Red bars represent the odd harmonics (pumps), whereas blue bars represent the even harmonics (signal and idlers).

Download Full Size | PDF

To test the theoretical predictions, we made scalar (one-polarization) and vector (two-polarization) simulations based on Eqs. (1) and (2), for β 3=0.03 ps3/Km, β 4=-3×10-4 ps4/Km and γs=10/Km-W. The pump powers Px=Py=0.5Wand the signal powers Psx=0.08 and Psy=0.02 mW (which correspond to the pump and signal amplitudes of Fig. 6). The average pump frequency was the ZDF, and the signal-pump frequency difference was 0.31 Tr/s (50 GHz). The simulation results are shown in Fig. 7. For amplification, the theoretical predictions and numerical results agree quantitatively [Figs. 6(a) and 7(a)], whereas for attenuation, they agree qualitatively [Figs. 6(b) and 7(b)]. Each harmonic in a cascade is continuously augmented or diminished by FWM processes. For the signal to be attenuated significantly, the processes that affect it must all diminish it simultaneously (cancel perfectly). Small amounts of dispersion (which are present in the simulations), modify these processes slightly and prevent this perfect cancellation. Overall, the simulation results confirm the theoretical predictions. Significant phase-sensitive amplification and attenuation was observed in a recent experiment [16], in which the inverse MI was used to phase-regenerate a signal.

 figure: Fig. 7.

Fig. 7. Simulated mode power [Eq. (1)] plotted as a function of mode frequency for cases in which P=0.5 W, P 0=0.1 mW and z=0.4 Km. (a) For ϕ 0=0.124 the signal is amplified by 10 dB. (b) For ϕ 0=1.69 the signal is attenuated by 20 dB.

Download Full Size | PDF

In Fig. 8 the simulated signal gain is plotted as a function of (input) signal phase and polarization [tan-1(ρ 0y/ρ 0x)], for the same parameters as Fig. 7. The phase dependence is nearly independent of polarization [Fig. 8(a)], and the polarization dependence is weak for the phases associated with maximal amplification and attenuation [Fig. 8(b)]. Thus, the gain is nearly polarization independent, as predicted.

 figure: Fig. 8.

Fig. 8. (a) Simulated signal gain [Eq. (1)] plotted as a function of signal phase for the signal polarizations θ=0 (solid curve) and θ=π/4 (dashed curve). (b) Simulated gain [Eqs. (1) and (2)] plotted as a function of polarization for the phases ϕ=0.124 (solid curve) and ϕ=1.69 (dashed curve). The other parameters are the same as those of Fig. 7.

Download Full Size | PDF

The ZDFs of real fibers vary with distance [21]. Phase-sensitive cascades should be affected only weakly by ZDF variations, because the frequencies of the harmonics are all close to the average ZDF. Although dispersion varies with distance, it remains weak at these frequencies.

4. Distant waves

Now suppose that the wave frequencies are not all close to the ZDF or each other. Then the effects of dispersion must be included. In this (distant-wave) case, one uses separate equations for each polarization and frequency component to study the wave interaction. Consider the interaction of four waves (1–4). By substituting the ansätze

X(t,z)=j=14Xj(z)exp(iωjt),
Y(t,z)=j=14Yj(z)exp(iωjt)

in the CSEs (1) and (2), and collecting terms of like frequency, one finds that

zX1=iβx(ω1)X1+iγs(X12+2X22+2X32+2X42)X1
+iγsX3*X22+i2γsX4*X2X3
+iγc(Y12+Y22+Y32+Y42)X1+iγcY1(Y2*X2+Y3*X3+Y4*X4)
+iγcY3*Y2X2+iγcY4*(Y2X3+Y3X2),
zX2=iβx(ω2)X2+iγs(X22+2X32+2X42+2X12)X2
+i2γsX2*X3X1+iγsX4*X32+i2γsX3*X4X1
+iγc(Y22+Y32+Y42+Y12)X2+iγcY2(Y3*X3+Y4*X4+Y1*X1)
+iγcY2*(Y3X1+Y1X3)+iγcY4*Y3X3+iγcY3*(Y4X1+Y1X4).

Equations (30) and (31) apply to the special case in which ω 2-ω 1=ω 3-ω 2=ω 4-ω 3. Similar equations apply to the general case in which the frequencies are arbitrary. One can deduce the equations for Y 1 and Y 2 from Eqs. (30) and (31), respectively, by interchanging X j and Yj. One can also deduce the equations for A 3 and A 4 from the equations for A 1 and A 2 by interchanging 1 and 4, and 2 and 3, respectively. (A is X or Y.) Collectively, these equations are called the coupled-mode equations (CMEs).

The first term on the right side of Eq. (30) represents a linear phase shift, caused by dispersion, whereas the second and fifth terms represent nonlinear phase shifts, caused by SPM and CPM. The third and fourth terms represent scalar FWM processes in which 2 π 2xπ 3x+π 1x and π 2x+π 3xπ 4x+π 1x, respectively, where πjk represents a photon with frequency ωj that is polarized parallel to the k-axis. The seventh term represents the vector FWM process in which π 2x+π 2yπ 3y+π 1x, and the eighth term represents the vector FWM processes in which π 2y+π 3xπ 4y+π 1x and π 2x+π3yπ 4y+π 1x. The sixth term represents a processes in which π 1x+π 2yπ 1y+π 2x, and similar processes that involve waves 3 and 4. Analyses show that these processes are cross-polarization rotations (CPRs). Similar statements apply to the terms in Eq. (31).

Define the total x- and y-polarized powers Px4 j=1|Xj|2 and Py4 j=1|Yj|2, respectively. Phase shifts do not change Px and Py. Scalar-FWM processes exchange energy between different x-polarized components or different y-polarized components. The photon-exchange equations show that, in each CPR and vector-FWM process, x- and y-polarized photons are destroyed and other x- and y-polarized photons are created. Hence, the CMEs conserve Px and Py, as do the coupled-component equations (CCEs) of Sec. 3. Power conservation is discussed in detail in the Appendix.

For parameters that are typical of many experiments, β(ω)=Σn≥1 βnωn/n!≈β 1 ω, where ω is the frequency of any wave (measured relative to a reference frequency, such as the ZDF). Let Xj(z)=X̄j(z)exp( 1x ωjz) and Yj(z)=Ȳj(z)exp( 1y ωjz), where β 1x=βd and β 1y=-βd. By making these substitutions in Eqs. (30) and (31), one obtains modified equations for X̄j and Ȳ j. The SPM and CPM terms in Eq. (30) are unaffected by this change of variables. The scalar FWM terms are multiplied by the phase factors exp[d(2 ω 2-ω 3-ω 1)] and exp[d(ω 2+ω 3-ω 4-ω 1)], both of which equal unity (because the frequency-matching conditions are satisfied by assumption). The CPR terms are multiplied by the phase factors exp[i 2 βd(ωj-ω 1)], the spatial average of which is zero. The vector FWM terms are multiplied by the phase factors exp[-d(ω 3-ω 1)], exp[-i 2 βd(ω 3-ω 1)] and exp[-i 2 β d(ω 2-ω 1)], the averages of which are zero. Similar statements apply to the terms in the modified version of Eq. (31). Hence, one can replace Eqs. (30) and (31) by the reduced equations

zX1iβx(ω1)X1+iγs(X12+2X22+2X32+2X42)X1
iγsX3*X22+i2γsX4*X2X3
+iγc(Y12+Y22+Y32+Y42)X1,
zX2iβx(ω2)X2+iγs(X22+2X32+2X42+2X12)X2
i2γsX2*X3X1+iγsX4*X32+i2γsX3*X4X1
+iγc(Y22+Y32+Y42+Y12)X2.

where the dispersion functions β(ω)=Σn≥2 βnωn/n! contain only second- and higher-order dispersion terms, and the bars were omitted for simplicity.

As stated above, the total powers Px and Py are conserved. Let Xj(z)=X̂j(z)exp[i(2 γsPx+γcPy)z] and Yj(z)=Ŷ j(z)exp[i(γcPx+2 γsPy)z]. For each polarization, every frequency component has the same phase shift. By making these substitutions in Eqs. (32) and (33), one removes the terms that produce this common shift, and obtains the reduced equations

zX1iβx(ω1)X1iγsX12X1+iγsX3*X22+i2γsX4*X2X3,
zX2iβx(ω2)X2iγsX22X2+i2γsX2*X3X1+iγsX4*X32+i2γsX3*X4X1,

where the carets were omitted for simplicity. Equations (34) and (35) describe scalar FWM processes in which 2 π 2xπ 3x+π 1x, π 2x+π 3xπ 4x+π 1x and 2 π 3xπ 4x+π 2x. Similar equations apply to the y-components of waves 1 and 2, and the x- and y-components of waves 3 and 4. These reduced CMEs imply that the x- and y-components of the waves evolve independently. Similar (reduced) CMEs can be derived for larger collections of waves (harmonics).

Equations (34) and (35) do not describe the vector FWM process associated with the seventh term on the right side of Eq. (30) and the eighth term on the right side of Eq. (31). In this process, which is often called the CPM instability [22, 23], a pump polarized at 45 ° to the birefringence axes drives an x- polarized sideband and a y-polarized sideband. Let ω denote the frequency difference between the x-polarized sideband and the pump. Then the dispersive contribution to the wavenumber mismatch is (β 1x-β 1y)ω+(β 2x+β 2y)ω 2/2. Phase matching only occurs for large frequency differences, for which the first two terms in the dispersion functions have the same magnitude. We omitted such frequencies from our analysis.

It is instructive to analyze the initial evolution of the FWM processes described by Eqs. (34) and (35). First, consider the MI 2 π 2π 3+π 1, in which two photons from a strong pump (2) are destroyed and two sideband, or signal and idler, photons (3 and 1) are created (πj is an abbreviation for π jx). By linearizing Eqs. (34) and (35), one obtains the input-output equations

X1(z)=μ(z)X1(0)+ν(z)X3*(0),
X3*(z)=ν*(z)X1(0)+μ*(z)X3*(0).

In Eqs. (36) and (37), the transfer functions

μ(z)=cos(kz)+i(δk)sin(kz),
ν(z)=i(γ̅k)sin(kz),

where the wavenumber-mismatch parameter δ=[β(ω 3)+β(ω 1)-2 β(ω 2)]/2+γsP 2, the coupling parameter γ̄=γsP 2 and the MI wavenumber k=(δ 2-γ̄2)1/2. The MI is unstable (amplifies the signal) whenever γ̄>|δ|. If the modulation frequency ω 2-ω 1 is low, only second-order dispersion is important and the instability criterion requires the pump frequency to be in the anomalous-dispersion regime [β 2(ω 2)<0]. The effects of higher-order dispersion on the existence and frequency dependence of MI gain are discussed in [17].

Second, consider the inverse MI π 3+π 12 π 2, in which photons from two strong pumps (3 and 1) are destroyed and two signal photons (2) are created. The inverse MI is characterized by the input-output equation

X2(z)=μ(z)X2(0)+ν(z)X2*(0),

where the transfer functions were defined in Eqs. (38) and (39), δ=[2 β 2(ω 2)-β(ω 3)-β(ω 1)]/2+γs(P 1+P 3)/2 and γ̄=2 γs(P 1 P 3)1/2. The inverse MI is also unstable whenever γ̄>|δ|. When one compares the formulas for the MI and inverse MI, one finds that the dispersive contributions to the mismatch parameters have opposite signs and the coupling parameters differ by a factor of 2. Hence, if the modulation frequency is low, the inverse MI exists when the pump frequency is in the anomalous- and normal-dispersion regimes, but is most unstable when the pump frequency is in the latter [β 2(ω 2)>0].

Third, consider the PC process π 2+π 3π 4+π 1, in which photons from two strong pumps (2 and 3) are destroyed, and signal and idler photons (4 and 1) are created. PC is characterized by the input-output equations

X1(z)=μ(z)X1(0)+ν(z)X4*(0),
X4*(z)=ν*(z)X1(0)+μ*(z)X4*(0),

where the transfer functions were defined in Eqs. (38) and (39), δ=[β(ω 4)+β(ω 1)-β 2(ω 2)-β(ω 3)]/2+γs(P 2+P 3)/2 and γ̄=2 γs(P 2 P 3)1/2. PC is also unstable whenever γ̄>|δ|. It is most unstable when the average of the pump frequencies is in the anomalous-dispersion regime [β 2(ωa)>0, where ωa=(ω 2+ω 3)/2]. Frequency tuning in the presence of higher-order dispersion is discussed in [17, 24].

Fourth, consider the BS process π 2+π 3π 4+π 1, in which photons from a strong pump (3) and a weak signal (2) are destroyed, and pump (1) and idler (4) photons are created. Although the photon-exchange equations for BS and PC have the same form, the identities of the pumps and sidebands are different. BS is characterized by the input-output equations

X2(z)=μ(z)X2(0)+ν(z)X4(0),
X4(z)=ν*(z)X2(0)+μ*(z)X4(0),

where the transfer functions were defined in Eqs. (38) and (39), δ=[β(ω 2)+β(ω 3)-β 4(ω 2)-β(ω 1)]/2+γs(P 1-P 3)/2, γ̄=2 γs(P 1 P 3)1/2 and k=(δ 2+γ̄2)1/2. BS is always stable: It provides coupling, but no amplification. If the pump powers are equal (P 1=P 3), a complete transfer of power from the signal to the idler occurs when the average of the signal frequency and the higher pump frequency equals the ZDF [ωa=ω 0, where ωa=(ω 2+ω 3)/2]. If the pump powers are not equal, a complete power transfer is possible when the average frequency is in either dispersion regime. Frequency tuning in the presence of higher-order dispersion is discussed in [25].

4.1. Phase-sensitive amplification

It was shown in Sec. 3.2 that the inverse MI of pumps whose frequencies are near the ZDF provides polarization-independent phase-sensitive amplification. However, because dispersion is weak near the ZDF, secondary FWM processes produce secondary pumps and idlers, which consume bandwidth and deplete the primary pumps. One can suppress the generation of secondary products by using primary pumps whose frequencies are far from the ZDF [17] (and it was for such a frequency configuration that inverse MI was proposed [12]). Equations (38)(40) were used to model the inverse MI for the fiber parameters β 3=0.03 ps3/Km, β 4=-3×10-4 ps4/Km, γ=10/Km-W and l=0.46 Km, and the pump powers P 1=P 3=0.25 W. (The slowness parameter β 1 has no effect on the inverse MI when the waves are co-polarized.) The dependence of the inverse MI on the pump and signal frequencies was studied. In Fig. 9, the signal gain is plotted as a function of the signal phase, for the pump frequencies ω 1=-8.25 and ω 3=11.75 Tr/s, and the signal frequency ω 2=1.75 Tr/s. These frequencies are measured relative to the ZDF. (Relative to the signal, the pump frequencies are ±ω d, where the frequency difference ωd=10 Tr/s.) For reference, the dispersion coefficient β 2(ω 2)=0.052. The inverse-MI equation (40) predicts that the signal gain P 2(z)/P 2(0) will attain a maximum of 20 dB when the signal phase ϕ 2(0)=0.78, and a minimum of 20 dB when ϕ 2=2.32.

 figure: Fig. 9.

Fig. 9. Signal gain plotted as a function of signal phase, for the fiber length l=0.46 Km, the pump powers P 1=P 3=0.25 Wand the signal power P 2=0.1 mW. The pump frequencies are -1.31 and 1.87 THz, and the signal frequency is 0.28 THz. The dashed curve denotes the theoretical prediction [Eq. (40)], whereas the solid curve denotes the simulation result [Eq. (1)].

Download Full Size | PDF

To test these predictions, we made scalar simulations based on Eq. (1), for the aforementioned parameters and the signal power P 2(0)=0.1 mW. In these simulations, a gain maximum of 19.3 dB was attained for ϕ 2=0.78 and a minimum of -19.3 dB was attained for ϕ 2=2.35, as shown in Fig. 9. The theoretical predictions agree quite well with the simulation results.

The small discrepancies between the theoretical and numerical results are due to the presence of weak FWM cascades. In Fig. 10, the simulated mode power [Eq. (1)] is plotted as a function of frequency (measured relative to the signal frequency), for the parameters of Fig. 9. In Fig. 10(a), for which ϕ 2=0.78, the amplified signal (9.3 dBmW) is stronger than its idlers (about -5 dBmW), whereas in Fig. 10(b), for which ϕ 2=2.32, the attenuated signal (-29.3 dBmW) is weaker than its idlers (about -25 dBmW). Pump-pump FWM does not involve the signal (and is intrinsically phase-insensitive), so the secondary pump powers are the same in both figures (about -10 dB).

 figure: Fig. 10.

Fig. 10. Simulated mode power [Eq. (1)] plotted as a function of mode frequency (measured relative to the signal frequency), for the parameters of Fig. 9. (a) When the signal phase ϕ 2=0.78, the signal is amplified by 19.3 dB. (b) When ϕ 2=2.35, the signal is attenuated by 19.3 dB. Red bars denote odd harmonics (pumps), whereas blue bars denote even harmonics (signal and idlers).

Download Full Size | PDF

To study the polarization dependence of the inverse MI, we made vector simulations based on Eqs. (1) and (2), for the parameters of Fig. 9. The signal gain is plotted as a function of the signal polarization in Fig. 11. For ϕ 2=0.78 the amplification varies by ±0.04 dB, whereas for ϕ 2=2.35 the attenuation varies by ±0.03 dB. Hence, polarization-independent (phase-sensitive) amplification is possible in a SBF, if the pump frequencies are dissimilar.

 figure: Fig. 11.

Fig. 11. Simulated signal gain [Eqs. (1) and (2)] plotted as a function of signal polarization, for the parameters of Fig. 9. The solid curve represents the signal phase ϕ 2=0.78, whereas the dashed curve represents ϕ 2=2.35.

Download Full Size | PDF

One can also suppress the generation of secondary pumps by using primary pumps and a signal, whose frequencies are all far from the ZDF (but not necessarily each other). The predictions of the three-mode model of inverse MI [Eq. (40)] are shown in Fig. 12(a), for the parameters β 2=5.0 ps2/Km, β 3=0.03 ps3/Km, β 4=-3×10-4 ps4/Km, γs=10/Km-W, l=0.46 Km, Px=0.25W (pumps) and Px=0.1 mW (signal). When the frequency difference ωd=1.0 Tr/s, dispersion cancels CPM exactly and the growth rate -ki=2 γsP is maximal. The inverse MI produces 20 dB of signal gain in a distance of 0.46 Km. When ωd=1.42, dispersion is stronger than CPM, but the signal still grows exponentially, with growth rate 31/2 γsP. Although the gain is reduced, the inverse MI still produces 18 dB of gain in 0.46 Km. When ωd=1.73, the total mismatch δ=2 γsP: The signal power increases quadratically with distance, rather than exponentially. In this case, the gain is reduced significantly, to only 13 dB. However, by increasing the distance to 0.99 Km, one can restore the gain to 20 dB (not shown). For reference, ωd=1 Tr/s corresponds to fd=159 GHz, which is comparable to the channel spacings of 40-Gb/s communication systems (100 and 200 GHz).

 figure: Fig. 12.

Fig. 12. Signal gain plotted as a function of the signal phase, for the fiber length l=0.46 Km, the pump powers P 1=P 3=0.25 W and the signal power P 2=0.1 mW. The solid, dot-dashed and dashed curves represent frequency differences of 0.16, 0.23 and 0.28 THz, respectively. (a) Theoretical predictions based on Eq. (40). (b) Simulation results based on Eq. (1).

Download Full Size | PDF

To test these theoretical predictions, we made scalar simulations based on Eq. (1), for the aforementioned parameters and the signal power P 2(0)=0.1 mW. When a frequency filter was used to restrict the simulations to three modes, the results obtained (not shown) agreed with the theoretical predictions. The simulation results obtained without a filter are shown in Fig. 12(b). In all three cases, the simulated gain depends periodically on the input phase, in qualitative agreement with the theoretical predictions. However, the simulated gain levels are much lower than their theoretical counterparts.

FWM cascades associated with the inverse MI are shown in Fig. 13, for cases in which ωd=1 Tr/s. (The other parameters are the same as those of Fig. 12.) For the signal phase ϕ 2(0)=0, the signal experiences little gain, whereas for ϕ 2=1.3, the signal experiences some gain (7 dB). These results are consistent with Fig. 12(b). However, in both cases there are significant secondary pumps (odd harmonics) and secondary signals, or idlers (even harmonics). The presence of these harmonics affects the signal evolution significantly, as shown in Fig. 12.

 figure: Fig. 13.

Fig. 13. Simulated mode power [Eq. (1)] plotted as a function of mode frequency, for the parameters of Fig. 12 and a frequency difference of 0.16 THz. (a) When the signal phase ϕ 2=0, the signal experiences little gain. (b) When ϕ 2=1.3, the signal experiences some gain. Red bars denote odd harmonics (pumps), whereas blue bars denote even harmonics (signal and idlers).

Download Full Size | PDF

It is well known that FWM is driven by nonlinearity and suppressed by dispersion (β 2 ω 2 d/2). Hence, the powers of the pump harmonics should depend inversely on the frequency difference (channel spacing). In Fig. 14, the power of the first pump harmonic (which has a frequency of 0.48 THz) is plotted as a function of distance, for three different channel spacings. In each case, the harmonic power is an oscillatory function of distance, so the measured power depends sensitively on the point at which it is measured. However, the peak power decreases as the channel spacing increases, as predicted.

 figure: Fig. 14.

Fig. 14. Simulated power [Eq. (1)] of a secondary pump plotted as a function of distance, for the parameters of Fig. 12. The solid, dot-dashed and dashed curves represent frequency differences of 0.16, 0.23 and 0.28 THz, respectively.

Download Full Size | PDF

The preceding results show that the potential signal gain is highest for a channel spacing of 1 Tr/s, but so also are the powers of the pump and signal harmonics. One can reduce the harmonic powers significantly by increasing the channel spacing to 1.73 Tr/s, but the benefits of weaker harmonics are offset by the significant reduction in gain. The best performance was obtained for the intermediate spacing of 1.42 Tr/s, for which moderate harmonic generation is accompanied by moderate gain. Consequently, we made additional simulations for this channel spacing, to determine how the gain depends on length. In Fig. 15, the signal gain is plotted as a function of the signal phase, for the fiber lengths l=0.2, 0.4, 0.6 and 0.8 Km. (The other parameters are the same as those of Fig. 4.) As the length increases, the gain increases, remains constant and decreases. The maximal gain of about 7 dB is attained for a length of about 0.5 Km [Figs. 12(b) and 13].

 figure: Fig. 15.

Fig. 15. Simulated signal gain [Eq. (1)] plotted as a function of signal phase, for the parameters of Fig. 12 and a frequency difference of 0.23 THz. The dotted, dashed, dot-dashed and solid curves represent lengths of 0.2, 0.4, 0.6 and 0.8 Km, respectively.

Download Full Size | PDF

FWM cascades associated with the inverse MI are shown in Fig. 16, for cases in which ωd=1.42 Tr/s. (The other parameters are the same as those for Fig. 12.) Figures 16(a)–16(d) correspond to the lengths l=0.2, 0.4, 0.6 and 0.8, and the signal phases ϕ 2=1.2, 1.4, 1.7 and 2.0 (maximal growth), respectively. In all cases, the secondary pumps are weaker than the primary pumps, and the idlers are weaker than the signal. (The secondary pump powers vary with distance, as shown in Fig. 13.) Nonetheless, the presence of weak FWM products is enough to limit the signal gain to about 7 dB. Further study of this frequency configuration (in which all frequencies are far from the ZDF) is required.

4.2. Phase-insensitive amplification and frequency conversion

Although the focus of this report was the inverse MI, which provides phase-sensitive amplification, the discussion of parametric interactions in Sec. 3 included PC [Eqs. (41) and (42)] and BS [Eqs. (43) and (44)]. PC provides phase-insensitive amplification, whereas BS provides phase-insensitive frequency conversion (provided that there are no input idlers). It was shown in [26] that strong birefringence suppresses vector PC, which involves x- and y-polarized pumps, and x- and y-polarized sidebands. This result, which was explained in Sec. 3, implies that scalar x-polarized PC and y-polarized PC are independent processes. It was stated in [25] that birefringence inhibits vector BS, which involves x-polarized pumps and y-polarized sidebands, or x- and y-polarized pumps, and x- and y-polarized sidebands [6]. This statement, which was also explained in the Sec. 3, implies that scalar x-polarized BS and y-polarized BS are independent processes.

 figure: Fig. 16.

Fig. 16. Simulated mode power [Eq. (1)] plotted as a function of mode frequency, for the parameters of Fig. 12 and a frequency difference of 0.23 THz. Figures 16(a)–16(d) correspond to lengths of 0.2, 0.4, 0.6 and 0.8, respectively. Red bars denote odd harmonics (pumps), whereas blue bars denote even harmonics (signal and idlers).

Download Full Size | PDF

Equations (38) and (41) were used to model PC for the fiber parameters β 3=0.03 ps3/Km, β 4=-3×10-4 ps4/Km, γ=10/Km-W and l=1.0 Km, and the pump powers P 2=P 3=0.25 W. (The slowness parameter β 1 has no effect on PC when the waves are co-polarized.) The dependence of the signal gain P 1(z)/P 1(0) on the pump and signal frequencieswas studied [17]. In Fig. 17, the gain is plotted as a function of the signal frequency, for the pump frequencies ω 2=-21.1 and ω 3=22.9 Tr/s (measured relative to the ZDF). For reference, β 2(ωa)=0.028, where ωa=(ω 2+ω 3)/2 is the average of the pump frequencies. The PC equations predict that signals with a wide range of frequencies will be amplified by about 36 dB. Outer gain bands (ω 1<ω 2<ω 3<ω 4) and inner bands (ω 2<ω 1<ω 4<ω 3) exist simultaneously because β 4<0 [24].

 figure: Fig. 17.

Fig. 17. Signal gain [Eqs. (38) and (41)] plotted as a function of signal frequency, for a fiber length of 1.0 Km. The pump frequencies are -3.35 and 3.65 THz (solid lines), and the (common) pump power is 0.25 W.

Download Full Size | PDF

To test the theoretical predictions, we made scalar and vector simulations based on the Eqs. (1) and (2), for the aforementioned fiber parameters (in addition to which β 1x=50 and β 1y=-50 ps/Km) and pump frequencies. In the first two simulations, the pumps were either x-polarized (P 2x=P 3x=0.25 and P 2y=P 3y=0.0 W) or y-polarized (P 2x=P 3x=0.0 and P 2y=P 3y=0.25 W) and x- or y-polarized noise sources (with powers of 10-5 mW per frequency bin) were used to seed PC. The spectra are shown in Fig. 18. Their main features (gain levels and locations of the gain maxima) are consistent with the theoretical predictions. For frequencies that are comparable to the pump frequencies, the predicted gains are higher than the simulated gains, because the PC equations omit the effects of sidebands produced by simultaneous MI and BS [17].

 figure: Fig. 18.

Fig. 18. Simulated sideband power [Eqs. (1) and (2)] plotted as a function of frequency. The pump frequencies are -3.35 and 3.65 THz (solid lines), the (common) pump power is 0.25 W and the noise power is -50 dBmW. (a) Pumps and sidebands are x-polarized. (b) Pumps and sidebands are y-polarized.

Download Full Size | PDF

In the third simulation both pumps were polarized at 45 ° to the fiber axes (P 2x=P 2y=0.25 and P 3x=P 3y=0.25 W). The spectra, which are shown in Fig. 19, are almost identical to the spectra of Fig. 18. (Differences exist between the simulated gains for individual frequencies, because different noise seeds were used.) This result confirms that polarization-independent (phase-insensitive) amplification is possible in a SBF, for a wide range of signal frequencies.

 figure: Fig. 19.

Fig. 19. Simulated sideband power [Eqs. (1) and (2)] plotted as a function of frequency. The pump frequencies are -3.35 and 3.65 THz (solid lines), the (common) pump power is 0.25 W in each polarization and the noise power is -50 dBmW in each polarization. (a) x-polarized power. (b) y-polarized power.

Download Full Size | PDF

Equations (39) and (44) were used to model BS for the fiber parameters β 3=0.03 ps3/Km, β 4=-3×10-4 ps4/Km, γ=10/Km-W and l=0.31 Km, and the pump powers P 1=P 3=0.25 W. (β 1 has no effect on BS when the waves are co-polarized.) The dependence of the signal-to-idler conversion efficiency P 4(z)/P 2(0) on the pump and signal frequencies was studied [25]. In Fig. 20, the conversion efficiency is plotted as a function of the lower pump frequency ω 1, for the signal frequency ω 2=-9.26 Tr/s and higher pump frequency ω 3=10.74 Tr/s (measured relative to the ZDF). For reference, β 2(ωa)=0.022, where ωa=(ω 2+ω 3)/2 is the average of the signal frequency and and the higher pump frequency. The BS equations predict that the conversion efficiency will be higher than -1 dB for a wide range of lower pump and idler frequencies (so the idler frequency can be tuned). The hollow circles denote lower pump and idler frequencies of -19.26 and 20.74 Tr/s, respectively, for which the predicted conversion efficiency is -0.36 dB, whereas the solid circles denote frequencies of -29.26 and 30.74 Tr/s, for which the conversion efficiency is -0.25 dB.

 figure: Fig. 20.

Fig. 20. Signal-to-idler conversion efficiency [Eqs. (39) and (44)] plotted as a function of the lower pump frequency, for a fiber length of 0.31 Km. The signal frequency is -1.47 THz (dashed line), the higher pump frequency is 1.71 THz (solid line) and the (common) pump power is 0.25 W. The hollow circles denote lower pump and idler frequencies of -3.07 and 3.3 THz, respectively, whereas the solid circles denote frequencies of -4.66 and 4.89 THz.

Download Full Size | PDF

To test the theoretical predictions, we made scalar and vector simulations based on Eqs. (1) and (2), for the aforementioned fiber parameters (in addition to which β 1x=50 and β 1y=-50 ps/Km), and pump and signal frequencies. The spectra shown in Fig. 21 were obtained for x-polarized pumps and sidebands (P 1x=P 3x=0.25 and P 1y=P 3y=0.0 W). For the lower pump frequency ω 1=-19.26 Tr/s and the idler frequency ω 4=10.74 Tr/s, a conversion efficiency of -0.50 dB was obtained [Fig. 21(a)]. This conversion efficiency is 0.14 dB lower than the predicted value, because of weak coupling to secondary modes. The figure shows two (pump) modes produced by pump-pump FWM, and two (idler) modes produced by MI and PC. For the frequencies ω 1=-29.26 Tr/s and ω 4=20.74 Tr/s, a conversion efficiency of 0.11 dB was obtained [Fig. 21(b)]. This conversion efficiency is 0.36 dB higher than predicted, because of weak (phase-sensitive) coupling to secondary modes. The figure shows two (pump) modes produced by pump-pump FWM. In both cases, the signal is attenuated by about 10 dB, so the corrections to the conversion efficiency are relatively small. Larger frequency differences between the lower pump and signal, and the higher pump (stronger dispersion), correspond(s) to weaker secondary modes.

 figure: Fig. 21.

Fig. 21. Simulated mode power [Eq. (1)] plotted as a function of mode frequency. The signal frequency is -1.47 THz, the higher pump frequency is 1.71 THz and the (common) pump power is 0.25 W. (a) The lower pump (idler) frequency is -3.07 (3.30) THz. (b) The lower pump (idler) frequency is -4.66 (4.89) THz. Red bars denote pumps, whereas blue bars denote signal and idlers.

Download Full Size | PDF

In Fig. 22, the conversion efficiency is plotted as a function of the signal polarization for pumps polarized at 45° to the fiber axes (P 1x=P 1y=0.25 and P 3x=P 3y=0.25 W). For both idler frequencies (3.30 and 4.89 THz), the polarization dependence of the conversion efficiency is very weak (≤0.01 dB). These results confirm that polarization-independent frequency conversion is possible in a SBF, for a wide range of idler frequencies.

 figure: Fig. 22.

Fig. 22. Simulated conversion efficiency [Eqs. (1) and (2)] plotted as a function of signal polarization. The signal frequency is -1.47 THz, the higher pump frequency is 1.71 THz and the (common) pump power is 0.25Win each polarization. The dashed curve represents a lower pump (idler) frequency of -3.07 (3.30) THz, whereas the solid curve represents a lower pump (idler) frequency of -4.66 (4.89) THz.

Download Full Size | PDF

The preceding results were based on the assumption that the second- and higher-order dispersion coefficients are the same for both axes. However, the ZDFs of some fibers differ by amounts of order 1 Tr/s. Such differences reduce the bandwidths over which amplification and frequency conversion are polarization independent (because the gain and conversion curves are offset in frequency). ZDF variations decrease the gain levels and conversion efficiencies [21], but should not affect the polarization dependence of PC and BS (unless the ZDFs vary independently).

5. Summary

The inverse modulation interaction (MI) is a degenerate four-wave mixing process in which two strong pumps drive a weak signal, whose frequency is the average of the pump frequencies. In this report, detailed studies were made of the phase and polarization dependence of the inverse MI in a strongly-birefringent fiber. In the first study (Sec. 3.2), all of the wave frequencies were comparable to the zero-dispersion frequency (ZDF) of the fiber, so the effects of dispersion were unimportant. For parameters that are typical of current experiments, birefringence decouples the evolution of the x and y components of the waves, where x and y denote the birefringence axes. Hence, if two pumps are linearly polarized at 45 ° to these axes, the phase-sensitive amplification (and attenuation) experienced by a linearly-polarized signal is independent of its polarization angle. This behavior was explained theoretically, in terms of the walk-off between time-dependent contributions to different components [Eqs. (7) and (8)], and verified numerically, by simulations based on the coupled Schrödinger equations [Eqs. (1) and (2)]. However, because dispersion is weak near the ZDF, secondary FWM processes produce a cascade of weak secondary pumps and signals (idlers), which consumes bandwidth and depletes the pumps.

In the second study (Sec. 4.1), most of the wave frequencies were far from the ZDF, so the effects of dispersion were important. For parameters that are typical of current experiments, dispersion reduces, but does not eliminate completely, the cascade of secondary waves produced by other FWM processes. Two alternative configurations were considered. In the first configuration, the pump frequencies were far from the ZDF, but the signal frequency was near the ZDF. The presence of a weak cascade did not affect the signal evolution significantly. If the pumps are linearly polarized at 45° to the birefringence axes, the signal still experiences polarization-independent phase-sensitive amplification. In the second configuration, the pump and signal frequencies were all far from the ZDF, but close to each other. Although the signal evolved in a phase-sensitive manner, the presence of a moderate cascade limited the signal amplification significantly.

This report focused on the inverse MI, in which two strong pumps drive one weak signal. However, similar results are obtained for a variant of the MI, in which one strong pump drives a pair of signals, whose frequencies are equally spaced about the pump frequency [27, 28]: If the pump is linearly polarized at 45° to the birefringence axes, a linearly-polarized signal pair experiences phase-sensitive amplification and attenuation that is independent of their (common) polarization angle.

Studies were also made of nondegenerate FWM processes in which two strong pumps drive weak (signal and idler) sidebands (Sec. 4.2). In phase conjugation (PC), the sum of the sideband frequencies equals the sum of the pump frequencies, whereas in Bragg scattering (BS), the difference between the sideband frequencies equals the difference between the pump frequencies. PC provides phase-insensitive amplification, whereas BS provides phase-insensitive frequency conversion. Dispersion affects both processes. If the pumps are polarized at 45 ° to the birefringence axes, signals experience polarization-independent amplification and frequency conversion. The weak cascades that occur do not degrade significantly the properties of PC and BS. However, the bandwidths of these polarization-independent processes are limited by the small difference between the zero-dispersion frequencies of the two axes.

Appendix: Power conservation by distant-wave interactions

The interaction of waves 1–4 is governed by CMEs (30) and (31) for X 1 and X 2, respectively, and their analogs for X 3, X 4 and Y 1-Y 4. The purpose of this appendix is to show that these CMEs conserve the total powers of the x- and y-components (Px and Py). Linear wavenumber shifting (caused by dispersion) and nonlinear wavenumber shifting (caused by SPM and CPM) do not change the powers of the x- and y-components of waves 1–4 (Pjx and Pjy), so these processes need not be considered explicitly.

First, consider the scalar FWM processes in which 2 π 2xπ 3x+π 1x, π 2x+π 3xπ 4x+π 1x and 2 π 3xπ 4x+π 2x. The first and third processes are degenerate, whereas the second is nondegenerate. By combining the associated terms in the CMEs, one finds that

dzX12=iγsX1*X22X3*+i2γsX1*X2X3X4*+c.c.,
dzX22=i2γsX1(X2*)2X3+i2γsX1X2*X3*X4+iγsX2*X32X4*+c.c.,
dzX32=iγsX1*X22X3*+i2γsX1X2*X3*X4+i2γsX2(X3*)2X4+c.c.,
dzX42=i2γsX1*X2X3X4*+iγSX2*X32X4*+c.c.

Equations (45)(48) imply that

dz(P1x+P2x+P3x+P4x)=0.

Although scalar FWM allows the x-components to exchange power, it conserves their total power. The properties of scalar FWM among the y-components are similar.

Second, consider the CPR process in which π 1x+π 2yπ 1y+π 2x. By combining the associated terms in the CMEs, one finds that

dzX12=iγcX1*Y1X2Y2*+c.c.,
dzY12=iγcX1Y1*X2*Y2+c.c.,
dzX22=iγcX1Y1*X2*Y2+c.c.,
dzY22=iγcX1*Y1X2Y2*+c.c.

Equations (50)(53) imply that

dz(P1x+P2x)=0=dz(P1y+P2y),
dz(P1x+P1y)=0=dz(P2x+P2y).

Not only does CPR conserve the total powers of the x- and y-components, it also conserves the total powers of waves 1 and 2. The properties of the other five CPR processes, in which πjx+πkyπjy+πkx, are similar.

Third, consider the CPM instability in which π 2x+π 2yπ 3y+π 1x. This FWM process is frequency degenerate, but polarization nondegenerate. By combining the associated terms in the CMEs, one finds that

dzX12=iγcX1*X2Y2Y3*+c.c.,
dzX22=iγcX1X2*Y2*Y3+c.c.,
dzY22=iγcX1X2*Y2*Y3+c.c.,
dzY32=iγcX1*X2Y2Y3*+c.c.

Equations (56)(59) imply that

dz(P1x+P2x)=0=dz(P2y+P3y).

Because one x-component is a pump and the other is a sideband, and one y-component is a pump and the other is a sideband, the total powers of the x- and y-components are conserved. The properties of the other CPM instabilities, in which π 2x+π 2yπ 3x+π 1y, π 3x+π 3yπ 4y+π 2x and π 3x+π 3yπ 4x+π 2y, are similar.

Fourth, consider the vector FWM process in which π 2y+π 3xπ 4y+π 1x. This process is frequency and polarization nondegenerate. By combining the associated terms in the CMEs, one finds that

dzX12=iγcX1*Y2X3Y4*+c.c.,
dzY22=iγcX1Y2*X3*Y4+c.c.,
dzX32=iγcX1Y2*X3*Y4+c.c.,
dzY42=iγcX1*Y2X3Y4*+c.c.

Equations (61)(64) imply that

dz(P1x+P3x)=0=dz(P2y+P4y).

Because one x-component is a pump and the other is a sideband, and one y-component is a pump and the other is a sideband, the total powers of the x- and y-components are conserved. The properties of the other vectorFWM processes, in which π 2y+π 3xπ 4x+π 1y, π 2x+π 3yπ 4y+π 1x and π 2x+π 3yπ 4x+π 1y, are similar. Thus, all the processes in which the waves participate conserve the total powers of the x- and y-components, as stated in Sec. 3.

References and links

1. J. Hansryd, P. A. Andrekson, M. Westlund, J. Li, and P. O. Hedekvist, “Fiber-based optical parametric amplifiers and their applications,” J. Sel. Top. Quantum Electron. 8, 506–520 (2002). [CrossRef]  

2. S. Radic and C. J. McKinstrie, “Two-pump fiber parametric amplifiers,” Opt. Fiber Technol. 9, 7–23 (2003). [CrossRef]  

3. C. J. McKinstrie, S. Radic, and A. H. Gnauck, “All-optical signal processing by fiber-based parametric devices,” Opt. Photon. News 18 (3), 34–40 (2007). [CrossRef]  

4. C. J. McKinstrie, S. Radic, and C. Xie, “Parametric instabilities driven by orthogonal pump waves in birefringent fibers,” Opt. Express 11, 2619–2633 (2003). [CrossRef]   [PubMed]  

5. C. J. McKinstrie, H. Kogelnik, R. M. Jopson, S. Radic, and A. V. Kanaev, “Four-wave mixing in fibers with random birefringence,” Opt. Express 12, 2033–2055 (2004). [CrossRef]   [PubMed]  

6. C. J. McKinstrie, H. Kogelnik, and L. Schenato, “Four-wave mixing in a rapidly-spun fiber,” Opt. Express 14, 8516–8534 (2006). [CrossRef]   [PubMed]  

7. T. Hasegawa, K. Inoue, and K. Oda, “Polarization independent frequency conversion by fiber four-wave mixing with a polarization diversity technique,” IEEE Photon. Technol. Lett. 5, 947–949 (1993). [CrossRef]  

8. K. K. Chow, C. Shu, C. Lin, and A. Bjarklev, “Polarization-insensitive widely tunable wavelength converter based on four-wave mixing in a dispersion-flattened nonlinear photonic crystal fiber,” IEEE Photon. Technol. Lett. 17, 624–626 (2005). [CrossRef]  

9. Z. Wang, N. Deng, C. Lin, and C. K. Chan, “Polarization-insensitive widely tunable wavelength conversion based on four-wave mixing using dispersion-flatened high-nonlinearity photonic crystal fiber with residual birefringence,” ECOC 2006, paper We3.P.18.

10. A. S. Lenihan and G. M. Carter, “Polarization-insensitive wavelength conversion at 40 Gb/s using birefringent nonlinear fiber,” CLEO 2007, paper CThAA2.

11. R. H. Stolen, M. A. Bösch, and C. Lin, “Phase matching in birefringent fibers,” Opt. Lett. 6, 213–215 (1981). [CrossRef]   [PubMed]  

12. C. J. McKinstrie and S. Radic, “Phase-sensitive amplification in a fiber,” Opt. Express 12, 4973–4979 (2004). [CrossRef]   [PubMed]  

13. C. J. McKinstrie, R. O. Moore, S. Radic, and R. Jiang, “Phase-sensitive amplification of chirped optical pulses in fibers,” Opt. Express 15, 3737–3758 (2007). [CrossRef]   [PubMed]  

14. C. J. McKinstrie and M. G. Raymer, “Four-wave-mixing cascades near the zero-dispersion frequency,” Opt. Express 14, 9600–9610 (2006). [CrossRef]   [PubMed]  

15. C. J. McKinstrie, S. Radic, M. G. Raymer, and L. Schenato, “Unimpaired phase-sensitive amplification by vector four-wave mixing near the zero-dispersion frequency,” Opt. Express 15, 2178–2189 (2007). [CrossRef]   [PubMed]  

16. K. Croussore and G. Li, “Phase regeneration of NRZ-DPSK signals based on symmetric-pump phase-sensitive amplification,” IEEE Photon. Technol. Lett. 19, 864–866 (2007). [CrossRef]  

17. C. J. McKinstrie, S. Radic, and A. R. Chraplyvy, “Parametric amplifiers driven by two pump waves,” IEEE J. Sel. Top. Quantum Electron. 8, 538–547 and 956 (2002). [CrossRef]  

18. C. R. Menyuk, “Nonlinear pulse propagation in birefringent optical fibers,” IEEE J. Quantum Electron. 23, 174–176 (1987). [CrossRef]  

19. C. J. McKinstrie, H. Kogelnik, G. G. Luther, and L. Schenato, “Stokes-space derivations of generalized Schrödinger equations for wave propagation in various fibers,” Opt. Express 15, 10964–10983 (2007). [CrossRef]   [PubMed]  

20. G. B. Whitham, Linear and Nonlinear Waves (Wiley, 1974).

21. M. Karlsson, “Four-wave mixing in fibers with randomly varying zero-dispersion wavelength,” J. Opt. Soc. Am. B 15, 2269–2275 (1998). [CrossRef]  

22. J. Rothenberg, “Modulational instability for normal dispersion,” Phys. Rev. A 42, 682–685 (1990). [CrossRef]   [PubMed]  

23. P. D. Drummond, T. A. B. Kennedy, J. M. Dudley, R. Leonhardt, and J. D. Harvey, “Cross-phase modulational instability in high-birefringence fibers,” Opt. Commun. 78, 137–142 (1990). [CrossRef]  

24. C. J. McKinstrie and S. Radic, “Parametric amplifiers driven by two pumps with dissimilar frequencies,” Opt. Lett. 27, 1138–1140 (2002). [CrossRef]  

25. C. J. McKinstrie, J. D. Harvey, S. Radic, and M. G. Raymer, “Translation of quantum states by four-wave mixing in fibers,” Opt. Express 13, 9131–9142 (2005). [CrossRef]   [PubMed]  

26. C. J. McKinstrie, S. Radic, and C. Xie, “Phase conjugation driven by orthogonal pump waves in birefringent fibers,” J. Opt. Soc. Am. B 20, 1437–1446 (2003). [CrossRef]  

27. R. Tang, P. Devgan, P. L. Voss, V. S. Grigoryan, and P. Kumar, “In-line frequency-nondegenerate phase-sensitive fibre parametric amplifer for fibre-optic communication,” IEEE Photon. Technol. Lett. 17, 1845–1847 (2005). [CrossRef]  

28. R. Tang, 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,” Opt. Express 13, 10483–10493 (2005). [CrossRef]   [PubMed]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (22)

Fig. 1.
Fig. 1. Normalized mode powers [Eq. (18)] plotted as functions of mode number for the case in which ρ0 =22.4, ρ 1x =0.283 and ρ 1y =0.141: (a) γsρ 2 0 z=0.0 and (b) γsρ 2 0 z=10.0.
Fig. 2.
Fig. 2. (a) Simulated mode powers [Eqs. (1) and (2)] plotted as functions of frequency for the case in which the pump powers Px =Py =0.5 W and the signal powers Psx =0.08 and Psy =0.02 mW: (a) z=0 Km and (b) z=2.0 Km. Red bars denote x components, whereas blue bars denote y components.
Fig. 3.
Fig. 3. Simulated gain [Eqs. (1) and (2)] of the signal (solid curve) and idler (dashed curve) plotted as functions of the signal polarization, for the same parameters as Fig. 2.
Fig. 4.
Fig. 4. (a) Signal gain [Eq. (26)] plotted as a function of signal phase and distance (γsρ 2 z). Light and dark regions correspond to high and low gains, respectively. The contour spacing is 0.5. (b) Gain plotted as a function of phase for the case in which γsρ 2 z=2.0.
Fig. 5.
Fig. 5. Input power plotted as a function of mode number (frequency) for the case in which there are two pumps and one signal. (a) ρ=22.4 and ρ 0=0.316. (b) P=0.5 W and P 0=0.1 mW. Red bars represent the odd harmonics (pumps), whereas blue bars represent the even harmonic (signal).
Fig. 6.
Fig. 6. Normalized mode power [Eq. (24)] plotted as a function of mode number for cases in which ρ=22.4, ρ 0=0.316 and γsρ 2 z=2.0. (a) For ϕ 0=0.124 the signal is amplified by 10 dB. (b) For ϕ 0=1.69 the signal is attenuated by 35 dB. Red bars represent the odd harmonics (pumps), whereas blue bars represent the even harmonics (signal and idlers).
Fig. 7.
Fig. 7. Simulated mode power [Eq. (1)] plotted as a function of mode frequency for cases in which P=0.5 W, P 0=0.1 mW and z=0.4 Km. (a) For ϕ 0=0.124 the signal is amplified by 10 dB. (b) For ϕ 0=1.69 the signal is attenuated by 20 dB.
Fig. 8.
Fig. 8. (a) Simulated signal gain [Eq. (1)] plotted as a function of signal phase for the signal polarizations θ=0 (solid curve) and θ=π/4 (dashed curve). (b) Simulated gain [Eqs. (1) and (2)] plotted as a function of polarization for the phases ϕ=0.124 (solid curve) and ϕ=1.69 (dashed curve). The other parameters are the same as those of Fig. 7.
Fig. 9.
Fig. 9. Signal gain plotted as a function of signal phase, for the fiber length l=0.46 Km, the pump powers P 1=P 3=0.25 Wand the signal power P 2=0.1 mW. The pump frequencies are -1.31 and 1.87 THz, and the signal frequency is 0.28 THz. The dashed curve denotes the theoretical prediction [Eq. (40)], whereas the solid curve denotes the simulation result [Eq. (1)].
Fig. 10.
Fig. 10. Simulated mode power [Eq. (1)] plotted as a function of mode frequency (measured relative to the signal frequency), for the parameters of Fig. 9. (a) When the signal phase ϕ 2=0.78, the signal is amplified by 19.3 dB. (b) When ϕ 2=2.35, the signal is attenuated by 19.3 dB. Red bars denote odd harmonics (pumps), whereas blue bars denote even harmonics (signal and idlers).
Fig. 11.
Fig. 11. Simulated signal gain [Eqs. (1) and (2)] plotted as a function of signal polarization, for the parameters of Fig. 9. The solid curve represents the signal phase ϕ 2=0.78, whereas the dashed curve represents ϕ 2=2.35.
Fig. 12.
Fig. 12. Signal gain plotted as a function of the signal phase, for the fiber length l=0.46 Km, the pump powers P 1=P 3=0.25 W and the signal power P 2=0.1 mW. The solid, dot-dashed and dashed curves represent frequency differences of 0.16, 0.23 and 0.28 THz, respectively. (a) Theoretical predictions based on Eq. (40). (b) Simulation results based on Eq. (1).
Fig. 13.
Fig. 13. Simulated mode power [Eq. (1)] plotted as a function of mode frequency, for the parameters of Fig. 12 and a frequency difference of 0.16 THz. (a) When the signal phase ϕ 2=0, the signal experiences little gain. (b) When ϕ 2=1.3, the signal experiences some gain. Red bars denote odd harmonics (pumps), whereas blue bars denote even harmonics (signal and idlers).
Fig. 14.
Fig. 14. Simulated power [Eq. (1)] of a secondary pump plotted as a function of distance, for the parameters of Fig. 12. The solid, dot-dashed and dashed curves represent frequency differences of 0.16, 0.23 and 0.28 THz, respectively.
Fig. 15.
Fig. 15. Simulated signal gain [Eq. (1)] plotted as a function of signal phase, for the parameters of Fig. 12 and a frequency difference of 0.23 THz. The dotted, dashed, dot-dashed and solid curves represent lengths of 0.2, 0.4, 0.6 and 0.8 Km, respectively.
Fig. 16.
Fig. 16. Simulated mode power [Eq. (1)] plotted as a function of mode frequency, for the parameters of Fig. 12 and a frequency difference of 0.23 THz. Figures 16(a)–16(d) correspond to lengths of 0.2, 0.4, 0.6 and 0.8, respectively. Red bars denote odd harmonics (pumps), whereas blue bars denote even harmonics (signal and idlers).
Fig. 17.
Fig. 17. Signal gain [Eqs. (38) and (41)] plotted as a function of signal frequency, for a fiber length of 1.0 Km. The pump frequencies are -3.35 and 3.65 THz (solid lines), and the (common) pump power is 0.25 W.
Fig. 18.
Fig. 18. Simulated sideband power [Eqs. (1) and (2)] plotted as a function of frequency. The pump frequencies are -3.35 and 3.65 THz (solid lines), the (common) pump power is 0.25 W and the noise power is -50 dBmW. (a) Pumps and sidebands are x-polarized. (b) Pumps and sidebands are y-polarized.
Fig. 19.
Fig. 19. Simulated sideband power [Eqs. (1) and (2)] plotted as a function of frequency. The pump frequencies are -3.35 and 3.65 THz (solid lines), the (common) pump power is 0.25 W in each polarization and the noise power is -50 dBmW in each polarization. (a) x-polarized power. (b) y-polarized power.
Fig. 20.
Fig. 20. Signal-to-idler conversion efficiency [Eqs. (39) and (44)] plotted as a function of the lower pump frequency, for a fiber length of 0.31 Km. The signal frequency is -1.47 THz (dashed line), the higher pump frequency is 1.71 THz (solid line) and the (common) pump power is 0.25 W. The hollow circles denote lower pump and idler frequencies of -3.07 and 3.3 THz, respectively, whereas the solid circles denote frequencies of -4.66 and 4.89 THz.
Fig. 21.
Fig. 21. Simulated mode power [Eq. (1)] plotted as a function of mode frequency. The signal frequency is -1.47 THz, the higher pump frequency is 1.71 THz and the (common) pump power is 0.25 W. (a) The lower pump (idler) frequency is -3.07 (3.30) THz. (b) The lower pump (idler) frequency is -4.66 (4.89) THz. Red bars denote pumps, whereas blue bars denote signal and idlers.
Fig. 22.
Fig. 22. Simulated conversion efficiency [Eqs. (1) and (2)] plotted as a function of signal polarization. The signal frequency is -1.47 THz, the higher pump frequency is 1.71 THz and the (common) pump power is 0.25Win each polarization. The dashed curve represents a lower pump (idler) frequency of -3.07 (3.30) THz, whereas the solid curve represents a lower pump (idler) frequency of -4.66 (4.89) THz.

Equations (77)

Equations on this page are rendered with MathJax. Learn more.

z X = i β x ( i τ ) X + i ( γ s X 2 + γ c Y 2 ) X ,
z Y = i β y ( i τ ) Y + i ( γ c X 2 + γ s Y 2 ) Y ,
( z + β d t ) X = i ( γ s X 2 + γ c Y 2 ) X ,
( z β d t ) Y = i ( γ c X 2 + γ s Y 2 ) Y ,
( z + β d t ) P x = 0 ,
( z β d t ) P y = 0 ,
X ( τ x , z ) = X ( τ x , 0 ) exp { i γ s P x ( τ x , 0 ) z + i γ c 0 z P y ( τ x + 2 β d z , 0 ) d z } ,
Y ( τ y , z ) = Y ( τ y , 0 ) exp { i γ c 0 z P x ( τ y 2 β d z , 0 ) d z + i γ s P y ( τ y , 0 ) z } .
X ( τ , 0 ) = ρ 0 x exp ( i ϕ 0 x ) + ρ 1 x exp ( i ϕ 1 x ) ,
Y ( τ , 0 ) = ρ 0 y exp ( i ϕ 0 y ) + ρ 1 y exp ( i ϕ 1 y ) ,
P x ( τ x , 0 ) = ρ 0 x 2 + ρ 1 x 2 + 2 ρ 0 x ρ 1 x cos ( ϕ 0 x ϕ 1 x ) .
P y ( τ y , 0 ) = ρ 0 y 2 + ρ 1 y 2 + 2 ρ 0 y ρ 1 y cos ( ϕ 0 y ϕ 1 y ) .
0 z P x ( τ y 2 β d z , 0 ) d z = ( ρ 0 x 2 + ρ 1 x 2 ) z ρ 0 x ρ 1 x { sin [ ω d ( τ y 2 β d z ) ϕ d x ]
sin [ ω d τ y ϕ d x ] } / β d ω d .
0 z P y ( τ x + 2 β d z , 0 ) d z = ( ρ 0 y 2 + ρ 1 y 2 ) z + ρ 0 y ρ 1 y { sin [ ω d ( τ x + 2 β d z ) ϕ dy ]
sin [ ω d τ x ϕ dy ] } / β d ω d ,
X ( τ x , z ) = n X n ( z ) exp ( i ϕ n x + i ψ x ) ,
ϕ n x ( τ x ) = ϕ 0 x ( τ x ) + n [ ϕ 1 x ( τ x ) ϕ 0 x ( τ x ) ]
ψ x ( z ) = γ s ( ρ 0 x 2 + ρ 1 x 2 ) z + γ c ( ρ 0 y 2 + ρ 1 y 2 ) z .
X n ( z ) = i n ρ 0 x J n ( ζ x ) + i n 1 ρ 1 x J n 1 ( ζ x ) ,
X 1 ( 1 + i γ s ρ 0 2 z ) ρ 1 x , X 1 ( i γ s ρ 0 2 z ) ρ 1 x .
X ( τ , 0 ) = ρ exp ( i ϕ 1 ) + ρ 0 x exp ( i ϕ 0 ) + ρ exp ( i ϕ 1 ) .
X ( τ , 0 ) = ρ exp ( i ω d τ x ) + ρ 0 x exp ( i ϕ 0 ) + ρ exp ( i ω d τ x ) ,
X ( τ x , z ) = n X n ( z ) exp ( i n ω d τ x + i ψ x ) ,
ψ x = 2 ( γ s + γ c ) ρ 2 z + ( γ s ρ 0 x 2 + γ c ρ 0 y 2 ) z ,
X n ( ζ ) = m i n m J m ( ζ ) [ ρ 0 x exp ( i ϕ 0 ) J n 2 m ( ε x ζ ) 2 i ρ J n 2 m ( ε x ζ ) ] ,
X n ( ζ ) i n 2 ρ 0 x { J n 2 ( ζ ) [ exp ( i ϕ 0 ) + 2 i ζ cos ϕ 0 ] + J n 2 ( ζ ) 2 ζ cos ϕ 0 ] ,
X 0 ( ζ ) = ρ 0 [ μ ( ζ ) exp ( i ϕ 0 ) + ν ( ζ ) exp ( i ϕ 0 ) ] ,
2 ϕ m ( ζ ) = tan 1 { J 0 2 ( ζ ) ζ [ J 0 2 ( ζ ) + J 1 2 ( ζ ) ] J 0 ( ζ ) J 1 ( ζ ) } .
X ( t , z ) = j = 1 4 X j ( z ) exp ( i ω j t ) ,
Y ( t , z ) = j = 1 4 Y j ( z ) exp ( i ω j t )
z X 1 = i β x ( ω 1 ) X 1 + i γ s ( X 1 2 + 2 X 2 2 + 2 X 3 2 + 2 X 4 2 ) X 1
+ i γ s X 3 * X 2 2 + i 2 γ s X 4 * X 2 X 3
+ i γ c ( Y 1 2 + Y 2 2 + Y 3 2 + Y 4 2 ) X 1 + i γ c Y 1 ( Y 2 * X 2 + Y 3 * X 3 + Y 4 * X 4 )
+ i γ c Y 3 * Y 2 X 2 + i γ c Y 4 * ( Y 2 X 3 + Y 3 X 2 ) ,
z X 2 = i β x ( ω 2 ) X 2 + i γ s ( X 2 2 + 2 X 3 2 + 2 X 4 2 + 2 X 1 2 ) X 2
+ i 2 γ s X 2 * X 3 X 1 + i γ s X 4 * X 3 2 + i 2 γ s X 3 * X 4 X 1
+ i γ c ( Y 2 2 + Y 3 2 + Y 4 2 + Y 1 2 ) X 2 + i γ c Y 2 ( Y 3 * X 3 + Y 4 * X 4 + Y 1 * X 1 )
+ i γ c Y 2 * ( Y 3 X 1 + Y 1 X 3 ) + i γ c Y 4 * Y 3 X 3 + i γ c Y 3 * ( Y 4 X 1 + Y 1 X 4 ) .
z X 1 i β x ( ω 1 ) X 1 + i γ s ( X 1 2 + 2 X 2 2 + 2 X 3 2 + 2 X 4 2 ) X 1
i γ s X 3 * X 2 2 + i 2 γ s X 4 * X 2 X 3
+ i γ c ( Y 1 2 + Y 2 2 + Y 3 2 + Y 4 2 ) X 1 ,
z X 2 i β x ( ω 2 ) X 2 + i γ s ( X 2 2 + 2 X 3 2 + 2 X 4 2 + 2 X 1 2 ) X 2
i 2 γ s X 2 * X 3 X 1 + i γ s X 4 * X 3 2 + i 2 γ s X 3 * X 4 X 1
+ i γ c ( Y 2 2 + Y 3 2 + Y 4 2 + Y 1 2 ) X 2 .
z X 1 i β x ( ω 1 ) X 1 i γ s X 1 2 X 1 + i γ s X 3 * X 2 2 + i 2 γ s X 4 * X 2 X 3 ,
z X 2 i β x ( ω 2 ) X 2 i γ s X 2 2 X 2 + i 2 γ s X 2 * X 3 X 1 + i γ s X 4 * X 3 2 + i 2 γ s X 3 * X 4 X 1 ,
X 1 ( z ) = μ ( z ) X 1 ( 0 ) + ν ( z ) X 3 * ( 0 ) ,
X 3 * ( z ) = ν * ( z ) X 1 ( 0 ) + μ * ( z ) X 3 * ( 0 ) .
μ ( z ) = cos ( kz ) + i ( δ k ) sin ( kz ) ,
ν ( z ) = i ( γ ̅ k ) sin ( kz ) ,
X 2 ( z ) = μ ( z ) X 2 ( 0 ) + ν ( z ) X 2 * ( 0 ) ,
X 1 ( z ) = μ ( z ) X 1 ( 0 ) + ν ( z ) X 4 * ( 0 ) ,
X 4 * ( z ) = ν * ( z ) X 1 ( 0 ) + μ * ( z ) X 4 * ( 0 ) ,
X 2 ( z ) = μ ( z ) X 2 ( 0 ) + ν ( z ) X 4 ( 0 ) ,
X 4 ( z ) = ν * ( z ) X 2 ( 0 ) + μ * ( z ) X 4 ( 0 ) ,
d z X 1 2 = i γ s X 1 * X 2 2 X 3 * + i 2 γ s X 1 * X 2 X 3 X 4 * + c. c . ,
d z X 2 2 = i 2 γ s X 1 ( X 2 * ) 2 X 3 + i 2 γ s X 1 X 2 * X 3 * X 4 + i γ s X 2 * X 3 2 X 4 * + c . c . ,
d z X 3 2 = i γ s X 1 * X 2 2 X 3 * + i 2 γ s X 1 X 2 * X 3 * X 4 + i 2 γ s X 2 ( X 3 * ) 2 X 4 + c . c . ,
d z X 4 2 = i 2 γ s X 1 * X 2 X 3 X 4 * + i γ S X 2 * X 3 2 X 4 * + c. c .
d z ( P 1 x + P 2 x + P 3 x + P 4 x ) = 0 .
d z X 1 2 = i γ c X 1 * Y 1 X 2 Y 2 * + c. c . ,
d z Y 1 2 = i γ c X 1 Y 1 * X 2 * Y 2 + c. c . ,
d z X 2 2 = i γ c X 1 Y 1 * X 2 * Y 2 + c. c . ,
d z Y 2 2 = i γ c X 1 * Y 1 X 2 Y 2 * + c. c .
d z ( P 1 x + P 2 x ) = 0 = d z ( P 1 y + P 2 y ) ,
d z ( P 1 x + P 1 y ) = 0 = d z ( P 2 x + P 2 y ) .
d z X 1 2 = i γ c X 1 * X 2 Y 2 Y 3 * + c. c . ,
d z X 2 2 = i γ c X 1 X 2 * Y 2 * Y 3 + c. c . ,
d z Y 2 2 = i γ c X 1 X 2 * Y 2 * Y 3 + c. c . ,
d z Y 3 2 = i γ c X 1 * X 2 Y 2 Y 3 * + c. c .
d z ( P 1 x + P 2 x ) = 0 = d z ( P 2 y + P 3 y ) .
d z X 1 2 = i γ c X 1 * Y 2 X 3 Y 4 * + c. c . ,
d z Y 2 2 = i γ c X 1 Y 2 * X 3 * Y 4 + c. c . ,
d z X 3 2 = i γ c X 1 Y 2 * X 3 * Y 4 + c. c . ,
d z Y 4 2 = i γ c X 1 * Y 2 X 3 Y 4 * + c. c .
d z ( P 1 x + P 3 x ) = 0 = d z ( P 2 y + P 4 y ) .
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.