Atomic-Clock-Based Synchronization

Updated 4 July 2026

Atomic-clock-based synchronization is a technique that aligns clocks by referencing atomic time scales, using GNSS, calibrated fiber links, and quantum protocols to overcome network delays and oscillator drift.
It employs detailed clock models, including linear and polynomial drift approximations and ensemble mean estimations with Kalman filtering, to accurately estimate and correct time offsets.
Quantum and entanglement-based methods further enhance precision by facilitating picosecond-level adjustments, crucial for secure communications and global time transfer.

Atomic-clock-based synchronization denotes the family of methods that align clocks, timestamps, oscillator frequencies, or symbol phases by referencing physically realized atomic time scales and by steering local clocks or ensemble means under measurement, propagation, and calibration constraints. In the recent literature, the topic spans atomic and quartz clock models driven by GNSS-derived references and explicit network propagation paths, global UTC/TAI comparison by satellite and optical fiber, real-time correction of free-running rubidium and cesium clocks against GPS Time, ensemble time-scale generation by observable canonical decomposition, and quantum protocols that use entangled photons for picosecond-class syntonization or sub-nanosecond global time transfer (Dai et al., 2024, Piester et al., 2011, Ishizaki et al., 22 Apr 2025, Pelet et al., 28 Jan 2025).

1. Time scales, reference systems, and operational notions

Atomic-clock-based synchronization begins from the distinction between atomic time as a physical realization, operational time scales such as TAI and UTC, and practical dissemination references such as GPS Time and other GNSS time scales. One line of work traces atomic time to early cesium standards and highlights NIST’s NIST-F2 cesium primary standard, which “will not gain or lose a second in at least 300 million years,” while the metrological infrastructure for International Atomic Time and Coordinated Universal Time is organized through BIPM, UTC(k) realizations, and interlaboratory comparisons using GPS and TWSTFT (Dai et al., 2024, Piester et al., 2011). GPS Time is a continuous timescale started at 1980-01-06 and does not include leap seconds, whereas UTC includes leap seconds; as of mid-2026, $\Delta(\mathrm{GPS}-\mathrm{UTC}) = 18\ \mathrm{s}$ (Dalmazzone et al., 2024).

A second distinction is between synchronization and syntonization. In the entanglement-based QKD literature, synchronization is defined as sharing both frequency and phase/date so two clocks tick at the same rate and read the same time, whereas syntonization is sharing only frequency so two clocks tick at the same rate but may have a fixed offset. For many elementary entanglement-based QKD links, syntonization suffices; full synchronization is needed for advanced multi-user quantum networking (Pelet et al., 28 Jan 2025).

A third distinction concerns what is being synchronized. Myers and Madjid formulate “logical synchronization” not as the establishment of a globally consistent time coordinate, but as phase-meshed symbol transfer between open machines stepped by local clocks. Each reading of a clock of an open machine $A$ has the form $m.\phi_m$ with $m$ the cycle count and $\phi_m$ the phase within the cycle, with $-1/2 < \phi_m \le 1/2$ . Successful symbol acceptance requires that, for some positive $\eta$ , any arrival phase $\phi_n$ satisfy

$|\phi_n| < (1-\eta)/2.$

This phase-window condition defines synchronization as local write-phase compatibility rather than global simultaneity (Myers et al., 2016).

These operational notions are not interchangeable. One literature centers epoch alignment to UTC or GPS Time; another centers rate alignment; another centers symbol acceptance within write/read windows. Atomic-clock-based synchronization therefore includes both metrological time transfer and local control of phase-sensitive computation and communication.

2. Clock models, drift processes, and ensemble means

At the device level, atomic-clock-based synchronization is built on explicit clock models. P-TimeSync defines wall-clock time $t$ , a software clock $A$ 0, and instantaneous offset

$A$ 1

It supports a linear model,

$A$ 2

and a second-order polynomial model,

$A$ 3

where $A$ 4 is time offset, $A$ 5 is frequency offset, $A$ 6 captures frequency drift, and $A$ 7 is random noise (Dai et al., 2024). The same work shows that fixed second-order drift can be physically misleading: with $A$ 8 and $A$ 9, the model reaches a local maximum at $m.\phi_m$ 0, which it describes as not realistic for physical clocks that do not self-correct after a peak. This motivates user-defined stochastic noise processes instead of fixed polynomial drift (Dai et al., 2024).

The same framework uses the conceptual clock and oscillator relations

$m.\phi_m$ 1

with $m.\phi_m$ 2 the initial offset, $m.\phi_m$ 3 the fractional frequency offset, and $m.\phi_m$ 4 and $m.\phi_m$ 5 noise processes. Atomic clocks can be configured with smaller $m.\phi_m$ 6 and lower noise than quartz oscillators, while environmental sensitivities can be injected through $m.\phi_m$ 7 (Dai et al., 2024). The practical motivation is explicit: miniaturized rubidium clocks and quartz oscillators are widely used, but MRC stability can degrade by at least an order of magnitude under high dynamics of about $m.\phi_m$ 8, and quartz oscillators are sensitive to temperature, humidity, pressure, acceleration, vibration, electric or magnetic fields, load, and radiation. Representative MRC performance cited in that work is approximately $m.\phi_m$ 9 for the AccuBeat AR133 series and about $m$ 0 for the Microsemi SA.3Xm (Dai et al., 2024).

At the ensemble level, synchronization becomes a state-estimation and control problem. “Explicit Ensemble Mean Clock Synchronization for Optimal Atomic Time Scale Generation” models a single atomic clock by phase and fractional frequency states, with white and random-walk frequency noise, and gives the free-running Allan variance

$m$ 1

The framework then decomposes an $m$ 2-clock ensemble into an observable synchronization-deviation component and an unobservable synchronization-destination component, identified as the ensemble mean. Within that decomposition, standard Kalman filtering is proved to be a special case corresponding to long-term Allan-variance optimization (Ishizaki et al., 22 Apr 2025).

The mixed-ensemble extension to cesium-type and hydrogen-maser-type clocks replaces Allan-variance weighting by Hadamard-variance weighting. For the unobservable mean, it gives

$m$ 3

with optimal weights

$m$ 4

Its short-term limit emphasizes the lowest white-frequency-noise clocks, whereas the long-term limit assigns zero weight to masers and weights only the cesium clocks (Dey et al., 23 May 2026). In this body of work, synchronization is not only pairwise clock matching; it is also explicit construction of a generated time scale.

3. Propagation asymmetry, two-way estimation, and correction loops

Clock models alone do not determine synchronization performance. P-TimeSync makes network propagation an explicit part of the timing model by decomposing total delay into router delay, bandwidth-dependent serialization delay, and distance-dependent propagation delay, and by dynamically selecting routes with Dijkstra’s algorithm. It states that “the path for sending data and the path for receiving data may differ” as router states change, and defines the total end-to-end delay as the sum of the active-path components. The simulator combines link bandwidths from Kbps to Gbps, distances from meters to thousands of kilometers, media including Wi‑Fi, Ethernet or fiber, and satellite or LEO paths, and per-hop delays and loss, all at nanosecond-level resolution (Dai et al., 2024).

Within that network model, P-TimeSync implements Cristian’s and Berkeley algorithms. Cristian’s algorithm estimates client time from a server reference using measured RTT, conceptually $m$ 5, while Berkeley performs distributed averaging among peers and allows weights to reflect clock quality, such as atomic versus quartz. The same work places these software algorithms in the context of NTP and PTP, noting that NTP is software-centric and millisecond-level, whereas PTP is hardware-centric and nanosecond-level, and that hardware timestamping is often required to achieve true nanosecond performance in practice (Dai et al., 2024).

Two-way reciprocity remains central in metrological links. For stations $m$ 6 and $m$ 7, the PTB treatment of TWSTFT and two-way fiber transfer writes

$m$ 8

$m$ 9

and, under reciprocity $\phi_m$ 0,

$\phi_m$ 1

This formulation isolates propagation cancellation from internal transmit–receive delay calibration (Piester et al., 2011).

A different correction-loop architecture appears in experimental particle physics. For Hyper-Kamiokande, the required target is synchronization with UTC and between sites with a precision better than $\phi_m$ 2. The proposed method compares a free-running rubidium clock to GPS Time using a GNSS receiver, fits the measured time difference $\phi_m$ 3 over a window with a polynomial

$\phi_m$ 4

and corrects timestamps by

$\phi_m$ 5

With an integration time window around $\phi_m$ 6 seconds, the corrected time difference stays within a $\phi_m$ 7 range in both offline and online modes (Dalmazzone et al., 2024).

A closely related real-time implementation was later applied to a low-cost rubidium clock and a magnetic cesium clock. It uses GNSS observations every $\phi_m$ 8 minutes, rolling least-squares fits, and digital correction of timestamps rather than oscillator disciplining. The reported residual difference to UTC(OP) remains within $\phi_m$ 9 with no apparent residual drift, with standard deviations of $-1/2 < \phi_m \le 1/2$ 0 for rubidium and $-1/2 < \phi_m \le 1/2$ 1 for cesium (Dalmazzone et al., 28 Oct 2025). In these systems, atomic-clock-based synchronization is realized as a software correction layer over a free-running atomic time base rather than as continuous hardware steering.

4. Metrological infrastructures: GNSS, TWSTFT, and calibrated fiber

At global scale, atomic-clock-based synchronization is embedded in the infrastructure of TAI, UTC, UTC(k), and interlaboratory links. PTB’s survey describes GPS time transfer using the C/A code, the P code, and operational modes including “Single Channel,” “Multi Channel,” “P3,” and “Precise Point Positioning (PPP).” It reports approximate statistical uncertainties $-1/2 < \phi_m \le 1/2$ 2 of $-1/2 < \phi_m \le 1/2$ 3 for Single Channel, $-1/2 < \phi_m \le 1/2$ 4 for Multi Channel, $-1/2 < \phi_m \le 1/2$ 5 for P3, and $-1/2 < \phi_m \le 1/2$ 6 for PPP (Piester et al., 2011).

The same survey reports that TWSTFT operates in the Ku-band and benefits from reciprocity, so ionosphere and troposphere effects cancel to first order. It gives time noise at $-1/2 < \phi_m \le 1/2$ 7 averaging of about $-1/2 < \phi_m \le 1/2$ 8, a fractional frequency instability approximated by

$-1/2 < \phi_m \le 1/2$ 9

and time-comparison uncertainties around $\eta$ 0 once internal delays are calibrated. Calibration campaigns are reported down to $\eta$ 1 in the U.S. and in the range about $\eta$ 2 to about $\eta$ 3 across Europe (Piester et al., 2011).

Optical fiber is used where local distribution must not dominate the error budget. PTB’s calibrated two-way fiber system applies the same two-way equations and a common-clock calibration constant

$\eta$ 4

leading to the operational relation

$\eta$ 5

Over a $\eta$ 6 campus link it demonstrated uncertainty below $\eta$ 7, with mean CCD $\eta$ 8 and standard deviation of single-second measurements $\eta$ 9. For a $\phi_n$ 0 link with $\phi_n$ 1 bidirectional erbium-doped fiber amplifiers, the projected uncertainty is about $\phi_n$ 2 for total-link calibration and about $\phi_n$ 3 when calibrating each amplifier individually (Piester et al., 2011).

This metrological layer directly conditions application-specific systems. The Hyper-Kamiokande timing chain, for example, uses a White Rabbit path to transport signals and validates against UTC(OP), while GNSS measurements are used for long-term anchoring and polynomial correction of a free-running rubidium clock (Dalmazzone et al., 2024). The broader implication is that atomic-clock-based synchronization depends as much on calibrated dissemination and validation links as on the intrinsic stability of the oscillator.

5. Quantum and entanglement-based synchronization

Quantum versions of atomic-clock-based synchronization use photon correlations as the measurement primitive while retaining atomic clocks as local references or holdover oscillators. In a deployed entanglement-based QKD link over $\phi_n$ 4 of optical fibers across the Métropole Côte d’Azur, paired-photon time correlations are used to syntonize two rubidium clocks. The system maintains a time offset under $\phi_n$ 5 at all times over multiple $\phi_n$ 6-hour runs, with one representative run showing a mean offset of $\phi_n$ 7 and a variance of $\phi_n$ 8. The coincidence peak has FWHM about $\phi_n$ 9, the coincidence window is set to $|\phi_n| < (1-\eta)/2.$ 0, and the BBM92 link runs continuously with an average secret key rate of about $|\phi_n| < (1-\eta)/2.$ 1 (Pelet et al., 28 Jan 2025).

Satellite proposals generalize the same idea to global distribution. One study of nanosatellite QCS networks uses SPDC sources, APDs, and CSAC-class clocks and concludes that establishing a global network of ground clocks synchronized to sub-nanosecond level, up to a few picoseconds, would be feasible. Its two-way estimator is path-independent:

$|\phi_n| < (1-\eta)/2.$ 2

A related “master clock in the sky” architecture studies a constellation of $|\phi_n| < (1-\eta)/2.$ 7 satellites in $|\phi_n| < (1-\eta)/2.$ 8 polar orbits at $|\phi_n| < (1-\eta)/2.$ 9 altitude. There, satellites reinforce each other’s sync capabilities by intra-orbit and inter-orbit synchronization, and the reported result is continuous global synchronization at sub-nanosecond precision. The required holdover for continuous inter-orbit synchronization is about $t$ 0 minutes at the $t$ 1 level, and $t$ 2 minutes yields continuous $t$ 3 shadows per orbit and robust city-to-city coverage (Ducoing et al., 2023).

A different quantum line considers direct synchronization of optical atomic clocks in LEO using two-mode squeezed vacuum. For a measurement window with $t$ 4 photons, the paper contrasts SQL scaling $t$ 5 and Heisenberg scaling $t$ 6, and reports quantum advantage over the SQL for symmetric channels with $t$ 7 even at $t$ 8 squeezing. For asymmetric channels, with $t$ 9, advantage persists when each path has $A$ 00. An explicit numerical example at $A$ 01, $A$ 02, and $A$ 03 yields $A$ 04, compared to an SQL value of about $A$ 05, described as about $A$ 06 improvement (Gosalia et al., 2023).

System	Architecture	Reported result
Entanglement-based QKD syntonization	48 km field QKD link with two rubidium clocks	time offset under 12 ps at all times
Satellite-based QCS network	nanosatellites with SPDC, APDs, and CSACs	sub-nanosecond level, up to a few picoseconds
Quantum-assisted master clock	50 satellites, 5 orbits, 500 km altitude	continuous global synchronization at sub-nanosecond precision

Across these quantum protocols, the atomic clock remains central: it provides the local time base, the holdover resource, or the optical reference being synchronized. The quantum resource changes the time-transfer observable and, in some cases, the achievable precision regime.

6. Constraints, applications, and open directions

Atomic-clock-based synchronization is limited not only by oscillator noise and network asymmetry, but also by physical and architectural constraints. In logical synchronization, curvature can force non-null phases on some channels. For a five-machine cluster in Schwarzschild geometry, Myers and Madjid derive a curvature-dependent phase and, from the phase-window condition, a lower bound on clock period

$A$ 07

For $A$ 08 and descent to $A$ 09 around Earth, the phase-window constraint implies

$A$ 10

so if an alphabet conveys $A$ 11 bits per character, the maximum bit rate satisfies $A$ 12 (Myers et al., 2016). This is not a clock-instability bound; it is a communication constraint induced by curvature through phase acceptance.

Application requirements differ sharply. Hyper-Kamiokande requires precision better than $A$ 13 between sites and to UTC; the rubidium-plus-GNSS correction method meets that requirement with margin by keeping residuals within $A$ 14 (Dalmazzone et al., 2024). Entanglement-based QKD requires timing that keeps the coincidence peak inside a narrow window; the deployed syntonization system maintains under- $A$ 15 offset to preserve operation of a link with about $A$ 16 secret key rate (Pelet et al., 28 Jan 2025). Global metrology seeks $A$ 17 in the ACES context, while quantum satellite constellations target sub-nanosecond continuous worldwide coverage (Piester et al., 2011, Ducoing et al., 2023).

Simulation and control frameworks suggest concrete engineering practices. P-TimeSync emphasizes higher-stability references for grandmasters, explicit modeling of per-hop delays, minimization of path asymmetry, and the use of hardware timestamping and two-way transfer for nanosecond requirements. It also notes that dynamic shortest-path routing can itself inject jitter and that future work includes ML-based delay prediction and security hardening (Dai et al., 2024). The real-time GNSS correction literature similarly favors tunable digital correction of timestamps over direct disciplining when auditability, robustness to intermittent GNSS reception, or preservation of short-term oscillator stability is required (Dalmazzone et al., 28 Oct 2025).

Several limitations remain explicit in the cited literature. Logical synchronization does not provide fully specified steering update laws (Myers et al., 2016). The QCS satellite simulations omit some environmental and relativistic effects from their numerics, even when they discuss the relevant correction terms (Haldar et al., 2022, Ducoing et al., 2023). Conventional Kalman filtering in ensemble time-scale generation exhibits divergence in the unobservable covariance subspace, which motivated observable-canonical-decomposition filters and explicit ensemble mean synchronization (Ishizaki et al., 22 Apr 2025, Dey et al., 23 May 2026).

Taken together, these lines of work show that atomic-clock-based synchronization is not a single protocol class. It is a layered field that includes time-scale realization, oscillator steering, delay calibration, network-aware correction, ensemble control, phase-window communication, and quantum time transfer. The unifying theme is the same: atomic clocks provide the physically realized temporal reference, but synchronization performance is determined jointly by the clock, the transfer channel, the estimator, and the control architecture.