The dataset contains the replication data used for generating plots and tables in the associated article. It contains both the data itself, stored in netcdf4 and pickle format. The data is best analyzed using a Python based stack for data analysis. The following packages are recommended for its viewing and running the examples below:

• numpy
• matplotlib
• xarray
• xrft
• arc-alkali-rydberg-calculator
• netcdf4
• pint
• tabulate

Due to the complexity of the data, analyzing them is best done based on the custom scripts below, developed while conducting the associated research. The structured nature of the netcdf format, as well as high level and quality code presented in the examples, serve as a documentation of its own. The code replicates all of the plots and values presented in the article.

All photon counts were measured using a single photon avalanche photodiode, while the probe spectra were measured using a regular avalanche photodiode, where the probe laser was scanned in frequency.

Contents

1. examples.py contains a Python script to be run in interactive mode, for viewing all the plots that are part of the linked article
2. comb_scan.nc contains broadband scans of the frequency comb showing photon counts with respect to the decoupling laser detuning, for 4 different atomic transitions. Used in Fig. 3-4
3. comb_theory.pkl theoretical predictions of the aforementioned scan for the 126 GHz transition, obtained with the theoretical model from appendix A, used in Fig. 4
4. map.nc conversion gain map, obtained from the aforementioned theoretical model. Used in Fig. 8-9 and to derive the theory in Fig. 5
5. photon_counts.nc measured photon counts for maximum comb signal intensity, thermal background and dark counts. Each measurement is a 1s aggregation and there are 30 copies of each level measurement. Used in Fig. 2 and sensitivity derivations shown in Table 1
6. rabi_fits.pkl experimental data showing Autler-Townes splitting, registered in the probe spectrum with a photodiode. Additionally, contains the resulting plots of the fitted theoretical model and estimated Rabi frequencies. Used in Fig. 7
7. rabis.pkl fitted Rabi frequencies with respect to the half-waveplate angles used for attenuating the signal. Used in Fig. 2 and the calibration procedure
8. tds.nc single frequency comb pulse spectrum, measured with a standard Time Domain Spectroscopy setup. Used in Fig. 3
9. tra_da.nc measured photon counts of the converted TRA chip signal with respect to the attenuation induced by the half-waveplate angles and the selected chip power mode. Used in Fig. 2
10. bandwidth.nc measured photon counts of the converted TRA chip signal with respect to its detuning, driven by a phase locked loop. Used as the conversion bandwidth for frequency comb mode number derivations presented in Table 2
11. beat_histogram.nc measured photon counts with respect to time, used for observing beat notes and shown in Fig. 5

This research was funded in whole or in part by National Science Centre, Poland grants No. 2021/43/D/ST2/03114 and No. 2024/53/N/ST7/02730. The "Quantum Optical Technologies" (FENG.02.01-IP.05-0017/23) project is carried out within the Measure 2.1 International Research Agendas programme of the Foundation for Polish Science co-financed by the European Union under the European Funds for Smart Economy 2021-2027 (FENG).