A new method is proposed and analyzed for measuring the timing jitter of the transmitted pulse relative to the reference pulse using two type II phase-matched nonlinear crystals for second harmonic generation (SHG). The polarizations of the two pulses are exchanged in two crystals and the difference between two detected second harmonic signals can reflect the transmitted jitter. This new method provides a high sensitivity and timing resolution compared with the conventional RF (radio frequency) method. Since the overlapping levels in the two crystals are the same, the final output is zero when there is no time delay between the two pulses. Thus no offset is necessary to be subtracted from the final output and no time delay adjustment is required between the two pulses, compared with the previous optical method using one crystal and two dichroic beamsplitters. The jitter measuring performance is studied theoretically using non-stationary nonlinear wave-coupled equations for type II SHG of two pulses. The theoretical computation and analysis show that the sensitivity and the dynamic range of this new method depend on pulse width, crystal pulses and group velocity difference between two fundamental pulses.
©2009 Optical Society of America
The reliable, accurate and precise measurement for timing jitter between the signal after transmitting through an optical link and the original reference is necessary in many areas such as frequency transfer and timing distribution [1–6]. In the application of frequency transfer and timing distribution through fiber link, a local reference is transmitted to the remote end by the fiber link. However, the performance of the transmitted frequency/timing signal at the remote end is degraded by the fiber-induced phase noise (timing jitter). In order to achieve a reliable and precise frequency signal at the remote end, an active compensation loop is used to cancel the phase noise. In this way, the timing jitter between the transmitted signal and the reference must be measured for the action of the loop.
In the case of the frequency comb as a reference , the timing jitter of transmitted pulse relative to the reference pulse is required. One conventional RF method of jitter measuring uses two photodetectors and a mixer to obtain the phase noise and timing jitter [3–6]. Recently, a balanced optical cross correlator based on a type II SHG crystal and two dichroic beamsplitters is used for jitter measurement in timing distribution system . The sensitivity of this optical method is higher than that of the conventional one. However, since the output of this method is not zero when the two pulses are overlapped entirely (no time delay between them) before entering into the crystal, the transmitted pulse is required to be adjusted to have a fixed time delay relative to the reference before measuring, or an offset signal should be subtracted from the output.
In this paper, a new method is proposed for measuring the jitter between two pulses (one is the transmitted pulse, while the other is the original pulse as a reference) using two type II phase-matched nonlinear crystals for SHG. Compared with the method in Ref. 8, when the two pulses are overlapped entirely (no time delay) before entering into the crystals, the overlapping levels in the two crystals are the same and the final output is zero for this new measuring method, so there is no need to subtract an offset from the final output and no time delay adjustment is required between the two pulses. Furthermore, this new method provides a high sensitivity compared with the conventional RF method.
The principle of this new method is demonstrated in section 2 to show how the final output reflects the jitter. The mathematical model and computation originated from non-stationary nonlinear wave-coupled equations are given in section 3. The performance of jitter measurement including sensitivity and dynamic range are analyzed theoretically in section 4. Conclusions are given finally in section 5.
2. The measuring principle
Figure 1 shows the schematic diagram for the proposed method of jitter measurement between the transmitted pulse and the reference pulse using type II phase-matched SHG. Two fundamental pulses in one nonlinear crystal are required to be perpendicular polarized in the case of type II SHG . The polarization of the reference pulse is adjusted by a half wave plate (HWP) and split by a polarizing beam splitter (PBS) to obtain two perpendicular polarized parts with equal intensity. In this diagram, we assume the group velocity of the part transmitted by PBS is larger than the one reflected by PBS when they are propagated in the crystals. Therefore, the transmitted part can be called slow part and the reflected one can be called fast part. The transmitted pulse is adjusted and split in the same situation with the reference, except that the incident direction on the PBS is orthogonal to that of the reference. The fast part of the reference pulse and the slow part of the transmitted pulse are collinearly propagated in the first crystal (crystal-1) to generate a second harmonic signal (SH signal-1); meanwhile the slow part of the reference pulse and the fast part of the transmitted pulse is collinearly propagated in the second crystal (crystal 2) to generate another second harmonic signal (SH signal-2). Each second harmonic signal is detected by a photomultiplier detector (PMD). The final output is the difference between the two converted electric signals.
The intensity of the second harmonic field at the crystal end depends on the intensities of the two perpendicular polarized fundamental pulses and the time overlapping level between them. Because of equal splitting for both the reference and the transmitted pulse, the intensity difference between SH signal-1 and SH signal-2 depends only on pulse overlapping level in the two crystals.
The transmitted pulse and the reference will be overlapped entirely all the time before entering into the crystal if there is no time delay (no jitter) between them. In this case, when the two pulses are split by PBS and propagated in the two crystals, the fast part becomes ahead of the slow part gradually in the crystal-1, and the same situation happens in crystal-2. Since the transmitted pulse has the same wavelength with the reference, their slow parts have the same group velocity and so do the fast parts. In this case, the pulse overlapping levels keep the same in two crystals, so there is no difference between two second harmonic signals and the final output is zero.
The transmitted pulse will have a variable delay time relative to the reference before their incidence to the crystal if there is jitter between them. When the transmitted pulse is delayed relative to the reference, the average overlapping level in crystal-2 is larger than in crystal-1, so the SH signal-2 is larger than the SH signal-1 and the final output becomes positive. On the other hand, when the transmitted pulse is advanced relative to the reference, the final output becomes negative. Furthermore, the larger the delayed or advanced time is, the bigger the difference between the two second harmonic signals becomes. Therefore, the final output increases as the delay time increases. The sign of the final output indicates the jitter direction (the transmitted pulse is delayed or advanced relative to the reference), and the amplitude of the final output represents the jitter value. Accordingly the final output can reflect the timing jitter of the transmitted pulse relative to the reference.
In the application, the reference pulse is generated by a mode-locked laser which is stabilized to a frequency reference, and the transmitted pulse is the signal after propagating through the fiber link. In practice, the timing jitter cannot be measured directly between the transmitted pulse at the remote end and the reference at the local, so a partial reflector is used at the remote end to return part of the transmitted signal to the local through a round trip of the fiber link and the jitter measurement can be performed between this round-trip signal and the reference. In this case, the reference’s coherence time must be longer than the round-trip time of the pulse propagation through the fiber link in order to achieve an accurate jitter measurement.
3. Mathematical model and computation
The non-stationary nonlinear wave-coupled equations are used to study the proposed jitter measuring method. The variation of field amplitude with time must be considered in the nonlinear equations in the case of pulse SHG . First the intensity of second harmonic field of type II SHG for two fundamental pulses with no jitter is derived, and then the relation between the final output and the timing jitter for this measuring method is given.
3.1The intensity of the second harmonic field
The two fundamental pulses are required to be perpendicular polarized in a nonlinear crystal for type II phase-matched SHG . One polarized pulse propagates slowly and the other propagates fast in the crystal. The generated second harmonic field should have the same polarization with the fast fundamental field so as to satisfy phase matching [9,11]. The two fundamental slow polarized field, fast polarized field and the second harmonic field propagated along the z direction, denoted by E 1 s(z,t), E 1 f(z,t) and E2(z,t) can be expressed in the complex form
Considering the case of no absorption in the crystal, if the phase is totally matched (Δk = k 2-k 1 s-k 1 f), the nonlinear wave-coupled equations of the two fundamental fields and the second harmonic field can be expressed as [9,10]
In the situation of the fundamental fields being low depleted and approximately unchanged through the whole crystal, the amplitudes of the fundamental fields can be defined asEq. (2)c) becomesEq. (4) being integrated in the range of the whole crystal length, the amplitude of the output second harmonic field at the crystal end can be obtained as
By substituting the variation t′ for the original variation t at the crystal end z = L, we can express Eq. (5) asEq. (6) can be removed from the integral and replaced by its value at midpoint t 1 = (β + αs)L/2, then Eq. (6) can be solved asEquation (9) shows how the intensity of second harmonic field at the crystal end varies with time.
3.2 The relation between the final output and the timing jitter
When the transmitted pulse has timing jitter relative to the reference pulse, delay time τ between the two pulses varies with time continually. We assume τ > 0 when the transmitted pulse is delayed relative to the reference and τ < 0 when the transmitted pulse is advanced relative to the reference. With consideration of delay time τ between the two fundamental pulses, the intensities of the two second harmonic signals at the end of crystal-1 and crystal-2, denoted by I 2-1 and I 2-2, can be written as
We assume that the pulse widths of the fundamental pulses are identical. We investigate two cases of the pulse type: Gaussian pulses and hyperbolic secant pulses. Since both reference pulse and transmitted pulse are equally split by the PBS as we can see in Fig. 1, the pulse intensities of the Gaussian type can be expressed asEq. (12)a) and Eq. (12)b) into Eq. (11)a) and Eq. (11)b), by using the method of Fourier transform, we can obtain the expressions for SPD1 and SPD2 for Gaussian pulse type as followsEq. (12)c) and Eq. (12)d) into Eq. (11)a) and Eq. (11)b), by using the method of Fourier transform, we can obtain the expressions for SPD1 and SPD2 for hyperbolic secant pulse type as follows
It is worth noticing that ξ is the time difference of the fundamental slow and fast pulse propagating through the whole crystal length, which describes the walk-off between the two fundamental pulses owing to the different group velocities.
The differential value Sdiff is the final output. Equation (14) shows how the final output varies with the delay time of the transmitted pulse relative to the reference. This expression reflects the relation between the final output and the timing jitter of the two pulses.
4. Theoretical analysis of the performance for timing jitter measurement
Figure 2 shows the final output signal Sdiff varying with the delay time τ of the transmitted pulse relative to the reference in the case of ξ/TP = 2. The two cases of Gaussian pulse type and hyperbolic secant pulse type are given in Fig. 2(a) and Fig. 2(b) respectively. The approximate linear region blocked by the dashed line in the middle of the curve is applicable to the measurement. Since final output is proportional to the delay time, the final output can reflect the timing jitter in this linear region. The gradient of this linear region KS represents sensitivity of jitter measurement and the width Δτ represents dynamic range of jitter measurement. The sensitivity and the dynamic range indicate jitter measuring performance. As we can see, the properties of the two pulse types are similar, with only the numerical difference, therefore we can investigate the performance of the jitter measurement using one of the two pulse type, and the Gaussian pulse type is used in the following analysis.
In the linear region, the expression of the sensitivity KS can be obtained approximately as following through calculating the derivative of Eq. (14)a) at τ = 0
Figure 3 shows KS and Δτ vary as a function of ξ/TP. For a certain value of L, Ks increases with ξ/TP in the range of ξ/TP ≤ 2 and decreases in the range of 2 ≤ ξ/TP ≤ 4. Ks arrives at its maximum value at ξ/TP = 2, which can be obtained from Eq. (15). To enhance the sensitivity, the value of ξ/TP is expected to be close to 2. Furthermore, the sensitivity is proportional to the square of the crystal length L with the value of ξ/TP fixed. From this angle we hope the crystal is long and the value of αs-αf is small. For a given value of TP, the dynamic range Δτ has an increasing trend in the whole extent of ξ/TP ≤ 4. The maximum value of Δτ is 4TP at ξ/TP = 4.
The parameter ξ/TP is the normalized time difference between two fundamental pulses propagating through the crystal length. If the group velocities of the two pulses are the same and so there is no time difference (ξ/TP = 0), the final output is always zero and the delay time cannot be measured. In this case, the sensitivity and the dynamic range are both zero. With the increasing of this normalized time difference, the original time delay between the two pulses before entering into the crystal can increase without the loss of pulse overlapping level in the crystal; hence the larger value for ξ/TP leads to the bigger value for the dynamic range. However, when ξ/TP is fixed, the dynamic range is proportional to the pulse half-width TP. As regard to a mode-locked laser with pulse width at hundreds of femtoseconds level, the dynamic range of this jitter measurement method can achieve hundreds of femtoseconds to several picoseconds.
It is worth noticing that all the above curves and results are given in the case of T P≥(αs-αf)L/4, which means L≤4T P/(αs-αf). Further analysis for this measuring method shows that the sensitivity and the dynamic range for L>4T P/(αs-αf) are the same as L = 4T P/(αs-αf), therefore increasing the crystal length further will not improve the jitter measuring performance any more.
As an example, we choose BBO crystal for second harmonic generation at the fundamental wavelength of 800nm. In this case, υg 1 s and υg 1 f are 1.780 × 108m/s and 1.843 × 108m/s respectively , from which we can obtain that the value of αs-αf is about 192fs/mm. With the pulse half-width TP of a mode-locked laser about 100fs, the crystal length L needs to be about 2mm or longer in order to arrive at a maximal dynamic range of 384fs. The value of the sensitivity KS is not only related to the crystal and the two fundamental pulses, but also to the performance of photodetectors. In general, the sensitivity can achieve the level of ~10mV/fs.
In the process of the pulse transmitting through the fiber link, both phase noise and amplitude noise can be introduced. In this paper, only phase noise (timing jitter) is required to be measured. From the above analysis, we can see that in the dynamic range of this jitter measurement, the final output, which represents the intensity difference of the two second harmonic signals, is nearly proportional to the delay time between the two pulses. In the vicinity of zero timing offset, the amplitude noise makes the same contribution to the two second harmonic signals in Fig. 1, so the intensity difference of the two signals (the final output) is not affected by the amplitude noise and only affected by phase noise. Hence, the timing jitter measurement can be robust against the amplitude noise.
In the application of jitter measuring, the level of timing jitter can be characterized by the root-mean-square (rms) timing jitter, ΔT rms, which is specified over a frequency bandwidth range from f 1 (lower frequency of the bandwidth) to f 2 (upper frequency of the bandwidth). We define the power spectral density (PSD) of the timing jitter, ST(f), as the mean squared timing fluctuation at Fourier frequency f in a 1Hz bandwidth, then the rms timing jitter ΔT rms in the bandwidth (f 1, f 2) can be expressed as
The timing resolution of this rms timing jitter in the measuring process can reflect the measuring precision. The timing resolution of conventional RF method for timing jitter measurement is limited by the use of microwave components and there is drift timing error because of the thermal drifts of these microwave components. These limitations make it hard to measure the timing jitter better than the resolution of 100fs . The proposed method uses optical crystals instead of microwave components such as mixers to avoid these limitations so as to obtain more precise jitter measurement of two pulses in the same wavelength, and this method is not affected by the temperature drifts. The timing resolution of the proposed method can achieve a level of better than fs. In addition, the RF method can be affected by amplitude noise of the pulse introduced by fiber link transmission, while the proposed method can be robust against the amplitude noise. However, we should note that the dynamic range is limited by the pulse width in this method. Furthermore, comparing this proposed method with the previous optical method in Ref. 8, there is no need to subtract an offset from the final output and no time delay adjustment is required between the two pulses before entering into the crystals as we have mentioned in section 1.
In summary, a new measuring method for timing jitter of the transmitted pulse relative to the reference pulse has been proposed. The method is based on type II phase-matched SHG. The polarizations of the measuring pulses are exchanged in two nonlinear crystals and the final output can reflect the timing jitter between the two pulses. The non-stationary nonlinear wave-coupled equations have been used to obtain the theoretical result. The effect of the important parameter ξ/TP on the measuring sensitivity and dynamic range is analyzed. For a given crystal length L, the sensitivity increases in the range of ξ/TP ≤ 2 and decreases in the range of 2 ≤ ξ/TP ≤ 4 as ξ/TP increases. The sensitivity has a maximum value of 1.2π1/2 ηAL 2 I 01 I 02 at ξ/TP = 2. With ξ/TP fixed, the sensitivity is proportional to L 2. For a given mode-locked laser with certain pulse width, the dynamic range has an increasing trend at ξ/TP ≤ 4 and has a maximum value of 4TP at ξ/TP = 4. With ξ/TP fixed, the dynamic range is proportional to TP. The proposed method has a higher sensitivity compared with the conventional RF method while the dynamic range is limited to the pulse width. The overlapping levels in the two crystals are the same and the final output is zero when there is no time delay between the transmitted pulse and the reference. There is no need to subtract an offset from the final output and no time delay adjustment is required between the two pulses before entering into the crystals compared with the previous optical method. The proposed method can be applied in timing distribution, frequency transfer, and other areas requiring reliable and precise jitter measurement.
References and links
1. S. M. Foreman, K. W. Holman, D. D. Hudson, D. J. Jones, and J. Ye, “Remote transfer of ultrastable frequency references via fiber networks,” Rev. Sci. Instrum. 78(2), 021101 (2007). [CrossRef] [PubMed]
2. F. Narbonneau, M. Lours, S. Bize, A. Clairon, G. Santarelli, O. Lopez, C. Daussy, A. Amy-Klein, and C. Chardonnet, “High resolution frequency standard dissemination via optical fiber metropolitan network,” Rev. Sci. Instrum. 77(6), 064701 (2006). [CrossRef]
3. F. O. Ilday, A. Winter, J. W. Kim, J. Chen, P. Schmuser, H. Schlarb, and F. X. Kartner, “Ultra-low timing-jitter passively mode-locked fiber lasers for long-distance timing synchronization,” Proc. SPIE 6389, 63890L (2006). [CrossRef]
4. M. Calhoun, S. Huang, and R. L. Tjoelker, “Stable photonic links for frequency and time transfer in the deep-space network and antenna arrays,” Proc. IEEE 95(10), 1931–1946 (2007). [CrossRef]
5. K. W. Holman, D. J. Jones, D. D. Hudson, and J. Ye, “Precise frequency transfer through a fiber network by use of 1.5-microm mode-locked sources,” Opt. Lett. 29(13), 1554–1556 (2004). [CrossRef] [PubMed]
7. S. N. Bagayev, S. V. Chepurov, V. I. Denisov, A. K. Dmitriyev, A. S. Dychkov, V. M. Klementyev, D. B. Kolker, I. I. Korel, S. A. Kuznetsov, Y. A. Matyugin, M. V. Okhapkin, V. S. Pivtsov, M. N. Skvortsov, V. F. Zakharyash, T. A. Birks, W. J. Wadsworth, and P. S. J. Russell, “Femtosecond optical clock with the use of a frequency comb,” Proc. SPIE 4900, 125–131 (2002). [CrossRef]
8. J. Kim, J. Chen, Z. Zhang, F. N. C. Wong, F. X. Kärtner, F. Loehl, and H. Schlarb, “Long-term femtosecond timing link stabilization using a single-crystal balanced cross correlator,” Opt. Lett. 32(9), 1044–1046 (2007). [CrossRef] [PubMed]
9. S. Qian and R. Zhu, Nonlinear Optics (Fudan University Press, Shanghai, 2005), Chap. 3.
10. J.-C. Diels and W. Rudolph, Ultrashort Laser Pulse Phenomena: Fundamentals, Techniques, and Applications on a Femtosecond Time Scale, 2nd ed. (Academic Press, Burlington, 2006), Chap. 3.
11. M. V. Hobden, “Phase-matched second-harmonic generation in biaxial crystals,” J. Appl. Phys. 38(11), 4365–4372 (1967). [CrossRef]