The role of master clock stability in quantum information processing

Experimentalists seeking to improve the coherent lifetimes of quantum bits have generally focused on mitigating decoherence mechanisms through, for example, improvements to qubit designs and materials, and system isolation from environmental perturbations. In the case of the phase degree of freedom in a quantum superposition, however, the coherence that must be preserved is not solely internal to the qubit, but rather necessarily includes that of the qubit relative to the ‘master clock’ (e.g., a local oscillator) that governs its control system. In this manuscript, we articulate the impact of instabilities in the master clock on qubit phase coherence and provide tools to calculate the contributions to qubit error arising from these processes. We first connect standard oscillator phase-noise metrics to their corresponding qubit dephasing spectral densities. We then use representative lab-grade and performance-grade oscillator specifications to calculate operational fidelity bounds on trapped-ion and superconducting qubits with relatively slow and fast operation times. We discuss the relevance of these bounds for quantum error correction in contemporary experiments and future large-scale quantum information systems, and consider potential means to improve master clock stability.


INTRODUCTION
A fundamental challenge to the broad application of quantum information science is the management of error in fragile quantum hardware. 1,2 The need for higher fidelity performance motivates research at all architectural levels, 3 from theoretical studies of fault tolerance and analyses of quantum error correction (QEC) implementations down to experimental improvements in the operational fidelity of elemental devices. A familiar aspect of this challenge is decoherence, a process by which even idle qubits undergoing free evolution (i.e., the identity operator) will gradually lose their stored quantum information, rendering them useless in subsequent computation.
A prevalent component of decoherence is dephasing, a randomisation of the relative phase between the basis states that form a coherent superposition state. Dephasing in qubit systems is commonly attributed to environmental fluctuations of a qubit bias or control parameter, e.g., an external magnetic field that modulates the qubit-state energy splitting and hence its dynamic phase evolution. In turn, experimentalists have primarily focused their efforts on reducing both the level of environmental noise and the qubit sensitivity to that which remains (c.f., refs 4-6). As a result, many qubit technologies 7 have witnessed remarkable performance improvements, with dephasing times approaching milliseconds or even seconds (depending on the qubit modality), 6,[8][9][10][11][12][13] and operational fidelities reaching 99.999%. [14][15][16] Achieving targets such as these-once thought all but impossible-has led many to recognise that scalable, error-corrected quantum computation may become a reality. In essence, by improving qubit coherence, error rates may be reduced below fault-tolerant thresholds, and, in principle, QEC may be employed to suppress errors at arbitrary system scales.
This prescription, however, belies issues that only arise as qubit coherence and operational fidelity improve. A fundamental example, one widely studied in the context of precision metrology, is that the control and measurement of highly coherent systems must be compared with a suitably coherent reference. 17 Qubit coherence, likewise, is inferred relative to a reference, generally in the form of a local oscillator (LO) used to control and interrogate the system, for instance by inducing Ramsey or Rabi oscillations. Here the evolution of the qubit's phase degree of freedom in a quantum superposition is effectively being compared against the accumulated ticks of a clock defined by the LO. 17 Although the importance of phase coherence is known to any experimentalist who has failed to lock an oscillator to a stable reference, in the past, qubit coherence times and operational fidelities have generally been limited by environmental noise. At today's performance levels, however, noise in the master clock is beginning to emerge as a contributor to the overall error rate. And, as qubits continue to improve, mitigating master clock instabilities will become material to high-fidelity operation.
This manuscript addresses the topic of master clock instability and its impact on qubit coherence and operational fidelities for quantum information applications. We begin by accounting for the master clock in a Hamiltonian treatment of qubit dephasing. From this starting point, we engage in an analysis of the relevance of master clock instabilities to qubit coherence and dephasing-induced error rates. We identify master clock phase fluctuations as an emerging source of error in today's qubit systems and one that will certainly become more prominent as qubit error rates continue to improve. Our presentation serves to unify concepts familiar to quantum information, 1 quantum control, 9,18,19 engineering 20,21 and precision frequency metrology 22 in order to allow the translation of LO phase noise specifications-as presented in experimental data sheets-into gate error probabilities for a variety of canonical single-qubit operations. We consider gates applied to superconducting and trapped-ion qubits, prominent qubit modalities that represent systems with relatively fast (10 ns) and slow (100 μs) gate times yielding sensitivity to different frequency regions of the oscillator noise spectrum. Calculations presented in this analysis highlight the fact that qubit dephasing errors arising from the phase fluctuations in commercial precision LOs, although not explicitly limiting experiments today, will likely become a significant consideration in the context of large-scale quantum information. They also reveal that many lab-grade oscillators are beginning to limit achievable fidelities in contemporary systems. We address the relevance of different frequency regimes of LO phase noise, and highlight 'far-from-carrier' phase noise as contributing to important upper bounds on qubit operational fidelities. The materials we present both demonstrate a path to mesoscale quantum information systems using existing master clocks and provide new motivation for investment in LO hardware and precision frequency metrology research in order to underpin large-scale quantum information processing. We augment this review with detailed Supplementary Information, aggregating a comprehensive theoretical foundation to understand clockinduced errors in quantum systems.

QUBIT DEPHASING INDUCED BY THE MASTER CLOCK
In a semiclassical picture, qubit phase coherence corresponds to maintaining the dynamical phase of the qubit relative to a reference (master) clock for the system. Implicit in the Bloch sphere representation of a qubit state is, in general, a transformation to a frame co-rotating with the nominal qubit Larmor frequency, o 0 . In a dephasing process, the relative phase difference between the qubit state and its reference frame evolves stochastically in time, introducing a degree of randomness to the qubit state. Typical formulations assign this stocastic evolution to instability in qubit frequency, induced, for instance, by fluctuations in external fields. Over time, these fluctuations cause the qubit to randomly advance or retreat relative to the co-rotating frame. However, what about the stability of the rotating frame being used as a phase reference? The simple observation that the qubit phase is defined relative to the rotating frame indicates that an experimentalist must consider not only the stability of the qubit but also the stability of the rotating frame itself, realised through the use of a local oscillator in the experiment. Accordingly, the master clock has a role that is more fundamental than the synchronisation of scheduled operations; the master clock determines, in part, the coherence of the underlying qubits (see Fig. 1). The impact of generic qubit dephasing may be formally captured through a Hamiltonian formulation, written as the sum of an ideal control component and a randomly fluctuating noise component in the three-dimensional Pauli basis.
where we have restricted attention to longitudinal dephasing noise by setting δh x,y (t) = 0. Terms from environmental dephasing δh env ð Þ z t ð Þ and LO-induced dephasing, δh LO ð Þ z t ð Þ, appear on an equal footing in this formulation. This Hamiltonian considers an ensemble of identically prepared noninteracting qubits subjected to a noisy environment and an external controller. As we assume an interaction picture co-rotating with a resonant driving field, the freely evolving qubit Hamiltonian H 0 1 2 o 0σz does not appear explicitly (see Equation (6) in Methods; for full derivation see Supplementary Information). Once a relative phase relationship between the qubit and LO is established via an initial interaction (control pulse), the presence of these two noise terms accounts for qubit dephasing at all times t.
Within this framework, we can build a connection between the LO-induced dephasing term and the noise properties of local oscillators. A local oscillator with amplitude A(t) and carrier frequency o LO possesses a phase that is easily partitioned into a 'control' component (desired, deterministic phase evolution), ϕ C (t), and a 'noise' component (unwanted, stochastic phase evolution), ϕ N (t), in the expression for the oscillator signal: ). Setting the LO frequency to be on resonance with the qubit, o LO = o 0 , the stability of the rotating frame is determined by the noisy phase evolution of the LO, ϕ N (t), as it produces a time-dependent frequency detuning of the LO relative to the qubit through its time derivative, δo LO t ð Þ _ ϕ N t ð Þ. With this qualitiative understanding in hand, we are able to associate the Hamiltonian term δh LO ð Þ, explicitly linking LO phase fluctuations to a dephasing Hamiltonian with an impact that is indistiguishable from other environmental dephasing sources.
In Figure 1, the equivalence of the two dephasing terms described above is represented on the Bloch sphere in the frame co-rotating with the resonant driving field. The term δh env ð Þ z t ð Þ induces rotations of the qubit Bloch vector about theẑ axis of the Bloch sphere, relative to a fixed coordinate frame set by an ideal LO. Conversely, LO instability captured by δh LO ð Þ z t ð Þ induces rotations of the coordinate frame relative to a fixed Bloch vector representing an ideal qubit. From the perspective of the qubit, both terms are sources of dephasing. Hereafter, without loss of generality, we will ignore environmental dephasing by setting δh env ð Þ z t ð Þ = 0 in order to restrict attention to LO-induced dephasing.

CALCULATING DEPHASING DUE TO PHASE NOISE IN THE LOCAL OSCILLATOR
We employ the filter transfer function formalism to incorporate noise power spectral densities for the Hamiltonian noise terms into calculations of operational fidelity and qubit coherence. 21 This approach treats the problem of calculating the impact of noise on a quantum system in terms of the overlap of the noise and an effective frequency-domain filter describing the action of the control. 23,24 It adopts concepts that are well known in the engineering literature 25 and has recently been developed and applied to the control of qubits. [26][27][28][29] In essence, the filter transfer function shapes the spectrum of the underlying noise as it couples to the qubit, passing certain bands and rejecting others. We see this quantitatively for a system Hamiltonian in the form of Equation (1) by writing the fidelity of an arbitrary qubit operation in terms of the overlap integral of S 1 ð Þ z o ð Þ, the unilateral (one-sided) power spectral density (PSD) (see Supplementary Information for discussion of the use of unilateral and bilateral power spectral densities) of the dephasing field δh z (t), and the filter transfer function G z,l (o), which captures the action of the control Hamiltonian h(t)σ. 21,30 We express the average fidelity The fact that the error integral, χ τ ð Þ, is expressed as a product of the noise power spectral density and filter transfer functions demonstrates how the frequency-domain shaping of the noise determines the ultimate contribution of the noise to operational infidelity. The transfer functions for the control, G z,l (o), may be calculated analytically and have contributions along all Cartesian directions (indexed by l), as dephasing noise present during a non-commuting control operation (e.g., pσ x ) induces both dephasing and amplitude damping. 21,30 Alternatively, the G z,l (o) corresponding to a Ramsey pulse sequence would determine dephasing during free evolution.
Our objective is to link noisy fluctuations in the LO phase to the quantity S 1 ð Þ z o ð Þ relevant to the Hamiltonian noise terms. As we are restricting attention to the case δh z (t) = δh LO ð Þ z t ð Þ (omitting environmental noise for convenience), the PSD S 1 ð Þ z o ð Þ will solely describe fluctuations in the LO phase. For this, we reference a large base of research from the field of precision metrology, where the characterisation of time and frequency signals is a key objective. 17 Conveniently, in this discipline the stability of a signal at a notionally fixed frequency is quantified by characterising temporal fluctuations in the phase of that signal, represented in the Fourier domain. Similar analyses are found in the context of communications engineering and stochastic signal processing.
The quantity of metrological significance 22 , the unilateral power spectral density of the LO phase fluctuations over a measurement time T. In this expression, is the complex amplitude of the harmonic Fourier component at frequency o for the time-gated signal ϕ N (t), defined for times |t|o T/2 and zero otherwise. The quantity |Φ T (o)| 2 is the energy density of this harmonic with units of energy per Hertz; angle brackets indicate an expectation value, and stationarity has been assumed to take the limit T → ∞. The phase fluctuations generate sidebands on the carrier-the unadulterated LO signal-at frequencies o LO ± o, effectively broadening the observed linewidth of the LO in frequency space (Figure 1). The Fourier harmonic o of the phase instability is therefore measured as an offset from the carrier, and the power of these sidebands is measured over a 1 Hz bandwidth.
The time derivative of the phase fluctuations is equivalent to a time-dependent detuning of the LO from resonance, as described above, and we may therefore relate their noise spectra to the dephasing power spectrum used in filter-function calculations as The factor 1/4 arises from omitting the counter-rotating term in the rotating wave approximation.
ð Þ is used in the metrological community, 22 most LO manufacturers use a metric for the single-sideband phase noise,L o ð Þ ¼ 10 log 10

expressed in logarithmic units
of dBc/Hz. The ultimate relationship therefore allows tabulated single-sideband phase-noise specifications to be converted to unilateral dephasing power spectral densities as This correspondence enables a connection between the metrics provided by oscillator manufacturers and those used for both precision oscillator characterisation and for qubit dephasing calcuations, 22 and it is represented graphically in the difference between curves in Figure 2a,b. A full accounting of all relevant factors and the formal derivation of these quantities is presented in the Supplementary Information in a single consistent notation. We have numerically confirmed that, for evolution times τ o100 ms, the presence of a sharp low-frequency cutoff below o/2πo1 Hz contributes a negligible correction to our calculated results based on extrapolatingL o ð Þ as described above.

CLOCK ERRORS INDUCED BY STATE-OF-THE-ART LOCAL OSCILLATORS
With Equation (5) in hand, it is possible to use LO phase-noise specifications found on instrumentation data sheets for calculations of qubit coherence and operational error rates. We calculate the operational infidelity 1 -F av τ ð Þ using published L o ð Þ for representative synthesisers and analytic filter transfer functions as inputs for Equation (4). We consider two standard operations: the identity operation,Î, and theX gate, defined as a driven rotation through angle π about thex axis of the Bloch sphere. These choices are representative of the single-qubit operations comprising a universal gate set for gate-based quantum computation. The former allows insight into the coherent lifetime of the qubit under free evolution (T 2 ), whereas the latter informs how gate fidelity is reduced during driven evolution by the presence of LO-induced dephasing noise. For each of these 'primitive' operations, the corresponding filter transfer function G z,l (o) is calculated following the techniques presented in refs 19,21,29-33. In addition to the primitive forms of these operations, we also employ specific dynamic error suppression (DES) strategies 24,[34][35][36][37][38][39][40][41][42][43][44][45][46][47][48] designed to reduce errors due to dephasing noise. Such controls generally constitute time-dependent modulation of the system dynamics with the aim of coherently averaging out slow fluctuations. Although these protocols have typically been associated with the mitigation of environmental decoherence, as the reader might expect based on the discussion above in 'Qubit dephasing induced by the master clock', a growing body of literature has shown that errors induced by imperfect controlincluding LO phase noise-can also be mitigated using these same techniques. 9,26,29 Here for the identity operation (free evolution), we use the simple example of a spin echo, whereas for the drivenX gate, we employ a dynamically corrected version of this gate, 44-46 also described as a Walsh amplitude modulated filter. 29,49 A description and analytic expressions for the filter functions used here may be found in the Supplementary Information.
For this analysis, we have selected two distinct grades of LO to demonstrate the significance of LO choice in quantum control experiments. Figure 2a shows the single-sideband phase noise at 10 GHz for both a 'lab-grade' synthesiser (Vaunix LMS-123) and a 'precision' synthesiser (Keysight/Agilent 8267D OPT-UNY). Across much of the band, as expected, the phase noise for the precision synthesiser is lower-by~40 dB-than the lab-grade unit. For the LOs studied here,L o ð Þ declines rapidly with offset frequency from the carrier, exhibiting various power-law-dependences over the entire frequency range.
Generically, these dependences are captured analytically by the form L o ð Þpo p (i.e., S z (o)∝o p+2 , see Eqs. 80-82 in Supplementary  Information). Moving away from the carrier, dominant processes common in LOs include random-walk frequency noise (p = − 4), flicker frequency noise (p = − 3), white frequency noise (p = − 2), flicker phase noise (p = − 1) and white phase noise (p = 0). Note that, as expected, the dephasing noise mechanisms carry an exponent p+2 with respect to the corresponding phase noise mechanisms with exponent p (Figure 2b). Details of the origins of the underlying processes may be found in ref. 20.
Calculations of infidelity for the four operations outlined above are presented in Figure 2c,d and summarised in Table 1. For the identity operator, we find an approximate improvement of 10 4 in residual error rate due to the use of a high-precision frequency source, and the infidelity remains below 10 − 6 out to 100 ms evolution time. In contrast, the lab-grade synthesiser induces an error exceeding~0.1% beyond a few microseconds of evolution time. The plateau-like behaviour in infidelity with increasing evolution time is due to the interplay of the dephasing power spectrum and filter-function, 50 and is similar to phase-error saturation phenomenology observed in precision oscillator characterisation. 20 For both free evolution and driven operations, the deleterious impact of using a LO increases as the duration of the operation grows, with relevance for all major qubit modalities, including trapped-ion, 51 superconducting 52 and semiconductor 53 qubits. For instance, for driven operations as long as 100 μs, particularly relevant for atomic qubits, the precision synthesiser only induces an error~10 − 7 , whereas the lab-grade synthesiser induces an error more than 20,000 × larger, reaching~0.2%. The primary reason for this behaviour is the superior far-from-carrier phasenoise performance in the precision synthesiser, especially over the offset frequency range 1-1,000 kHz. We note that the selected extrapolation of the high-frequency phase noise beyond o/2π = 10 MHz (see Figure 2) is particularly favourable to the lab-grade LO, meaning that calculated error rates at the shortest times-τ o100 ns-may underestimate the actual infidelities.
An interesting observation is that in the presence of these noise power spectral densities and for the realistic evolution times selected for our calculations, DES protocols have a minimal impact on gate performance. In Figure 2c,d, over the entire range of times-10 ns…100 ms-the use of these protocols offers only a small, sporadic improvement and, over several time spans, can even degrade the operational fidelity (note: these protocols may still give a substantial net improvement in the presence of more dominant environmental noise). This performance is explained by considering that, whereasL o ð Þ declines rapidly with offset frequency, the o 2 factor that transforms this spectrum to S z (o) reveals a high-frequency dominance of the resulting dephasing noise, in particular, for the precision LO, which exhibits an approximately Ohmic spectrum (S z (o) ∝ o; see Figure 2b). 54,55 DES is known to perform poorly for spectra with strong high-frequency content, where the noise evolves rapidly compared to the control, and the physics of coherent averaging fails-a violation of the socalled decoupling limit.
These results make clear that reducing far-from-carrier phase noise has the potential to provide augmented fidelities for quantum operations. However, it remains a question how well one might do by improving the quality of the LO. We consider achievable gains in operational fidelity by calculating a lower-bound to the error imparted by the LO's thermal noise floor; even if LO hardware were improved, we could do no better than saturating the thermal noise floor across the control bandwidth. This is generically quoted as − 174 dBm/Hz for a Main results calculated assuming an integration bandwidth of 10 MHz. The temperature-dependent noise-floor for different cutoff frequencies is also shown-these data are largely independent of operation time and type. Noise floor calculations strictly valid only for cutoff frequencies high compared with the inverse operation time.
Master clock stability in quantum information H Ball et al matched load at 290 K. This is an absolute noise power, meaning that for LOs with power 0 dBm the single-sideband phase noise floor will take a valueL min ¼ -174 dBc/Hz, and it will worsen (improve) with decreasing (increasing) LO power. We assume an otherwise ideal LO subject to this minimum noise floor by setting the phase noise to the constant valuẽ L o ð Þ ¼L min over the effective control bandwidth o ∈ [0, o c ]. Using this spectrum we calculate the resulting upper bound on operational fidelity in the presence of this noise using χ min κo c =2π ð Þ 10L min =10 . In this expression κ is a characteristic scaling factor depending on the control protocol, and χ min is approximately independent of τ in the limit o c /2π ≫ τ − 1 (Higherpower LOs can in principle have lower phase-noise floors, but will generally suffer from additional amplifier-induced phase noise. We therefore consider the thermal noise floor for a LO at 0 dBm output as an approximate lower-bound for the phase noise. We assume that the broadband thermal noise is actually manifested as temperature-dependent phase fluctuations in the LO rather than a background observed in the measurement process ( Supplementary Information)).
This straightforward calculation reveals that broadband LO phase noise due to thermal effects at room temperature imposes a non-negligible upper bound on gate fidelity, as shown in Figure 3. For instance, considering a typical bandwidth~20 GHz, the thermal noise floor induces qubit infidelity in excess of 10 − 8 . Reducing the thermal noise floor to that associated with a 4 K bath (−192 dBc/Hz) improves the fidelity by two orders of magnitude. Similarly, restricting the control bandwidth to 10 MHz sets the infidelity o10 − 11 .

DISCUSSION
The calculations we have presented demonstrate that master clock phase fluctuations are an emerging consideration for quantum information applications. Indeed, lab-grade oscillators may already limit the performance of today's qubits with gate times on microsecond scales. On the other hand, using high-performance precision LOs, the calculated error rates forl andX operations with both the fast and slow gates considered here are at several orders of magnitude smaller than the current state-of-art (Table 1). Therefore, although significantly more expensive, the use of such precision LOs is adequate for near-term, proof-of-concept experimental demonstrations with contemporary qubits.
Our results also enhance existing arguments about the merits of qubit modalities that accommodate fast control pulses. These arguments are typically based on a practical 'clock-speed' assessment, wherein a technology with faster gates will simply execute an algorithm more quickly. Here in considering LO phase noise, we find that reducing operation times from 10 μs to 10 ns generally also reduces infidelity (Figure 2c,d). The improvement can be substantial (3-4 orders of magnitude) for lab-grade LOs, and it remains about an order of magnitude for the precision LO considered here. The intuition is that shorter pulses are subject to phase fluctuations for less time and, therefore, suffer less dephasing-induced errors.
However, arguments in favour of shorter control operations must be made with a cognizance of the actual LO noise spectrum and the control bandwidth. First, shorter operations access higher frequencies in the LO spectrum, and, as illustrated in Figure 2, far-from-carrier phase noise dominates residual errors for commercial sources. Therefore, there may in principle be cases where decreasing the pulse duration actually increases the integrated noise level. Second, bandwidth-dependent thermal-noise floors will eventually pose upper bounds on operational fidelities that are more strict when using short control pulses possessing higher bandwidths. And, when they do, suppressing the effective thermal noise floor by further engineering the performance of room-temperature LOs and/or developing cryo-compatible LOs 56 embedded in cryogenic control architectures 57 may have an important role in the development of future, high-fidelity quantum information systems.
Similarly, whereas DES protocols are reasonably good at suppressing dephasing error due to low-frequency environmental fluctuations, they are less efficient at suppressing the impact of LO phase noise (Figure 2a) due to the high-frequency weight in the dephasing power spectrum (Figure 2b). In fact, over certain frequency bands in Figure 2c,d, the DES protocols (dashed lines) can actually enhance error due to LO phase noise. Noise at higher frequencies evolves too rapidly to be coherently averaged by the control, resulting in the filter transfer functions for the DES protocols passing noise in this band. The detailed spectrum ultimately determines whether shorter pulses and DES protocols yield a net win despite the higher noise far from the carrier.
Moving to mesoscale quantum information, such as quantum simulations and/or prototype logical-qubit demonstrations with a few hundred qubits, 58-60 we believe these results will serve to inform error budgets and motivate careful selection of the precision LO in use. Looking forward, however, the bounds on error rates calculated here for both free and driven evolution with precision oscillators are, in the view of these authors, remarkably high and will likely limit performance in larger scale systems, particularly in the context of quantum error correction.
It is widely viewed that quantum error correction, QEC, will be necessary to achieve sustained operation in faulty qubit systems. Error bounds due to LO phase instability (of the type calculated here) place architectural constraints on the error correction protocol: the QEC cycle must be completed within a time corresponding to a fixed phase-noise-induced error rate. The tabulated maximum QEC cycle times appearing in Table 2 show differences of 410 5 in the permissible QEC cycle times for the two oscillators considered here.
The use of a precision LO in this context readily supports phase noise infidelities at the level 10 − 7 , an order of magnitude below  These times may be viewed as an upper bound on the allowable QEC cycle period. Achievable cycle periods will be reduced because of other error sources. Error rates are derived from free-evolution calculations presented in Figure 2b. Driven operation error rates (see Figure 2c) yield similar results.
Master clock stability in quantum information H Ball et al the most pessimistic fault-tolerance thresholds, over a time period 600 μs. This provides strong evidence that, in principle, we should be able to achieve fault tolerance for QEC with existing oscillator technology. However, scalability requires consideration of corrected fault-tolerant error rates that are sufficiently low to permit algorithmic execution with manageable resource levels, 1 a target that grows increasingly challenging for larger scale applications. 2 Broadly, one would aim to have a logical error rate that scales as the inverse problem size, with only a small chance of a single logical error during the algorithmic execution. This results in calculated logical-qubit error rates reaching O 10 -15 À Á in the context of Shor's algorithm for factorisation of medium-sized keys. 61 Achieving such extraordinary logical-error-rate targets for largescale machines will involve a consideration of both the physical qubit error rate and the resource levels required to implement QEC. As error rates approach the fault-tolerance threshold, the encoding resource requirements diverge in time and qubit numbers. Stated in an alternate way, the efficiency of QEC grows as the error rate is suppressed relative to the fault-tolerance threshold. Accordingly, the error floors imposed by LO phase instability must be traded off at a system level against the overheads associated with QEC encoding, for instance, in the number of required physical qubits encoding a single logical bit. Calculations suggest that for the Bacon-Shor code, the presence of clock-induced errors near~10 − 7 would require a factor 100 − 1,000 times more resources relative to physical qubit error rates near 10 − 10 . Although the LO is only part of the story-errors due to environmental qubit dephasing must also be reducedone will also need to improve LOs to achieve these levels. 61 On the basis of current understanding of the relative challenge associated with QEC encoding in either large surface codes or concatenated schemes, and the continuing, anticipated improvements in physical qubit error rates, we believe that the clock-induced error rates we have identified here motivate investment in ultralow-phase noise LO development.
Another major consideration in the context of QEC is that clock-induced errors are highly correlated in space and time, posing a challenge to existing QEC analyses that have focused primarily on independent and identically distributed stochastic error models. As experimental teams incorporate clock distribution in quantum information architectures, 57,62,63 the complexity of QEC analyses will need to be augmented in order to handle the effects of such correlated errors.
Although we have focused on phase fluctuations far from the carrier-i.e., on short time scales such as the physical qubit gate time or the QEC cycle period-attention to long-term LO stability is also required. Slow phase diffusion of the LO-i.e., close to carrier phase noise-causes substantial error accumulation over long times. In the context of a single LO, such errors appear adiabatic and are generally correctable, provided the timescale of the phase diffusion is much longer than the QEC cycle period. By contrast, in systems with multiple control generators, long-time instabilities may be detrimental if the generators exhibit phase diffusion with respect to one another. In such cases, the diffusion represents a temporal, stochastic analogue to the spatial, deterministic 'clock skew' observed in classical semiconductor chips, and it is not generally accommodated by standard quantum error correction protocols. To give a concrete example, whereas generators A and B may individually exhibit correctable slow phase diffusion, the axis for anX gate on generator A may adiabatically skew to be (in an extreme case) in quadrature with the equivalentX gate on generator B. A logical controller directs generator A to apply anX gate to a qubit A, but the rotation is â Y rotation relative to generator B and the qubits it drives.
Temporally correlated low-frequency LO phase fluctuations can also introduce challenges in evaluating relevant gate errors for QEC because of the presence of potential biases in common quantum verification processes such as randomised benchmarking. 64 In evaluating the impact of such errors, it is important to distinguish the contribution of slow LO instabilities to single-operation error rates from their contribution to empirical estimates of qubit fidelities drawn from tomographic measurements over repeated trials. The process of data acquisition and averaging over many individual experiments can contribute error due to slow drifts that dominate the actual error rate experienced in any individual operation or evolution period (see ref. 33 for details).
In summary, we have presented a comprehensive review of critical issues relating to master-clock-induced errors for quantum information applications. We have unified established concepts from frequency metrology and quantum information, permitting experimentalists to translate common LO hardware specifications into estimates of qubit coherence and operational fidelities. Analyses employing information available on representative LO data sheets have revealed that far-from-carrier phase noise poses a performance-limiting upper-bound on operational fidelities. As a result we expect to see a growing emphasis on high-performance LO synthesis chains across all technology platforms, and we foresee a growing importance of LO-induced errors in the design of large-scale quantum information systems. We encourage future studies to examine not only clock synthesis but also clock distribution, with an eye towards architectural impacts of skewed clock distribution in quantum systems with single and multiple LOs.

METHODS
In the appropriate interaction picture, environmental dephasing processes and those induced by LO phase fluctuations are formally equivalent. We begin with the physical system Hamiltonian (ħ = 1) in the laboratory frame for a qubit possessing a nominal transition (angular) frequency o 0 , driven by a local oscillator with carrier frequency o LO , Appearing in this expression are the operatorsσ i representing the Pauli matrices in the Cartesian coordinate basis. The first term corresponds to the Hamiltonian of the qubit under free evolution, whereas the second term captures environmental noise (e.g., magnetic field fluctuations) that effectively changes the nominal qubit transition frequency by an amount δo(t) as The third term corresponds to the qubit-field interaction with amplitude Ω(t) driven by a local oscillator with nominal carrier frequency o LO . The time-dependent phase of the LO is partitioned into a 'control' component (desired phase evolution), ϕ C (t), and a 'noise' component (unwanted phase fluctuations), ϕ N (t).
Moving to the interaction picture co-rotating with the carrier frequency and making the rotating wave approximation we obtain where we define the qubit ladder operatorsσ 7 ¼σ x 7 iσ y . The LO phase fluctuation produces a time-dependent frequency detuning through its time derivative, δo LO t ð Þ _ ϕ N t ð Þ, effectively transforming the LO frequency as In order to make this phenomenology explicitly comparable to environmental dephasing, we perform a second interactionpicture transformation