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

Long distance measurement by dynamic optical frequency comb

Open Access Open Access

Abstract

In this paper, we propose a method aiming to measure the absolute distance via the slope of the inter-mode beat phase by sweeping the repetition frequency of the frequency comb. The presented approach breaks the inertial thinking of the extremely stable comb spacing, and the bulky phase-locking circuit of the repetition frequency is not required. In particular, the non-ambiguity range can be expanded to be infinite. To verify the performance of presented method, a series of distance experiments have been devised in different scenarios. Compared with the reference values, the experimental results show the differences within 25 µm at 65 m range in the laboratory, and within 100 µm at 219 m range out of the lab.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

During the past decades, optical frequency combs (OFCs), honored by the Nobel Prize, have attracted the attention of the whole world [13], and have seen a variety of applications as a versatile tool, e.g., frequency metrology [4,5], precision spectroscopy [6,7], time/frequency transfer [8,9], astrocombs [10,11], and absolute distance measurement [1214], etc. In the optical domain, OFC is composed of a series of equally-spaced frequency modes covering across a board spectral band, and these coherent lines can be expressed as n×frep+fceo clearly, where n is an integer, frep is the comb spacing, and fceo is the carrier-envelope-offset frequency. As long as frep and fceo are tightly locked to a RF frequency reference by a phase-locking loop, the bridge is constructed between the optical frequency and the RF frequency, and all the spectral lines will share the same frequency stability as the external reference. This inherent advantage has given rise to the substantial progress in the field of absolute distance measurement, which is of fundamental importance in both science and technology.

In the time domain, OFC is in the form of a pulse train with stable time interval, and the pulse period can be expressed simply as 1/frep. By means of the time-domain characteristics of frequency comb, researchers have made great progress during the past years. The way how to extract the accurate distance cues is the key point, and the measuring schemes including the pulse-to-pulse alignment [15], peak extraction [16], intensity detection [17] and the slope of the unwrapped phase [18] are well analyzed and conducted. However, the scanning of mechanical stage is indispensable to shift one pulse with regard to the other, to generate the interferograms. This long moveable component directly results in that the optical alignment in the free space is not easy, and the measuring speed is low. To overcome the above defects, researchers proposed an ingenious method called dual-comb ranging, which relies on two mode-locked lasers with slightly different repetition frequencies, and the pulses can automatically scan each other at a rapid pace to generate the interferograms. The limitation is that two well-stabilized mode-locked lasers are needed, making the system expensive and bulky [1921]. In the frequency domain, methods of dispersive interferometry [2224], chirped pulse interferometry, multi-heterodyne interferometry [25,26], and multi-wavelength interferometry [27,28] are proposed from diverse perspectives. Dispersive interferometry and chirped pulse interferometry can measure the distances via the appropriate processing of the spectral interferograms. Nevertheless, the measurement performances are limited by the resolution of the optical spectrum analyzer. The phases of the inter-mode beat can be used to determine the distances, and in general the high harmonic is needed to finely measure the distances with a smaller synthetic wavelength [29,30].

The frequency comb offers an intriguing solution to achieve the higher accuracy and the greater non-ambiguity range simultaneously, and of which the attraction of greater measuring range has prompted various efforts to improve the measuring scheme. In particular, the non-ambiguity range can be significantly extended to kilometer level by slightly changing the repetition frequency of the signal laser, for all the methods mentioned above (with one single comb or dual combs) in the time and frequency domains. However, in this case, the additional measurement step of changing the repetition frequency is required to unequally determine the distances. In this work, we propose a method enabling long distance measurement by using one single laser frequency comb. The distance can be measured by the phase slope of the inter-mode beats, while the repetition frequency is swept linearly. Compared with the general pulsed laser, the dynamic optical frequency comb can be achieved by linearly changing the cavity length, and the repetition rate can be controlled precisely through a high-precision linear stage, which greatly maintains the high coherence of frequency comb. Furthermore, the ultra-short pulse duration of optical frequency comb corresponds to a broadband spectrum, which can achieve the extension of the stable higher harmonics, and greatly improve the measurement accuracy. Surprisingly, the method can achieve arbitrary distance measurement in contrast to conventional time-of-flight methods, and the non-ambiguity range can be extended to be infinite in principle, which is very valuable to the large-scale measurement. The complex locking circuit for the stabilization of the repetition frequency is not needed any more, which can well satisfy the field-base experiments. The experimental results show that our method can measure the distances with high accuracy and precision.

2. Measuring principle based on dynamic comb spacing

Theoretically, the measuring system is devised as shown in Fig. 1. The output of the laser source is divided into three channels. One is the monitoring channel (green), which is used to monitor the real-time repetition frequency. One is the reference channel (blue), and the third one is the measurement channel (black). The repetition frequency of the laser source is linearly swept, and accordingly the phase difference between the measurement and reference channels would change linearly. The different distance L will correspond to a different phase slope, which means that the distance can be measured.

 figure: Fig. 1.

Fig. 1. Structure of measuring system

Download Full Size | PDF

2.1 Dynamic comb spacing based on fundamental inter-mode beat

The phases of the fundamental inter-mode beat for the reference and measurement channels (i.e., ϕref and ϕmeas) can be expressed respectively as:

$${\phi _{ref}} = \frac{{2\pi r}}{\lambda }$$
$${\phi _{meas}} = \frac{{2\pi (r + 2L)}}{\lambda }$$
where λ is the pulse-to-pulse interval, λ=c/(Ng·frep), c is the light speed in vacuum, Ng is the group refractive index, r is the length of the reference channel, and L is the measured distance. Equation (2) can be rewritten as:
$${\phi _{meas}} = \frac{{2\pi r}}{\lambda } + \frac{{4\pi L{N_g}}}{c} \cdot {f_{rep}}$$
Hence, the phase difference ϕ can be calculated via Eq. (1) and Eq. (3) as:
$$\phi = \frac{{4\pi L{N_g}}}{c} \cdot {f_{rep}}$$
With the sweeping of frep, the difference of phase ϕ will generate the linear modulation. Therefore, the measured distance L can be calculated from the slope of phase change as:
$$L = \frac{{\partial \phi }}{{\partial {f_{rep}}}} \cdot \frac{c}{{4\pi {N_g}}}$$
The phase difference can be greater than 2π as the measured distance increased, and it will cause the nominal ambiguous range for distance measurement. Hence, the non-ambiguity range LNAR needs to satisfy the following condition as:
$$\frac{{4\pi L{N_g}}}{c} \cdot {f_{rep1}} - \frac{{4\pi L{N_g}}}{c} \cdot {f_{rep2}} < 2\pi $$
where frep1 and frep2 are adjacent values in the process of repetition frequency variation, and LNAR can be collated as:
$${L_{NAR}} = \frac{c}{{2{N_g}\Delta {f_{rep}}}}$$
where Δfrep is the difference of contiguous change for the comb spacing, it also can be considered as a resolution of variation. We can observe that Δfrep is the decisive factor of non-ambiguity range from the above formula. Fortunately, the Δfrep could be tuned by external trigger of the signal generator. In theory, the Δfrep can be adjusted to infinitesimal if the trigger frequency ft is high enough.

2.2 Dynamic comb spacing based on higher harmonic

From the Eq. (5), with the sweeping of repetition frequency, the small phase variation will be the challenge we may encounter. In response to this problem, we utilize the higher harmonic as an “amplification factor” acting on the phase slope. Then, the phase difference ϕm between the reference and measurement channels can be redefined by:

$${\phi _m} = \frac{{4\pi L{N_g}m}}{c} \cdot {f_{rep}}$$
where m is an integer (i.e., the order of the higher harmonic). The slope of phase change can be derived as:
$$\frac{{\partial {\phi _m}}}{{\partial {f_{rep}}}} = \frac{{4\pi L{N_g}m}}{c}$$
Consequently, the measured distance L can be calculated as:
$$L = \frac{{\partial {\phi _m}}}{{\partial {f_{rep}}}} \cdot \frac{1}{m} \cdot \frac{c}{{4\pi {N_g}}}$$
Compared with Eq. (5), Eq. (10) involves the amplification factor m. Due to the higher harmonic (i.e., the smaller synthetic wavelength), the phase change can be improved to be more sensitive, which means that the distances can be finely measured.

3. Demonstration of the Yb-doped mode-locked laser

The mode-locked laser is generated based on the mechanism of nonlinear polarization rotation (NPR), using Yb-doped gain fiber. The high-gain efficient of the Yb-doped fiber makes the system easy to achieve the mode-locked state, and can maintain the high output power. The laser configuration is shown in Fig. 2. The fiber-ring cavity consists of the Yb-doped fiber (25 cm), single mode fiber, wavelength division multiplexing WDM, two collimators, optical isolator (Thorlabs IO-5-1030-VLP), and the NPR module. One collimator is installed at a mechanical stage (50 mm moving range), which is used to tune the repetition frequency slowly and largely. Part of the single mode fiber is stuck on a piezoelectric transducer (Thorlabs PK2FSP2), so that the repetition frequency can be fast and finely adjusted. In our experiments, the repetition frequency is about 100 MHz, with tuning range of about 1 MHz. The output power of the laser source is about 100 mW.

 figure: Fig. 2.

Fig. 2. Schematic of homemade frequency comb. LD, laser diode; WDM, wavelength division multiplexer; PZT, piezoelectric transducer; QWP, quarter-wave plate; HWP, half-wave plate; PBS, polarization beam splitter; ISO, isolator.

Download Full Size | PDF

Figure 3(a) shows the generated pulse train, with the period of 9.65 ns. The spectrum of the laser source is indicated in Fig. 3(b). We find that, the center wavelength is about 1030 nm, with about 15 nm spectral width (FWHM). The repetition frequency can be locked to the Rb clock by using the phase-locking loop [14]. Figure 4(a) shows the scatters of the stabilized repetition frequency for 5 hours, and the standard deviation is 0.2 mHz. The Allan deviation, shown in Fig. 4(b), can achieve 1.8 × 10−12 at 1 s, and 5.9 × 10−14 at 1000-s averaging time.

 figure: Fig. 3.

Fig. 3. (a) Pulse train with 103.6 MHz repetition frequency. (b) Spectrum of the laser source.

Download Full Size | PDF

 figure: Fig. 4.

Fig. 4. (a) Repetition rate variation. (b) Allan deviation of the relative repetition rate. The counter (53230A, Keysight) is synchronized to the Rb clock, whose gate time is 0.1 s.

Download Full Size | PDF

As shown in Fig. 5, the input voltage to the PZT (piezoelectric transducer) has a minor periodic fluctuation, which leads to a nonlinear phase error. The actual level of fluctuation of input voltage is 8.9 mV, and the standard deviation is 1.8 mV. The displacement of PZT (Thorlabs, PK4GA7P2) is lower than 0.67 µm/V at 25 °C. Hence, the impact on the cavity length is less than 5.9 nm, and the destabilization of repetition rate Dfrep is 0.15 Hz with regard to the input voltage of PZT. When the measured distance L is 65 m, the nonlinear phase error ϕerror can be calculated as:

$${\phi _{error}} = \frac{{4\pi L{N_g}}}{c} \cdot {D_{{f_{rep}}}}$$
where c is the light speed in vacuum, Ng is the air group refractive index. The level of periodic nonlinear phase error is below 4.0 × 10−7 rad, the standard deviation is 8.1 × 10−8 rad and the resulting error for the distance measurement is 9.4 × 10−8 m. We could find that the repetition rate fluctuation in Fig. 4 (1 mHz) is much less than the predicted value and consistent with the frequency stability of Rb clock. Due to the excessively marginal impact, the nonlinear phase error with respect to the input voltage can be ignored.

 figure: Fig. 5.

Fig. 5. Nonlinear phase error with the input voltage to the PZT. The left Y axis is the input voltage to the PZT, and the right Y axis is the nonlinear phase error.

Download Full Size | PDF

4. Long distance measurement using phase slope of the inter-mode beat

The experimental setup is shown in Fig. 6. The output of the laser source is split into three channels, as described in Sec. 2. To meet the requirements of ultralong distance experiments, the laser is amplified by Yb-doped fiber amplifier YDFA. Corresponding to the principle of the method, three photodetectors (MenloSystems FPD310-F) are used to measure the repetition frequency, the phase of the reference arm, and the phase of the measurement arm, respectively. As depicted in Fig. 6, two frequency counters (Keysight 53230A), which are synchronized by a signal generator (Tektronix AFG31000), are used as a tool to record repetition rates and phase signals. In the process of data acquisition, the measuring system shares the common time base of the rubidium clock (Microsemi 8040) as an external reference.

 figure: Fig. 6.

Fig. 6. Schematic overview of the measuring setup

Download Full Size | PDF

As shown in Fig. 6, the measurement and reference signals pass through a low-pass filter with a cut-off frequency of 100 MHz. After being amplified by the amplifier, they are connected to the frequency counter, and the fundamental frequency is electronically obtained. For the way of determining high inter-mode beats, the measurement and reference signals pass through a band-pass filter with a desired frequency band. After being amplified by the amplifier, they are connected to a mixer meeting the bandwidth requirement, and are mixed with the reference signal to obtain a low-frequency signal that retains the original phase difference information. Considering the change of the mixing signal caused by dynamic higher harmonics, a low-pass filter with a cut-off frequency of 5 MHz is subsequently connected, and finally they are input into the frequency counter after the secondary amplification of the amplifier.

To acquire the sweeping repetition rate accurately, both the frequency counter and the sampling signals output by the signal generator use the Rb clock as an external reference. At the same time, the frequency counter, which detects the phase change, is also locked to the same external reference.

4.1 Absolute distance measurements on the long optical rail

Based on the principle of approach, we performed a series of experiments on the 80 m long optical rail, at National Institute of Metrology, which has the stable environment conditions (T: 24.08 °C, P: 991.93 hPa, RH: 43.5%), and the group refractive index is 1.0002624 based on the Ciddor formula [31]. The experimental photograph is shown in the Fig. 7.

 figure: Fig. 7.

Fig. 7. Experimental photograph on the long rail.

Download Full Size | PDF

The incremental distance meter (Agilent 5519B) is used as the reference, and a series of experiments have been devised for different measured distances with 5 m step. In preparation for the measurement, the measurement arm is ensured to be parallel to the optical rail. By means of the long optical rail, the distance measurements have been performed, covering 65 m. With the dynamic comb spacing of fundamental frequency, the original data of phase changes are illustrated in Fig. 8.

 figure: Fig. 8.

Fig. 8. Experimental results for the fundamental frequency. (a) Sweeping of repetition rate for the experiment. (b) Measured data of unwrapped phase corresponding to 5 m, 35 m, 65 m.

Download Full Size | PDF

Figure 8 reveals that, the experimental results show a good agreement with theoretical derivation. We could observe the phase slope increases continuously with the increase of measured distance. It implies the measured distance has certain correspondence with the phase slope as the theoretical derivation. As shown in Fig. 8, the phase change is more distinct with the increase of measured distance. Nevertheless, the phase jitter is the problem we must solve for the calculation of phase slope. Due to the environmental disturbances, photoelectric conversion and the interference from electronics, etc., the acquired phase signal accompanying phase jitter is unavoidable. By means of the verification of original phase data from repeated experiments, the actual level of phase jitter is 8.6 mrad, and the standard deviation attains 6.7 mrad.

Based on the above factors, the acquired phase signals will generate phase jitter caused by random errors and non-linear periodic errors. The phase jitter has the characteristics of high frequency, hence, the frequency of phase jitter will be much higher than the frequency of phase change generated by the dynamic comb spacing. Compared with traditional Fourier transform, The Empirical Mode Decomposition (EMD) can analyze the local time-frequency characteristics of the signal [32], instead of performing a global transformation on the signal to obtain the overall frequency spectrum. Based on the time scale characteristics of the data itself, the method can decompose the signal into several Intrinsic Mode Functions (IMFs), and reconstruct the phase signals from component IMFs (C1, C2, C3…), as shown in Fig. 9. Because of time-frequency localization characteristics of EMD, the method is suitable for non-stable signals. For the situation of phase signal instability, phase jitter can be solved by taking advantage of EMD. Due to the local characteristics of the phase change, the IMF (C12) that matches the phase change is approximately linear. In the subsequent processing, the linear change of phase is obtained by linear fitting. We fit phase slope at different regions of the phase variation and average them to measure the distance, thereby eliminating the effect of phase jitter.

 figure: Fig. 9.

Fig. 9. Data processing by using EMD. It is clear that the phase jitter is greatly suppressed after EMD.

Download Full Size | PDF

Therefore, based on the phase slope of the inter-mode beat, the distances can be measured. In our experiments, each distance is fast measured for seven times. The gate time of the frequency counter is 1 ms, i.e., 1 kHz sampling rate. Compared with the reference distances, the residuals are shown in Fig. 10. Among that the cross marks represent the scatters of the measurements, and the solid mid points indicate the average of seven measurements, and the error bars show the standard deviation. We find, the residuals can be below 3 mm at 65 m range.

 figure: Fig. 10.

Fig. 10. Experimental results obtained by the fundamental frequency signal.

Download Full Size | PDF

Next, we use the 9th harmonic of repetition frequency (about 932.4 MHz) to finely measure the distances. Correspondingly, the gate time of the frequency counter is updated to 0.11 ms, i.e., 9 kHz sampling rate. The phase changes that are caused by the dynamic comb spacing have been depicted in Fig. 11.

 figure: Fig. 11.

Fig. 11. Experimental results of the 9th harmonic. (a) 9th harmonic of repetition rate sweeping. (b) Experimental data of unwrapped phase change for 5 m, 35 m, 65 m based on the 9th harmonic.

Download Full Size | PDF

In contrast to Fig. 8, Fig. 11 clearly indicates the fine step using the 9th harmonic of the repetition frequency shows more sensitive phase changes, which is actually magnified by 9 times. Each distance is measured for 21 times. The experimental results are depicted in Fig. 12. In the figure, the cross marks represent the scatters of the measurements, the solid mid points indicate the average of 21 measurements, and the error bars show the standard deviation. The residuals can be well below 25 µm at 65 m range, which is 3.8×10−7 in relative.

 figure: Fig. 12.

Fig. 12. Experimental results based on the 9th harmonic signal at 65 m range.

Download Full Size | PDF

4.2 Long distance measurement on the base line out of the lab

To further validate the performance of the presented approach, and verify the practical application of frequency comb outside the lab, open-air experiments are performed on the base line. The experimental scheme is shown in Fig. 13.

 figure: Fig. 13.

Fig. 13. (a) Experimental photograph on the base line at 1:00 am. (b) Top view of the measurement path.

Download Full Size | PDF

The experiments are performed at early morning with relatively stable environment (T: 20.06 °C, P: 988.28 hPa, RH: 47.4%) corresponding to the group refractive index of 1.0002651. We measure the distances at the initial position, 9, 27, 75, and 219 m, respectively. The phase changes at 219 m position for the fundamental and the 9th harmonic are indicated in Fig. 14, and the distance can be measured via the phase slope.

 figure: Fig. 14.

Fig. 14. Experimental results corresponding to 219 m. (a) Phase changes based on the fundamental frequency. (b) Phase changes based on the 9th harmonic.

Download Full Size | PDF

Figure 15 shows the experimental results of distance measurement. The reference values are provided by National Institute of Metrology. In the step of coarse measurement using fundamental frequency, the residuals can be within 4 mm. In the case of fine measurement using 9th harmonic, the residuals can be greatly improved to be below 100 µm.

 figure: Fig. 15.

Fig. 15. Distance measurement results over a travel of 219 m. (a) Experimental results based on the fundamental frequency. (b) Experimental results based on the 9th harmonic.

Download Full Size | PDF

We performed long-term measurements at the distance of 219 m, based on the fundamental frequency and the higher harmonic respectively. The scatters of residuals and Allan deviation are depicted in Fig. 16, and of which the data rates are 1 Hz and 5 Hz for the fundamental and the higher harmonic respectively. The standard deviations, for the fundamental and the higher harmonic, are 2.2 mm and 51.6 µm, respectively. The Allan deviation for the fundamental beat can reach 2.2 mm at 1 s, and 103.7 µm at 100 s averaging time. In the case of higher harmonic, the Allan deviation can achieve 20.4 µm at 1 s, and 5.6 µm at 100 s averaging time.

 figure: Fig. 16.

Fig. 16. Long-term experiments of the distance measurement at 219 m. (a) Measured residuals for fundamental frequency. (b) Measured residuals based on the higher harmonic. (c) Allan deviation of the distance measurement versus averaging time.

Download Full Size | PDF

4.3 Uncertainty evaluation

The uncertainty of measured distance is a tool that we appraise our method. Based on Eq. (10), the main influence factors of measurement uncertainty are related to the slope of phase S=∂ϕ/∂frep and the group refractive index Ng, hence it can be calculated as:

$${u_L} = \sqrt {{{(\frac{{\partial L}}{{\partial S}} \cdot {u_S})}^2} + {{(\frac{{\partial L}}{{\partial {N_g}}} \cdot {u_{{N_g}}})}^2}} = \sqrt {{{(\frac{1}{m} \cdot \frac{c}{{4\pi {N_g}}} \cdot {u_S})}^2} + {{(\frac{L}{{{N_g}}} \cdot {u_{{N_g}}})}^2}} $$
For the fundamental frequency, the uncertainty of the term related to the phase slope could be considered to be 1.3×10−10 rad/Hz (standard deviation). With the increase of measured distance, the standard deviation of phase slope is continuous decrease due to the slope gets larger. For the term related to the air refractive index, we used the Ciddor equation to correct the air refractive index. Thereinto, the inherent uncertainty of empirical equation is 2.0 × 10−8. In our measuring system, the inherent uncertainty of the sensors, the environmental stability, and the environmental inhomogeneity within the coverage of 65 m range should be considered, the uncertainties of temperature, air pressure and humidity are 35.2 mK, 10.3 Pa, and 2.1%, respectively. The uncertainties are corresponding to 3.2 × 10−8, 2.7 × 10−8, and 2.3 × 10−8 for the air refractive index, respectively. The combined uncertainty related to refractive index can attain 5.2 × 10−8. The second term in Eq. (11) corresponds to 5.2 × 10−8·L. For measurements in the long optical rail, the uncertainty caused by the stability of the rail at a 65 m range is below 0.43 µm. The Abbe error is generally better than 2.0 µm along the 65 m travel. The alignment deviation between the distance meter and the measuring system is well within 5.4 mm in the position of 65 m. Hence the cosine error is below 3.5×10−9 in relative. The combined uncertainty can be expressed as [(3 × 10−3 m)2 + (5.2 × 10−8·L)2 + (0.43 ×10−6 m)2 + (2 ×10−6 m)2 + (3.5 ×10−9·L)2]1/2, which is [(3 mm)2 + (5.2 × 10−8·L)2]1/2. The result could indicate the uncertainty of refractive index is continuously improving with the increase of measured distance. Nevertheless, the contribution of increased distance to the stability of phase slope is more significant in the first term. Experimental results demonstrate the introduction of higher harmonic can improve the uncertainty about phase slope to 9.3 × 10−12 rad/Hz (standard deviation), and the combined uncertainty is also reduced accordingly, i.e., [(25 µm)2 + (5.2 × 10−8·L)2]1/2.

Considering the above-mentioned factors in the open-air environment, the uncertainties of temperature, air pressure and humidity are 39.3 mK, 16.5 Pa and 5.1%, corresponding to the uncertainty of air refractive index of 3.6 × 10−8, 4.5 × 10−8, and 4.4 × 10−8, respectively. The combined uncertainty is 7.5 × 10−8. The uncertainty related to the phase slope is 1.68 × 10−8 rad/Hz (standard deviation) based on fundamental frequency, and the combined uncertainty can be expressed as [(4 mm)2 + (7.5 × 10−8·L)2]1/2. With the help of higher harmonic, the uncertainty is 3.7 × 10−11 rad/Hz (standard deviation) for the first term, and the combined uncertainty is [(99 µm)2 + (7.5 × 10−8·L)2]1/2.

5. Summary

Since the advent of optical frequency comb, researchers in many fields have been attracted by its charm. This paper proposes a simple method, and the distance can be measured via the phase slope of the inter-mode beats. It makes the method break the constraint of periodic ambiguity, and benefit simplification of measuring facilities. The introduction of higher harmonic can also ensure the measuring accuracy effectively. A series of long distance experiments are designed to verify the performances. Besides, the dynamic comb spacing makes an attempt to explore the application potential of frequency comb. The open-air experimental results demonstrate the method has the better practicability and the commercial value. For the further research, our work will inspire the superiority of frequency comb in the fields of ultralong absolute distance measurement and dynamic parameters detection, such as underwater object positioning, vibration measurement and ultra-fast dynamic object measurement.

Funding

National Natural Science Foundation of China (61505140); Beijing Aerospace Xinfeng Machinery Equipment Co., Ltd. Aerospace Testing Technology and Application Laboratory Innovation Fund (QT-2018-004); Natural Science Foundation of Tianjin City (18JCYBJC17100); National Basic Research Program of China (973 Program) (2016YFC1401203); Tianjin Municipal Education Commission (JWK1616).

Acknowledgments

We thank National Institute of Metrology for the technical support.

Disclosures

The authors declare no conflicts of interest.

References

1. T. W. Hänsch, “Nobel Lecture: Passion for precision,” Rev. Mod. Phys. 78(4), 1297–1309 (2006). [CrossRef]  

2. N. R. Newbury, “Searching for applications with a fine-tooth comb,” Nat. Photonics 5(4), 186–188 (2011). [CrossRef]  

3. A. Schliesser, N. Picqué, and T. W. HäNsch, “Mid-infrared frequency combs,” Nat. Photonics 6(7), 440–449 (2012). [CrossRef]  

4. T. Udem, R. Holzwarth, and T. W. Hänsch, “Optical frequency metrology,” Nature 416(6877), 233–237 (2002). [CrossRef]  

5. J. L. Hall, “Optical frequency measurement: 40 years of technology revolutions,” IEEE J. Sel. Top. Quantum Electron. 6(6), 1136–1144 (2000). [CrossRef]  

6. N. Picqué and T. W. Hänsch, “Frequency comb spectroscopy,” Nat. Photonics 13(3), 146–157 (2019). [CrossRef]  

7. A. Marian, M. C. Stowe, J. R. Lawall, D. Felinto, and J. Ye, “United time-frequency spectroscopy for dynamics and global structure,” Science 306(5704), 2063–2068 (2004). [CrossRef]  

8. L. C. Sinclair, F. R. Giorgetta, W. C. Swann, E. Baumann, I. Coddington, and N. R. Newbury, “Optical phase noise from atmospheric fluctuations and its impact on optical time-frequency transfer,” Phys. Rev. A 89(2), 023805 (2014). [CrossRef]  

9. M. G. Zeitouny, M. Cui, A. J. E. M. Janssen, N. Bhattacharya, S. A. van den Berg, and H. P. Urbach, “Time-frequency distribution of interferograms from a frequency comb in dispersive media,” Opt. Express 19(4), 3406–3417 (2011). [CrossRef]  

10. R. A. McCracken, J. M. Charsley, and D. T. Reid, “A decade of astrocombs: recent advances in frequency combs for astronomy [Invited],” Opt. Express 25(13), 15058–15078 (2017). [CrossRef]  

11. P. Roztocki and R. Morandotti, “Astrocombs for extreme-precision spectroscopy,” Nat Astron 3(2), 135–136 (2019). [CrossRef]  

12. S. Kim, “Combs rule,” Nat. Photonics 3(6), 313–314 (2009). [CrossRef]  

13. Y. S. Jang and S. W. Kim, “Distance Measurements Using Mode-Locked Lasers: A Review,” Nanomanuf. Metrol. 1(3), 131–147 (2018). [CrossRef]  

14. S. A. Diddams, “The evolving optical frequency comb [Invited],” J. Opt. Soc. Am. 27(11), B51–B62 (2010). [CrossRef]  

15. J. Zhu, P. Cui, Y. Guo, L. Yang, and J. Lin, “Pulse-to-pulse alignment based on interference fringes and the second-order temporal coherence function of optical frequency combs for distance measurement,” Opt. Express 23(10), 13069–13081 (2015). [CrossRef]  

16. 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(13), 1982–1984 (2009). [CrossRef]  

17. H. Wu, F. Zhang, S. Cao, S. Xing, and X. Qu, “Absolute distance measurement by intensity detection using a mode-locked femtosecond pulse laser,” Opt. Express 22(9), 10380–10397 (2014). [CrossRef]  

18. J. Lee, S. Han, K. Lee, E. Bae, S. Kim, and S. Lee, “Absolute distance measurement by dual-comb interferometry with adjustable synthetic wavelength,” Meas. Sci. Technol. 24(4), 045201 (2013). [CrossRef]  

19. I. Coddington, N. Newbury, and W. Swann, “Dual-comb spectroscopy,” Optica 3(4), 414–426 (2016). [CrossRef]  

20. I. Coddington, W. C. Swann, L. Nenadovic, and N. R. Newbury, “Rapid and precise absolute distance measurements at long range,” Nat. Photonics 3(6), 351–356 (2009). [CrossRef]  

21. M. Yu, Y. Okawachi, A. G. Griffith, N. Picqué, M. Lipson, and A. L. Gaeta, “Silicon-chip-based mid-infrared dual-comb spectroscopy,” Nat. Commun. 9(1), 1869 (2018). [CrossRef]  

22. K. E. Webb, J. K. Jang, J. Anthony, S. Coen, M. Erkintalo, and S. G. Murdoch, “Measurement of microresonator frequency comb coherence by spectral interferometry,” Opt. Lett. 41(2), 277–280 (2016). [CrossRef]  

23. K. Joo and S. Kim, “Absolute distance measurement by dispersive interferometry using a femtosecond pulse laser,” Opt. Express 14(13), 5954–5960 (2006). [CrossRef]  

24. X. Zhai, Z. Meng, H. Zhang, X. Xu, Z. Qian, B. Xue, and H. Wu, “Underwater distance measurement using frequency comb laser,” Opt. Express 27(5), 6757–6769 (2019). [CrossRef]  

25. L. Liao, W. Yi, Z. Yang, and G. Wu, “Synthetic-wavelength based absolute distance measurement using heterodyne interferometry of a femtosecond laser,” Acta Phys. Sin. 65(14), 140601 (2016). [CrossRef]  

26. S. Kray, F. Spöler, M. Först, and H. Kurz, “Dual femtosecond laser multiheterodyne optical coherence tomography,” Opt. Lett. 33(18), 2092–2094 (2008). [CrossRef]  

27. J. Jin, Y. Kim, Y. Kim, S. Kim, and C. Kang, “Absolute length calibration of gauge blocks using optical comb of a femtosecond pulse laser,” Opt. Express 14(13), 5968–5974 (2006). [CrossRef]  

28. G. Wang, Y. Jang, S. Hyun, B. Chun, H. Kang, S. Yan, S. Kim, and Y. Kim, “Absolute positioning by multi-wavelength interferometry referenced to the frequency comb of a femtosecond laser,” Opt. Express 23(7), 9121–9129 (2015). [CrossRef]  

29. N. R. Doloca, K. Meiners-Hagen, M. Wedde, F. Pollinger, and A. Abou-Zeid, “Absolute distance measurement system using a femtosecond laser as a modulator,” Meas. Sci. Technol. 21(11), 115302 (2010). [CrossRef]  

30. X. Zhang, W. Yi, M. Hu, Z. Yang, and G. Wu, “Large-scale absolute distance measurement using inter-mode beat of a femtosecond laser,” Acta Phys. Sin. 65(8), 080602 (2016). [CrossRef]  

31. P. E. Ciddor, “Refractive index of air: new equations for the visible and near infrared,” Appl. Opt. 35(9), 1566–1573 (1996). [CrossRef]  

32. X. Jia, Z. Liu, Z. Deng, W. Deng, Z. Wang, and Z. Zhen, “Dynamic absolute distance measurement by frequency sweeping interferometry based Doppler beat frequency tracking model,” Opt. Commun. 430, 163–169 (2019). [CrossRef]  

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 (16)

Fig. 1.
Fig. 1. Structure of measuring system
Fig. 2.
Fig. 2. Schematic of homemade frequency comb. LD, laser diode; WDM, wavelength division multiplexer; PZT, piezoelectric transducer; QWP, quarter-wave plate; HWP, half-wave plate; PBS, polarization beam splitter; ISO, isolator.
Fig. 3.
Fig. 3. (a) Pulse train with 103.6 MHz repetition frequency. (b) Spectrum of the laser source.
Fig. 4.
Fig. 4. (a) Repetition rate variation. (b) Allan deviation of the relative repetition rate. The counter (53230A, Keysight) is synchronized to the Rb clock, whose gate time is 0.1 s.
Fig. 5.
Fig. 5. Nonlinear phase error with the input voltage to the PZT. The left Y axis is the input voltage to the PZT, and the right Y axis is the nonlinear phase error.
Fig. 6.
Fig. 6. Schematic overview of the measuring setup
Fig. 7.
Fig. 7. Experimental photograph on the long rail.
Fig. 8.
Fig. 8. Experimental results for the fundamental frequency. (a) Sweeping of repetition rate for the experiment. (b) Measured data of unwrapped phase corresponding to 5 m, 35 m, 65 m.
Fig. 9.
Fig. 9. Data processing by using EMD. It is clear that the phase jitter is greatly suppressed after EMD.
Fig. 10.
Fig. 10. Experimental results obtained by the fundamental frequency signal.
Fig. 11.
Fig. 11. Experimental results of the 9th harmonic. (a) 9th harmonic of repetition rate sweeping. (b) Experimental data of unwrapped phase change for 5 m, 35 m, 65 m based on the 9th harmonic.
Fig. 12.
Fig. 12. Experimental results based on the 9th harmonic signal at 65 m range.
Fig. 13.
Fig. 13. (a) Experimental photograph on the base line at 1:00 am. (b) Top view of the measurement path.
Fig. 14.
Fig. 14. Experimental results corresponding to 219 m. (a) Phase changes based on the fundamental frequency. (b) Phase changes based on the 9th harmonic.
Fig. 15.
Fig. 15. Distance measurement results over a travel of 219 m. (a) Experimental results based on the fundamental frequency. (b) Experimental results based on the 9th harmonic.
Fig. 16.
Fig. 16. Long-term experiments of the distance measurement at 219 m. (a) Measured residuals for fundamental frequency. (b) Measured residuals based on the higher harmonic. (c) Allan deviation of the distance measurement versus averaging time.

Equations (12)

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

ϕ r e f = 2 π r λ
ϕ m e a s = 2 π ( r + 2 L ) λ
ϕ m e a s = 2 π r λ + 4 π L N g c f r e p
ϕ = 4 π L N g c f r e p
L = ϕ f r e p c 4 π N g
4 π L N g c f r e p 1 4 π L N g c f r e p 2 < 2 π
L N A R = c 2 N g Δ f r e p
ϕ m = 4 π L N g m c f r e p
ϕ m f r e p = 4 π L N g m c
L = ϕ m f r e p 1 m c 4 π N g
ϕ e r r o r = 4 π L N g c D f r e p
u L = ( L S u S ) 2 + ( L N g u N g ) 2 = ( 1 m c 4 π N g u S ) 2 + ( L N g u N g ) 2
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.