Abstract

In a nonlinear optical fiber communication (OFC) system with signal power much stronger than noise power, the noise field in the fiber can be described by linearized noise equation (LNE). In this case, the noise impact on the system performance can be evaluated by moment-generating function (MGF) method. Many published MGF calculations were based on the LNE using continuous wave (CW) approximation, where the modulated signal needs to be artificially simplified as an unmodulated signal. Results thus obtained should be treated carefully. Reliable results can be obtained by replacing the CW-based LNE with the accurate LNE proposed by Holzlöhner et al in Ref. [1]. In this work we show that, for the case of linearized noise amplified by EDFAs, its MGF can be obtained by calculating the noise propagator directly from the accurate LNE. Our results agree well with the experimental data of multi-span DPSK systems.

© 2014 Optical Society of America

1. Introduction

The amplified spontaneous emission (ASE) noise from optical amplifiers, e.g., Erbium Doped Fiber Amplifiers (EDFAs), is one of the fundamental reasons for the bit-error-rate (BER) in an optical fiber communication (OFC) system. For an OFC system with non-negligible Kerr nonlinearity, the ASE impact evaluation is complicated due to the nonlinear interaction between signal and ASE noise. In the case of signal power much stronger than the noise power, the noise-noise beating is relatively small so that the noise field in the fiber can be approximately described by linearized noise equation (LNE), which was proposed separately in Ref. [1] and Refs. [28].

Noise propagator is a matrix used to show noise field propagation in the fiber. In the case of linear perturbation, it is independent of the specific noise realizations, which makes it possible to calculate the moment-generating function (MGF) of the filtered photoelectric current at the receiver [1, 68].

MGF method is an approach making use of MGF to evaluate the noise impact on the system performance. Now, it is well known that this method is accurate for various linear OFC systems. For nonlinear OFC systems, the MGF method can also be applied, provided their noise fields obey LNE. One can see this from Doob’s theorem, which means that, in a linearizable system driven by Gaussian-distributed noise, each of the independent random variable keeps Gaussian (cf. P. 35 of Ref. [9]). Because of this, MGF of the received photoelectric current can be calculated like those in linear OFC systems.

Common form of LNE [28] is based on the continuous wave (CW) assumption, i.e., the noise-free signal in the LNE is artificially simplified as a CW wave. As a result, a semi-analytical form of noise propagator and the noise power spectral density (PSD), so called parametric gain (PG), can be obtained. The drawback of this simplification is that the noise-free signal in this LNE neglects chromatic dispersion (CD) effect. As a result, the couplings between noise components (in frequency domain) cannot be taken into account.

A LNE beyond CW, named accurate LNE in this work, was first proposed and discussed in Ref. [1]. Dynamically taking into account the local CD and Kerr nonlinearity along the fiber, this LNE provides accurate noise information, with its computational cost being much higher than the CW approach. For example, given a noise-free signal obtained from nonlinear Schrödinger equation (NLSE), the computation required to update the accurate LNE has cubic complexity in the number of Fourier components [1]. To reduce the computational complexity, covariance matrix method (CMM) was proposed in Ref. [1], where the noise covariance matrix was obtained by processing large noise realizations. In Refs. [10,11], the computational cost of CMM was further reduced by a deterministic approach using perturbation solution. Since the raw covariance matrix obtained from NLSE via Monte Carlo noise realizations [1] or perturbation solution [10, 11] may contain nonlinear noise contribution, it needs to be separated from the nonlinearity-induced phase and timing jitter. Thus, the obtained pdfs of the receiver voltage agrees well with Monte Carlo simulation [1, 10, 12].

With the help of linear perturbation, the noise covariance matrix can also be obtained from its ordinary differential equation (ODE) proposed by [1,13]. The covariance matrix obtained by solving such linear ODE does not need jitter separation, although this ODE is more complicated than the accurate LNE [13]. So far there is little comparison between the approaches of Ref. [13] and Refs. [1, 10, 11].

In this work, we simplify the CMM by showing that the noise propagator matrix can be obtained directly from the accurate LNE. Therefore, there is no nonlinearity-induced jitter. To effectively reduce the computational complexity in updating the LNE, one can decompose the Kerr effect related matrix into a symmetric and an antisymmetric matrices [cf. the discussion after Eq. (30) in Appendix A]. Making use of the fourth-order Runge-Kutta in the interaction picture (RK4IP) method [14, 15], the accurate LNE can be solved with large step size, as detailed in Sec. 2. We evaluate the impacts of noise propagator on moment-generating function (MGF) and BER in Sec. 3. The accuracy of this new approach depends on how far the linearized noise deviates from the actual noise. To numerically verify this new approach, we consider the BERs in a 20-span DPSK system with nonlinear phase of Φ̄N = 0.2π [5] in Sec. 5. Our BER calculations agree well with the published CMM results. In Sec. 5, we also simulate the experiments of the multi-span DPSK systems discussed in Ref. [7] and show that, to fit the experimental data, one needs to take into account the nonlinearity induced phase difference between noise and noise-free signal, which will affect the signal-noise beating significantly.

2. Noise propagator obtained from accurate LNE

In an OFC system amplified by EDFAs, the noise propagator is a fundamental matrix that determines the noise impacts on MGF and BER. In this section, we show that the noise propagator matrix in a fiber of length L can be obtained from the accurate LNE [1]. For a multi-span OFC system, one needs to introduce an equivalent noise propagator which can be obtained from PG.

2.1. Noise propagator in a fiber of length L

The noise propagator in a fiber of length L can be obtained by extending the RK4IP in Refs. [14, 15] to the accurate LNE [1]. [In Appendix A, this LNE is rewritten as Eq. (32), where the linear operator is associated with CD effect, whereas the nonlinear operator is caused by Kerr nonlinearity.] By introducing ã = e(zz0)ãI and I = e(zz0)N̂e(zz0), the accurate LNE in the interaction picture (IP) has the form

da˜Idz=N^Ia˜I

Taking z0 = zn + h/2 with step size h = zn+1zn and denoting ãn = ã(zn), ãn+1 = ã(zn+1), N^nI=eN^h/2N^(zn)eN^h/2, N^n+1/2I=N^(zn+h/2), and N^n+1I=eL^h/2N^(zn+1)eL^h/2, one can use RK4IP [14, 15] to to solve Eq. (1) with

a˜n+1=eL^h/2[a˜nI+hk16+hk23+hk33+hk46]a˜nI=eL^h/2a˜nk1=N^nIa˜nI=eL^h/2N^(zn)a˜nk^1a˜nk2=N^n+1/2I[a˜nI+hk12]=N^(zn+h/2)[eL^h/2+hk^12]a˜nk^2a˜nk3=N^n+1/2I[a˜nI+hk22]=N^(zn+h/2)[eL^h/2+hk^22]a˜nk^3a˜nk4=N^n+1I[a˜nI+hk3]=eL^h/2N^(zn+1)eL^h/2[eL^h/2+hk^3]a˜nk^4a˜n
or
a˜n+1=(eL^h/2[eL^h/2+hk^16+hk^23+hk^33]+h6N^(zn+1)eL^h/2(eL^h/2+hk^3))a˜n,
which means the noise propagator for the fiber of length h = zn+1zn can be calculated as
H(zn+1,zn)=eL^h/2[eL^h/2+13(h2k^1+2h2k^2+hk^3)]+h6N^(zn+1)eL^h/2[eL^h/2+hk^3]
For the fiber of length L, the noise propagator has the form
pn(L,0)=H(L,LhL)H(h1,0).

Note that the RK4IP used here is different from the RK4IP in Ref. [15], where what to be solved was the noise-free signal (a 1D matrix), while here what we want is the noise propagator (a 2D matrix). The computational complexity for this 2D matrix is O(Nw3), due to that each i (i = 2, 3, 4) in Eq. (2) needs one dense matrix multiplication. Here Nw is the number of Fourier components used for signal representation.

2.2. Equivalent noise propagator of a multi-span system

The ASE from an EDFA can be modeled as additive white Gaussian noise (AWGN). Given the AWGN injected at the input of a fiber and the noise propagator obtained from Eqs. (2), (4), and (5), the noise PSD (or PG) at the output of a fiber of length L can be written as [1, 5, 7, 8]

PG1=pn(L,0)σ2IpnT(L,0)=σ2pn(L,0)pnT(L,0),
where I is a unit matrix and pnT is the transpose of pn. In Eq. (6), with G being the EDFA gain shown in Fig. 1, the variance of the real or imaginary part of input ASE can be expressed as
σ2=N0/(2T0)[N0=nsp(G1)h¯ω]
[cf. Ref. [16] or the discussions near Eqs. (22) and (33) in Appendix A].

 

Fig. 1 Low-pass equivalent optical model. The receiver consists of optical and electrical filters and DPSK balance detection. In this work, calculations, except those in Sec. 5.1, are based on the assumptions N0k = N0 and Gk = G [k = 1, 2,...(K + 1)]. Also, the fiber in each span is assumed to have same length (L km) and same loss (α dB/km). Nin is the ASE noise added at the transmitter. Changing Nin will change the OSNR at the receiver. In a typical balanced DPSK receiver, the delay in one branch of the interferometer is Tb = 1/Rb. Here, Rb = 20 Gb/s.

Download Full Size | PPT Slide | PDF

For a K-span system consisting of (K + 1) EDFAs, as shown in Fig. 1, its PG has the form

PG^=(Gσin2+σ2)Pn(K,0)PnT(K,0)+σ2k=1KPn(K,k)PnT(K,k),Pn(K,k)=pn(KL,(K1)L)pn((k+1)L,kL),(k=0,1,,K1)Pn(K,K)=I,
where σin2=Nin/(2T0). In Eq. (8), the real symmetric matrix PG^ is positive definite. It can be factorized as
PG^=σ2Pn,eqPn,eqT,
where the equivalent noise propagator Pn,eq can be obtained either by using Cholesky decomposition or symmetric (square root) decomposition [8, 9]. The latter yields Pn,eq=Pn,eqT.

3. MGF calculation

With the noise propagator matrix, one can evaluate the BER in the OFC system by calculating the MGF of the electrically filtered current I(ts) expressed using Karhunen-Loève series expansion (KLSE).

Due to the noise in the OFC system, the received (or filtered photoelectric) current fluctuates around its expectation value. The MGF of such current is a useful form to show the noise probability distribution. To get a simple form of MGF, the received current needs to be expressed using KLSE. For a nonlinear OFC system with its noise being linearizable, the KLSE form of the received current can be formulated. Averaging the “canonical” noise |i (i = 1, ···, 4Mn + 2) with formula [16, 17]

E[exp(s(λ˜c2+2cb˜))]=dc2πσ2exp(c22σ2)exp[s(λ˜c2+2cb˜)]=exp[2σ2s2b˜212σ2sλ˜]12σ2sλ˜,
the MGF of the filtered current, detailed near the end of Apendix B, can be written as
Ψts(s)=E[exp[sI(ts)]]=exp[sIss(ts)]i=14Mn+2exp[2σ2s2b˜i2(ts)1sβi](1sβi)ξ,(βi=2σ2λ˜i)
where I(ts) is the filtered photoelectric current at time ts. It consists of signal-signal beating (yss), noise-noise beating (ynn), and signal-noise beating (yns). In Eq. (11), i(ts) is the ith component of |(ts)〉 detailed at the end of Appendix B, while λ̃i is the power of ith component of the noise in Karhunen-Loève presentation. In this work, we take ξ = 1/2 for polarized noise.

With the help of Eq. (11) as well as Eqs. (7) and (8) in Ref. [17], the BER can be calculated.

4. OSNR at the receiver

For an optical system with ASE power being much larger than other noise sources, the OSNR with reference bandwidth Br (0.1nm) can be calculated as

OSNR0.1nm=P¯sPASE(Br).
In Eq. (12), s is the time-averaged (noise free) signal power, while PASE (Br) is the noise power within Br. To obtain s and PASE (Br), one needs to notice that the measurement bandwidth Bm [e.g., the bandwidth of the transfer function of an optical spectrum analyzer (OSA)] may not be the same as Br. Thus the s in Eq. (12) becomes the power of the signal filtered by Bm, while PASE (Br) becomes the ASE filtered by Bm and weighted by a factor Br/Bm [18].

In a linear optical system, the ASE noise along the fiber can be treated as AWGN. Thus, for the system of Fig. 1, its OSNR can be simply calculated as

OSNRL,0.1nm=P¯PASE(Bm)×BmBr(PASE(Bm)=[GNin+N0(K+1)]Bm),
where PASE (Bm) is the ASE power within Bm and N0 is given by Eq. (22). In Eq. (13), the filter (Bm) effect on the ASE has been neglected.

In a nonlinear optical system, the ASE noise “amplified by” PG cannot be treated as a white noise. Similar to the noise-noise beating calculation detailed at the end of Appendix B, the measured ASE power, which is only relate with the self-beating terms of the noise, has the form

OSNRNL,0.1nm=P¯sPASE(Bm)×BmBr(PASE(Bm)=Tr(OmTPG^Om)σ2=Tr(PG^OmOmT)σ2).
Here Om is the low-pass transfer function of the bandpass filter (bandwidth Bm). In Eq. (14), PG^ is given by Eq. (8).

In the case of traditional OSA-based out-of-band OSNR monitoring, the ASE power can be interpolated using

PASE(Bm,±Δλ)=Tr[PG^(Om(Δλ)OmT(Δλ)+Om(+Δλ)OmT(+Δλ))]σ22,
where Om(±Δλ) is the filter function centered at ±Δλ. When Δλ = 0, Eq. (15) returns to the ASE power in Eq. (14), where PASE (Bm) = PASE (Bm, 0).

5. Applications to DPSK systems

To show that the new approach to get the noise propagator is numerically applicable, we will compare our RK4IP results with the CMM results given by Ref. [5] and with the experimental data given by Ref. [7]. Both consider systems with Rb =20 Gb/s, using RZ-50% DPSK modulation. In the receiver, the optical filter is Gaussian type, while the electric filter is the fifth-order Bessel type. In the following calculations, we set T0 = NTb by changing μ in Eq. (36). This means, given the noise propagator matrix, the computational cost for BER is much higher than that in the linear case. For BER calculations, the length of the de Bruijn sequence is N = 25 [16]. Based on the relation between the RK4IP step and the fiber dispersion length (LD) or nonlinear length (LN) discussed in Ref. [15] as well as the detailed fiber parameters in the following discussion, we let the RK4IP step for the transmission fiber (htr) and that for the DCF fiber (hDCF) be related with htr : hDCF = (5 ∼ 6) : 1. The value of htr and the required computational time will be detailed below.

5.1. Comparison with CMM results

We consider the 20-span DPSK system discussed in Ref. [5], where BERs using the CMM and the (improved) CW approaches were plotted against the received OSNR in the Fig. 8 of Ref. [5]. In fact this system is basically the same as the one shown in Fig. 1, provided that one removes the pre- and postcompensating fibers and their amplifiers in the Fig. 2 of Ref. [5] and removes the first EDFA and Nin in our Fig. 1. Thus the first term of PG^ [in Eq. (8)] needs to be ignored. Like Refs. [5, 19], where OSNR was calculated in the absence of PG, we obtain OSNR from Eq. (13) with Nin = 0 and (K + 1) being replaced by K. According to Ref. [5], we change OSNR by changing the nsp in Eq. (22). As plotted in Fig. 1, each span contains a transmission fiber followed by a dispersion-compensating fiber (DCF). The transmission fiber is l =100 km long with its CD parameter Dtx = 8 ps/nm/km. Each span is fully compensated. The nonlinear phase accumulated in the fiber, defined as Φ¯NL=0zγ(ξ)Pineα(ξ)ξdξ with Pin being the time averaged signal power (at the input of the fiber), is 0.2π. The bandwidth of the optical (electrical) filter in the receiver is Bo = 1.8Rb (Be = 0.65Rb), respectively.

 

Fig. 2 BER versus received OSNR for a 20Gb/s 20-span RZ-DPSK system with Φ̄NL = 0.2π. Solid: obtained using CMM of Ref. [1]. Dashed (dotted): improved CW approach of Ref. [5] with CW power being peak power (average power), respectively. Dash-dotted: RK4IP approach. All curves, except the dash-dotted, are obtained from Fig. 8 of Ref. [5].

Download Full Size | PPT Slide | PDF

To let our results be reproducible, we provide, as detailed as possible, other related parameters below. The DCF in each span is 8 km long with DDCF = −100 ps/nm/km. Transmission fiber and DCF are assumed to have same fiber loss (α =0.2 dB/km) and same nonlinear coefficient (γ =2.0 /W/km). The EDFA in each span is used to compensate the total loss in the fiber of L = (100 + 8) km. Therefore the signal power at the input of each span (Pin) keeps constant. Ignoring the nonlinear phase contribution of DCF [20], we set Pin=0.7307 mW, obtained from Φ̄NL = KγPin(1 − eαl)/α = 0.2π with K = 20 and l = 100 km. The OSNR is obtained using Eq. (13) with Bm/Br = 1.35.

As shown in Fig. 2, the curve using the proposed RK4IP approach (dash-dotted) agrees very well with the CMM curve (solid) given by Ref. [5]. In Fig. 2, the curves using improved CW approach [5] with CW power being transmitted peak power (dashed) and average power (dotted) are plotted for comparison.

For each calculated point in Fig. 2, the CPU time for the noise propagator calculation is ∼1.8 hr with RK4IP step for transmission fiber (20×100 km long) being htr =3.5 km. The CPU time for each BER calculation is ∼ 0.5 hr. In fact, htr ranged within 0.3 ∼ 5 km yields almost the same curve.

5.2. Comparison with experimental data

The optical system discussed in Ref. [7] can be modelled by Fig. 1, except that each EDFA should be replaced by an EDFA followed by an optical filter (Gaussian) Olk with bandwidth Blk =5 nm. Also, the input noise Nin needs to be filtered by an optical filter Oin (Gaussian, Bin =3 nm). As a result, in Eq. (8), the noise propagator pn ((k + 1)L, kL) (k = 0,···, K − 1) should be replaced with Olk pn ((k + 1)L, kL), Pn(K, K) = 1 with Pn(K, K) = Olk, and Gσin2Pn(K,0)PnT(K,0) with Gσin2Pn(K,0)OinOinTPnT(K,0). Since we only consider the curves plotted in Fig. 7(b) of Ref. [7], each fiber (L km long) in Fig. 1 contains a SMF (42 km long) followed by a DCF (7 km long). In our calculation, all the related fiber parameters are same as those given in Table 1 of Ref. [7]. In the receiver, the bandwidth of the optical filter is 1.87Rb, while the bandwidth of the electrical filter is 0.75Rb.

We first consider the back-to-back case. Similar to Ref. [7], we modify the 20Gb/s RZ-50% signal at the transmitter by comparing its calculated spectrum [21] and its measured spectrum [7]. Due to that the input noise Nin is filtered by Oin, the OSNR is calculated using Eq. (14) with Bm/Br = 0.95, yielding the back-to-back RK4IP curve shown in Fig. 3.

 

Fig. 3 BER versus received OSNR for the multi-span RZ-50% DPSK systems with Φ̄NL = 0.9. Each span consists of a SMF fiber (42 km long) and a DCF file (7 km long). Other fiber parameters were detailed in Table 1 of Ref. [7]. According to Fig. 7(b) in Ref. [7], where the experimental curves for 5-span, 10-span, and 25-span systems were almost the same, here we replot these three curves using a thick solid curve.

Download Full Size | PPT Slide | PDF

For 5-span, 10-span, and 25-span systems, their accumulated nonlinear phases are calculated according to Eq. (48) of Ref. [7]. Because of the spectral modification of the input signal, the optical power at the input of each span Pin = PSMF is smaller than Eb/Tb, where Eb is the energy per bit before the spectral modification. For example, to get nonlinear phase Φ̄N = 0.9 for the 25-span system, the fiber input power Pin = PSMF should be 1.516 mW, which means Eb/Tb =0.127 mW or G(Eb/Tb) = 2.316 mW (G = 18.197). Different from the DPSK receiver shown in Fig. 1, where the delay is Tb = 1/Rb = 50 ps, the delay in the receiver of Ref. [7] was T′b=(24.84 GHz)−1=40.26 ps. Thus, the DPSK phase factors given in Eq. (35) should be modified as

Dllss=ej2πlN+ej2πlN2,Dmmnn=ej2πmTbT0+ej2πmTbT02,Dmlns=ej2πmTbT0jΔ+ej2πlN+jΔ2,
with N′ = N(Tb/T′b), T′b = Tb + ΔTb. In Eq. (16), Δ is introduced as
Δ=ϕ0<δϕ>,
where <δϕ>, given by Eq. (41) in Appendix C, is the nonlinear phase difference between noise and noise-free signal. As shown in Fig. 3, all RK4IP curves (ϕ0 = 0.31) agree very well with the experiment results. The ASE power is calculated using Eq. (15) with Δλ = 2Bm. In Eq. (17), ϕ0 is a calibration constant that basically shifts the RK4IP curves in the OSNR direction, while <δϕ> determines the slope of the RK4IP curves. To show this, we plot in Fig. 4 the RK4IP results for the 25-span system with Δ = 0.31− <δϕ> and Δ = 0. Also, we consider the RK4IP curves using Eq. (14) to calculate ASE power. Our results for the 5-span, 10-span, and 25-span systems confirm that there is almost no difference between the curve using Eq. (15) with ϕ0 = 0.31 and the curve using Eq. (14) with ϕ0 = 0.57.

 

Fig. 4 BER vs RX-OSNR for the 5-span RZ-50% DPSK system with Φ̄NL = 0.9. Other parameters are same as those used in Fig. 3. Solid: experimental results. Dotted (Dash-dotted): numerical calculation using RK4IP with (without) phase shift Δ.

Download Full Size | PPT Slide | PDF

At present, the calibration constant ϕ0 in Eq. (17) can be treated as a fitting parameter. As mentioned above, it basically shifts the BER vs OSNR curve in the OSNR direction and is related with the detailed OSNR monitoring technique used in Ref. [7]. As there is few information about its OSNR measurement, evaluation of ϕ0 is expected to be discussed elsewhere.

In the above calculation, our CPU time for the noise propagator calculation is ∼0.8 hr with RK4IP step for transmission fibers (25×40 km long) being htr =6.0 km. The CPU time for each BER calculation is ∼ 0.5 hr. The step size htr within 0.3 ∼ 10 km will result in almost the same curve.

6. Summary

For linear OFC systems, MGF method is useful for one to evaluate the noise impacts on BERs. This is true not only because it is computationally efficient for the cases with low BER (e.g., < 10−9) but also it can provide reliable information for the cases using coherent detection, which is now widely used in modern OFC systems. It is now well recognized that traditional Gaussian fitting Q-factor approximation is accurate for OOK detection, while MGF method is accurate for various linear OFC systems.

To extend the MGF method to a nonlinear OFC system, one needs to make sure its noise propagator varies within linear regime. This means the noise-noise interaction needs to be neglected, so that the noise field can follow LNE.

In the approaches of CMM [1, 11, 12], the noise propagator was obtained from NLSE. It may contain the nonlinearity-induced jitter, which is beyond the linear regime and should be removed. In this work we simplify the CMM by directly solving the accurate LNE [Eq. (1)] with RK4IP. Like the CW approximations (cf. Refs. [28] and many others) as well as the approach of Ref. [13], where the noise propagation information obtained from the related linearized noise equation is automatically free from nonlinearity-induced jitter, our noise propagator obtained from accurate LNE also varies within linear regime.

To numerically verify our new approach, we consider a 20-span RZ-DPSK system discussed in Ref. [5]. The BERs obtained using this new RK4IP agree well with those using CMM in Ref. [5]. Taking account the phase difference between noise and noise-free signal leads to quantitative matching between numerical evaluation and experimental result [7].

Appendix A: Accurate LNE in the EDFA-based systems

The optical field u(z, t) in a fiber satisfies

uz=jβωω22ut2+βωωω63ut3jγ|u|2uα2u,
where α is the fiber loss and βωω = 2β/∂ω2 relates to the CD parameter D(λ) (ps/nm/km) with βωω=λ2D(λ)2πc (c=3 × 108 m/s). The slope parameter βωωω=(λ2πc)2[2λD(λ)+λ2D(λ)] [ D(λ)=dD(λ)dλ (ps/nm2·km)] can be neglected if bit rate Rb satisfies Rb > |βωωωωω| [16].

Introducing the transformation u(z, t) = v(z, t)eαz/2, Eq. (18) can be reduced as

vz=jβωω22vt2+βωωω63vt3jeαzγ|v|2v.
In a K-span system amplified by (K + 1) EDFAs (cf. Fig. 1), Eq. (19) can be modified as
vz=jβωω22vt2+βωωω63vt3jeαzγ|v|2vjw(z,t),
where w(z, t) is the ASE forcing modeled as the complex AWGN with correlation
r(z,z,t,t)=E{w(z,t)w*(z,t)}=δ(tt)δ(zz)k=1K+1N0kδ(z(k1)L).
In Eq. (21), the fiber length in each span is assumed to be L (km) long, According to Wiener-Khintchine theorem [22], N0k in Eq. (21) is the ASE PSD (in one polarization direction) at the output of the kth EDFA. Suppose each EDFA has the same gain G and spontaneous-emission parameter nsp, we have [16]
N0k=N0=nsp(G1)h¯ω(k=1,2,,K,(K+1))

Decomposing optical field v(z, t) in the fiber into noise-free field v0(z, t) and its perturbation δv(z, t) [i.e., v(z, t) = v0(z, t) + δv(z, t)] and assuming that |v0| >> |δv| (so that the nonlinear terms of δv can be neglected), Eq. (20) can be decomposed as [1]

v0z=jβωω22v0t2+βωωω63v0t3jeαzγ|v0|2v0
δvz=jβωω22δvt2+βωωω63δvt3j2eαzγ|v0|2δvjeαzγv02δv*jw(z,t).
Noise equation Eq. (24) differs from common equations in that the real and imaginary parts of the complex noise field need to be treated separately [1].

Denoting al = a(ωl) = ∫ δvelt dt and the circulant matrices [Mν]lm = νlm, [Mμ]lm = μl+m with

νl=ν(ωl)=eαz|v0|2ejωltdt,μl=μ(ωl)=eαzv02ejωltdt,
in frequency domain, Eq. (24) has the form
daldz=jβωω2ωl2aljβωωω6ωl3alj2γ(z)[Mν]lmamjγ(z)[Mμ]lmam*jWl,
where Wl = W (z, ωl) is the Fourier component of the forcing term w(z, t) in Eq. (20). As indicated in Ref. [1], since |v0|2 in (25) is real, νl=νl*. So Mν in (26) is Hermitian, or, its real part MνR is symmetric, while its imaginary part MνI is anti-symmetric. Also, as [Mμ]km = μk+m, both the real ( MμR) and imaginary parts ( MμI) of Mμ are symmetric.

The matrix form of Eq. (26) is

dadz=L¯a+νa+μa*jW
L¯=jLCD,ν=2jγ(z)(MνR+jMνI),μ=jγ(z)(MμR+jMμI)
Introducing ã = (aR, aI)T (for a = aR + jaI) and = (WI,WR)T (for −jW = WIjWR), Eq. (27) is equivalent to
da˜dz=(L^+ν^+μ^)a˜+W˜
L^=(0LCDLCD0),ν^=(νAAνSSνSSνAA),μ^=(μRμIμIμR),
with (LCD)ij=[βωω2ω2+βωωω6ω3]δij, νAA=2γMνI, νSS=2γMνR, μR=γMμI, and μI=γMμR. According to the discussion given below Eq. (26), ν̂ (μ̂) in Eq. (30) is antisymmetric (symmetric), respectively. Calculation of the Kerr term (ν̂ + μ̂) according to Eq. (30) has the computational complexity much less than O(NW3), where NW is the number of Fourier components used for signal representation. In fact, the computational cost of this way is basically determined by the FFTs in Eq. (25), which has the computational complexity of O(NW logNW).

In frequency domain, ASE correlation relation (21) has its matrix form ( W˜(ωl)W˜lΔf

E{W˜l(z)W˜l*(z)}=δz,zδl,lk=0KN02T0δz,kL(l=1,,4Mn+2),
where Eq. (22) has been used. In Eq. (31),T0 = 1/Δf and Mn are given by Eq. (36). Eq. (31) means that Eq. (29) can be equivalently replaced by
da˜dz=(L^+N^)a˜;(N^=ν^+μ^)
with boundary condition [23]
E{W˜(zf)W˜*(zf)}=N02T0I=σ2I,
where I is a (4Mn + 2) × (4Mn + 2) unit matrix.

Appendix B: Filtered photoelectric current expressed using KLSE

Given a linear optical system, based on the discussions in Ref. [16] and the notations introduced in Ref. [17], the filtered photoelectric current I(t) can be expressed in the form of I(t) = [〈so(t + Tb) + no(t + Tb)|so(t) + no(t)〉 + c.c.]/2. Here so(t) (no(t)) represents the signal (noise) field at the input of the optical filter. Dirac bra 〈x| is the conjugate transpose (or Hermitian transpose) of Dirac ket |x〉 [x = so(t), no(t), so(t) + no(t), etc.]. The Dirac ket differs from usual complex vector in that the ith element of the latter is just the ith Fourier coefficient of the (signal or noise) field, while the ith element of the former is the product of the ith Fourier coefficient and its base function (cf. Eq. (17) in Ref. [17]). According to Refs. [16, 17], the filtered current in a linear system can be formally expressed as I(t) = yss + ynn + yns with (l = −Ls, ··· Ls; m = −Mn ··· Mn)

yss(ts)=[so(ts+Tb)|Rss|so(ts)+c.c.]/2=so(ts)|RssD|so(ts)ynn(ts)=[no(ts+Tb)|Rnn|no(ts)+c.c.]/2=No|RnnD|No=Z|Λ|Zyns(ts)=[no(ts+Tb)|Rns|so(ts)+no(ts)|Rns|so(ts+Tb)+c.c.]/2=[Nin|OnnRnsD|so(ts)+c.c.]=[Z|bD(ts)+c.c.]
where ΛUOnnRnnDOnnU=diag{λ1,,λ2Mn+1}, |bD(ts)=UOnnRnsD|so(ts), and
Dllss=ej2πlN+ej2πlN2,Dmmnn=ej2πmTbT0+ej2πmTbT02,Dmlns=ej2πmTbT0+ej2πlN2.
In Eqs. (34)(35), |Z〉 represents the decoupled Gaussian random variables with zero mean and real part and imaginary part variance of σ2. The effects of the optical and electrical filters in the receiver are represented by matrices with their elements being (Onn)mm=δm,mHo(mT0), (Oss)ll=δl,lHo(lNTb) and (Rss)ll=Hr(llNTb), (Rnn)mm=Hr(mmT0), (Rns)ml=Hr(lNTbmT0). Due to the optical and electrical filters, signal (noise) components outside ±LsMn) can be neglected. Here [16]
Ls=ηNTbBo,Mn=ηBoT0,T0=μ(1Bo+1Br).

For a nonlinear optical system, to get the noise propagator from the accurate LNE (29), one needs to separate complex numbers into their real and imaginary parts. Denoting the Re-Im form of a complex matrix x as x˜=(Re{x}Im{x}Im{x}Re{x}) [where x can be any complex matrix in Eq. (34)] and introducing |no〉 = Pn,eq|a0〉 with |a0〉 being the AWGN from EDFA,

|s˜o=[Re{|so}Im{|so}],|a˜0=[Re{|a0}Im{|a0}],|Z˜=[Re{|Z}Im{|Z}],
it is easy to generalize the noise related currents in Eq. (34) as [1, 5, 7, 8, 16]
ynn(ts)=a˜0|U˜Λ˜U˜T|a˜0=Z˜|Λ˜|Z˜(Λ˜=U˜TPn,eqTO˜nmTPn,eqU˜=diag{λ˜1,,λ˜4Mn+2})yns(ts)=Z˜|U˜TP˜n,eqTO˜nmTR˜nsDO˜ss|s˜o(ts)Z˜|B˜|s˜o(ts)Z˜|b˜(ts)(B˜=U˜TPn,eqTO˜nnTR˜nsDO˜ss),
where Pn,eq can be obtained from Eqs. (2)(5) and (8)(9).

Appendix C: Nonlinearity induced phase difference between noise and noise-free signal

I: The phase difference caused by Nin

In this part, we assume that, in Fig. 1, the external noise injected at the transmitter (Nin) is much larger than the ASE noise from the EDFAs, which is true for the experiments discussed in Ref. [7]. Thus, one can only consider the phase difference caused by Nin and ignores the effect of N0k (k = 1, ···, K).

It is well known that, for the noise-free signal with its path average power being , its nonlinear phase accumulated at the fiber output is Φ̄NL = P̄γKL, where L is the fiber length of each span, as denoted in Fig. 1. Due to the optical power fluctuation δP̄, the actual nonlinear phase becomes

ϕNL=(P¯+δP¯)γKL.
Relative to the noise-free signal, the noise-induced phase change, ϕNL − Φ̄NL, varies randomly. The average variance of such phase noise can be calculated as
<δϕ2>=<ϕNL2Φ¯NL2>2P¯<δP¯>(γKL)2.
For the experiments in Ref. [7], Nin is filtered with bandwidth of Bin=3nm. Thus, we have <δP̄>= GNinBin. Eq. (40) yields
<δϕ>2P¯γKL/P¯/<δP¯>=2Φ¯NL/OSNR,
where OSNRP̄/(GNinBin) is the input OSNR.

Considering that δP̄ ≥ 0, we have ϕNL − Φ̄NL = δP̄γKL > 0. This means ϕNL rotates faster than Φ̄NL and there is a phase difference between the actual optical field and the noise-free signal. As part of the actual field, the noise field also has the same phase shift relative to the noise-free signal. Note that this phase shift will not affect signal-signal and noise-noise beatings. But it will affect the signal-noise beating. In fact, when calculating the signal-noise beating, the noise and signal should be treated consistently. Or, they should be considered within the same coordinate system. Thus the phase shift between noise field and and noise-free signal should be taken into account. In this work, <δϕ> given by Eq. (41) is approximated as the average of such phase shift. Our numerical results plotted in Figs. 3 and 4 confirm the validity of this approximation.

For the experiments in Ref. [7], we have Φ̄NL = 0.9. In this case, Eq. (41) yields

<δϕ1/OSNRπ/2arctan(OSNR),
where arctan(x) + arctan(1/x) = π/2 and arctan(x) ≈ x (for x → 0) have been used. Obviously, <δϕ> in Eq. (42) relates ϕGM in Eq. (23) of Ref. [24] with <δϕ> +ϕGM = π/2. Also, ϕ0 in Ref. [24] now becomes ϕ0 + π/2 → ϕ0 in Eq. (17), while ΔGM in Ref. [24] is named as Δ in this work.

II: The phase difference caused by Nin and N0k (k = 1, ···, K)

For the K-span system of Fig. 1 with Nin ≠ 0 and K not being large enough and with the ASE from each EDFA being filtered by Olk and the ASE injected at the transmitter being filtered by Oin, Eq. (40) can be generalized as

<δϕ2>=2[GNinBin+N0BlkK3(1+1K)(1+12K)]Φ¯NL2/P¯,
where k=1Kk2=K33(1+1K)(1+12K) has been used. Obviously, in the case of GNinBin >> KN0Blk [7], Eq. (43) yields Eq. (41).

Acknowledgments

The authors acknowledge the financial support from Canada Research Chair program. The first author sincerely thanks Paolo Serena and Leonardo D. Coelho for providing their detailed calculations in Refs. [5, 7]. The authors wish to thank the anonymous reviewers for their valuable comments and suggestions.

References and links

1. R. Holzlöhner, V. S. Grigoryan, C. R. Menyuk, and W. L. Kath, “Accurate calculation of eye diagrams and bit error rates in optical transmission systems using linearization,” J. Lightwave Technol. 20, 389–400 (2002). [CrossRef]  

2. K. Kikuchi, “Enhancement of optical-amplifier noise by nonlinear refractive index and group-velocity dispersion of optical fibers,” IEEE Photon. Technol. Lett. 5, 1079–1081 (1993).

3. A. Carena, V. Curri, R. Gaudino, P. Poggiolini, and S. Benedetto, “New analytical results on fiber parametric gain and its effects on ASE noise,” IEEE Photon. Technol. Lett. 9, 535–537 (1997). [CrossRef]  

4. R. Hui, M. O’Sullivan, A. Robinson, and M. Taylor, “Modulation instability and its impact in multispan optical amplified IMDD systems: theory and experiments,” J. Lightwave Technol. 15, 1071–1082 (1997). [CrossRef]  

5. P. Serena, A. Orlandini, and A. Bononi, “Parametric-gain approach to the analysis of single-channel DPSK/DQPSK systems with nonlinear phase noise,” J. Lightwave Technol. 24, 2026–2037 (2006). [CrossRef]  

6. A. Demir, “Nonlinear phase noise in optical-fiber-communication systems,” J. Lightwave Technol. 25, 2002–2032 (2007). [CrossRef]  

7. L. D. Coelho, L. Molle, D. Gross, N. Hanik, R. Freund, C. Caspar, E.-D. Schmidt, and B. Spinnler, “Modeling nonlinear phase noise in differentially phase-modulated optical communication systems,” Opt. Express 17, 3226–3241 (2009). [CrossRef]   [PubMed]  

8. M. Secondini, E. Forestieri, and C. R. Menyuk, “A combined regular-logarithmic perturbation method for signal-noise interaction in amplified optical systems,” J. Lightwave Technol. 27, 3358–3369 (2009). [CrossRef]  

9. R. Holzlöhner, “A covariance matrix method for the computation of bit errors in optical transmission systems,” Ph.D. thesis, University of Maryland Baltimore County. Baltimore, Maryland, USA (2003).

10. R. Holzlöhner, C. R. Menyuk, and W. L. Kath, “Efficient and accurate computation of eye diagrams and bit error rates in a single-channel CRZ system,” IEEE Photon. Technol. Lett. 14, 1079–1081 (2002). [CrossRef]  

11. R. Holzlöhner, C. R. Menyuk, and W. L. Kath, “A covariance matrix method to compute bit error rates in a highly nonlinear dispersion-managed soliton system,” IEEE Photon. Technol. Lett. 15, 688–690 (2003). [CrossRef]  

12. R. Holzlöhner and C. R. Menyuk, “Use of multicanonical Monte Carlo simulations to obtain accurate bit error rates in optical communications systems,” Opt. Lett. 28, 1894–1896 (2003). [CrossRef]   [PubMed]  

13. A. Demir, “Non-Monte Carlo formulations and computational techniques for the stochastic nonlinear Schrodinger equation,” J. Comput. Phys. 201, 148–171 (2004). [CrossRef]  

14. J. Hult, “A fourth-order Runge-Kutta in the interaction picture method for simulating supercontinuum generation in optical fibers,” J. Lightwave Technol. 25, 3770–3775 (2007). [CrossRef]  

15. Z. Zhang, L. Chen, and X. Bao, “A fourth-order Runge-Kutta in the interaction picture method for numerically solving the coupled nonlinear Schrödinger equation,” Opt. Express 18, 8261–8276 (2010). [CrossRef]   [PubMed]  

16. E. Forestieri, “Evaluating the error probability in lightwave systems with chromatic dispersion, arbitrary pulse shape and pre- and postdetection filtering,” J. Lightwave Technol. 18, 1493–1503 (2000). [CrossRef]  

17. Z. Zhang, L. Chen, and X. Bao, “Accurate BER evaluation for lumped DPSK and OOK systems with PMD and PDL,” Opt. Express 15, 9418–9433 (2007). [CrossRef]   [PubMed]  

18. D. Gariepy and G. He, “Measuring OSNR in WDM systemsEffects of resolution bandwidth and optical rejection ratio,” White paper, EXFO Inc. (2009).

19. P. Serena, A. Bononi, J. C. Antona, and S. Bigo, “Parametric gain in the strongly nonlinear regime and its impact on 10-Gb/s NRZ systems with forward-error correction,” J. Lightwave Technol. 23, 2352–2363 (2006). [CrossRef]  

20. A. Bononi, P. Serena, and A. Orlandini, “A unified design framework for single-channel dispersion-managed terrestrial systems,” J. Lightwave Technol. 26, 3617–3631 (2008). [CrossRef]  

21. E. Ip and J. M. Kahn, “Power spectra of return-to-zero optical signals,” J. Lightwave Technol. 24, 1610–1618 (2006). [CrossRef]  

22. A. Papoulis, Probability, Random Variables, and Stochastic Processes (McGraw-Hill, 1991).

23. E. Forestieri and G. Prati, “Exact analytical evaluation of second-order PMD impact on the outage probability for a compensated system,” J. Lightwave Technol. 22, 988–996 (2004). [CrossRef]  

24. Z. Zhang, L. Chen, and X. Bao, “The noise propagator in an optical system using EDFAs and its effect on system performance: accurate evaluation based on linear perturbation,” arXiv:physics.optics/1207.3362v1.

References

  • View by:
  • |
  • |
  • |

  1. R. Holzlöhner, V. S. Grigoryan, C. R. Menyuk, and W. L. Kath, “Accurate calculation of eye diagrams and bit error rates in optical transmission systems using linearization,” J. Lightwave Technol. 20, 389–400 (2002).
    [Crossref]
  2. K. Kikuchi, “Enhancement of optical-amplifier noise by nonlinear refractive index and group-velocity dispersion of optical fibers,” IEEE Photon. Technol. Lett. 5, 1079–1081 (1993).
  3. A. Carena, V. Curri, R. Gaudino, P. Poggiolini, and S. Benedetto, “New analytical results on fiber parametric gain and its effects on ASE noise,” IEEE Photon. Technol. Lett. 9, 535–537 (1997).
    [Crossref]
  4. R. Hui, M. O’Sullivan, A. Robinson, and M. Taylor, “Modulation instability and its impact in multispan optical amplified IMDD systems: theory and experiments,” J. Lightwave Technol. 15, 1071–1082 (1997).
    [Crossref]
  5. P. Serena, A. Orlandini, and A. Bononi, “Parametric-gain approach to the analysis of single-channel DPSK/DQPSK systems with nonlinear phase noise,” J. Lightwave Technol. 24, 2026–2037 (2006).
    [Crossref]
  6. A. Demir, “Nonlinear phase noise in optical-fiber-communication systems,” J. Lightwave Technol. 25, 2002–2032 (2007).
    [Crossref]
  7. L. D. Coelho, L. Molle, D. Gross, N. Hanik, R. Freund, C. Caspar, E.-D. Schmidt, and B. Spinnler, “Modeling nonlinear phase noise in differentially phase-modulated optical communication systems,” Opt. Express 17, 3226–3241 (2009).
    [Crossref] [PubMed]
  8. M. Secondini, E. Forestieri, and C. R. Menyuk, “A combined regular-logarithmic perturbation method for signal-noise interaction in amplified optical systems,” J. Lightwave Technol. 27, 3358–3369 (2009).
    [Crossref]
  9. R. Holzlöhner, “A covariance matrix method for the computation of bit errors in optical transmission systems,” Ph.D. thesis, University of Maryland Baltimore County. Baltimore, Maryland, USA (2003).
  10. R. Holzlöhner, C. R. Menyuk, and W. L. Kath, “Efficient and accurate computation of eye diagrams and bit error rates in a single-channel CRZ system,” IEEE Photon. Technol. Lett. 14, 1079–1081 (2002).
    [Crossref]
  11. R. Holzlöhner, C. R. Menyuk, and W. L. Kath, “A covariance matrix method to compute bit error rates in a highly nonlinear dispersion-managed soliton system,” IEEE Photon. Technol. Lett. 15, 688–690 (2003).
    [Crossref]
  12. R. Holzlöhner and C. R. Menyuk, “Use of multicanonical Monte Carlo simulations to obtain accurate bit error rates in optical communications systems,” Opt. Lett. 28, 1894–1896 (2003).
    [Crossref] [PubMed]
  13. A. Demir, “Non-Monte Carlo formulations and computational techniques for the stochastic nonlinear Schrodinger equation,” J. Comput. Phys. 201, 148–171 (2004).
    [Crossref]
  14. J. Hult, “A fourth-order Runge-Kutta in the interaction picture method for simulating supercontinuum generation in optical fibers,” J. Lightwave Technol. 25, 3770–3775 (2007).
    [Crossref]
  15. Z. Zhang, L. Chen, and X. Bao, “A fourth-order Runge-Kutta in the interaction picture method for numerically solving the coupled nonlinear Schrödinger equation,” Opt. Express 18, 8261–8276 (2010).
    [Crossref] [PubMed]
  16. E. Forestieri, “Evaluating the error probability in lightwave systems with chromatic dispersion, arbitrary pulse shape and pre- and postdetection filtering,” J. Lightwave Technol. 18, 1493–1503 (2000).
    [Crossref]
  17. Z. Zhang, L. Chen, and X. Bao, “Accurate BER evaluation for lumped DPSK and OOK systems with PMD and PDL,” Opt. Express 15, 9418–9433 (2007).
    [Crossref] [PubMed]
  18. D. Gariepy and G. He, “Measuring OSNR in WDM systemsEffects of resolution bandwidth and optical rejection ratio,” White paper, EXFO Inc. (2009).
  19. P. Serena, A. Bononi, J. C. Antona, and S. Bigo, “Parametric gain in the strongly nonlinear regime and its impact on 10-Gb/s NRZ systems with forward-error correction,” J. Lightwave Technol. 23, 2352–2363 (2006).
    [Crossref]
  20. A. Bononi, P. Serena, and A. Orlandini, “A unified design framework for single-channel dispersion-managed terrestrial systems,” J. Lightwave Technol. 26, 3617–3631 (2008).
    [Crossref]
  21. E. Ip and J. M. Kahn, “Power spectra of return-to-zero optical signals,” J. Lightwave Technol. 24, 1610–1618 (2006).
    [Crossref]
  22. A. Papoulis, Probability, Random Variables, and Stochastic Processes (McGraw-Hill, 1991).
  23. E. Forestieri and G. Prati, “Exact analytical evaluation of second-order PMD impact on the outage probability for a compensated system,” J. Lightwave Technol. 22, 988–996 (2004).
    [Crossref]
  24. Z. Zhang, L. Chen, and X. Bao, “The noise propagator in an optical system using EDFAs and its effect on system performance: accurate evaluation based on linear perturbation,” arXiv:physics.optics/1207.3362v1.

2010 (1)

2009 (2)

2008 (1)

2007 (3)

2006 (3)

2004 (2)

E. Forestieri and G. Prati, “Exact analytical evaluation of second-order PMD impact on the outage probability for a compensated system,” J. Lightwave Technol. 22, 988–996 (2004).
[Crossref]

A. Demir, “Non-Monte Carlo formulations and computational techniques for the stochastic nonlinear Schrodinger equation,” J. Comput. Phys. 201, 148–171 (2004).
[Crossref]

2003 (2)

R. Holzlöhner, C. R. Menyuk, and W. L. Kath, “A covariance matrix method to compute bit error rates in a highly nonlinear dispersion-managed soliton system,” IEEE Photon. Technol. Lett. 15, 688–690 (2003).
[Crossref]

R. Holzlöhner and C. R. Menyuk, “Use of multicanonical Monte Carlo simulations to obtain accurate bit error rates in optical communications systems,” Opt. Lett. 28, 1894–1896 (2003).
[Crossref] [PubMed]

2002 (2)

R. Holzlöhner, C. R. Menyuk, and W. L. Kath, “Efficient and accurate computation of eye diagrams and bit error rates in a single-channel CRZ system,” IEEE Photon. Technol. Lett. 14, 1079–1081 (2002).
[Crossref]

R. Holzlöhner, V. S. Grigoryan, C. R. Menyuk, and W. L. Kath, “Accurate calculation of eye diagrams and bit error rates in optical transmission systems using linearization,” J. Lightwave Technol. 20, 389–400 (2002).
[Crossref]

2000 (1)

1997 (2)

A. Carena, V. Curri, R. Gaudino, P. Poggiolini, and S. Benedetto, “New analytical results on fiber parametric gain and its effects on ASE noise,” IEEE Photon. Technol. Lett. 9, 535–537 (1997).
[Crossref]

R. Hui, M. O’Sullivan, A. Robinson, and M. Taylor, “Modulation instability and its impact in multispan optical amplified IMDD systems: theory and experiments,” J. Lightwave Technol. 15, 1071–1082 (1997).
[Crossref]

1993 (1)

K. Kikuchi, “Enhancement of optical-amplifier noise by nonlinear refractive index and group-velocity dispersion of optical fibers,” IEEE Photon. Technol. Lett. 5, 1079–1081 (1993).

Antona, J. C.

Bao, X.

Benedetto, S.

A. Carena, V. Curri, R. Gaudino, P. Poggiolini, and S. Benedetto, “New analytical results on fiber parametric gain and its effects on ASE noise,” IEEE Photon. Technol. Lett. 9, 535–537 (1997).
[Crossref]

Bigo, S.

Bononi, A.

Carena, A.

A. Carena, V. Curri, R. Gaudino, P. Poggiolini, and S. Benedetto, “New analytical results on fiber parametric gain and its effects on ASE noise,” IEEE Photon. Technol. Lett. 9, 535–537 (1997).
[Crossref]

Caspar, C.

Chen, L.

Coelho, L. D.

Curri, V.

A. Carena, V. Curri, R. Gaudino, P. Poggiolini, and S. Benedetto, “New analytical results on fiber parametric gain and its effects on ASE noise,” IEEE Photon. Technol. Lett. 9, 535–537 (1997).
[Crossref]

Demir, A.

A. Demir, “Nonlinear phase noise in optical-fiber-communication systems,” J. Lightwave Technol. 25, 2002–2032 (2007).
[Crossref]

A. Demir, “Non-Monte Carlo formulations and computational techniques for the stochastic nonlinear Schrodinger equation,” J. Comput. Phys. 201, 148–171 (2004).
[Crossref]

Forestieri, E.

Freund, R.

Gariepy, D.

D. Gariepy and G. He, “Measuring OSNR in WDM systemsEffects of resolution bandwidth and optical rejection ratio,” White paper, EXFO Inc. (2009).

Gaudino, R.

A. Carena, V. Curri, R. Gaudino, P. Poggiolini, and S. Benedetto, “New analytical results on fiber parametric gain and its effects on ASE noise,” IEEE Photon. Technol. Lett. 9, 535–537 (1997).
[Crossref]

Grigoryan, V. S.

Gross, D.

Hanik, N.

He, G.

D. Gariepy and G. He, “Measuring OSNR in WDM systemsEffects of resolution bandwidth and optical rejection ratio,” White paper, EXFO Inc. (2009).

Holzlöhner, R.

R. Holzlöhner and C. R. Menyuk, “Use of multicanonical Monte Carlo simulations to obtain accurate bit error rates in optical communications systems,” Opt. Lett. 28, 1894–1896 (2003).
[Crossref] [PubMed]

R. Holzlöhner, C. R. Menyuk, and W. L. Kath, “A covariance matrix method to compute bit error rates in a highly nonlinear dispersion-managed soliton system,” IEEE Photon. Technol. Lett. 15, 688–690 (2003).
[Crossref]

R. Holzlöhner, C. R. Menyuk, and W. L. Kath, “Efficient and accurate computation of eye diagrams and bit error rates in a single-channel CRZ system,” IEEE Photon. Technol. Lett. 14, 1079–1081 (2002).
[Crossref]

R. Holzlöhner, V. S. Grigoryan, C. R. Menyuk, and W. L. Kath, “Accurate calculation of eye diagrams and bit error rates in optical transmission systems using linearization,” J. Lightwave Technol. 20, 389–400 (2002).
[Crossref]

R. Holzlöhner, “A covariance matrix method for the computation of bit errors in optical transmission systems,” Ph.D. thesis, University of Maryland Baltimore County. Baltimore, Maryland, USA (2003).

Hui, R.

R. Hui, M. O’Sullivan, A. Robinson, and M. Taylor, “Modulation instability and its impact in multispan optical amplified IMDD systems: theory and experiments,” J. Lightwave Technol. 15, 1071–1082 (1997).
[Crossref]

Hult, J.

Ip, E.

Kahn, J. M.

Kath, W. L.

R. Holzlöhner, C. R. Menyuk, and W. L. Kath, “A covariance matrix method to compute bit error rates in a highly nonlinear dispersion-managed soliton system,” IEEE Photon. Technol. Lett. 15, 688–690 (2003).
[Crossref]

R. Holzlöhner, C. R. Menyuk, and W. L. Kath, “Efficient and accurate computation of eye diagrams and bit error rates in a single-channel CRZ system,” IEEE Photon. Technol. Lett. 14, 1079–1081 (2002).
[Crossref]

R. Holzlöhner, V. S. Grigoryan, C. R. Menyuk, and W. L. Kath, “Accurate calculation of eye diagrams and bit error rates in optical transmission systems using linearization,” J. Lightwave Technol. 20, 389–400 (2002).
[Crossref]

Kikuchi, K.

K. Kikuchi, “Enhancement of optical-amplifier noise by nonlinear refractive index and group-velocity dispersion of optical fibers,” IEEE Photon. Technol. Lett. 5, 1079–1081 (1993).

Menyuk, C. R.

Molle, L.

O’Sullivan, M.

R. Hui, M. O’Sullivan, A. Robinson, and M. Taylor, “Modulation instability and its impact in multispan optical amplified IMDD systems: theory and experiments,” J. Lightwave Technol. 15, 1071–1082 (1997).
[Crossref]

Orlandini, A.

Papoulis, A.

A. Papoulis, Probability, Random Variables, and Stochastic Processes (McGraw-Hill, 1991).

Poggiolini, P.

A. Carena, V. Curri, R. Gaudino, P. Poggiolini, and S. Benedetto, “New analytical results on fiber parametric gain and its effects on ASE noise,” IEEE Photon. Technol. Lett. 9, 535–537 (1997).
[Crossref]

Prati, G.

Robinson, A.

R. Hui, M. O’Sullivan, A. Robinson, and M. Taylor, “Modulation instability and its impact in multispan optical amplified IMDD systems: theory and experiments,” J. Lightwave Technol. 15, 1071–1082 (1997).
[Crossref]

Schmidt, E.-D.

Secondini, M.

Serena, P.

Spinnler, B.

Taylor, M.

R. Hui, M. O’Sullivan, A. Robinson, and M. Taylor, “Modulation instability and its impact in multispan optical amplified IMDD systems: theory and experiments,” J. Lightwave Technol. 15, 1071–1082 (1997).
[Crossref]

Zhang, Z.

IEEE Photon. Technol. Lett. (4)

K. Kikuchi, “Enhancement of optical-amplifier noise by nonlinear refractive index and group-velocity dispersion of optical fibers,” IEEE Photon. Technol. Lett. 5, 1079–1081 (1993).

A. Carena, V. Curri, R. Gaudino, P. Poggiolini, and S. Benedetto, “New analytical results on fiber parametric gain and its effects on ASE noise,” IEEE Photon. Technol. Lett. 9, 535–537 (1997).
[Crossref]

R. Holzlöhner, C. R. Menyuk, and W. L. Kath, “Efficient and accurate computation of eye diagrams and bit error rates in a single-channel CRZ system,” IEEE Photon. Technol. Lett. 14, 1079–1081 (2002).
[Crossref]

R. Holzlöhner, C. R. Menyuk, and W. L. Kath, “A covariance matrix method to compute bit error rates in a highly nonlinear dispersion-managed soliton system,” IEEE Photon. Technol. Lett. 15, 688–690 (2003).
[Crossref]

J. Comput. Phys. (1)

A. Demir, “Non-Monte Carlo formulations and computational techniques for the stochastic nonlinear Schrodinger equation,” J. Comput. Phys. 201, 148–171 (2004).
[Crossref]

J. Lightwave Technol. (11)

J. Hult, “A fourth-order Runge-Kutta in the interaction picture method for simulating supercontinuum generation in optical fibers,” J. Lightwave Technol. 25, 3770–3775 (2007).
[Crossref]

E. Forestieri, “Evaluating the error probability in lightwave systems with chromatic dispersion, arbitrary pulse shape and pre- and postdetection filtering,” J. Lightwave Technol. 18, 1493–1503 (2000).
[Crossref]

P. Serena, A. Bononi, J. C. Antona, and S. Bigo, “Parametric gain in the strongly nonlinear regime and its impact on 10-Gb/s NRZ systems with forward-error correction,” J. Lightwave Technol. 23, 2352–2363 (2006).
[Crossref]

A. Bononi, P. Serena, and A. Orlandini, “A unified design framework for single-channel dispersion-managed terrestrial systems,” J. Lightwave Technol. 26, 3617–3631 (2008).
[Crossref]

E. Ip and J. M. Kahn, “Power spectra of return-to-zero optical signals,” J. Lightwave Technol. 24, 1610–1618 (2006).
[Crossref]

M. Secondini, E. Forestieri, and C. R. Menyuk, “A combined regular-logarithmic perturbation method for signal-noise interaction in amplified optical systems,” J. Lightwave Technol. 27, 3358–3369 (2009).
[Crossref]

R. Hui, M. O’Sullivan, A. Robinson, and M. Taylor, “Modulation instability and its impact in multispan optical amplified IMDD systems: theory and experiments,” J. Lightwave Technol. 15, 1071–1082 (1997).
[Crossref]

P. Serena, A. Orlandini, and A. Bononi, “Parametric-gain approach to the analysis of single-channel DPSK/DQPSK systems with nonlinear phase noise,” J. Lightwave Technol. 24, 2026–2037 (2006).
[Crossref]

A. Demir, “Nonlinear phase noise in optical-fiber-communication systems,” J. Lightwave Technol. 25, 2002–2032 (2007).
[Crossref]

E. Forestieri and G. Prati, “Exact analytical evaluation of second-order PMD impact on the outage probability for a compensated system,” J. Lightwave Technol. 22, 988–996 (2004).
[Crossref]

R. Holzlöhner, V. S. Grigoryan, C. R. Menyuk, and W. L. Kath, “Accurate calculation of eye diagrams and bit error rates in optical transmission systems using linearization,” J. Lightwave Technol. 20, 389–400 (2002).
[Crossref]

Opt. Express (3)

Opt. Lett. (1)

Other (4)

R. Holzlöhner, “A covariance matrix method for the computation of bit errors in optical transmission systems,” Ph.D. thesis, University of Maryland Baltimore County. Baltimore, Maryland, USA (2003).

D. Gariepy and G. He, “Measuring OSNR in WDM systemsEffects of resolution bandwidth and optical rejection ratio,” White paper, EXFO Inc. (2009).

A. Papoulis, Probability, Random Variables, and Stochastic Processes (McGraw-Hill, 1991).

Z. Zhang, L. Chen, and X. Bao, “The noise propagator in an optical system using EDFAs and its effect on system performance: accurate evaluation based on linear perturbation,” arXiv:physics.optics/1207.3362v1.

Cited By

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

Alert me when this article is cited.


Figures (4)

Fig. 1
Fig. 1 Low-pass equivalent optical model. The receiver consists of optical and electrical filters and DPSK balance detection. In this work, calculations, except those in Sec. 5.1, are based on the assumptions N0k = N0 and Gk = G [k = 1, 2,...(K + 1)]. Also, the fiber in each span is assumed to have same length (L km) and same loss (α dB/km). Nin is the ASE noise added at the transmitter. Changing Nin will change the OSNR at the receiver. In a typical balanced DPSK receiver, the delay in one branch of the interferometer is Tb = 1/Rb. Here, Rb = 20 Gb/s.
Fig. 2
Fig. 2 BER versus received OSNR for a 20Gb/s 20-span RZ-DPSK system with Φ̄NL = 0.2π. Solid: obtained using CMM of Ref. [1]. Dashed (dotted): improved CW approach of Ref. [5] with CW power being peak power (average power), respectively. Dash-dotted: RK4IP approach. All curves, except the dash-dotted, are obtained from Fig. 8 of Ref. [5].
Fig. 3
Fig. 3 BER versus received OSNR for the multi-span RZ-50% DPSK systems with Φ̄NL = 0.9. Each span consists of a SMF fiber (42 km long) and a DCF file (7 km long). Other fiber parameters were detailed in Table 1 of Ref. [7]. According to Fig. 7(b) in Ref. [7], where the experimental curves for 5-span, 10-span, and 25-span systems were almost the same, here we replot these three curves using a thick solid curve.
Fig. 4
Fig. 4 BER vs RX-OSNR for the 5-span RZ-50% DPSK system with Φ̄NL = 0.9. Other parameters are same as those used in Fig. 3. Solid: experimental results. Dotted (Dash-dotted): numerical calculation using RK4IP with (without) phase shift Δ.

Equations (43)

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

d a ˜ I d z = N ^ I a ˜ I
a ˜ n + 1 = e L ^ h / 2 [ a ˜ n I + h k 1 6 + h k 2 3 + h k 3 3 + h k 4 6 ] a ˜ n I = e L ^ h / 2 a ˜ n k 1 = N ^ n I a ˜ n I = e L ^ h / 2 N ^ ( z n ) a ˜ n k ^ 1 a ˜ n k 2 = N ^ n + 1 / 2 I [ a ˜ n I + h k 1 2 ] = N ^ ( z n + h / 2 ) [ e L ^ h / 2 + h k ^ 1 2 ] a ˜ n k ^ 2 a ˜ n k 3 = N ^ n + 1 / 2 I [ a ˜ n I + h k 2 2 ] = N ^ ( z n + h / 2 ) [ e L ^ h / 2 + h k ^ 2 2 ] a ˜ n k ^ 3 a ˜ n k 4 = N ^ n + 1 I [ a ˜ n I + h k 3 ] = e L ^ h / 2 N ^ ( z n + 1 ) e L ^ h / 2 [ e L ^ h / 2 + h k ^ 3 ] a ˜ n k ^ 4 a ˜ n
a ˜ n + 1 = ( e L ^ h / 2 [ e L ^ h / 2 + h k ^ 1 6 + h k ^ 2 3 + h k ^ 3 3 ] + h 6 N ^ ( z n + 1 ) e L ^ h / 2 ( e L ^ h / 2 + h k ^ 3 ) ) a ˜ n ,
H ( z n + 1 , z n ) = e L ^ h / 2 [ e L ^ h / 2 + 1 3 ( h 2 k ^ 1 + 2 h 2 k ^ 2 + h k ^ 3 ) ] + h 6 N ^ ( z n + 1 ) e L ^ h / 2 [ e L ^ h / 2 + h k ^ 3 ]
p n ( L , 0 ) = H ( L , L h L ) H ( h 1 , 0 ) .
PG 1 = p n ( L , 0 ) σ 2 I p n T ( L , 0 ) = σ 2 p n ( L , 0 ) p n T ( L , 0 ) ,
σ 2 = N 0 / ( 2 T 0 ) [ N 0 = n s p ( G 1 ) h ¯ ω ]
P G ^ = ( G σ in 2 + σ 2 ) P n ( K , 0 ) P n T ( K , 0 ) + σ 2 k = 1 K P n ( K , k ) P n T ( K , k ) , P n ( K , k ) = p n ( K L , ( K 1 ) L ) p n ( ( k + 1 ) L , k L ) , ( k = 0 , 1 , , K 1 ) P n ( K , K ) = I ,
P G ^ = σ 2 P n , eq P n , eq T ,
E [ exp ( s ( λ ˜ c 2 + 2 c b ˜ ) ) ] = d c 2 π σ 2 exp ( c 2 2 σ 2 ) exp [ s ( λ ˜ c 2 + 2 c b ˜ ) ] = exp [ 2 σ 2 s 2 b ˜ 2 1 2 σ 2 s λ ˜ ] 1 2 σ 2 s λ ˜ ,
Ψ t s ( s ) = E [ exp [ s I ( t s ) ] ] = exp [ s I ss ( t s ) ] i = 1 4 M n + 2 exp [ 2 σ 2 s 2 b ˜ i 2 ( t s ) 1 s β i ] ( 1 s β i ) ξ , ( β i = 2 σ 2 λ ˜ i )
OSNR 0.1 nm = P ¯ s P ASE ( B r ) .
OSNR L , 0.1 nm = P ¯ P ASE ( B m ) × B m B r ( P ASE ( B m ) = [ GN in + N 0 ( K + 1 ) ] B m ) ,
OSNR N L , 0.1 nm = P ¯ s P ASE ( B m ) × B m B r ( P ASE ( B m ) = T r ( O m T P G ^ O m ) σ 2 = T r ( P G ^ O m O m T ) σ 2 ) .
P ASE ( B m , ± Δ λ ) = T r [ P G ^ ( O m ( Δ λ ) O m T ( Δ λ ) + O m ( + Δ λ ) O m T ( + Δ λ ) ) ] σ 2 2 ,
D l l s s = e j 2 π l N + e j 2 π l N 2 , D m m n n = e j 2 π m T b T 0 + e j 2 π m T b T 0 2 , D m l n s = e j 2 π m T b T 0 j Δ + e j 2 π l N + j Δ 2 ,
Δ = ϕ 0 < δ ϕ > ,
u z = j β ω ω 2 2 u t 2 + β ω ω ω 6 3 u t 3 j γ | u | 2 u α 2 u ,
v z = j β ω ω 2 2 v t 2 + β ω ω ω 6 3 v t 3 j e α z γ | v | 2 v .
v z = j β ω ω 2 2 v t 2 + β ω ω ω 6 3 v t 3 j e α z γ | v | 2 v j w ( z , t ) ,
r ( z , z , t , t ) = E { w ( z , t ) w * ( z , t ) } = δ ( t t ) δ ( z z ) k = 1 K + 1 N 0 k δ ( z ( k 1 ) L ) .
N 0 k = N 0 = n s p ( G 1 ) h ¯ ω ( k = 1 , 2 , , K , ( K + 1 ) )
v 0 z = j β ω ω 2 2 v 0 t 2 + β ω ω ω 6 3 v 0 t 3 j e α z γ | v 0 | 2 v 0
δ v z = j β ω ω 2 2 δ v t 2 + β ω ω ω 6 3 δ v t 3 j 2 e α z γ | v 0 | 2 δ v j e α z γ v 0 2 δ v * j w ( z , t ) .
ν l = ν ( ω l ) = e α z | v 0 | 2 e j ω l t d t , μ l = μ ( ω l ) = e α z v 0 2 e j ω l t d t ,
d a l d z = j β ω ω 2 ω l 2 a l j β ω ω ω 6 ω l 3 a l j 2 γ ( z ) [ M ν ] l m a m j γ ( z ) [ M μ ] l m a m * j W l ,
d a d z = L ¯ a + ν a + μ a * j W
L ¯ = j L C D , ν = 2 j γ ( z ) ( M ν R + j M ν I ) , μ = j γ ( z ) ( M μ R + j M μ I )
d a ˜ d z = ( L ^ + ν ^ + μ ^ ) a ˜ + W ˜
L ^ = ( 0 L C D L C D 0 ) , ν ^ = ( ν A A ν S S ν S S ν A A ) , μ ^ = ( μ R μ I μ I μ R ) ,
E { W ˜ l ( z ) W ˜ l * ( z ) } = δ z , z δ l , l k = 0 K N 0 2 T 0 δ z , k L ( l = 1 , , 4 M n + 2 ) ,
d a ˜ d z = ( L ^ + N ^ ) a ˜ ; ( N ^ = ν ^ + μ ^ )
E { W ˜ ( z f ) W ˜ * ( z f ) } = N 0 2 T 0 I = σ 2 I ,
y s s ( t s ) = [ s o ( t s + T b ) | R s s | s o ( t s ) + c . c . ] / 2 = s o ( t s ) | R s s D | s o ( t s ) y n n ( t s ) = [ n o ( t s + T b ) | R n n | n o ( t s ) + c . c . ] / 2 = N o | R n n D | N o = Z | Λ | Z y n s ( t s ) = [ n o ( t s + T b ) | R n s | s o ( t s ) + n o ( t s ) | R n s | s o ( t s + T b ) + c . c . ] / 2 = [ N in | O n n R n s D | s o ( t s ) + c . c . ] = [ Z | b D ( t s ) + c . c . ]
D l l s s = e j 2 π l N + e j 2 π l N 2 , D m m n n = e j 2 π m T b T 0 + e j 2 π m T b T 0 2 , D m l n s = e j 2 π m T b T 0 + e j 2 π l N 2 .
L s = η N T b B o , M n = η B o T 0 , T 0 = μ ( 1 B o + 1 B r ) .
| s ˜ o = [ Re { | s o } Im { | s o } ] , | a ˜ 0 = [ Re { | a 0 } Im { | a 0 } ] , | Z ˜ = [ Re { | Z } Im { | Z } ] ,
y n n ( t s ) = a ˜ 0 | U ˜ Λ ˜ U ˜ T | a ˜ 0 = Z ˜ | Λ ˜ | Z ˜ ( Λ ˜ = U ˜ T P n , eq T O ˜ n m T P n , eq U ˜ = diag { λ ˜ 1 , , λ ˜ 4 M n + 2 } ) y n s ( t s ) = Z ˜ | U ˜ T P ˜ n , eq T O ˜ n m T R ˜ n s D O ˜ s s | s ˜ o ( t s ) Z ˜ | B ˜ | s ˜ o ( t s ) Z ˜ | b ˜ ( t s ) ( B ˜ = U ˜ T P n , eq T O ˜ n n T R ˜ n s D O ˜ s s ) ,
ϕ N L = ( P ¯ + δ P ¯ ) γ K L .
< δ ϕ 2 > = < ϕ N L 2 Φ ¯ N L 2 > 2 P ¯ < δ P ¯ > ( γ K L ) 2 .
< δ ϕ > 2 P ¯ γ K L / P ¯ / < δ P ¯ > = 2 Φ ¯ N L / OSNR ,
< δ ϕ 1 / OSNR π / 2 arctan ( OSNR ) ,
< δ ϕ 2 > = 2 [ G N in B in + N 0 B l k K 3 ( 1 + 1 K ) ( 1 + 1 2 K ) ] Φ ¯ N L 2 / P ¯ ,

Metrics