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

Analysis and compensation of dispersion-induced bit loss in a photonic A/D converter using time-wavelength interweaved sampling clock

Open Access Open Access

Abstract

In this paper, the timing jitter induced by the fiber dispersion in photonic A/D converters using time-wavelength interweaved sampling clocks generated by optical time-division-multiplexing (OTDM) with fiber delay lines is analyzed and effective bit loss is calculated. A compensation method is proposed to decrease the dispersion-induced jitter. Simulations are performed and the results show the validity of the proposed compensation method. An experimental demonstration is carried out to verify the theoretical expression derived.

©2009 Optical Society of America

1. Introduction

In order to generate time-wavelength interweaved sampling clock of tens even hundreds GHz for photonic analog-to-digital converter (ADC), a laser source with wide spectrum and ultra-short pulse output is needed [13]. Passive mode lock femto-second fiber laser is one of the appropriate candidates. The spectrum of a commercially available femto-second fiber laser may cover 50nm in the 1550 nm optical communications band. Such a spectrum can be sliced into 20-30 wavelengths. The wavelength-division-multiplexing (WDM) technique is desirable not only for increasing the pulse rate, but also for providing the possibility of parallel signal processing [1, 2]. However, the repetition rate of such lasers is usually very low. Therefore, optical time-division-multiplexer (OTDM) should be used to further aggrandize the pulse rate (i.e. the sampling rate). It is relatively convenient and cheap to use fiber delay lines (FDLs) to build the OTDM [2]. The drawback is that the fiber dispersion may aggravate timing jitter, which is one of the factors [4] that limit the effective bits of AD converters. In this paper, we analyze the effect of fiber dispersion with simulation verifications and propose a compensation scheme.

2. Theory

Figure 1 is the principle diagram that generates ultrahigh-speed time-wavelength interweaved sampling clock from a passive mode lock femto-second fiber laser (MLFL) with repetition rate of RS. The WDM part has M wavelength channels, which are arranged at different time slot so as to form a pulse group by using M fiber delay lines. The OTDM part has N paths, which replicate the pulse group by N times. The repetition rate of the output pulse train is thus MNRS.

 figure: Fig. 1

Fig. 1 The principle diagram generating time-wavelength interweaved sampling clock (FRM: Faraday Rotator Mirror)

Download Full Size | PDF

In such a configuration, timing jitter, amplitude jitter and pulse width variation are three main factors that limit the effective bits [4, 5]. According to our experiment observation, both the WDM part and the OTDM part will induce power fluctuation and delay variation between different channels. In the WDM part, each channel uses a separated fiber. So the power fluctuation can be balanced through variable optical attenuator and the delay difference can be eliminated through careful adjustment. In the OTDM part, however, the fiber delay lines are shared by different channels. Though the amplitude fluctuation can be eliminated by using highly uniform couplers, the timing error caused by dispersion in the OTDM paths is hard to compensate.

Suppose the wavelength interval is Δλ and the fiber length difference between the adjacent delay units isΔL, then the jth wavelength and its delay, δj, in the WDM part, and the fiber length of ith delay unit, Li, can be expressed as:

λj=λ1+(j1)Δλ,1jM
δj=δ1+(j1)MNRs,1jM
Li=L1+(i1)ΔL,1iN
where λ1 is the lowest wavelength, δ1 is the shortest delay of the WDM, and L1 is the length of the shortest path of the OTDM, respectively. The ideal delay caused by ΔL should be 1/NRS and the ideal time delay for the wavelength j through path i should be:

t¯ij=iNRs+δj,1iN,1jM

However, the group velocity, v(λj), j=1,2…M, differs from each other due to fiber dispersion. The actual time delay of the wavelength j through path i is:

tij=Liv(λj)+δj',1iN,1jM
where δj' is the actual fiber delay of the jth channel. Suppose there is a wavelength, λs (1≤s≤M), whose delay after transmission through the fiber section of ΔL is 1/NRS (in fact, what we concern is the relative delay, so the assumption is valid which will be seen in Section 3). We regard this channel as the ideal channel and rewrite the jth wavelength using λs as:

λj=λs+(js)Δλ,1jM,1sM

Therefore the delay difference between λjandλsfor the fiber section of ΔL is given by [6]:

Δτ=ΔLv(λj)ΔLv(λs)ΔLD(λ(j+s)/2)Δλ(js)
where D(λ(j+s)/2) is the dispersion at the center wavelength between λj and λs. So the delay deviation, i.e. the timing jitter, of the jth wavelength through ith path, τij, is given by:

τij=tijtij¯LiD(λ(j+s)/2)Δλ(js)+(δj'δj)=[L1+(i1)ΔL]D(λ(j+s)/2)Δλ(js)+(δj'δj),1iN,1jM

Equation (8) shows the sampling time error of the pulse Pij relative to the ideal case, which is caused by fiber dispersion in the WDM part and the OTDM part. As mentioned above, the timing error in the WDM part can be emendated by careful adjustment of the related fiber length. Therefore we will ignore the (δjδj') term in the following analysis. According to the analysis by G. C. Valley, etc., the effective bits of photonic AD converters limited by the timing jitter is given by [4,5, 7]:

Neff=log2(1/3πfσj)
where f is the highest frequency of the signal sampled [4]. In Eq. (9), σj is the RMS of the timing jitter. It contains both the laser source jitter and jitter induced by other factors such as dispersion-induced jitter given by Eq. (8). Taking into consideration the fact that τijis periodic and can be equivalently treated as discrete uniform distribution, one can calculate the RMS of the dispersion induced timing jitter:

σjd=i=1Nj=1M|τij|2MN

Since the laser source timing jitter and dispersion induced timing jitter are independent, one can get the total equivalent timing jitter as follows:

σj=(σjd)2+(σjs)2
where σjs is the RMS timing jitter of the laser source.

3. Compensation method and analysis

We simulated the impact of dispersion induced timing jitter on the effective bits. The parameters used are as follows. The output pulse width of the MLFL is 70fs with optical spectrum of 1560±22.5nm, repetition period of 27.44ns and σjs=150fs. The WDM part has 16 channels with channel spacing of 200GHz (1.6nm) and λ1=1540.56nm. The OTDM part has 64 paths with ΔL≈8.575cm and L1=5m. The highest frequency of the sampled microwave signal is f=10GHz. The dispersion coefficient, D(λ), is given by [8]:

D(λ)=2πcλ2β2
where c is the speed of light in vacuum, β220ps2/km in the 1550nm range.

First we calculated the equivalent timing jitter σjdcaused by dispersion for different λs as shown in Fig. 2 . Obviously, the dispersion induced timing jitter is the lowest when s=M/2. It means that when we make an OTDM, the “ideal wavelength” should be the central wavelength of the spectrum used. Even though, the jitter caused by dispersion is still large and can cause significant effective bit loss. So it must be compensated. It is hard to compensate the dispersion for each wavelength in OTDM because OTDM paths are shared. One method to alleviate the dispersion effect is to use low dispersion fiber, such as dispersion shifted fiber (DSF) or dispersion flattened fiber, provided that devices such as optical couplers based on such fiber with precise coupling ration and ultra flatness over the interested bandwidth are available. The other way to decrease the jitter is to compensate the dispersion in the WDM part.

 figure: Fig. 2

Fig. 2 Dispersion-induced timing jitter

Download Full Size | PDF

Set s=M/2, then stretch the jth (j=1,2,…,M) fiber in WDM part so that it produces an extra compensation Cjtgiven by:

Cjt=CCj
where C is a constant. The relative delay to the “ideal wavelength channel”, λs, may be positive or negative. It is more convenient to produce extra delay by stretching the fiber. C is used to avoid the case that the total extra delay, Cjt, is negative. Because C is effectively equal to each other, we will not include it in the equations of the following analysis. The term “effectively equal” indicates that the dispersion induced delay for each wavelength, i.e., (δjδj'), is considered in C. After compensation, the sampling error becomes:

τij'=LiD(λj/2+M/4)Δλ(jM2)Cj

Different choice of Cjadds different effects onσjd. In this paper we use a simple strategy of compensation, i.e., Cj is given by:

Cj=LkD(λj/2+M/4)Δλ(jM2),1kN

Hence the setting of Cjbecomes the problem of choosing optimal k. Figure 3 shows the relationship between the effective bits and k when s=M/2. We can find that when k=N/2, the impact of the timing jitter is the lowest. Therefore, Eq. (15) gives the amount of compensation of the jth wavelength in the WDM part for the compensation strategy mentioned above, where Lk is given by Eq. (3) with k=N/2.

 figure: Fig. 3

Fig. 3 Relationship between the effective bits and k

Download Full Size | PDF

In order to prove the validity of the compensation scheme, we simulated the effective bits for different number of wavelengths. Figure 4 shows the relationship for the parameters given above. It can be clearly seen that the fiber dispersion induced timing jitter has strong influence on the effective bits in the system employing fiber delay lines to generate time-wavelength interweaved sampling clocks as shown in Fig. 1, especially when more wavelengths are used for generating ultrahigh rate sampling clocks. The simple compensation scheme we proposed can improve the performance of the photonic ADC by increasing 2-3 effective bits. We also simulated the extreme cases that the laser source is very stable (σjs=10fs) and very noisy (σjs =300fs). The results still show effectiveness of the compensation scheme.

 figure: Fig. 4

Fig. 4 Relationship between effective bits and number of channels

Download Full Size | PDF

Figure 5 shows the effective bits versus channel number when DSF is used to implement the OTDM. The DSF dispersion is modeled by formulas D(λ)=(λλ0)S0 recommended by EIA and ITU-T [6], where λ0 is the zero dispersion wavelength around 1550 nm and S0 is the dispersion slope. It can be seen that DSF can provide higher effective bits in comparison with the conventional single mode fiber. However, the performance can still be further improved by employing the compensation method proposed, especially when a stable laser source and more channels are used.

 figure: Fig. 5

Fig. 5 Relationship between effective bits and number of channels for DSF-based OTDM

Download Full Size | PDF

In addition to the theoretical calculation, the right Cj value can be found experimentally as follows. Firstly, measure the time interval between the input pulses of jth wavelength and the reference wavelength. Then measure the time interval again after the pulse train passes through the N/2th path of the OTDM as shown in Fig. 6 (a) . Cjcan be obtained by calculating the difference of the two time intervals.

 figure: Fig. 6

Fig. 6 Measurement of pulse intervals to determine Cj. (a) Output from the N/2th path of the OTDM before compensation. (b) Output from the N/2th path of the OTDM after compensation

Download Full Size | PDF

Figure 6 shows the measured results. In our experiment, the reference pulse is with wavelength of 1554.13nm and the two pulses used as example are at 1546.12nm and 1565.5 nm. Their intervals with the reference pulse before input into the OTDM are 66.4ps and 62.7ps, respectively. After passing through the OTDM, these values changed to 66.8ps and 63.2ps, respectively, as shown in Fig. 6(a). Hence C values are equal to −0.4ps and 0.5ps, respectively. From Eq. (15), one can calculate the theoretic values are −0.36ps and 0.51ps, respectively, which are in good agreement with the test result.

The compensation can be implemented by high precise fiber stretcher [9, 10] or fused micro drawing [11]. Figure 6(b) is the measured results after compensation using manually tunable mechanical stretchers we designed and implemented. It can work stably, as shown in Fig. 7 , with resolution sufficient for such applications.

 figure: Fig. 7

Fig. 7 The measured stability of the designed mechanical fiber stretcher

Download Full Size | PDF

4. Conclusion

In this paper, the timing jitter caused by dispersion in the shared OTDM fiber channels in a photonic A/D converter is analyzed. Expressions to calculate the equivalent RMS timing jitter are presented. A compensation strategy is proposed by simply stretching the fiber in the WDM part. The simulation results show that, by using the strategy, the performance of the photonic A/D converter can be effectively improved by 2-3 bits. Simulation shows that it is effective even if the OTDM is configured with dispersion shifted fiber. Experimental measurement verified the validity of the theoretical expressions derived.

Acknowledgments

This work was supported in part by 863 Project under Grant 2007AA01Z275 and 2006AA01Z242, National Science Foundation of China (NSFC) under Grant 90704002 and 60507013, 60877012, STCSM Project under Grant 09JC1408100 and Dawn Program for Excellent Scholars by the Shanghai Municipal Education Commission.

References and links

1. J. U. Kang and R. D. Esman, “Demonstration of time interweaved photonic four-channel WDM sampler for hybrid analogue-digital converter,” Electron. Lett. 35(1), 60–61 (1999). [CrossRef]  

2. G. Wu, M. Li, and B. Wang, at al., “An Approach to Generate Multi-wavelength Sampling Clock for Photonic A/D Converters,” The 14th OptoElectronics and Communications Conference (OECC 2009), Hong Kong, July 13–17, 2009.

3. T. R. Clark, J. U. Kang, and R. D. Esman, “Performance of a Time- and Wavelength-Interleaved Photonic Sampler for Analog–Digital Conversion,” IEEE Photon. Technol. Lett. 11(9), 1168–1170 (1999). [CrossRef]  

4. G. C. Valley, J. P. Hurrell, and G. A. Sefler, “Photonic analog-to-digital converters: fundamental and practical limits,” Proc. SPIE 5618, 96–106 (2004). [CrossRef]  

5. H. F. Taylor, “An optical analog-to-digital converter—design and analysis,” IEEE J. Quantum Electron. 15(4), 210–216 (1979). [CrossRef]  

6. Gerd Kerser, Optical Fiber Communications, (Third Edition). (MCGraw-Hill Companies, New York, 2000).

7. R. H. Walden, “Analog-to-digital converter survey and analysis,” IEEE J. Sel. Areas Comm. 17(4), 539–550 (1999). [CrossRef]  

8. P. G Agrawal, Nonlinear Fiber Optics (Third Edition). (Academic Press, San Diego, 2001).

9. www.generalphotonics.com, Fiber Stretcher-FST-001-B, General Photonics Corp.

10. www.optiphase.com, PZ1-Low-profile Fiber Stretcher, Optiphase, Inc.

11. X. Li, L. Peng, S. Wang, Y.-C. Kim, and J. Chen, “A novel kind of programmable 3(n) feed-forward optical fiber true delay line based on SOA,” Opt. Express 15(25), 16760–16766 (2007). [CrossRef]   [PubMed]  

Cited By

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

Alert me when this article is cited.


Figures (7)

Fig. 1
Fig. 1 The principle diagram generating time-wavelength interweaved sampling clock (FRM: Faraday Rotator Mirror)
Fig. 2
Fig. 2 Dispersion-induced timing jitter
Fig. 3
Fig. 3 Relationship between the effective bits and k
Fig. 4
Fig. 4 Relationship between effective bits and number of channels
Fig. 5
Fig. 5 Relationship between effective bits and number of channels for DSF-based OTDM
Fig. 6
Fig. 6 Measurement of pulse intervals to determine Cj . (a) Output from the N/2th path of the OTDM before compensation. (b) Output from the N/2th path of the OTDM after compensation
Fig. 7
Fig. 7 The measured stability of the designed mechanical fiber stretcher

Equations (15)

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

λj=λ1+(j1)Δλ,1jM
δj=δ1+(j1)MNRs,1jM
Li=L1+(i1)ΔL,1iN
t¯ij=iNRs+δj,1iN,1jM
tij=Liv(λj)+δj',1iN,1jM
λj=λs+(js)Δλ,1jM,1sM
Δτ=ΔLv(λj)ΔLv(λs)ΔLD(λ(j+s)/2)Δλ(js)
τij=tijtij¯LiD(λ(j+s)/2)Δλ(js)+(δj'δj)=[L1+(i1)ΔL]D(λ(j+s)/2)Δλ(js)+(δj'δj),1iN,1jM
Neff=log2(1/3πfσj)
σjd=i=1Nj=1M|τij|2MN
σj=(σjd)2+(σjs)2
D(λ)=2πcλ2β2
Cjt=CCj
τij'=LiD(λj/2+M/4)Δλ(jM2)Cj
Cj=LkD(λj/2+M/4)Δλ(jM2),1kN
Select as filters


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