PMD compensation in multilevel codedmodulation schemes with coherent detection using BLAST algorithm and iterative polarization cancellation

Transcription

1 PMD compensation in multilevel codedmodulation schemes with coherent detection using BLAST algorithm and iterative polarization cancellation Ivan B Djordjevic, Lei Xu*, and Ting Wang* University of Arizona, Department of Electrical and Computer Engineering, Tucson, AZ 85721, USA ivan@ecearizonaedu * NEC Laboratories America, Princeton, NJ 08540, USA Abstract: We present two PMD compensation schemes suitable for use in multilevel (M 2) bloc-coded modulation schemes with coherent detection The first scheme is based on a BLAST-type polarization-interference cancellation scheme, and the second scheme is based on iterative polarization cancellation Both schemes use the LDPC codes as channel codes The proposed PMD compensations schemes are evaluated by employing coded-ofdm and coherent detection When used in combination with girth-10 LDPC codes those schemes outperform polarization-time coding based OFDM by 1 db at BER of 10-9, and provide two times higher spectral efficiency The proposed schemes perform comparable and are able to compensate even 1200 ps of differential group delay with negligible penalty 2008 Optical Society of America OCIS codes: ( ) Optical communications; ( ) Polarization mode dispersion (PMD); ( ) Modulation; ( ) Multiplexing; ( ) Orthogonal frequency division multiplexing; ( ) Low-density parity-chec (LDPC) codes, ( ) Bell Laboratories layered space-time architecture (BLAST); ( ) Iterative polarization interference cancellation References and Lins 1 I B Djordjevic, M Cvijetic, L Xu, and T Wang, Using LDPC-coded modulation and coherent detection for ultra high-speed optical transmission, J Lightwave Technol 25, (2007) 2 L L Minov, I B Djordjevic, H G Batshon, L Xu, T Wang, M Cvijetic, and F Kueppers, Demonstration of PMD compensation by LDPC-coded turbo equalization and channel capacity loss characterization due to PMD and quantization, IEEE Photon Technol Lett 19, (2007) 3 W Shieh, X Yi, Y Ma, and Y Tang, Theoretical and experimental study on PMD-supported transmission using polarization diversity in coherent optical OFDM systems, Opt Express 15, (2007) 4 H Sun, K -T Wu, and K Roberts, Real-time measurements of a 40 Gb/s coherent system, Opt Express 16, (2008) 5 S Alamouti, A simple transmit diversity technique for wireless communications, IEEE J Sel Areas Commun 16, (1998) 6 I B Djordjevic, L Xu, and T Wang, PMD compensation in multilevel coded-modulation schemes with coherent detection using Alamouti-type polarization-time coding, in Proc IEEE LEOS Summer Topicals 2008, pp , July G J Foschini, Layered space-time architecture for wireless communication in a fading environment when using multi-element antennas, Bell Labs Tech J 1, (1996) 8 E Biglieri, R Calderban, A Constantinides, A Goldsmith, A Paulraj, and H V Poor, MIMO Wireless Communications (Cambridge University Press, Cambridge 2007) 9 I B Djordjevic, S Denic, J Anguita, B Vasic, and M A Neifeld, LDPC-coded MIMO optical communication over the atmospheric turbulence channel, J Lightwave Technol 26, (2008) 10 D Pennincx, and V Morenás, Jones matrix of polarization mode dispersion, Opt Lett 24, (1999) 11 I B Djordjevic, L Xu, T Wang, and M Cvijetic, Large girth low-density parity-chec codes for long-haul high-speed optical communications, in Proc OFC/NFOEC 2008, Paper no JWA53 (C) 2008 OSA 15 September 2008 / Vol 16, No 19 / OPTICS EXPRESS 14845

2 1 Introduction The bit-error ratio (BER) performance of fiber-optic communication systems operating at high data rates is degraded by intrachannel and interchannel fiber nonlinearities, polarization mode dispersion (PMD), and chromatic dispersion [1] To deal with PMD a number of methods have been proposed recently, four of them seem to be able successfully to tacle the PMD effects: (i) turbo equalization [2], (ii) polarization diversity orthogonal frequency division multiplexing (OFDM) [3], (iii) the channel equalization scheme described in [4], and Alamouti-type [5] polarization-time (PT) coding scheme introduced by authors in [6] In this paper we propose two alternative schemes suitable for PMD compensation, which do not require the increase of complexity as the differential group delay (DGD) increases The first scheme is based on the Bell Laboratories layered space-time architecture (BLAST) [7], originally proposed to deal with spatial interference in wireless communications We consider two versions of this scheme [8]: (a) the zero-forcing vertical-blast scheme (ZF V-BLAST), and the minimum-mean-square-error vertical-blast (MMSE V-BLAST) scheme Because the ZF V-BLAST scheme is derived by ignoring the influence of amplified spontaneous emission (ASE) noise, we proposed the second scheme that uses the output of ZF V-BLAST scheme as starting point and removes the remaining polarization interference in an iterative fashion This approach also leads to reducing the influence of ASE noise We evaluate the performance of those schemes when used in combination with coherent detection based OFDM We describe how to use those schemes together with multilevel modulation and forward error correction (FEC) The arbitrary FEC scheme can be used with proposed PMD compensation schemes, however, the use of low-density parity-chec (LDPC) codes leads to near channel capacity achieving performance [9] The proposed schemes outperform the polarization-diversity OFDM scheme [3], and PT-based OFDM [6] in terms of both BER and spectral efficiency 2 Description of proposed PMD compensation schemes For the first-order PMD study the Jones matrix, neglecting the polarization dependent loss and depolarization effects, can be represented by [10] h h e H = = RP R, P = jωτ /2 xx xy -1 ( ω) ( ω) jωτ /2 hyx h yy 0 e 0, (1) where τ denotes DGD, ω is the angular frequency, and R=R(θ,ε) is the rotational matrix [10] θ jε /2 θ jε/2 cos e sin e 2 2 R =, θ jε /2 θ jε /2 sin e cos e 2 2 with θ being the polar angle, and ε being the azimuth angle For the OFDM with coherent detection, the received symbol vector of th subcarrier in ith OFDM symbol r i, =[r x,i, r y,i, ] T can be represented by j CD ( ) T LO i, i, e φ + φ φ = + i, r H s n, (2) where s i, =[s x,i, s y,i, ] T denotes the transmitted symbol vector of th subcarrier in ith OFDM symbol, for both polarizations, n i, =[n x,i, n y,i, ] T denotes the noise vector dominantly determined by the amplified spontaneous emission (ASE) noise; φ T and φ LO denote the laser phase noise processes of transmitting and local lasers, φ CD () denotes the phase distortion of th subcarrier due to chromatic dispersion (CD), and the Jones matrix of th subcarrier H is already introduced in (1) The transmitted/received symbols are complex-valued, with real part corresponding to the in-phase coordinate and imaginary part corresponding to the (C) 2008 OSA 15 September 2008 / Vol 16, No 19 / OPTICS EXPRESS 14846

3 quadrature coordinate Figure 1 shows the magnitude responses of h xx and h xy coefficients of the Jones channel matrix against a normalized frequency fτ (the frequency is normalized with DGD τ so that the conclusions are independent on the data rate) for two different cases: (a) θ=π/2 and ε=0, and (b) θ=π/3 and ε=0 In the first case channel coefficient h xx completely fades away for certain frequencies, while in the second case it never completely fades away; suggesting that the first case represents the worst case scenario To avoid this problem, in direct detection OFDM systems someone can redistribute the transmitted power among subcarriers not being under fading, or use the polarization-diversity coherent detection OFDM [3] We propose two alternative approaches instead, which can be used for a number of modulation formats including M-ary phase-shift eying (PSK), M-ary quadrature-amplitude modulation (QAM) and OFDM 10 h xx h xy θ=π/2, ε=0 08 Magnitude Normalized frequency, f τ (a) h h xy xx θ=π/3, ε= Magnitude Normalized frequency, f τ (b) Fig 1 Magnitude response of h xx and h xy Jones matrix coefficients against the normalized frequency for: (a) θ=π/2 and ε=0, and (b) θ=π/3 and ε=0 The polarization interference cancellation scheme based on V-BLAST algorithm, which uses an LDPC code as channel code, is shown in Fig 2 The bit streams originating from m different information sources are encoded using different (n, i ) LDPC codes of code rate r i = i /n i denotes the number of information bits of ith (i=1,2,,m) component LDPC code, and n denotes the codeword length, which is the same for all LDPC codes The use of different LDPC codes allows us to optimally allocate the code rates The bit-interleaved coded modulation (BICM) scheme can be considered as a special multilevel coding (MLC) scheme in which all of the component codes are of the same rate [1] The outputs of m LDPC encoders are written row-wise into a bloc-interleaver bloc The mapper accepts m bits at (C) 2008 OSA 15 September 2008 / Vol 16, No 19 / OPTICS EXPRESS 14847

4 time instance i from the (mxn) interleaver column-wise and determines the corresponding M- ary (M=2 m ) signal constellation point (φ I,i, φ Q,i ) in a two-dimensional (2D) constellation diagram such as M-ary PSK or M-ary QAM (The coordinates correspond to in-phase and quadrature components of M-ary 2D constellation) The 2D signal constellation points are split into two streams for OFDM transmitters (see Fig 2(b)) corresponding to the x- and y- polarizations The QAM constellation points are considered to be the values of the fast Fourier transform (FFT) of a multi-carrier OFDM signal The OFDM symbol is generated as follows: N QAM input QAM symbols are zero-padded to obtain N FFT input samples for inverse FFT (IFFT), N G non-zero samples are inserted to create the guard interval, and the OFDM symbol is multiplied by the window function For efficient chromatic dispersion and PMD compensation, the length of cyclically extended guard interval should be linger than the total spread due to chromatic dispersion and DGD The cyclic extension is accomplished by repeating the last N G /2 samples of the effective OFDM symbol part (N FFT samples) as a prefix, and repeating the first N G /2 samples as a suffix After D/A conversion (DAC), the RF OFDM signal is converted into the optical domain using the dual-drive Mach-Zehnder modulator (MZM) Two MZMs are needed, one for each polarization The outputs of MZMs are combined using the polarization beam combiner (PBC) One DFB laser is used as a CW source, with x- and y-polarization separated by the polarization beam splitter (PBS) Source channels 1 LDPC encoder r 1 = 1 /n m LDPC encoder r m = m /n m Interleaver mxn Mapper OFDM transmitters s OFDM,x DFB s OFDM,y PBS MZM MZM to fiber PBC (a) QAM symbols S/P converter and Subcarrier mapper IFFT Cyclic extension insertion DAC DAC LPF LPF I Q (b) From SMF From local laser PBS PBS Coherent detector Coherent detector OFDM receivers + BLAST detector APP Demapper Bit LLRs Calculation LDPC Decoder 1 1 LDPC Decoder m m (c) Coherent detector π/2 v I v Q (d) ADC ADC LPF LPF FFT Symbol detection by BLAST detector P/S converter (e) Fig 2 The architecture of polarization interference cancelation scheme in combination with LDPC-coded OFDM: (a) transmitter architecture, (b) OFDM transmitter configuration, (c) receiver architecture, (d) coherent detector configuration, and (e) OFDM receiver configuration DFB: distributed feedbac laser, PBS(C): polarization beam splitter (combiner), MZM: dual-drive Mach-Zehnder modulator, APP: a posteriory probability, LLRs: loglielihood ratios (C) 2008 OSA 15 September 2008 / Vol 16, No 19 / OPTICS EXPRESS 14848

5 r i, C C r i, + - r = C r D s i, i, i, D s i, detector (a) r ( l ) i, + ( l + 1 ) ( l ) ( l ) r =r E s i, i, i, - E =C H -diag(c H ) s ( l ) i, detector (b) Fig 3 The configurations of polarization interference cancelation schemes: (a) BLAST-type polarization interference cancelation scheme, and (b) iterative polarization cancelation scheme On the receiver side, in PT-coded OFDM, we have the option to use only one polarization or to use both polarizations The polarization diversity OFDM [3], and polarization interference cancellation schemes proposed here require the use of both polarizations The receiver architecture employing both polarizations is shown in Fig 2(c) The configuration of polarization interference cancellation scheme by a BLAST-algorithm is shown in Fig 3(a) The received symbol vector in th subcarrier of ith OFDM symbol in both polarization (see Eq (2)) is linearly processed, the processing is described by matrix C related to channel matrix H as shown below, and the estimate of polarization interference obtained from preliminary decisions s i,, denoted as Ds i,, is removed from received symbol r The Euclidean detector can be used to create the preliminary decisions When the i, presence of ASE noise is ignored, the zero-forcing V-BLAST polarization interference cancellation scheme results The matrices C and D can be determined from QR-factorization of channel matrix H =Q R, as follows ( ) ( ) C R Q D R R I 1 1 = diag, = diag, (3) where I is the identity matrix, and with diag() we denoted the diagonal elements of R Notice that elements at the main diagonal in D are zero in order to have only polarization interference be removed (We use to denote the simultaneous transposition and complexconjugation) In the presence of ASE noise, the matrices C and D can be determined by minimizing the MSE, which leads to 1 ( )( ) ( ) C S S H D S S I 1 1 = diag, = diag, (4) where S is the upper triangular matrix obtained by Cholesy factorization of H H + I / SNR = SS, where SNR denotes corresponding electrical SNR, and I is the identity matrix The derivation of (3),(4) is equivalent to that for wireless communications [8], and as such is omitted here Because the ZF V-BLAST is derived by ignoring the influence of ASE noise, we propose to use the ZF V-BLAST as starting point, and perform the polarization interference cancellation in an iterative fashion as shown in Fig 3(b) If () l r denotes the i, processed received symbol of th subcarrier in ith OFDM symbol (for both polarizations) in lth iteration, then corresponding received symbol in (l+1)th iteration can be found by (C) 2008 OSA 15 September 2008 / Vol 16, No 19 / OPTICS EXPRESS 14849

6 ( l+ 1) ( l) [ ] () l r = r C H -diag( C H ) s, (5) i, i, i, where () l s denotes the transmitted symbol (of th subcarrier in ith OFDM symbol (for both i, polarizations)) estimate in lth iteration The matrices C and D are already introduced in (3) Notice that different matrix operations applied in (3)-(5) are trivial because the dimensionality of matrices is small, 2x2 The BLAST-detector soft estimates of symbols carried by th subcarrier in ith OFDM symbol, s i, xy, ( ), are forwarded to the a posteriori probability (APP) demapper, which determines the symbol log-lielihood ratios (LLRs) λ x(y) (q) (q=0,1,,2 b -1) of x- (y-) polarization by λ ( )( q) x y = ( Re s ( ) Re QAM,, ( map ix y ( q) ) ) 2 2σ ( Im s ix,, ( y) Im QAM ( map ( q) ) ) b ; q = 0,1,, σ where Re[] and Im[] denote the real and imaginary part of a complex number, QAM denotes the QAM-constellation diagram, σ 2 denotes the variance of an equivalent Gaussian noise process originating from ASE noise, and map(q) denotes a corresponding mapping rule (Gray mapping is applied here) (b denotes the number of bits per constellation point) Let us denote by v j,x(y) the jth bit in an observed symbol q binary representation v=(v 1,v 2,,v b ) for x- (y-) polarization The bit LLRs needed for LDPC decoding are calculated from symbol LLRs by : 0 exp qv λx( y)( q) j = L vˆ log (7) jx, ( y) = qv : 1 exp λx( y)( q) j = Therefore, the jth bit reliability is calculated as the logarithm of the ratio of a probability that v j =0 and probability that v j =1 In the nominator, the summation is done over all symbols q having 0 at the position j, while in the denominator over all symbols q having 1 at the position j The extrinsic LLRs are iterated bacward and forward until convergence or pre-determined number of iterations has been reached The LDPC code used in this paper belongs to the class of quasi-cyclic (array) codes of large girth (g 10) [11], so that the corresponding decoder complexity is low compared to random LDPC codes, and do not exhibit the error floor phenomena in the region of interest in fiber-optics communications ( ) 3 Evaluation of proposed PMD compensation schemes We are turning our attention to the BER performance evaluation of the proposed schemes The results of simulation for uncoded OFDM for different PMD compensation schemes are shown in Fig 4 The OFDM system parameters are chosen as follows: N QAM =512, oversampling is two times, OFDM signal bandwidth is set to 10 GHz, and N G =256 samples The MMSE V-BLAST and iterative polarization cancellation schemes (with ZF V-BLAST as starting point) perform identically (only MMSE curve is shown because the curves overlap each other), while ZF V-BLAST is slightly worse Polarization diversity OFDM outperforms MMSE V-BLAST at low BERs, but performs comparable at BERs above 10-2, which is the threshold region of girth-10 LDPC codes employed here Moreover, the spectral efficiency of MMSE V-BLAST is twice higher because in polarization diversity OFDM the same symbol is transmitted twice over both polarizations The MMSE V-BLAST OFDM outperforms the PTcoding based OFDM at both low and high BERs, and has two times higher spectral efficiency The MMSE V-BLAST OFDM scheme is able to compensate even 1200 ps of DGD with negligible penalty Notice that for corresponding turbo equalization [2] or maximum- 2 2 (6) (C) 2008 OSA 15 September 2008 / Vol 16, No 19 / OPTICS EXPRESS 14850

7 lielihood sequence estimation schemes, the detector complexity grows exponentially as DGD increases, and normalized DGD of 800 ps it would require the trellis description (see [2]) with 2 17 states, which is too high for practical implementation Our schemes, although of lower complexity, are able to compensate up to 1200 ps of DGD with negligible penalty The proposed schemes also outperform the scheme implemented by Nortel Networs researchers [4], capable of compensating the rapidly varying first order PMD with pea DGD of 150 ps Uncoded OFDM: M=4, Pol Div, τ= 800 ps M=4, ZF V-BLAST, τ=800 ps M=4, MMSE V-BLAST, τ=400 ps, 800 ps, 1200 ps M=4, MMSE V-BLAST, B2B M=2, MMSE V-BLAST, B2B and τ=800 ps M=2, PT-coding, τ=800 ps 10-1 Bit-error ratio, BER Optical SNR, OSNR [db/01 nm] Fig 4 BER performance of proposed schemes against polarization diversity OFDM, and PTcoding based OFDM B2B: bac-to-bac The results of simulations for LDPC-coded OFDM when MMSE V-BLAST polarization cancellation scheme is used are shown in Fig 5 The girth-10 LDPC(16935,13550) code of rate 08 and column weight 3 is used in simulations This code does not exhibit error floor phenomena for the region of interest in optical communications (see [11]) At BER of 10-9 the LDPC-coded OFDM with MMSE V-BLAST polarization interference cancellation scheme outperforms PT-coding based OFDM by about 1 db, and has the spectral efficiency twice higher In simulations shown in Figs 4-5 the average launch symbol power was set to -3 dbm (C) 2008 OSA 15 September 2008 / Vol 16, No 19 / OPTICS EXPRESS 14851

8 10-1 LDPC(16935,13550)-coded OFDM: M=2, PT-coding, τ = 800 ps M=2, MMSE V-BLAST, τ = 800 ps M=4, MMSE V-BLAST, τ = 400 ps, 800 ps, 1200 ps Bit-error ratio, BER M=2 M=4 Uncoded MMSE V-BLAST τ= 800 ps Optical SNR, OSNR [db/01 nm] Fig 5 BER performance of proposed polarization interference cancelation schemes for LDPCcoded OFDM 5 Conclusion We proposed two alternative PMD compensation schemes to turbo equalization, channel equalization scheme from [4], polarization diversity OFDM [3], and PT-coding based OFDM The proposed schemes are suitable for use in multilevel (M 2) bloc-coded modulation schemes with coherent detection In contrast to the PMD turbo equalization scheme whose complexity grows exponentially as DGD increases, the complexity of the proposed schemes stays the same The spectral efficiency of proposed schemes is two times higher than that of polarization diversity OFDM and PT-coding based OFDM The first scheme is based on MMSE V-BLAST algorithm, used in MIMO wireless communications to deal with spatial interference The second scheme is based on iterative polarization interference cancellation Those two schemes perform comparable, and are able to compensate up to 1200 ps of DGD with negligible penalty When used in combination with girth-10 LDPC codes, those schemes outperform PT-coding based OFDM by 1 db at BER of 10-9 Acnowledgments This wor was supported in part by the National Science Foundation (NSF) under Grant IHCS (C) 2008 OSA 15 September 2008 / Vol 16, No 19 / OPTICS EXPRESS 14852