# Studying fundamental limit of optical fiber links to $10^{-21}$ level

Dan Xu, Won-Kyu Lee, Fabio Stefani, Olivier Lopez, Anne Amy-Klein, And Paul-Eric Pottie

► To cite this version:

Dan Xu, Won-Kyu Lee, Fabio Stefani, Olivier Lopez, Anne Amy-Klein, et al., Studying fundamental limit of optical fiber links to $10^{-21}$ level. Optics Express, 2018, 26 (8), pp.9515-9527. 10.1364/OE.26.009515. hal-02283571

HAL Id: hal-02283571

https://hal.science/hal-02283571

Submitted on 8 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Studying fundamental limit of optical fiber links to $10^{-21}$ level

Dan Xu, Won-Kyu Lee, Fabio Stefani, Olivier Lopez, Anne Amy-Klein, And Paul-Eric Pottie

► To cite this version:

Dan Xu, Won-Kyu Lee, Fabio Stefani, Olivier Lopez, Anne Amy-Klein, et al. Studying fundamental limit of optical fiber links to $10^{-21}$ level. Opt.Express, 2018, 26 (8), pp.9515-9527. 10.1364/OE.26.009515. hal-02283571v2

HAL Id: hal-02283571

https://hal.science/hal-02283571v2

Submitted on 8 Jun 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Studying the fundamental limit of optical fiber links to the 10⁻²¹ level

DAN XU,<fnref n="1" />,* WON-KYU LEE,<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> FABIO STEFANI,<fnref n="1" /> OLIVIER LOPEZ,<fnref n="4" /> ANNE AMY-KLEIN,<fnref n="4" /> AND PAUL-ERIC POTTIE<fnref n="1" />

**Abstract:** We present a hybrid fiber link combining effective optical frequency transfer and evaluation of performances with a self-synchronized two-way comparison. It enables us to detect the round-trip fiber noise and each of the forward and backward one-way fiber noises simultaneously. The various signals acquired with this setup allow us to study quantitatively several properties of optical fiber links. We check the reciprocity of the accumulated noise forth and back over a bi-directional fiber to the level of 3.1(±3.9) × 10⁻²⁰ based on a 160000s continuous data. We also analyze the noise correlation between two adjacent fibers and show the first experimental evidence of interferometric noise at very low Fourier frequency. We estimate redundantly and consistently the stability and accuracy of the transferred optical frequency over 43 km at 4 × 10⁻²¹ level after 16 days of integration and demonstrate that a frequency comparison with instability as low as 8 × 10⁻¹⁸ would be achievable with uni-directional fibers in urban area.

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

OCIS codes: (060.2360) Fiber optics links and subsystems; (120.3930) Metrological instrumentation; (120.3940) Metrology; (120.5050) Phase measurement.

## References and links

1. N. R. Newbury, P. A. Williams, and W. C. Swann, "Coherent transfer of an optical carrier over 251km," Opt. Lett. 32(21), 3056-3058 (2007).

2. 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).

3. T. L. Nicholson, S. L. Campbell, R. B. Hutson, G. E. Marti, B. J. Bloom, R. L. McNally, W. Zhang, M. D. Barrett, M. S. Safronova, G. F. Strouse, W. L. Tew, and J. Ye, "Systematic evaluation of an atomic clock at 2 × 10⁻¹⁸ total uncertainty," Nat. Commun. 6, 6896 (2015).

4. N. Huntemann, C. Sanner, B. Lipphardt, Chr. Tamm, and E. Peik, "Single-Ion atomic clock with 3 × 10⁻¹⁸ systematic uncertainty," Phys. Rev. Lett. 116(6), 063001 (2016).

5. I. Ushima, M. Takamoto, M. Das, T. Ohkubo, and H. Katori, "Cryogenic optical lattice clocks," Nat. Photon. 9(3), 185-189 (2015).

6. M. Schioppo, R. C. Brown, W. F. McGrew, N. Hinkley, R. J. Fasano, K. Beloy, T. H. Yoon, G. Milani, D. Nicolodi, J. A. Sherman, N. B. Phillips, C. W. Oates, and A. D. Ludlow, "Ultra-stable optical clock with two cold-atom ensembles," Nat. Photon. 11(1), 48-52 (2017).

7. S. Falke, N. Lemke, C. Grebing, B. Lipphardt, S. Weyers, V. Gerginov, N. Huntemann, C. Hagemann, A. Al-Masoudi, S. Häfner, S. Vogt, U. Sterr, and C. Lisdat, "A strontium lattice clock with 3 × 10⁻¹⁷ inaccuracy and its frequency," New J. Phys. 16(7), 073023 (2014).

8. N. Nemitz, T. Ohkubo, M. Takamoto, I. Ushijima, M. Das, N. Ohmae, and H. Katori, "Frequency ratio of Yb and Sr clocks with 5 × 10⁻¹⁷ uncertainty at 150 seconds averaging time," Nat. Photon. 10(4), 258-261 (2016).

9. F. Riehle, "Towards a redefinition of the second based on optical atomic clocks," C. R. Phys. 16(5), 506-515 (2015).

10. S. M. F. Raupach, A. Koczwara, and G. Grosche, "Brillouin amplification supports 1 × 10⁻²⁰ accuracy in optical frequency transfer over 1400 km of underground fibres," Phys. Rev. A 92(2), 021801 (2015).

11. N. Chiodo, N. Quintin, F. Stefani, F. Wiotte, E. Camisard, C. Chardonnet, G. Santarelli, A. Amy-Klein, P.-E. Pottie, and O. Lopez, "Cascaded optical fiber link using the internet network for remote clocks comparison," Opt. Express 23(26), 33927-33937 (2015).

<footnote n="1">¹LNE-SYRTE, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités, UPMC, 61 Avenue de l'Observatoire, 75014 Paris, France</footnote>

<footnote n="2">²Korea Research Institute of Standards and Science, Daejeon 34113, South Korea</footnote>

<footnote n="3">³University of Science and Technology, Yuseong, Daejeon 34113, South Korea</footnote>

<footnote n="4">⁴Laboratoire de Physique des Lasers, Université Paris 13, CNRS, 99 Avenue Jean-Baptiste Clément, 93430 Villetaneuse, France</footnote>

<footnote>*Dan.Xu@obspm.fr</footnote>
# Studying the fundamental limit of optical fiber links to the 10⁻²¹ level

DAN XU,<fnref n="1" />,* WON-KYU LEE,<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> FABIO STEFANI,<fnref n="1" /> OLIVIER LOPEZ,<fnref n="4" /> ANNE AMY-KLEIN,<fnref n="4" /> AND PAUL-ERIC POTTIE<fnref n="1" />

**Abstract:** We present a hybrid fiber link combining effective optical frequency transfer and evaluation of performances with a self-synchronized two-way comparison. It enables us to detect the round-trip fiber noise and each of the forward and backward one-way fiber noises simultaneously. The various signals acquired with this setup allow us to study quantitatively several properties of optical fiber links. We check the reciprocity of the accumulated noise forth and back over a bi-directional fiber to the level of 3.1(±3.9) × 10⁻²⁰ based on a 160000s continuous data. We also analyze the noise correlation between two adjacent fibers and show the first experimental evidence of interferometric noise at very low Fourier frequency. We estimate redundantly and consistently the stability and accuracy of the transferred optical frequency over 43 km at 4 × 10⁻²¹ level after 16 days of integration and demonstrate that a frequency comparison with instability as low as 8 × 10⁻¹⁸ would be achievable with uni-directional fibers in urban area.

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

OCIS codes: (060.2360) Fiber optics links and subsystems; (120.3930) Metrological instrumentation; (120.3940) Metrology; (120.5050) Phase measurement.

## References and links

1. N. R. Newbury, P. A. Williams, and W. C. Swann, "Coherent transfer of an optical carrier over 251km," Opt. Lett. 32(21), 3056-3058 (2007).

2. 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).

3. T. L. Nicholson, S. L. Campbell, R. B. Hutson, G. E. Marti, B. J. Bloom, R. L. McNally, W. Zhang, M. D. Barrett, M. S. Safronova, G. F. Strouse, W. L. Tew, and J. Ye, "Systematic evaluation of an atomic clock at 2 × 10⁻¹⁸ total uncertainty," Nat. Commun. 6, 6896 (2015).

4. N. Huntemann, C. Sanner, B. Lipphardt, Chr. Tamm, and E. Peik, "Single-Ion atomic clock with 3 × 10⁻¹⁸ systematic uncertainty," Phys. Rev. Lett. 116(6), 063001 (2016).

5. I. Ushima, M. Takamoto, M. Das, T. Ohkubo, and H. Katori, "Cryogenic optical lattice clocks," Nat. Photon. 9(3), 185-189 (2015).

6. M. Schioppo, R. C. Brown, W. F. McGrew, N. Hinkley, R. J. Fasano, K. Beloy, T. H. Yoon, G. Milani, D. Nicolodi, J. A. Sherman, N. B. Phillips, C. W. Oates, and A. D. Ludlow, "Ultra-stable optical clock with two cold-atom ensembles," Nat. Photon. 11(1), 48-52 (2017).

7. S. Falke, N. Lemke, C. Grebing, B. Lipphardt, S. Weyers, V. Gerginov, N. Huntemann, C. Hagemann, A. Al-Masoudi, S. Häfner, S. Vogt, U. Sterr, and C. Lisdat, "A strontium lattice clock with 3 × 10⁻¹⁷ inaccuracy and its frequency," New J. Phys. 16(7), 073023 (2014).

8. N. Nemitz, T. Ohkubo, M. Takamoto, I. Ushijima, M. Das, N. Ohmae, and H. Katori, "Frequency ratio of Yb and Sr clocks with 5 × 10⁻¹⁷ uncertainty at 150 seconds averaging time," Nat. Photon. 10(4), 258-261 (2016).

9. F. Riehle, "Towards a redefinition of the second based on optical atomic clocks," C. R. Phys. 16(5), 506-515 (2015).

10. S. M. F. Raupach, A. Koczwara, and G. Grosche, "Brillouin amplification supports 1 × 10⁻²⁰ accuracy in optical frequency transfer over 1400 km of underground fibres," Phys. Rev. A 92(2), 021801 (2015).

11. N. Chiodo, N. Quintin, F. Stefani, F. Wiotte, E. Camisard, C. Chardonnet, G. Santarelli, A. Amy-Klein, P.-E. Pottie, and O. Lopez, "Cascaded optical fiber link using the internet network for remote clocks comparison," Opt. Express 23(26), 33927-33937 (2015).

<footnote n="1">¹LNE-SYRTE, Observatoire de Paris, PSL Research University, CNRS, Sorbonne Universités, UPMC, 61 Avenue de l'Observatoire, 75014 Paris, France</footnote>

<footnote n="2">²Korea Research Institute of Standards and Science, Daejeon 34113, South Korea</footnote>

<footnote n="3">³University of Science and Technology, Yuseong, Daejeon 34113, South Korea</footnote>

<footnote n="4">⁴Laboratoire de Physique des Lasers, Université Paris 13, CNRS, 99 Avenue Jean-Baptiste Clément, 93430 Villetaneuse, France</footnote>

<footnote>*Dan.Xu@obspm.fr</footnote>
12. S. Droste, F. Ozimek, Th. Udem, K. Predeh, T. W. Hänsch, H. Schnatz, G. Grosche, and R. Holzwarth, "Optical-frequency transfer over a single-span 1840 km fiber link," Phys. Rev. Lett. **111**(11), 110801 (2013).

13. L. S. Ma, P. Jungner, J. Ye, and J. L. Hall, "Delivering the same optical frequency at two places: accurate cancellation of phase noise introduced by an optical fiber or other time-varying path," Opt. Lett. **19**(21), 1777-1779 (1994).

14. C. E. Calosso, E. Bertacco, D. Calonico, C. Clivati, G. A. Costanzo, M. Frittelli, F. Levi, A. Mura, and A. Godone, "Frequency transfer via a two-way optical phase comparison on a multiplexed fiber network," Opt. Lett. **39**(5), 1177-1180 (2014).

15. O. Lopez, A. Haboucha, F. Kéfélian, H. Jiang, B. Chanteau, V. Roncin, C. Chardonnet, A. Amy-Klein, and G. Santarelli, "Cascaded multiplexed optical link on a telecommunication network for frequency dissemination," Opt. Express **18**(16), 16849-16857 (2010).

16. A. Bercy, F. Stefani, O. Lopez, Ch. Chardonnet, P. -E. Pottie, and A. Amy-Klein, "Two-way optical frequency comparison over 100 km telecommunication network fibers," Phys. Rev. A **90**(6), 061802 (2014).

17. F. Stefani, O. Lopez, A. Bercy, W. -K. Lee, C. Chardonnet, G. Santarelli, P. -E. Pottie, and A. Amy-Klein, "Tackling the limits of optical fiber links," J. Opt. Soc. Am. B **32**(5), 787-797 (2015).

18. W. -K. Lee, F. Stefani, A. Bercy, O. Lopez, A. Amy-Klein, and P. -E. Pottie, "Hybrid fiber link for accurate optical frequency comparisons," App. Phy. B **123**(5), 161 (2017).

19. J. Grotti, S. Koller, S. Vogt, S. Häfner, U. Sterr, C. Lisdat, H. Denker, C. Voigt, L. Timmen, A. Rolland, F. N. Baynes, H. S. Margolis, M. Zampaolo, P. Thoumany, M. Pizzocaro, B. Rauf, F. Bregolin, A. Tampellini, P. Barbieri, M. Zucco, G. A. Costanzo, C. Clivati, F. Levi, and D. Calonico, "Geodesy and metrology with a transportable optical clock," Nat. Phys. (to be published).

20. A. Tampellini, C. Clivati, F. Levi, A. Mura, and D. Calonico, "Effect of a timebase mismatch in two-way optical frequency transfer," Metrologia **54**(6), 805-809 (2017).

21. P. Delva, J. Lodewyck, S. Bilicki, E. Bookjans, G. Vallet, R. L. Targat, P. -E. Pottie, C. Guerlin, F. Meynadier, C. L. Poncin-Lafitte, O. Lopez, A. Amy-Klein, W. -K. Lee, N. Quintin, C. Lisdat, A. Al-Masoudi, S. Dörscher, C. Grebing, G. Grosche, A. Kuhl, S. Raupach, U. Sterr, I. R. Hill, R. Hobson, W. Bowden, J. Kronjäger, G. Marra, A. Rolland, F. N. Baynes, H. S. Margolis, and P. Gill, "Test of special relativity using a fiber network of optical clocks," Phys. Rev. Lett. **118**(22), 221102 (2017).

22. B. Argence, B. Chanteau, O. Lopez, D. Nicolodi, M. Abgrall, Ch. Chardonnet, C. Dossy, B. Darquié, Y. L. Coq, and A. Amy-Klein, "Quantum cascade laser frequency stabilization at the sub-Hz level," Nat. Photon. **9**(7), 456-460 (2015).

23. P. A. Williams, W. C. Swann, and N. R. Newbury, "High-stability transfer of an optical frequency over long fiber-optic links," J. Opt. Soc. Am. B **25**(8), 1284-1293 (2008).

24. K. Turza, A. Binczewski, W. Bogacki, P. Krehlik, and L. Śliwczyński, "Time and frequency transfer in modern DWDM telecommunication networks," in *Proceedings of IEEE Joint Conference on European Frequency and Time Forum and International Frequency Control Symposium* (IEEE, 2017), pp. 368-370.

25. S. T. Dawkins, J. J. McFerran, and A. N. Luiten, "Considerations on the measurement of the stability of oscillators with frequency counters," IEEE Trans. Ultrason. Ferroelect. Freq. Control **54**(5), 918-925 (2007).

26. W. -K. Lee, D. -H. Yu, C. Y. Park, and J. C. Mun, "The uncertainty associated with the weighted mean frequency of a phase-stabilized signal with white phase noise," Metrologia **47**(1), 24-32 (2010).

27. E. Benkler, C. Lisdat, and U. Sterr, "On the relation between uncertainties of weighted frequency averages and the various types of Allan deviations," Metrologia **52**(4), 565-574 (2015).

28. I. Sesia, E. Cantoni, A. Cernigliaro, G. Signorile, G. Fantino, and P. Tavella, "An efficient and configurable preprocessing algorithm to improve stability analysis," IEEE Trans. Ultrason. Ferroelect. Freq. Control **63**(4), 575-581 (2016).

29. P. Defraigne, N. Guyennon, and C. Bruyninx, "GPS time and frequency transfer: PPP and phase-only analysis," Int. J. Navig. Obs. **2008**, 175468 (2008).

30. S. M. F. Raupach, A. Koczwara, and G. Grosche, "Optical frequency transfer via a 660 km underground fiber link using a remote Brillouin amplifier," Opt. Express **22**(22), 26537-26547 (2014).

31. OFTEN, "Optical frequency transfer - a European network," http://www.ptb.de/emrp/often_home.html.

# 1. Introduction

Optical fiber links have been developed rapidly for time dissemination and frequency transfer over the last decade [1, 2]. They are the most efficient technique for the distant comparison of the rapidly progressing optical atomic clocks [3–8]. These impressive progresses in optical frequency metrology pave the way for a future redefinition of the second [9].

Optical frequency transfer via fiber networks has been performed with ever-increasing performance and fiber lengths. It demonstrates outstanding capabilities to transfer optical frequency with instability below $10^{-15}$ to the low $10^{-20}$ level, even on long-haul fiber links up to 1840 km [10–12]. Meantime, different schemes of optical fiber links have been investigated: the
12. S. Droste, F. Ozimek, Th. Udem, K. Predeh, T. W. Hänsch, H. Schnatz, G. Grosche, and R. Holzwarth, "Optical-frequency transfer over a single-span 1840 km fiber link," Phys. Rev. Lett. **111**(11), 110801 (2013).

13. L. S. Ma, P. Jungner, J. Ye, and J. L. Hall, "Delivering the same optical frequency at two places: accurate cancellation of phase noise introduced by an optical fiber or other time-varying path," Opt. Lett. **19**(21), 1777-1779 (1994).

14. C. E. Calosso, E. Bertacco, D. Calonico, C. Clivati, G. A. Costanzo, M. Frittelli, F. Levi, A. Mura, and A. Godone, "Frequency transfer via a two-way optical phase comparison on a multiplexed fiber network," Opt. Lett. **39**(5), 1177-1180 (2014).

15. O. Lopez, A. Haboucha, F. Kéfélian, H. Jiang, B. Chanteau, V. Roncin, C. Chardonnet, A. Amy-Klein, and G. Santarelli, "Cascaded multiplexed optical link on a telecommunication network for frequency dissemination," Opt. Express **18**(16), 16849-16857 (2010).

16. A. Bercy, F. Stefani, O. Lopez, Ch. Chardonnet, P. -E. Pottie, and A. Amy-Klein, "Two-way optical frequency comparison over 100 km telecommunication network fibers," Phys. Rev. A **90**(6), 061802 (2014).

17. F. Stefani, O. Lopez, A. Bercy, W. -K. Lee, C. Chardonnet, G. Santarelli, P. -E. Pottie, and A. Amy-Klein, "Tackling the limits of optical fiber links," J. Opt. Soc. Am. B **32**(5), 787-797 (2015).

18. W. -K. Lee, F. Stefani, A. Bercy, O. Lopez, A. Amy-Klein, and P. -E. Pottie, "Hybrid fiber link for accurate optical frequency comparisons," App. Phy. B **123**(5), 161 (2017).

19. J. Grotti, S. Koller, S. Vogt, S. Häfner, U. Sterr, C. Lisdat, H. Denker, C. Voigt, L. Timmen, A. Rolland, F. N. Baynes, H. S. Margolis, M. Zampaolo, P. Thoumany, M. Pizzocaro, B. Rauf, F. Bregolin, A. Tampellini, P. Barbieri, M. Zucco, G. A. Costanzo, C. Clivati, F. Levi, and D. Calonico, "Geodesy and metrology with a transportable optical clock," Nat. Phys. (to be published).

20. A. Tampellini, C. Clivati, F. Levi, A. Mura, and D. Calonico, "Effect of a timebase mismatch in two-way optical frequency transfer," Metrologia **54**(6), 805-809 (2017).

21. P. Delva, J. Lodewyck, S. Bilicki, E. Bookjans, G. Vallet, R. L. Targat, P. -E. Pottie, C. Guerlin, F. Meynadier, C. L. Poncin-Lafitte, O. Lopez, A. Amy-Klein, W. -K. Lee, N. Quintin, C. Lisdat, A. Al-Masoudi, S. Dörscher, C. Grebing, G. Grosche, A. Kuhl, S. Raupach, U. Sterr, I. R. Hill, R. Hobson, W. Bowden, J. Kronjäger, G. Marra, A. Rolland, F. N. Baynes, H. S. Margolis, and P. Gill, "Test of special relativity using a fiber network of optical clocks," Phys. Rev. Lett. **118**(22), 221102 (2017).

22. B. Argence, B. Chanteau, O. Lopez, D. Nicolodi, M. Abgrall, Ch. Chardonnet, C. Dossy, B. Darquié, Y. L. Coq, and A. Amy-Klein, "Quantum cascade laser frequency stabilization at the sub-Hz level," Nat. Photon. **9**(7), 456-460 (2015).

23. P. A. Williams, W. C. Swann, and N. R. Newbury, "High-stability transfer of an optical frequency over long fiber-optic links," J. Opt. Soc. Am. B **25**(8), 1284-1293 (2008).

24. K. Turza, A. Binczewski, W. Bogacki, P. Krehlik, and L. Śliwczyński, "Time and frequency transfer in modern DWDM telecommunication networks," in *Proceedings of IEEE Joint Conference on European Frequency and Time Forum and International Frequency Control Symposium* (IEEE, 2017), pp. 368-370.

25. S. T. Dawkins, J. J. McFerran, and A. N. Luiten, "Considerations on the measurement of the stability of oscillators with frequency counters," IEEE Trans. Ultrason. Ferroelect. Freq. Control **54**(5), 918-925 (2007).

26. W. -K. Lee, D. -H. Yu, C. Y. Park, and J. C. Mun, "The uncertainty associated with the weighted mean frequency of a phase-stabilized signal with white phase noise," Metrologia **47**(1), 24-32 (2010).

27. E. Benkler, C. Lisdat, and U. Sterr, "On the relation between uncertainties of weighted frequency averages and the various types of Allan deviations," Metrologia **52**(4), 565-574 (2015).

28. I. Sesia, E. Cantoni, A. Cernigliaro, G. Signorile, G. Fantino, and P. Tavella, "An efficient and configurable preprocessing algorithm to improve stability analysis," IEEE Trans. Ultrason. Ferroelect. Freq. Control **63**(4), 575-581 (2016).

29. P. Defraigne, N. Guyennon, and C. Bruyninx, "GPS time and frequency transfer: PPP and phase-only analysis," Int. J. Navig. Obs. **2008**, 175468 (2008).

30. S. M. F. Raupach, A. Koczwara, and G. Grosche, "Optical frequency transfer via a 660 km underground fiber link using a remote Brillouin amplifier," Opt. Express **22**(22), 26537-26547 (2014).

31. OFTEN, "Optical frequency transfer - a European network," http://www.ptb.de/emrp/often_home.html.

# 1. Introduction

Optical fiber links have been developed rapidly for time dissemination and frequency transfer over the last decade [1, 2]. They are the most efficient technique for the distant comparison of the rapidly progressing optical atomic clocks [3–8]. These impressive progresses in optical frequency metrology pave the way for a future redefinition of the second [9].

Optical frequency transfer via fiber networks has been performed with ever-increasing performance and fiber lengths. It demonstrates outstanding capabilities to transfer optical frequency with instability below $10^{-15}$ to the low $10^{-20}$ level, even on long-haul fiber links up to 1840 km [10–12]. Meantime, different schemes of optical fiber links have been investigated: the
active noise compensation (ANC) method [13] and the passive noise rejection methods called as Two-Way (TW) [14]. In ANC method, the propagation noise is measured after a round-trip and actively compensated with an opto-electronic setup [15]. On the contrary, TW comparison uses passive noise rejection after post-processing.

The TW comparison method can be divided further into two types: classical two-way (CTW) and local two-way (LTW) [16–18]. In CTW method, optical frequency comparisons are obtained by post-processing the data acquired from two physically separated locations, and efforts need to be paid to synchronize the acquisition of the data at two ends of the link [14, 19, 20]. In LTW method, two beat notes are recorded in either of the laboratories: in addition to the "classical" data, one can record another optical beat note between the laser and itself after a round-trip, which is the same as the one detected in an ANC setup. In this case, by processing the data acquired locally, either only at the local site or only at the remote site, the propagation noise can be rejected in real time. The advantage is that the data are self-synchronized and optimal fiber noise rejection is achieved without extra efforts.

We investigate here a hybrid architecture, using an ANC setup on one fiber and using an LTW setup on the other parallel fiber. Such a topology has been reported in [18] and used to analyze the contribution of laser and interferometric noise to stabilities. It is also used for optical frequency comparison between SYRTE (France) and NPL (UK) [21]. We present here an improved hybrid setup with the motivation to use the parallel fibers for effective frequency transfer for direct application in optical metrology and served as a test bed for studying fundamental limits of fiber links. It enables us to detect round-trip and one-way fiber noise, from which we can infer several perfectly synchronous observables and show greater reliability and robustness of the setup. It enables us to assess the performance of optical frequency transfer in real time, while allowing extremely accurate frequency comparison of the primary and secondary frequency standards between National Metrology Institutes (NMIs), and providing effective outputs for advanced molecular spectroscopy experiment [21, 22]. Moreover, for the first time to the best of our knowledge, we can check the reciprocity of the accumulated noise forth and back over a bidirectional fiber, and we can measure the noise floor of an effective frequency transfer setup in situ. With this hybrid setup, we can also evaluate the noise correlation between two adjacent fibers or between forth and back propagation in one fiber.

## 2. Experimental scheme

The hybrid experimental scheme we implement for optical frequency transfer and comparison is shown in Fig. 1. The two fibers, denoted by F1 and F2, are a pair of urban fibers of 43 km connecting two laboratories, SYRTE and LPL. An ANC setup is built on F1; an LTW setup

<figure><img src="image_3.png" /><figcaption>Fig. 1. Experimental scheme for the hybrid fiber links. AOM, acousto-optic modulator; PD, photodiode; PLL, phase-locked loop; bi-EDFA, bi-directional erbium-doped fiber amplifier; p-FM, partial Faraday mirror; PC, polarization controller.</figcaption></figure>
active noise compensation (ANC) method [13] and the passive noise rejection methods called as Two-Way (TW) [14]. In ANC method, the propagation noise is measured after a round-trip and actively compensated with an opto-electronic setup [15]. On the contrary, TW comparison uses passive noise rejection after post-processing.

The TW comparison method can be divided further into two types: classical two-way (CTW) and local two-way (LTW) [16–18]. In CTW method, optical frequency comparisons are obtained by post-processing the data acquired from two physically separated locations, and efforts need to be paid to synchronize the acquisition of the data at two ends of the link [14, 19, 20]. In LTW method, two beat notes are recorded in either of the laboratories: in addition to the "classical" data, one can record another optical beat note between the laser and itself after a round-trip, which is the same as the one detected in an ANC setup. In this case, by processing the data acquired locally, either only at the local site or only at the remote site, the propagation noise can be rejected in real time. The advantage is that the data are self-synchronized and optimal fiber noise rejection is achieved without extra efforts.

We investigate here a hybrid architecture, using an ANC setup on one fiber and using an LTW setup on the other parallel fiber. Such a topology has been reported in [18] and used to analyze the contribution of laser and interferometric noise to stabilities. It is also used for optical frequency comparison between SYRTE (France) and NPL (UK) [21]. We present here an improved hybrid setup with the motivation to use the parallel fibers for effective frequency transfer for direct application in optical metrology and served as a test bed for studying fundamental limits of fiber links. It enables us to detect round-trip and one-way fiber noise, from which we can infer several perfectly synchronous observables and show greater reliability and robustness of the setup. It enables us to assess the performance of optical frequency transfer in real time, while allowing extremely accurate frequency comparison of the primary and secondary frequency standards between National Metrology Institutes (NMIs), and providing effective outputs for advanced molecular spectroscopy experiment [21, 22]. Moreover, for the first time to the best of our knowledge, we can check the reciprocity of the accumulated noise forth and back over a bidirectional fiber, and we can measure the noise floor of an effective frequency transfer setup in situ. With this hybrid setup, we can also evaluate the noise correlation between two adjacent fibers or between forth and back propagation in one fiber.

## 2. Experimental scheme

The hybrid experimental scheme we implement for optical frequency transfer and comparison is shown in Fig. 1. The two fibers, denoted by F1 and F2, are a pair of urban fibers of 43 km connecting two laboratories, SYRTE and LPL. An ANC setup is built on F1; an LTW setup

<figure><img src="image_2.png" /><figcaption>Fig. 1. Experimental scheme for the hybrid fiber links. AOM, acousto-optic modulator; PD, photodiode; PLL, phase-locked loop; bi-EDFA, bi-directional erbium-doped fiber amplifier; p-FM, partial Faraday mirror; PC, polarization controller.</figcaption></figure>
is built on F2 as introduced in [18]. The fibers are fed with one ultra-stable laser, de-drifted by referencing it to an H-Maser using an optical frequency comb. It is split into two parts: one part is injected into F1 (see upper part of Fig. 1); the other is injected into F2 (see lower part of Fig. 1). Two couplers at the local site are used to build the strongly unbalanced Michelson interferometers required for noise compensation. The Michelson interferometers are housed in an aluminum box enclosed in a polyurethane foam box for thermal and acoustic isolation and passive stabilization. Three bi-directional erbium-doped fiber amplifiers (bi-EDFAs) are used to compensate the transmission loss of the ultra-stable laser. A partial Faraday mirror (p-FM) is used at LPL site, which makes a perfect fiber-length match in the remote setup as the uncommon fiber length is zero [18]. It enables the implementation of a simple but robust and temperature-insensitive interferometer at the end-user site. The 50/50 coupler and PD3 represented in dashed line in Fig. 1 lead to 6 dB additional loss. They were used in earlier work, but not effectively used in this experiment. The beat notes both on PD1 and PD2 at SYRTE, after detection, electronic amplification and filtering, are simultaneously recorded by two dead-time free frequency counters (K+K FXE) operated in Π-type and Λ-type with a gate time of 1 s.

Based on the above hybrid scheme, we demonstrate two approaches of two-way frequency comparison. The first one uses two parallel fibers from the same cable, F1 and F2, each fiber transmitting the light in a single direction, referred to as two-way uni-directional (TWU). This approach enables us to estimate the noise correlation between the two adjacent fibers. It opens the way to frequency comparison over telecommunication network with minimal modification of the network backbone. In addition, we realize a second approach which uses a single fiber, F2, through which the light propagates in both directions, referred to as two-way bi-directional (TWB). The TWB approach enables us to detect round-trip and one-way fiber noises and also check the reciprocity of the accumulated noise forth and back on F2.

### 3. Experimental results

At SYRTE, we obtain four beat notes: one on PD1 and three on PD2. On PD1, one records the beat note between the laser beam reflected by FM labeled A and by p-FM labeled D, shown in Fig. 1 as green filled circles. This signal exhibits the round-trip propagation noise on F1. It is detected, tracked and divided to apply an active noise compensation on AOM1. The correction of

<figure><img src="image_4.png" /><figcaption>Fig. 2. Measurements of frequency data acquired only at SYRTE: half the round-trip noise of F1, round-trip noise, forward and backward one-way noise of F2 (Λ-counting).</figcaption></figure>
is built on F2 as introduced in [18]. The fibers are fed with one ultra-stable laser, de-drifted by referencing it to an H-Maser using an optical frequency comb. It is split into two parts: one part is injected into F1 (see upper part of Fig. 1); the other is injected into F2 (see lower part of Fig. 1). Two couplers at the local site are used to build the strongly unbalanced Michelson interferometers required for noise compensation. The Michelson interferometers are housed in an aluminum box enclosed in a polyurethane foam box for thermal and acoustic isolation and passive stabilization. Three bi-directional erbium-doped fiber amplifiers (bi-EDFAs) are used to compensate the transmission loss of the ultra-stable laser. A partial Faraday mirror (p-FM) is used at LPL site, which makes a perfect fiber-length match in the remote setup as the uncommon fiber length is zero [18]. It enables the implementation of a simple but robust and temperature-insensitive interferometer at the end-user site. The 50/50 coupler and PD3 represented in dashed line in Fig. 1 lead to 6 dB additional loss. They were used in earlier work, but not effectively used in this experiment. The beat notes both on PD1 and PD2 at SYRTE, after detection, electronic amplification and filtering, are simultaneously recorded by two dead-time free frequency counters (K+K FXE) operated in Π-type and Λ-type with a gate time of 1 s.

Based on the above hybrid scheme, we demonstrate two approaches of two-way frequency comparison. The first one uses two parallel fibers from the same cable, F1 and F2, each fiber transmitting the light in a single direction, referred to as two-way uni-directional (TWU). This approach enables us to estimate the noise correlation between the two adjacent fibers. It opens the way to frequency comparison over telecommunication network with minimal modification of the network backbone. In addition, we realize a second approach which uses a single fiber, F2, through which the light propagates in both directions, referred to as two-way bi-directional (TWB). The TWB approach enables us to detect round-trip and one-way fiber noises and also check the reciprocity of the accumulated noise forth and back on F2.

### 3. Experimental results

At SYRTE, we obtain four beat notes: one on PD1 and three on PD2. On PD1, one records the beat note between the laser beam reflected by FM labeled A and by p-FM labeled D, shown in Fig. 1 as green filled circles. This signal exhibits the round-trip propagation noise on F1. It is detected, tracked and divided to apply an active noise compensation on AOM1. The correction of

<figure><img src="image_3.png" /><figcaption>Fig. 2. Measurements of frequency data acquired only at SYRTE: half the round-trip noise of F1, round-trip noise, forward and backward one-way noise of F2 (Λ-counting).</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>Fig. 3. Measurement of three sets of TWU and three sets of TWB by combining the four data sets: ANC, RT, OWB, OWF (Λ-counting).</figcaption></figure>

half the round-trip fiber noise on F1 is almost equal to one-way noise. We label this correction ANC (for active noise compensation signal). After active noise compensation, the optical phase at A is then copied to the remote p-FM labeled D. On PD2, we obtain three beat notes: 1) the beat note between the light and itself after a round-trip on F2, denoted as round-trip term (RT); 2) the beat note between the local light and the virtual ultra-stable laser (the point D on p-FM), which is transferred backward on F2. We denote it the one-way backward term (OWB); 3) the beat note between the local light transferred forward on F2, reflected by p-FM, then transferred backward on F2, and the virtual ultra-stable laser transferred backward on F2. This beat note should be detected at LPL, but due to our experimental arrangement, thanks to the p-FM, the beat note is actually detected at SYRTE. We denote this beat note here one-way forward term (OWF). This term is weaker by about 30 dB. By using 2 bi-EDFAs at SYRTE and 1 bi-EDFA at LPL, OWF can be detected at SYRTE. Using a double-down-conversion stage and narrow RF filtering with bandwidth less than 1 MHz, and tracking and dividing the beat notes, stable counting with almost no cycle slips is achieved. The four beat notes contain abundant information about the two parallel fibers' noise: the one-way noise of F1 (ANC), the round-trip noise of F2 (RT), the one-way forward noise of F2 (OWF), the one-way backward noise of F2 (OWB). Their relative frequency fluctuations with Λ-data are shown in Fig. 2. As expected, the one-way noises of F1 (ANC) and F2 (OWB, OWF), are almost at the same level. The round-trip noise on F2 (RT) is twice the other terms.

Based on the above four basic measurements acquired only at SYRTE, we can obtain 7 observables. Three observables which arise from two separated fibers in the same cable used for the uplink and downlink (F1 and F2), are given by

$$ \left\{ \begin{array}{l} \phi_{\text{TWU1}} = \phi_{\text{ANC}} - \phi_{\text{RT}}/2 \\ \phi_{\text{TWU2}} = \phi_{\text{ANC}} - \phi_{\text{OWB}} \\ \phi_{\text{TWU3}} = \phi_{\text{ANC}} - \phi_{\text{OWF}} \end{array} \right. \qquad (1) $$

These uni-directional observables are of the same nature as reported in [16, 18, 23], or similar architecture in DWDM network as reported in [24]. We interpret it as a virtual uni-directional LTW setup where the forward and backward signals propagate in two different (uni-directional)
<figure><img src="image_4.png" /><figcaption>Fig. 3. Measurement of three sets of TWU and three sets of TWB by combining the four data sets: ANC, RT, OWB, OWF (Λ-counting).</figcaption></figure>

half the round-trip fiber noise on F1 is almost equal to one-way noise. We label this correction ANC (for active noise compensation signal). After active noise compensation, the optical phase at A is then copied to the remote p-FM labeled D. On PD2, we obtain three beat notes: 1) the beat note between the light and itself after a round-trip on F2, denoted as round-trip term (RT); 2) the beat note between the local light and the virtual ultra-stable laser (the point D on p-FM), which is transferred backward on F2. We denote it the one-way backward term (OWB); 3) the beat note between the local light transferred forward on F2, reflected by p-FM, then transferred backward on F2, and the virtual ultra-stable laser transferred backward on F2. This beat note should be detected at LPL, but due to our experimental arrangement, thanks to the p-FM, the beat note is actually detected at SYRTE. We denote this beat note here one-way forward term (OWF). This term is weaker by about 30 dB. By using 2 bi-EDFAs at SYRTE and 1 bi-EDFA at LPL, OWF can be detected at SYRTE. Using a double-down-conversion stage and narrow RF filtering with bandwidth less than 1 MHz, and tracking and dividing the beat notes, stable counting with almost no cycle slips is achieved. The four beat notes contain abundant information about the two parallel fibers' noise: the one-way noise of F1 (ANC), the round-trip noise of F2 (RT), the one-way forward noise of F2 (OWF), the one-way backward noise of F2 (OWB). Their relative frequency fluctuations with Λ-data are shown in Fig. 2. As expected, the one-way noises of F1 (ANC) and F2 (OWB, OWF), are almost at the same level. The round-trip noise on F2 (RT) is twice the other terms.

Based on the above four basic measurements acquired only at SYRTE, we can obtain 7 observables. Three observables which arise from two separated fibers in the same cable used for the uplink and downlink (F1 and F2), are given by

$$ \left\{ \begin{array}{l} \phi_{\text{TWU1}} = \phi_{\text{ANC}} - \phi_{\text{RT}}/2 \\ \phi_{\text{TWU2}} = \phi_{\text{ANC}} - \phi_{\text{OWB}} \\ \phi_{\text{TWU3}} = \phi_{\text{ANC}} - \phi_{\text{OWF}} \end{array} \right. \qquad (1) $$

These uni-directional observables are of the same nature as reported in [16, 18, 23], or similar architecture in DWDM network as reported in [24]. We interpret it as a virtual uni-directional LTW setup where the forward and backward signals propagate in two different (uni-directional)
fibers, thus the TWU observables give the difference between the one-way fiber noise of F1 and half the round-trip, forward and backward one-way fiber noise of F2, respectively. As shown in Fig. 3, TWUs behave consistently which give a first indication of high noise correlation between two adjacent fibers. These correlations have been estimated and differ from 1 at the level of $10^{-4}$ with $\Lambda$-counting phase data. It indicates that the phase noises between two fibers are indeed highly correlated. The main limiting factors of TWUs arise from the remaining uncorrelated fiber noises.

The three following observables which arise from one fiber (F2), are given by

$$ \left\{ \begin{aligned} \phi_{\text{TWB1}} &= \phi_{\text{OWB}} - \phi_{\text{RT}}/2 \\ \phi_{\text{TWB2}} &= \phi_{\text{OWF}} - \phi_{\text{RT}}/2. \\ \phi_{\text{TWB3}} &= (\phi_{\text{OWB}} - \phi_{\text{OWF}})/2 \end{aligned} \right. \qquad (2) $$

We interpret it as a fully bi-directional LTW setup. In our previous setup [18], TWB1 was detected at SYRTE and TWB2 at LPL; TWB3 was obtained by post-processing the data between SYRTE and LPL. In the experimental arrangement here, TWB1, TWB2 and TWB3 are obtained at SYRTE, and all of them are self-synchronized. Note that TWB3 is a classical two-way term although no data exchange is needed. The correlations between half the round-trip, forward and backward one-way fiber noise of F2 were also evaluated and differ from 1 at the level of around $10^{-10}$ with $\Lambda$-counting phase data. The rejection of fiber noise can be almost perfect after some integration, and thus the dominating noise for TWBs is interferometric noise.

In addition, we introduce here a new two-way noise floor (TWNF) observable:

$$ \phi_{\text{TWNF}} = (\phi_{\text{OWB}} + \phi_{\text{OWF}} - \phi_{\text{RT}})/2. \qquad (3) $$

With this definition, the observable should converge to zero. The interferometric noise of OWB, OWF and RT exactly compensate themselves (see Appendix for details), then we expect that the TWNF exhibits interferometric noise free behavior. Such an interferometric noise here specially arise from the fibers of the local interferometers used for the beat notes detection. As it is the sum of three independent signals, it contains also the quadratic sum of the noise floors of the three signals. To the best of our knowledge, this is the first time that the noise floor of an effective frequency transfer setup can be measured in situ.

In the following section, we present first a detailed analysis of the phase evolution and frequency stability for the bi-directional observables. The accuracy and long-term behavior of uni-directional and bi-directional observables are presented in section 5.

# 4. Analysis and discussion of the bi-directional observables

## 4.1. Phase evolution and Fourier analysis

In order to provide a better insight of the differences between these TWBs and TWNF observables, we present their phase evolution in Fig. 4(a). Phase error is actually converted into time error and expressed in fs. The TWBs show the same phase evolution of about 3 fs over one day. Obviously, TWB1 and TWB2 have opposite sign, TWB1 and TWB3 coincide, which are expected from our definitions. Knowing the temperature evolution of the experimental room, the temperature sensitivity of optical fibers and the length mismatch of our interferometric setup from previous work, the phase evolution of TWBs can be attributed to interferometric noise [18].

By contrast, TWNF has a phase evolution less than 0.2 fs, as expected. Figure 4(a) illustrates that the phase wandering seen with TWBs are almost erased in TWNF. This is an evidence that TWNF is interferometric-noise free. This is a new situation compared to the noise floor measurements reported elsewhere previously. Usually noise floor is estimated by optically short-cutting the link with a short fiber patch cord, and setting the detection parameter nearly
fibers, thus the TWU observables give the difference between the one-way fiber noise of F1 and half the round-trip, forward and backward one-way fiber noise of F2, respectively. As shown in Fig. 3, TWUs behave consistently which give a first indication of high noise correlation between two adjacent fibers. These correlations have been estimated and differ from 1 at the level of $10^{-4}$ with $\Lambda$-counting phase data. It indicates that the phase noises between two fibers are indeed highly correlated. The main limiting factors of TWUs arise from the remaining uncorrelated fiber noises.

The three following observables which arise from one fiber (F2), are given by

$$ \left\{ \begin{aligned} \phi_{\text{TWB1}} &= \phi_{\text{OWB}} - \phi_{\text{RT}}/2 \\ \phi_{\text{TWB2}} &= \phi_{\text{OWF}} - \phi_{\text{RT}}/2. \\ \phi_{\text{TWB3}} &= (\phi_{\text{OWB}} - \phi_{\text{OWF}})/2 \end{aligned} \right. \qquad (2) $$

We interpret it as a fully bi-directional LTW setup. In our previous setup [18], TWB1 was detected at SYRTE and TWB2 at LPL; TWB3 was obtained by post-processing the data between SYRTE and LPL. In the experimental arrangement here, TWB1, TWB2 and TWB3 are obtained at SYRTE, and all of them are self-synchronized. Note that TWB3 is a classical two-way term although no data exchange is needed. The correlations between half the round-trip, forward and backward one-way fiber noise of F2 were also evaluated and differ from 1 at the level of around $10^{-10}$ with $\Lambda$-counting phase data. The rejection of fiber noise can be almost perfect after some integration, and thus the dominating noise for TWBs is interferometric noise.

In addition, we introduce here a new two-way noise floor (TWNF) observable:

$$ \phi_{\text{TWNF}} = (\phi_{\text{OWB}} + \phi_{\text{OWF}} - \phi_{\text{RT}})/2. \qquad (3) $$

With this definition, the observable should converge to zero. The interferometric noise of OWB, OWF and RT exactly compensate themselves (see Appendix for details), then we expect that the TWNF exhibits interferometric noise free behavior. Such an interferometric noise here specially arise from the fibers of the local interferometers used for the beat notes detection. As it is the sum of three independent signals, it contains also the quadratic sum of the noise floors of the three signals. To the best of our knowledge, this is the first time that the noise floor of an effective frequency transfer setup can be measured in situ.

In the following section, we present first a detailed analysis of the phase evolution and frequency stability for the bi-directional observables. The accuracy and long-term behavior of uni-directional and bi-directional observables are presented in section 5.

# 4. Analysis and discussion of the bi-directional observables

## 4.1. Phase evolution and Fourier analysis

In order to provide a better insight of the differences between these TWBs and TWNF observables, we present their phase evolution in Fig. 4(a). Phase error is actually converted into time error and expressed in fs. The TWBs show the same phase evolution of about 3 fs over one day. Obviously, TWB1 and TWB2 have opposite sign, TWB1 and TWB3 coincide, which are expected from our definitions. Knowing the temperature evolution of the experimental room, the temperature sensitivity of optical fibers and the length mismatch of our interferometric setup from previous work, the phase evolution of TWBs can be attributed to interferometric noise [18].

By contrast, TWNF has a phase evolution less than 0.2 fs, as expected. Figure 4(a) illustrates that the phase wandering seen with TWBs are almost erased in TWNF. This is an evidence that TWNF is interferometric-noise free. This is a new situation compared to the noise floor measurements reported elsewhere previously. Usually noise floor is estimated by optically short-cutting the link with a short fiber patch cord, and setting the detection parameter nearly
<figure><img src="image_6.png" /><figcaption>Fig. 4. (a) Phase evolution of the bi-directional LTW observables: TWBs and TWNF, (b) phase noise power spectral densities(PSD) of the roundtrip noise on F2 and the residual phase noise by bi-directional LTW method. Note that three PSDs corresponding to TWB1, TWB2 and TWB3 are overlaid on each other.</figcaption></figure>

the same as those in real experimental conditions. In that case, the bandwidth of the servo is inherently changed, as the propagation time delay is modified. The interferometric noise of the setup is kept unmodified. However, in this new situation, it is the opposite: the detection parameters are the very same as the ones of the running link : the propagation time delay is not modified, power levels are the same, but the interferometric noise is erased from the data set combination.

To study further the properties of the TWBs and TWNF, we derive the power spectral density (PSD) by calculating the fast Fourier transform of the auto-correlation of their phase noise [17], as shown in Fig. 4(b). We find out that in our hybrid fiber link, from 0.5 mHz to about 80 mHz, the PSD of RT falls off as a flicker frequency noise $f^{-3}$, where f is the Fourier frequency. Then the RT noise changes to white frequency noise ($f^{-2}$) above 80 mHz. It is interesting to note that the PSD of TWBs also fall off as $f^{-3}$ from 0.5 mHz to 5 mHz. It indicates that the TWBs residual noise are originating from fiber noise arising from non-compensated fiber segments, thus the interferometers’ noises dominate the TWBs noise at low frequencies. For higher Fourier frequency, approximately white phase noises are achieved for TWBs at about $10^{-3}$ rad$^2$/Hz, indicating that the phase noise is canceled within the bandwidth. The PSD of TWNF shows only white phase noise signature over almost the full frequency range. The constant slope even at very low Fourier frequency confirms that TWNF is interferometric-noise free. More informations are detailed in Appendix.

## 4.2. Frequency stability

The fractional frequency instabilities of TWBs and TWNF are plotted in Fig. 5 as modified Allan deviation (MDEV) calculated from the Λ-data [25]. The stability of TWBs is as low as $5 \times 10^{-17}$ at 1 s integration time. The stability of TWB3 is lower than that of TWB1 and TWB2, which could arise from the different detection noises. We expect indeed that the detection noise contribution is higher for TWB1 and TWB2 than for TWB3, for which the noise contribution of the backward-F2 and forward-F2 beat notes are both divided by two, following Eq. (7) in Appendix. At long integration time, the measurements reach the interferometric noise floor, and the three curves of TWBs perfectly merge themselves. The TWBs reach relative frequency stabilities of $1.5 \times 10^{-20}$ at 10 000s integration time. At the longer integration time, the MDEV of TWNF deviates from TWBs and keeps decreasing with a slope of $-\frac{3}{2}$. TWNF reaches $2.6 \times 10^{-22}$
<figure><img src="image_5.png" /><figcaption>Fig. 4. (a) Phase evolution of the bi-directional LTW observables: TWBs and TWNF, (b) phase noise power spectral densities(PSD) of the roundtrip noise on F2 and the residual phase noise by bi-directional LTW method. Note that three PSDs corresponding to TWB1, TWB2 and TWB3 are overlaid on each other.</figcaption></figure>

the same as those in real experimental conditions. In that case, the bandwidth of the servo is inherently changed, as the propagation time delay is modified. The interferometric noise of the setup is kept unmodified. However, in this new situation, it is the opposite: the detection parameters are the very same as the ones of the running link : the propagation time delay is not modified, power levels are the same, but the interferometric noise is erased from the data set combination.

To study further the properties of the TWBs and TWNF, we derive the power spectral density (PSD) by calculating the fast Fourier transform of the auto-correlation of their phase noise [17], as shown in Fig. 4(b). We find out that in our hybrid fiber link, from 0.5 mHz to about 80 mHz, the PSD of RT falls off as a flicker frequency noise $f^{-3}$, where f is the Fourier frequency. Then the RT noise changes to white frequency noise ($f^{-2}$) above 80 mHz. It is interesting to note that the PSD of TWBs also fall off as $f^{-3}$ from 0.5 mHz to 5 mHz. It indicates that the TWBs residual noise are originating from fiber noise arising from non-compensated fiber segments, thus the interferometers’ noises dominate the TWBs noise at low frequencies. For higher Fourier frequency, approximately white phase noises are achieved for TWBs at about $10^{-3}$ rad$^2$/Hz, indicating that the phase noise is canceled within the bandwidth. The PSD of TWNF shows only white phase noise signature over almost the full frequency range. The constant slope even at very low Fourier frequency confirms that TWNF is interferometric-noise free. More informations are detailed in Appendix.

## 4.2. Frequency stability

The fractional frequency instabilities of TWBs and TWNF are plotted in Fig. 5 as modified Allan deviation (MDEV) calculated from the Λ-data [25]. The stability of TWBs is as low as $5 \times 10^{-17}$ at 1 s integration time. The stability of TWB3 is lower than that of TWB1 and TWB2, which could arise from the different detection noises. We expect indeed that the detection noise contribution is higher for TWB1 and TWB2 than for TWB3, for which the noise contribution of the backward-F2 and forward-F2 beat notes are both divided by two, following Eq. (7) in Appendix. At long integration time, the measurements reach the interferometric noise floor, and the three curves of TWBs perfectly merge themselves. The TWBs reach relative frequency stabilities of $1.5 \times 10^{-20}$ at 10 000s integration time. At the longer integration time, the MDEV of TWNF deviates from TWBs and keeps decreasing with a slope of $-3/2$. TWNF reaches $2.6 \times 10^{-22}$
<figure><img src="image_7.png" /><figcaption>Fig. 5. Modified Allan deviation (MDEV) of 3 bi-directional LTW observables TWB1, TWB2, TWB3, as well as TWNF (Λ-counting).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Fig. 6. Fractional frequency stabilities of TWB1 in terms of overlapping Allan deviation (ADEV) and modified Allan deviation (MDEV) with Π-type counter and Λ-type counter.</figcaption></figure>

at 10 000 s, which further indicates that TWNF is insensitive to diurnal noise bump due to temperature variations.

Figure 6 shows the fractional frequency instabilities of bi-directional LTW in terms of overlapping Allan deviation (ADEV) and modified Allan deviation (MDEV). The free running round-trip fiber noise of F2 is shown as a black line. We show only TWB1 for the sake of clarity, but similar data were obtained for TWB2 and TWB3. The ADEV and MDEV of TWB1 that are calculated from the Π-data start from $6.3 \times 10^{-16}$ at 1 s integration time and decrease with the slope of $-1$ and $-\frac{3}{2}$, respectively. The ADEV and MDEV of TWB1 from the Λ-data start from $5.9 \times 10^{-17}$ at 1 s and decrease to $2.1 \times 10^{-20}$ and $1.5 \times 10^{-20}$ at 10 000s, respectively. The MDEVs of the Π-data and Λ-data converge to the same value and show good agreement. The
<figure><img src="image_6.png" /><figcaption>Fig. 5. Modified Allan deviation (MDEV) of 3 bi-directional LTW observables TWB1, TWB2, TWB3, as well as TWNF (Λ-counting).</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Fig. 6. Fractional frequency stabilities of TWB1 in terms of overlapping Allan deviation (ADEV) and modified Allan deviation (MDEV) with Π-type counter and Λ-type counter.</figcaption></figure>

at 10 000 s, which further indicates that TWNF is insensitive to diurnal noise bump due to temperature variations.

Figure 6 shows the fractional frequency instabilities of bi-directional LTW in terms of overlapping Allan deviation (ADEV) and modified Allan deviation (MDEV). The free running round-trip fiber noise of F2 is shown as a black line. We show only TWB1 for the sake of clarity, but similar data were obtained for TWB2 and TWB3. The ADEV and MDEV of TWB1 that are calculated from the Π-data start from $6.3 \times 10^{-16}$ at 1 s integration time and decrease with the slope of $-1$ and $-\frac{3}{2}$, respectively. The ADEV and MDEV of TWB1 from the Λ-data start from $5.9 \times 10^{-17}$ at 1 s and decrease to $2.1 \times 10^{-20}$ and $1.5 \times 10^{-20}$ at 10 000s, respectively. The MDEVs of the Π-data and Λ-data converge to the same value and show good agreement. The
ADEVs of the Π-data and Λ-data show a difference compatible with the statistical error bars. Actually we check that ADEVs converge to the same value at longer integration time, as seen in Table 1. Compared to the results shown in [18], we find excellent agreement and consistent behavior over 2 years of operation. The ratio between the Λ-data and Π-data at 1 s integration time is 10.6. Using K+K counter, the maximum expected ratio for pure white phase noise is approximately 31. The $τ^{-1}$ behavior of the MDEV using Λ-data and the noise floor measurement given by TWNF show that we are limited by the noise floor of the setup. We check this assumption by replacing the fiber link by a short patch fiber and setting an equivalent remote end (AOMs and p-FM) at SYRTE. We observe the same noise floor limitation from 1 s to 100 s integration time.

# 5. Uni-directional and bi-directional data over days

## 5.1. Accuracy and reciprocity

We evaluate the accuracy of the frequency comparison of the hybrid setup. We calculate the arithmetic mean of the difference between the expected and the measured beat frequency of the relevant observables TWUs and TWBs introduced above. The statistical uncertainty is estimated by the instability of the data set at a well-considered averaging time [10, 18, 26, 27]. The effects of missing data in comparisons were explored in the case of GPS and satellite two-ways, but not in the context of coherent fiber links [28, 29]. The question is not trivial, and out of the scope of this paper. So we use here a subset of 160 000 s (2 days) data without any cycle slip, recorded with 1 s gate time. The results for this cycle-slip-free 160 000 s data are presented in Table 1, with the mean value, ADEV and MDEV at 40 000 s of integration time, with Π-data and Λ-data.

The three TWU and three TWB are consistent, which show the reliability of the setup and give a greater confidence into the quality of the data set. We observe the convergence of ADEVs using Π-data and Λ-data. The MDEVs converge also to identical values. No systematic shifts are observed for all the observables within the statistical uncertainty at $1-\sigma$ level.

For uni-directional LTW, we obtain mean value of these relative frequencies of $2.2 \times 10^{-17}$, with a statistical uncertainty of $1.0 \times 10^{-17}$ and $7.0 \times 10^{-18}$ given by the long-term ADEV and MDEV at 40 000 s (Π-data and Λ-data), respectively. For bi-directional LTW, the mean values are ~ $3.0 \times 10^{-20}$, with a statistical uncertainty of $3.9 \times 10^{-20}$ and $2.8 \times 10^{-20}$ given by the long-term ADEV and MDEV at 40 000 s (Π-data and Λ-data), respectively. As a conclusion, the accuracy of our setup is limited by the statistical uncertainty contribution. We consider here then conservatively the worse case, to be $3.9 \times 10^{-20}$.

In particular, the TWB3 can be interpreted as a test of reciprocity using in field fibers between the forward and backward paths, at the level of $3.1(\pm3.9) \times 10^{-20}$. It is limited by the interferometric noise at SYRTE. The effect of a periodic perturbation can be also a problem at such high level of accuracy. Improvement of the setup and a better data processing should allow us to improve the accuracy of this first test of reciprocity.

## 5.2. Behavior over 16 days

The hybrid setup runs continuously for one month, as a part of an international clock comparison between NPL (UK), SYRTE (France) and PTB (Germany). We show in Fig. 7 the instability and phase evolution of the observables of the hybrid fiber link over 16 days. The interruptions seen on the plot at MJD = 57904 + 5 and +6 are not due to the setup, but to an interruption of the ultra-stable light (no input signal). Except for the interruptions, the uptime over 14 days is estimated conservatively to 96.5%. RT has a phase evolution less than 4 ns; TWU2 is at the level of ps; TWB3 is at the level of fs. The MDEV of TWU2 starts from $7.0 \times 10^{-16}$ at 1 s and decreases to $1.3 \times 10^{-17}$ at 400 000s integration time. The MDEV of TWB3 is $3.9 \times 10^{-17}$ at 1 s and reaches as low as $8.5 \times 10^{-22}$ at 400 000s integration time. It shows the reliability of the hybrid setup in the long-term. The stability curve of TWB3 shows a plateau from $10^4$ s to $4 \times 10^4$ s.
ADEVs of the Π-data and Λ-data show a difference compatible with the statistical error bars. Actually we check that ADEVs converge to the same value at longer integration time, as seen in Table 1. Compared to the results shown in [18], we find excellent agreement and consistent behavior over 2 years of operation. The ratio between the Λ-data and Π-data at 1 s integration time is 10.6. Using K+K counter, the maximum expected ratio for pure white phase noise is approximately 31. The $τ^{-1}$ behavior of the MDEV using Λ-data and the noise floor measurement given by TWNF show that we are limited by the noise floor of the setup. We check this assumption by replacing the fiber link by a short patch fiber and setting an equivalent remote end (AOMs and p-FM) at SYRTE. We observe the same noise floor limitation from 1 s to 100 s integration time.

# 5. Uni-directional and bi-directional data over days

## 5.1. Accuracy and reciprocity

We evaluate the accuracy of the frequency comparison of the hybrid setup. We calculate the arithmetic mean of the difference between the expected and the measured beat frequency of the relevant observables TWUs and TWBs introduced above. The statistical uncertainty is estimated by the instability of the data set at a well-considered averaging time [10, 18, 26, 27]. The effects of missing data in comparisons were explored in the case of GPS and satellite two-ways, but not in the context of coherent fiber links [28, 29]. The question is not trivial, and out of the scope of this paper. So we use here a subset of 160 000 s (2 days) data without any cycle slip, recorded with 1 s gate time. The results for this cycle-slip-free 160 000 s data are presented in Table 1, with the mean value, ADEV and MDEV at 40 000 s of integration time, with Π-data and Λ-data.

The three TWU and three TWB are consistent, which show the reliability of the setup and give a greater confidence into the quality of the data set. We observe the convergence of ADEVs using Π-data and Λ-data. The MDEVs converge also to identical values. No systematic shifts are observed for all the observables within the statistical uncertainty at $1-\sigma$ level.

For uni-directional LTW, we obtain mean value of these relative frequencies of $2.2 \times 10^{-17}$, with a statistical uncertainty of $1.0 \times 10^{-17}$ and $7.0 \times 10^{-18}$ given by the long-term ADEV and MDEV at 40 000 s (Π-data and Λ-data), respectively. For bi-directional LTW, the mean values are ~ $3.0 \times 10^{-20}$, with a statistical uncertainty of $3.9 \times 10^{-20}$ and $2.8 \times 10^{-20}$ given by the long-term ADEV and MDEV at 40 000 s (Π-data and Λ-data), respectively. As a conclusion, the accuracy of our setup is limited by the statistical uncertainty contribution. We consider here then conservatively the worse case, to be $3.9 \times 10^{-20}$.

In particular, the TWB3 can be interpreted as a test of reciprocity using in field fibers between the forward and backward paths, at the level of $3.1(\pm3.9) \times 10^{-20}$. It is limited by the interferometric noise at SYRTE. The effect of a periodic perturbation can be also a problem at such high level of accuracy. Improvement of the setup and a better data processing should allow us to improve the accuracy of this first test of reciprocity.

## 5.2. Behavior over 16 days

The hybrid setup runs continuously for one month, as a part of an international clock comparison between NPL (UK), SYRTE (France) and PTB (Germany). We show in Fig. 7 the instability and phase evolution of the observables of the hybrid fiber link over 16 days. The interruptions seen on the plot at MJD = 57904 + 5 and +6 are not due to the setup, but to an interruption of the ultra-stable light (no input signal). Except for the interruptions, the uptime over 14 days is estimated conservatively to 96.5%. RT has a phase evolution less than 4 ns; TWU2 is at the level of ps; TWB3 is at the level of fs. The MDEV of TWU2 starts from $7.0 \times 10^{-16}$ at 1 s and decreases to $1.3 \times 10^{-17}$ at 400 000s integration time. The MDEV of TWB3 is $3.9 \times 10^{-17}$ at 1 s and reaches as low as $8.5 \times 10^{-22}$ at 400 000s integration time. It shows the reliability of the hybrid setup in the long-term. The stability curve of TWB3 shows a plateau from $10^4$ s to $4 \times 10^4$ s.
<table><caption>Table 1. Accuracy of the frequency comparisons using the hybrid architecture</caption><thead><tr><th>Quantity</th><th></th><th>TWU1</th><th>TWU2<br/>×10<sup>-17</sup></th><th>TWU3</th><th>TWB1</th><th>TWB2<br/>×10<sup>-20</sup></th><th>TWB3</th></tr></thead><tbody><tr><td rowspan="3">Π</td><td>Mean</td><td>2.2</td><td>2.2</td><td>2.2</td><td>2.9</td><td>3.2</td><td>3.1</td></tr><tr><td>ADEV at 40 000s</td><td>1.0</td><td>1.0</td><td>1.0</td><td>4.0</td><td>4.0</td><td>3.9</td></tr><tr><td>MDEV at 40 000s</td><td>0.7</td><td>0.7</td><td>0.7</td><td>2.8</td><td>2.8</td><td>2.8</td></tr><tr><td rowspan="3">Λ</td><td>Mean</td><td>2.2</td><td>2.2</td><td>2.2</td><td>3.0</td><td>2.9</td><td>2.9</td></tr><tr><td>ADEV at 40 000s</td><td>1.0</td><td>1.0</td><td>1.0</td><td>3.9</td><td>3.9</td><td>3.9</td></tr><tr><td>MDEV at 40 000s</td><td>0.7</td><td>0.7</td><td>0.7</td><td>2.8</td><td>2.8</td><td>2.8</td></tr></tbody></table>

s integration time, and then a rapid decay with $τ^{-2}$, indicating a periodic perturbation with a time constant of approximately 40 000 s, which we attribute to the diurnal temperature fluctuations of the fiber, which is also evidenced by the instability curve of RT [blue dots in Fig. 7(a)] [30].

The typical stability of Cs Fountain Clock is shown as blue dashed line, taken as $4 \times 10^{-14}/\sqrt{\tau}$; the red dashed line shows the stability of one of the best optical clocks, taken as $8 \times 10^{-17}/\sqrt{\tau}$, as reported for instance in [6]. For achieving optical frequency comparison with the best accuracy, it is mandatory to use fully bi-directional setup. Moreover, the data acquired with our hybrid setup confirm the first results reported in [14], and show the ability of uni-directional setup to compare the best microwave cold atoms frequency standards, at least at a regional scale.

We use the data set of 16 continuous days for accuracy analysis. For TWU2, the mean offset frequency of $8.8 \times 10^{-18}$ shows nearly no deviation within the statistical uncertainty of $1.3 \times 10^{-17}$. For TWB3, we obtain a mean offset frequency of $4.2 \times 10^{-21}$, that deviates from its statistical uncertainty of $8.5 \times 10^{-22}$. This frequency shift could arise from subtle technical issues or from the limited time of observation (as compared to the 1-day periodic perturbation presented in this data set). The effect is not explained so far with simple models and need detailed analysis over very long run, and will be the object of future studies. Anyway, although a conclusion cannot be made about the origin of the effect, it is not a limiting factor for the accuracy of optical clock comparison by *at least* 2 orders of magnitude.

<figure><img src="image_9.png" /><figcaption>Fig. 7. (a) Long-term fractional frequency instabilities and (b) phase evolution of a uni-directional and a bi-directional two-way scheme(Λ-counting).</figcaption></figure>
<table><caption>Table 1. Accuracy of the frequency comparisons using the hybrid architecture</caption><thead><tr><th>Quantity</th><th></th><th>TWU1</th><th>TWU2<br/>×10<sup>-17</sup></th><th>TWU3</th><th>TWB1</th><th>TWB2<br/>×10<sup>-20</sup></th><th>TWB3</th></tr></thead><tbody><tr><td rowspan="3">Π</td><td>Mean</td><td>2.2</td><td>2.2</td><td>2.2</td><td>2.9</td><td>3.2</td><td>3.1</td></tr><tr><td>ADEV at 40 000s</td><td>1.0</td><td>1.0</td><td>1.0</td><td>4.0</td><td>4.0</td><td>3.9</td></tr><tr><td>MDEV at 40 000s</td><td>0.7</td><td>0.7</td><td>0.7</td><td>2.8</td><td>2.8</td><td>2.8</td></tr><tr><td rowspan="3">Λ</td><td>Mean</td><td>2.2</td><td>2.2</td><td>2.2</td><td>3.0</td><td>2.9</td><td>2.9</td></tr><tr><td>ADEV at 40 000s</td><td>1.0</td><td>1.0</td><td>1.0</td><td>3.9</td><td>3.9</td><td>3.9</td></tr><tr><td>MDEV at 40 000s</td><td>0.7</td><td>0.7</td><td>0.7</td><td>2.8</td><td>2.8</td><td>2.8</td></tr></tbody></table>

s integration time, and then a rapid decay with $τ^{-2}$, indicating a periodic perturbation with a time constant of approximately 40 000 s, which we attribute to the diurnal temperature fluctuations of the fiber, which is also evidenced by the instability curve of RT [blue dots in Fig. 7(a)] [30].

The typical stability of Cs Fountain Clock is shown as blue dashed line, taken as $4 \times 10^{-14}/\sqrt{\tau}$; the red dashed line shows the stability of one of the best optical clocks, taken as $8 \times 10^{-17}/\sqrt{\tau}$, as reported for instance in [6]. For achieving optical frequency comparison with the best accuracy, it is mandatory to use fully bi-directional setup. Moreover, the data acquired with our hybrid setup confirm the first results reported in [14], and show the ability of uni-directional setup to compare the best microwave cold atoms frequency standards, at least at a regional scale.

We use the data set of 16 continuous days for accuracy analysis. For TWU2, the mean offset frequency of $8.8 \times 10^{-18}$ shows nearly no deviation within the statistical uncertainty of $1.3 \times 10^{-17}$. For TWB3, we obtain a mean offset frequency of $4.2 \times 10^{-21}$, that deviates from its statistical uncertainty of $8.5 \times 10^{-22}$. This frequency shift could arise from subtle technical issues or from the limited time of observation (as compared to the 1-day periodic perturbation presented in this data set). The effect is not explained so far with simple models and need detailed analysis over very long run, and will be the object of future studies. Anyway, although a conclusion cannot be made about the origin of the effect, it is not a limiting factor for the accuracy of optical clock comparison by *at least* 2 orders of magnitude.

<figure><img src="image_8.png" /><figcaption>Fig. 7. (a) Long-term fractional frequency instabilities and (b) phase evolution of a uni-directional and a bi-directional two-way scheme(Λ-counting).</figcaption></figure>
## 6. Conclusion

We present fully new data and data analysis for the 2 × 43 km hybrid fiber link. We demonstrate a hybrid setup suitable for the study of fundamental limits of fiber links, while being a very effective link for optical frequency dissemination and optical frequency comparison. We show that the hybrid setup evaluates uni-directional LTW abilities to reach instability below 10⁻¹⁷ over 43 km of an urban network. The hybrid setup we have implemented enables us to transfer frequency with a relative frequency stability as low as 3.9 × 10⁻¹⁷ at 1 s integration time and below 10⁻²¹ over 16 days. Combined with another hybrid links between LPL and NPL, this hybrid setup is already a part of an "optical fiber links chain" linking NPL, SYRTE and PTB [31]. With this work, we are showing that far from being competitive, the combination of two techniques can even lead to improve efficiency, redundancy and overall performances.

## Appendix

Here, we analyze in detail the phase noise in the hybrid setup. According to the noise sources, the total phase noise of the beat note can be given by the sum of four terms splitting the phase wandering into the laser term, the RF term, the interferometric term and the fiber term:

$$ \phi = \phi_{\text{laser}} + \phi_{\text{RF}} + \phi_{\text{inter}} + \phi_{\text{fiber}}. \quad (4) $$

Since the same laser is used to inject the two fibers, and its frequency is almost constant, the first term can be neglected and we state $\phi_{\text{laser}} = 0$. By phase-locking the PLL reference synthesizers and the frequency counters to a common 10 MHz reference signal, the RF contribution to phase can be neglected, and we state $\phi_{\text{RF}} = 0$. The interferometric noise is due to the length mismatch in the local interferometers (a perfect match of the fiber length in the remote interferometer is accomplished by the use of a partial Faraday mirror at the remote setup) [18]. As the interferometers are isolated passively in a box with time constant which is much greater than the light propagation time in the fiber, we write the interferometric term as $\phi_{\text{inter}}(t)$. Thus for the four basic beat notes, the two dominant phase noises are the fiber link noise and interferometers' noise

$$ \left\{ \begin{aligned} \phi_{\text{ANC}}(t) &= \phi_{\text{inter-ANC}}(t) + \frac{1}{2}\phi_{\text{RT-F1}}(t) \\ \phi_{\text{RT}}(t) &= \phi_{\text{inter-RT}}(t) + \phi_{\text{RT-F2}}(t) \\ \phi_{\text{OWB}}(t) &= \phi_{\text{inter-OWB}}(t) + \phi_{\text{backward-F2}}(t) \\ \phi_{\text{OWF}}(t) &= \phi_{\text{inter-OWF}}(t) + \phi_{\text{forward-F2}}(t - \tau). \end{aligned} \right. \quad (5) $$

where $\tau$ is the light propagation delay between the local site and the remote site, $\phi_{\text{inter-ANC}}$ is the interferometric noise in the ANC setup, $\phi_{\text{RT-F1}}$ is the roundtrip propagation noise on F1, and other notations have the similar definitions. Then the TWUs can be expressed as

$$ \left\{ \begin{aligned} \phi_{\text{TWU1}}(t) &= \phi_{\text{inter-TWU1}}(t) + \frac{1}{2}\phi_{\text{RT-F1}}(t) - \frac{1}{2}\phi_{\text{RT-F2}}(t) \\ \phi_{\text{TWU2}}(t) &= \phi_{\text{inter-TWU2}}(t) + \frac{1}{2}\phi_{\text{RT-F1}}(t) - \phi_{\text{backward-F2}}(t). \\ \phi_{\text{TWU3}}(t) &= \phi_{\text{inter-TWU3}}(t) + \frac{1}{2}\phi_{\text{RT-F1}}(t) - \phi_{\text{forward-F2}}(t - \tau) \end{aligned} \right. \quad (6) $$

The TWBs can be expressed as
## 6. Conclusion

We present fully new data and data analysis for the 2 × 43 km hybrid fiber link. We demonstrate a hybrid setup suitable for the study of fundamental limits of fiber links, while being a very effective link for optical frequency dissemination and optical frequency comparison. We show that the hybrid setup evaluates uni-directional LTW abilities to reach instability below 10⁻¹⁷ over 43 km of an urban network. The hybrid setup we have implemented enables us to transfer frequency with a relative frequency stability as low as 3.9 × 10⁻¹⁷ at 1 s integration time and below 10⁻²¹ over 16 days. Combined with another hybrid links between LPL and NPL, this hybrid setup is already a part of an "optical fiber links chain" linking NPL, SYRTE and PTB [31]. With this work, we are showing that far from being competitive, the combination of two techniques can even lead to improve efficiency, redundancy and overall performances.

## Appendix

Here, we analyze in detail the phase noise in the hybrid setup. According to the noise sources, the total phase noise of the beat note can be given by the sum of four terms splitting the phase wandering into the laser term, the RF term, the interferometric term and the fiber term:

$$ \phi = \phi_{\text{laser}} + \phi_{\text{RF}} + \phi_{\text{inter}} + \phi_{\text{fiber}}. \qquad (4) $$

Since the same laser is used to inject the two fibers, and its frequency is almost constant, the first term can be neglected and we state $\phi_{\text{laser}} = 0$. By phase-locking the PLL reference synthesizers and the frequency counters to a common 10 MHz reference signal, the RF contribution to phase can be neglected, and we state $\phi_{\text{RF}} = 0$. The interferometric noise is due to the length mismatch in the local interferometers (a perfect match of the fiber length in the remote interferometer is accomplished by the use of a partial Faraday mirror at the remote setup) [18]. As the interferometers are isolated passively in a box with time constant which is much greater than the light propagation time in the fiber, we write the interferometric term as $\phi_{\text{inter}}(t)$. Thus for the four basic beat notes, the two dominant phase noises are the fiber link noise and interferometers' noise

$$ \left\{ \begin{aligned} \phi_{\text{ANC}}(t) &= \phi_{\text{inter-ANC}}(t) + \frac{1}{2}\phi_{\text{RT-F1}}(t) \\ \phi_{\text{RT}}(t) &= \phi_{\text{inter-RT}}(t) + \phi_{\text{RT-F2}}(t) \\ \phi_{\text{OWB}}(t) &= \phi_{\text{inter-OWB}}(t) + \phi_{\text{backward-F2}}(t) \\ \phi_{\text{OWF}}(t) &= \phi_{\text{inter-OWF}}(t) + \phi_{\text{forward-F2}}(t - \tau). \end{aligned} \right. \qquad (5) $$

where $\tau$ is the light propagation delay between the local site and the remote site, $\phi_{\text{inter-ANC}}$ is the interferometric noise in the ANC setup, $\phi_{\text{RT-F1}}$ is the roundtrip propagation noise on F1, and other notations have the similar definitions. Then the TWUs can be expressed as

$$ \left\{ \begin{aligned} \phi_{\text{TWU1}}(t) &= \phi_{\text{inter-TWU1}}(t) + \frac{1}{2}\phi_{\text{RT-F1}}(t) - \frac{1}{2}\phi_{\text{RT-F2}}(t) \\ \phi_{\text{TWU2}}(t) &= \phi_{\text{inter-TWU2}}(t) + \frac{1}{2}\phi_{\text{RT-F1}}(t) - \phi_{\text{backward-F2}}(t). \\ \phi_{\text{TWU3}}(t) &= \phi_{\text{inter-TWU3}}(t) + \frac{1}{2}\phi_{\text{RT-F1}}(t) - \phi_{\text{forward-F2}}(t - \tau) \end{aligned} \right. \qquad (6) $$

The TWBs can be expressed as
$$
\left\{
\begin{aligned}
\phi_{\text{TWB1}}(t) &= \phi_{\text{inter-TWB1}}(t) + \phi_{\text{backward-F2}}(t) - \frac{1}{2}\phi_{\text{RT-F2}}(t) \\
\phi_{\text{TWB2}}(t) &= \phi_{\text{inter-TWB2}}(t) + \phi_{\text{forward-F2}}(t - \tau) - \frac{1}{2}\phi_{\text{RT-F2}}(t). \\
\phi_{\text{TWB3}}(t) &= \phi_{\text{inter-TWB3}}(t) + \frac{1}{2}[\phi_{\text{backward-F2}}(t) - \phi_{\text{forward-F2}}(t - \tau)]
\end{aligned}
\right.
\qquad (7)
$$

It is demonstrated that the interferometric noises arise from the fiber sensitivity to external perturbations. Since acoustic and mechanic noises can be made negligible thanks to the well-designed isolation box, the so-called interferometric noise is mainly the result of long-term thermal effects. Inside the interferometers enclosure, the temperature is approximately homogeneous and we have

$$ \phi_{\text{inter}} = 2\pi\nu\gamma\Delta T\delta L. \qquad (8) $$

where $\gamma$ is a phase-temperature coefficient of silica fiber, which has a value of 37 fs/(K·m) for an optical carrier at 194.4 THz and at 298 K, and $\Delta T$ is the temperature variation in the local interferometers. We omit here the time dependency for sake of clarity. For OWB, OWF and RT, the temperature variation is the same, but the fiber length mismatch is different. The local fiber length mismatch consists of the uncommon paths from the ultra-stable laser to two Michelson-like interferometers and the two short arms of the two interferometers. Following the approach developed in [18], the RT term will contain an interferometric noise proportional to the length $2L_{BC}$; the OWB term contains an interferometric noise proportional to the optical length difference $2L_{BC} + L_{OA} - L_{OB}$; the OWF term contains an interferometric noise only proportional to the optical length difference $L_{OB} - L_{OA}$ (see Fig. 1 for the definition of point O, A, B and C). Thus the interferometric noises are given by

$$
\left\{
\begin{aligned}
\phi_{\text{inter-RT}} &= 2\pi\nu\gamma\Delta T \cdot 2L_{BC} \\
\phi_{\text{inter-OWB}} &= 2\pi\nu\gamma\Delta T \cdot (2L_{BC} + L_{OA} - L_{OB}). \\
\phi_{\text{inter-OWF}} &= 2\pi\nu\gamma\Delta T \cdot (L_{OB} - L_{OA})
\end{aligned}
\right.
\qquad (9)
$$

And the interferometric noise term for TWBs can be expressed as

$$
\left\{
\begin{aligned}
\phi_{\text{inter-TWB1}} &= 2\pi\nu\gamma\Delta T \cdot (L_{BC} + L_{OA} - L_{OB}) \\
\phi_{\text{inter-TWB2}} &= 2\pi\nu\gamma\Delta T \cdot (-L_{BC} - L_{OA} + L_{OB}). \\
\phi_{\text{inter-TWB3}} &= 2\pi\nu\gamma\Delta T \cdot (L_{BC} + L_{OA} - L_{OB})
\end{aligned}
\right.
\qquad (10)
$$

The phase evolutions for TWB1 and TWB2 are opposite, and those for TWB1 and TWB3 are the same, as shown in Fig. 4(a). We measured in [18] that the length mismatch $L_{BC} + L_{OA} - L_{OB} \approx 15$ cm. The phase evolution shown in Fig. 7(b) is compatible with a temperature variation of 1.5 K, a temperature variation of the local interferometers that we actually confirm experimentally during the acquisition period.

For the definition of TWNF, under the combination of OWB, OWF and RT, the interferometric noises compensate themselves, so

$$ \phi_{\text{inter-TWNF}} = 0. \qquad (11) $$

Showing that TWNF is interferometric-noise free.

## Funding

Agence Nationale de la Recherche (ANR) (ANR-10-LABX-48-01, ANR-11-EQPX-0039); H2020 (SIB02, 15SIB05); Korea Research Institute of Standards and Science (KRISS) (16011007); National Research Council of Science and Technology (NST) (CAP-15-08-KRISS).
$$ egin{cases} 
\phi_{\text{TWB1}}(t) = \phi_{\text{inter-TWB1}}(t) + \phi_{\text{backward-F2}}(t) - \frac{1}{2}\phi_{\text{RT-F2}}(t) \\ 
\phi_{\text{TWB2}}(t) = \phi_{\text{inter-TWB2}}(t) + \phi_{\text{forward-F2}}(t - \tau) - \frac{1}{2}\phi_{\text{RT-F2}}(t). \\ 
\phi_{\text{TWB3}}(t) = \phi_{\text{inter-TWB3}}(t) + \frac{1}{2}[\phi_{\text{backward-F2}}(t) - \phi_{\text{forward-F2}}(t - \tau)] 
\end{cases} 
\quad (7) $$

It is demonstrated that the interferometric noises arise from the fiber sensitivity to external perturbations. Since acoustic and mechanic noises can be made negligible thanks to the well-designed isolation box, the so-called interferometric noise is mainly the result of long-term thermal effects. Inside the interferometers enclosure, the temperature is approximately homogeneous and we have

$$ \phi_{\text{inter}} = 2\pi\nu\gamma\Delta T\delta L. \quad (8) $$

where $\gamma$ is a phase-temperature coefficient of silica fiber, which has a value of 37 fs/(K·m) for an optical carrier at 194.4 THz and at 298 K, and $\Delta T$ is the temperature variation in the local interferometers. We omit here the time dependency for sake of clarity. For OWB, OWF and RT, the temperature variation is the same, but the fiber length mismatch is different. The local fiber length mismatch consists of the uncommon paths from the ultra-stable laser to two Michelson-like interferometers and the two short arms of the two interferometers. Following the approach developed in [18], the RT term will contain an interferometric noise proportional to the length $2L_{BC}$; the OWB term contains an interferometric noise proportional to the optical length difference $2L_{BC} + L_{OA} - L_{OB}$; the OWF term contains an interferometric noise only proportional to the optical length difference $L_{OB} - L_{OA}$ (see Fig. 1 for the definition of point O, A, B and C). Thus the interferometric noises are given by

$$ egin{cases} 
\phi_{\text{inter-RT}} = 2\pi\nu\gamma\Delta T \cdot 2L_{BC} \\ 
\phi_{\text{inter-OWB}} = 2\pi\nu\gamma\Delta T \cdot (2L_{BC} + L_{OA} - L_{OB}). \\ 
\phi_{\text{inter-OWF}} = 2\pi\nu\gamma\Delta T \cdot (L_{OB} - L_{OA}) 
\end{cases} 
\quad (9) $$

And the interferometric noise term for TWBs can be expressed as

$$ egin{cases} 
\phi_{\text{inter-TWB1}} = 2\pi\nu\gamma\Delta T \cdot (L_{BC} + L_{OA} - L_{OB}) \\ 
\phi_{\text{inter-TWB2}} = 2\pi\nu\gamma\Delta T \cdot (-L_{BC} - L_{OA} + L_{OB}). \\ 
\phi_{\text{inter-TWB3}} = 2\pi\nu\gamma\Delta T \cdot (L_{BC} + L_{OA} - L_{OB}) 
\end{cases} 
\quad (10) $$

The phase evolutions for TWB1 and TWB2 are opposite, and those for TWB1 and TWB3 are the same, as shown in Fig. 4(a). We measured in [18] that the length mismatch $L_{BC} + L_{OA} - L_{OB} \approx 15$ cm. The phase evolution shown in Fig. 7(b) is compatible with a temperature variation of 1.5 K, a temperature variation of the local interferometers that we actually confirm experimentally during the acquisition period.

For the definition of TWNF, under the combination of OWB, OWF and RT, the interferometric noises compensate themselves, so

$$ \phi_{\text{inter-TWNF}} = 0. \quad (11) $$

Showing that TWNF is interferometric-noise free.

## Funding

Agence Nationale de la Recherche (ANR) (ANR-10-LABX-48-01, ANR-11-EQPX-0039); H2020 (SIB02, 15SIB05); Korea Research Institute of Standards and Science (KRISS) (16011007); National Research Council of Science and Technology (NST) (CAP-15-08-KRISS).
## Acknowledgments

The authors would like to thank Florian Frank for providing assistance and help in many experimental aspects. We acknowledge funding support from the Agence Nationale de la Recherche (Labex First-TF ANR-10-LABX-48-01, Equipex REFIMEVE+ANR-11-EQPX-0039), Conseil Régional Ile-de-France (DIM Nano'K), CNRS with Action Spécifique GRAM, the European Metrology Research Programme (EMRP) in the Joint Research Projects SIB02 (NEAT-FT) and the European Metrology Programme for Innovation and Research (EMPIR) in project 15SIB05 (OFTEN). The EMRP and EMPIR are jointly funded by the EMRP/EMPIR participating countries within EURAMET and the European Union.
## Acknowledgments

The authors would like to thank Florian Frank for providing assistance and help in many experimental aspects. We acknowledge funding support from the Agence Nationale de la Recherche (Labex First-TF ANR-10-LABX-48-01, Equipex REFIMEVE+ANR-11-EQPX-0039), Conseil Régional Ile-de-France (DIM Nano'K), CNRS with Action Spécifique GRAM, the European Metrology Research Programme (EMRP) in the Joint Research Projects SIB02 (NEAT-FT) and the European Metrology Programme for Innovation and Research (EMPIR) in project 15SIB05 (OFTEN). The EMRP and EMPIR are jointly funded by the EMRP/EMPIR participating countries within EURAMET and the European Union.