We investigate general properties of the interferograms from a frequency comb laser in a non-linear dispersive medium. The focus is on interferograms at large delay distances and in particular on their broadening, the fringe formation and shape. It is observed that at large delay distances the interferograms spread linearly and that its shape is determined by the source spectral profile. It is also shown that each intensity point of the interferogram is formed by the contribution of one dominant stationary frequency. This stationary frequency is seen to vary as a function of the path length difference even within the interferogram. We also show that the contributing stationary frequency remains constant if the evolution of a particular fringe is followed in the successive interferograms found periodically at different path length differences. This can be used to measure very large distances in dispersive media.
© 2011 Optical Society of America
Laser sources with carrier-envelope-offset stabilized femtosecond pulses offer the unique advantage of a fixed phase relationship between the pulses emitted by the laser. The long term stability of these sources opens up a whole range of applications in several precision metrology application areas such as precision optical frequency metrology, spectroscopy, and absolute distance measurements [1–7], to name a few. Measurement of interferograms, or correlations, between pulses gives an insight of the coherence properties of the light source. A normalized correlation from a balanced Michelson interferometer, termed as the autocorrelation, exhibits a peak intensity of unity at exactly equal arms and decaying oscillations as a function of the position of the scanning arm. The autocorrelation is a symmetric function and its Fourier transform yields the power spectral density (PSD) of the laser source according to the Wiener-Khintchin theorem. If the interferometer is unbalanced and the measurements are carried out in vacuum, for a frequency comb source, the correlation repeats itself periodically with a maximum coherence when the path length difference between both arms is equal to an integer multiple of the laser cavity length. In this case the correlation patterns remain undistorted and the position of the maximum fringe visibility can be easily and accurately determined.
In the case of distance metrology [8–15] in air a dispersive medium is placed in the measurement arm. The observed interference patterns, which are cross-correlation patterns, will show distortion when compared to the autocorrelation. Also, the position of the coherence maximum or the brightest fringe will vary as a function of the delay distance. We have shown that , even if the medium is homogeneous, the position of the maximum coherence varies non-linearly for short delay distances and linearly for larger delay. In addition, our simulations also show that the correlation patterns are non-linearly broadened at short delay distances and that linear broadening replaces the non-linear broadening at larger distances.
In this work our main focus is on the formation of the fringes and the shape of the cross-correlation patterns for large delay distances. The terms large delay and short delay depend on the extent of distortion in the cross-correlation patterns which in turn is dependent on the properties of the light source and the dispersive medium. We describe the cross-correlation from a mode-locked laser using a discrete model and then extend this to a continuous model, allowing us to derive some dispersion properties of the correlation patterns. In particular, we have shown that the cross-correlation exhibits a non-linear broadening for short delay distances and linear broadening for larger delays. The theoretical results obtained from the developed continuous model are in full agreement with the standard discrete model of the cross-correlation as well as with the measured data. The depth of non-linearity for pulses in dispersive media has been shown to quantitatively depend on the temporal and spectral properties of the emitted light and the dispersion properties of the medium. In this work, using the continuous model of the cross-correlation and the method of stationary phase, we study the behavior of cross-correlation patterns and the frequency content of the constituting optical fringes at large delay distances.
This paper is organized as follows. In Section II, a general description of the cross-correlation using the continuous model is presented. In Section III, we use the method of stationary phase to study the cross-correlation functions and compare it to numerical simulations. Finally, in Section IV a new proposal for distance measurement using the method of stationary phase is described.
2. Cross-correlation in dispersive media
The starting point of our analysis is provided by a model describing the propagation of plane waves from a pulsed laser along the x-direction. A detailed description of our work can be found in . The frequency spectrum emitted by a mode-locked laser consists of a comb of regularly spaced frequencies ωm = mωr + ω0, m = 0, 1, 2,... where ω0 is the common offset frequency and ωr is the repetition frequency fr expressed in angular notation , . Here lpp is the cavity length, c is the velocity of light in the medium of the cavity and Tr is the time distance between the pulses. The offset frequency ω0 is caused by the difference between the group velocity and the phase velocity inside the laser cavity. Both ω0 and ωr are stabilized to an atomic clock in most laboratories.
The correlation function can be readily measured experimentally by placing a slow detector at the output of the Michelson interferometer, which can be balanced or unbalanced. A detailed analysis of this function in case of an unbalanced interferometer in the presence of a homogeneous dispersive medium leads to the following equation Equation (1), which is a discrete representation of the correlation function, shows good agreement with the experimental results but is unable to provide a physical explanation of some properties of the cross-correlation function such as the the shift of the position of the maximum coherence, broadening effects and shape of the cross-correlation function. A continuous model was therefore developed for a better understanding of the problem.
We can rewrite the cross-correlation as Eq. (2), for Γ(X) reduces to at most a single term when hX(t) has a support length ≤ Tr. The limitation posed by this for the case of propagation in air is discussed in a subsequent section. From a physical point of view, the integer ℓ denotes the multiple of the laser cavity length at a given delay distance X and ℓTr is the propagation time of a pulse in “vacuum”. In the case where X > 0 the integer ℓ must be negative. Defining the time variable (t) as 16], Eq. (4) we see that for any arbitrary delay distance X, independent of the laser cavity length 2πc/ωr, a cross-correlation pattern hX (t) can be obtained by varying the time delay (t) where −Tr/2 ≤ t ≤ Tr/2. In practice, this time (t) is obtained by setting up a scanning short arm of one laser cavity length. For numerical comparison, both hX (t) and Γ(X), are normalised to unity.
Using the Wiener-Khintchin theorem and assuming a quadratic dispersive medium k(ω + ω0) = αω2 + βω + γ for large delay distances, we can write16]. These results have allowed us to define a particular function ζ(X) by Equation (7) gives an indication of the effective distance of non-linear dispersion effects on cross-correlations obtained from a light source having a coherence time τc and a carrier frequency ωc propagating in a refractive medium with a group delay dispersion α at ωc. Short delay and large delay distances can be defined using 𝒟.
2.1. Typical characteristics of cross-correlations at large delay distances
In this section, we numerically analyze cross-correlation patterns for delay distances larger than 𝒟. For standard air we have previously shown  that 𝒟 = 30 m. Our simulations were performed for standard air using the measured laser source spectrum . The typical parameters for our Ti:Sapphire based frequency comb are as follows. The central frequency of the comb is ωc = 2.3254 × 1015 rad/s, corresponding to a wavelength of 810 nm in vacuum, the bandwidth is typically Δω ≈ 5 × 1014 rad/s, which correspond to a pulse width of Δx ≈ 12 μm and a pulse duration of 40 fs. The frequency offset is typically ω0 ≈ 113 × 107 rad/s and the repetition frequency is ωr = 6.28 × 109 rad/s, corresponding to a cavity length lpp = 30 cm and period Tr ≈ 1ns. The frequencies ω0 and ωr are synchronized to a cesium clock. The spectral content of the initial pulse is the main input for the numerical model. Using the corrected updated Edlén’s equation  for the refractive index of air, 8 × 104 spectral lines fitted to the profile of the spectrum of the laser, are propagated and then recombined to form correlation patterns.
In this numerical example the delay distance X has been taken to be long enough, 120 meters, so that the correlation patterns are observed to be linearly broadened. Three fringes are picked up from this correlation pattern and further analysed. The cross-correlation is shown in Fig. 1.a where the vertical lines indicate the positions of the fringes we investigate. The fringe pattern is obtained by varying the time delay t around a delay distance X = 120 m. At each small delay the frequency content of the interfering fields will be constantly in phase or out of phase depending on t. We have chosen to analyze the cosine of the phase factor from Eq. (4). For each small delay scan one fringe the cosine of the phase function has been plotted. A 3-D plot can be obtained for each analyzed fringe showing the value of the cosine of the phase as a function of the time delay (t) and the frequencies (or the wavelength) of the optical field. These plots are shown in Fig. 1.(b,c,d). The plots show that within one fringe the cosine of the phase function has either fast or slow oscillations as a function of the frequency. Using this oscillating function in Eq. (4), where S(ω + ω0) is a slowly varying term, the integral transform is small. The integral will approach zero as the number of oscillations increases, Riemann-Lebesgue lemma. Thus this relatively high oscillatory part will have a minor contribution to the formation of the interference fringe. Only when the cosine of the phase function is slowly varying the contribution to the integral will be important. The plots analyzing the fringes show that the slowly varying part is always at a specific frequency for a particular fringe. In Fig. 1.(e,f,g) we plot the cosine of the phase function for each fringe as a function of the frequency. Five such samples at different intensity points of the fringe are plotted on top of each other. It is clearly seen that, for a given fringe, the position of the stationary frequency remains constant on the frequency-axis and changes only in value when t changes. Between different fringes, the stationary frequency varies, depending on the fringe position in the correlation pattern.
We study the evolution of the cosine of the phase function from Eq. (4) as a function of the distance in Fig. 2. We pay particular attention to the behaviour of this function at the brightest fringe of the cross-correlations. For short distances the cosine of the phase function is slowly oscillating. In this case we have a large contribution from the whole frequency band of the source. As the distance increases, the cosine of the phase function becomes highly oscillating and the formation of the fringe depends on a small frequency bandwidth which is the slowly oscillating part. The width of the slowly oscillating part converges as the distance increases. The central frequency of the slowly oscillating part is constant and becomes a dominant frequency after reaching a delay distance larger than 𝒟 which is more than 30 m in standard air.
2.2. Limit of the validity of the continuous model in standard air
In this section, we make an estimation of the maximum pulse propagation distance for which the correlations can be analyzed using the continuous model. For this, it is essential that the series in Eq. (2) for the cross-correlation Γ(X) contains one significant term only. This is the case when the length of the interval outside which hX is negligible does not exceed Tr ≈ 10−9s. We rewrite Eq. (4) asEq. (11) into Eq. (8), we get Eq. (8), the integration range [0, ∞) for ω can be replaced by the essential integration range I which is, in accordance with the values in Section 2.1 for the bandwidth Δω and center frequency ωc, contained in the interval [2 × 1015,2.6 × 1015]. Hence, instead of gX above, we consider Eq. (12). Thus hX is obtained as the real part of the convolution of hX=0 and . Now, hX=0 is concentrated on a time interval of length, see Section 2.1, Δx/c ≈ 4 × 10−14 s ≪ 10−9 s = Tr. So, the length of the time interval where hX is concentrated is chiefly determined by the set of points s for which gX,I(s) is non-negligible.
A rough indication of the set of points s where gX,I(s) of Eq. 13 is non-negligible can be obtained by appealing to the stationary phase principle. Accordingly, the integral in Eq. (13) is non-negligible when the phase of the integrand has a stationary point ω inside the integration range I, and is expected to be small otherwise. The condition for ω to be a stationary point isEq. (14) can hold for an ω ∈ I, we need that |s| is at most of the order , where .
Now standard air is a mildly refractive medium, with refractive index n, see Eq. (9), given by Edlén’s equation . We estimate that the number Mξ,I has order of magnitude 10−7. With n̄ = 1.00027, c ≈ 3 × 108 m/s, we see that the length of the interval of points s where Eq. (14) holds is less than Tr = 10−9 s when , i.e., X ≤ 1.5 × 106 m.
This analysis gives an approximate indication of the size of the set where gX,I is nonnegligible, but for the relevant cases that X ≤ 1000 km it seems that we can safely assume that the series in Eq. (2), that uses samples of hX at distances Tr apart, has only one significant term. In any case, these distances are beyond the coherence length of present day laser systems but maybe are relevant in the future.
In the next section we will show that for large delay distances, particularly larger than 𝒟, an asymptotic method can be used to derive the equation of the cross-correlation function. Results from this asymptotic method and Eq. (5), which uses the continuous model, will be compared to each other. Using this, a simple method for absolute distance metrology is proposed.
3. Stationary phase approximation of the cross-correlation function
To study the correlation functions and their properties for large delay distances, we define large delay as the distance beyond which the correlation patterns are only linearly broadening. We have shown by numerical simulations and experimental measurements that the shape of the cross-correlation converges to a particular profile for large delay distances. In air, we were able to study this for correlations patterns up to 200 m, where the numerical results agreed reasonably well with the experiments .
To simulate the correlation formation in air under experimental conditions we used the standard discrete model. In spite of the reasonable agreement that has been obtained between numerical and experimental results, we were not able to explain the physical origin of the shape convergence and disappearance of the effect of chirp at large delay distances. To understand the formation of the cross-correlation patterns at large distances we extended the discrete model to a continuous model. This enabled us to write the cross-correlation as an integral, Eq. (4) instead of a series. Using this integral we now try to study the quasi-asymptotic formation of cross-correlation patterns.
The cross-correlation in Eq. (4) is expressed as an integral of a slowly varying function, S(ω + ω0), multiplied by the cosine of a phase factor. The stationary phase method is an appropriate mathematical tool to study the asymptotic behaviour of such integrals . We rewrite Eq. (4) asEq. (15). This general class of oscillatory integrals, with polynomially growing phase functions, also has a solution where higher-order stationary phase approximation can be calculated. The asymptotic solution of Eq. (15) using only the first and the second order stationary phase approximation is given by 18]
Then, the first order asymptotic approximation of the cross-correlation function is given as19]. In this case, terms arising due to higher order stationary phase approximations yield minor corrections to the leading order approximations of the cross-correlation function. In the subsequent analysis we use the first and the second order terms and neglect the rest. We see from Eq. (19) that the enveloppe of the first order approximation is determined by the spectral function S(ω + ω0). This is in agreement with the previous experimental and numerical simulation , where we observe shape convergence of cross-correlation functions at large delay distances. The equation also shows that the term inside the cosine depends linearly on ωdom. Thus, with increasing delay distance the cross-correlation functions will spread linearly.
We compare various simulated correlations using Eq. (1) and their corresponding asymptotic functions using the stationary phase approximations from Eq. (16) in Fig. 3. We refer to the cross-correlation from Eq. (1) and the one from Eq. (16) as the exact and the asymptotic cross-correlations respectively. We start our comparison using the first order stationary phase method. Numerical results for 100 m delay distance in air are shown in Fig. 3.a. At the center of the interferogram a good match between the exact and the asymptotic correlation can be observed. An analysis of the wings of the interferogram show a fringe mismatch. A better match is obtained for a delay distance of 200 m as is shown in Fig. 3.c. This problem of shape matching between the exact and the asymptotic correlations can be alleviated using higher-order stationary phase approximation. Here we only show the improvements using the second-order stationary phase approximation. It is clearly seen that the second order stationary phase approximation yields a better agreement with the exact correlation. More importantly, for distance measurement purposes, the position of the central fringe of the asymptotic correlations has not shown any shift in comparison to the exact simulations.
We found that for larger distances, fringe patterns can be described by particular frequencies. In Fig. 4 we show the frequency content of a correlation pattern simulated for a delay of 200 m in air. The top figure shows the spectral distribution of the stationary frequencies after 200 m propagation in air. The bottom figure shows the cross-correlation function simulated using the first-order stationary phase method for the same propagation distance. The color indicates the frequency content ranging from 780 to 850 nm. It is clearly seen that the stationary frequencies change within a given cross correlation depending on the fringe order. In addition, the shape of the envelope of the cross-correlation follows the shape of the spectral distribution as was demonstrated in this section. Using this approximation greatly reduces the computational time when simulating the cross-correlation functions as can be seen when Eq. (19) and Eq. (1) are compared. As an example, the spectrum of the Ti:Saphirre frequency comb contains 105 frequencies. By using a scanning short arm with 104 steps, one has to compute a grid of 105 × 104 elements in order to obtain the final results. This can be simply reduced to a vector of 104 components by using the cross-correlation equations from the stationary phase method since we can associate one dominant frequency to each scanning step.
4. Stationary Phase Absolute Distance Metrology
In the previous section we observed that the analysis of cross-correlation patterns becomes simpler for large delay distances, because the stationary phase approximation can be applied. Each fringe in the cross-correlation is formed by its unique dominant frequency and the cross-correlation broadens only linearly. In this section we will show that by accurate knowledge of the spacing between two fringes of different frequencies we can determine the absolute distance of the measurement arm in dispersive media at large delays.
We would like to compare Eq. (5) and Eq. (19). Both equations are only valid for large delay distances and are calculated for a quadratic dispersive medium. In the result obtained from the first order stationary phase approximation as shown in Eq. (19), the cross-correlation has the envelope of the power spectral density of the laser source, S(ωdom + ω0), at the dominant frequencies ωdom. The envelope of the cross-correlation from Eq. (5) is given by the following integral
The exponential can be written as exp[iωp(X,t)t0], where ωp(X,t) = (t +βX − n̄X/c) / (2αX) is a particular angular frequency depending on the delay distance and the scanning time t. Here, α and β are the related dispersion properties of the medium as mentioned in Section 2. Therefore we can replace the above integral with S[ωp(X,t)] where we have used the Wiener-Khintchin theorem. Here S[ωp(X,t)] is the power spectral density at the particular angular frequencies ωp(X,t). Thus we can see that the envelope of the cross-correlation patterns formed after long propagation in quadratic dispersive media is determined by the shape of the PSD of the laser source. We observe a linear dependence between the stationary frequencies, ωp(X,t), and the scanning time t. In the regime where the extent of the cross-correlation function is smaller than the laser cavity length, the value of t is given by Eq. (3).
In the case where the output of the unbalanced Michelson interferometer is analyzed with a spectrometer, we would measure modulated spectra instead of cross-correlation patterns. Assume that we perform two independent measurements at the same distance X but at two different scanning positions of the short scanning arm, t1 and t2, respectively. From the expression for ωp(X,t) we obtain
The continuous model of cross-correlations allows us to understand many aspects of the formation of correlation patterns after propagation in dispersive media. The cross-correlation patterns are subject to non-linear broadening at short path length differences and linear broadening at larger path length differences as has been observed in measurements and by simulations. For short path length differences a large part of the spectral content of the pulse contributes in the formation of the correlation patterns and in particular the brightest fringe. The phase accumulation of the different frequencies in the spectral content contributing to the correlation pattern is non-linear because the change in refractive index varies non-linearly for different frequencies as seen from the Edlén equation. Therefore, for the initial part of the propagation the position of the maximum coherence varies non-linearly and the broadening of the correlation is also non-linear. For large path length differences very few frequencies of the spectra contribute to the fringes of the correlation pattern as can be seen from the analysis presented in this work and thus the non-linear phase accumulation plays a smaller role. Therefore both the variation in the position of the maximum coherence and the broadening are seen to be linear. For path length differences beyond the non-linear dispersion depth, the method of stationary phase was used to study the correlation patterns. It shows that the shape of the cross-correlation patterns at large delay distances is determined by the source spectral profile, though the cross-correlation will continue spreading linearly. It was also observed that each intensity point of the correlation pattern is formed by the contribution of one dominant stationary frequency. This stationary frequency is seen to vary as a function of the path length difference within the correlation pattern. We also show that the contributing stationary frequency remains constant if the evolution of a particular fringe is followed in the successive cross-correlation patterns found periodically at different delay distances for the long arm. Using this property a method of measuring very large distances has been proposed.
The authors would like to thank Omar El Gawhary for extensive discussions. The authors gratefully acknowledge J. J. M. Braat for his thoughtful comments on the manuscript. We also thank the peer reviewers for their helpful comments. This work is part of the Industrial Partnership Programme (IPP) Metrology with Frequency Comb Laser (MFCL) of the Stichting voor Fundamenteel Onderzoek der Materie (FOM), which is supported financially by Nederlandse Organisatie voor Wetenschappelijk Onderzoek (NWO). The IPP MFCL is co-financed by VSL, TNO and ASML. The research within this work is also a part of EURAMET joint research project and has received funding from the European Community’s Seventh Framework Programme, ERA-NET Plus, under Grant Agreement No. 217257.
References and links
1. D. J. Jones, S. A. Diddams, J. K. Ranka, A. Stentz, R. S. Windeler, J. L. Hall, and S. T. Cundiff,“Carrier-Envelope Phase Control of Femtosecond Mode-Locked Lasers and Direct Optical Frequency Synthesis,” Science 288, 635–639 (2000). [CrossRef] [PubMed]
2. R. Holzwarth, Th. Udem, T.W. Hänsch, J. C. Knight, W. J. Wadsworth, and P. St. J. Russell, “Optical Frequency Synthesizer for Precision Spectroscopy,” Phys. Rev. Lett.85, 2264–2267 (2000). [CrossRef] [PubMed]
3. S. T. Cundiff and J. Ye, “Colloquium: Femtosecond optical frequency combs,” Rev. Mod. Phys. 75, 325–342 (2003). [CrossRef]
6. L. Hollberg, C. W. Oates, E. A. Curtis, E. N. Ivanov, S. A. Diddams, T. Udem, H. G. Robinson, J. C. Bergquist, R. J. Rafac, W. M. Itano, R. E. Drullinger, and D. J. Wineland, “Optical frequency standards and measurements,” IEEE J. Quantum Electron. 37, 1502–1513 (2001). [CrossRef]
7. J. Ye, H. Schnatz, and L. W. Hollberg, “Optical frequency combs: from frequency metrology to optical phase control,” IEEE J. Sel. Top. Quantum Electron. 9, 1041–1058 (2003). [CrossRef]
9. M. Cui, R. N. Schouten, N. Bhattacharya, and S. A. van den Berg, “Experimental demonstration of distance measurement with a femtosecond frequency comb laser,” J. Eur. Opt. Soc. Rapid Publ.308003 (2008) [CrossRef]
10. Y. Salvade, N. Schuhler, S. Leveque, and S. Le Floch, “High-accuracy absolute distance measurement using frequency comb referenced multiwavelength source,” Appl. Opt. 47, 2715–2720 (2008) and references therein. [CrossRef] [PubMed]
13. M. Cui, M. G. Zeitouny, N. Bhattacharya, S. A. van den Berg, H. P. Urbach, and J. J. M. Braat, “High-accuracy long-distance measurements in air with a frequency comb laser,” Opt. Lett. 34, 1982–1984 (2009). [CrossRef] [PubMed]
14. I. Coddington, W. C. Swann, L. Nenadovic, and N. R. Newbury, “Rapid and precise absolute distance measurements at long range,” Nature Photon. 3, 351–356 (2009). [CrossRef]
15. J. Lee, Y.-J. Kim, K. Lee, S. Lee, and S. Kim, “Time-of-flight measurement with femtosecond light pulses,” Nature Photon. 4, 716–720 (2010). [CrossRef]
16. M. G. Zeitouny, M. Cui, N. Bhattacharya, S.A. van den Berg, A. J. E. M. Janssen, and H. P. Urbach, “From a discrete to a continuous model for interpulse interference with a frequency-comb laser,” Phys. Rev. A 82. 023808 (2010). [CrossRef]
17. K. P. Birch and M. J. Downs, “Correction to the Updated Edlén Equation for the Refractive Index of Air,” Metrologia 31, 315–316 (1994). [CrossRef]
18. V. A. Borovikov, Uniform Stationary Phase Method, IEE Electromagnetic Wave Series (1994).
19. K. E. Oughstun and N. A. Cartwright, “Physical significance of the group velocity in dispersive, ultrashort gaussian pulse dynamics,” J. Mod. Opt. 52, 1089–1104 (2005 )and references therein. [CrossRef]