Skip to content

Optical Clocks

An optical clock compares an optical local oscillator with a narrow atomic transition and steers the oscillator so that its frequency tracks the unperturbed transition frequency. The clock transition is typically in the range 101410^{14} to 1015 Hz10^{15}\ \mathrm{Hz}, about five orders of magnitude above the caesium microwave transition that currently defines the SI second.

The higher carrier frequency is valuable because a fixed interrogation time can resolve a similar absolute linewidth at a much smaller fractional linewidth. If

ν0=Ee−Egh,Q=ν0Δν,\nu_0 = \frac{E_e-E_g}{h}, \qquad Q = \frac{\nu_0}{\Delta\nu},

then a 1 Hz1\ \mathrm{Hz} line at 4.3×1014 Hz4.3\times10^{14}\ \mathrm{Hz} has Q≃4.3×1014Q\simeq4.3\times10^{14}. This is an opportunity, not a complete performance guarantee. A practical optical clock also requires:

  • an ultrastable laser that preserves phase during interrogation;
  • trapped atoms or ions whose motion and environment are controlled;
  • state preparation, coherent interrogation, and efficient readout;
  • a servo that converts atomic excitation into frequency correction;
  • an optical frequency comb or coherent optical link for usable output;
  • a correction model for Stark, Zeeman, motional, collisional, relativistic, and servo effects; and
  • independent comparisons capable of exposing unrecognized bias.

The useful output is therefore a disciplined optical oscillator plus its measurement model, not an isolated atom and not merely a narrow spectral line.

Atomic Clocks owns the generic passive clock loop, microwave architectures, clock stability, quantum projection noise, dead time, the Dick effect, accumulated time error, and the present caesium definition of the second. Ramsey Interferometry derives the separated-pulse fringe, discriminator, sensitivity function, and finite-pulse corrections. Frequency Combs derives comb-tooth frequencies, self-referencing, transfer-oscillator combinations, and optical frequency division.

Ion Traps owns Paul-trap and Penning-trap confinement, secular motion, micromotion, stability parameters, cooling, and heating. Optical Lattices owns periodic optical potentials, recoil scales, bands, tunneling, and many-body lattice physics. AC Stark Shift and Dynamic Polarizability own the general light-shift and polarizability derivations.

This page owns the optical-clock-specific synthesis:

  1. selection of optical clock transitions and species;
  2. the different roles of the atomic line, clock laser, servo, comb, and comparison link;
  3. trapped-ion, quantum-logic, and neutral-atom lattice architectures;
  4. the magic-wavelength condition and residual lattice shifts;
  5. optical-clock systematic ledgers, especially blackbody radiation, micromotion, time dilation, collisions, and electric-quadrupole shifts;
  6. absolute optical-frequency measurements and optical frequency ratios;
  7. clock networks, relativistic geodesy, and potential measurements;
  8. frequency-ratio tests of constants, relativity, and selected dark-matter models; and
  9. the current relation of optical clocks to the SI second and its possible future redefinition.

Frequency Standards owns international traceability, primary and secondary standards, time scales, and dissemination in detail. Tests of Fundamental Symmetries, Fundamental Constants, and Variation of Constants Searches own the broader theory and experimental landscape of those subjects.

Separate three frequencies:

  1. νatom\nu_{\mathrm{atom}} is the perturbed atomic transition sampled in a particular cycle.
  2. νLO\nu_{\mathrm{LO}} is the optical local-oscillator frequency delivered to the atoms after all controlled frequency shifts.
  3. νout\nu_{\mathrm{out}} is the reported or distributed clock frequency after steering and systematic corrections.

For a cycle indexed by kk, a simplified error signal is

ek≃D[νLO,k−νatom,k]+nk,e_k \simeq D \left[ \nu_{\mathrm{LO},k} - \nu_{\mathrm{atom},k} \right] + n_k,

where DD is the local discriminator slope and nkn_k includes projection, detection, and technical noise. A digital servo updates a control word, while a separate correction ledger connects the observed resonance to the specified unperturbed transition:

ν^0=ν‾LO−∑iΔνi^.\widehat{\nu}_0 = \overline{\nu}_{\mathrm{LO}} - \sum_i \widehat{\Delta\nu_i}.

This expression defines Δνi\Delta\nu_i as a physical shift of the observed transition. Some clock groups instead publish corrections to be added to the measured frequency. Either convention is valid; mixing the two changes the answer’s sign.

Comparison of trapped-ion and optical-lattice clocks, the coherent optical readout chain, and the principal systematic and relativistic coordinates

Modern optical clocks combine several logically distinct subsystems. A: a trapped-ion clock interrogates one or a few strongly confined particles, whereas a lattice clock interrogates many neutral atoms near the Lamb–Dicke regime at a magic trapping frequency. B: a cavity-stabilized laser is steered by atomic interrogation; a self-referenced frequency comb connects the optical output to another optical reference or to microwave electronics. C: ion and lattice clocks have different leading systematic coordinates, while every remote comparison also requires coherent transfer and a gravitational-potential correction.

A representative cycle contains:

  1. Load and cool. Capture an ion or neutral-atom ensemble and reduce motion to the regime required by the uncertainty target.
  2. Prepare. Pump into a selected Zeeman or hyperfine clock state.
  3. Interrogate. Apply a phase-coherent Rabi, Ramsey, hyper-Ramsey, or related sequence with the clock laser.
  4. Detect. Infer the final state, often by electron shelving or state-dependent fluorescence.
  5. Discriminate. Probe opposite sides of the resonance, or use a phase-stepped sequence, to form an error signal.
  6. Steer. Update the laser frequency through an acousto-optic or electro-optic actuator.
  7. Evaluate. Interleave changed experimental conditions to estimate shifts and monitor auxiliary fields.
  8. Compare. Count a comb-derived beat, optical ratio, or coherent link observable with defined gates and signs.

The atoms need not be present continuously. During loading, cooling, and detection, the pre-stabilized laser carries the phase. The clock’s short-term stability can therefore be limited by local-oscillator noise and dead-time aliasing even when the atoms themselves are nearly ideal.

For uncorrelated atoms, a useful order-of-magnitude instability expression is

σy(τ)≃12πν0TCTcNτ,\sigma_y(\tau) \simeq \frac{1}{ 2\pi\nu_0 T C } \sqrt{ \frac{T_c}{ N\tau } },

where TT is the coherent interrogation time, CC is fringe contrast, TcT_c is the cycle time, and NN is the number of detected atoms per cycle. The exact numerical factor depends on the sequence, operating point, and detection model. The important optical-clock advantages are visible:

  • large ν0\nu_0 lowers the fractional error corresponding to a fixed phase error;
  • narrow transitions permit long coherent interrogation;
  • laser cooling and trapping suppress first-order Doppler broadening;
  • carefully chosen clock states reduce sensitivity to fields; and
  • frequency combs preserve phase while connecting optical and electronic frequencies.

The formula also exposes the platform trade-off. A single ion has N=1N=1 but excellent isolation and readout. A lattice clock may interrogate 10210^2 to 10410^4 atoms, improving statistical precision, but must control collisions, inhomogeneity, and collective effects.

A clock transition is judged as part of an instrument, not solely by its natural lifetime. Desirable properties include:

  • a long-lived upper state and a linewidth narrower than the achievable interrogation linewidth;
  • a wavelength for which coherent lasers and optical components are available;
  • cooling, trapping, state-preparation, and detection transitions that can be operated reproducibly;
  • clock states with small or characterizable differential electric and magnetic response;
  • suppression of first-order Zeeman and electric-quadrupole shifts by angular-momentum choice or averaging;
  • low sensitivity to blackbody radiation or accurately calculable polarizability;
  • manageable motion, recoil, collisions, and probe shifts;
  • independent clocks and ratio links for reproducibility; and
  • atomic-structure sensitivity appropriate to any proposed fundamental-physics measurement.

The natural quality factor

Qnat=ν0Γ/(2π)Q_{\mathrm{nat}} = \frac{\nu_0}{\Gamma/(2\pi)}

can be enormous, but a finite pulse of duration TT produces a linewidth of order 1/T1/T. The observed line can instead be limited by laser coherence, residual motion, magnetic noise, inhomogeneous shifts, or decoherence. Quoting QnatQ_{\mathrm{nat}} as the achieved clock quality factor is therefore misleading unless lifetime-limited interrogation has actually been demonstrated.

The following table gives rounded frequencies and wavelengths to identify important platforms. It is not a substitute for the isotope, hyperfine, Zeeman, polarization, and averaging conventions in a clock report.

SystemRepresentative clock transitionFrequencyWavelengthTypical architecture
27Al+^{27}\mathrm{Al}^{+}1S0↔3P0^1S_0\leftrightarrow{}^3P_01121 THz1121\ \mathrm{THz}267 nm267\ \mathrm{nm}quantum-logic ion
199Hg+^{199}\mathrm{Hg}^{+}2S1/2↔2D5/2^2S_{1/2}\leftrightarrow{}^2D_{5/2}1065 THz1065\ \mathrm{THz}282 nm282\ \mathrm{nm}single ion
171Yb+^{171}\mathrm{Yb}^{+}electric-octupole line642 THz642\ \mathrm{THz}467 nm467\ \mathrm{nm}single ion
171Yb+^{171}\mathrm{Yb}^{+}electric-quadrupole line688 THz688\ \mathrm{THz}436 nm436\ \mathrm{nm}single ion
88Sr+^{88}\mathrm{Sr}^{+}2S1/2↔2D5/2^2S_{1/2}\leftrightarrow{}^2D_{5/2}445 THz445\ \mathrm{THz}674 nm674\ \mathrm{nm}single ion
40Ca+^{40}\mathrm{Ca}^{+}2S1/2↔2D5/2^2S_{1/2}\leftrightarrow{}^2D_{5/2}411 THz411\ \mathrm{THz}729 nm729\ \mathrm{nm}single ion
199Hg^{199}\mathrm{Hg}1S0↔3P0^1S_0\leftrightarrow{}^3P_01129 THz1129\ \mathrm{THz}265 nm265\ \mathrm{nm}optical lattice
171Yb^{171}\mathrm{Yb}1S0↔3P0^1S_0\leftrightarrow{}^3P_0518 THz518\ \mathrm{THz}578 nm578\ \mathrm{nm}optical lattice
87Sr^{87}\mathrm{Sr}1S0↔3P0^1S_0\leftrightarrow{}^3P_0429 THz429\ \mathrm{THz}698 nm698\ \mathrm{nm}optical lattice
88Sr^{88}\mathrm{Sr}magnetically enabled 1S0↔3P0^1S_0\leftrightarrow{}^3P_0429 THz429\ \mathrm{THz}698 nm698\ \mathrm{nm}optical lattice

Several of these radiations are secondary representations of the second in the BIPM list. That status means they can contribute to practical realizations and calibrations under specified procedures; it does not mean that each frequency is exact or that the SI second already has an optical definition.

As of 25 July 2026, the SI second remains defined by assigning the exact value

ΔνCs=9 192 631 770 Hz\Delta\nu_{\mathrm{Cs}} = 9\,192\,631\,770\ \mathrm{Hz}

to the unperturbed ground-state hyperfine transition of 133Cs^{133}\mathrm{Cs}. Optical standards are measured against that definition directly or through an international least-squares network of absolute frequencies and optical ratios.

The CCTF approved an updated list of recommended standard frequencies in 2025, supported by published values and a covariance-aware adjustment. The BIPM roadmap discusses a possible redefinition in 2030 or later, conditional on agreed criteria for clock performance, independent reproducibility, frequency ratios, time and frequency transfer, operational reliability, and international access. A candidate transition, an ensemble of transitions, and the continuity strategy are matters of metrological decision, not consequences of a single laboratory record.

Strong electric-dipole transitions are excellent for cooling and fluorescence, but their short excited-state lifetimes make them too broad for the highest-resolution clocks. Clock lines are instead weak because a leading transition amplitude is forbidden or strongly suppressed by angular momentum, parity, spin, or multipole selection rules.

Examples include:

  • J=0↔J=0J=0\leftrightarrow J=0 intercombination transitions between 1S0^1S_0 and 3P0^3P_0 states in group-II-like atoms and ions;
  • electric-quadrupole transitions such as S1/2↔D5/2S_{1/2}\leftrightarrow D_{5/2};
  • an electric-octupole transition in 171Yb+^{171}\mathrm{Yb}^{+}; and
  • transitions made weakly allowed by hyperfine or magnetic mixing.

The same weak coupling that creates a narrow natural line makes coherent excitation technically demanding. Probe intensity, polarization, magnetic field, pulse area, phase noise, and off-resonant couplings must be known well enough that the interrogation itself does not dominate the clock.

For fermionic isotopes such as 87Sr^{87}\mathrm{Sr} and 171Yb^{171}\mathrm{Yb}, hyperfine interaction mixes a small amount of J=1J=1 character into the nominal 3P0^3P_0 state. The 1S0↔3P0^1S_0\leftrightarrow{}^3P_0 line is then weakly electric-dipole allowed. The nuclear spin also creates multiple Zeeman components. Clock operation usually prepares selected mFm_F states and averages opposite projections to cancel the first-order Zeeman shift while retaining a magnetic-field diagnostic.

Even isotopes such as 88Sr^{88}\mathrm{Sr} have nuclear spin zero. Their J=0↔J=0J=0\leftrightarrow J=0 line lacks hyperfine-induced coupling and can be enabled by applying a static magnetic field that mixes the excited state with a nearby J=1J=1 state. This offers simple level structure but introduces probe and magnetic-field dependencies that must be calibrated jointly.

These statements are architecture-level rules, not universal rankings. Fermions can suffer density-dependent pp-wave interactions and imperfect spin polarization; bosons can require larger enabling fields and can experience ss-wave collisions. The optimum isotope depends on atom number, geometry, temperature, interrogation protocol, and target uncertainty.

During interrogation the local oscillator must preserve phase relative to the atomic superposition. A common chain is:

  1. lock a laser to a high-finesse reference cavity;
  2. suppress vibration, temperature, residual-amplitude, and optical-path noise;
  3. transfer that coherence to the atoms through a phase-stabilized path;
  4. steer the slowly varying laser frequency using the atomic error signal;
  5. compare or divide the disciplined light with a comb.

The cavity is a short-term flywheel. It does not provide the long-term accuracy of the clock transition. Thermal noise, drift, vibration, optical power, residual gas, and etalon effects move its resonance. Conversely, the atoms cannot correct laser phase between interrogations. Hybrid operation is essential.

For Rabi interrogation, the Fourier-limited line shape depends on pulse area and detuning. For Ramsey interrogation, the central-fringe scale is set mainly by the dark time. Rabi Oscillations and Ramsey Interferometry give the canonical derivations. Optical-clock reports should state the actual pulse duration, sequence, duty factor, contrast, atom number, and laser-noise conditions rather than infer them from the natural lifetime.

Weak clock lines can require non-negligible probe intensity. Off-resonant coupling then produces an AC Stark shift. A sequence may reduce sensitivity by extrapolating to zero intensity, alternating intensities, using hyper-Ramsey phases, or adding a calibrated frequency step during the pulses.

Composite interrogation does not make the shift disappear automatically. Its cancellation relies on pulse area, phase, decoherence, frequency-step calibration, and servo symmetry. A valid evaluation varies those coordinates and tests the residual with an out-of-loop comparison.

A radio-frequency Paul trap or static-field Penning trap can confine an ion for long periods. Laser cooling places its secular motion near the Lamb–Dicke regime, so first-order Doppler broadening is suppressed. Strong cycling transitions provide nearly unit-efficiency state detection through electron shelving.

A representative single-ion sequence is:

  1. Doppler cool and, when required, sideband cool;
  2. optically pump into a clock-state Zeeman sublevel;
  3. turn off or control perturbing cooling light;
  4. interrogate the clock transition;
  5. detect whether fluorescence is bright or dark;
  6. repump and repeat at the opposite discriminator point.

The single particle avoids density shifts and makes state control clean. Its statistical limitation is also clear: one binary outcome is obtained per cycle. Improved laser coherence, longer interrogation, rapid cooling, high duty factor, and correlation with another clock can reduce the averaging time.

Some excellent clock ions lack a convenient cycling transition. In quantum-logic spectroscopy, the clock ion is co-trapped with a logic ion. Shared motional modes provide a bus:

  1. the logic ion sympathetically cools the pair;
  2. the clock laser maps clock-state information to shared motion;
  3. logic operations map that motion to the logic ion;
  4. fluorescence of the logic ion reads out the clock state.

The canonical example is 27Al+^{27}\mathrm{Al}^{+} with a co-trapped 25Mg+^{25}\mathrm{Mg}^{+} logic ion. This architecture separates the species best suited to the reference transition from the species best suited to cooling and detection. It adds mode-frequency drift, heating, state-mapping fidelity, differential motion, collisions, and logic-pulse errors to the validation problem.

Quantum logic is not merely an indirect fluorescence detector. The mapping sequence and shared motion are part of the clock measurement model. A reported uncertainty must show that state-dependent mapping error does not move the lock point.

Important ion-specific effects include:

  • excess micromotion: displacement from the radio-frequency null creates driven motion and both time-dilation and rf Stark shifts;
  • secular-motion time dilation: residual thermal and driven motion gives a relativistic second-order Doppler shift;
  • electric-quadrupole shift: a state with nonzero quadrupole moment couples to static electric-field gradients;
  • trap-induced AC Zeeman shift: rf currents and fields can perturb magnetic sublevels;
  • background-gas collisions: rare events may shift phase, reorder an ion pair, or eject the ion;
  • heating and mode drift: motion can change between cooling and interrogation; and
  • logic-ion perturbations: cooling light, motion, and state-mapping operations can affect the clock ion.

For nonrelativistic speed v≪cv\ll c, the time-dilation contribution is

ΔνTDν0=−⟨v2⟩2c2.\frac{\Delta\nu_{\mathrm{TD}}}{\nu_0} = - \frac{ \langle v^2\rangle }{ 2c^2 }.

The relevant ⟨v2⟩\langle v^2\rangle includes secular motion, intrinsic micromotion, and excess micromotion under a stated model. It cannot be replaced by the Doppler-cooling temperature alone.

An electric-quadrupole shift has the angular structure

ΔνQ∝Θ[3mJ2−J(J+1)][3cos⁡2β−1]∇E,\Delta\nu_Q \propto \Theta \left[ 3m_J^2-J(J+1) \right] \left[ 3\cos^2\beta-1 \right] \nabla E,

where Θ\Theta is the state quadrupole moment, β\beta relates the quantization axis to a principal gradient axis, and ∇E\nabla E represents the relevant traceless field-gradient component. Exact prefactors depend on tensor and gradient conventions. Averaging suitable Zeeman components or three orthogonal quantization axes can cancel a static rank-2 contribution, provided field directions, gradients, and weighting remain controlled.

Neutral atoms cannot be confined as strongly by static electric fields. An optical lattice instead supplies a standing-wave AC Stark potential. Atoms are cooled and loaded into many lattice sites, where confinement suppresses recoil and first-order Doppler effects. Hundreds or thousands of atoms can be interrogated simultaneously.

The central design condition is that the ground and excited clock states experience equal leading-order light shifts. For lattice angular frequency ωL\omega_L, the electric-dipole magic condition is

ΔαE1(ωm)=αe(ωm)−αg(ωm)=0.\Delta\alpha_{\mathrm{E1}}(\omega_m) = \alpha_e(\omega_m) - \alpha_g(\omega_m) = 0.

At the magic frequency ωm\omega_m, the leading scalar E1 differential shift vanishes. The lattice still traps both states because their common polarizability need not vanish.

Magic operation does not imply zero lattice uncertainty. Residual terms can come from:

  • detuning from the operational magic frequency;
  • hyperpolarizability, proportional to higher powers of intensity;
  • electric-quadrupole and magnetic-dipole lattice couplings;
  • vector and tensor polarizabilities;
  • imperfect polarization and magnetic-field alignment;
  • motional-state dependence and tunneling;
  • intensity inhomogeneity and longitudinal or transverse temperature;
  • multiphoton or near-resonant structure; and
  • correlations among fitted lattice-shift coefficients.

With the convention

E(t)=E0cos⁡(ωLt),E(t) = E_0\cos(\omega_L t),

a schematic differential shift is

Δνlatt=−ΔαE1(ωL)4hE02−Δβ(ωL)64hE04+Δνqm+Δνmot+⋯ .\Delta\nu_{\mathrm{latt}} = - \frac{ \Delta\alpha_{\mathrm{E1}}(\omega_L) }{ 4h } E_0^2 - \frac{ \Delta\beta(\omega_L) }{ 64h } E_0^4 + \Delta\nu_{\mathrm{qm}} + \Delta\nu_{\mathrm{mot}} + \cdots.

Here Δβ\Delta\beta denotes a convention-dependent differential hyperpolarizability, Δνqm\Delta\nu_{\mathrm{qm}} collects multipolar terms, and Δνmot\Delta\nu_{\mathrm{mot}} collects motion-dependent corrections. The numerical coefficients change if intensity or rms electric field is used instead of peak field. A clock paper must define the convention rather than transport coefficients blindly.

Large atom number improves stability but introduces interactions. The observed shift can depend on:

  • site occupancy and density;
  • spin polarization and isotope statistics;
  • excitation fraction;
  • temperature and motional-state population;
  • tunneling and lattice dimensionality;
  • inhomogeneous Rabi coupling; and
  • the interrogation sequence.

For identical spin-polarized fermions, ss-wave collisions are suppressed, but imperfect indistinguishability and pp-wave interactions remain. Bosonic isotopes permit ss-wave interactions. A common evaluation interleaves high and low atom number and fits a density coefficient, but a single linear slope is insufficient if excitation fraction or temperature co-varies with atom number.

Three-dimensional lattices can isolate atoms site by site and suppress collisions, at the cost of added lattice fields and more complex loading. One-dimensional lattices offer simpler optical geometry and high atom number, but may require stronger density and interaction controls.

CoordinateTrapped-ion clockOptical-lattice clock
Particle numberusually one clock iontypically many neutral atoms
Statistical strengthone high-fidelity binary result per cycleensemble signal and fast averaging
Confinementelectromagnetic ion trapstate-insensitive optical lattice
Density shiftabsent for a single clock ionrequires interaction evaluation
Motionsecular motion and micromotionband, tunneling, temperature, recoil
Trap-light shiftusually no lattice fieldmagic condition plus residual lattice terms
Important tensor shiftelectric quadrupole for suitable statesvector/tensor lattice and Zeeman terms
Readoutelectron shelving; sometimes quantum logicstate-selective population detection
Duty factorcooling and readout can be substantialloading losses and destructive readout can be substantial
Strengthclean control and low interactionshigh atom number and rapid stability
Main scaling challengesingle-particle statistics and trap fieldsinteractions, inhomogeneity, and lattice model

No column is intrinsically more accurate. Uncertainty belongs to an evaluated realization under specified operating conditions. A clock with more favorable atomic coefficients can still perform worse if its fields, temperature, motion, or comparisons are inadequately controlled.

For independent clocks, local-oscillator noise can dominate a comparison. Several strategies reduce it:

  • increase laser coherence and interrogation duty factor;
  • interrogate two ensembles synchronously with a common laser;
  • use a second ensemble to estimate laser phase;
  • operate zero-dead-time interleaved ensembles;
  • transfer stability from a remote cryogenic or ultralow-noise cavity; or
  • use entanglement or spin squeezing while accounting for loss and decoherence.

Synchronous comparison rejects laser noise only to the extent that the two systems share the same optical phase and sensitivity function. Differential path noise, unmatched pulse timing, unequal transfer functions, cycle slips, and asynchronous dead time reintroduce noise. Common-mode rejection can also hide a shared systematic, so absolute evaluations and independent links remain necessary.

A self-referenced comb has teeth

νn=nfrep+fCEO.\nu_n = n f_{\mathrm{rep}} + f_{\mathrm{CEO}}.

Let the clock laser beat with tooth nn at a signed beat frequency fbf_b, and let acousto-optic modulators contribute signed shifts ∑jsjfAOM,j\sum_j s_j f_{\mathrm{AOM},j}. One explicit convention is

νclock=nfrep+fCEO+sbfb+∑jsjfAOM,j,sb,sj∈{−1,+1}.\nu_{\mathrm{clock}} = n f_{\mathrm{rep}} + f_{\mathrm{CEO}} + s_b f_b + \sum_j s_j f_{\mathrm{AOM},j}, \qquad s_b,s_j\in\{-1,+1\}.

The sign is determined by the physical optical-frequency ordering and diffraction path, not by whether a counter displays a positive number. A complete frequency record stores:

  • the integer tooth index nn;
  • the definitions and signs of all beat notes;
  • every fixed and time-dependent frequency shift;
  • counter mode and gate timing;
  • the comb locks and reference chain;
  • cycle-slip diagnostics; and
  • the optical reference plane to which the equation applies.

At νclock∼5×1014 Hz\nu_{\mathrm{clock}}\sim5\times10^{14}\ \mathrm{Hz}, a one-tooth error is of order frepf_{\mathrm{rep}}, often hundreds of megahertz. A sign error in a 40 MHz40\ \mathrm{MHz} AOM is also enormous by clock standards. These are not small systematic corrections.

An absolute frequency measurement expresses an optical transition in SI hertz. It must ultimately connect to the current caesium definition, usually through a microwave reference, a local time scale, satellite or fiber transfer, and primary or secondary standard calibrations. The comparison uncertainty and dead time can be much larger than the intrinsic optical-clock uncertainty.

An optical frequency ratio

RAB=νAνBR_{AB} = \frac{\nu_A}{\nu_B}

can be measured directly with a comb while avoiding much of the microwave chain. Because a ratio of frequencies is dimensionless, it is also the natural observable for comparing atomic structure and searching for changes in dimensionless constants.

Direct ratios do not eliminate all conventions. Both frequencies need systematic and gravitational corrections to declared reference conditions. Comb transfer noise, uncompensated paths, counter synchronization, data selection, and covariance among clocks must be evaluated.

Two lasers at νA\nu_A and νB\nu_B produce beat signals against teeth nAn_A and nBn_B. Suitable electronic or digital combinations can cancel comb fluctuations in frepf_{\mathrm{rep}} and fCEOf_{\mathrm{CEO}} without requiring the comb itself to be quieter than either optical source. In schematic form, construct a transfer observable

ft=(sAfb,A+fCEO)−nAnB(sBfb,B+fCEO),f_t = \left( s_A f_{b,A} + f_{\mathrm{CEO}} \right) - \frac{n_A}{n_B} \left( s_B f_{b,B} + f_{\mathrm{CEO}} \right),

with signs and AOM terms adapted to the actual optical layout. Substitution of the tooth equations removes the repetition-rate noise. The exact implementation must account for electronic delays, finite bandwidth, integer ratios, synthesizer phase noise, and sampled data timing.

Frequency Combs gives the full derivation. The operational lesson for clocks is that a comb is a coherent gearbox: it transfers a reference and its noise. It does not independently certify either clock.

For three co-located, simultaneously defined frequencies,

RABRBCRCA=1.R_{AB} R_{BC} R_{CA} = 1.

A measured logarithmic closure residual is

ϵcl=ln⁡R^AB+ln⁡R^BC+ln⁡R^CA.\epsilon_{\mathrm{cl}} = \ln\widehat R_{AB} + \ln\widehat R_{BC} + \ln\widehat R_{CA}.

Closure is a powerful check on comb signs, transfer chains, and ratio consistency. It is not a proof that all clocks are unbiased. A shared systematic shift can cancel from closure, and ratios measured at different times require a model for drift and covariance.

Remote clocks can be compared through actively stabilized optical fiber, free-space optical time-frequency transfer, transportable clocks, or satellite-based methods. The comparison observable includes:

  1. the two local clock corrections;
  2. local comb and synthesis chains;
  3. link phase-noise cancellation and residual asymmetry;
  4. cycle-slip and data-validity monitors;
  5. synchronized counter gates or a model for asynchronous sampling; and
  6. the gravitational-potential difference between reference locations.

At the 10−1810^{-18} level, the optical path through a laboratory, a patch cord exposed to temperature, or an incorrectly defined endpoint can matter. The phrase “fiber link uncertainty” is incomplete unless the physical reference planes and both local branches are included.

Write the fractional observation model as

yobs=νobs−νrefνref=y0+∑iδi+ϵ,y_{\mathrm{obs}} = \frac{ \nu_{\mathrm{obs}}-\nu_{\mathrm{ref}} }{ \nu_{\mathrm{ref}} } = y_0 + \sum_i \delta_i + \epsilon,

where y0y_0 is the target fractional frequency, δi\delta_i are systematic shifts, and ϵ\epsilon is statistical fluctuation under a specified averaging protocol. If the estimated corrections are collected in δ^\widehat{\boldsymbol\delta}, then

y^0=y‾obs−1Tδ^.\widehat y_0 = \overline y_{\mathrm{obs}} - \boldsymbol 1^{\mathsf T} \widehat{\boldsymbol\delta}.

For covariance matrix Σδ\Sigma_\delta, the systematic variance of the sum is

usys2=1TΣδ1.u_{\mathrm{sys}}^2 = \boldsymbol 1^{\mathsf T} \Sigma_\delta \boldsymbol 1.

Optical-clock coefficients are often shared across days, apparatuses, or species. Treating all entries as independent can understate or overstate uncertainty. The covariance model should distinguish measured environmental coordinates, fitted coefficients, atomic-theory inputs, and common calibrations.

EffectRepresentative coordinateTypical validation
Blackbody radiationtemperature field, static polarizability, dynamic correctionthermal map, sensor calibration, heated or cryogenic test
DC Starkstray electric field, electrode chargefield reversal, compensation scan
Lattice Starkdepth, frequency, polarization, motional stateinterleaved depth/frequency scan and model fit
Probe Starkpulse intensity and timingintensity extrapolation or composite-sequence test
Zeemanfield magnitude, direction, sublevelopposite-mm averaging and field scan
Electric quadrupolefield gradient and quantization axissublevel or three-axis averaging
Time dilationsecular motion, micromotion, temperaturesideband thermometry and micromotion diagnostics
Collisionsdensity, excitation, background gasdensity interleave, collision monitoring
Servodiscriminator asymmetry, drift, gainmodulation reversal and synthetic injection
Line pullingneighboring transitions, unwanted populationsstate purification and spectrum scan
Gravitationalgeopotential at the atomic reference pointgeodetic survey and potential model
Transfer chaincomb, AOM, fiber, countersredundant beats, closure, synchronized gates

A correction can be small while its uncertainty is large, and a large correction can be known accurately. Ranking effects only by correction magnitude is therefore unsound.

Thermal electromagnetic radiation produces a differential AC Stark shift. For an approximately isotropic environment at temperature TT, a common parameterization is

ΔνBBR(T)=−Δα(0)2h⟨E2(T)⟩[1+η(T)],\Delta\nu_{\mathrm{BBR}}(T) = - \frac{ \Delta\alpha(0) }{ 2h } \left\langle E^2(T)\right\rangle \left[ 1+\eta(T) \right],

where Δα(0)\Delta\alpha(0) is the differential static polarizability, η(T)\eta(T) is a dynamic correction, and

⟨E2(T)⟩∝T4.\left\langle E^2(T)\right\rangle \propto T^4.

The proportionality is exact for an ideal blackbody field when expressed through the radiation energy density; apertures, emissivity, view factors, temperature gradients, and non-equilibrium surfaces make the apparatus model more complicated.

If the T4T^4 term dominates and coefficient uncertainty is neglected, then

u ⁣(∣ΔνBBR∣)∣ΔνBBR∣≃4u(T)T.\frac{ u\!\left( \left|\Delta\nu_{\mathrm{BBR}}\right| \right) }{ \left|\Delta\nu_{\mathrm{BBR}}\right| } \simeq 4 \frac{u(T)}{T}.

This derivative is useful for planning but not a complete BBR budget. One must also propagate uncertainties in polarizability, dynamic correction, effective radiation temperature, gradients, emissivity, and the atom’s view of hot ovens, windows, electrodes, or apertures.

Three broad strategies are used:

  • operate near room temperature and measure the full thermal environment;
  • surround the atoms with a well-characterized radiation shield; or
  • operate cryogenically to reduce both the shift and its temperature sensitivity.

Temperature sensors report their own locations. The atoms sample a radiative field assembled from all visible surfaces. The two quantities coincide only under a validated thermal and radiative model.

A static field produces, to leading scalar order,

ΔνDC=−Δα(0)2hEDC2.\Delta\nu_{\mathrm{DC}} = - \frac{ \Delta\alpha(0) }{ 2h } E_{\mathrm{DC}}^2.

Patch charge on dielectrics, ultraviolet illumination, electrode contamination, and photoemission can create slowly varying fields. Because the leading shift is quadratic, simply reversing one electrode voltage does not necessarily reverse the shift. A compensation scan in several directions, together with a field model, is more informative.

Near the magic frequency, expand the differential E1 polarizability:

ΔαE1(ωL)≃∂ΔαE1∂ω∣ωm(ωL−ωm).\Delta\alpha_{\mathrm{E1}}(\omega_L) \simeq \left. \frac{ \partial\Delta\alpha_{\mathrm{E1}} }{ \partial\omega } \right|_{\omega_m} \left( \omega_L-\omega_m \right).

The leading residual shift is therefore bilinear in lattice intensity and magic-frequency detuning. At high depth, hyperpolarizability and multipolar terms can be comparable to the target uncertainty. An operational magic condition may deliberately choose a detuning and depth where several terms have reduced first derivative, rather than identifying one universal zero independent of motional state and intensity.

A robust evaluation samples a grid of depths and frequencies, not only two depths at one nominal wavelength. It records polarization, atom temperature, axial and radial occupation, and correlations among fitted coefficients. Extrapolating a polynomial outside the measured operating region is especially hazardous.

For a clock component labeled by mm, a schematic expansion is

ΔνZ(m,B)=k1mB+k2B2+⋯ .\Delta\nu_Z(m,B) = k_1 m B + k_2 B^2 + \cdots.

Averaging symmetric components gives

ν(+m)+ν(−m)2=ν0+k2B2+⋯ ,\frac{ \nu(+m)+\nu(-m) }{2} = \nu_0 + k_2B^2 + \cdots,

while their difference monitors the field:

ν(+m)−ν(−m)2≃k1mB.\frac{ \nu(+m)-\nu(-m) }{2} \simeq k_1mB.

This cancellation assumes the two components sample the same field, interrogation, and time. Magnetic drift between sequential probes, unequal line pulling, polarization error, tensor shifts, and servo weighting can leave a residual. The quadratic shift remains and must be measured or calculated.

The Lamb–Dicke regime suppresses first-order Doppler effects because the particle explores a region small compared with the probe wavelength. It does not eliminate:

  • relativistic time dilation;
  • recoil and motional-state dependence;
  • phase shifts from wave-front curvature or beam misalignment;
  • micromotion in an rf ion trap;
  • lattice tunneling and band effects; or
  • Doppler shifts from uncompensated optical-path motion.

For an optical phase delivered along path length L(t)L(t), mirror or fiber motion changes the phase at the atoms. Active path stabilization must sense the relevant endpoint. Stabilizing to a nearby mirror while leaving a moving segment between that mirror and the atoms can produce an interrogation-synchronous frequency error.

Neutral-atom density shifts are usually measured by interleaving populations and fitting frequency versus a density proxy. The proxy itself may depend on detection gain, cloud size, excitation fraction, temperature, and site occupancy. A trustworthy evaluation:

  1. defines the density observable;
  2. randomizes or balances the interleave order;
  3. checks linearity and excitation dependence;
  4. propagates uncertainty in the extrapolation point; and
  5. tests for correlations with lattice depth and loading.

Single-ion clocks avoid clock-ion density shifts but not background-gas collisions. Collision-induced phase steps can have non-Gaussian distributions, and energetic collisions can change the ion crystal without immediate loss. Event cuts require predeclared diagnostics and a bias test.

A clock servo can be biased by:

  • unequal sampling of the two line slopes;
  • cavity or laser drift during the modulation sequence;
  • different atom number or contrast at alternate points;
  • line-shape asymmetry and neighboring transitions;
  • digital quantization, saturation, or integrator leakage;
  • feed-forward error;
  • imperfect frequency steps during composite pulses; and
  • data rejection correlated with the measured detuning.

Useful tests include reversing probe order, changing servo gain and update law, injecting a known synthetic offset, locking with independent discriminators, and comparing against another clock. Agreement of an in-loop error signal with zero is expected for a working servo and is not an independent validation.

A typical optical-clock comparison may average approximately as

σy(τ)=aτ/s\sigma_y(\tau) = \frac{a}{ \sqrt{\tau/\mathrm{s}} }

over a range of τ\tau. This describes a noise process under the stated measurement conditions. It does not show that the mean is unbiased or that the uncertainty budget is complete.

Conversely, a systematic uncertainty of usysu_{\mathrm{sys}} is not a noise floor that every Allan-deviation curve must reach. A constant unknown offset does not contribute to Allan deviation, and correlated drift may average differently from white frequency noise. Reports should show stability analysis, corrected means, systematic ledgers, and comparison statistics separately.

The strongest evidence for an optical standard combines:

  • repeated evaluations over changed conditions and long times;
  • two clocks of the same species with substantially independent apparatus;
  • ratios between different species;
  • redundant combs or transfer chains;
  • remote comparisons through independently characterized links;
  • closure tests across a network;
  • blinded or precommitted analysis where feasible; and
  • uncertainty normalized residuals across laboratories.

Same-species agreement is sensitive to differential implementation errors but can miss shared atomic-theory or environmental assumptions. Different-species ratios probe different sensitivities but add a second clock’s systematic ledger. A mature evidence program uses both.

Performance numbers are moving targets and should always carry dates and conditions. Examples that illustrate architecture rather than establish a permanent ranking include:

  • optical lattice clocks have demonstrated total systematic uncertainties in the low-10−1810^{-18} range and below;
  • a 2025 27Al+^{27}\mathrm{Al}^{+} quantum-logic clock reported 5.5×10−195.5\times10^{-19} systematic uncertainty and instability 3.5×10−16/τ/s3.5\times10^{-16}/\sqrt{\tau/\mathrm{s}} using a one-second Rabi probe and remotely transferred laser stability; and
  • a three-species optical network reported Al+^{+}/Yb, Al+^{+}/Sr, and Yb/Sr ratios with fractional uncertainties at or below 8×10−188\times10^{-18} in 2021.

These results are not interchangeable. A systematic uncertainty, a comparison uncertainty, and a one-second instability answer different questions. They also depend on the reported coverage convention, correlations, data selection, and reference conditions.

In the weak-field, slow-motion limit, two stationary clocks at gravitational potentials differing by ΔU\Delta U have

Δνν≃ΔUc2.\frac{\Delta\nu}{\nu} \simeq \frac{\Delta U}{c^2}.

Near Earth’s surface, for a small vertical displacement in an approximately uniform field,

Δνν≃gΔhc2≃1.09×10−18(Δh1 cm).\frac{\Delta\nu}{\nu} \simeq \frac{ g\Delta h }{ c^2 } \simeq 1.09\times10^{-18} \left( \frac{\Delta h}{1\ \mathrm{cm}} \right).

The sign depends on the definition of ΔU\Delta U and which clock is in the numerator. With the usual convention, a clock at higher gravitational potential runs faster.

The centimeter conversion is a local approximation. A clock measures potential, not geometric height alone. Precise work includes:

  • Earth’s gravitational potential and its temporal variation;
  • centrifugal potential in the rotating terrestrial frame;
  • solid-Earth and ocean tides;
  • atmospheric and hydrological mass redistribution;
  • the reference potential used for conventional relativistic corrections;
  • coordinates of the atomic reference point; and
  • the potential difference along the comparison epoch.

Geodetic leveling and gravimetry remain necessary to interpret or validate the clock result. Optical clocks add a direct potential-difference observable and can connect separated height systems, monitor changes, or test geopotential models.

A local vertical comparison may use two clocks in one building connected by phase-stabilized fiber. A remote comparison can use:

  • fixed optical clocks linked through a fiber network;
  • a transportable clock brought to the field site;
  • free-space optical transfer;
  • a clock network tied through an intermediate transfer laser; or
  • hybrid optical and satellite techniques.

For corrected ratio RABR_{AB}, a schematic measurement model is

ln⁡RABobs=ln⁡RAB(0)+UA−UBc2+δlink+δclock,A−δclock,B.\ln R_{AB}^{\mathrm{obs}} = \ln R_{AB}^{(0)} + \frac{ U_A-U_B }{ c^2 } + \delta_{\mathrm{link}} + \delta_{\mathrm{clock},A} - \delta_{\mathrm{clock},B}.

RAB(0)R_{AB}^{(0)} is the ratio for a declared common reference potential. The link and both clock corrections must be estimated before interpreting the residual as geopotential.

A defensible clock-geodesy result states:

  1. which atomic reference points and epochs define the comparison;
  2. the frequency-ratio sign and gravitational-potential convention;
  3. the geodetic datum or reference potential;
  4. the clock, link, and geodetic uncertainty components and covariance;
  5. tidal and loading corrections;
  6. how transport, restart, or relocation shifts were tested; and
  7. whether the result estimates potential, geopotential number, or an equivalent height under a specified gravity model.

“The clock measured height” is shorthand only after these choices are declared.

Frequency ratios and dimensionless constants

Section titled “Frequency ratios and dimensionless constants”

Atomic transition frequencies depend differently on dimensionless constants because relativistic, hyperfine, nuclear, and many-electron contributions differ. For a ratio RAB=νA/νBR_{AB}=\nu_A/\nu_B, a linearized model is

δln⁡RAB=(Kα,A−Kα,B)δln⁡α+(Kμ,A−Kμ,B)δln⁡μ+⋯ .\delta\ln R_{AB} = \left( K_{\alpha,A} - K_{\alpha,B} \right) \delta\ln\alpha + \left( K_{\mu,A} - K_{\mu,B} \right) \delta\ln\mu + \cdots.

Here α\alpha is the fine-structure constant and, on this page,

μ≡mpme.\mu \equiv \frac{m_p}{m_e}.

Some literature defines μ=me/mp\mu=m_e/m_p, which reverses the associated sensitivity coefficient. Every analysis must state the convention.

The coefficients KK come from atomic, molecular, and nuclear-structure calculations. A measured ratio drift constrains a combination of parameters. It becomes a bound on one constant only after assuming the others are fixed or combining enough independent observables to estimate them.

A phenomenological redshift-violation model can assign each clock species a parameter βi\beta_i:

νi(U)−νi(0)νi(0)=(1+βi)Uc2.\frac{ \nu_i(U)-\nu_i(0) }{ \nu_i(0) } = \left( 1+\beta_i \right) \frac{U}{c^2}.

General relativity predicts βi=0\beta_i=0. A co-located ratio measured while the laboratory potential changes is sensitive to βi−βj\beta_i-\beta_j, not to a common universal term. Annual solar-potential modulation, terrestrial height comparisons, and spacecraft or tower experiments test different combinations of assumptions and systematics.

Environmental annual cycles can mimic a potential-correlated signal. Temperature, magnetic fields, humidity, electronics, link behavior, and data availability must be included in the nuisance model. Choosing the phase or period after inspecting the data also changes the statistical interpretation.

Some ultralight dark-matter models predict coherent oscillations of effective constants and therefore of clock ratios. A representative template is

δln⁡RAB(t)=Acos⁡(2πfϕt+φ),\delta\ln R_{AB}(t) = A \cos \left( 2\pi f_\phi t+\varphi \right),

where AA depends on a specified field density, coupling model, and sensitivity coefficients. Networks can also search for transient, correlated disturbances.

A null clock residual is not a model-independent exclusion of dark matter. A valid limit requires:

  • a Lagrangian or explicit coupling parameterization;
  • assumptions about field density, coherence time, and velocity distribution;
  • clock response and sampling transfer functions;
  • calibrated timing and network geometry;
  • treatment of gaps, colored noise, and systematics;
  • trials factors for a frequency or event-time scan; and
  • a declared confidence or credible-interval construction.

Clock networks can be exceptionally sensitive within such models because ratios are precise and geographically separated instruments provide correlation tests. The interpretation remains conditional.

Clock transitions can test Lorentz symmetry by searching for orientation- or boost-dependent frequency changes as Earth rotates and orbits. The observable may involve Zeeman sublevels with anisotropic electronic momentum sensitivity rather than the field-insensitive component optimized for timekeeping. Magnetic fields, tensor Stark shifts, quadrupole shifts, laboratory orientation, and sidereal-versus-solar aliases become central.

The result constrains coefficients in a specified effective theory and frame convention. It should not be summarized as a universal “test of Lorentz invariance” without naming the sector and coefficient combination.

Building a Trustworthy Optical-Clock Result

Section titled “Building a Trustworthy Optical-Clock Result”

A mature optical-clock campaign should preserve enough information to reconstruct the result:

  1. Define the measurand. State isotope, transition, Zeeman or hyperfine averaging, rest-frame convention, reference potential, and epoch.
  2. Freeze sign conventions. Document comb beats, AOM shifts, ratio order, correction signs, and relativistic potential signs.
  3. Record raw cycle data. Preserve excitation, atom number, servo state, environmental monitors, validity flags, and timing.
  4. Calibrate transfer functions. Measure discriminator slope, servo response, optical-path response, and counter behavior.
  5. Interleave systematics. Randomize or balance changed conditions and preserve the complete sequence.
  6. Propagate covariance. Keep common coefficients and calibrations correlated across days and clocks.
  7. Blind high-impact choices. Where practical, hide the final ratio or inject a secret offset while cuts and models are settled.
  8. Compare independently. Use another clock, comb, link, species, or analysis path.
  9. Audit residuals. Examine dependence on time, operating coordinates, data quality, and environmental variables.
  10. Publish evidence. Report enough metadata, code, covariance, and corrected data for an informed reanalysis.

The atomic transition is the reference. The running output comes from a laser, servo, synthesis chain, and counter. Noise or bias in any of them can move the delivered frequency.

Linewidth determines resolution and potential stability. Accuracy requires an evaluated shift model and comparison evidence.

Treating the magic wavelength as exact and universal

Section titled “Treating the magic wavelength as exact and universal”

The operational zero can depend on lattice depth, polarization, motional state, hyperpolarizability, multipolar terms, and fit convention.

Assuming symmetric-state averaging cancels everything

Section titled “Assuming symmetric-state averaging cancels everything”

Opposite Zeeman states cancel the leading odd-in-mm term under matched sampling. Quadratic Zeeman, tensor Stark, drift, line pulling, and unequal weights remain.

Using room temperature instead of radiation temperature

Section titled “Using room temperature instead of radiation temperature”

A sensor value is not automatically the effective temperature seen by the atoms. View factors, gradients, windows, ovens, and emissivity matter.

Counters display magnitudes. Frequency equations require signed optical ordering and every pass through every modulator.

Calling an optical ratio an absolute frequency

Section titled “Calling an optical ratio an absolute frequency”

A ratio is dimensionless and can avoid the microwave chain. An absolute frequency in SI hertz must connect to the caesium-defined second.

Two accurate clocks at different potentials should disagree before the relativistic correction. At 10−1810^{-18}, centimeters matter.

Reading a model-independent discovery into a ratio residual

Section titled “Reading a model-independent discovery into a ratio residual”

Variation-of-constants, redshift-violation, dark-matter, and Lorentz tests depend on sensitivity coefficients, effective-theory conventions, environmental models, and statistical trials.

  1. Bureau International des Poids et Mesures, The International System of Units (SI Brochure), 9th ed., version 4.01 (June 2026), doi:10.59161/AUEZ1291.
  2. Bureau International des Poids et Mesures, “Recommended values of standard frequencies: for the second,” updated following the 2025 adjustment, BIPM standard-frequency list.
  3. H. S. Margolis et al., “Data associated with the 2025 adjustment of standard frequencies” (2026), doi:10.59161/StdFreq2026.
  4. Bureau International des Poids et Mesures, “Roadmap to the redefinition of the second,” updated March 2025, BIPM redefinition portal.
  5. N. Dimarcq et al., “Roadmap towards the redefinition of the second,” Metrologia 61, 012001 (2024), doi:10.1088/1681-7575/ad17d2.
  6. T. M. Fortier, A. N. Luiten, and H. S. Margolis, “Optical atomic clocks: defining the future of time and frequency metrology,” Optica 13, 143–163 (2026), doi:10.1364/OPTICA.575770.
  7. A. D. Ludlow, M. M. Boyd, J. Ye, E. Peik, and P. O. Schmidt, “Optical atomic clocks,” Reviews of Modern Physics 87, 637–701 (2015), doi:10.1103/RevModPhys.87.637.
  8. H. Katori, M. Takamoto, V. G. Pal’chikov, and V. D. Ovsiannikov, “Ultrastable optical clock with neutral atoms in an engineered light shift trap,” Physical Review Letters 91, 173005 (2003), doi:10.1103/PhysRevLett.91.173005.
  9. A. Derevianko and H. Katori, “Colloquium: Physics of optical lattice clocks,” Reviews of Modern Physics 83, 331–347 (2011), doi:10.1103/RevModPhys.83.331.
  10. P. O. Schmidt et al., “Spectroscopy using quantum logic,” Science 309, 749–752 (2005), doi:10.1126/science.1114375.
  11. S. M. Brewer et al., “27Al+^{27}\mathrm{Al}^{+} quantum-logic clock with a systematic uncertainty below 10−1810^{-18},” Physical Review Letters 123, 033201 (2019), doi:10.1103/PhysRevLett.123.033201.
  12. M. C. Marshall et al., “High-stability single-ion clock with 5.5×10−195.5\times10^{-19} systematic uncertainty,” Physical Review Letters 135, 033201 (2025), doi:10.1103/hb3c-dk28.
  13. T. L. Nicholson et al., “Systematic evaluation of an atomic clock at 2×10−182\times10^{-18} total uncertainty,” Nature Communications 6, 6896 (2015), doi:10.1038/ncomms7896.
  14. W. F. McGrew et al., “Atomic clock performance enabling geodesy below the centimetre level,” Nature 564, 87–90 (2018), doi:10.1038/s41586-018-0738-2.
  15. E. Oelker et al., “Demonstration of 4.8×10−174.8\times10^{-17} stability at 1 s1\ \mathrm{s} for two independent optical clocks,” Nature Photonics 13, 714–719 (2019), doi:10.1038/s41566-019-0493-4.
  16. T. Bothwell et al., “Resolving the gravitational redshift across a millimetre-scale atomic sample,” Nature 602, 420–424 (2022), doi:10.1038/s41586-021-04349-7.
  17. M. Takamoto et al., “Test of general relativity by a pair of transportable optical lattice clocks,” Nature Photonics 14, 411–415 (2020), doi:10.1038/s41566-020-0619-8.
  18. J. Grotti et al., “Geodesy and metrology with a transportable optical clock,” Nature Physics 14, 437–441 (2018), doi:10.1038/s41567-017-0042-3.
  19. C. Lisdat et al., “A clock network for geodesy and fundamental science,” Nature Communications 7, 12443 (2016), doi:10.1038/ncomms12443.
  20. Boulder Atomic Clock Optical Network Collaboration, “Frequency ratio measurements at 18-digit accuracy using an optical clock network,” Nature 591, 564–569 (2021), doi:10.1038/s41586-021-03253-4.
  21. S. M. Middelmann, C. Lisdat, S. Falke, J. S. R. Vellore Winfred, F. Riehle, and U. Sterr, “Tackling the blackbody shift in a strontium optical lattice clock,” IEEE Transactions on Instrumentation and Measurement 60, 2550–2557 (2011), doi:10.1109/TIM.2010.2088470.
  22. S. G. Porsev and A. Derevianko, “Multipolar theory of blackbody radiation shift of atomic energy levels and its implications for optical lattice clocks,” Physical Review A 74, 020502(R) (2006), doi:10.1103/PhysRevA.74.020502.
  23. R. C. Brown et al., “Hyperpolarizability and operational magic wavelength in an optical lattice clock,” Physical Review Letters 119, 253001 (2017), doi:10.1103/PhysRevLett.119.253001.
  24. M. Bishof et al., “Inelastic collisions and density-dependent excitation suppression in a 87Sr^{87}\mathrm{Sr} optical lattice clock,” Physical Review A 84, 052716 (2011), doi:10.1103/PhysRevA.84.052716.
  25. T. Rosenband et al., “Frequency ratio of Al+^{+} and Hg+^{+} single-ion optical clocks; metrology at the 17th decimal place,” Science 319, 1808–1812 (2008), doi:10.1126/science.1154622.
  26. M. S. Safronova et al., “Search for new physics with atoms and molecules,” Reviews of Modern Physics 90, 025008 (2018), doi:10.1103/RevModPhys.90.025008.
  27. N. Huntemann, B. Lipphardt, C. Tamm, V. Gerginov, S. Weyers, and E. Peik, “Improved limit on a temporal variation of mp/mem_p/m_e from comparisons of Yb+^{+} and Cs atomic clocks,” Physical Review Letters 113, 210802 (2014), doi:10.1103/PhysRevLett.113.210802.
  28. F. Riehle, “Optical clock networks,” Nature Photonics 11, 25–31 (2017), doi:10.1038/nphoton.2016.235.
  29. T. Mehlstäubler, G. Grosche, C. Lisdat, P. O. Schmidt, and H. Denker, “Atomic clocks for geodesy,” Reports on Progress in Physics 81, 064401 (2018), doi:10.1088/1361-6633/aab409.
  30. C. Sanner et al., “Optical clock comparison for Lorentz symmetry testing,” Nature 567, 204–208 (2019), doi:10.1038/s41586-019-0972-2.
  31. M. Filzinger et al., “Improved limits on the coupling of ultralight bosonic dark matter to photons from optical atomic clock comparisons,” Physical Review Letters 130, 253001 (2023), doi:10.1103/PhysRevLett.130.253001.

1. Optical quality factor and interrogation

Section titled “1. Optical quality factor and interrogation”

A strontium lattice clock operates at ν0=4.29228×1014 Hz\nu_0=4.29228\times10^{14}\ \mathrm{Hz} and uses a one-second Rabi pulse. For an estimate, take the interrogation-limited linewidth to be Δν=0.80/T\Delta\nu=0.80/T.

  1. Estimate the linewidth.
  2. Find the corresponding observed quality factor.
  3. A metastable lifetime would imply a natural linewidth of 1 mHz1\ \mathrm{mHz}. Find the natural quality factor.
  4. Explain which quality factor should be used to describe the observed discriminator.
Solution

With T=1.00 sT=1.00\ \mathrm{s},

Δν=0.80T=0.80 Hz.\Delta\nu = \frac{0.80}{T} = 0.80\ \mathrm{Hz}.

The observed quality factor is

Qobs=4.29228×10140.80=5.37×1014.Q_{\mathrm{obs}} = \frac{ 4.29228\times10^{14} }{ 0.80 } = 5.37\times10^{14}.

For a natural linewidth of 10−3 Hz10^{-3}\ \mathrm{Hz},

Qnat=4.29228×101410−3=4.29×1017.Q_{\mathrm{nat}} = \frac{ 4.29228\times10^{14} }{ 10^{-3} } = 4.29\times10^{17}.

The observed discriminator is set by the actual pulse, laser coherence, contrast, and decoherence, so QobsQ_{\mathrm{obs}} is the relevant realized value. QnatQ_{\mathrm{nat}} describes the transition’s lifetime-limited potential and should not be presented as the measured linewidth.

A clock laser is above comb tooth n=2 145 678n=2\,145\,678. The comb has

frep=200.000000 MHz,fCEO=20.000000 MHz.f_{\mathrm{rep}} = 200.000000\ \mathrm{MHz}, \qquad f_{\mathrm{CEO}} = 20.000000\ \mathrm{MHz}.

The positive beat magnitude is fb=35.000000 MHzf_b=35.000000\ \mathrm{MHz}. The light then passes once through a +80.000000 MHz+80.000000\ \mathrm{MHz} AOM before reaching the atoms.

  1. Write the signed frequency equation at the atoms.
  2. Calculate the optical frequency.
  3. Find the error caused by using the wrong beat-note sign.
  4. Find the error caused by forgetting the AOM.
Solution

“Above the tooth” gives sb=+1s_b=+1, and the AOM is explicitly a positive shift:

νatom=nfrep+fCEO+fb+fAOM.\nu_{\mathrm{atom}} = n f_{\mathrm{rep}} + f_{\mathrm{CEO}} + f_b + f_{\mathrm{AOM}}.

The tooth contribution is

nfrep=2 145 678(200.000000 MHz)=429 135 600.000000 MHz.n f_{\mathrm{rep}} = 2\,145\,678 \left( 200.000000\ \mathrm{MHz} \right) = 429\,135\,600.000000\ \mathrm{MHz}.

Therefore

νatom=429 135 600.000000 MHz+20.000000 MHz+35.000000 MHz+80.000000 MHz=429 135 735.000000 MHz=429 135 735 000 000 Hz.\begin{aligned} \nu_{\mathrm{atom}} &= 429\,135\,600.000000\ \mathrm{MHz} \\ &\quad + 20.000000\ \mathrm{MHz} + 35.000000\ \mathrm{MHz} + 80.000000\ \mathrm{MHz} \\ &= 429\,135\,735.000000\ \mathrm{MHz} \\ &= 429\,135\,735\,000\,000\ \mathrm{Hz}. \end{aligned}

Changing +fb+f_b to −fb-f_b changes the answer by −2fb=−70 MHz-2f_b=-70\ \mathrm{MHz}. Forgetting the AOM changes it by −80 MHz-80\ \mathrm{MHz}. Both errors are vastly larger than any optical-clock uncertainty.

3. Residual lattice shift near the magic frequency

Section titled “3. Residual lattice shift near the magic frequency”

Near an operational point, suppose the measured fractional lattice shift is modeled by

δlatt=a(νL−νm)U+bU2,\delta_{\mathrm{latt}} = a \left( \nu_L-\nu_m \right) U + bU^2,

where νL−νm\nu_L-\nu_m is in megahertz and UU is the trap depth in units of ErE_r. Let

a=1.0×10−20(MHz Er)−1,b=−2.0×10−22Er−2.a = 1.0\times10^{-20} \left( \mathrm{MHz}\,E_r \right)^{-1}, \qquad b = -2.0\times10^{-22} E_r^{-2}.

Evaluate the shift at νL−νm=3.0 MHz\nu_L-\nu_m=3.0\ \mathrm{MHz} and U=50ErU=50E_r. Then find the frequency detuning that makes this two-term model zero at the same depth.

Solution

At the stated point,

δlatt=(1.0×10−20)(3.0)(50)−(2.0×10−22)(50)2=1.5×10−18−5.0×10−19=1.0×10−18.\begin{aligned} \delta_{\mathrm{latt}} &= \left( 1.0\times10^{-20} \right) (3.0)(50) \\ &\quad - \left( 2.0\times10^{-22} \right) (50)^2 \\ &= 1.5\times10^{-18} - 5.0\times10^{-19} \\ &= 1.0\times10^{-18}. \end{aligned}

For zero shift at fixed UU,

a(νL−νm)U+bU2=0.a \left( \nu_L-\nu_m \right) U + bU^2 = 0.

Thus

νL−νm=−bUa=−(−2.0×10−22)(50)1.0×10−20=1.0 MHz.\nu_L-\nu_m = - \frac{bU}{a} = - \frac{ (-2.0\times10^{-22})(50) }{ 1.0\times10^{-20} } = 1.0\ \mathrm{MHz}.

The result illustrates an operational magic point: cancellation at one depth need not identify the zero of the E1 polarizability alone.

At 300.0 K300.0\ \mathrm{K}, a clock has a fractional BBR correction of −5.0×10−15-5.0\times10^{-15}. Treat the shift as exactly proportional to T4T^4 and let the standard uncertainty of the effective radiation temperature be u(T)=0.050 Ku(T)=0.050\ \mathrm{K}.

  1. Find the fractional uncertainty contributed by temperature.
  2. How small must u(T)u(T) be for this contribution to be 1.0×10−181.0\times10^{-18}?
  3. Name three omitted contributions in a real BBR evaluation.
Solution

The relative sensitivity of the shift is

u(∣δBBR∣)∣δBBR∣=4u(T)T.\frac{u(|\delta_{\mathrm{BBR}}|)} {|\delta_{\mathrm{BBR}}|} = 4 \frac{u(T)}{T}.

Therefore

u(δBBR)=(5.0×10−15)40.050300.0=3.33×10−18.u(\delta_{\mathrm{BBR}}) = \left( 5.0\times10^{-15} \right) 4 \frac{0.050}{300.0} = 3.33\times10^{-18}.

For a target of 1.0×10−181.0\times10^{-18},

u(T)=(1.0×10−18)(300.0)4(5.0×10−15)=0.015 K.u(T) = \frac{ (1.0\times10^{-18})(300.0) }{ 4(5.0\times10^{-15}) } = 0.015\ \mathrm{K}.

Thus the effective radiation temperature would need a standard uncertainty of about 15 mK15\ \mathrm{mK} under this simplified model. Real evaluations also include differential polarizability, dynamic correction, gradients, view factors, surface emissivities, hot apertures, and sensor calibration and self-heating.

An ion’s total rms speed during interrogation is estimated as vrms=0.20 m s−1v_{\mathrm{rms}}=0.20\ \mathrm{m\,s^{-1}}.

  1. Find the fractional time-dilation shift.
  2. Find the absolute shift of a 1.121×1015 Hz1.121\times10^{15}\ \mathrm{Hz} transition.
  3. Why can the result not be assigned solely to the secular temperature without further evidence?
Solution

Using c=2.99792458×108 m s−1c=2.99792458\times10^8\ \mathrm{m\,s^{-1}},

Δνν=−vrms22c2=−(0.20)22(2.99792458×108)2=−2.23×10−19.\frac{\Delta\nu}{\nu} = - \frac{v_{\mathrm{rms}}^2}{2c^2} = - \frac{(0.20)^2} {2(2.99792458\times10^8)^2} = -2.23\times10^{-19}.

The absolute shift is

Δν=(−2.23×10−19)(1.121×1015 Hz)=−2.50×10−4 Hz.\Delta\nu = \left( -2.23\times10^{-19} \right) \left( 1.121\times10^{15}\ \mathrm{Hz} \right) = -2.50\times10^{-4}\ \mathrm{Hz}.

The velocity variance can include secular motion, intrinsic micromotion, excess micromotion, driven motion, and heating during the dark interval. A secular temperature alone does not determine all of these components.

Clock AA is 12.0 cm12.0\ \mathrm{cm} above clock BB in a laboratory where g=9.80 m s−2g=9.80\ \mathrm{m\,s^{-2}}. Ignore tides and use a uniform field.

  1. Find (νA−νB)/ν(\nu_A-\nu_B)/\nu.
  2. If the measured corrected ratio residual is (1.45±0.20)×10−17(1.45\pm0.20)\times10^{-17}, is it consistent with the height model at one standard uncertainty?
  3. What equivalent height corresponds to a fractional uncertainty of 3.0×10−193.0\times10^{-19}?
Solution

The potential difference is

ΔU=gΔh=(9.80)(0.120)=1.176 m2 s−2.\Delta U = g\Delta h = (9.80)(0.120) = 1.176\ \mathrm{m^2\,s^{-2}}.

Hence

Δνν=ΔUc2=1.176(2.99792458×108)2=1.31×10−17.\frac{\Delta\nu}{\nu} = \frac{\Delta U}{c^2} = \frac{1.176} {(2.99792458\times10^8)^2} = 1.31\times10^{-17}.

The residual difference from the model is

(1.45−1.31)×10−17=0.14×10−17=1.4×10−18.\left( 1.45-1.31 \right) \times10^{-17} = 0.14\times10^{-17} = 1.4\times10^{-18}.

This is 0.700.70 times the stated standard uncertainty 2.0×10−182.0\times10^{-18}, so it is consistent at one standard uncertainty.

The equivalent local height uncertainty is

u(h)=c2gu ⁣(Δνν)=(2.99792458×108)29.80(3.0×10−19)=2.75×10−3 m.u(h) = \frac{c^2}{g} u\!\left( \frac{\Delta\nu}{\nu} \right) = \frac{ (2.99792458\times10^8)^2 }{ 9.80 } (3.0\times10^{-19}) = 2.75\times10^{-3}\ \mathrm{m}.

Thus the simplified conversion is about 2.8 mm2.8\ \mathrm{mm}. A real result is an uncertainty in potential; conversion to height requires a gravity and geodetic model.

7. Sensitivity to variation of the fine-structure constant

Section titled “7. Sensitivity to variation of the fine-structure constant”

Two clock transitions have Kα,A=6.0K_{\alpha,A}=6.0 and Kα,B=0.0K_{\alpha,B}=0.0. Suppose their measured ratio drift is

ddtln⁡RAB=(1.2±1.8)×10−18 yr−1.\frac{d}{dt}\ln R_{AB} = \left( 1.2\pm1.8 \right) \times10^{-18}\ \mathrm{yr^{-1}}.

Assume for this exercise that only α\alpha varies.

  1. Infer dln⁡α/dtd\ln\alpha/dt and its standard uncertainty.
  2. Explain why this is not a model-independent statement about all constants.
  3. If an unmodeled annual temperature systematic exists, what additional evidence would be useful?
Solution

Under the one-parameter assumption,

ddtln⁡RAB=(Kα,A−Kα,B)ddtln⁡α.\frac{d}{dt}\ln R_{AB} = \left( K_{\alpha,A} - K_{\alpha,B} \right) \frac{d}{dt}\ln\alpha.

Therefore

ddtln⁡α=(0.20±0.30)×10−18 yr−1.\frac{d}{dt}\ln\alpha = \left( 0.20\pm0.30 \right) \times10^{-18}\ \mathrm{yr^{-1}}.

The inference assumes other dimensionless parameters are constant and the sensitivity coefficients are correct. More general models require additional independent clock and molecular ratios to separate parameter combinations.

Useful evidence includes measured temperature sensitivities, environmental regression fixed before looking at the final drift, multiple years of data, comparison with another laboratory and climate, residual spectra, and ratios with different KαK_\alpha. An annual signal alone can alias solar potential, temperature, humidity, and data-availability cycles.

8. Design and audit an optical-clock comparison

Section titled “8. Design and audit an optical-clock comparison”

Design one of the following:

  • a local Al+^{+}/Yb optical ratio;
  • a remote Sr/Sr relativistic-geodesy comparison;
  • a transportable Sr clock campaign;
  • a Yb+^{+}/Sr test of constant variation; or
  • a same-species comparison of two independent lattice clocks.

Specify:

  1. the measurand and ratio orientation;
  2. clock transition and interrogation for each endpoint;
  3. clock-laser and comb architecture;
  4. at least six platform-specific systematic effects per clock;
  5. path, link, and counter diagnostics;
  6. gravitational-potential treatment;
  7. covariance and data-selection model;
  8. at least three reversal or interleaving tests; and
  9. the evidence that would distinguish a clock discrepancy from a link or analysis discrepancy.
Solution

There is no unique solution. For a remote Sr/Sr geodesy comparison, a defensible design uses two independently evaluated 87Sr^{87}\mathrm{Sr} lattice clocks and defines

RAB=νAνBR_{AB} = \frac{\nu_A}{\nu_B}

at a common conventional reference potential. Each laboratory operates a cavity-prestabilized clock laser, symmetric mFm_F averaging, phase-stable delivery to the atoms, and an atomic servo. A coherent telecom transfer laser and self-referenced comb at each endpoint bridge the clock wavelength to a bidirectionally stabilized fiber.

Each clock ledger includes BBR, lattice E1 detuning, hyperpolarizability, multipolar and motional lattice shifts, quadratic Zeeman shift, density shift, probe shift, line pulling, optical-path phase, servo error, and background-gas collisions. The link ledger includes local uncompensated fibers, bidirectional asymmetry, interferometer offsets, amplifiers, cycle slips, comb transfer, AOM signs, counter dead time, and gate synchronization.

The geodetic model states both atomic reference points, comparison epochs, tidal and loading corrections, centrifugal potential, datum, and covariance between survey and gravity observations. Useful reversals include high/low lattice depth, high/low density, opposite Zeeman components, reversed probe order, alternate combs, and a deliberately changed link path.

To localize a discrepancy, preserve simultaneous local ratio beats, round-trip link diagnostics, independent comb channels, raw clock excitations, and common-view intervals. A transportable third clock or a loop with a second link can add closure. Agreement after swapping combs but not after bypassing a local fiber implicates the local transfer branch; agreement locally but not remotely implicates the link or geopotential model; disagreement between independently corrected local clocks implicates a clock ledger or shared optical reference. The conclusion must follow these diagnostics rather than the preferred application.