# Atomic Combination Clocks

###### Abstract

Atomic clocks use atomic transitions as frequency references. The susceptibility of the atomic transition to external fields limits clock stability and introduces systematic frequency shifts. Here, we propose to realize an atomic clock that utilizes an entangled superposition of states of multiple atomic species, where the reference frequency is a sum of the individual transition frequencies. The superposition is selected such that the susceptibilities of the respective transitions, in individual species, destructively interfere leading to improved stability and reduced systematic shifts. We present and analyze two examples of such combinations. The first uses the optical quadrupole transitions in a Ca - Yb two-ion crystal. The second is a superposition of optical quadrupole transitions in one Sr ion and three Hg ions. These combinations have reduced susceptibility to external magnetic fields and blackbody radiation.

Atomic clocks are composed of three main components. An electromagnetic oscillator, an atomic reference to which the oscillator locks; i.e. a narrow atomic transition, and a frequency counter which translates the number of oscillator periods to a time unit. The atomic reference is usually carefully chosen to have as little as possible coupling to varying conditions in the environment which can lead to systematic frequency shifts. These systematic shifts can be generally divided into two categories. The first is due to the susceptibility of the atomic transition to different kinds of environmental electromagnetic fields such as the ambient magnetic field or blackbody radiation Rosenband et al. (2006); Safronova and Safronova (2011); Dubé et al. (2013). The second kind of clock shifts are relativistic and emanate from the fact that the atoms are not investigated in an inertial frame. Here examples include second-order Doppler shifts or the gravitational red-shift Chou et al. (2010a). Relativistic shifts can be only mitigated by maintaining the atoms in an inertial frame. The first kind however can be reduced not only by trying to minimize background electromagnetic fields, but also by reducing the differential susceptibility of the clock transition to different electromagnetic perturbations. The blackbody radiation shift, for example, can be reduced by choosing a particular atomic transition with a small differential polarizability Rosenband et al. (2006). Alternatively, it was recently proposed to suppress blackbody radiation shift by using two different atomic clocks and a frequency comb to generate a “synthetic” frequency which is largely insensitive to the presence of electric fields Yudin et al. (2011).

Reducing the susceptibility of atomic superpositions to electromagnetic perturbations has been the focus of intensive research in the context of quantum computing. Here, the motivation is to improve the coherence of multi-qubit superpositions for the purpose of metrology. One of the most successful methods for prolonging the coherence times of superpositions is the use of decoherence-free subspaces Lidar et al. (1998). In this approach entangled superpositions of different states have reduced susceptibility to noise, leading to multi-second coherence times Roos et al. (2006); Kotler et al. (2014).

We propose to use several atomic transitions in different species, each ticking with it own environmentally-sensitive frequency. These transitions are chosen such that their susceptibilities compensate each other in an entangled superposition leading to significantly reduced sensitivity to environmental fields. The idea of combining different materials with opposite susceptibilities in clocks dates back to the work of John Harrison in the eighteenth century Usher (2013). Harrison revolutionized the fields of timekeeping and navigation by building a marine chronometer based on a bi-metallic balance spring to compensate for temperature variations of the spring constant.

The phase of a multi-species entangled state will evolve at a frequency which is a linear sum of the different frequencies of individual transitions. It is thus possible to construct an atomic clock superposition that will evolve at a frequency of an ultra-violate transition, however, no ultra-violate laser will be involved in the clock operation. Instead, the superposition phase would be compared to the linear sum of frequencies of the clock lasers of the different transitions using an optical frequency comb. Fig.1(A) shows a schematics of our method. Inquiry the two-specie superposition phase with an optical frequency comb serves as a ”mechanical differential” in analogy with the device used in vehicles which allows two wheels to rotate at different rotation frequencies while keeping the average rotation frequency fixed.

Maximally entangled multi-atom states can, in principle, be used to surpass the frequency estimation accuracy attainable with uncorrelated atoms (standard quantum limit) and reach the Heisenberg limit Wineland et al. (1994). However, as pointed by Huelga et.al Huelga et al. (1997), in practice, clock stability is ultimately limited by the finite coherence time of the system. For uncorrelated Markovian noise the two cases will give the same result. Nevertheless, in the some cases of correlated noise, one can sill gain by cleverly engineering the proper entangled state Roos et al. (2006); Kessler et al. (2014); Louchet-Chauvet et al. (2010).

## I The Combination clock

We begin by describing our general scheme and discuss in detail the quantum protocol that can be used in the case of a two-ion combination clock . Our scheme is a generalization of the Ramsey spectroscopy method similar to the case of identical ions in a maximally entangled state Bollinger et al. (1996); Leibfried et al. (2004).

An -ion, multi-species, atomic register is prepared in a superposition of two states . Where and are tensor products of states in the ground and excited state manifolds of the different ions. Here and can be either the ground or excited states of the clock transition in ion of species , and are the quantum numbers of the exact state in the ground or excited states manifolds of the transition (e.g. Zeeman or hyperfine state). The clock frequency is , where and is the energy of the state. Notice that with our definition can be negative if ion is prepared such that is the excited state and is the ground state of the clock transition. The total frequency is thus a linear sum of the clock transition frequencies of all the ions in the crystal.

After an interrogation time the phase between the atomic superposition and corresponding linear sum of local oscillators (i.e. clock lasers), , is obtained in the following way. A carrier pulse is applied to all ions in the crystal with their corresponding clock lasers followed by state detection and parity analysis. The probability of measuring even parity (i.e. an even number of ions in the excited state) will be given by . A following change to the frequency of one of the clock lasers can lock the linear sum of frequencies to the entangled superposition.

Comparing with the case of a single clock-ion species, the use of multi ion-species increases experimental complexity due to the need of having a set of all the lasers, including narrow-linewidth clock lasers, that are necessary for each of the ion-species. On the other hand, some of the technical aspects such as single-ion addressing and detection are simplified owing to the spectral distinguishably between species. Comparing with logic spectroscopy schemes Schmidt et al. (2005); Chou et al. (2010b), the combination clock does not involve motion-sensitive operations after the initial preparation stage and therefore continuous cooling or very-low heating rates are not a prerequisite for implementing this scheme.

The experimental sequence for the specific example of a two species two-ion crystal is depicted in Fig.1(B). The entanglement of two different ion species was recently demonstrated Tan et al. (2015); Ballance et al. (2015). As a simple example, we begin with the two ions in the ground state of their internal state as well as in the ground state of one of their normal modes of motion . After applying a blue-sideband (BSB) -pulse to ion A, we entangled the internal state of ion A with the motion of the two-ion crystal. Next we apply a red-sideband (RSB) -pulse on ion B which entangles the internal states of the ions while disentangling the internal states from motion. In the rotating-frame with respect to the sum of atomic transition frequencies , the resulting state is,

(1) | ||||

Here with as the detunings of the two lasers from their corresponding atomic transitions and is a constant phase for zero interrogation time resulting from finite pulse times, light shifts, etc. The two-ion entangled state acquires a phase which is the sum of the two transition frequencies times the interrogation time . This phase difference with respect to the sum of both clock laser frequencies is measured by a carrier pulse on each ion with its own clock laser, followed by state detection and parity analysis.

The combined two-ion clock is equivalent to a single atom with and with measured error signal . In the case of two uncorrelated lasers it is enough to feedback on only one of the lasers while making sure that each of the lasers is sufficiently close to resonance such that it can execute the necessary pulses for the protocol; i.e. to within their Rabi frequency. One way to experimentally retrieve the stable clock frequency is by sum frequency generation using a nonlinear crystal. A more practical way is to work with correlated lasers which may be derived from a frequency comb. In this case we can write the two lasers frequencies as and where and are integers and is the comb repetition rate. Here the comb offset frequency is assumed to be set to zero. The error signal from the experimental interrogation is then fed-back to stabilize to the value that locks . The extension of the locking of the comb to the linear sum in the case of more than two species is straightforward. .

Below we analyze two examples where the combination of two optical transitions with differential scalar polarizabilities of opposite sign such that the BBR shift is considerably reduced. This is one of the more difficult systematic shifts to mitigate because, for a given polarizability, it requires the evaluation of the BBR fields with high accuracy. As shown below, other differential shifts due to electromagnetic moments can be either reduced (e.g. first order Zeeman) or enhanced (e.g. electric quadrupole shifts). It should be mentioned that this method is not limited to the combination of optical transitions only. It is also possible to combine optical clock transitions with a microwave hyperfine transition in an entangled superposition. Here the resulting effective clock frequency is very similar to that of the optical component, however the magnetic susceptibility of the combination superposition can be controlled, especially when working at a high magnetic field in the Paschen-Back regime. In this regime the effective g-factor of the hyperfine transition can be smoothly tuned to null the superposition first order magnetic susceptibility.

We note that the use of entangled states is not always necessary. In many cases one can start with non-entangled states which have a significant overlap with the desired entangled state. If the part orthogonal to the desired clock superposition dephases quickly it will yield no parity signal. In this case, the measurement of parity will yield the same information as before, however at the cost of lower contrast and thus lower signal-to-noise ratio Chwalla et al. (2007). The advantage here would be the lower technical complexity as an entangling operation will not be necessary.

## Ii Ca - Yb combination clock

The first combination clock we examine is composed of the electric quadrupole transition in Ca and the same electric quadrupole transition; ; in Yb. This seems like a particularly nice combination since the two transitions are driven using fairly accessible wavelengths of and correspondingly; while the effective frequency at which this clock superposition rotates corresponds to an effective wavelength of 273 nm (1090 THz) and lifetime of ms limited by the Yb transition. This effective wavelength is very close to that of the Al optical clock (267 nm) with the advantage that no UV light has to be used in order to drive this transition. Furthermore the magnetic moments of these two transitions are the same to within a fractional difference of , so superpositions with states containing an equal number of states with opposite projections along the magnetic field direction can have almost zero magnetic susceptibility. But most importantly, the differential static polarizability of these two transitions are very similar in magnitude but have opposite signs Ludlow et al. (2015) and therefore a superposition of the two species can have a significantly smaller susceptibility to blackbody radiation. Up to normalization, the specific superposition we consider here is,

(2) |

Below is an analysis of the different systematic shifts for this combination clock.

### ii.1 First-order Zeeman

The first-order Zeeman shift due to a uniform magnetic field for this superposition is almost zero as both superposition parts have a total zero magnetic moment, up to a small () difference between the bound electron -factor in the two cases. Magnetic-field gradients, however, will still contribute a first order sensitivity to magnetic field. Here the magnetic field gradient can be measured, and subsequently nulled, by preparing the superposition,

(3) |

A similar superposition has been previously used to estimate and compensate magnetic field gradients to the 20 mHz level Kotler et al. (2014) which translates into a fractional uncertainty of . One can further improve on this uncertainty by interlacing experiments in which the spin direction of both ions is reversed, thus time-averaging out the gradient effect in the same manner that first order Zeeman shift is commonly averaged over Bernard et al. (1998).

### ii.2 Second order Zeeman

Other optical clocks such as Yb benefit from the possibility of a clock transition which is first order insensitive to the magnetic field. In most of these cases this happens due to the hyperfine interaction between the electron and the nucleus. This first-order insensitivity then comes at the expense of a relatively large second-order Zeeman shift (typically or more larger than in the first-order sensitive case). Here, since the first order insensitivity results from careful engineering of the clock superposition, the clock superposition amplitudes do not depend on the magnetic field value. A change in the magnetic field will not result in a change of the magnetic moment of the superposition states and hence no second order Zeeman shift. The small second order Zeeman shift will be, as in first order sensitive clocks, due to the fine-structure coupling in each of the transitions. The order of the second order Zeeman will be,

(4) |

where is the electron dimensionless magnetic moment, is Bohr magneton, is the magnetic field, is Plank’s constant and is the fine-structure (FS) splitting in the D manifold. The FS splitting in Ca is THz and in Yb it is 42 THz. This means that the second order Zeeman will be dominated by Ca. Working at a field of T we get,

(5) |

One can measure the average magnetic field by preparing the superposition,

(6) |

Here it is likely that the magnetic field will be measured with at least accuracy, allowing for a clock frequency error as low as 4 mHz which translates into a fractional uncertainty below .

### ii.3 Electric Quadrupole shift

The quadrupole shift (QS) of this clock superposition will be due to the quadrupole moment of the level of both ions since the orbital has no electric quadrupole moment. The quadrupole shift of level is given by Itano (2006),

(7) |

Here, is the level quadrupole moment and is the angle between the quantization axis set by the magnetic field and the quadrupole axis set by the trap geometry, where we assume full radial (cylindrical) symmetry. In the case of two ions in the trap, the equilibrium positions of the ions are , symmetrically positioned around the potential minimum along the trap axial direction. This results in an electric field gradient that is three times larger than the case of a single ion. Using and taking the case where as a worst case scenario, the QS for the will be

(8) |

This result is rather large due to the constructive addition of the entangled superposition and the position of the ions in the trap. In order to calibrate the QS one can place the ions in the superposition,

(9) |

Which, in the absence of magnetic moment evolves solely due to the QS Roos et al. (2006). In this case the QS will be factor of two larger as compared with that of the clock superposition above. Here a value of can be chosen to minimize (null) the QS. Moreover the QS can be eliminated by averaging over different magnetic sub-levels Dubé et al. (2013). For the case of the D level, only the m and m states are required.

### ii.4 Blackbody radiation shift

The blackbody shift is due to the scalar differential polarizability of the clock transition; . The blackbody radiation creates a varying electric field at the ion position . The shift of the clock transition is given by,

(10) |

Here is the differential static polarizability between the two clock transition states. The dynamic correction factor takes into account the effect of BBR fields not being static but rather peaked at around (30 THz) at 300 K. This dynamic correction factor is typically small () and we will therefore neglect it in the following. The effective differential static polarizability of the superposition above (for the purpose of frequency shift calculation) will therefore be the sum of the two polarizabilities of the two transitions involved.

The differential polarizability for the above transition in Ca is [] Arora et al. (2007); Safronova and Safronova (2011). This value is similar to the Ca clock transition which is investigated as an optical clock in several labs Chwalla et al. (2009); Wolf et al. (2011); Matsubara et al. (2012); Huang et al. (2016). The corresponding differential polarizability on the same clock transition in Yb is [] Schneider et al. (2005). The differential polarizability of the Ca-Yb combination clock is []. This is a significant reduction in the expected BBR shift of at least a factor of if not more. The expected shift to the clock transition at is corresponding to a fractional uncertainty of dominated by the uncertainty in the differential polarizability. With better measurement of the this polarizability and a accuracy in estimating the BBR field power, the BBR uncertainty of this clock will be in the range.

### ii.5 Second order Stark shift due to trap fields

Trap fields will induce frequency changes due to the differential polarizability of the transition. Trap fields are mostly due to uncompensated micromotion or inherent micromotion that results from ion thermal motion in the pseudo-potential. Here, since the trap fields have a well-defined direction both the scalar and the tensor part of the differential polarizability will play a role.

To the extent that the trap rf fields are identical on both ions then the scalar polarizability of the transition is suppressed by the same factor that is mentioned in the BBR section. A difference in the rf amplitude the two ions experience is likely. Firstly, any stray field orthogonal to the trap axis will tilt the two ion crystal and displace the ions into different rf regions. This is because the radial pseudo-potential depends on mass and therefore the two ions experience different radial restoring forces. Secondly, any inhomogeneous axial rf fields will lead to a difference too. However, with the minimization of micromotion this residual shift should be at least as small as in the case of each of the two clocks when operated on a single ion.

The tensor polarizability of Ca is Safronova and Safronova (2011). The tensor polarizability of is Ludlow et al. (2015). Therefore, the combined tensor polarizability of this clock transition is . The amplitude of rf fields that the ions experience can be minimized to the level. The order of magnitude of tensor shifts expected is therefore in the few corresponding to fractional uncertainty in the level. Furthermore, the tensor shift depend on the level and magnetic field direction, here the projection of the electric field, rather than the electric field gradient along the magnetic field direction is used. This means that the same average of measurements along different magnetic field directions or different levels that nulls the quadrupole shift will null the tensor polarizability contribution as well.

## Iii Sr - 3Hg combination clock

The second combination we look into here is motivated primarily by the desire to reduce BBR shifts as much as possible by reducing the combination differential scalar polarizability. The differential polarizability Ludlow et al. (2015) of the electric quadrupole transition at in Sr is about 3 times larger and opposite in sign than the differential polarizability Ludlow et al. (2015) of the same transition in Hg at . Thus, a superposition of a single Sr and three Hg ions would significantly reduce to .

(11) | |||

The effective frequency at which this superposition rotates is () - deep in the extreme UV. The expected shift to the clock transition at T=300 K is below 20 mHz corresponding to a fractional uncertainty of . With better measurement of the differential polarizability and 10% accuracy in estimating the BBR field power, the BBR uncertainty of this clock will be in the range.

The rest of the systematic shift analysis here is rather similar to the one presented above for the Ca - Yb combination clock. In short: The first-order Zeeman shift due to a uniform magnetic field for this superposition is zero as both superposition parts have a total zero magnetic moment. Moreover, with this particular spatial ordering of spin states, there would not be first order sensitivity to magnetic field gradients as well and only second order spatial magnetic field variation would cause a shift. The FS splitting in Hg is THz and in Sr it is THz and therefore second order Zeeman shifts will be dominated largely by Sr and can therefore be estimated using Eq. 4 to be Hz. Here as well, an entangled superposition can be used to measure the average magnetic field on the clock array. Evaluating it at the level, leads to a systematic shift of the order of . The quadrupole shift in this case will be rather large due to the large crystal and is estimate to be around 5 kHz for the same conditions considered above.

In conclusion, we propose and analyze the use of entangled states composed of multi ion-species for realizing an optical clock with reduced systematic shifts and improved stability. As concrete examples we proposed two different combinations. A Ca - Yb two ions clock and a Sr - 3Hg four-ion crystal. Both these examples are expected to operate with a significantly reduced BBR shift. Our approach opens a new dimension to choosing suitable atomic clock references. Instead of merely scanning the periodic table for a suitable transition in different atomic species, we propose to search for the suitable linear combination of transitions, in different species which would have the lowest susceptibilities to environmental conditions.

## References

- Rosenband et al. (2006) T. Rosenband, W. M. Itano, P. Schmidt, D. Hume, J. Koelemeij, J. C. Bergquist, and D. J. Wineland, in Frequency and Time Forum (EFTF), 2006 20th European (IEEE, 2006) pp. 289–292.
- Safronova and Safronova (2011) M. S. Safronova and U. I. Safronova, Phys. Rev. A 83, 012503 (2011).
- Dubé et al. (2013) P. Dubé, A. A. Madej, Z. Zhou, and J. E. Bernard, Phys. Rev. A 87, 023806 (2013).
- Chou et al. (2010a) C.-W. Chou, D. Hume, T. Rosenband, and D. Wineland, Science 329, 1630 (2010a).
- Yudin et al. (2011) V. Yudin, A. Taichenachev, M. Okhapkin, S. Bagayev, C. Tamm, E. Peik, N. Huntemann, T. Mehlstäubler, and F. Riehle, Phys. Rev. Lett. 107, 030801 (2011).
- Lidar et al. (1998) D. A. Lidar, I. L. Chuang, and K. B. Whaley, Phys. Rev. Lett. 81, 2594 (1998).
- Roos et al. (2006) C. F. Roos, M. Chwalla, K. Kim, M. Riebe, and R. Blatt, Nature 443, 316 (2006).
- Kotler et al. (2014) S. Kotler, N. Akerman, N. Navon, Y. Glickman, and R. Ozeri, Nature 510, 376 (2014).
- Usher (2013) A. P. Usher, A history of mechanical inventions: revised edition (Courier Corporation, 2013).
- Wineland et al. (1994) D. J. Wineland, J. J. Bollinger, W. M. Itano, and D. J. Heinzen, Phys. Rev. A 50, 67 (1994).
- Huelga et al. (1997) S. F. Huelga, C. Macchiavello, T. Pellizzari, A. K. Ekert, M. B. Plenio, and J. I. Cirac, Phys. Rev. Lett. 79, 3865 (1997).
- Kessler et al. (2014) E. M. Kessler, P. Komar, M. Bishof, L. Jiang, A. S. Sørensen, J. Ye, and M. D. Lukin, Phys. Rev. Lett. 112, 190403 (2014).
- Louchet-Chauvet et al. (2010) A. Louchet-Chauvet, J. Appel, J. J. Renema, D. Oblak, N. Kjaergaard, and E. S. Polzik, New J. Phys. 12, 065032 (2010).
- Bollinger et al. (1996) J. Bollinger, W. Itano, D. Wineland, and D. Heinzen, Phys. Rev. A 54, R4649 (1996).
- Leibfried et al. (2004) D. Leibfried, M. D. Barrett, T. Schaetz, J. Britton, J. Chiaverini, W. M. Itano, J. D. Jost, C. Langer, and D. J. Wineland, Science 304, 1476 (2004).
- Schmidt et al. (2005) P. O. Schmidt, T. Rosenband, C. Langer, W. M. Itano, J. C. Bergquist, and D. J. Wineland, Science 309 (2005).
- Chou et al. (2010b) C. W. Chou, D. B. Hume, J. C. J. Koelemeij, D. J. Wineland, and T. Rosenband, Phys. Rev. Lett. 104, 1 (2010b).
- Tan et al. (2015) T. R. Tan, J. P. Gaebler, Y. Lin, Y. Wan, R. Bowler, D. Leibfried, and D. J. Wineland, Nature 528, 380 (2015).
- Ballance et al. (2015) C. Ballance, V. Schäfer, J. Home, D. Szwer, S. Webster, D. Allcock, N. Linke, T. Harty, D. A. Craik, D. Stacey, et al., Nature 528, 384 (2015).
- Chwalla et al. (2007) M. Chwalla, K. Kim, T. Monz, P. Schindler, M. Riebe, C. Roos, and R. Blatt, Applied Physics B 89, 483 (2007).
- Ludlow et al. (2015) A. D. Ludlow, M. M. Boyd, J. Ye, E. Peik, and P. Schmidt, Rev. Mod. Phys. 87, 637 (2015).
- Bernard et al. (1998) J. Bernard, L. Marmet, and A. Madej, Opt. Comm. 150, 170 (1998).
- Itano (2006) W. Itano, Phys. Rev. A 73, 022510 (2006).
- Arora et al. (2007) B. Arora, M. S. Safronova, and C. W. Clark, Phys. Rev. A 76, 064501 (2007).
- Chwalla et al. (2009) M. Chwalla, J. Benhelm, K. Kim, G. Kirchmair, T. Monz, M. Riebe, P. Schindler, A. S. Villar, W. Hänsel, C. F. Roos, R. Blatt, M. Abgrall, G. Santarelli, G. D. Rovera, and P. Laurent, Phys. Rev. Lett. 102, 023002 (2009).
- Wolf et al. (2011) A. L. Wolf, J. Morgenweg, J. C. J. Koelemeij, S. A. van den Berg, W. Ubachs, and K. S. E. Eikema, Opt. Lett. 36, 49 (2011).
- Matsubara et al. (2012) K. Matsubara, H. Hachisu, Y. Li, S. Nagano, C. Locke, A. Nogami, M. Kajita, K. Hayasaka, T. Ido, and M. Hosokawa, Opt. Express 20, 22034 (2012).
- Huang et al. (2016) Y. Huang, H. Guan, P. Liu, W. Bian, L. Ma, K. Liang, T. Li, and K. Gao, Phys. Rev. Lett. 116, 013001 (2016).
- Schneider et al. (2005) T. Schneider, E. Peik, and C. Tamm, Phys. Rev. Lett. 94, 230801 (2005).