Electronic compensation of chromatic dispersion using a digital ...

5 downloads 4912 Views 2MB Size Report
[email protected]. Abstract: Digital signal processing (DSP) combined with a phase and ..... is the complex sign of the symbol i. X defined by. ⎪. ⎪. ⎩. ⎪. ⎪. ⎨.
Electronic compensation of chromatic dispersion using a digital coherent receiver Seb J. Savory, Giancarlo Gavioli, Robert I. Killey, Polina Bayvel Optical Networks Group, Dept. of Electronic & Electrical Engineering, University College London., Torrington Place, London WC1E 7JE,UK [email protected]

Abstract: Digital signal processing (DSP) combined with a phase and polarization diverse coherent receiver is a promising technology for future optical networks. Not only can the DSP be used to remove the need for dynamic polarization control, but also it may be utilized to compensate for nonlinear and linear transmission impairments. In this paper we present results of a 42.8Gbit/s nonlinear transmission experiment, using polarization multiplexed QPSK data at 10.7GBaud, with 4 bits per symbol. The digital coherent receiver allows 107,424 ps/nm of chromatic dispersion to be compensated digitally after transmission over 6400km of standard single mode fiber. ©2007 Optical Society of America OCIS codes: (060.1660) Coherent communications; (060.4510) Optical communications;

References and links 1. 2.

3. 4.

5. 6. 7. 8.

Y. Han and G. Li, "Coherent optical communication using polarization multiple-input-multiple-output," Opt. Express 13, 7527-7534 (2005). S. Tsukamoto, D.-S Ly-Gagnon, K. Katoh, and K. Kikuchi, “Coherent demodulation of 40-Gbit/s polarization-multiplexed QPSK signals with 16-GHz spacing after 200-km transmission”, in Proceedings of Optical Fiber Communications Conference 2005, paper PDP-29. M. G. Taylor, "Coherent detection method using DSP for demodulation of signal and subsequent equalization of propagation impairments,” IEEE Photon. Technol. Lett. 16, 674 – 676 (2004). D. McGhan, C. Laperle, A. Savchenkov, C. D. Li, G. Mak, and M. O'Sullivan, “5120 km RZ-DPSK transmission over G652 fiber at 10 Gb/s without optical dispersion compensation,” IEEE Photon. Technol. Lett. 18, 400 - 402 (2006). A. Hardy and W. Streifer, “Coupled mode solutions of multiwaveguide systems,” IEEE J. Quantum Electron. 22, 528 – 534 (1986). D. Godard, “Self-recovering equalization and carrier tracking in two-dimensional data communication systems,” IEEE Trans. Commun. 28, 1867 – 1875 (1980). J. G. Proakis, Digital Communications, (McGraw-Hill, 2001), Chap. 11. M. Jeruchim, “Techniques for estimating the Bit Error Rate in the simulation of Digital Communication Systems,” IEEE J. Sel. Areas Commun. 2, 153 – 170 (1984).

1. Introduction Phase and polarization diverse coherent detection combined with digital signal processing (DSP) is a powerful combination for next generation optical receivers. Not only does it remove the need for dynamic polarization control [1], but also enables the detection of polarization multiplexed spectrally efficient formats [2]. Furthermore since the receiver is linear in field, dispersion from the fiber can be compensated using digital finite impulse response (FIR) filters [3], in principle allowing for transmission over standard fiber without optical dispersion compensation. To date the longest distance reported for a 10Gbit/s system operating over standard fiber without optical dispersion compensation is 5120km using pre-distortion [4]. In this paper we demonstrate that a distance of 6400km can be achieved for a 42.8Gbit/s signal, by transmitting 10.7GBaud polarization multiplexed QPSK (PMQPSK) with 4 bits per symbol. This represents the longest reach achieved for transmission without optical dispersion #78394 - $15.00 USD

(C) 2007 OSA

Received 21 December 2006; revised 16 February 2007; accepted 18 February 2007

5 March 2007 / Vol. 15, No. 5 / OPTICS EXPRESS 2120

compensation, not only for a 40Gbit/s system, but also for any 10GBaud system operating over standard fiber. 2. Coherent receiver architecture 2.1 Overview The coherent receiver illustrated in Fig. 1 is based on an intradyne system, using a local oscillator which is nominally at the same frequency as the incoming signal. Ex i + E y j

Polarization beam splitter Photodiode

Re(E x )

Photodiode

Im(E x )

Input optical signal Polarization controller I/Q coupler Photodiode

Re(E y )

Photodiode

Im(E y )

Optical local oscillator Polarization controller

Fig. 1 Schematic of phase and polarization diverse receiver

The I/Q coupler acts as a 90 degree hybrid, separating the incoming optical signal into inphase and quadrature components, corresponding to the real and imaginary parts of a complex signal the analysis of which is presented in the next subsection. 2.2 Fiber I/Q coupler based hybrid In our receiver we used a 3x3 coupler as the 90 degree hybrid. The coupling between three identical fibers, may be modeled by the following matrix equation [5]

dE = jKE dz

(1)

Where E = (E1 , E 2 , E3 ) and K is a matrix in which all elements are equal to the coupling coefficient κ which determines the coupling between the three electric fields E1 , E 2 and E 3 . The solution to this equation may be expressed in term of the matrix T

exponential such that E( z ) = exp( jKz )E( z = 0) . If we choose a coupling distance L given by

κL = [2nπ + 2 arctan(3)] / 3 , where n is an integer, the coupler splits power in

(

)

the ratio 1:2:2, hence for input fields E( z = 0) = E sig ,0, E lo exp(− 3 jπ / 4 ) the output T

fields are given by

− 1 + j − 1 + j ⎞⎛ E sig ⎞ ⎛ 1 ⎟⎜ ⎟ 2+ j⎜ 1 0 − 1 + j ⎟⎜ E( z = L ) = ⎟ ⎜ −1 + j 5 ⎜ ⎟ ⎜ 1 ⎠⎝ Elo exp(− 3 jπ / 4)⎟⎠ ⎝ −1 + j −1 + j

(2)

These resulting electric fields are then detected using a photodiode, with the corresponding

#78394 - $15.00 USD

(C) 2007 OSA

Received 21 December 2006; revised 16 February 2007; accepted 18 February 2007

5 March 2007 / Vol. 15, No. 5 / OPTICS EXPRESS 2121

(

photocurrents (i1, i2 , i3 )T = E1 ( L) 2 , E2 ( L) 2 , E3 ( L) 2 ⎛ i1 ⎞ ⎜ ⎟ ⎜ i2 ⎟ = ⎜i ⎟ ⎝ 3⎠

2 2 5

⎛ ⎜ ⎜ − Re ⎜ ⎝

(

) ( )

Re E sig Elo* E sig E lo* − Im E sig E lo* Im E sig E lo*

(

)

⎞ ⎟ 1 ⎟+ ⎟ 5 ⎠

)

coherently( detected  terms

)

T

given by 2 2 ⎛ E + 2 E lo ⎞⎟ sig ⎜ 2 2⎟ ⎜ 2 E sig + 2 Elo ⎟ ⎜ 2 2 ⎜ 2E + E lo ⎟⎟ ⎜ sig ⎝ ⎠





(3)



directly detected terms

From Eq. (3) it can be seen that the coherently detected currents associated with i1 and i3 are in quadrature. Provided that E lo

2

2

>> E sig , if a DC block is employed, the directly detected

terms may be neglected, allowing phase diversity using just two photodiodes. 3. Digital signal processing algorithms 3.1 Overview DSP is used to equalize the channel and recover the phase of the signal. This allows for chromatic dispersion and polarization dependent effects to be compensated, in addition to tracking the phase difference between the incoming signal and the local oscillator. To achieve this, the samples pass through four distinct stages. Firstly the data is resampled to give ~2 samples per symbol, after which the signals are equalized using FIR filters. Following this the phase of the signal is recovered, before finally a threshold is applied to recover the data bits. 3.2 Details of the algorithms Most of the distortion occurs due to chromatic dispersion. By using an FIR filter with a quadratic phase response of value equal and opposite to that imposed by the chromatic dispersion, this may be compensated. The FIR filter coefficients (512 taps) are calculated from a truncated impulse response of the chromatic dispersion with a Kaiser window function. Polarization tracking and compensation of the remaining distortions are achieved using four shorter filters hxx, hxy, hyx and hyy, having just 5 taps each (corresponding to the elements in the inverse Jones matrix). If x and y represent the incoming electrical signal for the two polarizations, then the output signals for the two polarization are x ′ = h xx ⋅ x + h xy ⋅ y and

y ′ = h yx ⋅ x + h yy ⋅ y . These filters are initially adapted using the constant modulus algorithm[6], which exploits the fact that QPSK is “constant modulus” and so minimizes the errors

ε x = 1 − x′

2

, and

ε y = 1 − y′

2

using the stochastic gradient such that

h xy = h xy + με x x ′ ⋅ y

h xx = h xx + με x x ′ ⋅ x

(4)

h yx = h yx + με y y ′ ⋅ x h yy = h yy + με y y ′ ⋅ y where μ is the convergence parameter and x denotes the complex conjugate of x . Once the constant modulus algorithm has converged the offset frequency for the two polarizations are estimated using a 4th order nonlinearity to remove the QPSK data. The frequency estimated from the maximum of FFT ( x ′ ) and FFT ( y ′ ) respectively, with 4

4

the maximum frequency offset that can be resolved being ±2.675 GHz, a quarter of 10.7GHz Having removed the offset frequency the receiver moves into a decision directed mode. The phase-corrected output symbols are given by X = exp (− j φ x )x ′ and Y = exp − jφ y y ′ ,

(

)

where the phases are given by

#78394 - $15.00 USD

(C) 2007 OSA

Received 21 December 2006; revised 16 February 2007; accepted 18 February 2007

5 March 2007 / Vol. 15, No. 5 / OPTICS EXPRESS 2122

∑ N

φ x = arg

i =1

φ y = arg

csgn ( X i ) ⋅X i

∑ csgn (Y ) ⋅Y N

i

i =1

(5)

i

where csgn( X i ) is the complex sign of the symbol X i defined by

⎧ 1 + j [Re( x ) > 0, Im( x ) > 0] ⎪ 1− j Re( x ) > 0, Im( x ) < 0 ] ⎪ csgn( x ) = ⎨ ⎪− 1 + j Re( x ) < 0, Im( x ) > 0 ] ⎪⎩− 1 − j [Re( x ) < 0, Im( x ) < 0]

(6)

The tap weights of the FIR taps adapted using a decision directed mode such that

h xx = h xx + με x ⋅ x

h xy = h xy + με x ⋅ y

h yx = h yx + με y ⋅ x

h yy = h yy + με y ⋅ y

where

the

errors

for

(

the

decision

ε x = exp ( jφ x ) csgn ( X ) / 2 − X

) and ε

directed

y

(7)

algorithm

are

(

given

= exp ( jφ y ) csgn (Y ) / 2 − Y

)

by

4. Experimental investigation into the nonlinear propagation of PMQPSK The aim of the experiment was to determine, the maximum distance that can be achieved using 42.8Gbit/s polarization multiplexed QPSK, given that the fiber is nonlinear. 10.7Gb/s PPG PRBS length=11

42.8Gb/s PMQPSK

______ Delayed DATA DATA PC DQPSK Tx

PBS delay ~10ns

TX AOM LOOP AOM

PC

80km SMF

21.4Gb/s QPSK

10/90

Tekronix TDS6154C

PBS PIN PIN

50GHz AWG PC

PIN LO

PIN PC

Fig. 2. Experimental setup, with 42.8 Gbit/s obtained by polarization multiplexing two 21.4Gbit/s NRZ-QPSK signals. PBS: Polarization beam splitter; PC: Polarization controller; SMF: Single mode fiber; AOM: Acousto-optic modulator; PPG: Pulse pattern generator; AWG: Array waveguide filter.

#78394 - $15.00 USD

(C) 2007 OSA

Received 21 December 2006; revised 16 February 2007; accepted 18 February 2007

5 March 2007 / Vol. 15, No. 5 / OPTICS EXPRESS 2123

4.1 Experimental setup The experimental setup used (Fig. 2), consists of three distinct elements, the transmitter, the recirculating loop, and the coherent receiver. The DFB-based transmitter with wavelength of 1554.94nm and measured linewidth of 2MHz was operated at 10.7Gbaud, with a 211-1 PRBS, chosen to facilitate processing of the data in a realistic time frame. The data rate of 42.8Gbit/s was achieved using QPSK and polarization multiplexing. Within the recirculating loop the 80km fiber span generates 1342.8ps/nm of chromatic dispersion with a loss of 15.7dB. Two erbium doped fiber amplifiers, are used within the loop to overcome the loss of the fiber and loop components. At the receiver a DFB local oscillator with a linewidth of 100kHz was combined with the signal via two asymmetric 3 fiber couplers, with coupling ratios 1:2:2, to detect the in-phase and quadrature components of the two polarizations. On conversion into the electrical domain, the signal was digitized at 20GSa/s using a Tektronix TDS6154C digital storage oscilloscope with 8 bit resolution, with the waveforms then processed off-line using Matlab. 4.2 Initial characterization To determine the penalty due to transmission we performed back-to-back measurements which gave the required OSNR (ROSNR) at the FEC limit (BER=3x10-3) as 10.2dB (in 0.1nm). Having determined the ROSNR we next determined how noise accumulated with transmission distance. Operating the loop with a launch power of -5dBm we determined that the OSNR after N recirculations was given by (8) OSNR = Plaunch + 37.5 − 10 log 10 ( N ) With Plaunch being the launch power in dBm, with 37.5 arising due to the experimentally determined parameter that -37.5dBm of optical noise (in 0.1nm) is added per recirculation. 4.3 Characterizing the nonlinear transmission To obtain the optimal fiber launch power for maximum distance, the BER was measured at distances (in km) of {1280, 2480, 3840, 5120} as the power level was increased from -15 to 5 dBm, in steps of 2dB. To determine the point of convergence, additional measurements were taken after {5720, 6400, 7040} km as the power was increased from -9 to -1 dBm again in 2 dB steps. The DSP algorithms used the same parameter set and converged in all cases. The BER was obtained by counting errors over a total of 73692 bits, from which a contour plot of BER=3x10-3 based on interpolation of 20logQ was obtained, which closely followed that obtained in preliminary tests. From Fig. 3, we determine that the maximum reach of 6480km is obtained with a fiber launch power of -7dBm. At the optimum launch power after 6480km (81 recirculations), Eq. (8) gives the OSNR as 11.4dB, indicating that the OSNR penalty is 1.2 dB compared to the back-to-back performance. While these maximum distances are based on interpolation, the optimum point from Fig. 3 is close to a measurement point taken at 6400km with launch power of -7dBm. At this point, as illustrated in Fig. 4, the constellations are successfully recovered, with 176 errors counted over 73692 bits giving a BER of 2.4x10-3 with the 99% confidence interval being (2.1x10-3, 2.7x10-3) [8].

#78394 - $15.00 USD

(C) 2007 OSA

Received 21 December 2006; revised 16 February 2007; accepted 18 February 2007

5 March 2007 / Vol. 15, No. 5 / OPTICS EXPRESS 2124

contour of BER=3x10 -3

2 Launch power (dBm)

0

Nonlinearity limited

-2 -4 -6

-7dBm optimum power

Linear noise limited

-8 -10

maximum reach=6480km (81 recirculations)

-12 0

1000

2000

3000 4000 5000 Distance (km)

6000

7000

1.5

1.5

1

1

Quadrature (Y polarization)

Quadrature (X polarization)

Fig. 3. Contour plot of launch power versus distance for BER=3x10-3, indicating a maximum reach of 6480km

0.5 0 -0.5 -1 -1.5 -1.5

-1

-0.5 0 0.5 1 In-phase (X polarization)

1.5

0.5 0 -0.5 -1 -1.5 -1.5

-1

-0.5 0 0.5 1 In-phase (Y polarization)

1.5

Fig 4. Recovered constellation diagrams with -7dBm launch power after 6400km transmission, having compensated for 107,424 ps/nm of chromatic dispersion digitally with an overall BER=2.4x10-3

5. Conclusion In this paper we have demonstrated 42.8Gbit/s per wavelength transmission with a record distance of 6400km over standard fiber with no optical dispersion compensation. A record total dispersion of 107,424ps/nm was compensated using digital signal processing, with an OSNR penalty of 1.2dB. This is the first time that 40Gbit/s per wavelength data has been transmitted over such a distance on standard fiber, without the use of optical dispersion compensation and exceeds the distance previously reported for 10Gbit/s systems. Acknowledgments Financial support to SJ Savory and G Gavioli from the Leverhulme Trust is gratefully acknowledged, as is that received for the project from FP6 IST-IP NOBEL and EPSRC. The

#78394 - $15.00 USD

(C) 2007 OSA

Received 21 December 2006; revised 16 February 2007; accepted 18 February 2007

5 March 2007 / Vol. 15, No. 5 / OPTICS EXPRESS 2125

authors thank Dr R Griffin and Bookham for the loan of the DQPSK transmitter, AD Stewart and S Woods for their assistance and to Dr MG Taylor for valuable discussions.

#78394 - $15.00 USD

(C) 2007 OSA

Received 21 December 2006; revised 16 February 2007; accepted 18 February 2007

5 March 2007 / Vol. 15, No. 5 / OPTICS EXPRESS 2126