Advanced
Joint Probability Density Functions for Direct-Detection Optical Receivers
Joint Probability Density Functions for Direct-Detection Optical Receivers
Journal of the Optical Society of Korea. 2014. Apr, 18(2): 124-128
Copyright © 2014, Journal of the Optical Society of Korea
  • Received : February 02, 2014
  • Accepted : March 03, 2014
  • Published : April 25, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Share
Article
Author
Metrics
Cited by
TagCloud
About the Authors
Jae Seung Lee
jslee@kw.ac.kr
Abstract
We derive joint probability density functions (JPDFs) for two adjacent data from direct-detection optical receivers in dense wavelength-division multiplexing systems. We show that the decision using two data per bit can increase the receiver sensitivity compared with the conventional decision. Our JPDFs can be used for software-defined optical receivers enhancing the receiver sensitivities for intensity-modulated channels.
Keywords
I. INTRODUCTION
Recently, software-defined networks that can provide practical ways of handling multi-standard environments have been investigated extensively [1 , 2] . In a similar way, software-defined optical receivers are under development in optical fiber communication systems to satisfy many different kinds of modulation formats and baud rates [3 , 4] . Although coherent optical systems have become more practical than during their early stages, currently deployed optical transmission systems are mostly the intensity-modulated and direct-detection (IM/DD) type.
Conventional optical receivers for IM/DD systems use a single data per single bit for the decision [5] and require clock recovery circuits. If we could use multiple data for the decision, we do not need the clock recovery circuits [4] while enhancing the receiver sensitivities. However, there are no analyses, to our knowledge, about the correlations between two adjacent data from an optical receiver.
In this paper, we propose to use correlations of two adjacent data for the decision in direct-detection optical receivers. Using the receiver eigenmodes [6 - 12] , we derive the joint characteristic function (JCF) for two successive data from an optical receiver and evaluate the corresponding joint probability density function (JPDF) [13] . The receiver eigenmodes can describe accurately the effects of the amplified spontaneous emission (ASE), received optical waveforms, and shapes of optical and electrical filters within the receiver [7] . Recently, receiver eigenmode contributions have been analyzed as a function of time for the optical receiver output [12] . It has been found that, in conventional dense WDM systems [14 - 17] , where the channel spacing is comparable to the bit rate, the lowest-order (0-th) receiver eigenmode contributes dominantly. We will use this fact to find the correlations of two adjacent data and the threshold line for the decision to get higher sensitivities than conventional receivers.
II. JOINT CHARACTERISTIC FUNCTION
Before the derivation of the JPDF, we derive the JCF first. We will consider two adjacent samples at times t 1 and t 2 . Considering only the polarization component parallel to the received optical signal, we can find the voltage at t 1 , denoted as y 1 , as [7]
PPT Slide
Lager Image
where Sm ( t 1 ) and Nm ( t 1 ) are complex numbers that represent the signal and the noise amplitudes, respectively, for the m -th receiver eigenmode. λm is the m -th eigenvalue and k is a proportional constant. For all m , the real and the imaginary parts of Nm ( t 1 ) are mutually independent zero mean Gaussian random variables having an identical variance of
PPT Slide
Lager Image
In conventional dense WDM systems, the lowest-order receiver eigenmode contributes dominantly [12] and we can approximate
PPT Slide
Lager Image
We normalize y 1 such that it becomes a power ratio such that
PPT Slide
Lager Image
where the denominator is the noise power per receiver eigenmode per polarization. Similarly, we have the normalized voltage at t 2 as
PPT Slide
Lager Image
N 0 ( t 1 ) and N 0 ( t 2 ) are related by a correlation function as [12]
PPT Slide
Lager Image
The correlation function, C ( t 1 t 2 ), relates the 0 -th order eigenmode amplitudes at different times, which is given by
PPT Slide
Lager Image
where ø 0 ( ω ) is the 0 -th order eigenfunction. Let’s denote the real part of N 0 ( t 1 ) and N 0 ( t 2 ) as N 1r and N 2r , respectively, and the imaginary part of N 0 ( t 1 ) and N 0 ( t 2 ) as N 1i and N 2i , respectively. The covariance matrix for N 1r , N 1i , N 2r , and N 2i can be found as [13]
PPT Slide
Lager Image
where Cr and Ci are real and imaginary parts of C ( t 1 t 2 ), respectively. The inverse matrix of Σ is
PPT Slide
Lager Image
where
PPT Slide
Lager Image
For the Gaussian vector, X = (N 1r , N 1i , N 2r , N 2i ) T , where T means the transpose, its JPDF is [13]
PPT Slide
Lager Image
Thus the JPDF for N 1r , N 1i , N 2r , and N 2i can be written as
PPT Slide
Lager Image
The JCF for
PPT Slide
Lager Image
and
PPT Slide
Lager Image
can be written as
PPT Slide
Lager Image
PPT Slide
Lager Image
where E {⋅} is the ensemble average. The integrations can be done exactly using the Gaussian integration formula,
PPT Slide
Lager Image
where α and β are constants with Re α > 0 . Thus we find
PPT Slide
Lager Image
where
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
PPT Slide
Lager Image
S 1r and S 1i are real and imaginary parts of S 0 ( t 1 ), respectively. Also, S 2r and S 2i are real and imaginary parts of S 0 ( t 2 ), respectively. If we include the other polarization, where only the ASE components are present, we obtain the JCF exactly as follows:
PPT Slide
Lager Image
III. JOINT PROBABILITY DENSITY FUNCTION
The JPDF of
PPT Slide
Lager Image
and
PPT Slide
Lager Image
can be found from the twodimensional Fourier transform
PPT Slide
Lager Image
When the signal is absent, the JCF becomes
PPT Slide
Lager Image
which gives the JPDF of
PPT Slide
Lager Image
and
PPT Slide
Lager Image
in an exact form
PPT Slide
Lager Image
where I 1 is the modified Bessel function. We have used the following relation ( x , M > 0) :
PPT Slide
Lager Image
If the two data are independent, we find from (22)
PPT Slide
Lager Image
This JPDF is just a product of each sample’s PDF. To find (24), we have increased G satisfying the G < 1 condition. When the signal is present, the JPDF can be obtained numerically. Note that the integration over ξ 1 can be done exactly to yield
PPT Slide
Lager Image
where
PPT Slide
Lager Image
The remaining integral can be done using a fast-Fourier-transform algorithm.
We may use the asymptotic form for the Bessel function,
PPT Slide
Lager Image
, and evaluate the integral of (25) using the method of steepest descents [18], which gives
PPT Slide
Lager Image
The complicated expression of (27) is valid when the correlation is low ( G ≈1). We perform numerical evaluations of (25) to find BER throughout.
IV. RESULTS AND DISCUSSION
For our analysis, we use a Gaussian optical receiver [11 , 12] , where both optical and electrical filters are Gaussian. We assume A 1 = A 2 . We choose the 3-dB bandwidth of the optical filter two times larger that of the electrical filter. The 3-dB bandwidth of the electrical filter is 0.7 times the bit rate. Then the correlation function (6) becomes C (τ) = exp[−(0.7π) 2 2 q /{(1− q 2 )ln2}(τ/ T ) 2 ] [12] . We have q = 0.268 according to the bandwidth ratio between the optical and the electrical filters. T is the bit period. The two sampling points are T /4 apart and we obtain G = 0.396. We assume no phase changes of the signal between two sampling points, which gives
PPT Slide
Lager Image
Figure 1 shows two JPDFs evaluated numerically along the
PPT Slide
Lager Image
=
PPT Slide
Lager Image
line for A 1 = 40.65 and for A 1 = 0 simultaneously. The former corresponds to the mark and the latter corresponds to the space. The JPDF for the space has been scaled down by the factor of 10 compared with the JPDF for the mark.
PPT Slide
Lager Image
Joint probability density functions (JPDFs) evaluated numerically along the = line for the mark (A1,2= 40.65) and for the space (A1,2 = 0). The JPDF for the space has been scaled down by the factor of 10 compared with the JPDF for the mark. The two sampling points are T/4 apart, where T is the bit period. A Gaussian optical receiver is used with G = 0.396.
In Fig. 2 , we show the foregoing JPDFs in a 3-dimensional way. It has been obtained by adding the JPDFs’ functional values to
PPT Slide
Lager Image
after multiplying a constant 4×103. Only the points, where at least one of the JPDFs is larger than 10-10, are shown. The JPDF for the space has been scaled down by the factor of 10 here also. Since A 1 = A 2 , the JPDFs, including (22), are symmetric with respect to the =
PPT Slide
Lager Image
=
PPT Slide
Lager Image
line. Note that f (
PPT Slide
Lager Image
,
PPT Slide
Lager Image
) = 0 when at least one of
PPT Slide
Lager Image
and
PPT Slide
Lager Image
is zero. Along the
PPT Slide
Lager Image
=
PPT Slide
Lager Image
line, the width of each JPDF increases as
PPT Slide
Lager Image
increases. In conventional optical receivers, where a threshold voltage is used, we have BER = (2 A 1 /π) 1/4 exp(− A 1 /4)/4 [11] , which is 2.2×10 -5 . In our case, there is a threshold line as is shown in Fig. 2 where both JPDFs have the same value. The area of each JPDF beyond that line contributes to the BER [5] which is evaluated numerically as 6.8×10 -6 . This corresponds to about 0.5 dB enhancement in A 1 compared with the conventional BER of 2.2×10 -5 . A 1 and A 2 are actually signal-to-noise ratios [11] . Thus we can have higher system margins using the correlation between the data from the direct-detection optical receiver.
PPT Slide
Lager Image
Plot of the JPDFs of Fig. 1 in the (, ) plane. We have added the JPDF’s functional value to after multiplying a constant 4×103. We show only the points where at least one of JPDFs is larger than 10-10. The JPDF for the space has been scaled down by the factor of 10.
If the two sampling points are T /2 apart with all other parameters fixed, G increases to 0.866. The JPDFs in this case are shown in Fig. 3 , which become more circularly symmetric owing to the reduced correlation. The BER is 2.6×10 -7 and the enhancement in A 1 increases to about 1.6 dB. If the two sampling points are 3 T /4 apart with all other parameters fixed, G increases to 0.989. The BER is 1.3×10 -8 and the enhancement in A 1 is about 2.4 dB. If the received signal has a different phase between two sampled points, this will decrease H in (18). In this case, our numerical analyses show that the BER decreases as the phase difference increases. Thus our BER values can be reduced introducing the phase change within a bit.
PPT Slide
Lager Image
Plot of the JPDFs in (, ) plane, where the two sampling points are T/2 apart with all other parameters the same as Fig. 2. G = 0.866. The JPDFs become more circularly symmetric than those of Fig. 2 owing to the reduced correlation.
In order to use the threshold line for the decision, we need analog-to-digital converters (ADCs) and digital signal processing (DSP) circuits instead of conventional D-flip-flop type decision. The speed of ADC and DSP circuits has been increased remarkably up to coherent 200 Gb/s per channel during recent years [19] . Realizations of software-defined optical receivers using these devices are now technically feasible and our JPDFs can be used to upgrade their capabilities for IM/DD channels.
V. CONCLUSION
We have derived the JPDFs for two voltage data from direct-detection optical receivers in dense WDM systems. We have shown that, with our JPDFs, we can reduce the BER and enhance the receiver sensitivities by 0.5 dB ~ 2.4 dB. Our decision method can be used for software-defined optical receivers supporting both coherent and intensity-modulated channels simultaneously.
Acknowledgements
The work reported in this paper was conducted during the sabbatical year of Kwangwoon University in 2012 and this research was also supported partly by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science, and Technology (2013R1A1A2012918).
References
Channegowda M. , Nejabati R. , Simeonidou D. 2013 “Softwaredefined optical networks technology and infrastructure: Enabling software-defined optical network operations [invited],” IEEE/OSA Journal of Optical Communications and Networking 5 A274 - A282
Cui Y. , Xiao S. , Liao C. , Stojmenovic I. , Li M. 2013 “Data centers as software defined networks: Traffic redundancy elimination with wireless cards at routers,” IEEE Selected Areas in Communications 31 2658 - 2672
Ip E. , Ji P. , Mateo E. , Huang Y.-K. , Xu L. , Qian D. , Bai N. , Wang T. 2012 “100G and beyond transmission technologies for evolving optical networks and relevant physical-layer issues,” Proc. IEEE 100 1065 - 1078
Isautier P. , Stark A. , Mehta K. , de Salvo R. , Ralph S.E. 2013 “Autonomous software-defined coherent optical receivers,” Anaheim Convention Center Anaheim, California, USA in Proc. OFC OTh3B.4-1 - OTh3B.4-3
Personick S.D. 1973 “Receiver design for digital fiber optic communications systems, I,” Bell Syst. Tech. J. 52 843 - 874
Kac M. , Siegert A.J.F. 1947 “On the theory of noise in radio receivers with square law detectors,” J. Appl. Phys. 18 383 - 397
Lee J.S. , Shim C.S. 1994 “Bit-error-rate analysis of optically preamplified receivers using an eigenfunction expansion method in optical frequency domain,” J. Lightwave Technol. 12 1224 - 1229
Forestieri E. , Secondini M. 2009 “On the error probability evaluation in lightwave systems with optical amplification,” J. Lightwave Technol. 27 706 - 717
Winzer P.J. , Pfennigbauer M. , Strasser M.M. , Leeb W.R. 2001 “Optimum filter bandwidths for optically preamplified NRZ and RZ receivers,” J. Lightwave Technol. 19 1263 - 1273
Holzlohner R. , Grigoryan V.S. , Menyuk C.R. , Kath W.L. 2002 “Accurate calculation of eye diagrams and bit error rates in optical transmission systems using linearization,” J. Lightwave Technol. 20 389 - 400
Lee J.S. , Willner A.E. 2013 “Analysis of Gaussian optical receivers,” J. Lightwave Technol. 31 2987 - 2993
Seo K.H. , Lee J.S. , Willner A.E. 2013 “Time-dependent analysis of optical receivers using receiver eigenmodes,” J. Opt. Soc. Korea 17 305 - 311
Papoulis A. , Pillai S.U. 2002 Probability, Random Variables, and Stochastic Processes 4th ed. McGraw-Hill New York, USA
Hurh Y.S. , Hwang G.S. , Jeon J.Y. , Lee K.G. , Shin K.W. , Lee S.S. , Yi K.Y. , Lee J.S. 2005 “1-Tb/s (100×12.4 Gb/s) transmission of 12.5-GHz-spaced ultradense WDM channels over a standard single-mode fiber of 1200 km,” IEEE Photon. Technol. Lett. 17 696 - 698
Kim S.Y. , Seo K.H. , Lee J.S. 2012 “Spectral efficiencies of channel-interleaved bidirectional and unidirectional ultradense WDM for metro applications,” J. Lightwave Technol. 30 229 - 233
Gnauck A. , Tkach R. , Chraplyvy A. , Li T. 2008 “High-capacity optical transmission systems,” J. Lightwave Technol. 26 1032 - 1045
Batsuren B. , Kim H.H. , Eom C.Y. , Choi J.J. , Lee J.S. 2013 “Optical VSB filtering of 12.5-GHz spaced 64 × 12.4 Gb/s WDM channels,” J. Opt. Soc. Korea 17 63 - 67
Arfken G.B. , Weber H.J. 2005 Mathematical Methods for Physicists 6th ed. Elsevier Academic New York, USA
Laperle C. , O’Sullivan M. 2014 “Advances in high-speed DACs, ADCs, and DSP for optical coherent transceivers,” J. Lightwave Technol. 32 629 - 643